劉慶山，胡學(xué)平，來(lái)江濤

（安慶師范大學(xué) 數(shù)理學(xué)院，安徽 安慶 246133）

人口問題事關(guān)社會(huì)發(fā)展大局，對(duì)人口問題的研究和對(duì)人口數(shù)量的預(yù)測(cè)具有重要的意義.安徽省是我國(guó)中部人口大省，由于經(jīng)濟(jì)社會(huì)發(fā)展相對(duì)沿海城市有一定差距，一直是人口流出大省.隨著中部崛起戰(zhàn)略的推進(jìn)，安徽省經(jīng)濟(jì)發(fā)展也進(jìn)入了快車道，人口情況發(fā)生了一系列的變化，關(guān)于這些變化的研究對(duì)安徽省相關(guān)政策的制定有重要的意義.國(guó)內(nèi)學(xué)者在對(duì)人口進(jìn)行預(yù)測(cè)時(shí)使用了不同的模型和方法，江衛(wèi)華[1]等利用馬爾薩斯模型和羅杰斯蒂模型對(duì)西藏人口進(jìn)行了預(yù)測(cè)，得出人口發(fā)展要和資源環(huán)境相協(xié)調(diào)的結(jié)論；羅榮桂[2]等使用BP 神經(jīng)網(wǎng)絡(luò)對(duì)長(zhǎng)江流域人口進(jìn)行預(yù)測(cè)研究，發(fā)現(xiàn)相對(duì)于指數(shù)平滑法和Logistic 法，BP 神經(jīng)網(wǎng)絡(luò)法在人口數(shù)量預(yù)測(cè)中精確度最高；石敏軍[3]等使用優(yōu)化的 Leslie 人口預(yù)測(cè)模型對(duì)我國(guó)人口進(jìn)行了預(yù)測(cè)，得出中國(guó)人口老齡化程度將迅速加深的結(jié)論；張娟[4]等使用動(dòng)力學(xué)模型對(duì)人口數(shù)量進(jìn)行預(yù)測(cè)，認(rèn)為如果生育政策不放開，將會(huì)導(dǎo)致兒童撫養(yǎng)比降低，老齡化加深，這將不利于經(jīng)濟(jì)發(fā)展.但是Logistic 模型沒有從出生和死亡等細(xì)分角度考慮，只是籠統(tǒng)地考慮阻滯作用，BP 神經(jīng)網(wǎng)絡(luò)計(jì)算太復(fù)雜，Leslie 模型沒有考慮人口流動(dòng)對(duì)人口數(shù)量的影響，而宋健人口發(fā)展方程是分年齡計(jì)算，同時(shí)考慮了人口出生、死亡和遷移對(duì)人口數(shù)量的影響，對(duì)人口數(shù)量變化的原因概括得比較全面.丁研[5]使用人口發(fā)展方程對(duì)包頭市的人口數(shù)量進(jìn)行了預(yù)測(cè)，指出應(yīng)把控制出生人口數(shù)量、穩(wěn)定低生育水平、全面提升人口素質(zhì)放在人口發(fā)展的首位；屠小明[6]等使用人口發(fā)展方程對(duì)我國(guó)人口數(shù)量和人口老齡化進(jìn)行了研究，指出我國(guó)人口在短期到中期內(nèi)還將繼續(xù)增長(zhǎng).由于人口流動(dòng)對(duì)現(xiàn)在人口數(shù)量變動(dòng)的影響越來(lái)越大，因此有的學(xué)者在宋健人口發(fā)展方程基礎(chǔ)上加入了凈流動(dòng)人口數(shù)量，使用移動(dòng)平均線擬合凈遷移人口數(shù)量.但是由于用于預(yù)測(cè)的數(shù)據(jù)是常住人口的口徑，不是戶籍人口，常住人口數(shù)量已經(jīng)去除了一些流動(dòng)人口的人數(shù)，所以在常住人口基礎(chǔ)上減去流出人口并不符合實(shí)際.

基于已有的研究，本文從將人口數(shù)量分性別計(jì)算、將遷移人口改為流動(dòng)變動(dòng)人口等方面對(duì)宋健人口發(fā)展方程[7]進(jìn)行了修正，用修正方程對(duì)安徽省人口數(shù)量和人口結(jié)構(gòu)進(jìn)行了預(yù)測(cè)，根據(jù)分析結(jié)果，提出了相關(guān)促進(jìn)人口發(fā)展的建議..

