采用射线寻迹的井间层析成像

本文采用二维井间观测走时反演地震波速,用迭代射线寻迹走时成像方法进行层析成像.作者导出三次多项式插值基函数计算离散网格内的折射率分布,还给出了具有自己特点的二步扫描法,不但加快了正演计算时间,同时方便地建立走时残差与模型修正量的关系.考虑到实际应用时,震源与接收器常不位于同一平面,为此进行了校正,并给出了校正方法,提高了重建精度. 数值的与实际应用的结果表明,本文给出的层析成像算法是快速的、稳定的、适用的.     

采用射线寻迹的井间层析成像

本文采用二维井间观测走时反演地震波速,用迭代射线寻迹走时成像方法进行层析成像.作者导出三次多项式插值基函数计算离散网格内的折射率分布,还给出了具有自己特点的二步扫描法,不但加快了正演计算时间,同时方便地建立走时残差与模型修正量的关系.考虑到实际应用时,震源与接收器常不位于同一平面,为此进行了校正,并给出了校正方法,提高了重建精度. 数值的与实际应用的结果表明,本文给出的层析成像算法是快速的、稳定的、适用的.     

三维多层介质重力-地震同步联合反演

联合反演是地球物理勘探的重要解释手段,能够提高模型参数的反演精度.本文在归纳和分析重力与地震资料联合反演的研究和应用现状的基础上,利用三维多层介质模型的地震走时和重力正演公式,推导了地震走时和重力异常对界面深度的雅可比矩阵,实现了三维重力-地震同步联合反演界面成像.最后进行了数值理论模型模拟和实例计算,结果表明地震走时和重力同步联合反演很好的重建了三维多层介质界面.     

基于Chebyshev多项式的非对称走时Kirchhoff叠前时间偏移角道集求取

地震射线走时的求取方法是叠前时间偏移研究的核心问题之一,也是影响计算时间域角道集角度精确性的关键因素之一.本文基于Kirchhoff叠前时间偏移,应用第一类切比雪夫多项式,对弯曲射线对称走时加以改进,引进非对称项,优化后得到切比雪夫非对称走时方程,与高精度走时进行比较和误差分析,再将该走时求取方法应用于时间域角道集的求取中,得到地下较真实的入射角.通过模型计算和实际地震资料处理证明,此种非对称走时及其角道集的求取方法具有精度高、计算量少的优点.     

加权抛物Radon变换叠前地震数据重建

基于部分动校正（NMO）后反射同相轴在CMP道集上的抛物线走时近似,给出了加权抛物Radon变换叠前地震数据重建方法（WPRT）. WPRT通过在迭代过程中引入变化着的权系数,拓展和改进了传统抛物Radon变换方法,使其可同时完成不规则采样的规则化和空道及近偏移距道重建,且有更高的计算效率. 文中给出了应用WPRT进行近偏移距和中偏移距的空地震道重建及数据规则化的算法实现. 理论模型和实际地震资料的地震数据重建结果显示了本文算法的优点.     

VSP直达波和反射波波动方程走时联合反演

地面资料全波形反演采用低波数回折波敏感核和高波数偏移等时线恢复地下模型的长波长和短波长分量.当回折波的穿透深度有限时,很容易陷入局部极值.波动方程反射波走时反演采用反射波的透射敏感核更新速度模型.当浅层速度存在较大误差时,仍无法得到满意的反演结果.VSP资料中直达波是一种透射波,穿透深度大,可用于井旁背景速度场建模.本文发展了针对VSP观测方式的波动方程走时联合反演(直达波和反射波)方法.采用全局优化参数反演法实现VSP上下行波和纵横波分离,构建了基于下行直达波和上行反射波走时残差的混合目标函数,推导了相应的伴随状态方程和梯度公式,给出了背景速度和反射系数分步反演流程.理论和实际VSP资料应用表明：直达波和反射波走时联合反演可以得到运动学特征精确的速度模型,为后续全波形反演提供可靠的初始模型;联合反演比单独直达波/反射波走时反演的精度高,比直达波和反射波走时级联反演的耗时少.     

二维复杂层状介质中地震多波走时联合反演成像

采用新近提出的多次波射线追踪正演算法,结合共轭梯度法求解带约束的阻尼最小二乘最优化反演问题,分析讨论了利用多震相走时资料进行联合反演成像的方法及技术.考虑到不同震相走时的拾取误差不同,反演算法中引入了不同震相种类数据的权系数; 由于同时反演速度模型和反射界面起伏中不同模型参数变化对走时影响程度的不同, Jacobi偏导矩阵元素中引入了不同参数的归一化因子; 另外,为了克服射线密度过大(或过小)区域速度模型的过度(或欠)更新问题,反演算法中引入了等权射线密度的概念.几种数值模拟实例表明(含噪声敏感性试验): 多波走时的联合或同时反演成像技术是一种提高走时成像空间分辨率,进而降低重建模型失真度行之有效的方法.     

