水资源系统分析-第7章系统模拟研讨.ppt

第7章 系统模拟方法与应用 -7.1.1 功能和用途- 1、通过对仿真运行过程的观察和统计，得到仿真系统的输出参数和基本特性，以此来估计和推断实际系统的真实参数和性能。 2、分析系统在不同的输入或参数下，系统的行为变化。 3、通过模拟，可以获得近似的优化方案。 可用于：改造旧系统，设计新系统、发现 异常、寻找优化的运行策略。 仿真模拟的步骤 1、问题阐述 2、设置目标和项目计划 3、建立模型 4、收集数据 5、编制程序 6、验证 7、设计试验方案 8、模型运行 9、反馈与修正 10、实施 -水资源系统模型- 描述复杂的水资源系统 1、1955，Morrice和Allan模拟尼罗河规划问题，考虑了17个水库（电站）。 2、美国陆军工程师兵团，哥伦比亚河流域模拟，分析了25个水库，45条输水线路。 3、1974，MIT开发了多目标流域规划模拟模型。 系统模型 b.系统模型将所研究的流域视作一种动力系统，建立输入与输出资料的数学模型，这样就可以由新的输入预测输出。这种模型不牵涉物理实质，只要求系统响应。 物理模型-SWAT 简介 SWAT（soil and water assessment tool）是由美国农业部农业研究中心（USDA-ARS）开发的流域尺度模型，可以用于模拟地表水和地下水的水质、水量，预测土地管理措施对大尺度流域水文、泥沙和农业化学物质的影响。SWAT模型自20世纪90年代初开发以来，已经经历了不断的发展，先后经过了SWAT94.2、SWAT96.2、WAT98.1、SWAT99.2、SWAT2000、SWAT2005系列官方版本。 SWAT模型是一个具有很强物理机制的、长时段的流域水文模型，适用于具有不同的土壤类型、不同的土地利用方式和管理水平的复杂流域，并能在资料缺乏的地区建模。从模型结构上看（图1），SWAT属于第二类分布式水文模型，即在每一个子流域上应用传统概念性模型，再进行汇流演算，最后求得出口断面流量. 运行SWAT 模型时，首先要求输入流域的数字高程模型（DEM） 模型基于一定的原理把整个流域划分成不同的子流域或水文响应单元（Hydrologic response units，HRUs） 再分别输入流域的土地覆被栅格图和土壤栅格图 通过计算逐个水文响应单元的径流量，最后得到流域的总径流量。 在子流域内部生成水文响应单元（HRU）。 SWAT模型根据子流域内不同的土壤类型、土地利用类型组合来生成HRU。 在SWAT模型中有两种划分HRU的方法：第一种是把每个子流域视为一个HRU，其性质由子流域内的主要土地利用分类和土壤类型来确定， -蒙特卡罗(Monte Carlo)法- 用计算机来模拟统计的随机抽样过程。 初始的蒙特卡罗方法只用于生成某种随机数列，现在已经有了更多的用途. 只要能把问题转化为随机抽样，即可应用. 历史 早在17世纪，人们就知道用事件发生的“频率”来决定事件的“概率”。 19世纪人们用投针试验的方法来决定圆周率π。本世纪40年代电子计算机的出现，特别是近年来高速电子计算机的出现，使得用数学方法在计算机上大量、快速地模拟这样的试验成为可能。 比如:y=x2从0积到1。结果就是下图红色部分的面积： 用matlab模拟，做一百万次(即共取1000000个点)，结果为0.3328。 理论基础 理论基础是概率论中的大数定理。即在相同的条件下，对事件A进行n次独立的试验，当n无限增大时，事件A的n个观测值的平均值依机率收敛于其数学期望值 蒙特卡罗方法的近似值与真值的误差问题，概率论的中心极限定理对此作出了刻画。 起源： 第二次世界大战期间，美国论证出了制造原子弹的可能性，但要制出实际可用的核武器，逐项解决大量复杂的理论和技术问题，如中子轨道和辐射轨道等等问题。描述这些过程需要相当复杂的微分、积分的耦合方程组。 科学家们采用建立基础的物理模型，用随机抽样法在计算机上进行模拟的方法。这种方法与数学模拟法不同，不从物理模型建立数学模型（太复杂），而是在计算机上进行试验，用随机抽样法解决问题。 Monte Carlo是一个赌城，赌博的特点就是随机性，由此对本法进行命名。 Monte Carlo方法主要创始人： Stanislaw Marcin Ulam Enrico Fermi John von Neumann Nicholas Metropolis 1949年Metropolis和Ulam发表了第一篇论文，学术界以此作为Monte Carlo法诞生的标志。Monte Carlo法是在计算机发展起来后才出现的。根据基础的物理模型在计算机上改变条件、观测结果，作出结论。因此也称为“计算机实验”。 最优化问题的随机枚举法 -