1 人口預(yù)測(cè)方程

1.1 宋健人口發(fā)展方程

由于連續(xù)形式的人口發(fā)展方程不方便計(jì)算，所以本文主要介紹人口發(fā)展方程的離散形式.令xi(t)為第t年滿i周歲但不滿i+1周歲的人口總數(shù)，Ψ(t)為第t年出生的嬰兒總數(shù)，則

式中：β(t)為第t年的平均婦女生育率；［a1,a2］為第t年育齡婦女的育齡區(qū)間；ki(t)為第t年滿i周歲但不滿i+1周歲的年齡組中女性人口的比例；hi(t)為生育模式，而且滿足.可得人口發(fā)展離散方程組為

式中：g i(t)為第t年滿i周歲但不滿i+1周歲的遷移人口數(shù)量；u i(t)為第t年滿i周歲但不滿i+1周歲的人口死亡率（其中u∞表示嬰兒死亡率）；下標(biāo)n為人口存活的最大年齡.

令bi(t)=（1-u∞(t)）（1-u0(t)）k i(t)hi(t)，引入向量和矩陣，向量形式表示為

1.2 人口發(fā)展方程的改進(jìn)

對(duì)人口發(fā)展方程進(jìn)行2個(gè)方面的處理：（1）方程（1）沒有區(qū)分性別，是把男性和女性放在一起分析.由于性別比例一直在變動(dòng)，而且在死亡率等很多方面都有差異，所以本文對(duì)性別進(jìn)行區(qū)分單獨(dú)分析.（2）方程（1）中g(shù)(t) 是遷移人口，由于本文的數(shù)據(jù)是常住人口數(shù)據(jù)，常住人口數(shù)據(jù)是在戶籍人口的基礎(chǔ)上已經(jīng)去除了流動(dòng)出去的人口，所以將遷移函數(shù)選擇為流動(dòng)變動(dòng)函數(shù)m(t)，指的是以2010 年的凈流出人口數(shù)為基準(zhǔn)，如果人數(shù)增加則為正，否則為負(fù).

得到改進(jìn)后的人口發(fā)展方程為

式中：下標(biāo)n為人口存活的最大年齡；(t)為第t年滿i周歲但不滿i+1周歲的男性人口總數(shù)；(t)為第t年滿i周歲但不滿i+1周歲的女性人口總數(shù)；(t)為第t年滿i周歲但不滿i+1周歲的男性流動(dòng)變動(dòng)人口總數(shù)；(t)為第t年滿i周歲但不滿i+1周歲的女性流動(dòng)變動(dòng)人口總數(shù)；(t)為第t年滿i周歲但不滿i+1周歲的男性人口死亡率；(t)為第t年滿i周歲但不滿i+1周歲的女性人口死亡率.

2 基于修正人口發(fā)展方程的安徽省人口預(yù)測(cè)

基于《安徽省統(tǒng)計(jì)年鑒》收集安徽省2010—2019 人口數(shù)據(jù).由于2010 年數(shù)據(jù)是第六次人口普查數(shù)據(jù)，有詳細(xì)的分年齡、分性別的死亡率和流動(dòng)變動(dòng)人口數(shù)量，從2011 年開始的抽樣調(diào)查人口數(shù)據(jù)只有整年的死亡率和流動(dòng)變動(dòng)人口數(shù)量，因此在使用人口發(fā)展方程計(jì)算時(shí)，需要確定生育模式函數(shù)h(i)、年平均死亡率函數(shù)u(t)、年平均流動(dòng)人口變動(dòng)函數(shù)m(t) 和總和生育率函數(shù)β(t) .確定u(t)，m(t) 后，再將方程（2）中的分性別、分年齡的人口死亡率和流動(dòng)變動(dòng)人口總數(shù)按比例變動(dòng)，就得到2011 年以后的相關(guān)數(shù)據(jù).

2.1 修正的人口發(fā)展方程中的參數(shù)函數(shù)

2.1.1 生育模式函數(shù) 根據(jù)宋健[8]的測(cè)算，標(biāo)準(zhǔn)化的生育模式函數(shù)滿足卡方分布，如果i代表年齡，則可得到生育模式函數(shù)

