部分叠加资料纵横波速度测井约束反演

本文提出一种角度部分叠加资料同时反演纵波和横波速度的方法,以角度部分叠加资料为基础,利用地震波振幅随入射角变化与弹性参数间的数学关系,基于非线性最优化理论,通过将最小平方问题转化为大型带状矩阵的求解问题,采用逐道外推技术依次求得该角度剖面的每个点纵波和横波速度值,利用地震波速度在各向异性介质中与入射角的关系,进一步得到地震波零入射角纵、横波速度,从而可以求得纵横波速度比及泊松比,为地震资料的岩性及含油气性解释提供了丰富的参数信息,为保证反演的稳定性,降低多解性,采用测井、叠后构造层位解释结果的共同约束,并做好叠前去噪以及振幅保持处理.     

测井约束地震波形反演的同伦摄动法

测井数据和地震数据是地震勘探中两种最重要的资料. 测井约束地震波形反演是在非线性波形反演的基础上,利用已知测井资料详细的垂直分辨能力和地震资料均匀密集的水平采样特点, 通过迭代反演来求取一个具有较高分辨率的速度参数.本文建立了测井约束反演模型,研究了测井约束下地震波形反演的同伦摄动求解方法.同伦摄动法作为一种新的、求解数学物理中各种非线性问题的有效方法,具有计算速度快、计算精度高的优点.这对于提高反演的精度和效率是十分有益的. 为了表征该方法的有效性和稳定性,分别对水平层状介质模型和逆冲断层带模型进行了数值模拟,并与Landweber迭代法相对比,结果表明该算法具有更好的收敛性,能够取得更为满意的反演效果.     

测井约束地震波形反演的同伦摄动法

测井数据和地震数据是地震勘探中两种最重要的资料. 测井约束地震波形反演是在非线性波形反演的基础上,利用已知测井资料详细的垂直分辨能力和地震资料均匀密集的水平采样特点, 通过迭代反演来求取一个具有较高分辨率的速度参数.本文建立了测井约束反演模型,研究了测井约束下地震波形反演的同伦摄动求解方法.同伦摄动法作为一种新的、求解数学物理中各种非线性问题的有效方法,具有计算速度快、计算精度高的优点.这对于提高反演的精度和效率是十分有益的. 为了表征该方法的有效性和稳定性,分别对水平层状介质模型和逆冲断层带模型进行了数值模拟,并与Landweber迭代法相对比,结果表明该算法具有更好的收敛性,能够取得更为满意的反演效果.     

重力与地震资料的模拟退火约束联合反演

联合反演是综合地球物理研究的重要定量解释手段.本文在总结和分析重力与地震资料联合反演的研究现状基础上,利用改进的全局寻优的快速模拟退火算法,实现了重力和地震资料的约束同步联合反演.针对性地设计了密度和速度界面不完全一致的模型,理论模型的试验说明了方法的效果和适用性.结合最近完成的广东徐闻地区实际资料的处理和解释,表明该方法可准确确定复杂构造物性界面的密度和速度结构,在该地区的油气勘探中发挥了作用.在先验信息约束下,该联合反演方法要明显优于单独的重力反演.     

地震数据互相关驱动的多道反演方法

地震反演是储层定量描述和地震油气识别的关键技术,反演结果在复杂构造区域的横向连续性和保真性是影响地震资料定量解释精度的重要因素.基于此,本文发展了地震数据互相关驱动的多道反演方法.考虑地层反射系数与地震数据在结构上具有相似性的特点,基于地震数据互相关描述地层反射系数的结构特征,并将其作为多道地震反演的横向约束条件;此外,为改善地震数据本身横向连续性差对反演结果的影响,在目标泛函的惩罚项中引入局部优化算子,构建了一个易于求解的多道地震反演目标泛函.与常规多道地震反演方法相比,本文方法能够设计更合理、更符合实际情况的横向约束算子,提高反演结果的横向连续性,并且能有效降低地震资料质量对反演结果的影响.模型测试和实际应用验证了本方法的可靠性和稳定性.

     

