张新启， 龚爱民， 徐兴倩， 罗聪聪， 谢 非

（云南农业大学水利学院，昆明 650201）

边坡安全系数的计算是评价其稳定性的基本依据，极限平衡分析方法仍是边坡防治工程中获取安全系数的主流方法. 其中，瑞典圆弧法具有简单有效，易于编程的优点，被广泛应用于滑面为圆弧形的土质边坡滑动面分析［1］. 采用该方法计算边坡安全系数，最关键的问题是如何从大量可行的计算滑动面中展开精准搜索，并最终确定可靠的边坡临界滑动面. 因此，边坡临界滑面搜索方法的优化研究对于边坡稳定性分析及防治具有重要意义.

目前，主流数值计算软件对圆弧滑面的搜索模块均使用的区格法，如SLIDE、STABL、SLOPE/W（Geostudio）等. 该方法通过n×n点的正方形网格确定临界滑面中心点的位置，计算不同区格点及其半径增量，从而定位安全系数最小的圆弧面［2］. 区格法在边坡圆弧滑动面搜索中应用广泛，求解严格，但该方法搜索范围广、计算量大、受区格划分精度影响不能保证其得到的破坏面就是全局最小值［3］. 因此，搜索空间存在多个局部最小点是定位边坡临界滑面研究的难点问题.

针对此问题，许多学者借助优化算法开展边坡临界滑面的定位研究. 国内外最常用的方法有遗传算法［4-5］、模拟退火算法［6-7］、蚁群算法［8-9］、粒子群算法［10］、鱼群算法［11］等，以上方法均有其各自的优势，但同样也存在些许不足. 例如：遗传算法易于早熟陷入局部最优解，对新空间的探索能力有限；模拟退火算法全局搜索性能受初始值设置影响较大，收敛速度慢；……