层析成像数据和射线不均匀覆盖问题解决方法.pdfVIP

中国地球物理201I ·565· (17)油气田与煤田地球物理勘探 层析成像数据和射线不均匀覆盖问题的解决方法 tohandlenon-uniformdataand Approaches raycoverage in tomography problem 黄光南’刘洋 Liu HuangGuangnanYang 中国石油大学(北京)CNPC物探重点实验室北京102249 1．层析成像的数据不均匀覆盖问题 层析成像最早应用于医学，地球物理领域内的层析成像主要研究如何利用地震波、电磁波或其它场的 数据对地球内部结构进行成像…。层析成像在地震领域的应用受到诸多因素的限制，分辨率不足是主要的 限制因素之一，它是受观测系统的限制所导致的。这种由于观测系统的限制所引起的数据和射线不均匀覆 盖决定了地震层析方程组是一个混定问题，在这种情况下，合理加入各种正则化方法是非常必要的，在某 种程度上层析解严重依赖于正则化过程中加进先验信息和对解的某种预期估计。 2．数据和射线不均匀覆盖问题的原因 在实际射线层析问题中，数据和射线不均匀覆盖问题主要来源于两个方面：一个方面是与观测系统有 关，包括不均匀分布的炮点间隔、不均匀分布的检波点间隔、不同炮点对应有不同的偏移距、以及不同地 震道具有不同的数据质量等，这些因素均会对射线的不均匀覆盖产生影响，速度结构也影响射线的覆盖方 式，由于速度的非均匀性射线可能在某些地方聚集而在其它某些地方发散：另一方面原因是内在的，对一 炮而言，所有的射线都从炮点出发，炮点处射线是非常密集的，而远离炮点的地方射线则相对稀疏，特别 是模型的边界处射线的覆盖更为稀疏甚至根本没有覆盖。如果忽略不均匀覆盖问题，层析过程中模型校正 量会正比于射线的覆盖程度【2】：射线覆盖程度高的网格被过分校正，射线覆盖程度低的网格校正量不足， 而没有被射线穿过的网格校正量为零，这样会导致所得层析解不能较好地反映真实的速度分布情况。 3．数据和射线不均匀覆盖问题的解决方法 自回归外推方法利用观测数据对射线覆盖区域进行外推，扩大射线覆盖区域范围，从而得到更能反映真 实速度分布的层析模型【3】。自回归外推方法包括两个阶段：第一个阶段是寻找最佳预测误差滤波器；第二阶 段是用预测误差滤波器寻找缺失的数据值。对于变网格间距方法，如果扩大反演网格的间距，网格所含的射 线数将会增加。因此，我们可以对射线覆盖程度低的区域的网格间距进行扩大来加强对数据的控制，换句话 说，变网格间距方法是对数据信息的重新调整，使每个模型网格的覆盖程度更加均匀。变阻尼方法是根据射 线对网格的覆盖程度对不同网格参数加以不同的权重，从而达到缓解校正量正比于网格覆盖程度的目的，在 反演的正则化过程中单独使用变阻尼方法会带来不稳定性，为了降低模型参数的不确定性，通常在正则化过 程中同时加入平滑约束条件，层析过程中同时考虑变阻尼方法和平滑约束条件，我们可以得到如下表达式： 7’G1A+A(D7’D+￡G))一ArG’ 《1一(A 式中，D为平滑算子，G为变阻尼算子，A为平衡数据拟合度与约束条件之间的参数，F为平衡平滑和阻 尼约束的相对权重参数。 模型数据试算验证了上述方法的正确有效性。 本研究由国家科技重大专项课题(2008ZX05024．001)资助。 参考文献 【1】成谷．地震反射走时层析理论与应用研究【D】．同济大学海洋与地球科学学院博士学位论文，2004． L to data W．Effective non·unifor