地震数据互相关驱动的多道反演方法

地震反演是储层定量描述和地震油气识别的关键技术,反演结果在复杂构造区域的横向连续性和保真性是影响地震资料定量解释精度的重要因素.基于此,本文发展了地震数据互相关驱动的多道反演方法.考虑地层反射系数与地震数据在结构上具有相似性的特点,基于地震数据互相关描述地层反射系数的结构特征,并将其作为多道地震反演的横向约束条件;此外,为改善地震数据本身横向连续性差对反演结果的影响,在目标泛函的惩罚项中引入局部优化算子,构建了一个易于求解的多道地震反演目标泛函.与常规多道地震反演方法相比,本文方法能够设计更合理、更符合实际情况的横向约束算子,提高反演结果的横向连续性,并且能有效降低地震资料质量对反演结果的影响.模型测试和实际应用验证了本方法的可靠性和稳定性.     

基于局部相关性约束的三维大地电磁数据和重力数据的联合反演

为解决地球物理反演中多解性的问题,综合多种地球物理信息的联合反演受到了广泛的关注.本文依据不同地球物理响应可能由相同异常体引起,而不同地球物理分布参数之间存在相关性等特点,提出了一种基于局部Pearson相关系数约束的联合反演方法.该方法假设每个局部区域模型参数的分布具有线性相关特性,在拟合不同类型观测数据时,对局部模型参数施加相关性约束,进行联合反演以减少多解性.本文采用交替迭代联合反演流程,改善了同一目标函数下联合反演收敛性和速度问题.基于新的联合反演方法和流程,我们测试了三维大地电磁和重力仿真数据的联合反演.结果表明,本文提出的基于局部相关性约束的联合反演方法,能充分利用大地电磁和重力观测数据信息,有效改善单一地球物理反演收敛性和多解性的问题,反演效果得到明显提升.

     

大地电磁与地震资料仿真退火约束联合反演

在对仿真退火算法研究的基础上,将非常快速仿真退火的连续模型空间改为可变的离散化形式,理论模型试算表明,使用本文方法可提高工作效率5-10倍. 为适应联合反演的需要, 对非常快速仿真退火的误差"单判断"检验准则改为误差"双判断"检验准则,避免Metrop olis接受准则的判断失误. 在此基础上,实现了一维及二维起伏地形条件下大地电磁与地震资料多参量信息的联合反演. 理论模型的试算及实例说明了约束联合反演的有效性.     

层析成像反演速度建模方法研究

本文应用层析成像反演方法优化叠前深度偏移速度模型。阐述叠前深度偏移速度建模流程,详细介绍利用两种不同层析反演方法来优化深度域的层速度模型,对同一实际地震资料速度模型建立的应用效果表明,在资料信噪比较高的地区能提供准确解释层位信息的前提下,模型层析反演优化后的速度模型精度要比网格层析反演优化后的速度模型更高,且稳定性强。     

大地电磁与地震资料仿真退火约束联合反演

在对仿真退火算法研究的基础上,将非常快速仿真退火的连续模型空间改为可变的离散化形式,理论模型试算表明,使用本文方法可提高工作效率5-10倍. 为适应联合反演的需要, 对非常快速仿真退火的误差"单判断"检验准则改为误差"双判断"检验准则,避免Metrop olis接受准则的判断失误. 在此基础上,实现了一维及二维起伏地形条件下大地电磁与地震资料多参量信息的联合反演. 理论模型的试算及实例说明了约束联合反演的有效性.     

Two-step interface and velocity inversion

