主权项 |
1.一种滑坡主滑方向的识别方法,其特征在于按下列步骤进行:步骤一、获取滑体的工程勘察资料:对滑坡滑体进行工程勘察,获取该滑坡滑体的工程勘察资料,所述工程勘察资料分为:含所述滑体地形等高线的滑坡平面图、勘察剖面图、所有钻孔柱状图、滑体岩土重度指标γ,其中所述勘察剖面图与钻孔柱状图中包含滑面的位置资料;步骤二、将所述工程勘察资料发送给地理信息系统GIS进行空间分析,由GIS系统构建出所述滑体的相关基础数据地图主题,该基础数据地图主题分为地质剖面数据主题、钻孔资料数据主题和地表高程数据主题;步骤三、对所述滑体以边长L1×L2进行栅格单元的划分,可以划分为M×N个单元,栅格单元标记为(i,j),其中:1≤i≤M,1≤j≤N,每个栅格(i,j)滑体的柱体高度h(i,j);并对所述GIS系统内的地质剖面数据主题、钻孔资料数据主题和地表高程数据主题进行相应的栅格单元划分,分别获取地质剖面栅格主题、钻孔资料栅格主题和地表高程栅格主题;步骤四、通过GIS系统,对所述地质剖面栅格主题和所述钻孔资料栅格主题进行空间分析,构建出滑面高程地图栅格主题;步骤五、通过GIS系统,对所述地表高程栅格主题和所述滑面高程地图栅格主题进行空间分析,构建出滑体厚度栅格主题;步骤六、通过GIS系统的空间分析功能,结合所述滑体岩土重度指标γ,获得所述每个栅格(i,j)滑体的重力W(i,j):W(i,j)=L1×L2×h(i,j)×γ;步骤七、调用GIS系统空间分析功能中的生成坡度命令,建立所述滑面的坡度地图主题,得到所述每个栅格(i,j)滑体的坡角θ(i,j);步骤八、得出所述每个栅格(i,j)滑体的下滑力T(i,j):T(i,j)=W(i,j)×sinθ(i,j);步骤九、调用GIS系统空间分析功能中的生成坡向命令,建立所述滑面的倾向地图主题,得到每个栅格(i,j)滑体的滑动方位角Asp(i,j);步骤十、获得滑坡滑体的主滑方向:根据所述下滑力T(i,j)和滑动方位角Asp(i,j),获得每个栅格(i,j)滑体的下滑力矢量<img file="DEST_PATH_RE-FSB00000208359500011.GIF" wi="146" he="77" />通过矢量求和运算求得滑坡滑体整体的下滑力矢量<img file="DEST_PATH_RE-FSB00000208359500012.GIF" wi="80" he="86" /><img file="DEST_PATH_RE-FSB00000208359500013.GIF" wi="317" he="97" />该整体的下滑力矢量<img file="DEST_PATH_RE-FSB00000208359500014.GIF" wi="48" he="79" />对应的方位角Asp即为所述滑坡滑体的主滑方向;所述每个栅格(i,j)滑体的滑动方位角Asp(i,j),也是该栅格(i,j)滑体的倾向方位角;所述每个栅格(i,j)滑体的倾向方位角:是该栅格(i,j)滑体的下滑力T(i,i)的水平分量与正北方向之间的夹角;所述每个栅格(i,j)滑体的坡角0≤θ(i,j)≤90°,所述每个栅格(i,j)滑体的下滑力T(i,j)的滑动方位角0≤Asp(i,j)<360°,所述整体的下滑力矢量<img file="DEST_PATH_RE-FSB00000208359500015.GIF" wi="48" he="81" />对应的方位角0≤Asp<360°。 |