利用地震层析成象技术确定横向非均匀介质中的速度分布和界面位置

本文提出一种横向非均匀介质中的地震层析成象方法，由地震走时资料确定介质的速度分布和界面位置。计算中采用了参数分离、变阻尼、奇异值分解等技术，方程的解较稳定，且可得到解的分办。除反射波外，该方法还可同时利用分析波的资料。因此，适宜于折射及宽角度反射地震测深剖面的资料处理。数字试验结果表明，在模型内部，该方法能得到可靠的反演结果。在边界附近效果较差，受初始值影响较大。     

共接收点倾斜叠加波动方程偏移

共接收点倾斜叠加波动方程偏移,本质上是一种叠前偏移方法.每给定一个斜率P,对经过叠前（动校正前）常规处理的地震记录中的各共接收点道集,沿直线t=τ+px进行倾斜叠加,就形成一个共接收点倾斜叠加剖面.对之进行波动方程偏移,该偏移剖面将代表地下真实构造.对一系列的p,我们可以得到一系列这样的偏移剖面.对它们作共接收点叠加,偏移叠加剖面的信噪比将超过水平叠加剖面.本文导出了在均匀、水平层状及非均匀介质条件下的共接收点倾斜叠加波动方程偏移算法.     

琼东南盆地深水陆架坡折区地震速度建模与反演的研究与应用

无论是构造简单或者复杂的地方,高精度速度建模与反演一直是地震处理所追求的目标,为进一步提升复杂构造区域的速度建模精度,提升该区域的成像质量.文中以琼东南盆地陆架坡折区速度建模与反演为例,采用构造形态对速度反演过程进行约束,通过模型和实际数据验证了混合邻域构造约束建模与反演技术的有效性,解决了坡折大角度构造常规速度建模及反演技术存在速度异常的问题.结果表明:相对于常规速度反演方法,本文方法得到与地质更加吻合的速度场,能够更好地解决成像及岩性认识问题.结论认为:本方法在速度建模的过程中实现了自动提取构造张量并约束建模过程,提高陆坡下方速度精度,改善了坡折下方成像,该方法对于崎岖海底地区速度建模具有较大应用价值.     

腾冲火山及邻区速度结构的三维层析成象

利用１９８４～１９９８年新间期间发生在腾冲火山及邻区３０个台站记录的数知个地震的Ｐ波资料,采用地震层析成象法重建了腾冲火山及邻区地壳、上地幔不同深度的三维速度结构层析成象图。结果给出,在腾冲火山区下面地壳内约３～９ｋｍ存在低速带,扰动量达－２０％;１０～１５ｋｍ深度为高波速带,扰动量达７．３％;１６～２４ｋｍ深度为低波速带。扰动量达－９．１％;２５～４０ｋｍ为波速带,扰一达１０％;大于莫霍面深度（     

MSH法直接将地震勘控测线转化为深度、密度剖面

根据运动学中的动量守恒定律、分析地震勘探过程中地震波在地下介质中的传播过程,利用反射能量、透射能量与入射能量的计算公式和检波器转化的电信号与地下密度界面反射地震波波动能量的等效原理,推导出直接求层速度和层密度公式.利用某一反射界面直接测得的波速或依据钻孔资料求得的平均速度,在计算机上就可直接计算出各地震测线物性层的埋深和密度(速度、波阻抗),由此地质工作人员依据已知地质柱状剖面(钻孔资料)可进一步将深度、密度(速度、波阻抗)剖面解释成深度岩性地质剖面,为解释油气层、计算油气储量、各种地质构造提供更为可靠的物探依据.     

青藏高原地区地壳升降运动的分析

通过青藏高原及其邻区地壳、上地幔三维S波速度的分布, 换算出不同深度的压强分布规律;结合岩石类型, 地震及地壳相对运动进行综合分析, 探讨青藏高原地区地壳上地幔结构, 并针对人们关注的问题(青藏高原隆起的驱动力是什么？它还会上升多高？它会不会突然下落？)做了讨论。      

中国西部及邻区的地壳三维速度结构

用网格频散反演技术计算了155条瑞利面波混合频散数据,将中国大陆分为目前最小分格2°×2°,得到了中国西部及邻区深至80km地壳和上地幔顶部的三维剪切波SV速度结构。结果表明:(1)天山褶皱系地壳平均厚度为54km,上地壳平均速度为3.3km/s,下壳平均速度3.93km/s,上地幔盖层平均速度为4.64km/s,都偏大。北天山壳厚比天山薄,壳和上地幔盖层下均有低速层。(2)塔里木盆地中部地壳厚46km,到其边缘的昆仑山、天山地带壳厚增加到52～55km。壳内无低速层,壳平均速度很低,为3.49km/s,地表有9km的沉积层。(3)青藏高原莫霍面呈凹形,平均壳厚为67km,平均壳速度为3.58km/s。东西端存在壳内低速层,埋深为15～35km。(4)南北带为构造复杂的区域,莫霍面深度具有东高西低,南北两端高、中间低的特点,带内普遍存在壳内低速层。南北带从群速度、速度、莫霍面等深线都可看到明显的横向差异。(5)准噶尔盆地地壳厚47km,无壳内低速层,中地壳不明显。该区呈稳定的大陆壳结构。(6)柴达木盆地在正常速度的中壳顶部有一高速夹层,速度为3.72km/s,厚10km,上顶部埋深14km。     

含油岩石弹性特征及其与油气的关系

本文用准噶尔盆地东部油区58口井共2257块油气储层岩芯标本,在常温常压和模拟地下温度压力下对V_P/V_S及衰减系数进行了精确测量,结合密度、孔隙度及岩石学研究获得个同岩性地层含油、水、气不同流体的V_P/V_S比、泊松比及不同深度的速度校正曲线;用回归分析获得V_P-V_S,密度-速度,密度-孔隙度的相关关系式;用数学统计获得稠油地区岩芯砂岩的速度特征值;用AVO技术进行模型计算和实例验证.为油气碳氢检测,测井速度标定,声阻抗及演,模型正演等提供了有价值的数据.     

隧道地震反射法超前预报

利用隧道垂直地震剖面所得到的负视速度同相轴,用此作为识别来自隧道施工掌子面前方界面的反射波标志．然后利用这一标志推演出反射界面的空间位置与产状,达到超前预报的目的．     

一种适于近场理论地震图宽频域计算的矩阵分解算法

理论地震图方法在地震震源过程的研究中得到了广泛的应用。为了研究震源过程的细节,必须利用近场地震资料的高频信息,但是在应用Haskell矩阵法计算近场理论地震图时,格林函数的高频成分的数值不稳定性是计算宽频域理论地震图的一个基本困难。因此,目前仅有不超过10Hz的算例。本研究采用Haskell矩阵的一种新的分解组合形式,在数值计算时有效地避免了在计算过程中数值结果的溢出,实现了近场理论地震图的宽频域计算。同时,由于利用了矩阵运算中的解析关系,减少了运算的次数,从而提高了计算的速度。根据本算法建立的Fortran程序,在Univoc-1100计算机上进行了数值检验。结果表明,计算得到的格林函数至少在0—40Hz频率域范围内仍保持良好的数值稳定性。本文给出的这种Haskell矩阵分解组合形式的矩阵元素简单,且具有一定的对称性。对于这种分解所包含的物理意义,将有待进一步深入研究。