This paper studies the computation method of two-step inversion of interface and velocity in a region. The 3-D interface is described by a segmented incomplete polynomial; while the reconstruction of 3-D velocity is accomplished by the principle of least squares in functional space. The computation is carried out in two steps. The first step is to inverse the shape of 3-D interface; while the second step is to do 3-D velocity inversion by distributing the remaining residual errors of travel time in accordance with their weights. The data of seismic sounding in the Tangshan-Luanxian seismic region are processed, from which the 3-D structural form in depth of the Tangshan seismic region and the 3-D velocity distribution in the crust below the Tangshan-Luanxian seismic region are obtained. The result shows that the deep 3-D structure in the Tangshan seismic region trends NE on the whole and the structure sandwiched between the NE-trending Fengtai-Yejituo fault and the NE-trending Tangshan fault is an uplifted zone of the Moho. In the 3-D velocity structure of middle-lower crust below that region, there is an obvious belt of low-velocity anomaly to exist along the NE-trending Tangshan fault, the position of which tallies with that of the Tangshan seismicity belt. The larger block of low-velocity anomaly near Shaheyi corresponds to a denser earthquake distribution. In that region, there is an NW-trending belt of high-velocity anomaly, probably a buried fault zone. The lower crust below the epicentral region of the Tangshan M _S=7.8 earthquake is a place where the NE-trending belt of low-velocity anomaly meets the NW-trending belt of high-velocity anomaly. The two sets of structures had played an important role in controlling the preparation and occurrence of the M _S=7.8 Tangshan earthquake. Contribution RCEG97006, Research Center of Exploration Geophysics, China Seismological Bureau, China. This project is supported by the Chinese Joint Seismological Science Foundation.     

一种新的地震子波提取与层速度反演方法

粒子群优化算法是近十年发展起来的一种基于群智能的非线性全局最优化新方法.本文详细介绍了粒子群优化算法的基本原理,并将其应用到子波提取与层速度反演中.通过模拟数值算例,从不同角度研究了粒子群优化算法的可行性及其效率问题.试算结果表明,粒子群优化算法在不同分辨率、不同信噪比、不同相位子波合成的地震记录反演中效果明显.     

利用远震接收函数反演重庆地震台下方的速度结构

选取重庆地震台2010年至2012年记录的60个远震宽频带数字地震记录,采用频率域反褶积法获得台站的接收函数,采用H-Kappa叠加方法反演台站下方的地壳厚度和泊松比,作为台站下方波速反演的约束条件,以减少反演的非唯一性.计算结果显示,重庆地震台下方地壳厚度为42 km,与中国大陆中西部地区Moho面深度在38-45 km保持一致.该研究对增强该区的深部地质构造特征、分析孕震机制等具有积极意义.     

基于平均速度分析的近地表纵波速度反演

面波频散曲线对于横波速度、纵波速度、层厚等近地表地球物理参数的敏感度相差较大, 现阶段通过频散曲线可以获得较为精确的近地表横波速度与厚度信息, 但无法直接对纵波速度进行反演.研究表明, 泊松比对于波长(W)-探测深度(D)关系较为敏感.基于这一发现, 本文根据频散曲线与反演获取的横波速度结构计算平均速度, 获取W-D关系曲线.但合成数据测试证明, 因近地表浅层对应的W-D曲线变化较小, 且浅层纵波速度反演不准确会使误差累积, 直接反演W-D曲线无法获取可靠的纵波速度剖面.本文改进了基于平均速度分析的近地表纵波速度反演方法, 在目标函数中加入了对浅层信息较为敏感的W/D-D信息, 同时对W-D曲线与W/D-D曲线进行联合反演.合成数据测试证明可以获取较为准确地浅层与深层纵波速度.将该方法应用于实际地震数据中, 联合反演得到的纵波速度剖面与微测井数据较为吻合, 证明本文提出的方法可以不借助其他信息, 仅通过面波频散信息, 获取更为准确地近地表纵波信息.

     

河北省南部及邻区地壳P波速度结构三维反演

1980—2012年河北省及邻区测震台网地震记录,使用了河北省南部及邻区（34.0°—38.0°N,112.0°—118.0°E）63个固定地震台站和4 540个地震事件,得到27 709条P波到时数据,采用速度结构与地震位置联合反演的方法,获得研究区内地壳P波三维速度结构,重新确定中小地震震源位置。速度结构揭示：研究区域内地壳的P波速度结构存在明显的横向不均匀性,在10—25 km深度上横向不均匀性更加显著;大地震基本发生在速度异常体或高低速交界区域。地震重新定位结果显示：地震P波走时均方根残差（RMS）从1.68 s降到0.82 s;地震呈明显条带状分布,震源深度与地质构造年代具有一定负相关性。     

