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

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

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

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

简介：阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期（超过100秒）,逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。

简介：这篇文章重点研究改进的Gabor小波（improvedGaborwavelet，IGW）变）并讨论了它在地震信号处理和解释中的应用ThispaperintroducesanimprovedGaborwaveletanditscompletetransform，andmainlyanalysestheirpropertiesanddiscussesapplicationsofthesepropertiesinseismicsignalprocessandinterpretation。改进的Gabor小波变换具有以下特性：1）IGWT把时域信号映射到时间一频率域，而传统Gabor小波变换把时间信号映射到时间一尺度域；2）IGWT可用于信号分频，通过固定变换的主频参数dominantfrequency，并变换能提取相应的子带信号，且其主频部分的信息与原信号相应频率部分的信息一致，通过调节变换的分辨率因子，变换能有效控制子带信号的带宽；3）用IGWT和IGWIT构建的滤波器有良好的时一频局部性，在指定时一频范围内能实现针对性滤波。文章用仿真实验和实际用例验证IGWT的这些特性，并在提高地震信号分辨率、地震信号分频和识别小断层等地震信号处理和解释等方面的应用中取得良好效果。

简介：地震数据的精确解释和分析非常依赖于所使用算法的稳定性。我们着重于盐丘地震勘探的稳定检测。我们讨论一个基于地震成像中最佳结构特征属性排序分类的盐丘探测新模型。该算法克服了现有的基于结构属性技术的局限性，因为该技术非常依赖于盐丘固有的地质属性与盐丘检测所用的属性数量。该算法综合了灰度共生矩阵（GLCM）属性，Gabor滤波器，以及含有使用属性特征排序信息协方差矩阵的本征结构等属性。将排序前列的属性组合起来形成一组最优的特征集，以保证算法即使在沿盐丘边界没有强反射层的情况下也能有效。与现有的盐丘检测技术相比，本文的算法稳定和计算高效，并能处理小尺度特征集。我用荷兰F3地块评价该算法的性能。实验结果表明，本文提出的基于信息理论的工作流程用于检测盐丘，其精度优于现有盐丘检测技术。

简介：逆时偏移被认为是处理陡倾角断层、回转波和逆掩推覆构造最有效的成像方法。考虑到当今勘探对地下复杂构造体、薄互层和隐蔽储层的高分辨率成像需求，井间、井中地震资料成为地面地震资料的一个很重要的补充。但由于地面地震、井中地震和井间地震存在很大的频率差异，本文将相应的变网格算法应用到偏移中。另外，也将Lanczos滤波引入到波场外推计算中，从而可以很好地消除变网格引起的人为反射。试算表明，对比地面、井中和井间的偏移，可以看到微构造、陡反射层都得到了更好地刻画和描述。基于变网格的逆时偏移在储层预测和岩性解释方面有广阔的应用前景。

简介：随着并行计算技术的发展,非线性反演计算效率在不断提高,但对于基于单点搜索的非线性反演方法,其并行算法的实现则是一个难题。本文将群体搜索的思想引入到基于单点搜索的非线性反演方法,构建了并行算法,以量子蒙特卡罗方法为例进行了二维地震波速度反演及实际资料波阻抗反演,并测试了使用不同节点数进行计算的效率。计算结果表明：该并行算法在理论和实际资料反演中是可行的和有效的,具有很好的通用性;算法计算效率随着使用节点数的增加而提高,但算法计算效率的提高幅度随着使用节点数的增加逐渐减小。

简介：在井下微震监测中，各向异性模型用来获取事件位置和给出对介质实际描是十分有用的。在各向异性介质中，因为横波的分裂和奇异点的出现，使得波场的结构变得复杂。我们的研究展示了对强烈各向异性VTI介质利用速度模型校准和双偶极震源进行常规处理和运动学反演的结果。因小走时的拟合差，常用的质量评价标准，不是总正确的，人工加入各向同性介质层以降低拟合差，也可能产生非物理模型，负泊松比，事件位置的位移，

简介：本文描述了在全波实验中所见到的三种现象,即:1.纵横波初至的振幅(绝对值)具有相同的变化规律;2.纵横波初至的相位彼此总是反向的;3.纵横波初至的振幅变化周期为2π。经过充分的分析和解释指出这些现象应该出现在全波测井条件下,从而为利用相位差从全波中提取有用的信息奠定了基础。更多还原

简介：静冰压力是影响水工结构物安全运行的主要因素之一。现有静冰压力检测方法主要存在两个局限，一是实时性差，二是传统电阻应变式压力膜盒传感器因边壁效应影响测量结果的正确性。针对以上两个不足之处，作者基于反射式强度调制原理，研发了一种带杆伞状圆盘结构的压力膜盒静冰压力光纤传感器，实现了利用光纤传感技术检测静冰压力。我们用研制的新型光纤传感器在实验室完成了-30℃～5℃温度范围内冰生长过程中及5℃～-30℃温度范围内冰消融过程中静冰压力的测量试验。实验结果表明，与传统电阻应变式压力膜盒传感器相比，该传感器工作稳定、分辨率达0．02kPa、灵敏度达2．74×10-4／l（Pa，并可实现自动连续测量。

简介：静冰压力是影响水工结构物安全运行的主要因素之一。现有静冰压力检测方法主要存在两个局限,一是实时性差,二是传统电阻应变式压力膜盒传感器因边壁效应影响测量结果的正确性。针对以上两个不足之处,作者基于反射式强度调制原理,研发了一种带杆伞状圆盘结构的压力膜盒静冰压力光纤传感器,实现了利用光纤传感技术检测静冰压力。我们用研制的新型光纤传感器在实验室完成了-30℃~5℃温度范围内冰生长过程中及5℃~-30℃温度范围内冰消融过程中静冰压力的测量试验。实验结果表明,与传统电阻应变式压力膜盒传感器相比,该传感器工作稳定、分辨率达0.02kPa、灵敏度达2.74×10-4/kPa,并可实现自动连续测量。更多还原

简介：研究了基于Poynting矢量的角度域逆时偏移成像及成像幅度的校正方法。根据Poynting矢量进行波场角度分解，由此构建局部成像矩阵及局部照明矩阵。在局部成像矩阵中建立的角度域成像条件，有效地消除了低波数干扰，同时可在局部成像矩阵中进行角度域共成像点道集抽取、倾角估计等运算。利用局部照明矩阵进行了基于全波波动方程的时域照明分析，在局部照明矩阵中计算倾角域幅度校正因子。根据逆时偏移的像计算共倾角像，利用校正因子对各角度的像进行校正，进而实现对成像结果的校正，从而实现了一种高效的倾角域幅度校正方法。最后通过SEG／EAGE模型进行数值计算验证了文中所述的计算方法。

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

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

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

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

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

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