地震转换波对叠前同时反演的影响研究

叠前同时反演是油气探测的一种有效工具.其理论基础是平面P波Zoeppritz方程计算的反射系数的近似,是入射角的函数.叠前同时反演可以利用三项或两项Fatti方程进行反演分析.本文针对实际油田的测井数据,利用反射率法模拟了仅包含P波一次反射记录,包含P波一次反射和P波层间多次波记录以及全波场地震记录,再利用叠前同时反演对合成地震记录进行反演研究.研究结果表明,在大偏移距处P波主要反射受到其它模式波的污染,从而影响了叠前同时反演结果的精度.对于薄互层介质当转换波影响严重时,使用小角度数据的两项AVO反演比使用大角度数据的三项AVO反演更合理可靠.     

基于全区视电阻率的线源FCSEM二维正演模拟

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.     

复杂三维介质的大地电磁正演模拟

存在多个三维且不规则异常体的复杂模型,与实际地质模型更为接近,用有限单元法数值模拟计算无疑具有明显优势.但由于刚度系数矩阵为大型、不定且严重病态的,因此普遍存在求解效率低的问题.在引入不完全LU分解改善其条件数后,分析、比较了几种流行的Krylov子空间方法:CGS、BICG和BICGSTAB,得出BICGSTAB具有收敛平滑、求解速度快等优点,更适合于求解大地电磁三维有限单元法正演中的线性方程组问题.通过与一维模型解析计算结果的对比以及COMMEMI-3D模型的正演计算结果,验证了算法和编制程序的有效性,为后续构建大地电磁三维反演提供基础.     

地下隐蔽水库的高密度电法探测研究

地下空间开发过程中,要尽量避开地下已建成的建筑物和市政设施,然而由于各种原因造成地下已建成的建筑物、市政设施等具体位置不明确,形成为地下空间开发过程中的障碍物,对设计和施工造成不利影响.高密度电法具有成本低,效率高,信息丰富,解释方便等特点,威为解决这些地下障碍物问题的物探方法之一.为了探测在地铁选线过程中遇到的地下隐蔽水库,选用了高密度电法进行探测,并且依据调研到的资料建立了正演模型.模拟结果表明高密度电法可以圈定地下隐蔽水库的边界范围.最终的探测结果和正演结果基本一致.这说明在地下空间开发过程中采用高密度电法探测地下隐蔽水库具有可行性,探测的结果能为设计和施工提供参考依据.同时,在探测实施时根据已知条件进行正演模拟有助于解译高密度电法实测结果.     

可控源音频大地电磁测深法勘查深度研究

可控源音频大地电磁测深法的勘查深度问题在理论研究和实际工作中都得到了广泛的关注,本文分析了水平层状大地电磁正反演方法的基本原理,说明了可控源音频大地电磁测深法在深部勘查方面具有可靠性.通过建立两个三层地电模型,并对吉林某地区的CSAMT实测数据进行了带源反演成像,与已知钻井资料对比分析表明,可控源音频大地电磁测深法对地下1000~2000 m深度的反演较为准确,能够反映地下电性的垂向变化.     

下扬子黄桥地区构造特征及其与油气的关系

下扬子黄桥地区位于下扬子对冲构造带上,记录了多期构造变形,并且保存了较完整的中古生界海相地层,发育多 套烃源岩,具有良好的油气前景。文中采用三维地震工区构造精细解释及运动学正演模拟方法,分析了研究区的构造变形 特征及其变形过程,结果表明黄桥地区中上古生界和下中生界地层在自SE向NW的挤压应力作用下,以前展式发育了多个 断层相关褶皱构造。在构造前锋,断层相关褶皱反向发育,于工区中部偏北西侧形成明显的对冲构造,地层遭受强烈的褶 皱、隆升,构造总缩短率高于47.55%;之后,区域应力场转变为NW-SE向拉张,工区南东侧发育负反转断层,且自南西 向北东负反转程度逐渐变大,由上正下逆复合型转变为似正断层型;此次拉伸正断作用导致断层回滑,形成上白垩统填充 的半地堑构造,对冲构造也受到了一定的破坏。在此基础上,结合前人油气成藏的研究成果,探讨了黄桥构造变形与其油 气成藏的关系,并提出了对冲构造背斜、对冲构造下盘及南东侧的广陵构造背斜为黄桥地区的有利油气区块。研究结果对 进一步认识黄桥及下扬子区的构造变形特征和演化、以及该区的油气勘探提供一定的参考。     

基于谱元法的各向异性黏弹性介质地震波场模拟和分析

基于谱元法计算理论,通过坐标旋转算法实现任意各向异性黏弹性介质中的地震波场计算,研究复合介质中地震波场特征。根据计算结果,qSV波表现为剪切滞弹性形变,qP波表现为膨胀滞弹性形变;品质因子不仅影响地震波振幅,还影响其传播频率;黏弹性介质吸收高频快,在一定频率范围内,黏弹性介质没有频散现象。     

基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案

作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.     

基于COMSOL MULTIPHYSICS的直流电法正演模拟

以点电源和线电源为例,叙述了基于COMSOL MULTIPHYSICS的直流电法正演模拟过程。通过对COMSOL MULTIPHYSICS正演模拟解与理论值的对比分析,验证了应用COMSOL MULTIPHYSICS进行直流电法正演模拟的正确性与可行性。鉴于COMSOL MULTIPHYSICS强大的剖分求解功能及丰富的后处理操作,可以快速、准确实现正演模拟,其在地球物理正演模拟中的应用的研究意义重大。      

张性盆地裂后异常沉降的正反演数值模拟方法

世界上许多张性沉积盆地存在远大于McKenzie模型理论值的裂后异常沉降, 南海北部陆缘的沉积盆地也是如此; 确定裂后异常沉降的特征和分布是认识其成因机制和对油气成藏影响的前提.介绍了估算裂后异常沉降的3种方法: 古水深比较法、应变速率反演法和沉降过程二维正反演法, 并指出了各方法的应用前提和优缺点.对于由作者提出的后一种方法还结合在珠江口盆地的应用实例进行了较详细的讨论, 表明这种方法能在考虑岩石圈挠曲强度的基础上正演模拟出单幕或多幕盆地沉降及相应的岩石圈伸展系数, 从而计算出盆地理论热沉降, 与通过回剥反演得到的实测构造沉降进行对比; 还指出了该方法存在的问题和需进一步研究之处.