简介：传统Robinson褶积模型主要受缚于三种不合理的假设,即白噪反射系数、最小相位地震子波与稳态假设,而现代反射系数反演方法（如稀疏约束反褶积等）均在前两个假设上寻求突破的同时却忽视了一个重要事实：实际地震信号具有典型的非稳态特征,这直接冲击着反射系数反演中地震子波不随时间变化的这一基础性假设。本文首先通过实际反射系数测试证实,非稳态效应造成重要信息无法得到有效展现,且对深层影响尤为严重。为校正非稳态影响,本文从描述非稳态方面具有普适性的非稳态褶积模型出发,借助对数域的衰减曲线指导检测非稳态影响并以此实现对非稳态均衡与校正。与常规不同,本文利用对数域Gabor反褶积仅移除非稳态影响,而将分离震源子波和反射系数的任务交给具有更符合实际条件的稀疏约束反褶积处理,因此结合两种反褶积技术即可有效解决非稳态特征影响,又能避免反射系数和地震子波理想化假设的不利影响。海上地震资料的应用实际表明,校正非稳态影响有助于恢复更丰富的反射系数信息,使得与地质沉积和构造相关的细节特征得到更加清晰的展现。

简介：非局部均值滤波是一种基于图像信息冗余的去噪方法，其认为图像自身的有效结构具有一定的重复性，而随机噪声则不具备这一特点，通过利用图像本身的自相似性来达到压制随机噪声的目的，是一种全局的去噪方法。本文把这一思想引入地震数据随机噪声压制中，针对传统非局部均值滤波计算量过大的问题，文章采用分块非局部均值的方式来减少计算量；针对滤波参数选取会影响非局部均值滤波效果的问题，提出一种简单的自适应滤波参数地震数据分块非局部均值算法。模型和实际数据处理结果表明：相对于传统的去噪算法（如f-x反褶积），该方法在压制随机噪声的同时对有效信号保护地更好，具有更高的保真度，更有利于后续的处理和解释工作。

简介：常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI（ATI）条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的＂精确走时＂结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。

简介：通过最新的全球地磁模型——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）处衰减较快,在地表处则趋于稳定。

简介：以阿姆河盆地右岸地区碳酸盐岩台内滩气田为例,运用基于图论的多分辨率聚类分析法开展了以常规测井资料为基础的测井相聚类及岩相定量预测研究。该方法不需要分析数据体的结构及聚类数等先验知识为指导,能够自动优选最佳聚类个数,并允许按照实际需求控制聚类级别,进行聚类合并。依据岩芯资料岩相识别及测井相标定结果,本文最终建立了包含5个聚类的测井相划分模型及岩相定量解释图版,其中,聚类测井相1~5分别对应于泻湖泥、石膏坪、滩间、低能滩及高能滩,取芯段符合率达85%以上,能够较好的运用于非取芯段岩相预测研究。据此,我们进行了连续的聚类测井相划分及岩相预测,并对层序地层格架内岩相分布及物性特征进行了分析。

简介：多分量地震数据中低频缺失是弹性波全波形反演中的一大难题，低频的缺失导致全波形反演无法有效恢复介质的长波长成分进而使反演陷入局部极值。为此，本文提出了反演介质纵横波速度长波长分量的弹性波包络反演方法。该方法利用包络算子具有的解调多分量数据中隐含的低频信息的能力，构造多分量地震数据的包络目标函数进行反演，用以恢复地下介质纵横波速度的长波长成分。一系列数值试验表明，即使在多分量地震数据中缺失低频信息、并且初始模型缺少先验信息的情况下，这种弹性波包络反演方法能够有效降低波形反演的非线性，可以为后续的常规弹性波全波形反演或者深度偏移提供足够精确的初始模型，且该方法对横波速度长波长分量的重建尤为有效。Marmousi-2模型的高精度纵横波速度的反演结果表明，利用该方法反演的纵横波速度作为常规弹性波全波形反演的初始模型，可以显著提高反演结果的精度。此外，本文对弹性波包络反演方法的适用性也进行了初步的研究与讨论。

简介：本文利用x射线CT获取反映岩心微观结构的三维数字岩心，利用数学形态学中的开运算模拟了岩石的油驱水排驱过程中，不同含水饱和度下油和地层水在孔隙空间中的分布。利用有限元方法计算了岩石电阻率，进而得到岩石地层因素和电阻率指数，并考查了岩石润湿性对岩石电阻率指数的影响。数值模拟结果表明：基于数字岩心的水湿岩石地层因素和电阻率指数数值模拟结果与实验结果一致，拓展了岩石电阻率实验的能力；岩石润湿性对岩石电性有重要影响，在相同含水饱和度下，油湿岩石电阻率高于水湿岩石电阻率，油湿岩石饱和度指数远大于水湿岩石饱和度指数。

