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

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

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

简介：本文通过人工变换T2分布和建立管-球模型模拟法研究含水合物地层渗透率与水合物含量之间的关系。首先,在渗透率的模拟试验中,我们改变了束缚水与可动水的比例、总孔隙度以及与之关联的T2分布。试验结果表明,相对渗透率与水合物含量之间的关系受到这些因素的制约。随后,我们用管-球模型表示水合物生长的孔隙空间,并把水合物的生长过程看成是向孔隙空间随机扔小球的过程。在此过程中,采用两种方法计算渗透率,一是Schlumberger＇sT2公式（即SDR模型）,二是Darcy定律与Poiseuille流动方程相结合的方法。前人的实验研究表明,在一定的水合物含量范围内,渗透率基本保持不变。以此为参考,我们将计算结果与之进行比较。我们发现,采用SDR模型时,渗透率的数值模拟曲线与Masuda模型N=15时的结果相近。而采用Darcy定律时,渗透率模拟值较高,但与实验结果的趋势相一致,都会出现渗透率的平直阶段。尤其,当水合物晶体在孔隙体内优先生成时,优先的概率越高,渗透率的平直范围越大。

简介：渗透率评价是储层评价中的一个复杂问题，传统测井方法难以对储层渗透率参数进行直接和有效的评估。储层渗透率对钻井过程中的泥浆滤液侵入有着较大的影响，因此本文提出一种利用泥浆侵入效应来评价储层渗透率的数学模型和方法。首先构造含泥饼的泥浆侵入数值模型，然后以达西渗流理论为基础导出泥浆侵人深度和储层渗透率的近似数学关系，以此可以利用侵入深度测量值来估算储层渗透率。对孔隙度、渗透率和含水饱和度三个主要储层参数的敏感性分析，发现该方法适用于受到泥浆侵入影响的低孔隙度、低渗透率的油层或油水同层。采用现场测井和取芯数据进行数值模拟计算，结果表明估算出来的渗透率曲线和预设的渗透率曲线吻合较好，证明该方法具有一定的可行性。

简介：基于斑块饱和模型计算渗透率变化的地震反射特征，为流体流动性的地震描述提供依据。基于传播矩阵理论设计反射系数与合成记录算法，实现了频率域岩石物理模型与地震响应计算的无缝连接。斑块饱和储层地震响应包含如下动力学信息：分界面处波阻抗差异、储层内部波的频散与衰减，以及顶底界面波的调谐与干涉。模拟结果表明，渗透率的增加显著降低纵波速度，使其在高、低频弹性极限之间发生频散。储层速度频散与层状构造共同导致反射系数的频变现象。在储层与围岩波阻抗接近的情况下，地震响应对渗透率变化具有敏感性，对于不同储层厚度，当围岩为高速页岩时，反射波叠加振幅随渗透率增加而增加；当围岩为低速页岩时，叠加振幅随渗透率增加而降低。

简介：本文利用取自南海神狐海域钻井岩芯进行了核磁共振测量,研究了岩石中所含水合物对地层渗透率的影响,该岩芯样品为一＂孔隙充填型＂水合物样品。核磁共振测量结果表明岩石中所含水合物对地层渗透率影响很大,随着含水合物饱和度的增加,地层渗透率迅速降低。通过对测量结果的进一步分析可知＂孔隙充填型＂水合物对地层渗透率的影响-存在多方面的因素,包括水合物对地层孔隙度、孔隙迂曲度的影响以及对水合物对孔隙喉道的堵塞等。最后本文利用Masuda下降模型建立了适用于南海神狐海域的相对渗透率与含水合物饱和度的经验公式,下降指数N=7.9718。

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

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

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

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

简介：阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期（超过100秒）,逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。

简介：随着全张量重力梯度（FTG）测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU（GraphicsProcessingUnit）和预处理分解技术（Preconditioningmethods）在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛（SSOR）技术与不完全乔列斯基分解共轭梯度算法（ICCG）提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATeslaC2050GPU的并行SSOR-ICCG算法和在2.0GHzCPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。