基于改進(jìn)PDE的人口預(yù)測建模與仿真

廖松　付鵬霖, 王哲浩　周鋼

摘要：此文基于我國人口政策的變動及第五、六次人口普查數(shù)據(jù)，指出了其他人口預(yù)測模型的缺陷并采用改進(jìn)PDE人口預(yù)測模型對我國人口數(shù)量和結(jié)構(gòu)進(jìn)行了建模和仿真。選取了2010年后的20年作為仿真的時間范圍，分別在不同的人口政策導(dǎo)致相應(yīng)生育率變化的情況下建立模型，按照月份進(jìn)行迭代。得到所需人口數(shù)據(jù)并進(jìn)行分析，最終獲得三種政策下對應(yīng)的2010-2030年人口數(shù)量變化和人口結(jié)構(gòu)變化情況。

關(guān)鍵詞：PDE模型;人口預(yù)測;生育率;死亡率;差分迭代法

中圖分類號：TP24 ? ? 文獻(xiàn)標(biāo)識碼：A

文章編號：1009-3044（2019）15-0266-06

Abstract： According to the changes in Chinas population policy and the data from the fifth and sixth population censuses，In this paper， the shortcomings of traditional population forecasting models are analyzed. and uses the improved PDE population forecasting model to model and simulate the population size and structure of China. The 2010-2030 is selected as the simulation interval， and the population model under the one-child policy， the individual two-child policy， and the comprehensive two-child policy is iteratively iterated monthly to obtain the required population data and analyze it， and finally adopt the proper three policies. Population size and structure changes in 2010-2030. The results show that the current policy of “individual second child” gradually transitioning to”full second child” is an ideal choice.

Key words： PDE model; population policy; fertility rate; mortality rate; difference iteration

1 引言

人口政策作為我國的基本國策之一，必須符合我國的國情民情，與環(huán)境、社會發(fā)展?fàn)顩r相統(tǒng)一[1]。幾十年里，計劃生育政策有效地降低了人口增長率，提升了人口素質(zhì)。另一方面，與之同時而來的生育率持續(xù)走低、男女比例失調(diào)、人口老齡化嚴(yán)重等問題已開始影響到我國發(fā)展。為了解決上述問題，近年來我國人口政策先后經(jīng)歷了“獨生子女”“雙獨二孩”“單獨二孩”“全面二孩”階段。本文在每一種人口政策下，都抽象出對應(yīng)的數(shù)學(xué)模型對人口總量和年齡結(jié)構(gòu)展開預(yù)測及仿真，分析了不同政策對中國人口未來變化的影響，為今后人口政策的修訂提供了信息支持，具有重要研究意義。

2 人口預(yù)測模型建立

出于保證統(tǒng)計數(shù)據(jù)的真實性、有效性考慮，本文采用的數(shù)據(jù)都是最近兩次的人口普查數(shù)據(jù)。為便于研究和仿真，本文提出以下假設(shè)：

（1）本文所采用的數(shù)據(jù)均準(zhǔn)確、有效;

（2）本文只考慮出生、死亡對人口變化的影響，不考慮人口遷移因素;

（3）在本文提供預(yù)測的年限內(nèi)，不考慮不可抗力和人為因素的干涉。

2.1 模型選取

當(dāng)前，可供人口預(yù)測的模型包括線性回歸模型、聯(lián)合國模型、平均增長量法、指數(shù)平滑法、平均增長率法、灰色系統(tǒng)G-M法以及人口-發(fā)展-環(huán)境分析模型（PDE）。

上述模型和方法中，線性回歸模型假定人口數(shù)量的歷史趨勢不變，這不符合實際。聯(lián)合國模型相對于用于全國性人口的研究，更適合用于城鎮(zhèn)/鄉(xiāng)村人口比例的研究。平均增長率法、平均增長率法和指數(shù)平滑法的脫離實際應(yīng)用環(huán)境，精準(zhǔn)度度較低。灰色系統(tǒng)G-M法缺乏統(tǒng)計檢驗，同樣不適合我國人口的預(yù)測[2]。

