基于系统动力学的人口预测

3.2基于系统动力学的人口预测 21世纪是人类面临三大问题：第一是人口膨胀，第二是就业困难，第三是环境污染，这三大问题的焦点在于人口。因此，如何对未来的人口进行预测和控制，一直是人们关心的重要领域。 本课题是在宋健人口模型的基础上，考虑到上海作为一个开放城市，改良建立了双线性开放/动态人口模型。采用上述基于人口结构模型，预测上海2010—2050年的人口年龄、性别结构。为了更准确地研究人口系统，我们将人口按0-4岁、5-9岁、10-14岁、…、95-99岁、100岁及以上分群，分为21个群，并假设女性的生育时间以不同的概率分布在15-49岁之间。然后以政策系数和生育时间的分布概率为政策参数进行仿真分析和政策试验。 3.2.1系统模拟的一些基本假设 ● 人口分年龄数据 2000年人口普查的数据上海常住人口总数为1640万，而根据上海统计年鉴2000年上海常住人口总数为1608万。因为后续计算都是采用上海统计年鉴上的数据，所以按上海统计年鉴的常住人口总数1608万对2000年人口普查的数据进行了同比例调整。通过《上海统计年鉴》、《上海市2000年人口普查资料》、《2005年上海市1%人口抽样调查资料》等文献的搜索，得2000年上海市分年龄段的男、女人数数据见表1。 表1 上海市2000年分年段男、女人数（单位：万人）   第五次普查(2000年)   第五次普查(2000年) 年龄段 合计 男 女 年龄段 合计 男 女 0~4岁 49 25.7 23.3 55~59岁 66.61 34.3 32.31 5~9岁 58.67 30.53 28.14 60~64岁 57.73 28.54 29.19 10~14岁 93.42 47.95 45.48 65~69岁 68.31 32.25 36.06 15~19岁 129.51 64.9 64.62 70~74岁 53.13 24.37 28.76 20~24岁 135.3 69.72 65.57 75~79岁 36.71 15.94 20.77 25~29岁 135.4 72.09 63.31 80~84岁 19.39 7.67 11.73 30~34岁 139.48 76.16 63.32 85~89岁 7.9 2.71 5.19 35~39岁 150.65 81.27 69.38 90~94岁 2.16 0.61 1.55 40~44岁 170.02 88.98 81.04 95~99岁 0.42 0.11 0.32 45~49岁 159.15 82.89 76.26 100岁及以上 0.03 0.004 0.02 50~54岁 107.8   56.33 51.46                         ● 妇女生育时间 根据人口生育的一般规律可知，对出生有贡献的只有15-49岁的女性人口。出生率受人口政策的影响，如果严格实行“一对夫妇一个孩”的人口政策，那么任何一个女性在一生中只能生育一次。我们假设生育时间是在15-49岁之间均匀分布，于是有出生率=1/35≈2.9%。通过对统计资料和参考文献的整理和分析，可得妇女生育时间到俄分布规律如表所示。 表 2 妇女生育时间的分布规律 生育时间 15-19岁 20-24岁 25-29岁 30-34岁 35-39岁 40-44岁 45-49岁 概率 0.05 0.30 0.40 0.15 0.05 0.04 0.01                 ● 性别比 性别比是一个统计数据，是指新生婴儿中男性人口与女性人口的比例。新出生的人口可能是男性，也可能是女性。在自然出生的情况下，男性和女性的概率都是50%。但是根据前面的分析，新生婴儿中，男性与女性的平均性别比为105：100。 ● 政策系数 政策系数是一个政策参数，表明生育政策执行的严格程度。如果严格执行“一对夫妇一个孩”的人口政策，政策系数=1，随着执行程度的放松，其值增加。例如，如果实施“一对夫妇两个孩”的人口政策，政策系数=2。 ● 男、女性出生速率 根据政策系数，有 男性出生速率=“女性15-49”*出生率*（性别比）/（100+性别比）*政策系数； 女性出生速率=“女性15-49”*出生率*100/（100+性别比）*政策系数。 ● 死亡率 但不同年龄组死亡率存在差异。0-10岁组是少年儿童阶段，死亡率呈下降趋势，10-14岁组死亡率水平为最低，以后随着年龄的增长，死亡率逐步上升。由于上海市2008年男性预期寿命为79.06岁，女性预期寿命为83.50岁，人均寿命已经达到较高的水平，接近许多世界发达国家的水平，上升的空间已经不是很大，故在未来若干年中死亡率减低的速度必然逐步减弱。以2000年男性、女性死亡率为基期我们假设截止2050年上海人均死亡率每十年分别较上一个十年下降10%。 表3  上海市分年龄死亡率对比分析 1990年（‰） 1995年（‰） 2000年（‰） 2005年（‰） 0-4岁 2.88 0.939 1.1 0.98 5-9岁 0.32 0.298 0.24 0.07 10-14岁 0.33 0.375 0.21 0.23 15-19岁 0.45 0.608 0.25 0.05 20-24岁 0.54 0.59 0.36 0.08 25-29岁 0.61 0.927 0.38 0.21 30-34岁 0.8 0.938 0.51 0.42 35-39岁 1.18 0.88 0.87 0.61 40-44岁 1.52 2.303 1.43 1.05 45-49岁 2.36 2.293 1.99 1.66 50-54岁 4.25 3.641 3.01 2.05 55-59岁 7.04 6.254 4.73 3.43 60-64岁 11.71 10.874 8.21 6.03 65-69岁 19.89 17.217 14.66 11.76 70-74岁 34.96 29.4 27.13 17.96 75-79岁 62.95 54.411 48.68 37.12 80-84岁 106.48 94.775 89.11 63.61 85-89岁 166.33 150.899 155.35 124.83 90-94岁 276.33 215.847 262.48 187.72 95-99岁 282.609 354.52 358.33 100岁以上 142.857 384.62           ● 净迁入人口 按照前面的分析，以相同的速率进行。 3.2.2人口系统模型 根据以上的分析，建立人口系统模型如下： 图1 各年龄段人口系统模型 公式说明： (001)    女性出生速率=("女性15-19岁"*"P15-19"+"女性20-24岁"*"P20-24"+"女性25-29岁"*"P25-29"+"女性30-34岁"*"P30-34"+"女性35-39岁"*"P35-39"+"女性40-44岁"*"P40-44"+"女性45-49岁"*"P45-49"    )/5*100/(100+性别比)*政策 (002)    "女性0-4岁"= INTEG (女性出生速率-F1-FD1,23.3) (003)    "女性5-9岁"= INTEG (F1-F2-FD2,30.53) (004)    "女性10-14岁"= INTEG (F2-F3-FD3,47.95) (005)    "女性15-19岁"= INTEG (F3-F4-FD4,64.9) (006)    "女性20-24岁"= INTEG (F4-F5-FD5+"女性20-24岁净迁入",69.72) (007)    "女性25-29岁"= INTEG (F5-F6-FD6+"女性25-29岁净迁入",72.09) (008)    "女性30-34岁"= INTEG (F6-F7-FD7+"女性30-34岁净迁入",76.16) (009)    "女性35-39岁"= INTEG (F7-F8-FD8+"女性35-39岁净迁入",81.27) (010)    "女性40-44岁"= INTEG (F8-F9-FD9+"女性40-44岁净迁入",88.98) (011)    "女性45-49岁"= INTEG (+F9-F10-FD10+"女性45-49岁净迁入",82.89) (012)    "女性50-54岁"= INTEG (F10-F11-FD11+"女性50-54岁净迁入",56.33) (013)    "女性55-59岁"= INTEG (F11-F12-FD12,34.3) (014)    "女性60-64岁"= INTEG (F12-F13-FD13,28.54) (015)    "女性65-69岁"= INTEG (F13-F14-FD14,32.25) (016)    "女性70-74岁"= INTEG (F14-F15-FD15,24.37) (017)    "女性75-79岁"= INTEG (F15-F16-FD16,15.94) (018)    "女性80-84岁"= INTEG (F16-F17-FD17,7.67) (019)    "女性85-89岁"= INTEG (F17-F18-FD18,2.71) (020)    "女性90-94岁"= INTEG (F18-F19-FD19,0.61) (021)    "女性95-99岁"= INTEG (F19-F20-FD20,0.11) (022)    女性100岁及以上= INTEG (F20-FD21,0.004) (023)    男性出生速率=("女性15-19岁"*"P15-19"+"女性20-24岁"*"P20-24"+"女性25-29岁"*"P25-29"+"女性30-34岁"*"P30-34"+"女性35-39岁"*"P35-39"+"女性40-44岁"*"P40-44"+"女性45-49岁"*"P45-49"    )/5*性别比/(100+性别比)*政策 (024)    "男性0-4岁"= INTEG (男性出生速率-M1-MD1,25.7) (025)    "男性5-9岁"= INTEG (M1-M2-MD2,28.14) (026)    "男性10-14岁"= INTEG (M2-M3-MD3,45.48)    (027)    "男性15-19岁"= INTEG (M3-M4-MD4,64.62)    (028)    "男性20-24岁"= INTEG (M4-M5-MD5+"男性20-24岁净迁入",    65.57) (029)    "男性25-29岁"= INTEG (M5-M6-MD6+"男性25-29岁净迁入",    63.31) (030)    "男性30-34岁"= INTEG (M6-M7-MD7+"男性30-34岁净迁入",    63.32) (031)    "男性35-39岁"= INTEG (M7-M8-MD8+"男性35-39岁净迁入",    69.38) (032)    "男性40-44岁"= INTEG (M8-M9-MD9+"男性40-44岁净迁入",    81.04) (033)    "男性45-49岁"= INTEG (    +M9-M10-MD10+"男性45-49岁净迁入",    76.26) (034)    "男性50-54岁"= INTEG (M10-M11-MD11+"男性50-54岁净迁入",51.46)