多震相的层相成像与在介质中的多波型射线跟踪

本文介绍了多震相的层析成像的思路和算法，通过穿透和反射走时可以同时作出２维和３维慢度（速度的倒数）重建。我们分析了在穿透和反射数据中确定速度和深度的不确定性，并认识到深度扰动对反射走时异常比慢度扰动更敏感。由不同波类型所提供的对速度和深度的约束，这个算法实际上减少了在一般反射层析成像在速度和深度之间的不确定性，并且也避免了在穿透层析成像中的不确定问题。线性化反演是通过从反射界面深度由分离速度参数进     

基于图形结构的三维射线追踪方法

在地震层析成像研究中,为了克服最小走时射线路径追踪方法存在的问题,对该方法计算过程中的关键步骤进行了改进．在节点走时的计算中引入Ｂｒｅｓｅｎｈａｍ画线算法;在最小走时节点查寻中,结合使用快速排序算法与插入排序算法,替代以往方法中多采用的堆排序算法;所采用的节点设置方式,可以引入速度界面,还可以实现反射波射线追踪．模型计算证明,改进的最小走时射线路径方法具有精度高,速度快的特点,所提出的三维空间反射波射线追踪算法简便易行。     

多震相菲涅耳体射线走时同时反演成像

高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线附近一定范围(称其为第一菲涅耳体)内的速度异常分布有关.鉴于此,本文提出了计算多震相地震波菲涅耳体有限频射线的方法,并定义了走时敏感核函数,同时给出了利用多震相菲涅耳体有限频射线进行速度模型和反射界面同时反演成像的公式.利用多震相走时资料,使用传统射线层析成像方法与有限频射线层析成像方法进行了速度和界面的同时反演成像.结果表明,当射线密度较小时,无论是对速度模型的重建还是对反射界面几何形状的更新,有限频射线层析成像方法均优于传统射线层析成像方法, 而变频有限频射线层析成像则是实际地震层析成像的首选反演算法.      

反射地震层析成像在海洋天然气水合物储层识别中的应用

海底地层速度结构是识别海洋天然气水合物储层的直接依据,本文应用地震反射走时层析成像建立了海底地层速度模型.采用不规则网格对模型进行离散化,使速度单元与反射界面单元完全耦合;利用基于不规则单元波前扩展和走时插值的射线追踪方法,精确确定反射射线路径和反射波走时;在反演中同时使用先验约束、平滑约束、归一化和正则化技术,提高了层析反演的稳定性和结果的可靠性.对南海北部神狐海域SH2井附近的二维地震测线资料,利用多域人机交互法拾取了反射走时,用反射走时层析成像方法获得了SH2井附近含水合物目标区的速度结构.该反演结果与测井声波速度和钻探结果一致,其中的高速带对应水合物储层,表明反射走时层析成像能够有效地得到海洋天然气水合物储层的速度结构,为海洋天然气水合物储层识别提供依据.     

延庆─怀来地区地壳细结构──利用深地震反射剖面

１９９２年底在延庆－怀来地区实施了深地震反射剖面测量，目的是通过对该区地壳细结构的研究来探讨延庆－怀来盆地形成演化的机制和过程，以及与该区地震活动图像的关系．结果表明，研究区上部地壳总体上呈“透明”性质，厚２５－２６ｋｍ；下部地壳具有强烈的反射性质，厚约１４－１５ｋｍ．壳幔过渡带是由双程走时约２ｓ和厚约５－６ｋｍ的水平反射叠层所组成，莫霍界面具有深度深（约３９－４０ｋｍ）和平坦的特点。     

延庆─怀来地区地壳细结构──利用深地震反射剖面

１９９２年底在延庆－怀来地区实施了深地震反射剖面测量，目的是通过对该区地壳细结构的研究来探讨延庆－怀来盆地形成演化的机制和过程，以及与该区地震活动图像的关系．结果表明，研究区上部地壳总体上呈“透明”性质，厚２５－２６ｋｍ；下部地壳具有强烈的反射性质，厚约１４－１５ｋｍ．壳幔过渡带是由双程走时约２ｓ和厚约５－６ｋｍ的水平反射叠层所组成，莫霍界面具有深度深（约３９－４０ｋｍ）和平坦的特点。     

计算最小走时和射线路径的界面网全局方法

用慢度分块均匀正方形模型将介质参数化，仅在正方形单元的边界上设置计算结点，这些结点构成界面网．根据Ｈｕｖｓｅｎｓ和Ｆｅｒｍａｔ原理，由不断扩张、收缩的波前点扫描代替波前面搜索，在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似，通过比较确定最小走时和相应的次级源位置，记录在以界面网点位置为指针的３个一维数组中．借助这些数组通过向源搜索可计算任意点（包括界面网以外的点）上的全局最小走时和射线路径．这一方法不受介质慢度差异大小限制，占内存少，计算速度较快，适于走时反演和以Ｍａｓｌｏｖ射线理论为基础的波场计算．     

