《隧道系统动态稳定性的分析》

1、隧道系统动态稳定性的分析                隧道系统动态稳定性的分析张勇 华安增（中国矿业大学建工学院，江苏徐州221008）摘要：针对隧道系统演化发展的特点，采用信息扩散的方法对其动态稳定性进行分析，进而给出了判断其稳定状况的相关函数和指标。关键词：信息扩散；动态稳定；隧道岩土工程在施工过程中必然受到自然和人为等不确定性因素的影响，对于隧道工程则表现为一个开放性的系统，使得系统的稳定性的分析成为一个复杂的系统工程，并且是一

2、个由各种介质组成的高度复杂的动态系统，该系统的动态行为是各种因素共同作用的结果，其动态行为机制很难用公式进行准确描述，稳定性也表现出一个随时间演化发展的动态过程。现场技术人员及专家对隧道工程施工过程和运行期间的稳定性都有大致的评价，并有定量的认识。这类认识未用系统的理论进行深入的分析。利用系统动态稳定性理论，可以进一步探讨隧道系统的动态稳定性，分析隧道系统的动态行为，从而可进行信息化施工，以便更好地确保隧道支护结构体系在施工和运行过程中稳定性，更好地控制隧道开挖对周围岩体的影响，最小限度地减少工程事故的发生，确保工程的安全、高效。本文将根据信息扩散理论把这类经验与利用系统演化发展状态空间切向量

3、的特征指数，即lyapounov指数分析系统的动态稳定性，并与其动态稳定度相结合，对隧道系统的动态稳定性进行分析和研究。1 隧道系统动态稳定性的分类对于隧道系统单个状态变量实测数据的时间序列y1，y2，yn，若在时刻tj的实测值为yij，经过tj后（即在时刻tj+1时）为yi,j+1，根据文献1则有：yi,j+1=yijeleijtj(1)式中leij对应时刻tj的实测值为yij时的lyapounov指数。该状态变量在时刻tj的变化率ij为：ij=yi,j+1-yiyij=eleijtj-1(2)由此可得：leij=ln(ij+1)tj(3)上式即为lyapounov指数与状态变量变化率ij之

4、间的关系。可见，现场技术人员或专家经常依据状态变量的变化率ij评价隧道的稳定性，显然有理论根据。为了定量评价系统的动态稳定度，建议将开挖过程中隧道的动态稳定性按最大lyapounov指数lejmax的大小分为五类，即类，很稳定；类，稳定；类，较稳定；类，较不稳定；类，不稳定。对lejmax有：lej,max=maxle1j，le2j，,lenj为了由式（3）确定lej,max，必须先确定ij。对于ij的大小的取定及其与各类动态稳定性之间的界限，依据技术规范、诸多工程实践及有关专家学者的工程总结，给出如表1所示的关系。表1隧道系统动态稳定性分类类别类类类类类ij00102020808139139

5、lej,max000182 30182 30470470871 30871 3注：ij是单个状态变量的变化速率，lejmax是tj时刻的最大lyapounov指数。由表1可见，若某隧道系统的状态变量在某时刻的最大lyapounov指数为028，则隧道系统处于较稳定的状态；而在另外一时刻的最大lyapounov指数为014，则按表1既可判定隧道系统处于稳定状态，也可判定隧道系统处于较稳定状态。这类矛盾是由对隧道动态稳定性分类设置明确的区间界限引起的，如将分类模糊化，即对动态稳定性分类所依据的信息进行分配与扩散，这类现象就会消除。基于这个想法，下面先说明有关信息分配与扩散的部分理论，然后提出依据信

6、息扩散理论的隧道动态稳定性的分类。2 基于信息扩散理论的隧道系统动态稳定性的分类2.1 信息扩散原理设a=a1,a2,am为未知样本，x为一连续基础论域，aj(j=1,2,m）的观测值为xj。若a是非完备的，则存在函数（xj，x），使点xj获得的量值为1（即隶属度为最大）的信息可按（xj，x）的量值扩散到点x；，且扩散所得的信息分布q（x)=ni=1(xj，x）能更好地反映a的总体规律。由上述信息扩散原理，可构造如下信息扩散公式：q(xj,x)=0398 9exp-(x-xi)22h2(4)式中q(xj,x)观测值xj扩散到点x上的信息量；x信息吸收点，相当于信息分配中的信息控制点si；h窗宽

7、，表示信息扩散的控制范围，其值与样本a的容量有关。设样本a=a1,a2,am的最大、最小观测值为xmax、xmin，则窗宽h可按下式取值：h=1698 7(xmax-xmin)m，2m51445 6(xmax-xmin)m，6m71423 0（xmax-xmin)m,8m91420 8(xmax-xmin)m,m10(5)上述关于信息分配及信息扩散的理论，可用于对最大lyapounov指数进行信息处理，并由此对隧道系统的动态稳定性建立分类方法。2.2 隧道系统动态稳定性的分类lyapounov指数谱可视为与系统演化模式相应的样本观测值，并可据以对隧道系统的动态稳定性建立分类法。首先需要确定窗宽

8、，根据样本观测值a：le1，max、le2,max、lem,max的最大、最小观测值lemax,lemin，分别为：lemax=maxle1,max,le2，max, ,lem,maxlemin=minle1,max,le2，max, ,lem,max再由容量m,按式（5）可求得信息扩散的范围，即窗宽h其次计算扩散信息，按下式分别计算由最大lyapounov指数lej,max扩散给各动态稳定性信息控制点si的信息q(lej,max,si)：q(lej,max,si)=0398 9exp-（lej,max-si)22h2（6）式中si为信息控制点。按表1的定义，隧道系统动态稳定性分类的界限值分别

