简介：现阶段的叠前地震反演技术中用于描述反射系数与纵、横波速度和密度之间的关系几乎完全是Zoeppritz方程的近似式，由于这些近似公式在大角度和弹性参数变化剧烈时误差较大，这不仅降低了反演解的精度，而且增加了叠前反演的多解性。本文探索了直接利用Zoeppritz方程求解精确反射系数的理论方法，并基于广义线性反演理论详细推导了基于叠前大角度地震资料的纵、横波速度和密度三参数同步反演算法，同时在反演过程引入正则化约束阻尼因子和共轭梯度算法，有效降低了反演的不适定性和提高了反演收敛性。理论模型试算和实际工区应用表明，本文提出的反演方法能够有效利用大角度（一般入射角〉30°）的叠前地震数据，获得更精确的地震弹性参数反演结果，并且反演结果忠实于地震资料，与井吻合较好。

简介：本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。

简介：由于CRS叠加考虑了反射层的局部特征和第一菲涅耳带内的全部反射，从而更充分地利用了多次覆盖反射数据的信息。就目前的地震资料处理技术而言，它是最佳的零偏移距成像方式。本论文利用改进型的参数优化技术，得到高质量的CRS运动学参数剖面，并利用参数剖面计算出叠加孔径，实现了基于最优孔径的CRS叠加，使CRS参数的用途得到了充分利用。模型数据和实际资料的试算表明，基于最优孔径的CRS叠加的成像剖面与传统CRS叠加剖面相比，有着较高的信噪比和同相轴的连续性。

简介：阵列声波信号是典型的非线性、非平稳信号，Hilbert～Huang变换（HHT）是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解（EMD）和对瞬时频率的求解，可以获得声波信号的时一频谱。其关键技术就是进行经验模态分解，任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征，在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程，并列举了一些声波测井波列实例，说明了该方法的有效性。

简介：提高地震资料分辨率的一个有效途径就是衰减补偿,通过对地震波的衰减和频散效应进行校正,提高地震资料的分辨率。常规衰减补偿方法都是基于波场延拓的反Q滤波方法。本文利用Futterman衰减模型,导出了一种衰减介质中合成地震记录的计算方法,在此基础上将衰减补偿问题归结为一个Fredholm积分方程反问题,利用反演方法来实现衰减补偿。针对衰减补偿问题的不稳定性,利用Tikhonov正则化方法提高反演过程的稳定性,数值模拟资料和实际资料处理结果验证了方法的有效性。

简介：推导了二维TTI介质的相速度表达式,并且依据推导出来的相速度表达式,模拟并分析了二维TTI介质相速度的传播快照以及TI介质相速度的传播快照;对比并分析了TTI介质和TI介质模型的相速度理论计算值的X分量特征的差异。TTI介质的相速度研究具有较高的理论研究价值和实际应用价值．

简介：信噪比和分辨率是评价地震数据质量的两个重要的，但相互矛盾的指标。为了在提高分辨率的同时，保持地震数据的信噪比，引入信号纯度谱的概念并使用它来确定反褶积之后希望输出的振幅谱。由于实际反射系数序列具有蓝谱特征，并非白谱，在针对反射系数白噪假设对子波的影响进行实验分析的基础上，对谱白化之后的地震记录进行有色补偿处理。实际地震数据处理表明，基于信号纯度谱的谱白化与有色补偿的串联应用有效地改善了地震数据反映薄层细节的能力。

简介：阵列声波信号是典型的非线性、非平稳信号,Hilbert-Huang变换(HHT)是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解(EMD)和对瞬时频率的求解,可以获得声波信号的时-频谱。其关键技术就是进行经验模态分解,任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征,在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程,并列举了一些声波测井波列实例,说明了该方法的有效性。

