浮游生态系统分岔特性研究.docVIP

浮游生态系统分岔特性研究 葛龙 （天津大学） [摘要] 建立了新的藻类-浮游动物生态非线性动力学模型，该模型考虑了藻类与浮游动物之间的增长关系。研究了模型平衡点的稳定性,以及Hopf分岔现象。通过对选取的分岔参数进行分析，揭示了该系统的动力学特性。 [关键词] 浮游生态系统 分岔 Research on the Bifurcation Characters of the Plankton Ecosystem Ge Long Wang Hongli (Tianjin University) [Abstract] A simple algae-phytoplankton nonlinear dynamical model was built. We studied the model to analyze the relationship of the bloom recurrence. With one parameter changing continually, a hopf bifurcation was observed. To analyze the parameters, we can find the dynamic properties of the algae-phytoplankton nonlinear dynamical model [key word] Plankton ecosystem Bifurcation 自从上世纪90年代以来，海洋生态方面的研究日趋活跃，海洋生态系统动力学模型的研究成为本领域内的一个重要方向。本文通过查阅相关学术资料，建立海洋生态系统藻类-浮游动物的动力学模型，并分析了此模型的分岔特性，揭示了系统的复杂动力学行为。 模型 考虑藻类与浮游动物的捕食关系，建立二者的模型。 对于藻类密度变化率，考虑三方面因素的影响：首先是自身的增长，增长函数符合logistic模型，其中最大生长速率r和允许载荷量K；其次是浮游动物捕食的影响：藻类的损失率与浮游动物的捕食率和存量有关，它们之间遵循第三类功能反应（III）函数，浮游动物半饱和功能反应值hA和最大单位放牧量g；最后是水流因素的影响。得到方程如下式所示： (1) 其中，A(t)代表藻类浓度随时间的变化。 对于浮游动物密度变化率，考虑三方面因素的影响：首先是捕食藻类；其次是考虑浮游动物自身的死亡；最后考虑捕捞因素的影响。得到方程如下式: (2) 其中，Z(t)代表浮游动物密度随时间的变化；f表示捕捞强度。 综上所述，得到了藻类-浮游动物模型，该模型由两个等式组成。 (3) 其它各参数意义及其取值范围如下表1： 表1 参数意义及其取值范围 符号 意义 默认取值 单位 r 藻类的最大生长率 0.5 day-1 K 藻类的允许载荷量mg/L g 浮游动物的捕食率 0.4 day-1 hA 浮游动物半饱和功能反应 0.6 mg/L i 藻类扩散率 0.01 — e 浮游动物摄食转化率 0.6 — m 浮游动物死亡率 0.15 day-1 hZ 鱼的半饱和功能反应 0.5 mg/L 动力学分析 由于捕捞强度f作为人类对海洋生态系统的一种直接干扰，对海洋生态系统的动力学特性产生直接的影响，因此本文着重探讨捕捞强度对所建模型动力学特性的影响。 在藻类扩散率i=0.01时，分别研究捕捞强度f=0.007、0.1、0.3时，系统的动力学特性。 当捕捞强度f=0.007时，系统的相图如图1所示。 此时，系统产生极限环，分别取初始值为（0.1，0.1）和（1.1，2.1），系统地运动轨迹都趋向于同一个极限环。系统存在一个稳定的极限环。 图1 当f=0. 007，系统的相图 当捕捞强度f =0.1时，系统的相图如图2所示。 此时，系统分别取初始值为（0.1，0.1）和（1.1，2.1），系统的运动轨迹有极限环和稳态两种可能。 图2 f=0. 1，系统的相图 当捕捞强度f=0.3时，系统的相图如图3所示。此时在不同的条件下，系统均趋于稳定态。 图 3 f=0. 3，系统的相图 如上所示，在藻类扩散率i=0.01时，捕捞强度f=0.007、0.1、0.3三个离散值时，系统的动力学特性如图24、25、26所示。 保持藻类扩散率i=0.01不变，当捕捞强度f的值连续变化时，系统的动力学特性如图4所示。随着参数的变化，系统先后发生了三次分岔。首先，在（1.024859，2.159144，0.006245）处发生霍普分岔；其次，在（5.116208，3.319798，0.287300）处发生Fold分岔；接着，在（5.631772，3.244563，0.284144）处产生Neutral S