简介:针对不同矿物组分含量对盐膏岩速度影响大,盐膏岩速度复杂难题。本文把岩石物理模版的研究方法引入到盐膏岩的速度研究中,首次建立了盐膏岩四元组分(膏、盐、泥、孔隙地层水)含量与岩石速度关系的模版,在分析盐膏岩组分矿物骨架速度的基础上。基于多矿物测井解释模型,并以纵波模量比与孔隙度趋势作为约束条件,以组分含量为横坐标,岩石速度为纵坐标。在把不同尺度孔隙度值时,泥质、膏、盐矿物含量发生变化时,获得相应盐膏岩速度的变化值放在同一坐标系中,构成了盐膏岩四元矿物组分与岩石速度关系模版。并应用实测岩芯数据和测井数据,验证了该模版的可靠性和适应性,利用该模板能识别盐膏岩中关键结构和重要矿物的变化趋势,能直观的评估盐膏岩矿物组分蛮化对岩石涑席的影响.
简介:以阿姆河盆地右岸地区碳酸盐岩台内滩气田为例,运用基于图论的多分辨率聚类分析法开展了以常规测井资料为基础的测井相聚类及岩相定量预测研究。该方法不需要分析数据体的结构及聚类数等先验知识为指导,能够自动优选最佳聚类个数,并允许按照实际需求控制聚类级别,进行聚类合并。依据岩芯资料岩相识别及测井相标定结果,本文最终建立了包含5个聚类的测井相划分模型及岩相定量解释图版,其中,聚类测井相1~5分别对应于泻湖泥、石膏坪、滩间、低能滩及高能滩,取芯段符合率达85%以上,能够较好的运用于非取芯段岩相预测研究。据此,我们进行了连续的聚类测井相划分及岩相预测,并对层序地层格架内岩相分布及物性特征进行了分析。
简介:有限频走时层析成像是近年发展起来的一种新方法,这种新方法的一个主要过程是走时灵敏度算核的计算。求解灵敏度算核要多次用到同一散射点的走时,多次地求解同一走时是相当耗时的任务,如果介质为均匀或速度线性变化等简单模型,散射点的走时可以用解析公式快速地求出,从而灵敏度算核的计算耗时相对较少。然而各种地球模型中,介质速度大多为分层模型,从解析公式中得到走时信息就比较繁锁。为了提高计算效率,本文采用查表算法研究地球分层速度模型中的P震相有限频走时灵敏度算核的计算,选用的速度模型是地球AKl35模型,用查表算法求解走时,节省了约50%的计算时间。在相同的速度模型下,与已有结果的对比,本文所用的查表算法,能在兼顾精度的前提下,以较小的存储要求换取较高的计算效率,这对提高有限频走时层析成像算法的速度具有一定的参考价值。
简介:常规土地垂直地震介绍(VSP)探索通常为接收装置使用P波浪来源和三部件的地震检波器,强调P--并且变换S波浪。以前的研究证明从在表面的可控制的地震来源的两炸毁地上凿穿射击和垂直颤动将生产相对强壮的纯P波浪和更弱的纯S波浪。大泊松的比率差别的接口在强壮的播送变换S波浪的形成上有积极影响。由从来源的纯S波浪的比较分析并且在去的S波浪下面变换了,当时,我们相信纯S波浪的主要频率通常比纯P波浪低主要频率在变变换下面,S波浪接近P波浪的。我们学习了零偏移量并且从陆地P波浪的偏移量VSP数据采购原料。结果证明纯S波浪通常在波浪紧张与差别在这些数据存在。S波浪速度能从P波浪来源被获得零偏移量的VSP数据。最后,我们讨论VSPP的联合申请的明亮的未来--并且S波浪和在P波浪的S波浪的完整的使用采购VSP数据。
简介:噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。
简介:噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。更多还原
简介:储层物性参数的改变伴随着弹性参数某种程度的改变,而这种改变并非简单的线性关系。加上观测信息的缺失、重叠、噪声破坏及模型理想化等原因,致使这种改变关系的获取具有较大的不确定性。针对传统统计岩石物理物性参数反演方法的不足,以利用弹性参数反演储层物性参数为目的,依据贝叶斯反演框架,我们建立了新的储层物性参数目标反演函数。首先,采用兼具确定性与随机性特点的统计岩石物理模型,考虑到不同弹性参数间的精度存在差异,引入权重系数,建立起储层物性参数与弹性参数间的加权统计关系。其次,基于这种加权统计关系,结合马尔科夫链蒙特卡洛随机模拟技术产生储层物性参数、弹性参数随机联合样本空间作为目标函数求解样本空间。最后,建立解的快速求解准则,求取最大后验概率密度对应的储层物性参数取值作为最终解。实际应用表明,该方法具有较高的反演效率,应用前景较好。
简介:本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法(FastFourierTransformMovingAverage,FFT-MA)和逐渐变形算法(GradualDeformationMethod,GDM)得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。
简介:弹性参数在甜点区预测和页岩气的开发过程中扮演着重要的角色,因此研究等效弹性参数随页岩气储层属性的变化是一项很有意义的工作。研究中我们用x射线CT扫描技术获得了较为精确的页岩样品微观结构图像。从这些图像中,我们可以获得孔隙度和矿物的详细情况,据此,我们构建了三维数字岩心,并应用有限元法对弹性参数进行了数值模拟,其间深入考察了子样选取、网格划分、求解器类型以及边界条件等,该方法易于区别不同的矿物及其百分含量。本文重点研究孔隙度和干酪根含量对弹性参数的影响,计算结果表明,孔隙度和干酪根含量对弹性性质有较大的影响,当孔隙度和干酪根含量增加时,弹性模量降低,且当孔隙度小于0.75%左右、干酪根含量大于3%左右时弹性参数减小速率较缓。因为孔隙度仅仅为4.5%,孔隙中填充油或气对弹性参数的影响甚微。不同岩心样本具有不同的孔隙度和干酪根含量,传统岩石物理实验不仅昂贵而且费时,而数值模拟是基于数字岩心来计算弹性参数,更加经济、方便。本研究证实了将页岩样品的微观结构图像与弹性模量的计算相结合来预测页岩弹性参数的可行性。
简介:提高地震资料的分辨率是地震数据处理流程中的重要环节,对后续的精细构造解释起到重要作用。传统的提高分辨率方法大都假设地震资料是稳态的并且噪声水平不随空间发生变化,而实际情况不满足这一假设,导致提高分辨率处理后的效果达不到预期要求。针对这一问题,本文提出了一种基于时频二次谱的提高地震资料分辨率方法。首先,文中提出了基于S变换的时频二次谱,并结合模型论述了时变子波和反射系数在时频二次谱中的特征及其可分离性;其次,依据时变子波和反射系数在时频二次谱中的特征差异,构建了二维滤波器在地震记录的时频二次谱中提取时变子波的振幅谱;再次,文中研究了噪声环境中时变提高分辨率算子设计方法,并提出了依据时频谱能量强弱相对关系自适应确定频带拓宽范围的时变提高分辨率算子设计,进行提高分辨率处理;最后,文中对该方法进行了模型和实际数据的试处理,并与传统谱模拟方法和Q补偿方法的处理结果进行了对比分析,对比结果表明:本方法不需要估计Q值,提高分辨率能力不受震源子波频带的限制,在兼顾信噪比的前提下能够充分提高不同时间局部的地震数据的分辨率。
简介:全波形反演是一种高精度的地震成像方法,可以对地下介质物性参数模型进行准确的重构。然而在实际应用中,尤其是在三维复杂介质反演中,计算成本太大是该方法的一个重要缺陷。将混叠震源技术引入到频率域全波形反演中可以大幅度地降低计算成本,提高反演效率。但是使用震源编码技术也带来了两个问题:一方面,参与编码的各个震源之间会产生“串扰噪声”,导致反演结果中出现假象;另一方面,基于震源编码的频率域全波形反演方法周围噪声较为敏感,使该方法对含噪数据反演质量较差。本文引入一种频率组编码方法来压制“串扰噪声”,并基于震源编码技术提出一种频率域自适应全波形反演方法,通过一个与频率相关的自适应选择机制,将常规频率域全波形反演方法和基于震源编码的全波形反演方法联合起来,在保证反演质量的同时也最大程度地提高了反演效率。
简介:针对转换波各向异性叠前时间偏移实践中的四参数速度模型估计中参数间相互影响、难于准确的确定,本文将简化的转换波两参数动校正方程从叠加速度分析拓展到叠前时间偏移速度模型修正,形成了基于两参数动校正方程的新的四参数法转换波叠前时间偏移速度修正新方法和流程。在该方法和流程中,先进行转换波两参数叠加速度分析,获得转换波叠前时间偏移初始的速度和各向异性参数,然后通过分析偏移后的共成像点道集中的剩余时差修正速度和各向异性参数。其中叠前时间偏移速度模型的垂直速度比仍是需要利用纵波和转换波叠加剖面,层位标定后,利用相关的方法获得。有效速度比需结合纵波速度分析,利用Thomsen(1999)公式计算得到初始模型,最终依靠百分比扫描偏移处理获得最终的有效速度比模型。该方法简化了转换波高质量成像的速度估计方法,减小了多个参数估计的不确定性,在实际应用中也取得了较好的成像效果。
简介:无网格法形函数构造不依赖预定义的单元,具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法(EFGM)、点插值法(PIM)与径向基点插值法(RPIM)三种全域弱式无网格法的近似原理及特点;以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演,给出了RPIM形状参数的推荐值;分析了均匀介质模型大地电磁(MT)二维正演无网格法边界条件直接加载与罚函数法加载的精度差异,结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势,提出了EFG—PIM及EFG.RPIM耦合算法,数值计算结果验证了耦合算法的有效性。研究发现:无网格法及其耦合方法适用于电磁法数值模拟;支持域无量纲尺寸取1.0时无网格法精度与效率高,场节点与背景网格重合时计算效果佳;泊松方程求解PIM及RPIM精度较EFGM低,计算均匀介质MT响应精度较EFGM高;RPIM改善了PIM计算涉及的奇异性问题,对应支持域无量纲尺寸选择空间大。