一种评估模型科学严谨性的方法

  • 格式:doc
  • 大小:613.97 KB
  • 文档页数:24

以深圳市为例探讨洪灾损失预测研究的科学性与严谨性摘要本文针对2020年和2050年深圳市洪灾损失进行建模以对其进行科学的预测,从而探讨洪灾损失预测研究的科学性和严谨性。

我们通过参考资料确定了影响洪灾损失的主要因素:降雨量、海平面高度、地质沉降、深圳GDP水平。

由于未来的降雨量存在较大的随机性故在此忽略。

另外查找资料得知地质沉降大约平均为1mm/年,所以本文只需建立模型先对海平面高度和深圳GDP水平进行预测,进而由得到的结果预测洪灾损失。

模型一是基于logistic模型的预测模型。

我们根据历史数据运用matlab 编程得到GDP的散点图,其大致走势是呈曲线上升的,但增长到一定程度后GDP 增长率会逐年减小。

由此分析我们建立了logistic模型并结合实际对模型进行了改进,预测深圳2020年与2050年GDP分别为2.1961万亿与2.7582万亿。

然后运用残差分析对所得数据进行误差估计,验证了模型的可行性。

模型二是基于时间序列的海平面高度的预测模型。

首先通过matlab绘制的散点图可以看出海平面的大致走势是上升的,然后通过曲线拟合和实际分析我们确定运用趋势移动平均法来建立海平面高度预测模型。

通过matlab数据计算预测得到2020年和2050年海平面高度分别为143.28mm、322.53mm。

最后通过绘制预测值和实际值的走势曲线可以看出所建模型的预测结果有较高的可信度。

模型三是基于BP神经网络的洪灾损失预测模型。

由于洪灾损失数据存在较大的跳跃性,而BP神经网络又有较好地容错性和较高的学习能力,故我们运用BP神经网络建立模型,通过matlab计算得到的2020年和2050年洪灾损失的预测值分别为37866.7万元38796.6万元。

同时由原始数据与仿真数据的对比图可以看出该预测模型可以很好地对洪灾损失进行合理的预测。

关键词: logistic、时间序列、BP神经网络、曲线拟合、误差估计一:问题重述洪涝灾害是当今世界上主要的自然灾害,防治洪水灾害是世界各国普遍关注的问题,洪灾损失评估是防洪减灾领域的一项基础性工作。

据报载,国外有研究报告将广州预测为受洪灾损失最重的城市,也将深圳列为洪灾损失严重的城市(见附件1和附件2)。

我国有关专家和专业人士认为该报告结论与事实存在出入,因而怀疑其所用方法及支撑数据的正确性与准确性。

我们需要解决以下几个问题:1、收集深圳市的相关资料,通过数学建模的方法,分析经济合作与发展组织(OCED)研究报告(附件2)中可能存在的问题。

2、基于问题1对2020和2050深圳可能遭受的洪灾损失做出预测,同时对比评价我们的模型与研究报告所用模型的优缺点。

3、基于我们的研究结果,给普通百姓写一份不超过一页的建议书,说明研究报告和我们的结果是怎样得到的,并提出一些建议,使普通百姓能够正确对待信息时代所谓“科学结论”快速传播带来的问题,比如“预测”给人们带来的不确定性和焦虑感。

4、给深圳市政府写一份不超过一页的建议书,除了说明研究报告和我们的结果是怎样得到的、可信度如何以及市政府应该做什么等(包括后续研究应该做些什么)。

二:问题分析洪水灾害是我国发生频率高、危害范围广、对国民经济影响最为严重的自然灾害,亦是威胁人类生存的十大自然灾害之一。

本文我们先通过阅读经济合作与发展组织(OCED)研究报告(题给附件2),分析研究报告所用模型及方法,加上李平日教授对于该报告的反驳(题给附件3),我们对报告的合理性进行了评判,再预测洪灾给深圳市2020年和2050年带来的损失,从而对洪灾损失预测研究的科学性和严谨性进行探讨。

要对深圳2020年和2050年的洪灾损失进行预测,就必然联系到降雨量、海平面上升、地质沉降和GDP等一系列的问题,即洪灾损失是一个复杂多变量预测的问题。

故而我们通过对主要影响因素的筛选,可以进一步确定其所适用的预测理论和模型。