本文選取基于PDE的按月迭代，按年儲存的人口-發(fā)展-環(huán)境人口預(yù)測模型來預(yù)測我國人口數(shù)量和結(jié)構(gòu)的變化。

2.2 模型建立

PDE（人口-發(fā)展-環(huán)境）人口模型的核心思想是，第n年的人口是第n-1的人口經(jīng)過出生、死亡、遷入和遷出等變化過程后得到的。所以，如果已知某年人口數(shù)量和該年出生率、死亡率、遷入遷出率等數(shù)據(jù)，就可以通過迭代方法得到后面每一年的人口數(shù)量。

本文采用按月迭代和按年存儲的思路，即：本月人口數(shù)=上月人口數(shù)+上月出生數(shù)-上月死亡數(shù)+上月遷入-上月遷出，本月結(jié)果作為下個月的基數(shù)，再將一年12個月的數(shù)據(jù)匯總算作一年的變化，作為本年度的新出生人口。首先分別對模型內(nèi)部各個參數(shù)（如生育率、死亡率、性別比等）進(jìn)行合理預(yù)測分析得出數(shù)值，再根據(jù)往年人口總量進(jìn)行按月迭代運算，設(shè)置仿真的時間范圍，最終得到所需數(shù)據(jù)。

人口預(yù)測模型數(shù)學(xué)表達(dá)式如下：

2.3 模型參數(shù)分析

2.3.1生育率

事實上，由于出生人口的漏報，由初始始統(tǒng)計數(shù)據(jù)得到的總生育率嚴(yán)重偏低于真實水平。用人口漏報率作為一個因素，對總體的生育率進(jìn)行修正，修正后，可得到精確度較高的生育率。

研究表明，2010年第六次全國人口普查漏報率為0.12%，遠(yuǎn)低于第五次全國人口普查1.81%的漏報率[3]。計入漏報率的影響，可以將生育率數(shù)學(xué)表達(dá)式修正如下：

2.3.2 獨生子女政策因素

本方案假設(shè)我國維持獨生子女政策不變，時間為2010-2030年。該方案與近期的人口情況最為接近，同時該方案的預(yù)測結(jié)果很可能是嚴(yán)重的人口老齡化，人口數(shù)量在某一年份到達(dá)峰值后會迅速減少，對未來人口政策調(diào)整具有一定的參考價值。

2.3.3單獨二孩政策因素

本方案的提出背景為2013年底正式實行的“單獨二孩”政策。在此政策背景下，將人口增長率分為兩個階段，第一階段為2010-2013年，第二階段為2014-2030年。將“獨生子女”和“單獨二孩”政策下的增長率分別作為第一階段和第二階段人口增長率。第二階段的生育率以最近一次人口普查得到的生育率為參照，考慮“單獨二孩”政策的影響，得到單獨二孩政策下較為準(zhǔn)確的生育率公式：

2.3.4 全面二孩政策

2017年，生育第二個孩子的限制徹底解除。故本方案中，將人口增長率分為三個階段，第一階段為2010-2013年，第二階段為2014-2016年，第三階段為2017-2030年。將“獨生子女”“單獨二孩”“全面二孩”政策下的增長率分別作為第一階段、第二階段和第三階段的人口增長率。第二階段的生育率以第六次人口普查得到的生育率為基礎(chǔ)，第三階段的生育率以第二階段的生育率為基礎(chǔ)?？紤]“全面二孩”政策調(diào)整后的生育率，具體調(diào)整方案數(shù)學(xué)表達(dá)式如下：

2.3.5 死亡率

目前，死亡率預(yù)測模型包括隨機(jī)死亡模型，包括APC 模型、Cairns-Blake-Dowd 模型[4]、Lee-Carter 模型[5]及其擴(kuò)展模型。