琼西北地区的基底速度图象

利用人工地震测深纵剖面和非纵剖面的资料,对琼西北地区P_g波的走时进行了层析反演,对比了代数重建法和同时迭代重建法的效果.综合了纵剖面的模型后,在三维弯曲界面上重建了基底速度图象,探讨了用射线椭圆评价分辨的途径. 层析成象结果给出了许多有价值的信息,揭示琼西北的结晶基底存在三条北东向断裂和两条北西向断裂,后者是继承性的活动断裂,控制着新生代次级断块的发展与火山活动.地表所推断的东西向王五-文教深大断裂,在该处只呈断续隐伏状显示.该区的基底存在儋县、临高和光村三个高速区,主要反映着基性岩体或花岗侵入体.基底整体性最好的是在儋县-王五地区,破碎程度最严重的是在王五-中和-光村以西的地区.     

一种扇束扫描模式下的图像重建迭代算法

CT图像重建的扫描模式有平行束、扇束、锥束等,在扇束扫描模式下的图像重建算法大多基于图像的正方形网格剖分。本文建立了扇束扫描模式下新的图像重建离散化模型,并给出了基于新模型的代数迭代校正格式和重建算法。对新的模型下迭代算法几何意义进行了讨论,基于新模型的代数迭代重建算法有助于提高成像质量,启发新的图像重建算法。     

利用单程波算子进行叠前时间偏移走时计算(英文)

传统的叠前时间偏移公式都是基于水平层状模型的假设而推导而来,在速度横向变化剧烈的介质中则不能够达到理想的聚焦效果。本文基于李代数积分方程以及拟微分算子等理论推导了非对称走时公式,计算公式由于包含有速度横向导数项,在速度横向变化大的介质中,也能有较高的聚焦效果,而且也为走时计算提供了精确的相对振幅保持权系数。本文对推导的方法进行模型测试并进行实际数据的试算,其结果证明非对称走时方法的成像精度远高于对称走时计算方法。     

井间地震数据直达波走时层析成像

本文利用矢量射线追踪正演模拟技术计算地震波直达波传播的路径及走时,进而利用射线走时及路径的内插,发展了弯曲射线迭代反演技术.该方法可用来重建井间地层的速度图像.基于所发展的方法,我们对较为复杂的典型地质模型进行了井间速度重建.结果表明该方法是一项快速、高精度的走时层析成像技术.     

地震层析成象的-变换

由Radon变换出发给出了地震走时问题的——变换.定义模型坐标(x,z)与震源位置()和射线斜率()间的线性关系为x=+z.文中给出了-正反变换的理论公式,并分析了-变换与富氏变换的关系。理论表明,连续函数的-变换可得到走时反演问题的唯一解。 在地震学非完全层析成象中,欲提高图象质量必须压制假频和提高分辨.作者对重建图象的分辨进行了详细讨论,最后通过数值模拟讨论了观测系统和滤波处理的影响。      

适于Kirchhoff叠前深度偏移的地震走时李代数积分算法

本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对以往Kirchhoff叠前时间偏移中所用算法进行深化调整.文中数值算例对比了本文李代数积分表达与时间积分的区别,本算法计算结果与线性横向变速介质中的理论值相当吻合.通过走时多项式中各项对结果的影响分析,可知非对称项使计算精度得到了进一步提高.数值试验表明,本算法对横向变速介质中走时求取是可行的,且不需要存储海量走时表,有利于提高Kirchhof叠前深度偏移的精度和效率.     

物理模型超声CT图像代数法重建和结构识别

本文采用木板、有机玻璃板制作二维物理模型,并在模型中加入三个大小不一、位置各异、用石膏充填的方洞。对两个物理模型布置相同的观测系统,形成相同的射线方式,在等间距矩形像素条件下,每条射线的积分步长（射线穿过单个像素的射线段长度）皆为常数。由超声仪器沿指定射线采集各条透射波射线的走时数据（即Radon变换线积分的“观测值”）。应用“代数重建法”（ART加法修正迭代法）求解Radon变换的离散化方程组,重建二维物理模型的CT数字图像。在重建的模型CT数字图像中,能够有效地识别出“石膏洞”的存在和位置,验证了“代数重建法”CT成像技术及其对模型结构识别的有效性。本文是ART加法修正迭代的基础性应用工作,可为有关部门提供“代数重建法”CT成像实际应用参考。      

跨孔观测地震数据的速度重建

