简介:通过地震数据获取裂缝储藏中流体的性质并对流体类型进行识别,是地震勘探岩性反演的重要问题之一。由于地震波的速度、储层的密度等弹性参数对某些流体不具有很强的敏感性,使只依赖振幅信息进行流体识别的传统AVO方法面临困境。作为传统叠前振幅反演的一个拓展,频变AVO(FDAVO)技术进一步考虑了振幅对频率的依赖关系,将这种依赖关系与地下裂缝结构、流体填充对应起来,能带来更丰富的流体信息。利用该技术,本文提出了一种基于地震数据参数化Chapman模型的贝叶斯反演新方法(BIDCMP),它包含两步算法,即,FDAVO反演储层的非弹性属性和贝叶斯框架下的流体识别。首先,通过匹配观测数据和模型数据,构造差函数反演裂缝储层非弹性参数。随后,在贝叶斯框架下,使用马尔科夫随机场(MRF)作为先验模型,联合多参数场识别流体。本方法在计算过程中,除综合考虑了弹性参数场、测井资料等常规信息外,还特别地加人了第一步中反演得的非弹性参数的约束,从而充分利用了流体粘性差异,最后在最大后验概率(MAP)准则下输出最佳岩性一流体识别结果。分别对合成地震记录和模拟岩性—流体剖面验证本文方法的有效性,结果证明本文方法获得的流体识别结果准确可信。
简介:地震数据的精确解释和分析非常依赖于所使用算法的稳定性。我们着重于盐丘地震勘探的稳定检测。我们讨论一个基于地震成像中最佳结构特征属性排序分类的盐丘探测新模型。该算法克服了现有的基于结构属性技术的局限性,因为该技术非常依赖于盐丘固有的地质属性与盐丘检测所用的属性数量。该算法综合了灰度共生矩阵(GLCM)属性,Gabor滤波器,以及含有使用属性特征排序信息协方差矩阵的本征结构等属性。将排序前列的属性组合起来形成一组最优的特征集,以保证算法即使在沿盐丘边界没有强反射层的情况下也能有效。与现有的盐丘检测技术相比,本文的算法稳定和计算高效,并能处理小尺度特征集。我用荷兰F3地块评价该算法的性能。实验结果表明,本文提出的基于信息理论的工作流程用于检测盐丘,其精度优于现有盐丘检测技术。
简介:静冰压力是影响水工结构物安全运行的主要因素之一。现有静冰压力检测方法主要存在两个局限,一是实时性差,二是传统电阻应变式压力膜盒传感器因边壁效应影响测量结果的正确性。针对以上两个不足之处,作者基于反射式强度调制原理,研发了一种带杆伞状圆盘结构的压力膜盒静冰压力光纤传感器,实现了利用光纤传感技术检测静冰压力。我们用研制的新型光纤传感器在实验室完成了-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模型进行数值计算验证了文中所述的计算方法。
简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。
简介:提高地震资料分辨率的一个有效途径就是衰减补偿,通过对地震波的衰减和频散效应进行校正,提高地震资料的分辨率。常规衰减补偿方法都是基于波场延拓的反Q滤波方法。本文利用Futterman衰减模型,导出了一种衰减介质中合成地震记录的计算方法,在此基础上将衰减补偿问题归结为一个Fredholm积分方程反问题,利用反演方法来实现衰减补偿。针对衰减补偿问题的不稳定性,利用Tikhonov正则化方法提高反演过程的稳定性,数值模拟资料和实际资料处理结果验证了方法的有效性。
简介:本文将相干成像用于WVSP数据的一次反射波和下行地表多次波的成像中,通过对叠加次数分布的分析,发现一次反射波相干成像的成像范围小但在井位置处的叠加次数高;下行地表多次波相干成像的成像范围大而在井位置处的叠加次数少。基于此,本文尝试利用一种WVSP的联合相干成像方法来解决这一问题。该方法利用WVSP的直达波与下行地表多次波和直达波与一次反射波来进行联合相干成像。这种联合相干成像方法在有效拓宽成像面积的同时增加了覆盖次数(特别是在井旁),使干涉记录中的虚假同相轴得到了有效地压制,提高了成像剖面的信噪比,使成像结果更加可靠。正演模拟得到的WVSP合成记录和野外实际采集到的WVSP数据处理结果均表明这种方法的有效性。更多还原
简介:Tikhonov正则化(TR)方法在重磁数据处理中发挥了重要的作用,本文在研究如何利用Tikhonov正则化方法方法解决重力数据3D反演的同时,深入讨论了可进一步提高拟合误差的ExtrapolationTikhonov正则化方法(EXTR)的原理,并就其参数选择方法及各参数对拟合误差、迭代次数及反演结果的影响进行研究。常密度及变密度组合模型试算结果表明,与TR方法相比,EXTR方法不仅可以达到解释人员设定的先验拟合误差水平,在计算时间及迭代次数相应增加的前提下有更高的拟合精度;同时其反演结果也更加紧致,进一步改善了TR反演结果的发散性;并且其反演数据范围更贴近预设模型参数范围,模型特征与预设模型密度分布吻合较好。更多还原
简介:有限频走时层析成像是近年发展起来的一种新方法,这种新方法的一个主要过程是走时灵敏度算核的计算。求解灵敏度算核要多次用到同一散射点的走时,多次地求解同一走时是相当耗时的任务,如果介质为均匀或速度线性变化等简单模型,散射点的走时可以用解析公式快速地求出,从而灵敏度算核的计算耗时相对较少。然而各种地球模型中,介质速度大多为分层模型,从解析公式中得到走时信息就比较繁锁。为了提高计算效率,本文采用查表算法研究地球分层速度模型中的P震相有限频走时灵敏度算核的计算,选用的速度模型是地球AKl35模型,用查表算法求解走时,节省了约50%的计算时间。在相同的速度模型下,与已有结果的对比,本文所用的查表算法,能在兼顾精度的前提下,以较小的存储要求换取较高的计算效率,这对提高有限频走时层析成像算法的速度具有一定的参考价值。
简介:随着全张量重力梯度(FTG)测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU(GraphicsProcessingUnit)和预处理分解技术(Preconditioningmethods)在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛(SSOR)技术与不完全乔列斯基分解共轭梯度算法(ICCG)提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATeslaC2050GPU的并行SSOR-ICCG算法和在2.0GHzCPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。