简介：逆时偏移被认为是处理陡倾角断层、回转波和逆掩推覆构造最有效的成像方法。考虑到当今勘探对地下复杂构造体、薄互层和隐蔽储层的高分辨率成像需求，井间、井中地震资料成为地面地震资料的一个很重要的补充。但由于地面地震、井中地震和井间地震存在很大的频率差异，本文将相应的变网格算法应用到偏移中。另外，也将Lanczos滤波引入到波场外推计算中，从而可以很好地消除变网格引起的人为反射。试算表明，对比地面、井中和井间的偏移，可以看到微构造、陡反射层都得到了更好地刻画和描述。基于变网格的逆时偏移在储层预测和岩性解释方面有广阔的应用前景。

简介：在复杂地表条件的区域,地震数据的采集和处理是一项极大的挑战。虽然可以通过静校正来消除起伏地表的影响,然而当地表高程以及近地表速度剧烈变化时,简单的垂直时移对地震波场造成的畸变会严重降低偏移成像的质量。基于射线的偏移方法可以直接在起伏地表面进行波场的延拓成像,是解决上述问题的有效手段。本文针对复杂地表条件下的高斯波束叠前深度偏移进行研究,对倾斜叠加公式进行修改,使之包含地表高程以及速度的信息,通过直接在复杂地表面进行平面波的合成,得到了一种具有更高成像精度的改进方法。首先简单介绍常规高斯波束偏移的基本原理和计算流程,并以此为基础,给出复杂地表条件下高斯波束偏移原有的实现方法以及本文的改进方法,最后通过模型和实际资料的试算验证本文方法的有效性。

简介：复杂地层曲面构建可广泛应用于石油勘探、地质建模、地质构造分析等领域中,也是这些领域数据可视化和可视分析的重要基础。现有方法在多种类型断层处理、曲面光滑性等方面存在不足。为解决地层数据分布不规则性带来的曲面光滑性问题,本文引入几何偏微分方程的曲面造型方法开展地层曲面构建;针对多种类型复杂断层存在情况下的曲面构建问题,提出三维空间与二维平面的投影互换算法,给出了基于几何偏微分方程的多种形式断层的统一处理方法;针对复杂地层曲面构建问题,建立了相应几何偏微分方程,给出演化求解算法。实际地层数据的空间曲面构建验证了本文方法计算效率高、能处理不规则的数据分布,同时能对存在断层尤其是逆掩断层的地层面进行曲面恢复。

简介：初至拾取是起伏地表地震资料处理的关键环节,直接影响近地表建模的精度及静校正效果,但起伏地表资料极低的信噪比使传统的自动初至拾取算法几乎失效,而手动拾取要耗费巨大的人工和时间成本。本文研究了改进超虚干涉法,结合多道多域初至质量监控技术,实现起伏地表资料初至波自动拾取。改进超虚干涉法首次将近地表散射波纳入干涉法提高信噪比的范围,并通过折射波和散射波的线性组合使干涉法适用于起伏地表条件下任意复杂射线路径的初至波类型;提出反向干涉和多域干涉的概念,显著增强了所估计的虚源信号;使用波形反褶积滤波器较好地抑制了干涉导致的＂假事件＂的形成;采用多道多域初至质量监控技术,实现错误＂假事件＂初至的自动归位,提高了拾取初至的稳定性。本文研究的初至波自动拾取理论与技术具有突出的鲁棒性和稳定批量处理大量实际三维地震数据的能力,在中国西部山区三维地震勘探数据处理的应用效果显著,质量优于某常用商业软件。

简介：碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性，使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型，以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上，研究采用岩石网格化方法，将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体，根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量，并根据不同孔隙连通性进行流体替换，计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量，采用Hashin．Shtrikman—Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析，明确不同孔隙对岩石弹性参数的影响特征，模拟分析结果与实际资料认识一致。

简介：频率域全波形反演充分利用全波场的振幅、相位以及频率信息，采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础，对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中，针对截断边界问题，并考虑多频率联合反演中计算区域采用同一套剖分网格的需求，提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子，通过一定方式调节衰减因子使得声波在衰减层中充分衰减，达到压制截断边界影响的目的。根据指数衰减规律，文中推导出了一种新的衰减因子计算公式，并给出了不同频率条件下衰减层厚度计算公式；在反演算法中，采用共轭梯度法求解高斯牛顿反演迭代方程组，避免直接求解雅克比矩阵和Hessian矩阵带来的巨额计算量，并采用相同的反演模型，对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明，本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。

