基于分段常值水平集的抛物型方程的参数识别算法.pdfVIP

第31卷 第4期 数学理论与应用 Vo1．31No．4 2011年 12月 MATHEMAnCALTHE0RY AND APPUCATIONS DeC．2011 基于分段常值水平集的抛物型方程的参数识别算法 纪双西 李维国 同登科 (中国石油大学(华东)理学院，青岛，266555) 摘 要 研究了线性抛物型方程不连续参数的识别算法．根据原有算法对于加噪观测数据计算不收敛的问 题。本文基于分段常值水平集方法，根据水平集函数和优化过程的特点，修正原有Uzawa型算法中的带有总变 差(TV)正则化的极小化模型和对常值向量的极小化模型，并且利用分裂Bregman迭代算法处理TV范数的优 越性，构造一种新的参数识别算法格式．数值实验结果显示，新算法具有计算时间短、精度高、抗噪性强的优 点． 关键词 参数识别 增广Lagrange方法 分段常值水平集 分裂Bregman迭代 A New Algorithm SchemeforParameterIdentification BasedonPiecewiseConstantLevelSetM ethod JiShuangxi LiWeiguo TongDengke (SchoolofScience，ChinaUniversityofPe~oleum (EastChina)，Qingdao，266555，China) Abstract Apiecewiseconstantlevelsetmethodwasappliedtoparabolicinverseproblemwithdiscontinuousparame— ter．Accordingtothefeatureoflevelsetfunction andoptimization process，combinedwith splitBregmanmethod whichhashteadvantageofhandingthetotalvariationBonn，theexistingminimizationfunctipnsinUzawaalgorithm andoriginalalgorithm schemeWaSmodified．Andthenanew algorithm schemewasgiven．Numerical experiments showthatthenew algorihtm hasadvantagesofshortcomputingtime，highprecisionnadstrongperformnaceofna ti— noise． Keywords Praameteridentification AugmentedLagrnagianmethods PCLSM SplitBremganmethod 0 引言 参数识别问题作为反问题的一个重要方面，具有理论上和实践上的双重价值，长期以来其 研究均受到国内外众多学者的关注，如文献[1—6]所示；但参数重构却具有理论上和实际计 算上的双重困难，如何更高效、更精确的获得参数重构结果，显然具有重要的实际意义． 国家 自然科学基金40874044)；中央高校基本科研业务费专项资金(09CX04004A)；中国石油大学研究生 创新基金 (，CXYB11—16) 收稿 日期：2011年 10月27日 基于分段常值水平集的抛物型方程的参数识别算法 l5 本文考虑如下的抛物型方程，其物理意义可以看作是一个非稳态的扩散过程，工程应用意 义为油藏一维单相不可压缩流的渗流方程或地下水渗流方程的简化数学模型： 一 · (g() )=厂 ∈力，￡∈(o，) M(，0)=12,0() in ， (，t)=0 ∈an 其中，q