富春江电厂山体滑坡分析.pptVIP

富春江电厂山体滑坡分析 富春江水电站是一个低水头河床式电站，挡水建筑物自右至左有灌溉渠道、船闸、溢流坝段、鱼道和发电厂房等组成，全长554.35米。由于水工建筑物的受力情况复杂，受各种因素的影响，工作状况随时都在变化，为了掌握建筑物的原体变形，即监测大坝的稳定性，确保安全，对大坝轴线端点进行变形观测，从1977年起由同济大学测量系重建了工程监测平面控制网，用以校核定大坝视准线位移左、右两端点有否变动。电厂左边有一山体有滑动迹向，因而在其上设置了8个控制点，以大网的J1、J2为控制点来监测滑坡的变化。 二　数学模型 2.1　位移量的求法 　　滑坡的原始数据是各时期的点位坐标，而用于拟合的是各时期的位移差。因而需要把坐标数据转换成位移差。如果第一期的横坐标为 ，纵坐标为 , 第i期的横坐标为 ，纵坐标为 。那么第i期的总位移量可以用以下公式得　　　　　　　　　　　　　　　　　　　　　 2.2　对所求方向进行投影的处理 　　由于在X方向的位移和在Y方向的位移合成总位移后的方向是随机的，而一般只关心的是在滑坡方向的位移。那么就有必要将总的位移进行投影 滑坡的方向的方位角为a,总位移的方位角为c,那么其夹角b为 　　　　　　　　　　　　　　　　　　　 根据这个原理，可以求出总位移在任何方向的位移。 如果总位移为,在滑坡方向上的位移为，总位移方向与滑坡方向的夹角是b。那么，应用如下的公式： 　　　　　　　　　　 　　　　　　　　　　　　　　　　　　 可得到总位移在滑坡方向上的位移分量。 在本程序中提供了对话框，让用户输入想求的方向的方位角，然后再与总位移的方向求夹角。这样就可以投影过去了。 2.3　直线拟合 设直线方程为y=ax+b 第i点的误差方程为 　 代入直线方程得 　 然后，把a b 当成未知数写成向量的形式 这样就得出了一个分量的误差方程，把所有的观测值都写成误差方程的形式中，然后再写成矩阵形式为 　 　 再进行平差后，就可以得到未知数a,b的值。有了直线的参数后，就可以画出直线图，算出各点到直线的噪声和各点的信号。 2.4　曲线拟合 设第一段的曲线方程为 ，那么在第一段中的第i点的误差方程为 　　 把曲线方程代入可得 把参数提出来当未知数 　　 如果把第一段中所有的观测值都写在一起组成矩阵形式，那么 　　　　　　　　 V=A*X-L 样条函数要求在连接处连续，一阶导数连续，二阶导数连续。这个可以用条件式来进行约束。设第i点为连接点，左边段样条参数为，右边为，那么 由连续条件可得： 　　　　　　　 　　　　　　　 有几个连接点就写几组这样的方程。跟上面一样，最后也可以写成矩阵形式。 　　　　　 CX-W=0 其中C是关于x的矩阵，而且X是由参数a b c d 组成的矩阵。 未知数X根据样条段数的不同长度也不同。 假设有10个点，共分三段，那么有两个连接点，那么　　　　　　　　　　 三 程序设计 3.1　数据前期处理 程序设计成单一任务，一开始运行便读入数据。在读入数据的同时会计算每期坐标数据与第一期数据之差，便是每期的位移量.再把X和Y方向的位移量进行矢量合成，合位移的值为 ，方位角的正切为y/x，再具体看y、x的正负来确定方位角的象限。最后把总的位移量投影到指定方向。这个方向是由投影方向角α决定的。 3.2　数据格式： 观测的个数 观测年　观测月分　观测日期 第一点的X　第一点的Y 第二点的X　第二点的Y ...... ...... 最后一点的X　最后一点的Y 例如 00 01 20 835.7196 6403.2005 845.1479 6375.9992 866.1237 6403.4051 863.6678 6397.9559 867.1476 6382.4196 876.6244 6408.2660 883.0623 6394.8970 879.6478 6383.4402 四　实例计算分析 自1999年到2004年，沿江各点均向江里移动了几毫米，甚至更多。因而把大网中的J1、J2 这两点作为已知点平差滑坡上的8个点时，应该考虑到