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

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

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

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

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

简介：裂缝识别对碳酸盐岩储层测井评价至关重要，传统的测井仪器由于探测深度浅（小于3米）而无法对井壁外围（大于3米）裂缝发育情况进行评价，远探测声波测井仪器采用相控阵发射、同时加大源距，有效提升了测井仪器的探测深度。但由于缺少正演模拟研究，目前对于裂缝的解释往往是基于经验而缺乏理论依据，导致很多现象难以解释。本文利用高阶有限差分方法对远探测声反射波测井裂缝识别进行了正演模拟及叠前逆时偏移成像研究，首先在理论研究的基础上构建正演模型，重点研究不同裂缝离井壁距离、裂缝张开度和倾角的响应特征；其次在单因素变化基础上提取成像区域的能量强度，分析确定出在实际地层速度有波动变化时远探测声反射波测井方法能够识别裂缝的条件；最后通过对裂缝识别的影响因素定量化分析，确定裂缝识别的最大距离、最小张开度和最小倾角，降低了裂缝识别与评价中的多解性。研究成果对远探测声反射波测井仪器的发展、数据处理方法的改进，以及后续的测井解释工作都有一定的借鉴意义。

简介：基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。

简介：有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程，高阶显式有限差分方法通常用来提高求解精度，已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数，发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明：忽略这些小的差分系数能够近似维持有限差分的精度，但是显著减小了计算量。然后，引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性，采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后，给出了针对均匀和非均匀介质模型的弹性波模拟例子，表明了本文方法的优点。

简介：地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层（PML）边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件（HABC）具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。

简介：本文以中观孔隙结构的White模型为基础，构建了部分饱和孔隙介质模型，利用Biot方程的建立思路和Johnson推导的体变模量，推导了部分饱和孔隙介质中的纵波方程，并以平面波为例，求取了方程的衰减系数，分析了地震勘探频带范围内地震波的衰减特性。结果表明：在部分饱和孔隙介质中，地震波在低频段也会发生明显的衰减和频散现象，频率越大，衰减越大；且第二纵波的衰减比第一纵波更为明显；这一结论弥补了Biot理论在描述地震勘探频带范围内波的衰减现象的不足。文中还研究了孔隙度、饱和度和模型内径尺寸对第纵波衰减特性的影响机理，主要表现在在地震勘探频带范围内，波的衰减随孔隙度的增大而增大，随含油气饱和度的增大而减小，当孔隙内径尺寸小于二分之一外径尺寸时，波的衰减随内径尺寸的增大而增大，当内径尺寸大于二分之一外径尺寸时，波的衰减随内径尺寸增大而减小。

简介：分数阶S变换（FRST）具有较强的时频聚集性。利用FRST处理地震数据，通过合适的分数阶参数将频率轴旋转到适当位置，即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同，因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体，并利用共频率数据体进行了低频伴影分析，然后提出FRST和盲分离结合的方法，不需要对地震数据的最优分数阶参数进行估计，即可提取识别有效地质特征信息的独立频谱，提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据，并与已知井信息进行比对，验证了其有效性。

简介：分数阶S变换(FRST)具有较强的时频聚集性。利用FRST处理地震数据,通过合适的分数阶参数将频率轴旋转到适当位置,即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同,因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体,并利用共频率数据体进行了低频伴影分析,然后提出FRST和盲分离结合的方法,不需要对地震数据的最优分数阶参数进行估计,即可提取识别有效地质特征信息的独立频谱,提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据,并与已知井信息进行比对,验证了其有效性。

简介：目前对动电测井的研究主要集中于数值模拟和岩石物理实验。为研发出能下井测量的动电测井仪器,进行了一系列实验研究工作。文中介绍了一种融合了不同结构和材料要求的新型声电复合式动电测井探测器的结构及原理,对声、电信号接收器在仪器机械设计上无法布置在同一位置的难题也提出了解决方法。在水池中对其辐射声场特性进行了测试,计算了辐射声场声压及发射电压响应,分析了基于相控线阵的阵列式声波发射探头的指向性及实际应用效果。研究表明,在源距1500mm处产生的声压可达47.2kPa,且会随着激励声源主频增大而减小。随着相邻声波发射子阵激励信号延迟时间的增加,辐射声束主瓣会发生明显偏转,主瓣方向的能量也逐渐增强,可有效增强声电转换效率。本文为探测器研制及开展井下动电测井研究奠定了重要基础。

简介：特征值相干属性是一种高效、成熟的裂缝预测技术。但即使在连续、稳定的目的层段,部分裂缝隐没于相近高灰度值区域而难以识别;在破碎带发育区域,由于成像能量难以完美聚焦,中、小尺度的裂缝区域呈现相近低灰度值的云雾状模糊而无法区分。本文提出一种以直方图均衡化处理为核心技术的新型裂缝增强方法。该方法能够强化相干图像中不连续性信息与背景的差异性,突出裂缝的线性结构:采用相干图像阈值逐级调节的方式,实现对不同尺度裂缝的预测。该方法同时提升了显性、隐性裂缝识别能力和精度。更多还原