简介：本文将相干成像用于WVSP数据的一次反射波和下行地表多次波的成像中,通过对叠加次数分布的分析,发现一次反射波相干成像的成像范围小但在井位置处的叠加次数高;下行地表多次波相干成像的成像范围大而在井位置处的叠加次数少。基于此,本文尝试利用一种WVSP的联合相干成像方法来解决这一问题。该方法利用WVSP的直达波与下行地表多次波和直达波与一次反射波来进行联合相干成像。这种联合相干成像方法在有效拓宽成像面积的同时增加了覆盖次数(特别是在井旁),使干涉记录中的虚假同相轴得到了有效地压制,提高了成像剖面的信噪比,使成像结果更加可靠。正演模拟得到的WVSP合成记录和野外实际采集到的WVSP数据处理结果均表明这种方法的有效性。更多还原

简介：地震波场和电磁场耦合产生的震电效应与储层物性参数有关，含流体孔隙介质中震电效应的研究有助于更好的描述储层特性。本文我们对饱和砂岩样品中的震电效应进行了实验研究，构建了一套震电测量装置。利用此测量装置记录了在声波激励下砂岩样品中产生的震电转换信号，观测得到了砂岩界面产生的震电信号的衰减特性，在此基础上进一步研究了震电信号与砂岩物理参数之间的关系。在震电效应的实验中发现尽量减小参考电位与震电信号扰动区的基准电位之间的电位差是保证震电实验测量精确性的关键点，能够显著提高震电信号的可探测性。震电测量结果证实了地震波在含流体孔隙介质中传播时，能诱导震电耦合，同时给出了震电界面响应的特点。震电信号振幅随着波源与岩样距离的增加呈线性衰减，随着接收电极与岩样界面距离的增加呈指数衰减。不同渗透率砂岩样品中产生的震电响应结果表明震电响应的强弱与样品渗透率成正相关，震电效应可以作为研究储层渗透率的一种新方法。

简介：有限频走时层析成像是近年发展起来的一种新方法,这种新方法的一个主要过程是走时灵敏度算核的计算。求解灵敏度算核要多次用到同一散射点的走时,多次地求解同一走时是相当耗时的任务,如果介质为均匀或速度线性变化等简单模型,散射点的走时可以用解析公式快速地求出,从而灵敏度算核的计算耗时相对较少。然而各种地球模型中,介质速度大多为分层模型,从解析公式中得到走时信息就比较繁锁。为了提高计算效率,本文采用查表算法研究地球分层速度模型中的P震相有限频走时灵敏度算核的计算,选用的速度模型是地球AKl35模型,用查表算法求解走时,节省了约50%的计算时间。在相同的速度模型下,与已有结果的对比,本文所用的查表算法,能在兼顾精度的前提下,以较小的存储要求换取较高的计算效率,这对提高有限频走时层析成像算法的速度具有一定的参考价值。

简介：在储层发育带，通常地震波波形会变复杂，所以关联维会变小。但由于地震波是有一定带宽的信号，波形受这个带宽内所有频率成分的影响，如果直接进行关联维的计算，其结果往往不能很好地反映储层。鉴于这种情况，提出了经验模态分解（EMD）与关联维相结合的方法，利用EMD对地震波进行分解，对分解后的每一个平稳的固有模函数（IntrinsicModeFunction，IMF）分量进行关联维计算，将计算结果与有效的IMF分量进行储层预测。利用该方法对xx区三叠系中油组进行了试算，取得了较好的效果。

简介：野外空间采样密度的提高将增加室内数据分析的工作量，常规的基于地震数据的点与线分析方法有一定的局限性。本文简要说明空间子集的抽取方法，列举了正交子集和斜交子集的特点，并通过三维可视化展示了子集数据的空间特性。提出在数据处理中利用子集的时间切片：（1）分析地震道空间分布的均匀性和规则性；（2）研究面波及规则干扰的空间分布特点；（3）检测叠前数据中的异常信息；（4）监控叠前去噪的效果。实际数据的应用结果表明基于空间子集的分析方法是一种独特有效的地震数据分析方法，从另一种视角观察地震数据，可以发现其中某些新的特征，以提高处理人员对数据的洞察能力。

