一种基于面向对象思想正演反演统一框架设计.pdfVIP

·490· 中国地球物理2011 (15)信息技术与地球物理 一种基于面向对象思想的正演反演统一框架设计 A frameworkof andinversion program geophysicalModeling baseon method object-oriented 张，志勇+李曼谢尚平 LiManXie ZhangZhiyong Shangping 东华理工大学 江西抚州344000 地球物理勘探是通过观测物理场的空间与时间分布来确定地质单元物性分布的探测技术。正问题和反 问题是地球物理研究的两个重要问题。本文运用面向对象的思想分析了正演与反演的层次结构，建立有限 单元法正演、正则化非线性最小二乘反演为基础的地球物理正反演统一框架。这一框架提高了代码的重用 性与开发速度；对地质单元采用统一的空间剖分，有利于开展联合反演程序开发。 1．统一框架结构分析 有限单元由于具有剖分形式灵活、对地质体、地形的拟合能力好、计算精度高等特点，在地球物理正 演计算中被广泛采用；正则化非线性最小二乘法可以比较有效地避免了反问题的多解性和不适定性，因此 在地球物理反演方法中被广泛采用。所以我们考虑采用有限单元和正则化非线性最小二乘法设计统一的地 球物理正反演框架。 】)有限单元正演 有限单元基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元， 在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导 数的节点值与所选用的插值函数组成的表达式，借助于变分原理或加权余量法，将微分方程离散求解。 采用有限单元法进行地球物理正演计算要解决的问题：①正演问题等价的变分或加权余量方程；②计 算区域的剖分方法；③单元插值函数的选取：④线性稀疏系统的求解。对于不同的正演计算问题，除了问 题①基本方程不同外，其它三个问题都是要解决的共同问题。即使是对于不同的基本方程，最后离散到各 单元的形式也是有限几种，这正是设计统一框架的基础。 2)正则化梯度类反演 反演的实质是在给定正演算法下求取物性分布，使正演计算数据与实测数据在某种距离定义下达到最 小。正则化的梯度类最优化方法在地球反问题中广泛采用，如牛顿法及其变形、共轭梯度法等。反问题求 解过程中，有几个共同需要解决的问题：①正演计算：②Fnchet矩阵的求取：③正则化因子的选择。④反 演方法的选择与流程设计。其中前两个问题由具体问题决定：后两个问题具有很强的通用性。 2．类结构与继承关系 用于完成模型的离散、生成刚度矩阵，对各种系数矩阵中单元属性的部分采用纯虚函数设计；建立单元类 如线段单元、三角单元、四面体单元等，做为模型剖单元，存储剖分单元的拓普结构、单元的物理属性。 3．结论 通过分析地球物理正演与反演的特点，建立了以有限单元正演、正则化梯度法反演为基础的框架结构。该 结构将地球物理正演与反演计算中具有通用性的部分进行抽象，提高了代码的重用性，加快程序开发速度，同 时使研究工作重心移动到具体的问题上，另外对地质单元采用统一的空间剖分有利于开展联合反演程序开发。 参考文献 J andelectrical thefinite-element …1 CoggonH．Electromagnetic modelingby 【2】LiYDW inversionofinduced Oldenburg．3-D polarizationdata【J】．Geophysics，2000，65，6：193l~1945． M R reversionof 【3】LokeH，Barker