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