Lee-Carter模型中，死亡率變化和年齡因子x和時間因子t有關(guān)?；谶@樣的思路可以構(gòu)建函數(shù)表示時間為t年齡為x歲的人群的死亡率，其數(shù)學(xué)表達(dá)式如下：

2.3.6 性別比因素

本文對人口預(yù)測的研究涉及只出生率和死亡率，性別比是對人口數(shù)據(jù)分析得到的一個結(jié)果，對人口數(shù)量預(yù)測問題本身無影響。

2.3.7 人口遷入與遷出因素

遷移人口量相對于國內(nèi)人口總量顯得過于小了，可忽略不計。外國人口遷入與我國遷出不在本文考慮之列。

3 人口模型仿真

3.1 仿真設(shè)計

本文以Microsolt Visual Studio 2017為開發(fā)平臺，以上述建立的人口預(yù)測模型為基礎(chǔ)，對不同時期人口政策（包括獨生子女政策、單獨二孩政策以及全面二孩政策）對未來人口數(shù)量、年齡結(jié)構(gòu)的影響進(jìn)行仿真設(shè)計。

為了獲取較高的仿真精度，迭代方案采用有限差分算法進(jìn)行按月迭代。預(yù)設(shè)仿真開始年份為2010年，截止年份為2030年，后臺支持仿真至2050年，但仿真結(jié)果在2035年以前的年份精確度更高。

仿真界面軟件具備以下功能：

（1）選擇不同人口政策，對選定人口政策下2010年至2030年前任意年份時間區(qū)間內(nèi)的人口數(shù)量進(jìn)行預(yù)測，顯示預(yù)測后的統(tǒng)計數(shù)據(jù)和圖表;

（2）選擇不同人口政策，對選定人口政策下2010至2030年前任意年份時間區(qū)間的人口結(jié)構(gòu)進(jìn)行預(yù)測，顯示仿真時間范圍內(nèi)的人口年齡結(jié)構(gòu)圖以及仿真截止年份的人口年齡結(jié)構(gòu)的圖表。

圖1為仿真程序流程圖。

3.2 軟件功能

軟件加載完畢后，進(jìn)入軟件運行主界面，可以看到本仿真平臺的各個主要功能模塊和仿真說明。本次仿真預(yù)設(shè)區(qū)間為2010-2030年，可以得到三種政策下人口特征的變化。主界面包括以下：仿真截止年份設(shè)置、人口政策選擇與切換、預(yù)測數(shù)據(jù)及相應(yīng)圖表顯示。按照程序提示，輸入數(shù)據(jù)，選擇功能，即可推演不同人口政策下的人口特征的演化。

3.3 仿真結(jié)果分析

3.3.1 人口數(shù)量分析

仿真得到我國2010-2030年的人口數(shù)量變化如下方的圖3和表5。獨生子女政策下的推演結(jié)果表明，我國人口總量將于2021年前后達(dá)到人口最高峰13.760億，之后人口數(shù)量開始加速減少，將在2030年前后降至13.491億;在單獨二孩政策下的推演結(jié)果表明，我國的人口總量從2010-2016年緩慢增長至13.60億，2016年起增速加快，于2020年后增速開始放緩趨于穩(wěn)定增長，至2030年，將達(dá)到14.377億;全面二孩政策下的仿真結(jié)果顯示，我國的人口增長與單獨二孩政策下的增長模式幾乎一致，在2020年前后，人口增速繼續(xù)加快并逐漸放緩，人口數(shù)量持續(xù)增加，到了2030年，將達(dá)到14.831億。

3.3.2 人口結(jié)構(gòu)分析

通過模擬三種政策下的人口變化，還能夠得到我國2010-2030年在不同人口政策下的人口結(jié)構(gòu)變化情況。

