基于差分进化算法和降阶模型的渗流场反问题研究

来源 :西安理工大学 | 被引量 : 1次 | 上传用户:pisces_ww
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
中国水能资源的开发已逐渐向西部偏远地区推进,在建和拟建的大型水利水电工程坝址区多位于河谷深切、地质条件复杂的西部地区,查清库区渗流问题对水利水电工程的建设和安全管理十分重要。岩土体渗透系数是控制地下水渗流特征的关键参数之一,未知的渗透系数会对地下水模拟的可靠性产生严重影响。由于工程岩体的渗透性常具有空间变异性,仅依靠传统的现场试验法已不能满足工程需求。利用逆模型进行参数估计是地下水模拟的一个重要组成部分。作为一种典型的反演方法,数值模型和优化算法相结合的模拟-优化方法需多次调用数值模型,以对大量随机生成的候选解进行评估。即使使用高速处理器,参数反演也是一项非常耗时且计算量大的任务。本文针对模拟-优化方法的高耗时问题进行了研究,在尽可能减小引入误差的前提下,从优化算法、反演参数和数值模型三方面研究了减少模拟-优化方法时间成本的方法。由于水头为渗透系数的非线性函数,本文使用差分进化算法作为参数反演的优化算法。主要研究内容和成果如下:(1)阐述了模型校准与参数反演的关系,给出了有限元软件ADINA与优化算法结合的方法,建立了估计渗透系数的模拟-优化模型(ADINA-MMRDE)。通过一个算例阐述了参数灵敏度分析在参数反演中的重要性,研究了不同目标函数、测量误差和种群大小对ADINA-MMRDE模型性能的影响。结果表明,目标函数对ADINA-MMRDE影响甚微,ADINA-MMRDE对测量误差非常敏感。相比ADINA与其他优化算法结合的ADINA-DE和ADINA-PSO模型,ADINA-MMRDE模型反演精度更高,能更快、更稳定地搜索全局最优解。(2)针对经典差分进化算法的变异策略收敛速度慢、全局收敛性不佳及算法停滞等问题,提出了一种兼具局部与全局收敛性能的新型变异策略。基于该变异策略,进一步提出了一种基于轮盘赌选择的多种变异策略的差分进化算法(MMRDE)。经49个测试函数测试,结果表明,与一些改进的差分进化算法相比,MMRDE能在探索和开采之间取得更好的平衡。(3)为了在保证模拟精度的前提下减少模型的计算时间,阐述了基于投影法的降阶模型技术(本征正交分解法和贪心样本法)的降阶机理、构建步骤和误差估计方法。改进了贪心样本法的迭代终止条件,比较了本征正交分解法和贪心样本法的计算成本,以及二者在参数集、网格密度和参数数量方面的性能表现。结果表明,当样本规模较少时,不同的样本集生成方法对降阶精度影响较大;单元尺寸影响降阶模型的构建时间,但对降阶模型的精度影响不大;反演参数越多,降阶模型的省时优势越明显。(4)针对将模型降阶技术应用于参数反演中的一些关键性的程序设计难点,设计了一种集识别反演参数、矩阵分块技术以及边界处理于一身的渗透矩阵处理程序,设计了一套高效的内存存储方案以解决使用传统有限元的Skyline稀疏存储格式可能导致的内存不足问题。针对钻孔位置不在网格节点上时的水头计算问题,提出了基于本文提出的MMRDE算法的有限元插值程序插值计算钻孔处的水头。设计了基于降阶模型的参数反演程序,使用算例测试了其的反演精度、对观测误差的敏感性与时间成本。结果表明:推荐采用训练参数规模为500的贪心样本法用于参数反演;基于降阶模型与基于原始模型的参数反演程序对误差的敏感程度以及反演精度非常相近,但耗时差别较大;同等计算能力条件下采用算例中的三维模型时,使用降阶模型的参数反演程序的反演时长约为使用全阶模型的16.67%,因此能明显的节省时间成本。(5)将基于ADINA模型与基于降阶模型的反演程序共同用于估算某水电站坝基岩体的渗透系数,这两种反演程序都集成了本文提出的MMRDE算法。建立了初始渗流场分析模型(反演模型)来估算渗透系数,在反演模型的基础上建立了工程运行期模型以验证反演效果。共有20个勘探期钻孔水位和13个大坝监测孔水位数据,前者用作参数反演的观测数据,后者用于验证反演结果。结果表明,两反演模型的反演结果相差较小,但基于降阶模型的反演程序的时间成本远小于基于ADINA模型的反演程序(维数为6和13时,反演时长分别约节省19.1和21.4倍)。因此,使用本文提出的MMRDE算法作为优化算法时,降阶模型可替代原始模型用于执行大型工程的初始渗流场的反演任务。
其他文献
目的:观察分析西吡氯铵含漱液治疗牙龈炎的临床疗效。方法回顾性分析我院于2012年5月-2013年9月收治的90例牙龈炎患者的病理资料,将患者随机的分为对照组和治疗组,对照组和治
本文概述了松软煤层炮掘金梯支护破岩托顶煤施工方法,分析了顶煤容易冒落的原因,提出了防治措施及注意事项。
探讨了在Delphi中控制Excel电子表格的方法和技术,以及Excel表导入数据库、数据库内容导出为Excel电子表格的方法.
由于钱家营矿业分公司矿井地热及大功率设备应用,工作面空气温度有时超过规程规定,极大影响了员工的身心健康和制约矿井生产效率提高,为保证员工生命健康,提高生产率,针对工作面空
目的:观察疏肝降逆方治疗功能性消化不良的疗效。方法:180例随机分为治疗组90例,给予疏肝降逆方:柴胡、枳壳、党参、山药、茯苓、陈皮、石菖蒲各12g,白术、黄苓、郁金、法半夏、广
目的:观察痰饮膏穴位外敷对慢性支气管炎的治疗效果。方法:痰饮膏(干姜、花椒、细辛、肉桂、桂枝等)从初伏开始,分别敷于双侧肺俞穴及膻中穴处,48h后去掉药膏,间歇24h重复敷
目的:分析脑出血患者在急性期的相关并发症的预防措施和护理体会。方法选取我院近2年间收治的脑出血急性期患者共60例的临床资料进行研究,患者在入院后密切观察病情变化,进行基
为了给葡萄倍性选育工作提供理论依据和科学参考价值,对不同葡萄品种二、四倍体的形态特征、叶绿素含量、原花色素含量及POD、CAT、MDA含量进行了测定、分析。结果表明:(1)不同
逼仓事件是期货市场最典型、造成损失最大的风险事件。监控逼仓事件是期货市场风险管理的核心任务之一,建立期货逼仓监控指标体系是其中的首要环节。本文尝试建立了一套监控指
证明了如下抛物方程解的整体存在性和梯度估计ut-div{σ-(| △↓u(t)}+g(t,x,u,△↓u)=Ω×[0,∞)。初边值条件为u(x,0)=u0∈W0^1,p0(Ω),u(x,t)/aΩ=|0,0≤t<t+∞。其中:力是Rw的有界区域,σ(|△↓