由于中长期预测会受到很多不确定因素的影响(比如极端灾难性气候),所以本文我们忽略降雨量对洪灾损失的影响,并假设重大洪灾发生的概率和政府抗洪投入维持相对稳定,虽然在短期预测中可能出现较大误差,但是有利于得到中长期的统计性预测结果。

其次,地质沉降并不具有连续变化的特征,短期内也并不在总体上表现下降,但对于长期预测而言,可以将其看做连续变化的时间序列,满足一次函数的特征,取2003年为初始参考值1mm,用历史平均下降速率约1mm每年对2020和2050年的地质沉降做出预测,简化模型,以此可得到2020年和2050年的地质沉降数据(相对2003年)分别为18mm和48mm。

再者,我们收集到了大量海平面高度和深圳市GDP变化数据,由得到的数据然后联系实际我们分析得知时间序列和logistic模型可以很好地分别对二者进行合理的预测,最后由于洪灾损失数据存在较大的跳跃性,而BP神经网络又有较好地容错性和较高的学习能力,我们确定运用神经网络知识建立模型再预测2020年和2050年的洪灾损失。

在以上基础下,通过分析模型预测结果,可以科学且通俗易懂地给老百姓提适当建议以安抚人心;对有关政府部门提出相应建议,以帮助其制定合理的防灾政策,来最小化地减少洪灾造成的损失。

三:符号说明A...................................................................饱和水平C...............................................................速度增长因子 ...................................................................形状因子(1) M..................................................第t年第一次移动平均值t(2) M..................................................第t年第一次移动平均值ty...........................................................第t年海平面高度ta........................................................... 平滑系数(截距)tb.............................................................平滑系数(斜率)tt..................................................................当前时期数T.............................................为由t至预测期的时期的间隔数ˆy..................................................第t年海平面高度的预测值tˆy.............................................第t年海平面高度预测值的均值tˆy.............................................第t年海平面高度预测值的均值tQ..........................................调整激励函数形式的Sigmoid 参数L...................................................................学习次数Y ...................................................................样本个数 (1k Y)k x ≤≤....................................................... .样本数据 (1k Y)k d ≤≤.........................................................期望输出 P ...............................................................隐含节点个数 (1j )j net P ≤≤...................................................隐含节点输入 qk d ..............................................样本q 的节点k 的输出期望值 qk z ..............................................样本q 的节点k 的实际输出值 q E ................................................................误差函数值四:模型假设(1)假设在未来几十年,人民币不会贬值,汇率稳定;(2)假设深圳市未来经济增长符合正常规律,未出现类似于金融危机等的重大经济变革;(3)假设深圳地质沉降速率保持历史沉降速度,平均每年1mm 左右;(4)假设重大洪灾发生的概率和政府抗洪投入维持相对稳定;(5)假设深圳市相关防洪政策无重大调整;(6)假设未来无造成经济损失大幅度波动的极端灾难性气候。

五:模型的建立与求解5.1 分析报告中存在的问题洪灾损失体系如下图一所示:图一 洪灾损失预测体系 洪灾损失预测 致灾因素 环境因素 社会因素 防洪标准 地质沉降高度 损失估计根据附件1、附件2和附件3的分析,我们总结出研究报告中的问题如下:1、防洪标准:OECD 报告中指出的防洪标准为20年一遇,但事实上深圳的防洪标准早已达到200年一遇。

2、地质沉降高度:OECD 报告中指出2050年广州地质沉降40厘米,事实上,每年沉降只有一到两毫米,到2050年大约沉降40毫米。

3、水灾损失占GDP 比例:OECE 报告中指出水灾损失占GDP 比例较大。

报告中采用了不适合珠江三角洲的洪灾损失评估方法和计算参数如对AAL/GDP 采用线性关系,以GDP 及AAL/GDP 的比值反推AAL (年均洪灾经济损失),使得造成的灾害损失被高估。

5.2 生产总值(GDP )的预测模型—logistic 模型5.2.1 数据分析搜集深圳市各年GDP 总值(见附录【1】)对得到的数据画出散点图如下: 0123456789051015x 107x=(年份-1980)/4G D P 值/万图二 GDP 散点图由散点图可以看出GDP 呈曲线上升趋势,但由于GDP 不可能无限制地增长,所以后期增长速率会逐渐下降,类似于人口增长模型,故我们采用logistic 曲线建立模型来预测2020年及2050年GDP 。