anisotropic diffusion方法分析及其在提高地震资料品质方面的应用word格式论文.docxVIP

structure.Inthispaper,wepresentanisotropicdiffusionpioneeredbyWeickertandusethisedge-preservingfilteringfortheseismicimagetoimprovethequalityofseismicdata.Basedonthebasicprinciples,thenumericaldiscretization,parameteradjustmentanditerationstoppingcriteriaarediscussed.Byselectingtheappropriatemethodandadjustingeveryparameter,wetrytoachieveoptimalfiltering.IninterpretationsoftwareGeoScope,wedevelopatoolaccordingthealgorithmofanisotropicdiffusionusingtheC++programming.Weusethistooltoprocessseismicdatainsomeprojects.Theresultsshowtheremarkableeffectsofthisfilter,anditisveryhelpfulforfaultinterpretationandsedimentaryinterpretation.Keywords:AnisotropicDiffusionStructureTensorDiffusionTensorPartialDifferentialEquation独创性声明本人声明所呈交的学位论文是本人在导师指导下进行的研究工作及取得的研究成果。据我所知，除了文中特别加以标注和致谢的地方外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含为获得成都理工大学或其他教育机构的学位或证书而使用过的材料。与我一同工作的人员对本研究所做的任何贡献均已在论文中作了明确的说明并表示谢意。学位论文作者签名：年月日学位论文版权使用授权书本学位论文作者完全了解成都理工大学有关保留、使用学位论文的规定，有权保留并向国家有关部门或机构送交论文的复印件和磁盘，允许论文被查阅和借阅。本人授权成都理工大学可以将学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存、汇编学位论文。（保密的学位论文在解密后适用本授权书）学位论文作者签名：学位论文作者导师签名：年月日第1章引言1.1选题依据地震资料的品质是构造解释和岩性解释的基础。常规提高地震资料品质的手段往往在改善信噪比的同时，会损失地震的横向分辨率，如模糊断点、地层接触关系以及地质体边界，而这些信息对构造解释、地层及沉积解释非常重要，它们都是发现和描述油藏的重要依据。于是，我们在滤波的过程中引入了边界保护，以便保留地震图像中重要的地质特征。而要保护边界，就要先检测到边界的位置。所以，总体来说实现这种构造约束下的保边滤波（Structure-OrientedFilter-EdgePreserving）需要满足以下三个条件：1、方向分析：确定地震反射层的局部方向；2、边缘检测：确定可能的反射终止位置；3、保边平滑：沿着地震反射层的局部方向进行平滑，在检测到的反射终止处停止平滑。保边滤波的要点是做平行于地震反射层的平滑操作，在反射终止（断层）处不平滑。也就是说，一旦确定了倾角方位角，滤波就可以沿着反射轴进行，更重要是在所探测到的信号不连续处不再滤波，从而保护主要的断层和地层边缘。在传统的图像处理中，以上的每一步都可以用许多方法来实现，但它们都是相对独立的过程。近十多年来，随着基于偏微分方程的图像处理方法的迅速发展，非线性各向异性扩散算法不但满足以上条件，而且将这些过程有机地结合起来，与其它传统方法相比显示出了更加优越的性能。它从物理的扩散现象引申而来，将地震图像作为初始条件，通过解关于时间的偏微分方程得到扩散后的图像。在扩散方程中，通过引入“结构张量”提取图像局部的结构信息，根据这些结构信息设计“扩散张量”，在不同方向上选取不同的扩散策略，从而控制扩散过程，达到提高资料品质的同时保护地质体边缘信息的效果。1.2国内外发展与现状1.2.1 保边滤波作为图像处理的一部分，保边滤波从90年代才开始见诸文献，它很明显地反映在各向异性扩散的方法中，而在各向异性扩散方法提出以前，已经有一些方法实现了保边滤波。有一种算法从三维地震数据中提取二维地震数据，以及局部结构。然后对这个二维地震数据进行保边滤波，最后把结果再写回到三维数据体中。更好效果的边缘保护由Kuwahara-type方法实现（Kuwahara等人，1976；Nagao 和M