基于DEM的泥石流堆积物产量计算方法的探讨.pdfVIP

  其中： V ─ 任意两个相邻断面间体积 i Si 、Si+1 ─ 任意两条相邻断面面积 DEF、ABC ─ 泥石流发生前、后的地面线 图 2 断面法计算原理 h ─ 相邻断面间距 整个泥石流堆积物产量为： n V ∑ V i i 1, 2 , 3 ,L , n i 1 断面法计算简单，精度取决于横断面密度，其缺点是求面积工作量大，数据容易在计算时带来 较大的误差，特别是地形较复杂时，算出的面积值误差会很大。一般只适用于山地及高差变化比较 大、自然地面较复杂的狭长地带。 1.3 数字高程模型（DEM）法 数字高程模型（DEM）是数字地面模型（DTM）的一种， DEM以数字的形式按一定的结构组织在 一起，表示实际地形的高低起伏等特征。DEM 主要有规则格网（grid）与不规则格网两种表现形式， 而不规则三角网（TIN）是不规则格网中最重要的一种形式。TIN 通过由不规则分布的数据点生成的 连续三角面来逼近地形表面。从表达地形信息的角度而言，TIN 模型的优点是它能以不同层次的分 辨率来描述地形表面，即 TIN 的三角形随点集密度变化而变化，当点集密集时生成的三角形小而密， 稀疏时生成的三角形大而疏。这与实际的地形特征恰好一致，因其能够较好地反映实际地形信息而 得以广泛使用，其函数表示如下： Z i f (x i , y i ) i 1, 2 , 3,L ,n 其中： x , y ─ 地面上任意一点的平面坐标 i i Z ─ 点（ x , y ）对应的高程 i i i DEM 法计算泥石流堆积物产量的基本原理是基于二 重积分。如图 3 所示，泥石流堆积物产量计算范围相当 于 xoy 平面上的区域 D，P2、P1 分别对应于泥石流发生 前、后的地面。用均匀栅格将地面等分，设 DEM 分辨率 （栅格大小）为 S，P1对应的地面方程为： Z i f 1 (x i , y i ) i 1, 2 , 3, L ,n P2 对应的地面方程为： Z i f 2 (x i , y i ) i 1, 2, 3,L ,n 则图 3 中所示曲顶柱体体积为： 图 3 DEM 法计算原理 71   V { f ( x , y ) =− f ( x , y ) } d s i 1, 2 , 3 , L , n i ∫∫ 1 i i 2 i i S 当栅格很小时，由于 f ( x , y ) 连续，对于同一个小区域来说， f ( x , y ) 变化很小，这时小曲 i i i i 顶柱体可近似看作平顶柱体，此时体积公式变为： V ≈ [ f ( x , y ) − f