利用矢量射线追踪正演模拟技术 ,计算地震波传播的路径及走时 ,进而利用射线走时及路径的内插 ,发展了弯曲射线迭代反演技术 .该方法可用来重建井间地层的速度图像 .基于所发展的方法 ,我们对两种较为复杂的典型地质模型进行了井间速度重建 .结果表明该方法是一项快速、高精度的跨孔数据速度重建技术     

代数重建技术(ART)应用于跨孔地震资料——Ⅰ(方法)

本文将剖面成象技术应用于两组跨孔地震资料,来确定速度结构以及判明各种算法对买际资料的可靠性和分辨率。实验是在纽约Retsoff盐矿和瑞典Stripa矿地下放射性废料研究基地进行的。Retsoff的走时资料质量很高,并且在长度直到500m的整个射线路径上都能记到。Retsoff地区结构相当复杂,其速度差高达50％。Stripa场地位于速度差仅有百分之几的花岗岩上。实验范围不大,最大射线长度刚刚超过10m。该场地采集的资料精度很高,源和接收点位置的测量精度小于1.0mm走时的读数误差为0.001ms。本文将若干类似于医学成象的代数重建技术(ART)应用于这些资料。算法作了一些修改,例如,加权模式的应用,阻尼参数的应用,利用弯曲射线路径。所得到的速度场与已知场进行比较,相互对照以便确定一种最优方法。业已发现,这些算法是重建实际资料慢度场的一种速度可靠的方法,低速带的位置和速度值精确地被恢复原状。同时还发现,为保证利用适当的阻尼参数,取适当的叠代次数,应用这些技术时必须非常小心,否则重建技术的结果会很差。     

层析建模技术在库车坳陷东秋山地区静校正应用探索

地震走时层析建模技术是解决复杂近地表模型速度建模问题的重要技术.该方法是一种迭代反演方法,在地震反演过程中需反复计算地震射线走时.本文采用交互初至拾取方法,通过层析反演出近地表模型,对比分析了走时层析反演表层模型与野外表层模型差别,从而引入了控制点成果约束的层析反演综合建模技术.最终将野外调查成果同走时层析反演建模结合起来,运用非线性走时层析反演技术,通过可视图片、数据、曲线的比较,得出如下结论:在速度纵横向变化大的复杂地区,通过野外资料与室内处理信息交互约束,可以为建立符合实际情况的近地表模型奠定基础,进一步为地震采集资料静校正处理方法的确定和叠加剖面品质提高创造了条件.     

复杂介质下保真振幅Kirchhoff深度偏移

本文详细讨论地下复杂介质地震成像问题.当速度模型含有强横向变化时,射线场会折叠形成多次走时,在这种情况下,文中证明在保真振幅叠前深度偏移时,应当考虑所有的到时才能得到定量的像.通过实例证明,在多次走时情况下,常规的共炮集、共检波器集都存在强的假像;由此提出共角度偏移的概念.最后,在多次走时下,讨论多次偏移算子(走时、相移、振幅等)的高效重构,给出了三维叠前深度偏移的快速算法.     

地壳三维构造反演和速度层析成像

本文给出利用空间深地震测深反射资料重建地壳三维构造和速度分布的方法，适用于任意分层并且每层可由若干断块组成的地壳模型。通过采用适当的反演方法获得地壳界面的三维构造，并在此结果基础之上，进一步将剩余走时残差归因于地壳速度的不均匀性，采用东分块的反演方法，重建地壳三维速度分布，数值模拟的结果表明算法是有效的。     

基于交叉梯度耦合的大地电磁与地震走时资料三维联合反演

为了有效解决目前大地电磁和地震走时资料单方法反演结果一致性不好的问题,同时克服基于岩石不同物性参数间关系耦合约束联合反演的局限性,本文研究了基于交叉梯度耦合约束的大地电磁与地震走时资料的三维联合反演算法.以较为成熟的天然地震走时资料三维正反演和大地电磁三维正反演算法为基础,实现了具有共同的反演网格,以交叉梯度结构耦合约束,并能同时获得电阻率和速度模型的三维联合反演算法.分别利用单棱柱体模型和双棱柱体模型合成数据进行了联合反演试算.结果表明:无论是单棱柱体模型还是双棱柱体模型,联合反演结果比单独反演对异常体的空间形态都有更好的恢复,其中单棱柱体模型反演的异常体电阻率更接近于真实电阻率,双棱柱体模型的联合反演结果不仅消除了围岩的部分电阻率假异常,而且增强了对异常体深部速度结构特征的恢复程度.联合反演还能同时改善电阻率和速度反向变化异常体的单独反演结果,进一步证明交叉梯度耦合不依赖于岩石物性关系,而强调地下结构的相似性,具有更普遍的适用性.