利用井间地震层析成像方法探测隐伏断层

针对目前隐伏断层探测中存在的技术难题,提出了在地面反射地震勘探初步确定的断点附近,用井间地震层析成像方法探测隐伏断层及其上断点的埋深。首先介绍井间地震技术的应用现状和井间地震层析成像的基本原理,并结合实例讨论了井间初至走时层析成像的野外数据采集和室内资料处理方法。芦花台隐伏断层的井间地震层析成像结果表明,采用弯曲射线追踪方法反演得到的跨芦花台隐伏断层的井间纵波速度结构图像与钻孔联合剖面揭示的断层特征相吻合。      

近震资料确定三维地震速度结构的一种方法

我们研究了由近震得到的地震走时资料估计地震波速度三维结构的反演方法。这种方法由同时确定速度界面的二维深度分布、三维速度分布、台站校正和震源来表征。速度界面的深度分布和慢度扰动的分布以纬度、经度和深度的幂级数来拟合。这种三维地震结构的表示不仅有助于减少未知参数的个数,也有助于解析地求取界面深度、速度和走时关于未知参数偏导数的值。连接震源和台站之间的射线轨迹,是在具有确切界面深度分布的简化的速度分布条件下由称为试射法的射线追踪方法确定的。数值试验证明,我们的反演方法可成功地估计三维速度结构。     

地震射线辛几何算法初探

走时计算已广泛用于地震建模、成像及速度分析等诸方面．基于地震波场的Ｈａｍｉｌｔｏｎ力学性质,本文探索应用适于Ｈａｍｉｌｔｏｎ力学的计算方法──辛几何算法对地震射线 进行走时及路径的计算．利用辛几何算法和属于耗散算法的四阶 Ｒｕｎｇｅ－Ｋｕｔｔａ算法进行射线 路径和走时的计算对比,结果表明这两类算法数值精度相当,但辛几何算法的速度却快了３ 倍．本文还利用二阶Ｅｕｌｅｒ型辛差分格式对Ｍａｒｍｏｕｓｉ模型进行了射线追踪计算,结果显示所 得射线具有良好的光滑性和较强的阴影区穿透能力．     

地震立体层析成像的实现方法及效果分析

立体层析成像是一种新的地震反射波层析成像方法,能为叠前深度偏移提供较为精确的宏观速度模型。本文研究了立体层析成像的实现方法,包括斜率与走时数据的拾取、离散速度模型构建和初始化、射线参数的确定、斜率和走时及射线计算以及反演问题解法等,建立立体层析成像的算法流程。并通过对Marmousi模型试验,对立体层析成像运行所需的主要参数,如初始速度模型、拾取数据量、离散网格尺寸、速度平滑权重等进行测试和分析,总结这些不同参数对立体层析反演结果的影响规律,用以指导生产实践。     

地震层析成像反演中3种广义解的误差分析与评价

在地震波速度和地震波　值层析成像中,对于非线性反演问题解的可靠性研究非常重要．为了讨论ＡＲＴ（ＡｌｇｅｂｒａｉｃＲｅｃｏｎｓｔｒｕｃｔｉｏｎ　Ｔｅｃｈｎｉｑｕｅ）、ＳＶＤ（Ｓｉｎｇｕｌａｒ　ＶａｌｕｅＤｅｃｏｍｐｏｓｉｔｉｏｎ）、ＬＳＱＷ（ＬｅａｓｔＳｕａｒｅ）３种方法的算法稳定性以及反演效果,本研究将这３种算法用于求解非常病态问题和一般稳定问题,根据数值计算结果,分别给出了在非常病态问题和一般稳定问题中３种算法的相对误差及其算法稳定性评价．进而针对３种算法在同一地质模型上的速度成像结果进行了比较,为地震层析成像反演效果分析提供了定量参数．     

地震层析成像反演中3种广义解的误差分析与评价

在地震波速度和地震波　值层析成像中,对于非线性反演问题解的可靠性研究非常重要．为了讨论ＡＲＴ（ＡｌｇｅｂｒａｉｃＲｅｃｏｎｓｔｒｕｃｔｉｏｎ　Ｔｅｃｈｎｉｑｕｅ）、ＳＶＤ（Ｓｉｎｇｕｌａｒ　ＶａｌｕｅＤｅｃｏｍｐｏｓｉｔｉｏｎ）、ＬＳＱＷ（ＬｅａｓｔＳｕａｒｅ）３种方法的算法稳定性以及反演效果,本研究将这３种算法用于求解非常病态问题和一般稳定问题,根据数值计算结果,分别给出了在非常病态问题和一般稳定问题中３种算法的相对误差及其算法稳定性评价．进而针对３种算法在同一地质模型上的速度成像结果进行了比较,为地震层析成像反演效果分析提供了定量参数．     

