简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。
简介:阵列声波信号是典型的非线性、非平稳信号,Hilbert~Huang变换(HHT)是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解(EMD)和对瞬时频率的求解,可以获得声波信号的时一频谱。其关键技术就是进行经验模态分解,任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征,在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程,并列举了一些声波测井波列实例,说明了该方法的有效性。
简介:阵列声波信号是典型的非线性、非平稳信号,Hilbert-Huang变换(HHT)是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解(EMD)和对瞬时频率的求解,可以获得声波信号的时-频谱。其关键技术就是进行经验模态分解,任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征,在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程,并列举了一些声波测井波列实例,说明了该方法的有效性。
简介:阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以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结构区域内的电场参数估计发挥的重要作用。
简介:随着全张量重力梯度(FTG)测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU(GraphicsProcessingUnit)和预处理分解技术(Preconditioningmethods)在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛(SSOR)技术与不完全乔列斯基分解共轭梯度算法(ICCG)提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATeslaC2050GPU的并行SSOR-ICCG算法和在2.0GHzCPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。
简介:山区重力勘探中,重力外部校正是重力勘探料处理和解释的前提,其内容包括地形校正、中间层校正等.本文在分析常规外部校正方法的基础上,提出了一套适合于山区重力勘探的球面地形校正及有限球壳中间层校正方法,并严格推出了球面地形校正公式.对中国南方ZJJ地区重力资料处理结果表明,该方法提高了地形校正的精度,使山区重力资料品质达到或接近平原区的水平.
简介:从大量的地震属性中提取最能反映地质特征的综合属性是储层预测技术的关键,通常选用降维方法来优选属性。目前应用最为广泛的线性降维方法。但是,由于地震属性与地质特征的关系通常是非线性的,基于线性变换的地震属性降维优化方法不能充分地反映这种非线性关系,降低了储层预测的精度。流形学习是一种新的非线性学习方法,它是通过保持数据局部结构的方式将高维数据投影到低维空间,挖掘和发现隐藏在数据中的内在特征与规律性,开拓了地震属性降维优化研究的新领域。本文首次实现了3D地震数据的层问属性特征提取,讨论了LLE方法及其关键技术,并以奥陶系礁滩相储层实例说明LLE和PCA两种方法降维及聚类的不同效果。理论模型分析和实例应用表明:LLE较好地保持了数据本身的原始结构;提取的综合属性和聚类相图较好地刻画了沉积相带、储层和流体的特征。这说明流形学习具有更好的特征提取性能。
简介:井地电位成像是通过套管向井中供电或将电源放在井中,在地表观测电位异常的一项技术,其供电源有线源和点源两种类型。为了研究这两种电源对地下异常体产生的电位异常特征,本文针对不同激励源,采用有限差分方法进行数值模拟研究,在线性方程组求解电位时引入不完全Cholesky共轭梯度(ICCG)迭代方法,分别实现了点源和线源井地电位成像技术的三维正演。最后,基于阻尼最小二乘法实现了井地电位成像技术的电阻率三维反演。设计不同地电模型分别进行正演和反演试算,正演结果表明,供电电源的类型不同,异常体在地表的电位异常特征也不同;反演结果表明,低阻体的反演结果要好于高阻体,点源置于异常体下方时反演的电阻率对异常体边界的识别比线源更加准确。
简介:类中梯装置三维激电成像技术兼顾了激电测深与激电剖面二者的特点,能快捷高效的获得测区三维电阻率及极化率信息,空间分辨率高,探测深度大。本文从类中梯装置的数据采集方式入手,阐述了采用类中梯装置进行三维激电成像的工作方式.建立了一个三维地电模型,采用类中梯装置进行三维正演模拟,并进行反演及成像。数值模拟结果表明,采用类中梯装置进行三维激电成像能较好地刻画实际地电模型的特点。以甘肃省某多金属矿为例,利用阵列式电磁法综合测量系统,进行了三维大功率激电成像技术的应用研究。实际应用结果表明,采用三维激电成像技术,可以多角度、多细节显示测区地下介质电阻率、极化率的分布情况,明确电性突变界面的延展状态,快速圈定成矿有利区。该研究对多金属矿产勘探等领域具有一定的指导性意义。
简介:通过地震数据获取裂缝储藏中流体的性质并对流体类型进行识别,是地震勘探岩性反演的重要问题之一。由于地震波的速度、储层的密度等弹性参数对某些流体不具有很强的敏感性,使只依赖振幅信息进行流体识别的传统AVO方法面临困境。作为传统叠前振幅反演的一个拓展,频变AVO(FDAVO)技术进一步考虑了振幅对频率的依赖关系,将这种依赖关系与地下裂缝结构、流体填充对应起来,能带来更丰富的流体信息。利用该技术,本文提出了一种基于地震数据参数化Chapman模型的贝叶斯反演新方法(BIDCMP),它包含两步算法,即,FDAVO反演储层的非弹性属性和贝叶斯框架下的流体识别。首先,通过匹配观测数据和模型数据,构造差函数反演裂缝储层非弹性参数。随后,在贝叶斯框架下,使用马尔科夫随机场(MRF)作为先验模型,联合多参数场识别流体。本方法在计算过程中,除综合考虑了弹性参数场、测井资料等常规信息外,还特别地加人了第一步中反演得的非弹性参数的约束,从而充分利用了流体粘性差异,最后在最大后验概率(MAP)准则下输出最佳岩性一流体识别结果。分别对合成地震记录和模拟岩性—流体剖面验证本文方法的有效性,结果证明本文方法获得的流体识别结果准确可信。