CN101403798B - 一种滑坡主滑方向的识别方法 - Google Patents
一种滑坡主滑方向的识别方法 Download PDFInfo
- Publication number
- CN101403798B CN101403798B CN2008102330373A CN200810233037A CN101403798B CN 101403798 B CN101403798 B CN 101403798B CN 2008102330373 A CN2008102330373 A CN 2008102330373A CN 200810233037 A CN200810233037 A CN 200810233037A CN 101403798 B CN101403798 B CN 101403798B
- Authority
- CN
- China
- Prior art keywords
- grid
- theme
- gliding mass
- data
- sliding
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Abstract
一种滑坡主滑方向的识别方法,其特征在于按下列步骤进行:获取滑体的工程勘察资料,构建滑体的基础数据地图主题,对滑体进行栅格单元划分;并获取地质剖面栅格主题、钻孔资料栅格主题和地表高程栅格主题;通过GIS系统得出滑体厚度栅格主题;结合滑体岩土重度指标γ,得出每个栅格(i,j)滑体的重力W(i,j):调用坡度命令,得到每个栅格(i,j)的下滑力T(i,j):调用坡向命令,建立滑面的倾向地图主题,得到每个栅格(i,j)滑体的滑动方位角Asp(i,j);获得滑体的主滑方向。本发明的显著效果是:受外部因素的影响小,可靠性高,能为滑坡治理工程提供准确、科学的依据,操作简单,保证了优化设计,节约了投资成本。
Description
技术领域
本发明涉及山体滑坡预防和施工技术领域,具体地讲,是一种滑坡滑体主滑方向的识别方法。
背景技术
在滑坡治理工程设计中,支挡结构(如抗滑桩、重力式档墙等)的设计抗力应该与滑体下滑推力互为反作用力,而整个滑体下滑力的方向即为滑坡的主滑方向。滑坡主滑方向的确定是计算滑动推力、设计滑坡支挡抗力,并进一步设计支挡结构的重要参数之一。滑坡主滑方向判别不准确,可能导致滑动推力大小计算不正确,以及设计滑坡支挡抗力的方向错误,影响到支挡结构的布置及设计偏差,或者盲目增大治理工程投资,甚至会导致整个支挡工程的失效。工程实践中,基于滑坡主滑方向判别不准确而导致治理工程失败的案例并不少见,所以滑坡主滑方向的准确判别意义十分重大。
目前,对于滑坡主滑方向的确定主要是根据现场调研、勘察进行初步判断,或者根据工程钻孔揭露滑面上岩土擦痕的方向来确定滑坡滑体的主滑方向。现场踏勘观察的判断方法太宏观而又易受人为主观影响,钻孔揭露滑面上岩土擦痕的方法又太局部或者细观,难能代表滑体的整体主滑方向。
现有技术的主要问题:受现场实际条件、局部破坏情形及人为主观因素等多方面的影响较大,判别值存在较大的预测误差,可靠性低,且操作复杂。
发明内容
本发明的目的是提供一种滑坡主滑方向的识别方法,能克服外部因素的影响、预测的准确性好、可靠性高、且操作简单。
为达到上述目的,本发明所述的一种滑坡主滑方向的识别方法,其关键在于按下列步骤进行:
步骤一、获取滑体的工程勘察资料:
对滑坡滑体进行工程勘察,获取该滑坡滑体的工程勘察资料,所述工程勘察资料分为:含所述滑体地形等高线的滑坡平面图、勘察剖面图、所有钻孔柱状图、滑体岩土重度指标γ,其中所述勘察剖面图与钻孔柱状图中包含滑面的位置资料;
步骤二、将所述工程勘察资料发送给地理信息系统GIS进行空间分析,由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)为滑体所处滑面与水平面的夹角。
步骤八、得出所述每个栅格(i,j)滑体的下滑力T(i,j):
T(i,j)=W(i,j)×sinθ(i,j);
步骤九、调用GIS系统空间分析功能中的生成坡向命令,建立所述滑面的倾向地图主题,得到每个栅格(i,j)滑体的滑动方位角Asp(i,j);
步骤十、获得滑坡滑体的主滑方向:
所述每个栅格(i,j)滑体的滑动方位角Asp(i,j),也是该栅格(i,j)滑体的倾向方位角。
所述每个栅格(i,j)滑体的倾向方位角:是该栅格(i,j)滑体的下滑力T(i,j)的水平分量与正北方向之间的夹角。
所述每个栅格(i,j)滑体的坡角0≤θ(i,j)≤90°,所述每个栅格(i,j)滑体的下滑力T(i,j)的滑动方位角0≤Asp(i,j)<360°,所述整体的下滑力矢量对应的方位角0≤Asp<360°。
本发明的显著效果是:受外部因素的影响小,主滑方向预测值的可靠性高,能为滑坡治理工程提供准确、科学的依据,且操作简单,保证了优化设计,节约了投资成本。
附图说明
附图1为本发明的流程步骤图;
附图2为本发明滑体重力W(i,j)的分解图;
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明:
如图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);
如图2:步骤八、得出所述每个栅格(i,j)滑体的下滑力T(i,j):
T(i,j)=W(i,j)×sinθ(i,j);
步骤九、调用GIS系统空间分析功能中的生成坡向命令,建立所述滑面的倾向地图主题,得到每个栅格(i,j)滑体的滑动方位角Asp(i,j);
步骤十、获得滑坡滑体的主滑方向:
如图3所示:所述每个栅格(i,j)滑体的滑动方位角Asp(i,j),也是该栅格(i,j)滑体的倾向方位角。
所述每个栅格(i,j)滑体的倾向方位角:是该栅格(i,j)滑体的下滑力T(i,j)的水平分量与正北方向之间的夹角。
Claims (1)
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);
步骤十、获得滑坡滑体的主滑方向:
所述每个栅格(i,j)滑体的滑动方位角Asp(i,j),也是该栅格(i,j)滑体的倾向方位角;
所述每个栅格(i,j)滑体的倾向方位角:是该栅格(i,j)滑体的下滑力T(i,i)的水平分量与正北方向之间的夹角;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102330373A CN101403798B (zh) | 2008-11-12 | 2008-11-12 | 一种滑坡主滑方向的识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008102330373A CN101403798B (zh) | 2008-11-12 | 2008-11-12 | 一种滑坡主滑方向的识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101403798A CN101403798A (zh) | 2009-04-08 |
CN101403798B true CN101403798B (zh) | 2011-01-26 |
Family
ID=40537869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008102330373A Expired - Fee Related CN101403798B (zh) | 2008-11-12 | 2008-11-12 | 一种滑坡主滑方向的识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101403798B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102568011B (zh) * | 2011-09-02 | 2013-08-14 | 广西远长公路桥梁工程有限公司 | 根据等高线图自动生成数字航道图的方法 |
CN103455711B (zh) * | 2013-08-15 | 2017-07-11 | 广州地理研究所 | 一种基于机理分析的小流域地区滑坡灾害风险区划方法 |
CN105184050A (zh) * | 2015-08-10 | 2015-12-23 | 南京师范大学 | 一种gis与人工智能技术支持下的滑坡危险度评价方法 |
CN105572736B (zh) * | 2016-01-07 | 2017-12-05 | 西南石油大学 | 一种利用地震资料定量预测坡折带的方法 |
CN106777734A (zh) * | 2016-12-27 | 2017-05-31 | 湖南师范大学 | 一种山地滑坡预测方法 |
CN110781825B (zh) * | 2019-10-25 | 2023-05-23 | 云南电网有限责任公司电力科学研究院 | 一种电网滑坡区域识别系统及方法 |
CN113932791B (zh) * | 2021-09-15 | 2023-11-07 | 江苏徐工工程机械研究院有限公司 | 一种露天矿区装卸作业区地图采集方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1584542A (zh) * | 2004-05-28 | 2005-02-23 | 中国科学院力学研究所 | 一种水诱发滑坡模拟试验装置及坡面位移监测方法 |
CN101009024A (zh) * | 2007-02-09 | 2007-08-01 | 肖盛燮 | 一种滑坡灾害可视化的实现方法 |
-
2008
- 2008-11-12 CN CN2008102330373A patent/CN101403798B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1584542A (zh) * | 2004-05-28 | 2005-02-23 | 中国科学院力学研究所 | 一种水诱发滑坡模拟试验装置及坡面位移监测方法 |
CN101009024A (zh) * | 2007-02-09 | 2007-08-01 | 肖盛燮 | 一种滑坡灾害可视化的实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101403798A (zh) | 2009-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101403798B (zh) | 一种滑坡主滑方向的识别方法 | |
VandenBerg et al. | Late Mesozoic and Cenozoic movements of the Italian peninsula: further paleomagnetic data from the Umbrian sequence | |
CN102128625B (zh) | 重力辅助惯性导航系统中重力图匹配的初始匹配方法 | |
CN109891191B (zh) | 用于随钻自主导引的机会传感器融合算法 | |
CN105509628B (zh) | 一种磁测定位装置以及利用该装置进行滑坡深部位移监测的方法 | |
CN107991714A (zh) | 基于湖盆古地形恢复的定量方法 | |
WO2018145557A1 (zh) | 一种河流历史形态数字复原方法 | |
CN109826174A (zh) | 一种边坡加固深部区域范围确定方法 | |
Henton et al. | Crustal motion and deformation monitoring of the Canadian landmass | |
CN102736118A (zh) | 一种用于全球重力场测量的综合型卫星系统 | |
Leprêtre et al. | Vertical movements along the northern border of the West African Craton: the Reguibat Shield and adjacent basins | |
CN102393535A (zh) | 基于双星能量插值原理的卫星重力反演方法 | |
CN105911592B (zh) | 一种实钻约束的三维地震资料长波长静校正方法 | |
CN102183788B (zh) | 地层条件下横波测井曲线的合成方法 | |
Engel et al. | Using speleothems to constrain late Cenozoic uplift rates in karst terranes | |
Szczygieł | Quaternary faulting in the Tatra Mountains, evidence from cave morphology and fault-slip analysis | |
Haque et al. | A magnetostratigraphic age constraint for the proximal synorogenic conglomerates of the Late Cretaceous Cordilleran foreland basin, northeast Utah, USA | |
Teixell et al. | The geology of vertical movements of the lithosphere: An overview | |
Bender et al. | Gas age–ice age differences and the chronology of the Vostok ice core, 0–100 ka | |
CN106522924A (zh) | 一种随钻测量中方位角的获取方法 | |
Estrada et al. | Pilot phase of the Aguada Pichana Este block, gas window | |
Chen et al. | Active Faulting of Landforms along the Tuosuhu‐Maoniushan Fault and Its Seismotectonic Implications in Eastern Qaidam Basin, China | |
Armstrong et al. | Declining Basal Motion Dominates the Long‐Term Slowing of Athabasca Glacier, Canada | |
Liu et al. | Intelligent filter for accurate subsurface heading estimation using multiple integrated mems sensors | |
Anzidei et al. | Lake Albano: bathymetry and level changes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110126 Termination date: 20141112 |
|
EXPY | Termination of patent right or utility model |