9、为：v2=0，v3=0182 3，v4=047和v5=0871 3。为使分类区域具有均一性，取v1=-0212 4，v6=1106 3，并令：si=vi+vi+12,i=1,2,5(7)最后对信息q(lej,max,si)进行归一化处理，以上各点信息属离散论域，样点地位均不相同，需进行归一化处理。依照信息总量和为1的原则，可得归一化分布的计算式为：q(lej,max,si)=q(lej,max,si)5i=1q(lej,max,si)（8）由此，可得基于信息扩散理论的隧道系统动态稳定性的分类矩阵，信息分布矩阵的具体形式可表述如下：q(le1,max,s1)q(le1,max,s2)q(le1,

10、max,s5)q(le2,max,s1)q(le2,max,s2)q(le2,max,s5)q(le5,max,s1)q(le5,max,s2)q(le5,max,s5)(9)式（9）的矩阵有m行（m为时序资料长度），各行均为某时刻隧道系统动态稳定性分类的信息，每一行有5列，分别表示该时刻隧道系统动态稳定性归于类类的隶属度。按隶属度最大原则，可由其中数值最大的一列确定该时刻隧道系统动态稳定性的类别，由于这种分类既可考虑技术人员的现场经验，又可考虑专家经验的模糊性，因而具有良好的实用性。3 工程实例根据对京珠高速公路粤境北段的部分隧道工程的监控量测的数据所进行的计算分析，得到了比较满意的结论。以

11、其中的乌坑坝隧道为例，该隧道位于乳源县乳城镇北西约4 km处，为上下行左右分线的双洞单线行车隧道，在此进行了周边收敛、拱顶下沉、围岩内部位移、衬砌应力、围岩压力及层间支护压力和型钢支撑应力量测等，通过对测试数据的理论分析和计算，得到隧道系统的最大lyapounov指数谱（图1）。表2基于信息扩散原理的隧道系统动态稳定性分类情况时间动态稳定性分类类稳定性隶属度类稳定性隶属度类稳定性隶属度类稳定性隶属度类稳定性隶属度10000 00097 20902 80000 00000 020000 00981 10018 90000 00000 030000 00000 00859 90140 00000

12、040304 30000 00767 60232 40000 050198 40000 00526 40473 60000 060000 00000 00511 30488 70000 070000 00000 00957 50042 50000 080000 00000 00448 30551 70000 090000 00000 00363 30636 70000 0100000 00000 00176 30823 70000 0110000 00000 00296 20703 80000 0120000 00000 00655 30344 70000 0130000 00000 0056

13、8 80431 20000 0140000 00000 00336 30663 70000 0150000 00000 00368 80631 20000 0160000 00000 00334 50665 50000 0170000 00000 00230 30769 70000 0180000 00000 00446 60553 40000 0190000 00000 00519 70480 30000 0200000 00000 00502 60497 40000 0210000 00000 00667 50332 50000 0220000 00000 00778 10221 9000

14、0 0230000 00000 00806 20193 80000 0240000 00000 00972 60027 40000 0250000 00000 00967 00033 00000 0260000 00187 50812 50000 00000 0270000 00631 30368 70000 00000 0280000 00600 20399 80000 00000 0290000 00476 40523 60000 00000 0300000 00085 30914 70000 00000 0310000 00105 30894 70000 00000 0320195 80

15、804 20000 00000 00000 0330100 10899 90000 00000 00000 0由图1所示的最大lyapounov指数谱，根据基于信息扩散理论的隧道系统动态稳定性分析分类进行分类。分类的具体结果如表2所示。为了详细而又客观地了解分类情况，将分类的信息分布矩阵也列入了表2中。在表2中，将动态稳定度隶属度最大的值用黑体强调。从该表中也可以得到与图1所示的相似结论：开挖初期比较小，而后开始增大，表示开挖初期，隧道系统的动态稳定度较大，动态稳定性较好，其分类是类即一般稳定类；随着围岩状况的恶化，隧道系统动态稳定度减小，动态稳定性变差，在8 d18 d这段时间内动态稳定性类

16、属于类，即不稳定类；但随着初期支护的施工，在第18 d后，隧道系统的动态稳定度有所增大，动态稳定性随之改善，分类是类即一般类；而第26 d以后，隧道系统动态稳定度开始增大，其动态稳定性在加强，动态稳定性分类可归属类。图1最大lyapounov指数lej,max示意归属同一类的隧道系统动态稳定性的隶属度也是各不相同的，如第6 d、7 d的围岩类别划分属于类的隶属度分别为0511 3、0957 5，而它们归属类的隶属度分别为0488 7、0042 5，说明第7 d的动态稳定性归于类的把握比第6 d的大，这说明每一时刻隧道系统的稳定状况是不同的。这与其实际开挖工况是一致的。由此说明了隧道系统动态稳定性是一个随时间发展的非确定量。4 结论鉴于隧道工程的复杂性及使用要求的差异，加上工程岩体的非均质性等特点，隧道工程不同施工时间的稳定性是动态变化的，而影响隧道工程系统稳定性的状态变量往往不是单一的，且这些状态变量在隧道稳定性的动态