此函數(shù)中i為變量，n為未知參數(shù).為了確定未知參數(shù)，宋健同時(shí)通過函數(shù)計(jì)算，提出計(jì)算未知參數(shù)n的公式，即amax=a1+n-2，式中：amax為生育峰值年齡，a1為生育年齡區(qū)間最低年齡.根據(jù)2010 年安徽省統(tǒng)計(jì)年鑒的數(shù)據(jù)，峰值生育年齡為20 歲，生育最低年齡為15 歲，故取n=7.

2.1.2 年平均死亡率函數(shù) 以2010 年第六次全國(guó)人口普查數(shù)據(jù)為起始，再結(jié)合2011—2019 年人口抽樣調(diào)查數(shù)據(jù)，設(shè)t代表年份，分別采用非線性差值法和ARMA 模型擬合年平均死亡率函數(shù).

運(yùn)用Matlab 軟件，計(jì)算得到年平均死亡率插值函數(shù)為

由于年平均死亡率數(shù)據(jù)也是時(shí)間序列，所以另外選用ARMA 模型來(lái)模擬.R 語(yǔ)言程序運(yùn)行結(jié)果表明，此序列一階差分后的自相關(guān)圖是一階截尾（見圖1a），偏自相關(guān)圖是拖尾的（見圖1b），所以選擇ARIMA(0,1,1)來(lái)模擬.通過R 語(yǔ)言編程得出的擬合模型為

圖1 年平均死亡率一階差分序列的自相關(guān)和偏自相關(guān)

比較2種方法得到的預(yù)測(cè)值和實(shí)際死亡率之間的誤差，結(jié)果見表1.由表1可以看出，ARIMA(0,1,1)模型誤差更小，故本文選擇ARIMA(0,1,1)模型來(lái)預(yù)測(cè)死亡率.

表1 非線性差值法和ARIMA(0,1,1)模型死亡率函數(shù)的預(yù)測(cè)值和實(shí)際死亡率比較 （‰）

2.1.3 年平均流動(dòng)人口變動(dòng)函數(shù) 本次人口方程使用的數(shù)據(jù)是常住人口，雖然現(xiàn)在流動(dòng)人口數(shù)量對(duì)總?cè)丝谟绊懞艽?，但是常住人口?jì)算時(shí)已經(jīng)扣除了凈流出人口數(shù)，以安徽省2010 年人口普查流動(dòng)人口數(shù)為基準(zhǔn)，計(jì)算出相對(duì)2010 年的人口變動(dòng)數(shù)，運(yùn)用Matlab 軟件計(jì)算得出的插值函數(shù)為

使用R 語(yǔ)言編程繪制流動(dòng)人口變動(dòng)數(shù)曲線，圖像呈逐年增長(zhǎng)趨勢(shì)，可見序列并不平穩(wěn)，所以需要進(jìn)行差分處理.一階差分后圖像沒有明顯的趨勢(shì)，差分后序列平穩(wěn).其自相關(guān)圖是一階截尾（見圖2a），偏自相關(guān)圖是拖尾的（見圖2b），所以使用ARIMA(0,1,1)模型進(jìn)行擬合.通過R 語(yǔ)言編程得出的擬合模型為

圖2 流動(dòng)人口變動(dòng)數(shù)量一階差分序列的自相關(guān)和偏自相關(guān)

將2個(gè)不同擬合方法的預(yù)測(cè)值和實(shí)際的流動(dòng)人口變動(dòng)數(shù)量進(jìn)行比較，結(jié)果見表2.由表2可以看出，ARMA 的總體誤差更小，故本文選擇ARMA 模型預(yù)測(cè).

表2 非線性差值法和ARIMA(0,1,1)模型流動(dòng)人口變動(dòng)數(shù)量預(yù)測(cè)值和實(shí)際變動(dòng)人數(shù)比較（萬(wàn)）

2.1.4 總和生育率 總和生育率指某地區(qū)育齡婦女平均每人生育的子女?dāng)?shù).根據(jù)安徽省統(tǒng)計(jì)年鑒，從2010年開始，安徽省總和生育率一直在下降，從2010 年的1.94 下降到2015 年的1.24.在2015 年全面放開二胎后，總和生育率開始上升，從2015 年的1.24 上升到2017 年的1.54，2018—2019 年都一直保持在1.54，2020 年安徽省總和生育率還未公布，從全國(guó)總和生育率1.3 來(lái)看可能會(huì)降低.今年中央全面放開三胎，對(duì)總和生育率預(yù)計(jì)會(huì)有促進(jìn)作用.本文根據(jù)生育率的實(shí)際變化情況，按照低方案、中方案和高方案3種不同的總和生育率變化方案對(duì)人口數(shù)量進(jìn)行預(yù)測(cè).低方案認(rèn)為雖然中央放開了三胎，但是由于人們的生育意愿低，導(dǎo)致2020 年后總和生育率還是在緩慢地降低，假設(shè)總和生育率按照每年0.01 的趨勢(shì)逐年遞減；中方案認(rèn)為隨著三胎的放開，雖然人們生育意愿降低，但是部分三胎的出生使得總和生育率一直保持不變；高方案假設(shè)總和生育率保持每年0.01 的逐年遞增趨勢(shì).

