简介:Theattenuationfactororqualityfactor(Q-factororQ)hasbeenusedtomeasuretheenergyattenuationofseismicwavespropagatinginundergroundmedia.ManymethodsareusedtoestimatetheQ-factor.WeproposeamethodtocalculatetheQ-factorbasedontheprestackQ-factorinversionandthegeneralizedS-transform.TheproposedmethodspecifiesastandardprimarywaveletandcalculatesthecumulativeQ-factors;then,itfindstheinterlaminarQ-factorsusingtherelationbetweenQandoffset(QVO)andtheDixformula.TheproposedmethodisalternativetomethodsthatcalculateinterlaminarQ-factorsafterhorizonpicking.Becausethefrequencyspectrumofeachhorizoncanbeextractedcontinuouslyona2Dtime–frequencyspectrum,themethodiscalledthecontinuousspectralratioslope(CSRS)method.ComparedwiththeotherQ-inversionmethods,themethodoffersnearlyeffortlesscomputationsandstability,andhasmathematicalandphysicalsignificance.WeusenumericalmodelingtoverifythefeasibilityofthemethodandapplyittorealdatafromanoilfieldinAhdeb,Iraq.TheresultssuggestthattheresolutionandspatialstabilityoftheQ-profileareoptimalandcontainabundantinterlaminarinformationthatisextremelyhelpfulinmakinglithologyandfluidpredictions.
简介:本文给出一种既能有效衰减地震噪音又可保护地层及构造的不连续性的新方法。构造约束保边平滑技术需要已知反射局部方位和边界信息,通常这些信息由全频率地震资料估算获得,但在资料信噪比很低的情况下,噪音往往会降低估算的可靠度。对于信噪比极低的地震资料,其主频成分相对非主频成分信噪比高,所以由主频资料获取的方位和边界信息比由其它频率成分获取的更可靠。方位和边界信息通常用倾角和相干值差异来描述。由于不同频率所引起的倾角和相干值差异的变化均比地震记录的变化缓慢,所以由主频资料获取的倾角及边界信息能够近似代表所有频率成分的倾角及边界信息。Ricker子波广泛用于地震勘探,Marr小波与Ricker子波在时间和频率域均具有相同的形态,所以选用Marrl小波变换将地震数据按照倍频程分为几个分频体。扫描主频分频体,用不等权二次曲面拟合并求解极大值来获取视倾角,通过比较9个滑动窗口的相干值来确定反射边界。将这些信息用构造约束保边平滑技术可选择性地(selectively)对主频、低频、高频分频体做平滑处理,最后将平滑后的各频段地震记录合成为滤波去噪后的地震记录。理论模型和实际资料处理效果表明该方法能有效压制噪音,保护边界,保护同相轴的连续性,且灵活地保留地震记录中的有用信息。
简介:针对矩形网射线追踪存在的模型剖分灵活性差、速度界面描述精度差等问题,研究了复杂结构三角网最小走时射线追踪全局算法。(1)根据剖分区域点、线、面的结构关系,遵循Delaunay三角剖分的优化准则进行三角网格剖分;(2)定义三角单元射线追踪的拓扑关系;(3)波源点及某一时刻波到达的每一个节点点构成波行面,在波行面扩展过程中计算节点的最小走时和次级源位置,实际次级源检索采用双曲线近似算法;(4)利用各节点走时和次级源方向信息,通过最小走时搜索,拾取从接收点到源点的射线路径。数值模拟结果表明,三角网射线追踪方法模型剖分时灵活性强、速度间断面的描述精度高,追踪结果准确。
简介:无网格法形函数构造不依赖预定义的单元,具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法(EFGM)、点插值法(PIM)与径向基点插值法(RPIM)三种全域弱式无网格法的近似原理及特点;以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演,给出了RPIM形状参数的推荐值;分析了均匀介质模型大地电磁(MT)二维正演无网格法边界条件直接加载与罚函数法加载的精度差异,结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势,提出了EFG—PIM及EFG.RPIM耦合算法,数值计算结果验证了耦合算法的有效性。研究发现:无网格法及其耦合方法适用于电磁法数值模拟;支持域无量纲尺寸取1.0时无网格法精度与效率高,场节点与背景网格重合时计算效果佳;泊松方程求解PIM及RPIM精度较EFGM低,计算均匀介质MT响应精度较EFGM高;RPIM改善了PIM计算涉及的奇异性问题,对应支持域无量纲尺寸选择空间大。
简介:为了有效处理人工源的影响,本文开发了带源的CSAMT二维正反演算法,可用于全区(近区、过渡区和远区)资料的反演。引入正则化因子完成磁法二维反演,并且将模型参数调整为磁化率的对数,保证反演过程中磁化率始终为正值。本文基于交叉梯度原理,将CSAMT和磁法进行联合反演,通过搜索交叉梯度项权重的方法,避免了不同异常源引起的两种异常相互干扰的问题。理论模型算例表明基于交叉梯度的联合反演方法优于单独反演。本文开发的带源CSAMT二维正反演算法,有效处理了人工源的影响,保证了最终联合反演算法的可靠性。
简介:针对不同矿物组分含量对盐膏岩速度影响大,盐膏岩速度复杂难题。本文把岩石物理模版的研究方法引入到盐膏岩的速度研究中,首次建立了盐膏岩四元组分(膏、盐、泥、孔隙地层水)含量与岩石速度关系的模版,在分析盐膏岩组分矿物骨架速度的基础上。基于多矿物测井解释模型,并以纵波模量比与孔隙度趋势作为约束条件,以组分含量为横坐标,岩石速度为纵坐标。在把不同尺度孔隙度值时,泥质、膏、盐矿物含量发生变化时,获得相应盐膏岩速度的变化值放在同一坐标系中,构成了盐膏岩四元矿物组分与岩石速度关系模版。并应用实测岩芯数据和测井数据,验证了该模版的可靠性和适应性,利用该模板能识别盐膏岩中关键结构和重要矿物的变化趋势,能直观的评估盐膏岩矿物组分蛮化对岩石涑席的影响.
简介:在具有垂直对称轴横向各向同性介质中,利用四种参数来确定中间至远偏移距转换波(C-波)动校正.它们是C-波叠加速度VC2,垂直速度比和有效速度比γ0和γeff,以及各向异性参数χeff.我们将这四种参数作为C波叠加速度模型.C-波速度分析的目的就是确定这种叠加速度模型.C-波叠加速度模型VC2,γ0,γeff,和χeff可以由P-波和C-波反射动校正资料获得.然而错误的传播是C-波反射动校正反演中的严重问题.当前短排列叠加速度由于是从双曲线动校正推算而得,因而其精度不足以为各向异性参数提供有意义的反演值.中间偏移非双曲线动校正不再被人们所勿略,而是可以用一个背景γ加以量化.非双曲线分析通过中间偏移距的γ校正量可以产生VC2,若数据不含燥音,其误差小于1%.方法稳健,允许γ启始假定值的误差达20%.该方法也适用垂直非均匀各向异性介质.精度的提高使能够用4分量地震资料计算各向异性参数.为此提出了两种工作流程:双扫描和单扫描流程.理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效.
简介:本文在Oristaglio等(1984)和Adhidjaja等(1985)工作基础上,给出线源二维时间域瞬变电磁二次场的DuFort-Frankel有限差分数值解,有效避免了在总场求解法中场源附近的奇异问题,并对地-空边界电导率的处理、归一化感应电动势偏导数的计算、推进时间步的确定,提出了改进方法;吸取前人成就中二次场地-空边界向上延拓和零值边界处理技术,从而简化了计算方法;通过对均匀大地、水平层状大地模型的计算,二次场求解法与解析法的最大相对误差小于0.01%,计算速度比总场求解法提高了约3倍;模拟计算不同时刻瞬变电磁场在地下的分布形态,描绘出感应涡流向下向外的传播特征,以及与地下异常体相互作用的物理过程。
简介:以阿姆河盆地右岸地区碳酸盐岩台内滩气田为例,运用基于图论的多分辨率聚类分析法开展了以常规测井资料为基础的测井相聚类及岩相定量预测研究。该方法不需要分析数据体的结构及聚类数等先验知识为指导,能够自动优选最佳聚类个数,并允许按照实际需求控制聚类级别,进行聚类合并。依据岩芯资料岩相识别及测井相标定结果,本文最终建立了包含5个聚类的测井相划分模型及岩相定量解释图版,其中,聚类测井相1~5分别对应于泻湖泥、石膏坪、滩间、低能滩及高能滩,取芯段符合率达85%以上,能够较好的运用于非取芯段岩相预测研究。据此,我们进行了连续的聚类测井相划分及岩相预测,并对层序地层格架内岩相分布及物性特征进行了分析。
简介:有限频走时层析成像是近年发展起来的一种新方法,这种新方法的一个主要过程是走时灵敏度算核的计算。求解灵敏度算核要多次用到同一散射点的走时,多次地求解同一走时是相当耗时的任务,如果介质为均匀或速度线性变化等简单模型,散射点的走时可以用解析公式快速地求出,从而灵敏度算核的计算耗时相对较少。然而各种地球模型中,介质速度大多为分层模型,从解析公式中得到走时信息就比较繁锁。为了提高计算效率,本文采用查表算法研究地球分层速度模型中的P震相有限频走时灵敏度算核的计算,选用的速度模型是地球AKl35模型,用查表算法求解走时,节省了约50%的计算时间。在相同的速度模型下,与已有结果的对比,本文所用的查表算法,能在兼顾精度的前提下,以较小的存储要求换取较高的计算效率,这对提高有限频走时层析成像算法的速度具有一定的参考价值。
简介:储层物性参数的改变伴随着弹性参数某种程度的改变,而这种改变并非简单的线性关系。加上观测信息的缺失、重叠、噪声破坏及模型理想化等原因,致使这种改变关系的获取具有较大的不确定性。针对传统统计岩石物理物性参数反演方法的不足,以利用弹性参数反演储层物性参数为目的,依据贝叶斯反演框架,我们建立了新的储层物性参数目标反演函数。首先,采用兼具确定性与随机性特点的统计岩石物理模型,考虑到不同弹性参数间的精度存在差异,引入权重系数,建立起储层物性参数与弹性参数间的加权统计关系。其次,基于这种加权统计关系,结合马尔科夫链蒙特卡洛随机模拟技术产生储层物性参数、弹性参数随机联合样本空间作为目标函数求解样本空间。最后,建立解的快速求解准则,求取最大后验概率密度对应的储层物性参数取值作为最终解。实际应用表明,该方法具有较高的反演效率,应用前景较好。
简介:中国西部复杂山地是地震勘探的难点区域,地表和地下都复杂,从而导致地表激发困难,地震资料干扰严重、信噪比低。地下构造形态复杂,导致反射点分散严重,前人尝试过弯线采集,但很少考虑地下实际构造特征。本文章提出基于CRP的不规则采集设计方法,使地表避开难以激发的区域,地下考虑构造特征及CRP点分布位置,优化地表接收线位置,通过地表有目的布置接收点位置,最大程度的保证地下CRP点不分散,从而提高地震数据的信噪比.并以四川盆地某工地的实际资料通过模型数据验证了方法的合理性与可行性。表明该方法能够解决复杂地区地震采集存在的部分难题,为复杂地区的地震勘探采集提供一个新的思路。
简介:目前在研究地球物理参数的空间变异特性时常常是基于几何各向异性假设,然而实际中几何各向异性假设常常不成立。本文提出了一种求解地球物理参数各向异性的新方法,首先对各方向上的变异值作归一化处理,然后用拟合的方法得到各方向单位步长上的变异性强度,进而实现对全局各向异性的求解。以上海地区合成孔径雷达干涉测量(InSAR)中差分大气延迟样本数据为例,验证了新方法求解各向异性的优越性:新方法的偏差只有6.4%,而基于几何各向异性假设方法的偏差达到了21.2%,而且新方法的运算效率有非常显著的提高。进一步,通过克里金插值交叉验证实验,证明了基于新算法得到的结构函数的插值效果最好,从而从另外一个角度验证了新算法求解的各向异性更准确,更好地表征了区域化变量的空间结构特征。因此,与几何各向异性假设的方法相比,新方法不仅能更准确的描述地球物理参数的各向异性,而且有更高的运算效率,为更准确地估计出所需的地球物理参数奠定了基础。
简介:中国西部复杂山地是地震勘探的难点区域,地表和地下都复杂,从而导致地表激发困难,地震资料干扰严重、信噪比低。地下构造形态复杂,导致反射点分散严重,前人尝试过弯线采集,但很少考虑地下实际构造特征。本文章提出基于CRP的不规则采集设计方法,使地表避开难以激发的区域,地下考虑构造特征及CRP点分布位置,优化地表接收线位置,通过地表有目的布置接收点位置,最大程度的保证地下CRP点不分散,从而提高地震数据的信噪比.并以四川盆地某工地的实际资料通过模型数据验证了方法的合理性与可行性。表明该方法能够解决复杂地区地震采集存在的部分难题,为复杂地区的地震勘探采集提供一个新的思路。
简介:如何快速、精确地利用叠前深度偏移进行偏移速度分析是勘探地震学的一项重要研究内容,针对该问题,本文提出一种二阶精度广义非线性全局最优的偏移速度反演方法。我们将首先去掉速度模型修正量与成象深度差呈线性关系的假设,推导出具有二阶精度的速度模型修正量计算公式,使每一次迭代得到的速度模型尽可能地接近实际模型;然后采用广义非线性反演方法反演获得对所有道集的全局最优的速度模型修正量,不仅极大地加快了收敛速度,而且反演过程中陷入局部极小的可能性也减小了。理论模型和Marmousi模型的处理结果表明:本方法精度高、处理速度快,提高了偏移速度分析方法的实用性和对复杂构造成像的准确性。
简介:提高地震资料的分辨率是地震数据处理流程中的重要环节,对后续的精细构造解释起到重要作用。传统的提高分辨率方法大都假设地震资料是稳态的并且噪声水平不随空间发生变化,而实际情况不满足这一假设,导致提高分辨率处理后的效果达不到预期要求。针对这一问题,本文提出了一种基于时频二次谱的提高地震资料分辨率方法。首先,文中提出了基于S变换的时频二次谱,并结合模型论述了时变子波和反射系数在时频二次谱中的特征及其可分离性;其次,依据时变子波和反射系数在时频二次谱中的特征差异,构建了二维滤波器在地震记录的时频二次谱中提取时变子波的振幅谱;再次,文中研究了噪声环境中时变提高分辨率算子设计方法,并提出了依据时频谱能量强弱相对关系自适应确定频带拓宽范围的时变提高分辨率算子设计,进行提高分辨率处理;最后,文中对该方法进行了模型和实际数据的试处理,并与传统谱模拟方法和Q补偿方法的处理结果进行了对比分析,对比结果表明:本方法不需要估计Q值,提高分辨率能力不受震源子波频带的限制,在兼顾信噪比的前提下能够充分提高不同时间局部的地震数据的分辨率。
简介:转换波技术是解决气云区成像问题的有效手段,而二极化现象则是制约转换波构造成像精度提高的一个重要因素。解决该问题的一个有效办法是二极化校正:将转换波速度分解为基本速度和速度扰动,分两步对转换波地震数据进行叠加速度分析和叠前时间偏移速度分析。柴达木盆地三湖地区天然气资源丰富,气云广泛发育,通过对三湖地区二维三分量地震数据的研究分析,发现该地区转换波二极化现象明显,对转换波的准确成像有较大影响。通过对转换波二极化校正方法的应用,不论是在叠加成像,还是叠前偏移成像都取得了较好的效果,提高了气云区转换波构造成像的精度。二极化校正对解决气云区转换波剩余静校正有重要促进作用,二极化校正与转换波剩余静校正的结合是陆上气云区转换波高精度成像流程中必不可少的环节。
简介:频率域全波形反演充分利用全波场的振幅、相位以及频率信息,采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础,对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中,针对截断边界问题,并考虑多频率联合反演中计算区域采用同一套剖分网格的需求,提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子,通过一定方式调节衰减因子使得声波在衰减层中充分衰减,达到压制截断边界影响的目的。根据指数衰减规律,文中推导出了一种新的衰减因子计算公式,并给出了不同频率条件下衰减层厚度计算公式;在反演算法中,采用共轭梯度法求解高斯牛顿反演迭代方程组,避免直接求解雅克比矩阵和Hessian矩阵带来的巨额计算量,并采用相同的反演模型,对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明,本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。