裂缝介质中不同震源横波分裂现象的数值模拟

采用TTI介质中二维三分量一阶交错网格应力-速度弹性波方程,模拟胀缩源、垂直集中力源、剪切源在垂直裂缝介质中的井间地震波场,分析了不同方位角情况下波场的传播特征,并讨论了横波分裂时的能量分配情况。结果表明：地震波从各向同性介质进入各向异性介质时会观测到快慢横波;在一定方位角的情况下,垂直裂缝介质中会得到快慢横波,这为进一步分析横波分裂特征和地震实际数据提供了依据。     

对各向异性介质中横波分裂现象的研究

在本文中，我们系统研究了各向异性介质中的横波分裂现象，探讨了各向异性介质中横波分裂现象产生的机理，并给出了引起各种异性的薄层组合以及裂隙定向分布的数学模型。     

高阶旋转交错网格有限差分方法模拟TTI介质中横波分裂

笔者给出了一种能够模拟弹性波在任意各向异性介质中传播的二维三分量高阶有限差分算法.相对于常规交错网格有限差分方法,旋转交错网格有限差分方法在介质具有强差异性时能更精确地模拟地震波的传播,避免常规交错网格中因对弹性系数进行插值而带来的误差.采用高阶旋转交错网格有限差分方法模拟并分析了零偏移距横波分裂现象随裂缝介质方位角和倾角变化的响应特征.结果表明:结合完全匹配层(PML)吸收边界条件的高阶旋转交错网格有限差分方法能获得高精度的地震波场模拟数据,并且在边界具有良好的吸收效果;横波分裂现象主要受裂缝走向与波的极化方向之间的夹角影响,受裂缝倾角影响较小,且快慢横波的能量也跟裂缝走向与波极化方向间的夹角有关.具有倾斜对称轴的横向各向同性(TTI)介质倾角的变化可能会导致记录中波到达时的变化,影响快慢横波的时差.利用横波分裂的能量分布和方位各向异性特征,可以帮助检测裂缝的方位角和倾角.横波在多层TTI介质中传播时会发生多次分裂的现象.     

一种孔隙介质声波模型正演方法

Biot。的孑L隙介质理论描述了地震波动力学特征与油藏参数(孔隙度、 渗透率、粘度等)之间的关系。声波模型方法可用来研究碳氢油气藏特征。 本文基于Biot，声波方程提出一种二维有限差分方法，可用来模拟孔隙介 质低频声波传播。文中讨论了边界条件，网格频散等问题，模拟了一个孔 隙介质层状模型，并给出了其时间剖面和快照     

含直立裂隙介质的弹性波动方程正演模拟

从具有水平对称轴的横向各向同性（HTI）介质中的弹性波动方程出发,在交错网格空间中采用高阶差分算子对弹性波动方程进行差分离散,得到了HTI介质中地震波正演的高阶有限差分格式,研究并实现了PML吸收边界条件。在此基础上实现了HTI介质中弹性波方程的多波正演。数值算例表明,该方法能够精确模拟弹性波在复杂各向异性介质中的传播过程,得到高精度的正演记录。     

VTI介质中弹性波方程正演的一阶混合吸收边界

波动方程正演是研究地震波在地下介质中的传播机理以及提高地震资料采集、处理、解释精度和效果的重要工具,截断边界的处理是波动方程正演的重要研究内容。针对在地震数值模拟中常用的三种(吸收)边界处理方法存在的不足与问题,提出了适用于VTI介质中弹性波动方程正演的一阶混合吸收边界条件,并给出了该方法在过渡区域的线性叠加系数公式。数值模拟实验表明:相同条件下,VTI介质一阶混合吸收边界方法不仅在吸收效果方面优于常规的PML边界条件,而且在存储量和计算效率方面也具有较大优势。     

随机岩溶介质的电阻率正演模拟

受岩溶发育的影响,岩溶地区电阻率分布具有很大的不确定性。采用传统的正演方法在一定程度上不能满足岩溶探测的需要。为了更真实地反映岩溶介质的电性特点,本文提出了一种根据随机过程的谱分解理论来构造随机介质模型的方法,并以此建立了岩溶地区随机地电模型。对理论和实测地电模型的正演计算表明,随机介质能更有效模拟实际地电结构,正演结果更接近实测的异常形态。因此,岩溶地区的电阻率正演模拟计算采用随机介质模型是十分有必要的,一方面有助于对异常的认识,另一方面,随机介质模型可以更好地模拟岩溶地区的地电模型。      