简介：随着煤矿的深部开采，多层积水采空区探测成为煤田水文物探工作的新内容。中心回线瞬变电磁（简称TEM）法因场结构与地层的耦合关系，对低阻层探测有利。但低阻层的屏蔽作用不仅使得探测同样的深度需要更长的观测时间，而且还会减弱下伏地层的异常响应。本文通过直接时域数值模拟和水平分层大地的模型正演，估算了探测目的层所需要的时间长度，根据噪声对观测数据造成的影响，给出了多层积水采空区可分辨的标准。山西大同达子沟煤矿水文勘探中获得的TEM实测曲线，表现了对多层积水采空区的探测能力。在实测曲线的定量反演解释中，利用电测井资料作为初始参数进行定量反演，约束了煤系薄层的等值性。所推断的三层和两层积水采空区，为钻孔所证实。研究结果表明，当观测时间有足够的长度、下伏地层的异常显示大于观测误差，中心回线TEM法探测多层采空区积水是可行性的。

简介：已有的均匀磁化长方体磁场计算方法大都假设观测点位于上半无源空间，对于起伏地形条件而言，这些方法可能存在解析“奇点”。为此，本文基于地磁场基本理论，采用变量替换的积分方法，导出了改进的磁场表达式，详细讨论并有效解决了整个无源空间的所有奇点问题。相比前人的方法，其积分过程更自然、简单，最后的积分结果形式更加统一，并且不需要坐标变换即可求出无源区任意点处的磁场值，从而简化了正演过程。对比模型试验表明，新导出的磁场无解析“奇点”理论表达式是正确的，并能适应地形起伏的情况。

简介：礁滩储层是我国海相碳酸盐岩油气勘探的重要目标，也是我国油气产能的重要接替领域之一。由于礁滩储层储集空间类型多样，非均质性极强，依靠常规测井判别其沉积相和岩性非常困难。成像测井能够清晰地反映礁滩储层的结构组分和沉积构造，为沉积相和岩性的判别提供了可靠的依据。在对大量礁滩储层的成像和岩心对比观察的基础上，我们提出了9种典型的成像解释模式，建立了利用成像解释模式自动判别礁滩储层沉积相和岩性的方法并研制了相应的处理软件。该方法在塔中和川东北地区礁滩储层的实际应用中取得了良好的效果。

简介：在多个区块的测井评价工作发现,孔隙结构直接影响储集层的品质和油气层的电阻率,是测井准确评价流体性质的关键。岩石物理资料表明不同区块内影响储层孔隙结构的因素不同,但效果是一致的,即孔隙结构的复杂程度控制着储层的储集能力和渗透能力。孔隙结构的复杂程度影响了储层中导电流体的分布和含量,从而控制了储层的电阻率。储层出现低阻油气层的内因均为复杂的孔隙结构（骨架导电及工程原因除外）。测井储层评价在分析控制储层孔隙结构复杂程度的地质因素及储层分类的基础上,针对不同类型储层采用不同的模型、参数和标准,可以有效的认识储层品质和识别不同类型储层的流体性质。

简介：洪泽地区由于沉积的特点,储层横向变化快,油藏受构造、岩性、油源多因素控制.在对该区三维AVO属性体解释中,利用多元回归方法求取了横波曲线,分岩性和含油气性统计了纵、横波、泊松比参数分布规律,建立了本区的含油砂岩的正演模型,从而降低了AVO解释的多解性.通过井-震结合对四种AVO属性数据体进行了标定,并确定了各属性体应用范围,进而进行了储层和含油气检测.实践表明,该方法能有效地利用AVO属性数据体进行储层预测及油气检测,具有一定的推广价值.

简介：白云岩成岩收缩晶间孔、洞、缝与构造网状缝相互沟通可以组成良好的油气储层,但是这种复杂裂缝-孔隙型储层分布随机性强、发育程度和差异性大,储层预测难度大.本文以中国西部酒泉盆地青西油田下白垩统下沟组湖相白云岩裂缝-孔隙型储层为例,提出地球物理综合预测碳酸盐岩裂缝的方法.本文描述了综合地球物理方法预测碳酸盐岩裂缝储层的实施和应用效果.