简介:通过引入岩石骨架参数表达式到Gassmann近似流体方程中,实现了基于该方程的地震孔隙度反演。但是岩石骨架参数表达式有很多种,使得地震孔隙度反演公式也各不相同,实际应用不方便,反演结果也难以评价。作者在研究现有常用岩石骨架参数表达式如Esheby-Walsh,Pride,Geertsma,Nur,Keys-Xu以及Krief等的基础上,提出了含有调节参数的岩石骨架模型统一表示式,发展了基于Gassmann方程的地震孔隙度反演方法。该方法应用范围宽,参数调节方便、灵活。为验证所提公式和方法的实用性,作者们结合ZJ地区钻井岩心的岩石物理样品测试数据和测井数据,进行了地震孔隙度反演。由于该地区油气储层的存在与中等程度孔隙度的分布有关,因此地震孔隙度反演对预测和识别油气层、干层和含水层有重要作用。反演结果与工区内孔隙度测井瞌线吻合度高,说明本文所提公式与方法是有效的,可靠的。
简介:野外空间采样密度的提高将增加室内数据分析的工作量,常规的基于地震数据的点与线分析方法有一定的局限性。本文简要说明空间子集的抽取方法,列举了正交子集和斜交子集的特点,并通过三维可视化展示了子集数据的空间特性。提出在数据处理中利用子集的时间切片:(1)分析地震道空间分布的均匀性和规则性;(2)研究面波及规则干扰的空间分布特点;(3)检测叠前数据中的异常信息;(4)监控叠前去噪的效果。实际数据的应用结果表明基于空间子集的分析方法是一种独特有效的地震数据分析方法,从另一种视角观察地震数据,可以发现其中某些新的特征,以提高处理人员对数据的洞察能力。
简介:常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。
简介:边界识别是重磁数据解释中的常用方法之一,依据其结果可划分出地质体的水平范围。边界识别结果受地质体埋深及导数计算误差的影响所识别边界与真实边界之间存在一定的差距,且边界识别法无法直观地给出地质体的深度信息。为了获得异常体的水平位置和深度信息,本文提出空间归一化边界识别方法,其对不同深度的边界识别函数进行归一化计算,空间归一化边界识别法的最大值对应于异常体的水平位置和深度。常规边界识别结果的误差随理深的减小而减小,而空间归一化边界识别法是通过最大值来判断地质体的位置,最大值是在地质体处获得,因此归一化边界识别方法所获得的结果是准确的。通过理论模型试验证明归一化边界识别方法能有效地完成异常体的水平位置和深度的计算,所获得的水平位置和深度信息与理论值相一致,为下一步的勘探计划提供了更加可靠的依据。将其应用于实际航磁数据的解释,获得了断裂的具体分布形式。更多还原
简介:频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。
简介:由于孔隙型碳酸盐岩储层的复杂孔隙结构和强非均质性,岩石导电效率与含水孔隙度之间理想的线性关系相对复杂或并不存在,故前人基于岩石导电效率理论建立含水饱和度的计算公式并不完全适用于碳酸盐岩储层。基于岩石导电效率理论,推导了岩石导电效率的计算公式,阐明了岩石导电效率与含水孔隙度之间线性关系相对复杂或并不存在的根本原因,发现了岩石导电效率与电阻率之间的幂函数关系,分析了所建立含水饱和度计算公式误差的主控因素,得出了岩石导电效率的计算精度是该方法是否能推广应用的关键。与Archie公式相比,基于岩石导电效率理论建立的含水饱和度计算公式能更准确地计算研究靶区孔隙型碳酸盐岩储层的含水饱和度。在伊拉克某油田和印尼某气田3口井碳酸盐岩储层中的应用,表明当计算岩石导电效率的相对误差不大于0.1时,所计算储层含水饱和度的绝对误差不大于0.1,基本满足孔隙型碳酸盐岩储层精细评价的需求。
简介:有限频走时层析成像是近年发展起来的一种新方法,这种新方法的一个主要过程是走时灵敏度算核的计算。求解灵敏度算核要多次用到同一散射点的走时,多次地求解同一走时是相当耗时的任务,如果介质为均匀或速度线性变化等简单模型,散射点的走时可以用解析公式快速地求出,从而灵敏度算核的计算耗时相对较少。然而各种地球模型中,介质速度大多为分层模型,从解析公式中得到走时信息就比较繁锁。为了提高计算效率,本文采用查表算法研究地球分层速度模型中的P震相有限频走时灵敏度算核的计算,选用的速度模型是地球AKl35模型,用查表算法求解走时,节省了约50%的计算时间。在相同的速度模型下,与已有结果的对比,本文所用的查表算法,能在兼顾精度的前提下,以较小的存储要求换取较高的计算效率,这对提高有限频走时层析成像算法的速度具有一定的参考价值。
简介:传统的f-x域经验模态分解法(Empiricalmodedecomposition,EMD)能够有效地对主要由水平同相轴构成的地震记录进行随机噪声衰减。然而,当同相轴倾斜时,f-x域经验模态分解法在衰减随机噪声的同时去除大部分有效信号。本文提出了一种基于f-x域经验模态分解法的改进算法。我们通过局部相似度对所去除的噪声信号中的有效信号进行提取。局部相似度可以用来检测噪声信号中的有效信号点并用来构造一权重算子进行信号提取。新方法与f-x域经验模态分解法、f-x域预测滤波法以及f-x域经验模态分解预测滤波法相比能够在衰减随机噪声的同时保留更多的有用信号。数值模拟实验以及实际地震资料处理结果均表明该方法能更为有效地去噪。
简介:通过最新的全球地磁模型——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)处衰减较快,在地表处则趋于稳定。