2.2 預(yù)測(cè)結(jié)果

2.2.1 預(yù)測(cè)結(jié)果 按照總和生育率的高、中、低3個(gè)方案，利用R 語(yǔ)言編程預(yù)測(cè)死亡率和流動(dòng)人口變動(dòng)數(shù)，得到2018—2025 年人口數(shù)量、少兒系數(shù)和老年人口系數(shù)的預(yù)測(cè)數(shù)據(jù)（見圖3）

圖3 2018—2025 年人口數(shù)量、少兒系數(shù)和老年人口系數(shù)的預(yù)測(cè)

2.2.2 誤差比較 通過和實(shí)際人口總數(shù)的比較，計(jì)算出2018—2019 年預(yù)測(cè)人口總數(shù)的誤差率（見表3）.由表3可以看出，誤差率很低，本文修正的人口發(fā)展方程預(yù)測(cè)效果良好.

表3 2018—2019 年預(yù)測(cè)人口總數(shù)的誤差率

2.3 預(yù)測(cè)分析

通過誤差比較發(fā)現(xiàn)，本次構(gòu)建的修正人口發(fā)展方程預(yù)測(cè)效果良好.從預(yù)測(cè)的折線圖可以看出，按照低方案，總和生育率保持遞減，將會(huì)造成人口總量自2018 年開始下降，2020 年后下降速度加快，與此同時(shí)，老年人口系數(shù)快速上升，勞動(dòng)人口比例下降的同時(shí)養(yǎng)老負(fù)擔(dān)反而加重，不利于經(jīng)濟(jì)社會(huì)發(fā)展.按照中方案，總和生育率保持不變，總?cè)丝趶?018 年開始到2021 年緩慢下降，2021 年以后緩慢上升；少兒系數(shù)在2018—2019 年基本持平后開始上升，老年人口系數(shù)自2018 年開始一直較快上升，則若按照此方案，人口會(huì)緩慢增長(zhǎng)，但是老年人口系數(shù)增長(zhǎng)較快，勞動(dòng)力人口比例仍然處于下降態(tài)勢(shì).按照高方案，人口總量在2018—2019 年先緩慢增長(zhǎng)，經(jīng)歷了2019—2021 年3 年的緩慢下降后又較快增長(zhǎng)起來(lái)；老年人口系數(shù)只在2019—2020 年短暫增長(zhǎng)，其它時(shí)間均處于下降態(tài)勢(shì)，說(shuō)明了老年人口比重在減少，勞動(dòng)年齡人口比例在增長(zhǎng)，有利于經(jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展.

3 建議

為了緩解安徽省人口老齡化程度的快速加深，提高人口總和生育率是有效的辦法.2021 年中央放開了三胎，對(duì)提高人口總和生育率起到了積極的作用.但是由于經(jīng)濟(jì)發(fā)展帶來(lái)的育兒壓力以及觀念和生活方式的轉(zhuǎn)變，年輕人的生育意愿普遍偏低.首先，建議政府在執(zhí)行中央政策時(shí)多一些配套的政策，在減輕育兒壓力等方面出臺(tái)一些措施；其次，人口老齡化直接帶來(lái)老年人的生活和撫養(yǎng)等社會(huì)問題，而當(dāng)前老年事業(yè)的發(fā)展較緩慢，需要政府統(tǒng)籌規(guī)劃，建議政府在應(yīng)對(duì)人口老齡化問題時(shí)要及早布局，未雨綢繆；最后，建議政府改善人才環(huán)境，吸引人，留住人.安徽省仍然是人口流出大省，而且流動(dòng)人口逐年增加，人才流失嚴(yán)重.希望政府能出臺(tái)相關(guān)政策改善人才環(huán)境，留住人才.