与弦搜寻演算法在土石坝渗流破坏中应用

《与弦搜寻演算法在土石坝渗流破坏中应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、与弦搜寻演算法在土石坝渗流破坏中应用　　摘要：根据土石坝坝体破坏案例所累积的经验，除坝顶溢流溃决破坏外，坝体渗漏亦为造成管涌破坏主要原因之一。管涌会造成坝体突然、全面的破坏，后果非常严重，本研究最佳化演算法采用和弦搜寻法，分析中以MATLAB作为和弦搜寻最佳化流程的主控伺服，常用的地工分析软体FLAC则扮演计算引擎的角色。以评估本研究所研拟的和弦搜寻演算法应用于土石坝渗流问题诊断的功效，以说明此方法的可行性与有效性。关键词：土石坝；渗流；和弦搜寻演算法中图分类号：TV223.4文献标识码：A文章编号：1006-4311（2014）05-0070-020引言我国某些地区地形

2、陡峭、河川短促，降雨后雨水迅速流入海洋，又因降雨四季分配不均，以致常面临缺水危机，藉由建造水库以储蓄水源乃是最直接及不得已的水源开发手段。受限于地质条件，水库最常见的坝体型态为土石坝，以堆填的土石堤坝发挥阻水功能以蓄水。7过去国外不乏土石坝破坏的案例，土石坝常见的破坏主要原因包括库水溢顶溃决、管涌或边坡滑动破坏。溢顶溃决一般因为集水区带来流量超过本身设计容量、或因管理不当或坝体严重沉陷等因素，导致水流溢顶而直接冲刷坝体造成破坏，因此需透过水文分析、水库管理与操作来降低风险。边坡滑动破坏需检讨不同条件，包括刚完工、急泄降、稳态渗流、地震力影响等情境下的稳定状况，透过位移监测

3、与完整性检查及检测来评估其稳定性。管涌或边坡滑动破坏则系导因内部结构的破坏。土石坝不可能完全无渗漏水，但若渗漏水过多，一则造成库容的损失，更严重的是异常渗流水流入坝体内部、过大的水力坡降可能导致坝体内部冲蚀，逐渐扩大造成管涌终究可造成坝体全面的破坏，后果非常严重，日常需要透过渗漏量与水压监测数据来研判安全性。监测数据若有异常、其异常的真正原因有必要透过坝体内部问题诊断来寻求答案，对于水库安全是相当重要的课题，但是如何诊断坝体内部可能出现的问题却非常棘手。本研究以土石坝渗流相关问题为例，运用监测数据，配合数值模式模拟分析，建构坝体异常问题诊断与诠释的方法。1常见启发式最佳化

4、演算法启发式演算法与一般数值最佳化演算法最大不同在于求解方式为有系统的多次随机搜寻直到达目标条件为止，而非以计算导数值方式搜寻求解。启发概念即套用于随机搜寻模式中，并非毫无头绪的随机组合；而启发顾名思义就是对于人事物现象的观察有所启示，继而学以致用。7这类方法最大优点是：求解效果佳，且相对于传统求解最佳化法上，较不受限于初始参数值如何明确给定的问题，并以随机方式搜索得到的解较灵活、较能跳脱局部解情形，因此也不用记录求解过程大量函数导数资料，应用范围也较广。以下简单介绍一些常见的启发式最佳化演算法：遗传演算法、模拟退火法、禁忌搜寻法、蚂蚁族群演算法（AntColonyOpt

5、imization，ACO）、粒子群法（ParticleSwarmOptimization，PSO）、和弦搜寻演算法。2和弦搜寻演算法2001年由Geem等人所提出发表，为一种启发式演算法进化版，发展至今已有10年之久。此方法模拟音乐家们使用不同乐器一起即兴演奏下，每位音乐家各自记忆所弹奏的曲调，并藉由每次合奏后来调音，因此，经过数次的即兴演奏后所演奏的音乐会越来越和谐、越美妙，依此种方式产生的最佳演奏概念来求解，称和弦搜寻演算法。2.1和弦搜寻演算法步骤步骤1：问题公式化（problemformulation）步骤2：参数设定（algorithmparameterset

6、ting）步骤3：随机产生初始记忆（randomtuningformemoryinitialization）步骤4：改善和弦：随机选取、考虑记忆与调音（harmonyimprovisation：randomselection，memoryconsideration，andpitchadjustment）7步骤5：记忆更新（memoryupdate）步骤6：满足终止目标（performingtermination）步骤7：终曲（cadenza）一开始将问题以数学模式表示，即公式化；接着给定参数值或相关参数设定后以完全随机模式产生一组解（和弦）作为初始和弦记忆向量（HM），以

7、上（步骤1～3）归类为问题的初始化；接下来（步骤4～6）为问题的搜寻方式，考虑参数机率（HMCR、PAR、bw）关系后随机产生一个解（和弦）；此新产生的和弦若是比HM中任何一个好，则剔除最差者并取代之；接着一直重复步骤4、5直至达终止目标为止。步骤7类似乐曲接近结尾的一段装饰奏，将最佳解（和弦）再做一次演奏或修饰后来收尾。其中，和弦搜寻演算法所使用的主要参数有两个：和弦记忆比率（HarmonyMemoryConsideringRate，HMCR）与调整比率（PitchAdjustingRate，PAR），并以两者参数关系来改善