简介：针对波前构建法中非规则的波前四边形网格与规则的矩形网格节点之间的相对定位问题,对比分析了在计算机图形学中常用的四种网格定位方法(叉乘判断法、夹角之和检验法、射线法(交点计数检验法)以及符号比较法)在波前构建法中的具体实现问题,通过对定位方法稳健性分析及计算实例表明,矢量叉乘判断法在实现规则的矩形网格节点与非规则的波前四边形网格相对定位研究中,其稳健性和计算效率是最高的,从而为下一步做偏移或者反偏移提供了精确的网格点属性值。更多还原

简介：频率域全波形反演虽然克服了时间方向上的局部极小值问题,但是地下介质的复杂性使其在空间域仍然存在局部极小值缺陷。在优化梯度法基础上,本文采用预条件双共轭梯度稳定算法和多重网格方法计算反演中的波场传播和目标函数的梯度,在保证计算速度的同时,减小计算机内存的消耗。频率域波形反演和多重网格的多尺度性质有�

简介：通过最新的全球地磁模型——CALS10K.1b,结合CALS3K.4与IGRF11模型,计算并分析了10000BC~1990AD期间中国大陆及邻近地区非偶极子（ND）磁场Z分量的时空变化。为了深入了解ND场的变化,从场源的角度,对2n（n=2-10）极子ND场及其对应的能量进行了分析。结果显示在研究期间ND场的变化可分为3个阶段。在10000BC-2500BC期间,ND场以正值为主并持续了近7500a,在2500BC-1500AD期间强度转弱为以负值为主并持续至1500AD,自此快速增强为以正值为主。东亚地区ND场异常基本在截断阶数（n）为3时即形成,且该异常区已在大陆地区内形成了封闭的圆形区域,这意味着前3阶的ND场占据了总ND场强度的大部分。ND场在核幔边界（CMB）处衰减较快,在地表处则趋于稳定。

简介：随着全张量重力梯度（FTG）测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU（GraphicsProcessingUnit）和预处理分解技术（Preconditioningmethods）在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛（SSOR）技术与不完全乔列斯基分解共轭梯度算法（ICCG）提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATeslaC2050GPU的并行SSOR-ICCG算法和在2.0GHzCPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。

简介：在地震勘探数据采集中,随机噪声严重影响地震资料质量,给后期解释工作带来很大困难。如何在不损失剖面有效信息的前提下压制随机噪声,有效地提高地震资料的信噪比和保真度,是本文的研究目标。构造导向滤波技术的核心是构造方向表征的求取以及如何实现非平稳滤波,来达到提高地震数据信噪比和保真度的目的。本文首先通过分析函数二维导数与希尔伯特变换的频率响应关系,推导出了基于二维希尔伯特变换的非迭代地震同相轴倾角求取算子,进而达到了构造方向表征的求取;其次选取多项式拟合作为构造导向滤波中的非平稳滤波方法,扩展了非平稳多项式拟合的应用范围;最后沿构造倾角方向进行变振幅同相轴的非平稳多项式拟合,实现和构建了新的自适应构造导向滤波方法。理论模型和实际地震资料处理的结果表明,所提出的方法实现了既保护构造信息又有效地压制了随机噪声的目的。

简介：Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。

简介：相干体是一种能够有效反映断层和裂缝等地质上非连续性的地震资料解释工具。然而由于受到附近地层,尤其是具有较强反射系数的地层的影响,微小断层和裂缝在相干体的时间切片和纵剖面上仍然难以辨别。本文提出了一种基于局部直方图规定化的相干体增强方法,实现相干体中的微小断层和裂缝的增强。与传统直方图规定化方法不同,本文方法处理三维相干体时不需要对数据进行离散化以统计其分布情况。在将相干体划分为若干子块后,以相干体整体的分布函数作为每个子块内的目标分布函数进行直方图规定化。另外,相邻子块的部分区域需要重叠,以克服分块产生的边界效应。对实际数据的处理结果表明,该方法可以提供更多的微小断层和裂缝的细节信息。