层析成像读书报告.docxVIP

层析成像读书报告 姓名： 学号： 专业： 2013.11.11 层析成像是在物体外部发射物理信号，接受穿过物体且携带物体内部信息的信号，利用计算机图像重建方法，重修物体内部二维或三维清晰图像。 该技术最大的特点是在不损坏物体的条件下，探知物体内部结构的几何形态与物理参数的分布。 目前，层析成像有X射线层析、核磁共振（NMR）层析、正电子发射断层成像、微波层析、超声波层析、地震层析等到多个技术领域，而且还在继续发展中。 层析成像首先在诊断医学中得到应用，然后逐步推广到非医学领域中，诸如石油地震勘探、油气开发、地球物理学、射电天文学、超声学、微波、雷达、结晶学、电子显微技术等等。本课程针对我们地球物理系，主要讨论层析成像技术在井间地震的应用。 一、层析的物理原理： 地震层析成像(Seismic Tomography)，是利用地震观测资料信息来反演地质内部三维结构的一种方法。 对在地质中传播的地震波的走时和波形进行接收并分析计算，提出地质信息，借助计算机图象重建技术，重现地质内部二维或者三维清晰图象，进而揭示岩层性质、矿藏、地质构造等信息。 层析成像(Tomography)意思是把一个成像物体假想成一系列薄片，并显示出每个薄片上的物体图片。这种图像，医学CT成像中显示的是X射线穿透人体内脏的强度；在地震层析成像中，表示的是地震波穿透地下介质的慢度，也即速度的倒数。 层析的数学理论基础： 各种层析的数学理论基础是J.Radon的论文“论如何根据某些流形上的积分来确定函数”，现在人们称之为Radon变换或经典Radon变换，在经典Radon变换的基础上还发展了一种广义的Radon变换。 广义Radon变换是经典Radon变换的延伸，它由函数沿某种曲线簇的积分值，重建该函数或者恢复该函数的某种特性，以反演出地震地质结果。 三、层析成像的基本算法： 投影重建图像的基本算法，由于观测方式的不同、所应用的物理领域不同，对图像的要求不同，因而重建的算法也不同，但是它们都有共同的理论基础：Radon变换。 图像重建算法，就射线理论而言，主要有两大类：变换法和级数展开??，也有学者分成分析法和代数重建方法。分析方法中主要是Radon变换方法、Fourier方法、滤波反射投影方法或称褶积反射投影方法；在代数重建中，主要是ART方法、SIRT方法、LSQR方法；此外，还有级数展开方法。 四、层析的实现： 地震层析成像是典型的地球物理反演问题，大多数地震层析成像问题都涉及到以下几个方面： 型参数化； 正演（射线追踪）； 反演； 解的评价。 成像的目的就是要获得接近实际地下结构的模型，所以在成像前必需要建立模型来描述地层结构。 模型参数化可分为两类，一类是“不分块”的参数化，不对模型进行离散化，反演完全在泛函空间中进行，只是在最后计算想要的截面时采取离散化。由于反演在泛函空间中进行，理论上可以计算空间任何位置上的速度，结果不受离散化的影响，有利于成像的显示。 另一类是离散化的模型参数化，其优点是数学上容易处理，运算相对简单；缺点是在一般方法中出现的某些简化，在用离散时可能被掩盖掉。现在通用的大都是离散化的模型参数化。 一旦模型参数化的方式确定了，接下来要解决的问题就是怎么计算理论走时，这个过程就是正演。地震层析成像正演一般可以分为两种，一种是基于射线理论的射线追踪法，另一种是对波场进行数值模拟方法。 射线追踪方法是一种快速有效的波场近似计算方法，计算快速高效，所提供的射线路径非常直观。传统的射线方法可分为打靶法、弯曲法、及近似算法。 当构造比较复杂时，或者由于间断面的存在，可能会出现无法追踪到准确射线路径的情况，而全波场就不存在这样的问题。地震波场的正演数值计算是一个求解变系数偏微分方程的过程。目前，解偏微分方程比较常用的有：有限差分法、有限元法、边界元法、伪谱法、积分法等。 通过正演计算，我们可以得到关于震源和介质参数的观测方程式，反演的目的就是通过求解方程式来获得模型参数。地震层析成像反演方法可以分为两类：第一类是基于算子的线性或拟线性反演方法；另一类是基于模型的完全非线性反演方法，又称为“随机反演方法”。 地球内部成像结果反应的不仅仅是真实速度结构非均匀性，而且还有数据误差、有限的地震射线采样、模型参数化、线性化以及实施算法等因素带来的影响，因而经常导致最终图像的虚假异常。反演问题的解必须计算相应的分辨核和协方差矩阵，如果不计算这些参数来正确评估反演解的可靠性，这样反演的解几乎是没有任何价值。 五、井间地震层析??像发展展望： 如何提高图像分辨率和成像效果，仍是今后提高井间地震析成像技术实用化程度所研究的重点引起图像分辨率偏低、成像效果不好的主要原因是：有限观测角；介质的非均匀性和各