简介：Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。

简介：目前观测系统设计中接收道数的选择基本上是基于水平层状介质假设，这显然不能满足复杂地区采集的需求。本文在朱金平研究基础上，通过野外试验法获得最优接收道数问题，进行数值模拟实现，在固定道间距情况下，建立接收总能量与平均效率能量之间的数学模型，获得最优接收道数，并利用一个复杂工区的地震数据对方法的正确性进行验证，模型数据计算的结果与实际数据处理的结果表明，本文方法与实际处理结论一致，说明该方法具有较高的科学性，可以为复杂地区观测系统设计提供新的思路。

简介：常规地震观测系统设计方法基于地下水平层状介质的假设,通常不能适应复杂构造情况。我们从控制照明的思想出发,提出了一种面向目标成像的地震观测系统设计方法,该方法需要一个由初步地震解释得到的速度模型。利用单程傅立叶有限差分波场传播算子将目标层的平面源延拓到地表,通过分析从目标层延拓到地表的波场能量的分布,可以确定目标层成像所需要的炮点或者检波点的位置。利用SEG-EAGE盐丘模型数值试算结果表明,该方法用于设计面向目标成像的特定地震采集系统。

简介：本文提出了一种将高分辨率阵列侧向和方位电极系综合在一起的三维侧向测井电极系3D-LS,该电极系具有径向、纵向和周向探测能力。通过有限元数值模拟计算,考察了井眼尺寸、冲洗带电阻率、侵入深度、层厚及围岩电阻率对六种不同探测模式的影响,确定了电极系尺寸和探测特性。分析伪几何因子,低侵时电极系的探测深度最深可达1.5m,其值接近斯伦贝谢双侧向电极系深探测深度,而大于高分辨率方位侧向成像仪深探测深度,并且三维侧向测井电极系可提供多条径向不同深度曲线,可更好地描述地层侵入剖面。无限厚地层条件下,方位电极可识别出厚度0.1m的异常体,利用方位侧向曲线半幅点对应异常体厚度判断,对异常体纵向分层能力可达0.5m。高阻背景下,异常体的电阻率越低,越靠近井眼,方位越大于15度,越易被方位电极探测。数值模拟结果为后续三维侧向测井电极系的研究奠定了基础,对低阻异常评价具有一定的指导意义。

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

简介：这篇文章重点研究改进的Gabor小波（improvedGaborwavelet，IGW）变）并讨论了它在地震信号处理和解释中的应用ThispaperintroducesanimprovedGaborwaveletanditscompletetransform，andmainlyanalysestheirpropertiesanddiscussesapplicationsofthesepropertiesinseismicsignalprocessandinterpretation。改进的Gabor小波变换具有以下特性：1）IGWT把时域信号映射到时间一频率域，而传统Gabor小波变换把时间信号映射到时间一尺度域；2）IGWT可用于信号分频，通过固定变换的主频参数dominantfrequency，并变换能提取相应的子带信号，且其主频部分的信息与原信号相应频率部分的信息一致，通过调节变换的分辨率因子，变换能有效控制子带信号的带宽；3）用IGWT和IGWIT构建的滤波器有良好的时一频局部性，在指定时一频范围内能实现针对性滤波。文章用仿真实验和实际用例验证IGWT的这些特性，并在提高地震信号分辨率、地震信号分频和识别小断层等地震信号处理和解释等方面的应用中取得良好效果。

简介：地震数据的精确解释和分析非常依赖于所使用算法的稳定性。我们着重于盐丘地震勘探的稳定检测。我们讨论一个基于地震成像中最佳结构特征属性排序分类的盐丘探测新模型。该算法克服了现有的基于结构属性技术的局限性，因为该技术非常依赖于盐丘固有的地质属性与盐丘检测所用的属性数量。该算法综合了灰度共生矩阵（GLCM）属性，Gabor滤波器，以及含有使用属性特征排序信息协方差矩阵的本征结构等属性。将排序前列的属性组合起来形成一组最优的特征集，以保证算法即使在沿盐丘边界没有强反射层的情况下也能有效。与现有的盐丘检测技术相比，本文的算法稳定和计算高效，并能处理小尺度特征集。我用荷兰F3地块评价该算法的性能。实验结果表明，本文提出的基于信息理论的工作流程用于检测盐丘，其精度优于现有盐丘检测技术。

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

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

简介：在井下微震监测中，各向异性模型用来获取事件位置和给出对介质实际描是十分有用的。在各向异性介质中，因为横波的分裂和奇异点的出现，使得波场的结构变得复杂。我们的研究展示了对强烈各向异性VTI介质利用速度模型校准和双偶极震源进行常规处理和运动学反演的结果。因小走时的拟合差，常用的质量评价标准，不是总正确的，人工加入各向同性介质层以降低拟合差，也可能产生非物理模型，负泊松比，事件位置的位移，

简介：本文描述了在全波实验中所见到的三种现象,即:1.纵横波初至的振幅(绝对值)具有相同的变化规律;2.纵横波初至的相位彼此总是反向的;3.纵横波初至的振幅变化周期为2π。经过充分的分析和解释指出这些现象应该出现在全波测井条件下,从而为利用相位差从全波中提取有用的信息奠定了基础。更多还原

