简介：频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播，但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性，本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验，验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理，基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像，说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力，我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理，对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力，具有一定实际应用价值。

简介：随着并行计算技术的发展,非线性反演计算效率在不断提高,但对于基于单点搜索的非线性反演方法,其并行算法的实现则是一个难题。本文将群体搜索的思想引入到基于单点搜索的非线性反演方法,构建了并行算法,以量子蒙特卡罗方法为例进行了二维地震波速度反演及实际资料波阻抗反演,并测试了使用不同节点数进行计算的效率。计算结果表明：该并行算法在理论和实际资料反演中是可行的和有效的,具有很好的通用性;算法计算效率随着使用节点数的增加而提高,但算法计算效率的提高幅度随着使用节点数的增加逐渐减小。

简介：为了进一步提高大地电磁三维快速松弛反演的计算效率，在深入分析大地电磁三维快速松弛反演算法的基础上，结合MPI自身的优越性，确定了并行计算的思路，实现了三维快速松弛反演的并行计算。通过理论模型合成数据和实测数据对实现的三维快速松弛反演并行程序进行了试算，分析对比了在多种情况下程序的执行效率。测试结果表明，所实现的三维快速松弛反演并行程序运行结果正确，效率提高明显。此成果在普通微机上实现，推动了大地电磁三维反演技术的实用化，可为其它地球物理三维正反演研究所借鉴。

简介：针对矩形网射线追踪存在的模型剖分灵活性差、速度界面描述精度差等问题,研究了复杂结构三角网最小走时射线追踪全局算法。（1）根据剖分区域点、线、面的结构关系,遵循Delaunay三角剖分的优化准则进行三角网格剖分;（2）定义三角单元射线追踪的拓扑关系;（3）波源点及某一时刻波到达的每一个节点点构成波行面,在波行面扩展过程中计算节点的最小走时和次级源位置,实际次级源检索采用双曲线近似算法;（4）利用各节点走时和次级源方向信息,通过最小走时搜索,拾取从接收点到源点的射线路径。数值模拟结果表明,三角网射线追踪方法模型剖分时灵活性强、速度间断面的描述精度高,追踪结果准确。

简介：正则化因子的选择方法决定了正则化方法的稳定性和准确度。在分析改进的奇异值分解法与正则化方法的关系的基础上，给出一个正则因子计算公式。数值模拟试验表明改进的奇异值分解法和正则化方法适应低信噪比。当信噪比低于30时正则化方法比奇异值分解法结果更好，当信噪比高于30时奇异值分解法更优。采用本文提出的正则因子的正则化方法可以适用于实际核磁测井的低信噪比（5〈SNR）情况。数值模拟及实际资料的计算实验表明，改进的SVD反演算法与正则化方法都具有适应低信噪比、求解速度快等优点，可以满足核磁共振测井工作的需求。

简介：基于曾新吾和MacBeth提出的时域内双源累积旋转方法,建立了频率相关介质中分析多分量VSP数据的横波分裂参数提取算法（DCTF）。该算法可以在频域中针对单个频率提取横波各向异性参数（快横波的极化方向以及快、慢横波间的时间差）,从而避免了目前常用方法中使用窄带通滤波可能带来的误差。通过对地震合成记录的数值分析,确定了该算法的可行性和正确性,并与目前常用方法的应用结果进行了比较。结果显示,频率相关横波分裂参数可以利用DCTF从地震四分量数据中直接提取。在地震频率范围内,含较大尺度裂缝时各向异性参数将表现出频率相关性,这意味着在地震频率范围会出现频散。随着频率的增加,各向异性有降低的趋势。

简介：随着全张量重力梯度（FTG）测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU（GraphicsProcessingUnit）和预处理分解技术（Preconditioningmethods）在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛（SSOR）技术与不完全乔列斯基分解共轭梯度算法（ICCG）提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATeslaC2050GPU的并行SSOR-ICCG算法和在2.0GHzCPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。

简介：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.

简介：致密砂岩具有一定的各向异性。本文通过不同方向样品的声速测量及相关配套实验对致密砂岩进行了孔隙结构及各向异性分析，同时与页岩进行了对比。结果表明，致密砂岩样品所表现出的各向异性主要是由沉积环境差异形成的成分层、薄互层造成，具有典型的TI介质特征。不同方向杨氏模量随围压增加而增大，泊松比变化不明显。各向异性系数均随有效压力增加而减少，同时ε与γ和δ之间存在一定的线性关系。最后总结了不同地区的各向异性规律，为不同地区选择适合的评价模型奠定了基础。本次研究可为致密砂岩的测井评价、地震资料解释及压裂开发提供实验依据。