基于稀疏反演理论的自动叠加速度反演方法

叠加速度分析技术是常规地震资料处理中的重要环节,也是经典的时间域速度建模方法.叠加速度分析技术主要包括速度谱计算和拾取两个步骤.至今为止,多数研究工作通过提高速度谱的分辨率以及抗噪声能力,获得高质量的速度谱从而有利于拾取.本文的目标是将叠加速度分析技术转为一个全自动化的处理流程.从参数估计的角度出发,将叠加速度估计转化为稀疏反演框架下的模型参数估计问题,并通过稀疏反演算法自动反演叠加速度,进而提高叠加速度建模的效率.为实现这一目标,首先给出了正问题的定义,即层状介质中CMP道集的预测模型,利用叠加速度、垂向双程走时(t_0)以及反射子波以及CMP道集时距关系(如双曲时距关系)可以预测CMP道集.接着,速度分析反问题可以描述为已知观测的CMP道集,估计模型参数(叠加速度及t_0时间等).利用模型参数的稀疏性作为约束条件并用L_0范数作为模型稀疏性的度量准则,叠加速度分析可以转化为L_0范数约束下的稀疏反演问题.本文提出了一种基于预测校正思想的匹配追踪算法求解上述反问题,实现了自动叠加速度建模并为后续的高精度速度反演方法提供较好的初始模型.理论和实际资料的测试结果证明了本文方法的有效性.     

测井约束地震反演综述

阐述了测井约束地震反演研究的必要性、技术内涵及研究思路，回顾了测井约束地震反演技术的发晨历程，比较了三种非线性测井约束地震反演方法的优缺点。小波变换以其独特的变时窗特性，在测井约束地震反演技术发展中潜力较大，值得深入研究。同时对测井约束地震反演技术的研究现状做了小结，指出当前的研究对象已转向了复杂地质体。特别是海相碳酸盐岩复杂地质体，由于海相和陆相油气田存在许多不同之处，所以测井约束地震反演技术将面临巨大的机遇和挑战。最后，指出了测井约束地震反演技术的发展方向。     

随机反演在储层预测中的应用

针对隐蔽油气藏储层预测的需要,开展了地震反演研究,根据目前的实际应用将储层预测中的基于模型的地震反演分为三个实施阶段:即构造反演、声波波阻抗或弹性波阻抗反演以及岩性反演,并对每个阶段的目的、关键技术及其原理进行了详细描述,尤其是详细描述了基于马尔科夫链的蒙特卡罗随机模拟技术.最后给出了一个综合应用测井、地质、地震资料进行反演,从而进行储层预测的实例.     

叠前弹性参数一致性反演方法研究(英文)

常规叠前反演中,纵波速度、横波速度和密度三参数之间,在反演精度上存在明显的差异,"三参数"的一致性反演遂成为重要的研究目标。本文从导致它们精度差异的根源入手,提出了新的叠前反演算法和思路,通过合理的近似,构成参数间的互动和相互约束,使三个参数的反演精度得以同步提高。理论模型试算和实际资料应用表明,三个弹性参数均有较高的反演精度且保持了一致性,与理论模型和实际资料吻合。该方法具有较好的应用前景。     

测井约束的地震反演在鄂尔多斯盆地大牛地气田中的应用

有测井约束的地震反演技术是国内外油气储层预测的关键技术之一,本文介绍了储层预测中常用的反演方法,对地震反演技术在储层预测过程中的一些关键技术环节作了初步探讨,重点介绍了测井约束的地震反演方法在鄂尔多斯盆地储层预测中的应用情况,最后就反演结果作了初步的跟踪分析．关键词地震反演,地质统计反演,储层预测。