物联网 我不是很熟悉，个人建议，校园位置服务就是基于位置的服务（LBS）里面的应用。你可以 结合 校园三位景观的建模，来写基于位置的服务。
校园三维景观建模 就是属于三维GIS的范畴，LBS 又结合了GPS的内容。属于学科交叉。
既然你在这里问，那么很可能你是个本科生，因为研究生导师会帮着选好的。
校园三位景观建模+LBS，写个差不多，让你论文通过已经绰绰有余了。
浅谈GIS技术在水利工程中的应用展望论文
摘要： 近年来， GIS （地理信息系统） 技术在水利工程应用领域中发挥着技术先导的作用。文章通过分析目前GIS技术在水利规划、水资源管理等水利工程行业的具体应用形式， 结合地理信息学科的发展方向， 展望GIS在水利工程中的应用前景。
关键词： GIS； 水利工程；
1 GIS技术概述
20世纪60年代， 世界上第一个GIS （加拿大地理信息系统CGIS） 诞生， 其核心是用计算机来处理和分析具有空间属性地理信息。GIS技术能够有效地管理地理信息资源， 是其能够有效利用的核心技术。20世纪以来， 计算机和网络技术日的新月异极大地推进了GIS在我国的发展进程， 由于应用后能够明显地提高工作效率和经济效益， 因此， GIS技术已成为资源与环境各领域应用中不可或缺的前沿技术。
近20年， 我国经济水平持续增长， 人民生活质量稳步提高， 同时， 水资源需求与水资源利用率相对不足的矛盾逐渐加剧， 这促使水利工作者不断探索如何利用现代信息技术手段缓解这一矛盾。为有效地利用水资源， 在当前条件下最大限度地发挥水利工程的调节作用， 减少建设、管理人员的投入量， GIS技术作为信息化的体现之一， 在水利工程各环节中的应用范围不断扩大， 应用层次也逐渐深入。通过近年来的发展， GIS的应用形式已从最初单纯的可视化应用过渡为集合分析、模拟、预测等多位一体的复杂应用， 功能也提升为对多时期的地理信息变化进行动态监测和分析比较， 将数据收集、空间分析和决策过程统一打包的应用。实践表明， 通过使用GIS系统， 有利于设计人员对工程进行总体规划， 方便施工人员对实施过程进度和质量把控， 有效提高管理人员的工作效率， 从而提升水利行业整体信息化水平。
2 GIS技术在水利工程中的应用
水利工程为消除水害和开发利用水资源而修建， 重要性不言而喻， 而其本身具有一定的地理空间属性， 对工程的地理环境依赖较高。GIS的发展恰好为水利工程的空间属性提供了可行的表现途径和有效的模拟、分析方法。运用GIS技术， 矢量化地理信息数据， 搭建水利工程地理信息平台， 能够直观地展示工程环境， 结合水利资料， 实现环境分析模拟， 有利于管理人员决策[1]。
目前， GIS在水利工程方面的应用主要表现在以下几方面：
（1） 水利工程规划
GIS技术广泛应用于水库选址、复杂工程布局、工程测量[2]、库容量计算、开挖土石方量计算、工程建设监测[3]、工程变形监测等方面， 促使水利工程规划、管理科学发展。
（2） 水资源管理
运用GIS技术能够确定水资源分区， 建立科学有效的管理模式， 分析水资源量， 直观地展示水资源分布和数量的动态变化， 有助于实现水资源的合理利用。
（3） 防洪减灾
GIS应用主要表现在防汛决策支持系统[4]和洪灾损失评估， 它以GIS为基础， 实现了决策方案、防汛信息、损害范围的直观化和形象化表达， 为全国防汛决策提供有力的技术支撑。
（4） 水土保持
利用GIS技术能够进行水土流失预测和动态监测[4]， 查询、统计、分析土地分区和土地利用情况， 有利于重点区域水土流失综合治理措施的制定， 提升治理效果。
（5） 水质监测
建立GIS水质监测应用平台， 能够掌握水质的实时动态变化， 及时辨析污染源头， 有利于地表水、地下水储量分析， 模拟水量调度， 提高水资源利用水平。
（6） 水文预报
通过GIS技术， 整合水文基础数据， 实现基础背景数据管理， 对空间和属性数据的查询， 统计数据以及显示检索， 为水文预报提供基础数据支撑[5]。
3 GIS技术应用展望
随着计算机技术的发展， GIS技术在水利工程应用的内容已从结合地理要素的水利信息查询展示发展到利用GIS空间运算完成水利信息的分析、计算、模拟与统计。近年来， 结合遥感、GPS等前沿技术， 构建3D、4D地理信息系统平台， 实现多维度水工建筑仿真应用已成为水利信息化应用新趋势 （如三维可视化洪水淹没分析与灾情评估系统、4D水利施工管理系统等[6]） 。水利要素具有的空间地理属性促进了GIS在水利应用中的高速发展， 纵观发展历史， 水利GIS应用的发展趋势主要表现为应用内容的发散性扩展及应用技术的融合性集成。
3。1 应用内容扩展
利用GIS技术， 实现工程环境分析模拟的形象化， 并且在一定程度上提高了获取信息的时效性和准确性， 有利于管理和决策。因此， 今后GIS技术将逐渐覆盖水利行业的各个方面， 促使水利信息化进一步发展。
（1） 洪水模拟
利用洪水历史资料， 结合GIS平台， 建立合理的流域模型， 通过可视化模拟， 分析流域洪水成因、洪水特性及其规律；通过地理信息要素分类， 分析各主要河道的自然条件；运用GIS技术， 进行上下游洪水演进模拟分析， 促使建立一定洪水标准下合理的蓄、滞、泄关系及管理措施等。
（2） 供水预案
建立供水、需水、蓄水地理信息管理平台， 通过电子地图可视化展示水源地分布， 有助于研究可能的供水方案与主要工程措施及用水管理与供水水质保证措施。
（3） 水质分析
实现水质采集自动化， 分析地质、水质信息的地理属性， 应用GIS技术研究土壤侵蚀分区；分析各类可能污染源；利用历史资料， 模拟预测不同规划水平年的污染负荷量和水质变化， 以此为基础研究保护水资源应采取的措施。
（4） 工程测量
应用无人机搭载遥感、GPS设备， 采集、处理工程测量数据， 利用GIS技术实现测量数据可视化， 减少外勘人力成本， 提高工程测量的.精度与效率， 便于数据成果的合理应用。
3。2 应用技术深入
（1） GIS技术发掘
随着GIS技术的发展， 多种空间数据结构 （如“真三维”、“时空四维”） [7]应运而生， 面向对象的数据模型和实用的界面语言正处于高速发展阶段， 数据自动输入技术不断完善， 各行各业GIS应用模型开发力度不断加大， GIS的网络共享能力 （webGIS） 不断增强。在这样的大环境中， 水利与GIS的结合必将更加深入， 具体表现为水利数据信息的表达将趋于多维度、丰富性、立体性、共享性， 水利应用模型将趋于结构化、可扩展性， 形成整个水利行业的GIS综合管理平台， 逐渐打开GIS综合系统与水利专业系统共存共荣的局面。
（2） 相关技术集成
鉴于水利空间数据的时间性和复杂性， 单一的GIS技术很难满足水利要素的处理要求。因此， 在水利工程的应用中， “3S”技术 （全球定位系统GPS、遥感RS、地理信息系统GIS） 、无人机技术的集成已成为必然的发展趋势。GPS为GIS的快速定位和更新提供手段， 遥感技术的多谱段、多时相、多传感器和多分辨率的特点， 为GIS不断注入“燃料”， 反过来又可利用GIS支持从遥感影像数据中自动提取语义和非语义信息， 而无人机则可作为遥感设备提供载体[7]。利用无人机航拍， GPS和RS赋予了GIS实时、动态属性。3S技术整体结合所构成的水利地理信息系统是高度自动化、实时化的GIS系统。这种系统不仅具有自动、实时地采集、处理和更新数据的功能， 而且能够分析和运用数据， 为水利各学科应用提供科学的决策咨询。
多媒体地理信息系统 （MGIS） 将文字、图形 （图像） 、声音、色彩、动画等技术融为一体， 为GIS应用开拓了新的领域和广阔的前景[8]。它能够以最直观的方式表达和感知空间地理信息， 以形象化的、可触摸 （触屏） 的甚至声控对话的人机界面操纵空间地理信息处理的技术。应用MGIS的水利地理信息系统将对结构、功能及应用模式的设计产生极大的影响， 使得各类信息的表现形式更丰富， 更灵活， 更友好。
4 结语
毋庸置疑， 水利GIS应用的基础是地理空间数据和水利专题要素数据， 相关专业 （如测量专业） 的发展将对其产生明显的限制或促进作用。因此， GIS在水利工程中的应用水平应依据各专业的发展稳步提高， 防止闭门造车， 构建海市蜃楼。
复杂技术集成的GIS在水利工程上的应用将会成为今后一段时期内水利信息化的主要趋势。搭建集工程规划、建设、管理为一体的水利工程综合地理信息平台， 能够实时、有效地展示工程运行情况， 模拟防洪工况， 提高管理决策能力， 为水利信息化事业发展添砖加瓦。
参考文献
[1]刘林宁。地理信息系统在水利工程中的应用探讨[J]。农村经济与科技， 2016， 27 （18） ：65。
[2]鞠尊洲。3S技术在水利工程测量中的应用[J]。科技信息 （学术研究） ， 2008 （14） ：239—240。
[3]李蕾。浅谈在水利工程建设管理中信息技术的应用[J]。建设科技， 2016 （17） ：93。
[4]王芳， 卜倩倩。探析计算机地理信息系统在水利中的应用[J]。资源节约与环保， 2015 （03） ：154。
[5]赵晓敏， 李本怀。GIS技术在水文预报中的应用—以长春新立城水库为例[J]。水利技术监督， 2013， 21 （02） ：64—66。
[6]刘丹， 刘萍， 樊耔均。GIS在水利现代化中的应用和发展趋势[J]。硅谷， 2012， 5 （24） ：43+42。
[7]张会玲。GIS技术应用及未来发展趋势[J]。信息系统工程， 2012 （04） ：98—99。
[8]董凤服， 刘洋。GIS技术的发展趋势[J]。网络与信息， 2007 （08） ：65。
译自 Environment Geo1ogy,2003(43):503～512。
Mowen Xie1Tetsuro Esaki1Guoyun Zhou1Yasuhiro Mitani1著
张晓娟2译 罗靖筠2校 朱汝烈2复校
（1Environmental System Institute,Kyushu University,Hakozaki 6-10-1,Higashi Ku,Fukuoka,Japan；2中国地质调查局水文地质工程地质技术方法研究所，河北保定，071051）
【摘要】本文在传统的边坡稳定性三维分析模型的基础上，提出了一个全新的基于GIS的边坡稳定性三维栅格分析模型。在这个模型中，假定初始滑动面就是椭球底面，采用蒙特卡洛（Monte-Carlo）随机模拟方法，在求取最小安全系数法的同时，确定出最危险滑动面。运用GIS栅格模型和GIS数据模拟滑坡滑动过程时，滑坡体将沿主滑方向滑动，直到其安全系数上升到1为止。所有的计算均可通过一个称为三维边坡地理信息系统（3DSLOPGIS）的计算程序来完成，该程序主要利用GIS的空间数据处理分析功能。
【关键词】确定性模型 地理信息系统（GIS） 蒙特卡洛（Monte-Carlo）模拟 滑动模拟 三维边坡稳定性
1 引言
滑坡不稳定性和风险评价不但已成为地学家和工程专家们感兴趣的主要课题，同时也成了世界各地政府部门和管理者关注的焦点。据统计世界上每年约有600人葬身于滑坡灾害中。在许多发展中国家，自然灾害所带来的经济损失，占总国民生产总值的1%～2%。
近年来，由于地理信息系统具有强大的空间数据处理功能，被广泛运用于自然灾害评价领域。GIS是由硬件和软件组成的系统，它可以实现数据采集、输入、操作、转换、可视化、组合、质疑、分析、建模和输出等过程。GIS对空间数据具有强大的分析和处理功能。同时，基于GIS的地质技术分析模型，可以简便而有效地分析滑坡稳定性。目前它已经被广泛地用于土木工程和地质工程中，进行边坡稳定性的分析。
我们通常认为一个传统的模型无论是对均质滑坡还是非均质滑动都是适用的。稳定性指数是被广泛应用的、基于岩土工程模型和物理力学参数的安全系数。安全系数的计算需要几何数据、剪切强度数据及孔隙水压力数据，正确的结果取决于可靠的数据和恰当的模型。尽管输入的数据会较大程度地影响安全系数，但一个可靠的确定性模型对于取得可靠结果则更为重要。确定性计算可在GIS系统内执行，也可利用其他程序完成。若使用其他程序计算，则GIS只作为一个空间数据库用来存储、显示、更新输入数据。此方法主要优点是利用外部模型计算可以节约时间；而其缺陷是对从外部模型获得的数据进行转化时较为复杂。因为每一个程序都有其自己的数据格式和数据结构，数据转换成为一个主要的问题。有些程序的输入模块只允许人工输入数据。只有当这些程序所默认的数据格式都是 ASCII码时，数据转换才可直接进行。运用外部模型的另一个缺点是计算结果通常不是按GIS的空间分布模式来表达，而是以点或线的形式表述的。因此，改变这种计算结果的表达形式也是个主要的问题。
用来计算安全系数稳定性模型的边坡是二维或三维的。因为一个地区包括很多边坡，而且必须分别对每个边坡做分析，所以利用这些模型计算安全系数的空间分布非常花费时间。要克服数据转换的困难，可以利用GIS内部确定性计算模型来实现。然而这一方法也有缺点，那就是由于应用复杂算法、迭代过程及在常规二维 GIS中的三维体积等复杂局限性，使得只有简单的模型能较容易实现。当前，只有基于GIS的无限边坡模型能分别计算出每个像元的安全系数。研究表明，只有当越来越多的成熟的三维模型和GIS系统得到使用后，才能彻底解决这类问题。
从近来对 GIS用于边坡稳定性分析的调查中发现，大部分研究者潜心于运用统计学方法来确定边坡破坏与影响因素之间的关系。尽管GIS能对区域数据进行了准备和处理，但是只有极少量的研究者运用了GIS的集成功能和边坡稳定性的确定性模型。
即使在很短的距离范围内，边坡破坏在空间上都有其不同的几何结构。因而，运用三维模型分析边坡稳定性是合理的。从20世纪70年代中期以来，三维稳定性模型的发展和运用日益受到关注。在地质力学的著作中提到了几个三维分析方法。
上面提到的大部分方法都用到了柱状图法。这些方法将柱体之间的作用力，或者说作为三维安全系数计算的假定前提，都忽略不计。因为所有与斜坡相关的GIS数据都可转成栅格数据，所以这些基于三维模型的柱体，就可能借助于使用GIS栅格数据用来进行三维稳定性的计算。然而，长期以来大家习惯采用人尽皆知的“一维模型”——“无限斜坡”模型，来描述滑动面与地面平行的长期天然边坡的潜在危险性。这样的模型仅仅可以用于浅层斜坡失稳分析和一些存在深层滑坡的区域性研究。
由于算法复杂、步骤重复和三维数据在二维GIS中难于表达，早期的文献中并没有提及三维确定模型的应用。为了克服 GIS数据的外部转换和GIS内部算法复杂等困难，此次研究中，在GIS软件组件（a GIS component）中使用了Visual Basic程序。三维因子的计算和滑动过程的模拟由计算机内的三维边坡地理信息系统（3-DSLOPGIS）的计算程序完成。在这个系统中，GIS组件（ESRI公司生产的MapObjects2.1）可以完成所需的GIS功能，就像普通的GIS软件一样，它可以有效的管理和分析所有与滑动相关的数据。所有用来计算三维斜坡安全系数的数据都采用GIS的数据格式（例如矢量和栅格数据层），因此，没必要在GIS数据格式和其他程序的数据格式之间进行数据转换；同时，复杂算法和三维问题的交互程序也可以理想的实现。
在此次研究中，将基于GIS栅格数据和基于柱状图的三维边坡稳定性分析模型相结合（Hovland,1977），演绎了一个新的基于GIS栅格的三维确定性分析模型。
运用蒙特卡洛随机模拟方法求最小安全系数值，从而确定临界滑动条件。假定基本滑动面是一椭球体的较低部分，临界滑动则受不同地层受力情况和不连续界面状况的影响而变化。客观事物的这种变化引出最小三维安全系数。
如果滑坡的三维安全系数小于1，滑坡就有滑动的危险，那么评估滑坡灾害的规模和影响范围是非常重要的。因此，在此研究中，采用基于GIS三维栅格数据模型和GIS栅格数据来模拟滑坡滑动过程的目的，就是评估滑坡危险性和预测其影响范围。
2 基于GIS的三维模型
利用GIS的空间分析功能，所有与三维安全系数计算有关的输入数据（如高程、倾向、坡度、地下水、地层、滑动面和力学参数等）都有其对应的栅格元，而所有与斜坡相关的数据都是栅格化的。当这些数据输入到确定的边坡稳定性模型中时，就可计算出一个安全系数值。下面在Hovland模型的基础上，详细介绍基于GIS的三维模型。在这个模型中，考虑了孔隙地下水压力，所有输入数据都能简单地转换成栅格数据。
图1是具有潜在滑动面的滑体的三维几何示意图。滑坡的稳定性与地质岩层、地貌、地质力学参数和水动力条件有关。
图1 边坡坍塌三维景观
图2所示是土壤（或岩石）小柱状研究体物质的离散性。所有与滑坡相关的数据都可用如图2所示的柱状三维可视图来表示。假定每一个柱体单元的垂面均为无摩擦面（柱体单元的垂面不受其他边界影响，或其影响可忽略不计），三维安全系数可用公式（1）表示：
地质灾害调查与监测技术方法论文集
式中：F3-D为三维斜坡安全系数，W为一个柱体的重量，A为滑动面面积，c为内聚力，φ为内摩擦角，θ为滑动面的角度，而J、I为在斜坡破坏范围栅格内的行列数和柱体数。如果没有GIS，则基于柱体模型的三维安全系数的计算将是冗长且耗时的工作，数据的更新和增加也极其不便。然而，在GIS中，通过运用GIS空间数据处理与分析功能，整个研究区的边坡稳定性相关数据可用如图3所示的矢量图层来描述；而对于每一层，则可通过GIS空间数据处理与分析功能得到栅格数据，其像元大小可根据精度需要而定。
图2 滑动面和三维棚格柱状图
现在，将斜坡破坏划分为基于栅格数据的柱体。参考图2，诸如地表、地层、地下水、裂缝和滑动面之类的空间数据均可从栅格数据层中得到。因为与斜坡相关的数据量非常大，所以不能高效的管理所有的栅格数据集。因此，在三维边坡地理信息系统中，有一个专门储存这些栅格数据的点数据库，其中，有一个属性表用来链接所有与滑动相关的数据。每个栅格柱状图的中心点设置点类型，其他区域则设置与滑坡相关的一些数据（例如地面高程、地层和裂缝的高程、地下水、滑动面的深度等等）。表1所示即是属性表的一个实例。
图3 边坡稳定性分析GIS图层
表1 点数据库的实例描述
另一方面，为了控制滑坡边界和有效管理空间数据并进行分析，滑坡的边界线被定义为多边形类型文件。
基于这种点数据库，公式1可以改成基于GIS的方程。这里所有的阻力和滑力都是沿着滑动方向的，而不必如 Hovland的模型所用的Y轴方向。在本研究中，假定斜坡区域的主要倾斜方向为可能滑动方向。根据图4，滑动表面面积可由公式（2）得到。
地质灾害调查与监测技术方法论文集
从图4推导出如下公式：
地质灾害调查与监测技术方法论文集
地质灾害调查与监测技术方法论文集
接着，x和y轴的倾角推导如下：
地质灾害调查与监测技术方法论文集
记α＝cellsize/cosθxz和b=cellsize/cosθyz，则一个栅格柱状图的滑动面面积为：
地质灾害调查与监测技术方法论文集
滑坡范围主滑动方向的倾角计算公式如下：
地质灾害调查与监测技术方法论文集
至此，三维边坡水平滑动方向安全系数可以用下面的公式计算：
地质灾害调查与监测技术方法论文集
图4 三维安全因子推导公式的一个栅格柱状图
这里，对于每个栅格，Zji,zji分别为地表高程和滑动面高程，uji为在滑动面上的孔隙水压力，而 γ′为单位重量。
为了检验基于栅格的GIS三维稳定分析模型，我们运用这个模型做了一个实例计算。实例问题为一个均质的粘土滑坡，具有球形滑动面，其他各种参数如图5所示。在图5中，c为内聚力，φ为摩擦角，R为瞬时摩擦力，γ为土的单位重量。运用封闭式（closed-form）算法得出三维安全系数为1.402。运用CLARA模型算得安全系数为1.422。同样的问题运用三维边坡模型算得三维安全系数范围为1.386到1.472，它取决于用于被分离的边坡柱体的数量。
图5 实例问题验证
运用基于GIS栅格的三维稳定分析模型（图5），并将格网尺寸定为0.5m时，算得三维安全系数为1.386；而当格网尺寸为0.6m时，算得安全系数为1.388。很明显，与封闭式算法相比，基于栅格模型的GIS可有效的用于三维边坡稳定性评估。
3 确定临界滑动表面和蒙特卡洛模拟
滑动面只能通过岩土工程调查来确定，由于地质调查的费用比较昂贵，因此滑动面通常是很难确定的。因此，边坡稳定性评价对临界滑动面的确定是非常重要的。
为了判定三维临界滑动情况，利用蒙特卡洛随机模拟方法来计算三维安全系数最小值。假定最初的滑动面是一个椭球体的较低部分，边坡表面则根据不同地层受力情况和不连续界面条件而改变。最终得到危险滑动面，同时可得到相关三维安全系数的最小值。
4 椭圆坐标转换
假定最初的滑动面是一椭球体的较低部分，椭球体的倾斜方向设置为与研究区主要的倾斜方向一致；将椭圆的倾角基本上设定得与研究区起伏变化的倾角接近。其主倾向为α，主倾角为β，它们是由边坡破坏区域主要栅格像元的值确定的。假定倾向和倾角属正常分布，则将主倾向α和倾角β代入分布模型中：
地质灾害调查与监测技术方法论文集
运用公式（10）和（11）完成坐标转换。图6显示了坐标转换过程。
图6 坐标转换过程
地质灾害调查与监测技术方法论文集
地质灾害调查与监测技术方法论文集
式中：x、y、z为全球大地坐标， 为当地坐标，x0、y0、z0为椭球体中心点坐标。
5 Z值的确定和滑动面的倾斜度
滑动面上“B”点的Z值是根据直线 AB和椭圆，由公式（12）计算的结果确定的（见图7）。
地质灾害调查与监测技术方法论文集
对于每个栅格像元，滑动面的倾向和倾角可通过下面的公式计算得出，像元（j,i）的倾角可以通过图8中点1～4的Z值来确定。点1～4的值由公式（13)(14)(15）算出，滑动面的倾向和倾角由公式（16）算出。
图7 确定滑动面上的Z值
图8 滑动倾角的计算
地质灾害调查与监测技术方法论文集
地质灾害调查与监测技术方法论文集
地质灾害调查与监测技术方法论文集
地质灾害调查与监测技术方法论文集
这里，Z(j,i）为像元（j,i）的Z值，θ为倾角，β0是相对于X轴的倾向。在GIS中，倾向是与 Y轴之间的夹角。因此，当最高点是点3时，倾向是90－β0；当最高点是点4时，倾向是90＋β0；当最高点是点2时，倾向是270－β0；当最高点是点1时，倾向是270＋β0。
6 随机模拟
为了确定临界滑动面，蒙特卡洛模拟通常用于为三维边坡稳定性分析选择变量。这些变量是椭球体的中心点、几何参数和倾角。椭球体的中心点作为研究区的中心点需要首先确定，然后在一个确定的范围内随机选择。
椭球体的几何参数a、b、c是由用户在一定范围内随机设定的，确定范围如公式（17）：
地质灾害调查与监测技术方法论文集
假定a,b,c都均匀分布，则蒙特卡洛模拟的随机变量由公式（18）和（19）来算出。
在[0,1]范围内平均分布的随机变量可通过全等乘积方法得出：
地质灾害调查与监测技术方法论文集
地质灾害调查与监测技术方法论文集
式中：ri为在[0,1]范围内平均分布的随机变量。在[a,b]范围内平均分布的随机变量可由公式（19）计算得出。
地质灾害调查与监测技术方法论文集
式中：xi为在[a,b]范围内平均分布的随机变量。
椭球体的倾角设定为平均分布的一个随机变量。平均分布范围为主倾角及其在一个确定的波动范围之内变化的变量。
7 计算三维安全系数最小值的过程
整个研究区（或边坡破坏范围）可以被均分为若干小矩形栅网，如同基于栅格的GIS一样。关于基于栅格的三维边坡稳定性分析的数值计算，所有的计算过程都可以通过前面提到的Visual Basic（利用GIS组件）来完成。这个软件叫三维边坡地理信息系统，是运用 Visual Basic 6.0和ESRI公司生产的MapObjects 2.1开发的。MapObjects作为GIS的一个组件，用来对GIS数据进行组织和空间分析。计算三维安全系数的过程如图9所示。
图9 三维安全因子最小值计算过程
在这个过程中，数据模块的功能用来获得所有与边坡相关的地质、地貌、水动力学数据和地质力学参数；随机变量参数模块用来随机选择蒙特卡洛模拟的实验滑动面；三维边坡稳定性模块可用于计算三维安全系数；而危险滑动面及其安全系数可以通过一些实验计算得出。在图9中可以看到，关于GIS空间分析功能的所有模块可以通过GIS组件来实现。因为一个GIS组件是在三维边坡地理信息系统系统中完成的，所以可以有效地计算三维安全系数；同时利用与边坡相关的GIS数据，所有的相关数据和结果可以在三维边坡地理信息系统系统中实现可视化。
实例剖面如图10所示。在这个实例中考虑的因素有：4个地层、地下水和破坏面；其物理和力学参数如表2所示。
表2 研究实例的物理和地质力学参数
图10 断层面研究实例
图11 计算次数与最小三维安全因子实验
为确定临界滑动面，对蒙特卡洛随机计算次数进行了实验，总共计算次数达到了1000次。每次实验计算的三维安全系数最小值的结果如图11所示。图中明确显示在实验计算了300次后，得到的安全系数最小值。这300次实验的结果见图12，这些计算结果差别不太大，其最小值为1.34，最大值是1.68。这个临界滑动的研究程序是建立在最小安全系数的计算基础之上的。而最小安全系数的计算结果取决于参数的随机选择。有关这一临界滑动实例的三维可视图见图13。通过三维模型与二维模型结果的比较，用Janbu法确定临界滑动面时，使用的是图10所示的二维模型和表2所列的参数，通过这种二维模型计算出的安全系数为1.18，这要比用三维模型计算出结果的极小值（1.346）略小一点。
图12 三维安全因子分布曲线
8 滑坡滑动过程模拟
基于GIS栅格三维边坡稳定性分析模型和GIS栅格数据，对滑坡滑动过程进行了模拟，直到三维安全系数大于1为止。滑动方向按滑动面的主滑方向确定。图14中展示了由滑动面确定的八个滑动方向。例如，若滑面方向的倾角在22.5°～67.5°之间，则滑坡将要滑动的方向恰在该图的右上方（即“5”方向）。
图13 临界滑动面三维展视图
图14 滑动面的滑动主倾向
图15 滑坡滑动过程模拟流程方框图
滑坡滑动过程的模拟流程见图15。首先，要计算滑坡初始状态时的三维安全系数，以确定其滑动的可能性。若其安全系数小于1，则接着进行下一步滑动过程模拟。先沿着由滑面主倾向确定的滑动方向移动滑坡多边形；接着，在新的滑坡多边形范围内，分步（每一步等于一个栅格大小）计算每一个栅格的DEM和滑动的变化，并再次计算下一步滑动的新滑动方向。并在新的DEM数据和滑动多边形范围的基础上，计算出新的三维安全系数。如果三维安全系数仍然小于1，则进行以下的新滑动步骤模拟。
在这种滑动模拟模型中，假定滑动面内摩擦角不改变，但除了在初始三维边坡安全系数的计算过程之外，假定滑动面没有内聚力（即内聚力为零）。
仍然用同样的实例（如图5所示），用不同的两种动力学参数进行滑坡滑动过程模拟：
情况1:c=4kN/m2，φ＝110,y=23kN/m3
情况2∶c=6kN/m2，φ＝10.5°，γ＝23kN/m3
第一种情况下，初始边坡安全系数为0.82，在进行7步滑动之后，滑坡体开始趋于稳定，其安全系数是1.04。部分滑动步骤剖面及三维视图变化如图16所示。在此图中，DEM的改变及滑坡体移动过程一目了然。运用三维边坡地理信息系统，也可将可视滑动过程表现为GIS地图和剖面图的形式。滑坡体沿水平方向的最终滑动距离为3.0m。
图16 不同滑动阶段的地表和剖面三维视图
第二种情况下，滑坡体将一直向下滑动到平坦地区，水平方向滑动距离为14m。滑坡体最后停止滑动位置的三维展视图如图17所示。
图17 滑坡体最后停止位置
9 讨论和结论
在三维边坡稳定性柱状分析模型的基础上，开发了一个全新的基于GIS栅格的三维确定性模型，并且通过一个问题实例证实了其正确性。在三维边坡稳定性分析模型中，假定其初始滑面为一椭球面；其三维临界滑面，是利用蒙特卡洛随机模拟求取最小三维安全系数而确定的。基于GIS的栅格三维模型，滑坡滑动过程模拟用于判断滑坡灾害和预测滑动距离。已开发了作为计算程序软件的三维边坡地理信息系统，它足以完成一切有关三维边坡问题的计算，其中的GIS组件用于实现GIS的空间分析功能和有效数据的管理。因其具有空间分析、数据管理和与边坡相关的综合数据的GIS可视化等优点，所以三维边坡稳定性问题已经比较易于研究。自打全新的基于GIS栅格三维边坡稳定性分析模型问世，就为惯于使用传统数学方法研究边坡稳定性的工作者拓展了一个新的研究领域和数据库方法。
基于CAD和GIS的校园绿地信息管理的研究
基于GIS的城市土地利用与交通网络一体化模型
基于GIS和RS的应用分析系统集成研究
基于GIS技术的社区综合信息管理系统
基于组件式GIS的模拟导航电子地图设计
浅谈电信部门对电子地图的要求和制作
物流园区信息平台规划与设计
......
另外，可以看看《第X届ArcGIS暨ERDAS中国用户大会论文集》，上面有很多，
我这有第七届的。
一、国际发展状况 地理信息系统的存在与发展已历经30余年。用户的需要、技术的进步、应用方法论的提 高，以及有关组织机构的建立等因素，深深地影响着地理信息系统的发展。 综观GIS发展，尤其是北美地区的实际情况，可将地理信息系统发展分为以下几个阶段： （1）60年代为地理信息系统开拓期，注重于空间数据的地学处理。例如，处理人口统计 局数据（如美国人口调查局建立的DIME）、资源普查数据（如加拿大统计局的GRDSR）等 。许多大学研制了一些基于栅格系统的软件包，如哈佛的SYMAP、马里兰大学的MANS等。 综合来看，初期地理信息系统发展的动力来自于诸多方面，如学术探讨、新技术的应用 、大量空间数据处理的生产需求等。对于这个时期地理信息系统的发展来说，专家兴趣 以及政府的推动起着积极的引导作用，并且大多地理信息系统工作限于政府及大学的范 畴，国际交往甚少。 （2）70年代为地理信息系统的巩固发展期，注重于空间地理信息的管理。地理信息系统 的真正发展应是70年代的事情。这种发展应归结于以下几方面的原因：一是资源开发、 利用乃至环境保护问题成为政府首要解决之疑难，而这些都需要一种能有效地分析、处 理空间信息的技术、方法与系统。二是计算机技术迅速发展，数据处理加快，内存容量 增大，超小型、多用户系统的出现，尤其是计算机硬件价格下降，使得政府部门、学校 以及科研机构、私营公司也能够配置计算机系统；在软件方面，第一套利用关系数据库 管理系统的软件问世，新型的地理信息系统软件不断出现，据IGU调查，70年代就有 80 多个地理信息系统软件。第三，专业化人才不断增加，许多大学开始提供地理信息系统 培训，一些商业性的咨询服务公司开始从事地理信息系统工作，如美国环境系统研究所 （ESRI）成立于1969年。这个时期地理信息系统发展的总体特点是：地理信息系统在继 承60年代技术基础之上，充分利用了新的计算机技术，但系统的数据分析能力仍然很弱 ；在地理信息系统技术方面未有新的突破；系统的应用与开发多限于某个机构；专家个 人的影响削弱，而政府影响增强。 （3）80年代为地理信息系统大发展时期，注重于空间决策支持分析。地理信息系统的应 用领域迅速扩大，从资源管理、环境规划到应急反应，从商业服务区域划分到政治选举 分区等，涉及到了许多的学科与领域，如古人类学、景观生态规划、森林管理、土木工 程以及计算机科学等。许多国家制定了本国的地理信息发展规划，启动了若干科研项目 ，建立了一些政府性、学术性机构。如中国于1985年成立了资源与环境信息系统国家重 点实验室，美国于1987年成立了国家地理信息与分析中心（NCGIA），英国于1987年成立 了地理信息协会。同时，商业性的咨询公司，软件制造商大量涌现，并提供系列专业性 服务。这个时期地理信息系统发展最显著的特点是商业化实用系统进入市场。 （4）90年代为地理信息系统的用户时代。一方面，地理信息系统已成为许多机构必备的 工作系统，尤其是政府决策部门在一定程度上由于受地理信息系统影响而改变了现有机 构的运行方式、设置与工作计划等。另一方面，社会对地理信息系统认识普遍提高，需 求大幅度增加，从而导致地理信息系统应用的扩大与深化。国家级乃至全球性的地理信 息系统已成为公众关注的问题，例如地理信息系统已列入美国政府制定的“信息高速公 路”计划；同美国副总统戈尔提出的“数字地球”战略、我国的“21世纪议程”和“三 金工程”也包括地理信息系统。毫无疑问，地理信息系统将发展成为现代社会最基本的 服务系统。二、国内发展状况 我国地理信息系统方面的工作自80年代初开始。以1980年中国科学院遥感应用研究所成 立全国第一个地理信息系统研究室为标志，在几年的起步发展阶段中，我国地理信息系 统在理论探索、硬件配制、软件研制、规范制订、区域试验研究、局部系统建立、初步 应用试验和技术队伍培养等方面都取得了进步，积累了经验，为在全国范围内展开地理信息系统的研究和应用奠定了基础。 地理信息系统进入发展阶段的标志是第七个五年计划开始。地理信息系统研究作为政府 行为，正式列入国家科技攻关计划，开始了有计划、有组织、有目标的科学研究、应用 实验和工程建设工作。许多部门同时展开了地理信息系统研究与开发工作。如全国性地 理信息系统（或数据库）实体建设、区域地理信息系统研究和建设、城市地理信息系统 、地理信息系统基础软件或专题应用软件的研制和地理信息系统教育培训。通过近五年 的努力，在地理信息系统技术上的应用开创了新的局面，并在全国性应用、区域管理、 规划和决策中取得了实际的效益。
