用同伦方法反演流体饱和孔隙介质的参数.pdfVIP

第23 卷 第1 期 岩石力学与工程学报 23(1) ：129～136 2004 年1 月 Chinese Journal of Rock Mechanics and Engineering Jan. ，2004 用同伦方法反演流体饱和孔隙介质的参数* 1 1 2 韩 华 章梓茂 魏培君 (1 北方交通大学土建学院 北京 100044) (2 北京大学力学系 北京 100871) 摘要 从材料响应的理论合成应与实际测量数据相拟合这一出发点，将参数反演问题转化为非线性方程组的零点求 解问题，从而应用一种大范围收敛的同伦方法来求解非线性方程组，并把这种方法用于 Simon1984 年给出的具有解 析解的一维双相介质模型的数值模拟。数值模拟实例的结果表明了同伦方法的可行性和稳健性。 关键词 数值模拟，同伦法，双相介质，反演，大范围收敛 分类号 O 189.23 文献标识码 A 文章编号 1000-6915(2004)01-0129-08 HOMOTOPY METHOD FOR INVERSING PARAMETERS OF WAVE EQUATION IN POROUS MEDIA 1 1 2 Han Hua ，Zhang Zimao ，Wei Peijun 1 ( Department of Civil Engineering ，Northern Jiaotong Uuniversity ， Beijing 100044 China ) (2Department of Mechanics ，Peking University ， Beijing 100871 China ) Abstract According to the principle that the computed response and measured response should be fitted ，the parameter inversion problem is reduced to a problem of solving nonlinear equations’ zero. Then the homotopy method can be used and it is widely convergent to solve nonlinear equations. At last ，this method is used to solve a 1-D model in porous media which has an analytical solution given by Simon in 1984. The numerical results confirm that homotopy algorithms are feasible and robust. Key words numerical simulation ，homotopy method ，porous media ，inversion，widely convergent 为非线性方程组的零点求解问题。传统的求解非线 1 引 言 性方程组的 Newton 法及其变形算法，存在着局部 收敛的问题，即只有当初始值与真解比