一种基于阵列声波测井的胶结指数计算方法在饱和度评价中的应用

利用饱和度公式评价储层含油性时,胶结指数一般在某一层位取某一固定值.在复杂岩性储层中,储层非均质性强,孔隙结构复杂,连通性差,导致胶结指数变化较大,利用常数胶结指数计算的含油饱和度精度降低,给储层测井评价带来了极大的困难.基于弹性介质理论,利用阵列声波数据得到地层骨架的剪切模量和体积模量,并计算出骨架的纵波速度.分析发现,地层和骨架的纵波速度比与地层孔隙度和胶结指数具有一定的函数关系,通过地层和骨架的纵波速度比和孔隙度可以计算出连续的胶结指数,进而得到地层饱和度.实验结果表明,本文方法计算得到的地层饱和度与岩心分析饱和度吻合较高,方法一定程度上克服了传统的饱和度模型在复杂岩性储层中的不适用性,为复杂岩性储层含油性准确评价提供了依据.     

地震波数值模拟的非规则网格PML吸收边界

以格子法为基础,以声波方程为例研究非规则网格PML(Perfectly Matched Layer)方法.本方法的核心是建立局部坐标系下的分裂方程和基于积分近似的微分方程弱形式.该非规则网格模拟方法允许在计算域内设置任意形状的人工边界.对于二维半空间问题,与采用矩形人工边界相比,采用半圆形人工边界可减少计算量20%以上.采用光滑的曲边界,不仅可减少计算区域,还可避免常规的PML吸收边界在吸收带角点区域的特殊处理.本方法事先计算和存储边界单元的局部几何参数,在计算的每一时间步查表调用这些参数,与常规的直边界PML方法相比,不增加任何计算量.     

复杂地质构造波动方程数值模拟

利用全声波方程进行地震波场的数值模拟，虽然可以计算全声波波场值，但是也常常会出现层间多次波，给资料解释造成一定的困难。作者在本文中根据无反射声波方程，采用有限差分数值解法，实现了自激自收时间剖面的制作，避免了层间多次波的产生。本方法不但允许速度有纵向变化，而且允许速度有横向变化，这使得模拟复杂地质构造的地震响应成为可能。实际计算结果表明，该方法是一种行之有效的自激自收时间剖面的计算机数值模拟方法。     

一种优化的频率域三维声波有限差分模拟方法

为提高频率域有限差分（FD,finite-difference）正演模拟技术的计算精度和效率,基于旋转坐标系统的优化差分格式被广泛应用,但是只应用于正方形网格的情况.基于平均导数法（ADM）的优化差分格式,应用于正方形和长方形网格模拟.这些频率域有限差分算子,各自具有不同的差分格式和对应的优化系数求解表达式.本文基于三维声波方程发展了一种新的优化方法,只要给定FD模板形式,可直接构造频散方程,求取FD模板上各节点的优化系数.此方法的优点在于频率域FD算子的优化系数对应各个节点,可扩展优化其他格式.运用此优化方法,计算得到了不同空间采样间距比情况下27点和7点格式的优化系数.数值实验表明,优化27点格式与ADM 27点格式具有相同的精度,优化7点格式比经典的7点格式具有更小的数值频散.

     

一个近似的VTI介质声波方程

声学近似是提高各向异性介质准纵波数值模拟和偏移处理计算速度的有效方法之一.为了获得VTI(具有垂直对称轴的横向各向同性)介质声波方程,根据VTI介质的精确相速度表达式,用待定系数法简化其中的根号项,得到一个近似相速度公式,由此公式通过反Fourier变换推导出一个时间二阶、空间四阶偏微分方程.定量相速度计算表明,当Thomsen参数ε=δ(即相应介质为椭圆各向异性介质)时,由该方程所确定的相速度是精确的;当ε≠δ时,该方程所确定的相速度随ε和δ差值的增大逐渐增加.该方程所确定的最大相速度百分比相对误差当ε=0.1、δ=0.2时为0.13%,当ε=0.8、δ=0.3时为-1.65%.有限差分数值模拟算例表明该方程是一个纯准纵波方程,其数值模拟波场快照中没有准横波.     

基于多重互相关函数分析剪切波速

介绍了一种提高单孔法剪切波速测量精度的新方法-多重互相关函数法。考虑到判断场地剪切波到时差的困难,提出了用多重互相关函数来求得两剪切波的到时差。互相关分析可求取两相似信号的到时差,但当两信号相似性比较差时,所求得到时差的误差较大;自相关函数本身可进行多重自相关分析,互相关函数结合多重自相关函数后,也可以进行多重互相关分析来求取两相似信号的到时差。试验表明：相比互相关函数,多重互相关函数可改善触发信号和井中三分向探头输出信号的相似性及更好地抑制噪声,因此能更精确地计算出触发信号和井中三分向探头输出信号的的到时差。基于此原理设计了剪切波波速的测试系统,此系统能自动分析出测点的剪切波速和形成完整的报告。     

北山花岗岩热破裂室内模拟试验研究

岩石热破裂是高放废物地质处置工程中需深入研究的课题。对我国高放废物重点预选场址甘肃北山的花岗岩开展室内热破裂模拟试验研究,采用多通道温度测试仪、声发射、波速层析成像和数码显微镜等手段研究了该花岗岩热破裂过程。试验表明,（1）热破裂从试件端部开始产生,逐步向内缓慢扩展,表现出分段性和独立性;（2）根据声发射撞击率可将热破裂可分为稳定热损伤、宏观裂纹形成、宏观裂纹扩展、裂纹冷却闭合4个阶段,声发射定位的时空演化规律清楚地揭示了裂纹从试件上端部向内部扩展的规律;（3）波速层析成像指示了宏观裂纹位置及高温对岩石造成显著损伤的区域,热应力产生的损伤集中在试件边界,范围小,损伤严重,高温造成的损伤集中在钻孔附近高温区,范围较大,损伤略轻微;（4）监测多通道温度,获得了试件内的温度场并为数值模拟参数选取提供验证,采用有限元程序进行了热力耦合数值模拟,从机制上初步解释了热破裂现象,研究认为综合声发射实时监测热破裂过程和波速层析成像能实现对热损伤的量化的特性可实现岩石热破裂的动态监测和损伤量化,为今后地下实验室相关试验的开展和认识高放废物处置长期稳定性做了有意义的探索。     

声学法深海热液速度场测量重建算法研究

介绍了利用声学法测量深海热液速度场的基本原理.利用往返飞渡时间差与流场速度的关系,应用反问题求解技术重建测量区域的速度场.应用最小二乘法对两种典型的热液口速度场模型进行了仿真重建,分析了换能器数量、实验测量误差以及换能器的布放对重建结果的影响,并对重建结果进行了流量分析.仿真重建的绝对误差、相对误差以及流量分析结果表明...     

花山岩画开裂岩体声波测试与分析

阐述了声波测试的方法和过程,通过在岩体中声波的传播速度,建立开裂岩体厚度与劣化程度的数学模型。结果表明,开裂岩体厚度小于30 mm时,强度损失较大;开裂岩体厚度大于30 mm时,强度基本保持不变,为今后花山岩画的保护提供科学可靠的理论依据。