边坡变形破坏预测的混沌与分形研究.pdfVIP

第凹卷第3期 JOURNALOFHOHAI UN／VEP~ITY河 海 大 学 学 报 ．凹 No 3 1999年5月 M t999 边坡变形破坏预测的混沌与分形研究 曾开华 陆兆溱 7 (河海大学土木工程学院 南京 210098) 摘 要 运用非线性科学理论和方法，对边坡变形破坏预测方法进行了研究，得出了反映边坡系统 动态特性的关联维数D2和Kolmogomv熵 的计算公式．实例分析表明，滑坡变形失稳过程具有混 沌和分维特性，可以用分形理论来研究滑坡预测问题． 关键词 边坡；预测；非线性动力学；混沌；分形 中图号 TU457 毒矗‘ R I 以往研究边坡问题的方法，大多数都是确定论和概率论的方法．然而，边坡作为一个真实的地质系统，存 在着与外界物质和能量的交换，因而，本质上是一个开放、耗散和复杂的非线性动力学系统，它既受内部非线 性因素的作用，又受外部随机因素的影响；既有确定性，又有不确定性随机的一面，即混沌现象．因此，我们可 以用非线性科学中的有关理论，如混沌与分形理论来解决边坡预测问题、 1 边坡演变动力学方程及单变量时间序列的拓展 实际边坡监测通常是对某一时间变量进行预测，而有些变量难以直接观测．有人认为，单变量的时序提 供的信息对于边坡复杂系统是有限的；但用动力学系统方法可以发现单变量的时序蕴藏着参与动态过程的 全部变量的痕迹，是边坡演化过程的综合反映， 由于边坡本身的复杂性，迄今还不能正确写出其系统演化的动力学方程．但我们可以得到具有n个状 态变量缸随时间演变的非线性动力学方程组 J =JS(X 2m ，…， ) i=1 2一，n (1) 对上述 n维非线性动力学系统，通过消元，使其变换为一个n阶非线性微分方程 ( ) = ， ，… ， ( )) (2) 变换后，系统新的轨迹 (f)=( (t)， (t)，…， t 一 ’(‘)) (3) 描写了同样的动力学系统．代替上式，可用不连续的时序 (t)和它的(n一1)时滞位移 ( )=( (f)， (f+f)，…， (t+(n一1)r)) (4) 其中r为延滞参数．将r选作时间序列的长时间尺度，将会保证延滞坐标线性无关， 设t≈}(i=1，2，…，Ⅳ)为监测所得的位移时间序列，取一个固定的时问间隔r，将原有序列 (t)延拓成 维相空间的一个相型分布 1： (t．1)， ( 1+r)， ( 1+2r)，…， (tl+(n一1)r) 1 ： (t2)， (t2十r)，X(t2+2r)，…， (t2十( 一1)r) l i f ： ( )， ( +r)， (f +2r)，…， ( +( 一1)r)J 这里r=kAt为延滞时问，k为整数，△f为采样间隔，则相点数为m=N一(n一1)r 墨代表n维相空间矢 量，设置为其相点，则坐标为{ ( )， (t‘+r)，…， ( +(n一1)r)} 式(5)构成m个相型，相点间的连续描 述了系统在n维相空间的演化轨迹． 收稿日期：1997-12-31 第 一作者筒介：曾开华 男．障 研究生 岩 工程专业 上要从事趣蛙稳定与非饱和上研究 10 河 海