简介：静冰压力是影响水工结构物安全运行的主要因素之一。现有静冰压力检测方法主要存在两个局限，一是实时性差，二是传统电阻应变式压力膜盒传感器因边壁效应影响测量结果的正确性。针对以上两个不足之处，作者基于反射式强度调制原理，研发了一种带杆伞状圆盘结构的压力膜盒静冰压力光纤传感器，实现了利用光纤传感技术检测静冰压力。我们用研制的新型光纤传感器在实验室完成了-30℃～5℃温度范围内冰生长过程中及5℃～-30℃温度范围内冰消融过程中静冰压力的测量试验。实验结果表明，与传统电阻应变式压力膜盒传感器相比，该传感器工作稳定、分辨率达0．02kPa、灵敏度达2．74×10-4／l（Pa，并可实现自动连续测量。

简介：静冰压力是影响水工结构物安全运行的主要因素之一。现有静冰压力检测方法主要存在两个局限,一是实时性差,二是传统电阻应变式压力膜盒传感器因边壁效应影响测量结果的正确性。针对以上两个不足之处,作者基于反射式强度调制原理,研发了一种带杆伞状圆盘结构的压力膜盒静冰压力光纤传感器,实现了利用光纤传感技术检测静冰压力。我们用研制的新型光纤传感器在实验室完成了-30℃~5℃温度范围内冰生长过程中及5℃~-30℃温度范围内冰消融过程中静冰压力的测量试验。实验结果表明,与传统电阻应变式压力膜盒传感器相比,该传感器工作稳定、分辨率达0.02kPa、灵敏度达2.74×10-4/kPa,并可实现自动连续测量。更多还原

简介：研究了基于Poynting矢量的角度域逆时偏移成像及成像幅度的校正方法。根据Poynting矢量进行波场角度分解，由此构建局部成像矩阵及局部照明矩阵。在局部成像矩阵中建立的角度域成像条件，有效地消除了低波数干扰，同时可在局部成像矩阵中进行角度域共成像点道集抽取、倾角估计等运算。利用局部照明矩阵进行了基于全波波动方程的时域照明分析，在局部照明矩阵中计算倾角域幅度校正因子。根据逆时偏移的像计算共倾角像，利用校正因子对各角度的像进行校正，进而实现对成像结果的校正，从而实现了一种高效的倾角域幅度校正方法。最后通过SEG／EAGE模型进行数值计算验证了文中所述的计算方法。

简介：现阶段的叠前地震反演技术中用于描述反射系数与纵、横波速度和密度之间的关系几乎完全是Zoeppritz方程的近似式，由于这些近似公式在大角度和弹性参数变化剧烈时误差较大，这不仅降低了反演解的精度，而且增加了叠前反演的多解性。本文探索了直接利用Zoeppritz方程求解精确反射系数的理论方法，并基于广义线性反演理论详细推导了基于叠前大角度地震资料的纵、横波速度和密度三参数同步反演算法，同时在反演过程引入正则化约束阻尼因子和共轭梯度算法，有效降低了反演的不适定性和提高了反演收敛性。理论模型试算和实际工区应用表明，本文提出的反演方法能够有效利用大角度（一般入射角〉30°）的叠前地震数据，获得更精确的地震弹性参数反演结果，并且反演结果忠实于地震资料，与井吻合较好。

简介：本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。

简介：由于CRS叠加考虑了反射层的局部特征和第一菲涅耳带内的全部反射，从而更充分地利用了多次覆盖反射数据的信息。就目前的地震资料处理技术而言，它是最佳的零偏移距成像方式。本论文利用改进型的参数优化技术，得到高质量的CRS运动学参数剖面，并利用参数剖面计算出叠加孔径，实现了基于最优孔径的CRS叠加，使CRS参数的用途得到了充分利用。模型数据和实际资料的试算表明，基于最优孔径的CRS叠加的成像剖面与传统CRS叠加剖面相比，有着较高的信噪比和同相轴的连续性。

简介：阵列声波信号是典型的非线性、非平稳信号，Hilbert～Huang变换（HHT）是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解（EMD）和对瞬时频率的求解，可以获得声波信号的时一频谱。其关键技术就是进行经验模态分解，任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征，在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程，并列举了一些声波测井波列实例，说明了该方法的有效性。

简介：提高地震资料分辨率的一个有效途径就是衰减补偿,通过对地震波的衰减和频散效应进行校正,提高地震资料的分辨率。常规衰减补偿方法都是基于波场延拓的反Q滤波方法。本文利用Futterman衰减模型,导出了一种衰减介质中合成地震记录的计算方法,在此基础上将衰减补偿问题归结为一个Fredholm积分方程反问题,利用反演方法来实现衰减补偿。针对衰减补偿问题的不稳定性,利用Tikhonov正则化方法提高反演过程的稳定性,数值模拟资料和实际资料处理结果验证了方法的有效性。

简介：推导了二维TTI介质的相速度表达式,并且依据推导出来的相速度表达式,模拟并分析了二维TTI介质相速度的传播快照以及TI介质相速度的传播快照;对比并分析了TTI介质和TI介质模型的相速度理论计算值的X分量特征的差异。TTI介质的相速度研究具有较高的理论研究价值和实际应用价值．