超声波层析成像技术中图像重建方法的对比研究.pdfVIP

西部探矿工程 2011年第 1期 超声波层析成像技术中图像重建方法的对比研究 顾兆君 ，于长青 ，李成刚 ，王 勇 (1．中国石油天然气华东勘察设计研究院岩土工程公司，山东 青岛266071； 2．吉林大学建设工程学院，吉林长春 130026) 摘 要：地震层析成像是一种调查和认识地下地质结构的有效方法，被广泛应用于工程与地质灾害和 地球物理勘探等领域。它包括模型的参数化、正演计算、反演及图像重建、反演结果的评价四部分，其 中反演及图像重建是地震层析技术的关键部分。因此，对图像重建算法的研究就尤为重要。在建立 数值模型的基础上，通过对ART和 SIRT数字化 图像重建效果的分析与比较，得出SIRT方法重建 图像效果最好。 关键词：层析成像 ；观测系统；图像重建；超声波 中图分类号：P631．8 文献标识码 ：A 文章编号：10O4—5716(2O11)01一O034一O3 地震层析成像技术在地球物理领域的应用非常广 约束求△ 的L。模的极小解。由拉格朗日乘子法令 泛，它是用地震数据来反演地下结构的物质属性，并逐 目标函数为： 层剖析绘制其图像的技术，其主要目的是确定地球内部 Q一 (△ ～ )+尬 ri min (2) 的精细结构和局部的不均匀性。在许多科学技术、工程 其中 为拉格朗日乘子。由aQ／a(aj~)：O得： 等领域中，由于投影重建图像问题频频出现，多年来人 们做了大量的工作去研究和寻找各种算法，通过数字化 △ 一鲁 (3) 图像重建效果的分析与 比较，重建一幅数字化的图像， 代入 (1)式有 ： 假如它是根据测试模型的模拟投影数据重建的，那么， 将它与数字化测试模型进行比较，就可以判断重建图像 At-i 鲁口； (4) 的质量的好坏，以此对重建算法的优劣做出评价和改 即得 ： 进。 1 基本原理 舍：==△／ ； (5) 1．1 代数重建技术(ART) 因此再由(3)式便可写出对第 i条射线及第 个象 ART(AlgebraicReconstructionTechnique)算法 元求图像 向量修改增量的公式： 是基于物理上的考虑而发展起来的。它的基本思想是 △ 一 Av／∑a； (6) 先给被重建的区域一个初始值，然后将所得到的投影值 由于图像重建是不适定问题，为增加计算的稳定 残差一个一个沿其射线方向均匀地反投影回去，同时不 断地对重建图像进行校正，直到满足所需要求为止 。 性，记 为第q轮迭代时用第 i条射线对第 个象元 AR T算法在每轮迭代过程中，将每一条射线的走 的图像向量估计算，而第i条射线的理论投影数据为： 时误差分配给射线经过的所有单元，略去r=Afq-e中 一 (7) 的误差项e，令波慢修改增量为af，有 AI-i~_∑△ 口，i=1，2，…，j (1) 对下一条射线(+1)，如果它通过象元 ，则由(6) J 式中对射线经过的象元 求和。作为迭代算法