论文部分内容阅读
中国水能资源的开发已逐渐向西部偏远地区推进,在建和拟建的大型水利水电工程坝址区多位于河谷深切、地质条件复杂的西部地区,查清库区渗流问题对水利水电工程的建设和安全管理十分重要。岩土体渗透系数是控制地下水渗流特征的关键参数之一,未知的渗透系数会对地下水模拟的可靠性产生严重影响。由于工程岩体的渗透性常具有空间变异性,仅依靠传统的现场试验法已不能满足工程需求。利用逆模型进行参数估计是地下水模拟的一个重要组成部分。作为一种典型的反演方法,数值模型和优化算法相结合的模拟-优化方法需多次调用数值模型,以对大量随机生成的候选解进行评估。即使使用高速处理器,参数反演也是一项非常耗时且计算量大的任务。本文针对模拟-优化方法的高耗时问题进行了研究,在尽可能减小引入误差的前提下,从优化算法、反演参数和数值模型三方面研究了减少模拟-优化方法时间成本的方法。由于水头为渗透系数的非线性函数,本文使用差分进化算法作为参数反演的优化算法。主要研究内容和成果如下:(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算法作为优化算法时,降阶模型可替代原始模型用于执行大型工程的初始渗流场的反演任务。