简介:传统上,有限差分的差分系数一般可以通过泰勒级数展开法或优化方法来极小化频散误差得到。基于泰勒级数展开的差分法在有限的波数范围内精度较高,但在这个范围之外会产生较强的数值频散;基于最小二乘的优化有限差分法能在更大的波数范围内达到较高的精度,并可以在较小的计算需求内获得全局最优解。本文将基于最小二乘的优化有限差分法从二维正演模拟推广到三维,形成了计算效率高、高精度范围宽、适合并行计算的三维声波优化有限差分方法。频散分析及正演模拟表明本文发展的有限差分方法可以很好地压制数值频散。最后,将本文发展的有限差分方法应用到三维逆时偏移的震源波场延拓和检波点波场延拓中,并结合有效边界存储策略与checkpointing技术在GPU集群上实现三维逆时偏移以提高计算效率、减少存储量。三维逆时偏移试算结果表明本文三维优化有限差分方法与传统的有限差分法相比可以获得更高精度的偏移成像结果。
简介:频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播,但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性,本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验,验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理,基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像,说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力,我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理,对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力,具有一定实际应用价值。
简介:常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。
简介:频率空间域地震波数值模拟具有独特的优势:可以同时模拟多源的波传播、每个频率之间独立并行地计算、计算频带选择灵活、不存在累计误差、容易模拟粘弹性介质中地震波传播。但是该方法的最大瓶颈是对于计算机内存的巨大需求。我们使用压缩存储系数矩阵的方法,极大地减少了计算机内存的需求量。同时为了减少短差分算子的数值频散,引用了频率空间域25点弹性波波动方程的差分格式,并使用了最小二乘意义下求出的优化差分系数。为了克服边界反射,采用了最佳匹配层吸收边界条件。数值模拟试验证明:用压缩存储系数矩阵及优化差分系数的频率空间域25点差分格式进行弹性波正演模拟,可以减少数值频散,提高计算精度。使用较大的网格间距,降低计算机内存需求,并保持较高的计算效率。该正演方法为后续弹性波偏移和弹性参数反演提供较好的基础。
简介:High-qualityseismicgeometryisthekeytoobtainhigh-qualityseismicdata,andcanaffecttheaccuracyofdataprocessingandimaging.Basedontheanalysisoftherelationshipbetweenthequalityofthegeometryandthefouracquisitionparameters(thenumberoftraces,shotlinespacing,andthespaceandnumberofreceiverlines),aqualityevaluationmethodofthegeometrybasedoncomprehensivequalityfactor(CQF)isproposed,andtherelationshipbetweenthegeometryqualityandthefourparametersisgiven.WeusefielddatacollectedinanoilfieldinWesternChinawithcomplexgeology:Firstweuseawideazimuthgeometry.Then,wecalculatetherelationshipcurvebetweengeometryanddataqualitybyvaryingeachparameterwhilekeepingtherestfixed.andtheanalysisresultsaregivenbyusingtheCQFevaluationmethod.Theresultsshowthattheshot-linespacinghasthegreatesteffectonthequalityofthegeometry,andtheincreaseofthereceiverlinespacingcanappropriatelyimprovethequalityofthegeometry,andtheincreaseofthenumberofreceivingtracescanimprovethegeometryquality.Thedifferentacquisitionparametershavedifferenteffectsontheimagingqualityofshallowanddeepevents.Themodelforwardandprestackdepthmigrationareusedtogenerateprestackdepthmigrationprofileswithdifferentacquisitionparameters.Theimagingresultsareconsistentwiththeabovecalculatedresults.Accordingtothedepthofthetargetlayer,thequalityfactorevaluationmethodisappliedtoguidethedesignofthegeometryandoptimizetheacquisitionparameterstoimprovetheimagingaccuracyofseismicdata.
简介:在文献中,Kirchhoff型反偏移场的稳相分析主要是在下列两个条件下进行的:(1)等时面和目标反射面相切;(2)深度偏移像场信号的长度接近于零。对于与目标反射面不相切的等时面和长度远大于零的深度偏移像场子波,已有的结果将不再成立。为了在等时面和目标反射面不相切和深度偏移像场的子波长度远大于零的条件(一般条件)下对Kirchhoff型反偏移场进行稳相分析,我推导了出现在二维稳相分析公式中的诸因子的计算公式,并从中发现:(1)对于不同的等时面,距离差函数的稳相点具有不同的水平坐标;(2)Kirchhoff型真振幅反偏移的输出场由两部分(真振幅反偏移信号与振幅畸变因子)的乘积组成。由此得到下列两个结论:(1)一个给定的反偏移信号由多个深度偏移信号上的采样点组装而成,反偏移信号上的采样点个数等于对于这种组装有贡献的偏移信号的个数。(2)振幅畸变效应是Kirchhoff型反偏移中的固有效应,靠反偏移本身无法消除。如果一定要消除这种振幅畸变效应,必须对反偏移结果进行振幅校正。
简介:野外空间采样密度的提高将增加室内数据分析的工作量,常规的基于地震数据的点与线分析方法有一定的局限性。本文简要说明空间子集的抽取方法,列举了正交子集和斜交子集的特点,并通过三维可视化展示了子集数据的空间特性。提出在数据处理中利用子集的时间切片:(1)分析地震道空间分布的均匀性和规则性;(2)研究面波及规则干扰的空间分布特点;(3)检测叠前数据中的异常信息;(4)监控叠前去噪的效果。实际数据的应用结果表明基于空间子集的分析方法是一种独特有效的地震数据分析方法,从另一种视角观察地震数据,可以发现其中某些新的特征,以提高处理人员对数据的洞察能力。
简介:裂缝识别对碳酸盐岩储层测井评价至关重要,传统的测井仪器由于探测深度浅(小于3米)而无法对井壁外围(大于3米)裂缝发育情况进行评价,远探测声波测井仪器采用相控阵发射、同时加大源距,有效提升了测井仪器的探测深度。但由于缺少正演模拟研究,目前对于裂缝的解释往往是基于经验而缺乏理论依据,导致很多现象难以解释。本文利用高阶有限差分方法对远探测声反射波测井裂缝识别进行了正演模拟及叠前逆时偏移成像研究,首先在理论研究的基础上构建正演模型,重点研究不同裂缝离井壁距离、裂缝张开度和倾角的响应特征;其次在单因素变化基础上提取成像区域的能量强度,分析确定出在实际地层速度有波动变化时远探测声反射波测井方法能够识别裂缝的条件;最后通过对裂缝识别的影响因素定量化分析,确定裂缝识别的最大距离、最小张开度和最小倾角,降低了裂缝识别与评价中的多解性。研究成果对远探测声反射波测井仪器的发展、数据处理方法的改进,以及后续的测井解释工作都有一定的借鉴意义。
简介:为了探明介质的非完全弹性对瑞利型槽波传播的影响以及利用瑞利型槽波品质因子QR反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽波频散方程,计算了煤层瑞利型槽波的品质因子QR曲线;采用控制变量法,分析了瑞利型槽波品质因子QR随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽波品质因子QR曲线与群速度曲线变化趋势一致;埃里相(Airy-phase)频率附近瑞利型槽波的品质因子QR最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽波品质因子QR随煤层横波品质因子QS2的增大而增大。利用瑞利型槽波品质因子QR曲线可以实现对瑞利型槽波的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。更多还原
简介:ThevariationsintheEarth’srotationareimportanttospacedynamictheoryandnaturaldisastersbecausetheyaffectthelength-of-day(LOD)andconsequentlyhumanlives.Weusemaximumentropymethod(MEM),Lombmethod(LOMB),andphasedispersionminimization(PDM)todeterminethenaturalperiodsofequallyspacedLODtimeseries.WetransformtheobservationalmonthlyLODtimeseries(LODM)tounequallysampledseries(LODMD)byremovingeveryfourth,third,andhalfofthetotalsamples.WealsoapplysplineinterpolationtoLODMDtoyieldequallyspacedtimeseries(LODMDN).Theresultssuggestthatregardlessofthetimeseries,theMEMfrequencyis0.1660month^-1and0.0840month^-1,whereasLOMBandPDMyield0.166month^-1and0.083month^-1,respectively.Furthermore,missingdatathatarelessthanhalfofthetotaldataorsplineinterpolationdonotaffecttheanalysis.Fortheamplitude,neithermissingdatanorsplineinterpolationaffecttheanalysis.
简介:地震低频信息能够提高分辨率与成像精度,改善反演质量,甚至直接进行油气检测,需要对其进行有效保护与拓展。对于子波而言,缺失低频信息会导致主瓣幅度降低、第一旁瓣幅度增加,并出现次级旁瓣呈周期震荡衰减的现象;从合成地震记录和典型地质模型来看,低频缺失会产生假同相轴,造成分辨率提高的假象,且模型不同位置的特征存在一定差异;对缺失低频的模型数据进行波阻抗反演,会造成构造失真、岩性改变的假象,特别是高陡构造和薄互层。针对缺失低频的地震资料,本文还研究了基于压缩感知与稀疏约束的拓频方法,开发了相应的模块,并对实际CIP道集进行处理,取得了较好的应用效果。
简介:针对转换波各向异性叠前时间偏移实践中的四参数速度模型估计中参数间相互影响、难于准确的确定,本文将简化的转换波两参数动校正方程从叠加速度分析拓展到叠前时间偏移速度模型修正,形成了基于两参数动校正方程的新的四参数法转换波叠前时间偏移速度修正新方法和流程。在该方法和流程中,先进行转换波两参数叠加速度分析,获得转换波叠前时间偏移初始的速度和各向异性参数,然后通过分析偏移后的共成像点道集中的剩余时差修正速度和各向异性参数。其中叠前时间偏移速度模型的垂直速度比仍是需要利用纵波和转换波叠加剖面,层位标定后,利用相关的方法获得。有效速度比需结合纵波速度分析,利用Thomsen(1999)公式计算得到初始模型,最终依靠百分比扫描偏移处理获得最终的有效速度比模型。该方法简化了转换波高质量成像的速度估计方法,减小了多个参数估计的不确定性,在实际应用中也取得了较好的成像效果。
简介:大型、构造复杂的三维物理模型可用于模拟油气勘探。构造逼近实际地质状况的模拟具有制作技术难度大、质量控制严格等特点,可用于采集宽方位、多方位和全方位的地震数据,从而进行多种三维处理、解释方法验证。本文针对中国西部前陆盆地地表条件复杂地下构造复杂,导致成像不理想等问题,基于复杂的地下构造,设计制作了目前世界上模拟施工面积最大、构造最复杂的KS(塔里木盆地克深勘探工区)物理模型。本文的模型技术的进步主要涉及3个方面:模型的设计方法、模型的浇铸流程和数据采集,首次给出了物理模型的三维真实速度模型,定量分析了物理模型的制作精度,绝对误差小于3mm,可以满足方法试验的需要。该模型基于三维形态测量技术建立了三维真实速度模型,可作为方法试验的基础数据。因此,该模型可作为地震物理模拟技术的标准。
简介:地震反演是当今最广泛应用于含油气储层预测的技术之一,取得了很多很好的预测效果,但也有失败的例子,达不到区分岩性和识别流体的目的。而本文介绍的建立在岩石物理建模和分析基础上的地震弹性反演,可将含油气储层预测由定性向(半)定量推进一步。根据岩石物理建模和正演扰动分析,可深刻理解岩石物性参数与地震弹性参数之间的内在关系,进而寻找储层岩性、物性和含油气性的敏感地震弹性参数,建立起理论岩石物理解释图版。岩石物理分析结果和所建立的岩石物理解释图版分别用以指导地震反演和反演结果的解释,实现油气储层分布预测和流体检测的目的。文中的含气砂岩分布预测实例研究应用结果表明,这种方法较叠后地震反演储层预测技术具有无可比拟的优越性,效果更佳,效率更高。
简介:地震纹理属性与地震相和储层特征密切相关,广泛地应用在地震资料的解释中。传统的地震纹理属性基于叠后数据提取,受叠加作用的影响,易造成地层特征信息的损失,不利于复杂储层的描述。本文提出叠前纹理属性,其不仅可以精细地刻画不同反射点波形的横向连续性,也能体现AVO、各向异性和介质的均质性。基于叠前纹理属性丰富的地层特征表达能力,结合SOM聚类算法,形成了利用叠前数据进行地震相分析的方法。该方法应用于中国某工区宽方位地震资料,通过对比证实了叠前纹理属性描述地层横向变化的优越性,并能揭示各向异性特征及非均质性特征,基于叠前纹理的分类结果能有效区分不同的地震反射模式,为地震相分析提供了可靠的依据。
简介:基于岩石物理模型和广义各向异性Zoeppritz方法在频率域计算裂缝型反射体模型反射波方位地震AVO响应。反射体模型为粘弹各向异性、有限厚度的地质体,其地震反射波形序列包含如下动力学信息,即分界面处介质的波阻抗和非弹性差异、反射体内部波的各向异性传播、在传播路径上的频散与衰减,以及来自顶底界面的反射波的调谐与干涉等。计算表明,速度频散和衰减增顶界面反射波大入射角反射时的振幅,而减弱底界面反射振幅。对于固定人射角的方位地震方位地震响应,PP波反射特征表现为随方位角的增加反射波形序列延续时间变长,而PSV和PSH转换类型反射波的方位各向异性变化特征稳定且受储层厚度影响较小,表现为PSV波反射振幅随方位角增加而增加,PSH波在0°和90°方位无反射能量,在45°方位反射振幅最强。
简介:在水库现场试验了RISK2型探地雷达探测水库冰厚度的能力,试验时所用天线频率为600MHz;同步钻孔测量雷达探测处的冰厚度;以及在一个点上取样测试分析冰晶体、冰内气泡和冰密度。试验时冰面积雪厚度0.03-0.05m,冰层上部有0.24m粒状冰,其下均为柱状冰;冰内气泡含量呈表层高底层低分布;冰密度随气泡含量变化;冰厚度在平面内不均一。通过探测厚度和实测厚度的对比分析以及气泡含量对介电系数影响的理论分析,建立了积雪、粒状冰和柱状冰三层介质模型,获取雷达波在冰内的理论传递时间。结果发现:能够利用等效介电常数或等效传播速度评价雷达波传递时间,结冰期冰层1/3深度处的对应介电常数或传递速度可以作为等效值;另外因冰内大气泡造成的理论传递时间大于雷达探测时间,其差值随理论传递时间或冰厚的增加呈非线性增加。
简介:基于斑块饱和模型计算渗透率变化的地震反射特征,为流体流动性的地震描述提供依据。基于传播矩阵理论设计反射系数与合成记录算法,实现了频率域岩石物理模型与地震响应计算的无缝连接。斑块饱和储层地震响应包含如下动力学信息:分界面处波阻抗差异、储层内部波的频散与衰减,以及顶底界面波的调谐与干涉。模拟结果表明,渗透率的增加显著降低纵波速度,使其在高、低频弹性极限之间发生频散。储层速度频散与层状构造共同导致反射系数的频变现象。在储层与围岩波阻抗接近的情况下,地震响应对渗透率变化具有敏感性,对于不同储层厚度,当围岩为高速页岩时,反射波叠加振幅随渗透率增加而增加;当围岩为低速页岩时,叠加振幅随渗透率增加而降低。