天然地震折射与反射速度CT的模型实验

为提高天然地震折射与反射速度CT的精度,我们运用＂走时最小纵横波速度层析成像＂与＂波形基本解纵横波速度层析成像＂相结合的方法,反演各层界面形状及各层纵波、横波速度。反演结果图与＂正演地球模型剖面图＂基本一致。本文在数学上发展了＂层析成像＂的内容,提出＂波形基本解纵横波速度层析成像＂基本概念,是天然地震简单模型实验的初步...     

多震相走时联合三参数同时反演成像

采用新近研制的分区多步不规则最短路径多震相地震射线追踪正演技术,结合流行的子空间反演算法,提出了一种联合多震相走时资料进行地震三参数 （速度、反射界面和震源位置） 同时反演的方法技术.数值模拟反演实例、以及与双参数 （速度和反射界面或速度和震源位置） 同时反演的对比分析表明：三参数同时反演成像结果大体接近双参数同时反演成像的结果.另外,噪声敏感性试验表明：所提算法对到时数据中可容许的随机误差并不敏感,结果说明多震相走时的联合三参数同时反演成像方法技术不失为一种提高走时成像空间分辨率、进而降低重建模型参数失真度、行之有效的方法技术.     

垂向非均匀介质的走时反演

本文提出了一种利用反射波走时曲线计算垂向非均匀介质速度和反射界面深度的方法。当在地球表面没有获得来自地下某一深度范围内介质的任何信息时,可以认为这一深度范围内的地震波速度具有连续性。利用来自其底部反射界面的反射波走时曲线,可以计算出该深度范围的地震波速度结构。对三种模型进行了理论计算,所得反演解与真实值较为一致,其中计算出的反射界面深度最为精确。利用本文提出的方法可以计算两个相邻反射面之间的垂向非均匀速度结构,如低速层等。     

界面和速度反射联合成像──理论与方法

依据波逆行原理重新推导了由地震波反射走时资料反演界面深度的关系式──走时对界面偏导数关系，它不仅适用于反射波，也适用于透射波．还给出了任意多个复杂界面情况下走时对界面偏导数关系的离散形式．另外，正交算子投影法被推广用来解决速度与界面的同时成像问题．数值模拟的结果表明，本文的方法是有效的．     

长白山天池火山区三维地壳结构层析成像

利用长白山天池火山区三维空间深地震测观测系统所采集的反射P波走时资料,采用层析成像技术,重建了该区地壳界面构造形态和速度分析图像。界面成像结果表明：研究区地壳界面总体上显示了由北西向东南加深的趋势;马鞍山—三道白河断裂和富尔河—红旗河断裂是本区两条主要的深部构造,尤其是马鞍山—三道白河断裂,北北东向穿越天池火山口,其两侧的地壳界面存在明显的错断,预示了该处地壳厚度陡变或深大断裂带的存在,速度成像结果显示在10km深度,明显的P波低速异常分布在天池周围;15km深度上它表现为一个近南北向的P波低速异常条带,其延展尺度南北向为80~90km,东西向30~40km;随着深度的增加,P波低速异常分布在天池西侧,其尺度有明显的缩小,分布范围更加集中,而且低速扰动幅度更大。这种P波速度异常的变化图像在一定程度上反映了天池火山口下方壳内岩浆系统的空间分状况。     

滇西地区地壳浅部基底速度细结构的研究

本文论述滇西地区思茅－中甸剖面基底导速度细结构的研究结果。基底层项部Ｐｇ界面深度为０－３．５ｋｍ，而底部Ｐ１＾０界面深度为１１．０－１７．０ｋｍ。基底层速度在金河－洱海断裂以南为５．７０－６．３０ｋｍ／ｓ，断裂以北增至６．３０－６．５０ｋｍ／ｓ，其过渡带位于剑川附近，剖面上断裂附近，除界面深度变化外，Ｐｇ面速度横向变化也较明显。速度等值线较为稀疏的景云桥炮南侧、大仓炮附近与支梯炮南侧地区，估计发生     

甘肃天水地区地壳上部二维速度结构成像研究

利用有限差分和不规则网格最小二乘正反演算法,对青藏高原东北缘暨西秦岭天水－武都地区两条高分辨人工地震测深剖面的Pg波走时进行层析成像反演,结果显示该区基底界面起伏较大,构造复杂,断裂丰富、较陡且延伸较大,多数断层切割基底.礼县-西和之间的低速异常反映了沉积盆地的结构特征,结晶基底深度在5 km左右.武山以东速度等值线横向变化明显,表现为褶皱.成县盆地速度变化平缓且速度值相对较高,基底埋深较浅,深度不到3 km.礼县附近存在一条近南北断裂,5 km深度以下有明显的上涌低速异常,可能为青藏高原深部幔源物质上升的通道.