水葫芦生及治理模型.pptVIP

水葫芦生长及治理模型 周继顺 0018103 复旦大学 03.04.02 提要 提出问题 分析问题 提出解决问题的方案 利用数学方法解决问题 应用举例 模型改进 1.水葫芦污染的现状 水葫芦是50多年前从南美的巴西引进来的，是外来物种。它喜欢生长于温暖向阳及富含养分的水域中，无性繁殖能力特别强，每年的九、十月份是生长旺季。在适宜的条件下，每5天就能繁殖新株，也能开花结实产生种子而进行有性繁殖，一枝花大约结300粒种子。一公顷水面的水葫芦就能挤满200万株，重达300多吨。在上海，每年内河水域的水葫芦总量高达200万吨，覆盖了405.57平方公里中25%的水面。 2.水葫芦污染的形成 由于水葫芦在水中不及时打捞上岸，会腐烂变臭，从而造成环境污染。上海的水葫芦污染已日趋严重，水域整治已迫在眉睫，在上海水域的水葫芦污染形成有以下几个原因。 由长江上游及苏浙地区漂移而下。 水葫芦自身强大的繁殖能力。 水葫芦与周围的许多昆虫和病原菌建立了稳定的生态平衡关系。 3.水葫芦的危害 河道堵塞 航运受阻 水质污染 破坏生态 4.建立水葫芦生长方程 由于上游的水葫芦也在不断繁殖，所以每天从上游漂移进上海的水葫芦的速度也不断增加，于是我们设流速为f(t)，而已在上海水域的水葫芦自身繁殖加速度与水域中水葫芦增加速度成正比，我们设为k v。而这个k的大小取决于此水域生态系统的大小。我们设水葫芦增加速度为v，这样我们可排出如下方程。 5.计算水域中水葫芦的数量 由以上方程我们可以利用微分方程数值方法计算每一个t时刻，水葫芦增加量v的值，我们只需得到若干个v的值，就能利用数值积分方法计算水葫芦数量M。（见右边） 6.水葫芦数量增长分析 我们由假设可知流速函数f ( t )必为递增的，且其一阶导数也为递增的。而k必为正值，那么可知v也为递增的。 这样，水葫芦的增长必定越来越快，所以我们必须及时对它进行治理，否则其必将迅速蔓延，无法控制。以下我们将讨论如何治理。 7.水葫芦治理方法 截断上游漂移的水葫芦 大力打捞漂浮在水域中的水葫芦 运用生物（如象甲）防治以降低水葫芦繁殖速度 破坏水葫芦生态系统 8.各种方法的单独应用（1） 我们先使用打捞这一最简单的方法。设v(0)=10000吨/月，而流速f( t )为 吨/月,且k=0.5，每月可打捞4000 吨,那么我们可以得出方程如下： 8（2） 使用RK4方法可观察到前十月的速度曲线从中可知单纯打捞是无济于事的，水葫芦的增长速度很快会超过打捞速度。 8（3） 再看利用生物防治方法，设象甲幼虫定期投放入水域，而随着幼虫的生长，其对水葫芦的食用量也越来越大，我们设其食用量为 ，则方程可写为： 8（4） 使用RK4方法可观察到前十月的速度曲线从中可知单纯生物治理也是没有什么作用的。 8（5） 由于破坏水葫芦生态系统只是单纯减少之前的系数，速度仍增长很快，所以也是没有很大作用的。 从观察可知截断上游漂移的水葫芦可使方程中的 这一项没有，这样再使用以上方法就会有明显效果。 9.各种方法的综合使用（1） 我们在截断上游污染的前提下，使用8（1），8（3），8（5）三种方法。得到方程为： 9（2） 我们也观察10个月的速度变化（取值如前），可以看出速度得到了控制，污染不会再蔓延了。 10.具体例子分析（1） 几点假设： 黄浦江中某水域现有水葫芦10万吨 水葫芦繁殖系数0.5,生态系统破坏后为0.3 现今水葫芦增长速度为1万吨/月 每月可打捞4000吨 每批象甲虫食用水葫芦速度为500 吨/月 10（2） 在以上几点假设下，问二年内黄浦江中此水域能否完全消除水葫芦污染？ 我们先建立速度方程。 10（3） 我们通过 matlab命令ode45作出v在两年内的函数图象。（见右边） 10（4） 从图中可以看出当t=13时，水葫芦的增长速度几乎为0，即在第13个月时，我们就能控制水葫芦增长了。 接下来我们就要计算能不能在余下的几个月中清理完已在水域中的水葫芦。 10（5） 要计算水域中的水葫芦数量，我们可以利用已学过的数值积分的Matlab命令。（具体操作见右边） 10（6）