双相介质瑞雷面波有限差分正演模拟

为了研究双相介质瑞雷面波的形成机制及传播规律,促进瑞雷面波资料处理方法的发展。文章根据弹性波动方程,采用交错网格有限差分算法,对二维各向同性弹性介质做解析解与数值解的对比,在此基础上,将PML吸收边界条件,改进的镜像法应用于双相介质波动方程中,并作了稳定性分析,对双相介质水平层状、起伏分界面等典型模型瑞雷面波及体波在内的全波场进行研究。结果表明:基于弹性介质解析解与数值解的对比,在误差接受范围内,研究双相介质是可行的;把稍作改进的镜像法应用于双相介质中,能够有效地处理瑞雷面波自由边界问题;通过详细分析双相介质瑞雷面波及体波在内的全波场的信息,对以双相介质为基础的地震波勘探有一定的指导作用。     

裂隙各向异性介质横波双折射特性研究

借助于有限元方法,对裂隙各向异性介质一雏四分量弹性波场中横波双折射的有关问题进行了研究。它们是裂隙方位与相速度,双折射与震源主频,双折射与裂隙丰度,双折射与能量分配,双折射与延迟时,双折射及其振幅谱等。这些问题对深入认识横波双折射特性,对理论和实际工作都有一定帮助。     

基于ADIPI的单程波三维VSP正演模拟方法

三维VSP正演模拟在三维VSP观测系统设计、三维VSP波场分析及三维VSP偏移成像方法研究中均有重要的作用。以三维单程波动方程为基础,采用基于正交方向加插值(ADIPI)的傅里叶有限差分(FFD)算子,实现了三维VSP正演模拟。克服了对三维算子的双向正交分裂所导致的方位各向异性误差,既保证了模拟精度,又提高了计算效率,还可根据实际需要单独模拟上行波或下行波。模型试验结果证明了方法的有效性。     

放射性核素在裂隙介质中迁移模型研究综述

指出了研究核素迁移问题的意义,归纳了裂隙介质的几种简化形式,概括总结了国内外对放射性核素在裂隙介质中迁移模式研究中的进展.提出了迁移数学模式的分类方法,并指出了目前迁移问题研究中存在的不足.     

裂隙-管道介质泉流域水文地质模拟试验

为了深入认识岩溶泉流域的流量过程影响因素和作用机制,以西南岩溶泉流域水文地质概念模型为原型,建立了实验室尺度下的裂隙-管道介质物理模型和数值模型。选取补给雨强、泉口直径、落水洞与开采井的位置作为影响因子,模拟了泉流量过程。试验结果表明:单次降水试验可以观测得到泉流量增加、平稳波动和衰减共三个阶段;同一泉口直径条件下,补给强度对泉流量衰减系数的影响极小,直径为3,4,5,7 mm的泉口直径对应的衰减系数分别为0.003 6,0.006 7,0.011 5,0.012 9,泉口直径与衰减系数呈近似对数关系;落水洞汇流水量占泉流量的比例远大于裂隙;2~3 mm是微小裂隙与较大裂隙的过渡区域,两种类型的介质泉流量过程差异明显;开采井分别位于泉口附近、模型中间位置和落水洞时,泉流量衰减系数依次减小。     

Equivalent Permeability and Simulation of Two-Phase Flow in Heterogeneous Porous Media

The problem of calculating equivalent grid block permeability tensors for heterogeneous porous media is addressed. The homogenization method used involves solving Darcy's equation subject to linear boundary conditions with flux conservation in subregions of the reservoir and can be readily applied to unstructured grids. The resulting equivalent permeability tensor is stable as defined relative to G-convergence. It is proposed to use both conforming and mixed finite elements to solve the local problems and compute approximations from above and below of the equivalent permeability, respectively. Comparisons with results obtained using periodic, pressure and no-flux boundary conditions and the renormalization method are presented. A series of numerical examples demonstrates the effectiveness of the methodology for two-phase flow in heterogeneous reservoirs.     

多孔介质渗流力学理论研究现状及发展趋势

多孔介质广泛存在于地下岩土结构、工程材料、生物体中,由于多孔介质微结构的复杂性,多孔介质渗流问题受到广泛关注。总结介绍了多孔介质渗流理论在国内外的研究现状,对多孔介质的描述、渗流问题的研究方法以及多孔介质中溶质运移、热量运移以及数值模型等方面的研究进展作了评述,并就多孔介质渗流理论在环境、农业、生态等领域的应用与发展方向提出看法。     