简介：地震低频信息能够提高分辨率与成像精度，改善反演质量，甚至直接进行油气检测，需要对其进行有效保护与拓展。对于子波而言，缺失低频信息会导致主瓣幅度降低、第一旁瓣幅度增加，并出现次级旁瓣呈周期震荡衰减的现象；从合成地震记录和典型地质模型来看，低频缺失会产生假同相轴，造成分辨率提高的假象，且模型不同位置的特征存在一定差异；对缺失低频的模型数据进行波阻抗反演，会造成构造失真、岩性改变的假象，特别是高陡构造和薄互层。针对缺失低频的地震资料，本文还研究了基于压缩感知与稀疏约束的拓频方法，开发了相应的模块，并对实际CIP道集进行处理，取得了较好的应用效果。

简介：针对转换波各向异性叠前时间偏移实践中的四参数速度模型估计中参数间相互影响、难于准确的确定,本文将简化的转换波两参数动校正方程从叠加速度分析拓展到叠前时间偏移速度模型修正,形成了基于两参数动校正方程的新的四参数法转换波叠前时间偏移速度修正新方法和流程。在该方法和流程中,先进行转换波两参数叠加速度分析,获得转换波叠前时间偏移初始的速度和各向异性参数,然后通过分析偏移后的共成像点道集中的剩余时差修正速度和各向异性参数。其中叠前时间偏移速度模型的垂直速度比仍是需要利用纵波和转换波叠加剖面,层位标定后,利用相关的方法获得。有效速度比需结合纵波速度分析,利用Thomsen（1999）公式计算得到初始模型,最终依靠百分比扫描偏移处理获得最终的有效速度比模型。该方法简化了转换波高质量成像的速度估计方法,减小了多个参数估计的不确定性,在实际应用中也取得了较好的成像效果。

简介：大型、构造复杂的三维物理模型可用于模拟油气勘探。构造逼近实际地质状况的模拟具有制作技术难度大、质量控制严格等特点,可用于采集宽方位、多方位和全方位的地震数据,从而进行多种三维处理、解释方法验证。本文针对中国西部前陆盆地地表条件复杂地下构造复杂,导致成像不理想等问题,基于复杂的地下构造,设计制作了目前世界上模拟施工面积最大、构造最复杂的KS（塔里木盆地克深勘探工区）物理模型。本文的模型技术的进步主要涉及3个方面：模型的设计方法、模型的浇铸流程和数据采集,首次给出了物理模型的三维真实速度模型,定量分析了物理模型的制作精度,绝对误差小于3mm,可以满足方法试验的需要。该模型基于三维形态测量技术建立了三维真实速度模型,可作为方法试验的基础数据。因此,该模型可作为地震物理模拟技术的标准。

简介：地震反演是当今最广泛应用于含油气储层预测的技术之一，取得了很多很好的预测效果，但也有失败的例子，达不到区分岩性和识别流体的目的。而本文介绍的建立在岩石物理建模和分析基础上的地震弹性反演，可将含油气储层预测由定性向（半）定量推进一步。根据岩石物理建模和正演扰动分析，可深刻理解岩石物性参数与地震弹性参数之间的内在关系，进而寻找储层岩性、物性和含油气性的敏感地震弹性参数，建立起理论岩石物理解释图版。岩石物理分析结果和所建立的岩石物理解释图版分别用以指导地震反演和反演结果的解释，实现油气储层分布预测和流体检测的目的。文中的含气砂岩分布预测实例研究应用结果表明，这种方法较叠后地震反演储层预测技术具有无可比拟的优越性，效果更佳，效率更高。

简介：地震纹理属性与地震相和储层特征密切相关,广泛地应用在地震资料的解释中。传统的地震纹理属性基于叠后数据提取,受叠加作用的影响,易造成地层特征信息的损失,不利于复杂储层的描述。本文提出叠前纹理属性,其不仅可以精细地刻画不同反射点波形的横向连续性,也能体现AVO、各向异性和介质的均质性。基于叠前纹理属性丰富的地层特征表达能力,结合SOM聚类算法,形成了利用叠前数据进行地震相分析的方法。该方法应用于中国某工区宽方位地震资料,通过对比证实了叠前纹理属性描述地层横向变化的优越性,并能揭示各向异性特征及非均质性特征,基于叠前纹理的分类结果能有效区分不同的地震反射模式,为地震相分析提供了可靠的依据。

简介：基于岩石物理模型和广义各向异性Zoeppritz方法在频率域计算裂缝型反射体模型反射波方位地震AVO响应。反射体模型为粘弹各向异性、有限厚度的地质体，其地震反射波形序列包含如下动力学信息，即分界面处介质的波阻抗和非弹性差异、反射体内部波的各向异性传播、在传播路径上的频散与衰减，以及来自顶底界面的反射波的调谐与干涉等。计算表明，速度频散和衰减增顶界面反射波大入射角反射时的振幅，而减弱底界面反射振幅。对于固定人射角的方位地震方位地震响应，PP波反射特征表现为随方位角的增加反射波形序列延续时间变长，而PSV和PSH转换类型反射波的方位各向异性变化特征稳定且受储层厚度影响较小，表现为PSV波反射振幅随方位角增加而增加，PSH波在0°和90°方位无反射能量，在45°方位反射振幅最强。