基于计算岩石物理方法的页岩储层弹性参数提取

发展了应用数值计算方法获取页岩储层的速度、各向异性参数的计算岩石物理系列方法.该系列方法包括了大尺度精细地质模型数值建模、计算网格尺度的地球物理建模和地震波数值模拟提取岩石物理弹性参数.本文方法利用储层的统计数据而不是具体岩心的测量数据，可获得储层岩石物理弹性参数的变化规律.相比于基于岩心测试的岩石物理方法，本文方法可精细考虑实际储层的非均匀特征，可得到岩心测试难以求取的与尺寸效应高度相关的弹性参数，也避免了求取弹性参数变化规律时获取不同地质特征岩心的困难.本文发展了计算岩石物理方法，为计算岩石物理面临的大尺度地质建模和计算能力限制问题提供了有效的解决方案.文中以胜利罗家的页岩储层为例，求得了储层TOC含量从3%到21%变化情况下储层的P波、S波速度以及各向异性参数变化规律.     

基于岩石物理模型的页岩孔隙结构反演及横波速度预测

准确预测储层的等效孔隙纵横比对页岩储层岩石物理建模及横波速度预测具有重要意义。为分析页岩储层孔隙纵横比及预测横波速度,提出了基于岩石物理模型的页岩孔隙纵横比反演及横波速度预测方法。本文首先通过岩石物理模型建立岩石的纵、横波速度与孔隙纵横比、孔隙度和矿物组分等参数之间的定量关系,寻找最佳孔隙纵横比;然后通过使理论预测与实际测量的纵波速度之间误差达到最小的方式反演孔隙纵横比,并以此为约束预测横波速度。实际测井数据反演结果表明,龙马溪组页岩地层的孔隙纵横比稳定,而围岩的孔隙纵横比变化范围较大;说明与围岩相比,页岩的孔隙结构更为稳定。同时,预测得到的页岩横波速度与实测横波速度的误差较小,另外对于缺少矿物组分资料的页岩层段,用平均矿物组分预测得到的横波速度误差仍较小;说明与矿物组分相比,龙马溪组页岩的纵、横波速度对孔隙纵横比参数更敏感。综上所述,利用该方法可预测到较为准确的等效孔隙纵横比和横波速度。     

预条件共轭梯度法在地震数据重建方法中的应用

基于最小平方的Fourier地震数据重建方法最终转化为求解一个线性方程组, 其系数矩阵是Toeplitz矩阵,可以用共轭梯度法求解该线性方程组.共轭梯度法的迭代次数受系数矩阵病态程度的影响,地震数据的非规则采样程度越高,所形成的系数矩阵病态程度越高,就越难收敛和得到合理的计算结果.本文研究了基于Toeplitz矩阵的不同预条件的构造方法,以及对共轭梯度法收敛性的影响.通过预条件的使用,加快了共轭梯度法的迭代速度, 改进了共轭梯度算法的收敛性,提高了计算的效率.数值算例和实际地震数据重建试验证明了预条件共轭梯度法对计算效率有很大的提高.     

泥页岩岩石物理建模研究

泥页岩由于其复杂的岩石特性（主要是裂缝及有机质的存在）,目前还没有有效的岩石物理模型可以较为精确的模拟其性质.本文在自洽模型和微分等效介质模型的基础上,引入Berryman三维孔隙形态及Brown-Korringa固体替代技术,建立适用于富有机质泥页岩的新型岩石物理模型.在此基础上进行正演分析,讨论不同孔隙形态对于自洽模型的临界孔隙度以及岩石速度的影响.正演分析的结果表明即使将未知的混合岩石作为背景岩石,微分有效介质模型的引入使得固体相和流体相仍然不是对称的,临界孔隙度不一定要落在0.4到0.6之间.且不同的孔隙形状对于自洽模型的临界孔隙度以及岩石的速度具有明显的影响.此外,基于岩石物理模型,文章讨论了不同孔隙形态、不同泥质含量时有机质对于岩石弹性性质的影响.最后利用一口页岩气井对该模型进行验证,预测的纵横波速度与测井结果吻合的很好,证明了该模型对于富有机质泥页岩的适用性.     

地震数据重建方法综述

地震勘探的目的是为了获得地下构造的精确成像.由于人为因素和环境原因,地震数据在空间方向上往往是不规则采样或缺失采样的,这可能会在数据处理时产生假象,最终导致错误的解释.因此经常需要在空间方向对缺失的地震数据进行重建.重建问题可以看作是一个反演问题,即从不完整的地震数据中重建出完整的地震波场.本文重点研究了国内外比较成熟的地震数据重建方法,如基于滤波的重建方法、基于波场延拓算子的重建方法、基于变换域的重建方法以及相干倾角插值等方法,分析了这些方法的优缺点,以及目前地震数据重建所面临的挑战.