地质雷达数值模拟中有损耗介质吸收边界条件的实现

利用时间域有限差分(FDTD)法将麦克斯韦方程进行离散化,可以对地质雷达进行数值模拟。利用完全匹配层(PML)作为吸收边界条件可以有效地吸收向外的电磁波,从而大大提高了计算效率。对于有损耗介质的情况,采用扩张坐标系下改正的麦克斯韦方程,只要扩张变量满足一定的关系,同样可有效地吸收向外的电磁波, 这种吸收边界条件称为通用完全匹配层。     

TTI介质弹性波频率-空间域有限差分数值模拟

由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。     

利用高阶交错网格有限差分法数值模拟VTI介质井孔声场

横向各向同性（TI）介质是岩石地球物理中常见的一种现象,研究其井孔声场传播特征对声波测井理论以及为声波测井解释提供依据具有重要意义。针对具有垂直对称轴的横向各向同性（VTI）介质,根据柱坐标系条件下的弹性波波动方程,推导了速度-应力交错有限差分公式,采用时间二阶、空间十阶的交错有限差分算法对VTI介质中的井孔声场进行数值模拟。给出了在均匀介质中井孔声场不同时刻的波场快照,以及不同各向异性系数的VTI介质中的波场快照,计算了井轴上声源激发出的声波全波列波形。结果表明,在其他条件不变的条件下,VTI地层的各向异性系数的增大对横波的传播影响不大,但会使得纵波在纵向上的传播速度相对变小,径向上变化不大。各向异性系数的增大会使声波测井全波列首波信号时差变大,声波幅度略变小。     

水力旋转填料生物接触氧化处理城市污水的研究

通过对中试规模的新型水力旋转填料生物接触氧化法处理城市污水的试验研究,利用新型填料的独特构型,增强反应器的传质性能,探讨了接触氧化池对COD、NH+4-N的去除,及水力停留时间、气水比、温度对处理效果的影响。试验结果表明在水力停留时间为2~2.5h,气水比为(7~8)∶1,进水COD、NH4+-N浓度分别为130.1~366.5mg/L、11.3~19.2mg/L的条件下,两者的平均去除率为75.9%、65.4%,出水COD、NH4+-N浓度均已达到国家二级排放标准;气水比和水力停留时间对接触氧化池处理效果的影响较为明显;低温下接触氧化池对有机物、氨氮的去除仍具有较好的效果。     

基于CFP的岩溶管道流数值模拟研究——以桂林寨底地下河子系统为例

为了描述岩溶多重含水介质存在非达西流问题,研究管道流CFP模型在岩溶区数值模拟中的应用。在分析管道流CFP基本原理的基础上,通过建立概念模型算例(5个管道,上下临界雷诺数分别为2 000和4 000),探讨不同管道水文地质参数对出口流量及水流状态的影响,并对参数敏感性进行分析,最后将CFPM1模型应用于广西寨底岩溶地下河系统,探讨在实际应用中的管道流模型适用性。概念模型结果表明,参数敏感度从大到小依次为:管道直径、水力梯度、管道渗透系数、弯曲度及粗糙度,流量与雷诺数随管道参数增大而增大,水流状态从层流变为紊流。寨底岩溶管道流模型中3个观测井及地下河出口实测流量表明CFPM1能较好地模拟出岩溶地下水水位变化趋势。得出结论:CFPM1管道流模型允许岩溶管道与基岩含水层进行水流交换,能够较好地刻画岩溶区管道介质非达西流特征,但管道参数获取较为困难,且不能刻画岩溶管道形态变化。该方法具有一定通用性,可为实际应用提供指导。     

用随机介质模型方法描述孔洞型油气储层

孔洞型储集层具有强烈的非均质性,在空间分布上变化剧烈,规律性差,其地球物理参数(如速度、密度、弹性参数等)在空间上的变化很难用传统的层状介质模型进行刻画.本文借助随机模型方法,研究孔洞型储层中不同孔洞尺度、分布密度与模型参数的关系,建立反映实际孔洞油气储层介质地球物理参数空间分布统计特征的随机介质模型,以更方便地对孔洞型油气储层进行描述,并为碳酸盐岩孔洞型储层地震反射特征的研究提供有效的方法途径.