地下水流并行有限层方法及同伦反演研究  

Parallel finite layer method for groundwater flow and its homotopy inversion algorithm

作　　者：王少伟 徐进[1] 杨伟涛 WANG Shao-wei;XU Jin;YANG Wei-tao(School of Civil Engineering,Yantai University,Yantai 264005,China)

机构地区：[1]烟台大学土木工程学院,烟台264005

出　　处：《计算力学学报》2020年第6期756-762,共7页Chinese Journal of Computational Mechanics

基　　金：国家自然科学基金(41602284)资助项目.

摘　　要：根据有限层求解格式存在的解耦性,实现了地下水三维流问题的高效并行化计算。在此基础上,结合非线性同伦方法,提出了地下水参数反演分析的并行同伦算法,利用MATLAB编译了相应的正反演计算程序。与已有解析解和有限差分解的对比以及数值算例,验证了并行化正反演方法及程序的正确性,探讨了并行算法的计算效率。研究表明,并行方法可以有效提高计算速度,较串行方法具有明显优势,同时同伦反演方法具有大范围收敛的特点,不依赖于参数值的初始选取。According to the decouplization of the finite layer scheme,efficient parallel computation of the three-dimensional groundwater flow problem is presented.On this basis,and further combined with the nonlinear homotopy method,a parallelized homotopy algorithm for groundwater parameter inversion analysis is proposed.The corresponding forward and inversion compute programs are developed in MATLAB.Through the comparison of existing analytical solutions,finite difference solutions and the verification numerical examples,the validity and computational efficiency of the parallel forward and inverse methods are demonstrated.The research shows that the parallel method can effectively improve the calculation speed,and has obvious advantages over the serial method.Moreover,the homotopy parallel inversion method has the characteristics of large-scale convergence and does not depend on the initial selection of parameter values.

关 键 词：地下水流 有限层法 并行计算 参数反演 同伦方法 

分 类 号：P641[天文地球—地质矿产勘探] O35[天文地球—地质学]