仿真結(jié)果說明，如果繼續(xù)實行獨生子女政策，2030年我國0-16歲的少年兒童占總?cè)丝诘谋戎亟抵?2.9%，老年人（60歲以上）的比重升至26.3%，人口結(jié)構(gòu)老齡化嚴(yán)重。單獨二孩政策下仿真結(jié)果顯示，2010-2030年間我國兒童人口比重為16%-18%，比重略微下降;老年人人口比重為12%-22%，上升將近一成，老齡化現(xiàn)象依然凸顯，但較獨生子女政策下的預(yù)測比重有了明顯改善。全面二孩下，即2010年后的20年，經(jīng)過獨生子女、單獨二孩、全面二孩政策后的人口預(yù)測仿真結(jié)果顯示，2010年后的20年，我國兒童人口比重由17.9%降至16.85%，在三種政策下，下降幅度是最小的;老年人人口比重增至25.5%，介于前兩者。見圖4、表2所示。

3.3.3 綜合分析

由于本次仿真中不同政策僅會對出生率造成影響，即模型的人口變化主要反映在新出生人口上。因此，在2010-2030年間每隔5年選擇0-15歲的人口數(shù)據(jù)進(jìn)行對比分析，如圖5所示。

由圖5可知，按時間變化，單獨二孩政策到全面二孩政策確實對新出生人口數(shù)量的穩(wěn)定增加，從而改善人口年齡結(jié)構(gòu)起到了相當(dāng)重要的作用。三種政策下人口數(shù)目和構(gòu)成特點各有不同，對其進(jìn)行橫向?qū)Ρ龋梢缘玫饺缦陆Y(jié)論：

（1）若繼續(xù)幾十年的獨生子女政策，某一時間后（可以認(rèn)為是2021年）我國的新生人口將發(fā)生銳減，同時老年人比重也將逐年升高，加深我國社會老齡化水平;

（2）若實行單獨二孩政策，中國新生人口數(shù)量平穩(wěn)增長，我國老人老齡化能夠得到一定改善;

（3）若實行全面二胎政策，我國嬰幼兒、少年兒童數(shù)量將增加，青壯年占人口比例將增加;

綜上所述，本文認(rèn)為目前實施單獨二孩政策并逐步向全面二孩過渡的政策是目前較為理想的選擇。

4 總結(jié)及展望

4.1 仿真平臺優(yōu)勢

本仿真平臺設(shè)置大量功能接口以便于模塊化設(shè)計，平臺可操作性和可觀性較強(qiáng);與此同時，還有部分功能接口冗余，使本仿真平臺在未來具備了功能擴(kuò)展的巨大潛力。

4.2 ?不足與展望

由于時間倉促，未能將人口預(yù)測模型進(jìn)一步細(xì)化具體至各省市自治區(qū)，從而得到區(qū)域性人口密度等人口數(shù)據(jù)，提升人口預(yù)測分析的覆蓋面和精確度。

本次仿真更多的是基于現(xiàn)有數(shù)據(jù)，對于環(huán)境和其他各方面擾動因素考慮不夠充分，在未來會進(jìn)一步做出改進(jìn)。

最后衷心希望本次仿真能夠為未來人口政策的修訂提供一定幫助。

參考文獻(xiàn)：

[1] 趙艷芳. 論我國實施的全面二胎人口政策[J]. 山西青年， 2016（20）.

[2]孟令國，李超令，胡廣.基于PDE模型的中國人口結(jié)構(gòu)預(yù)測研究[J].中國人口·資源與環(huán)境，2014，24（2）：132-141.

[3]王金營.2000 年中國第五次人口普查漏報評估及年中人口估計[J].人口研究，2003，27（5）：53-62.

[4] Cairns A，Blake D，Dowd K. A Two-factor Model for Stochastic Mortality with Parameter Uncertainty： Theory and Calibration [J].Journal of Risk and Insurance，2006，73（4）：687 -718.

[5]李志生，劉恒甲.Lee-Carter死亡率模型的估計與應(yīng)用——基于中國人口數(shù)據(jù)的分析[J].中國人口科學(xué)，2010（03）：46-56，111.

【通聯(lián)編輯：梁書】