简介：在具有垂直对称轴横向各向同性介质中,利用四种参数来确定中间至远偏移距转换波(C-波)动校正.它们是C-波叠加速度VC2,垂直速度比和有效速度比γ0和γeff,以及各向异性参数χeff.我们将这四种参数作为C波叠加速度模型.C-波速度分析的目的就是确定这种叠加速度模型.C-波叠加速度模型VC2,γ0,γeff,和χeff可以由P-波和C-波反射动校正资料获得.然而错误的传播是C-波反射动校正反演中的严重问题.当前短排列叠加速度由于是从双曲线动校正推算而得,因而其精度不足以为各向异性参数提供有意义的反演值.中间偏移非双曲线动校正不再被人们所勿略,而是可以用一个背景γ加以量化.非双曲线分析通过中间偏移距的γ校正量可以产生VC2,若数据不含燥音,其误差小于1%.方法稳健,允许γ启始假定值的误差达20%.该方法也适用垂直非均匀各向异性介质.精度的提高使能够用4分量地震资料计算各向异性参数.为此提出了两种工作流程:双扫描和单扫描流程.理论数据和实际数据的应用表明这两种流程得出的结果其精度相似,但是单扫描流程比双扫描更有效.

简介：Tikhonov正则化(TR)方法在重磁数据处理中发挥了重要的作用,本文在研究如何利用Tikhonov正则化方法方法解决重力数据3D反演的同时,深入讨论了可进一步提高拟合误差的ExtrapolationTikhonov正则化方法(EXTR)的原理,并就其参数选择方法及各参数对拟合误差、迭代次数及反演结果的影响进行研究。常密度及变密度组合模型试算结果表明,与TR方法相比,EXTR方法不仅可以达到解释人员设定的先验拟合误差水平,在计算时间及迭代次数相应增加的前提下有更高的拟合精度;同时其反演结果也更加紧致,进一步改善了TR反演结果的发散性;并且其反演数据范围更贴近预设模型参数范围,模型特征与预设模型密度分布吻合较好。更多还原

简介：相干属性是检测地下介质不连续性特征的主要参数之一，但倾斜地层会对相干属性产生干扰，故在计算相干属性时需要消除地层倾角的影响。传统的相干属性计算一般使用时域倾角扫描法来消除地层倾角的影响，但这种方法的精度受到采样间隔等因素的限制，并不能很好地消除地层倾角的影响。根据振幅谱具有不受地层倾角影响的特点，本文提出将预先定义的子数据体内的地震数据变换为振幅谱，再利用振幅谱构造协方差矩阵，协方差矩阵最大特征根与矩阵迹的比值作为子数据体中心点的相干值。利用振幅谱计算相干属性的方法能够较好地消除地层倾角的影响，而且可以使用不同频段的振幅谱计算相干体，具备多尺度相干体算法的特性。模型数据和实际地震数据的应用结果表明，振幅谱相干属性可以更好地消除地层倾角对相干属性的影响。

简介：地震数据重构是地震数据处理的重要步骤之一,重构算法的精度、效率与抗噪性是地震数据重构技术的核心研究内容。研究针对傅里叶域凸集投影(POCS)算法,在定义的最优阈值评价标准基础上,提出了反比例阈值模型,该模型具有在大系数区间比指数模型更快下降速率、而在小系数区间比指数模型更慢下降速率,从而在保证弱反射信号重构精度的同时有效提高POCS地震数据重构算法计算效率。为提高反比例阈值对不同地震数据特点的适应性,在地震数据谱能量分布差异性特征分析基础上,研究提出了在反比例阈值模型分母上增加适应地震数据谱能量特征的因变参数,通过调节该因变参数获得适应不同地震数据特点的最佳阈值曲线,进一步提高算法的计算精度与计算效率。为了实现重构过程中随机噪音的自适应衰减,提高重构后地震数据信噪比,研究提出了数据驱动的加权回加系数计算策略,利用每次迭代对应数据驱动阈值占阈值区间的百分比获得加权回加系数。研究将新方法应用于模拟三维数据和实际三维地震数据,分析结果表明反比例阈值相对传统阈值在提高数据重构计算效率和精度方面具有明显的优越性,新提出的加权回加系数计算策略能有效提高重构数据的信噪比。

简介：在水库现场试验了RISK2型探地雷达探测水库冰厚度的能力,试验时所用天线频率为600MHz;同步钻孔测量雷达探测处的冰厚度;以及在一个点上取样测试分析冰晶体、冰内气泡和冰密度。试验时冰面积雪厚度0.03-0.05m,冰层上部有0.24m粒状冰,其下均为柱状冰;冰内气泡含量呈表层高底层低分布;冰密度随气泡含量变化;冰厚度在平面内不均一。通过探测厚度和实测厚度的对比分析以及气泡含量对介电系数影响的理论分析,建立了积雪、粒状冰和柱状冰三层介质模型,获取雷达波在冰内的理论传递时间。结果发现：能够利用等效介电常数或等效传播速度评价雷达波传递时间,结冰期冰层1/3深度处的对应介电常数或传递速度可以作为等效值;另外因冰内大气泡造成的理论传递时间大于雷达探测时间,其差值随理论传递时间或冰厚的增加呈非线性增加。

简介：目前在研究地球物理参数的空间变异特性时常常是基于几何各向异性假设,然而实际中几何各向异性假设常常不成立。本文提出了一种求解地球物理参数各向异性的新方法,首先对各方向上的变异值作归一化处理,然后用拟合的方法得到各方向单位步长上的变异性强度,进而实现对全局各向异性的求解。以上海地区合成孔径雷达干涉测量（InSAR）中差分大气延迟样本数据为例,验证了新方法求解各向异性的优越性：新方法的偏差只有6.4%,而基于几何各向异性假设方法的偏差达到了21.2%,而且新方法的运算效率有非常显著的提高。进一步,通过克里金插值交叉验证实验,证明了基于新算法得到的结构函数的插值效果最好,从而从另外一个角度验证了新算法求解的各向异性更准确,更好地表征了区域化变量的空间结构特征。因此,与几何各向异性假设的方法相比,新方法不仅能更准确的描述地球物理参数的各向异性,而且有更高的运算效率,为更准确地估计出所需的地球物理参数奠定了基础。