江湖顶托模型中分岔现象研究.pdfVIP

江湖顶托模型中分岔现象的研究 吴大曼唐云 (清华大学数学科学系，北京，100084) ■ 薹 本文在以前工作的基础上将原有江湖顶托模型非线性化，同时将顶托系数设为参 数．用Newton法求出非线性模型的正平衡点．然后根据H唧itz定理判别该平衡点的稳定性 与参数的关系，最后证明系统出现的是Hopf分岔．并初步分析了此现象的工程意义． 关■谰 非线性，顶托，Ne咖法，Hu唧itz定理，H叩f分岔 0前言 荆江与洞庭湖由于其中水、沙因素与边界条件互相依赖、此响彼应，形成一系列反 馈环节，从而构成一个复杂的自我适应系统．在此系统中，城陵矶顶托是影响江湖关系变 迁的两个关键因素之一．所谓顶托是指两条河流交汇时，其中一条河流对另一条入汇的阻 碍作用．洞庭湖接纳上荆江在松滋、太平、藕池三口处的分流，并在城陵矶处与下荆江交 汇，交汇点在分流口门和荆江三次裁弯点下游．监利和螺山分别是交汇点上游和下游的第 一个水文站．江湖交汇点与城陵矶一螺山河段均在荆江监利河段内，江湖交汇河势如图1 湖汇流区域内，荆江与湖口河段里垂直相交，这使得荆江与洞庭湖之间的项托作用达到最 大，不利于洞庭湖的出流． 文献【1堋统计方法对捆江裁 弯前后以城陵矾站为代表的洞庭 湖水位抬高问题作了定量分析， 以监利和城陵矶站的水位一流量 为参量一变量，用幂函数描述江 湖相互项托的规律．文献【2]运用 控制理论中动态反馈的观点及水 流、电流之间的类比关系，建立 了一个反映江湖顶托关系的线性 模型(如式(O．1))，并对模型 中系数矩阵AC作了递推最小 二乘估计： 图1． 江湖交汇与城一螺河段 H=A日+CQ． (O．1) 式中：日=(．111，^：，心)7，I(f=1，2，3)分别为监利、城陵矶、螺山三站的水位值； 百=(丘，‘：，忘)7，五．为^，对时间f的导数：A为常数矩阵，这里％o(f=l，2，3)， 223 n12’口2I0；c为常数对角矩阵，C=d妞g(cl，。2，c3)，q0(f=1，2，3) Q=(Q．，Q2，岛)。，Qj(f=1，互3)分别为监利、城陵矶、螺山三站的来流量． 但是，根据资料分析及长期数据观测II．3“，荆江一洞庭湖区各站水位、流量时常有 失稳现象，如水位、流量周期性波动．洪水暴发等，而这正是线性模型所无法刻画的．为 进一步描述复杂的实际水利背景，本文在文献【21的基础上，将原有模型非线性化．考虑到 相互项托的重要影响，我们引入顶托系数^和九．且设之为一组参数．在求出非线性模 型的正平衡水位后，利用HuⅢitz方法判定相应于不同参数的正平衡水位的稳定性，同时 对于临界状态作定性分析；最后可以看出当取^和厶菜些值时，系统会出现H0pf分 岔，即系统存在周期解．这与水利工程中观测和分析到的各测站水位有时出现周期性波动 的现象是一致的．… 1模型的改进 由文献【1冲的数据分析得出监利、城陵矶、螺山三站的流量与水位分别有如下近似 关系： Q=d．， l=1，23 式中：正O，％l(f=l，2，3)均为常数． 另外，在文献【2】的建模过程中我们引入顼托系数^和L，其中^的太小反映了城 陵矶一螺山河段出流吼对监利一螺山河段出流口。的顶托作用的程度，即当玑固定时， ^越太，则吼越小；反之，赴则反映旱l对吼的顶托情况．为了突出考虑顶托因素对平衡 水位的影响，我们将丑，丑设为一组参数，同时，将原模型(o．1)中的常数矩阵A作如下 处理： q】 ^q2 ^=A(^，t)=屯d2l 口22 ∽z， 日3】 口32 小Ⅷ㈣． 将式(1．1)和(1．2)代A式(0．i)得 矗=A(^，^2)日+口H“ (1．3