天然气水合物和游离气饱和度估算的影响因素

讨论了不同水合物胶结类型的流体饱和多孔隙固体中地震波的衰减情况,分析了估算天然气水合物和游离气饱和度影响因素.结果表明, 地层孔隙度、纵波速度模型和弹性模量的计算方法是影响反演水合物和游离气饱和度的关键因素.含水合物地层的吸收与水合物胶结类型密切相关,当水合物远离固体颗粒,像流体一样充填在孔隙时,品质因子出现负异常,而当水合物胶结固体颗粒影响骨架的弹性性质,其品质因子出现正异常.根据布莱克海台地区164航次995井的测井资料,分别应用低频和高频速度模型估算了水合物和游离气饱和度.由低频速度模型得到的水合物饱和度(占孔隙空间的)10%～20%,游离气饱和度(占孔隙空间的)05%～1%;而由高频速度模型得到的水合物饱和度(占孔隙空间的)5%～10%,游离气饱和度(占孔隙空间的)1%～2%.     

改进的等效半空间法及瑞雷波频散曲线反演

层状模型的理论频散曲线表明,层状模型表面瑞雷波同一频率对应的相速度具有多值特点,这就是所谓的瑞雷波多阶性. 面波的多阶性导致了实测频散曲线的复杂性,当地层中存在软弱夹层时这一问题尤为突出,从而给实测频散曲线的解释带来了困难. 本文从工程实用的角度出发,根据等效半空间理论,提出了一种计算理论频散曲线的新算法——改进的等效半空间法,避开了面波多阶性这一复杂问题. 由此方法计算得到的面波相速度并不同于某个具体阶数的面波,而是对多阶面波的综合反映. 按照改进的等效半空间法编制了拟合反演程序,并在工程中进行了运用,取得较好的效果.     

全球地幔垂直流动速度研究

用高分辨率地震体波速度成像以及相关的地球物理资料,计算地幔垂直流动形式及流动速度,得到全球地幔流垂直运动模式.从全球尺度来看,地幔流基本可划分为以下几个区域：欧亚大陆—澳大利亚、北美洲—南美洲为两个大规模下降流区域,西印度洋—非洲及大西洋、中南太平洋及东太平洋为两个大规模地幔上升流区域.地幔上升流起源于核幔边界,主要表现在地幔中部和上地幔下部.地幔垂直流动速度约每年1～4cm.地幔流动对地表板块运动、海洋中脊和中隆、俯冲带和碰撞带的分布起着控制作用.地幔上升流与地表现代热点有密切关系.从东亚尺度看,地幔流大体分为三个区域：东亚边缘裂谷系和西太平洋边缘海为上升流、西伯利亚地幔深度表现为物质下降流、青藏高原—缅甸—印度尼西亚特提斯俯冲带地幔下降流,这三个区域地幔流动与地表的西太平洋构造域、亚洲构造域和特提斯构造域相吻合.勾勒出南海地区构造特征：从上到下的大体结构是上部呈“工"字型、中间为圆柱型、底部呈盾形的地幔上升流.     

无单元法用于地震波波动方程模拟与成像

波动方程方法是解决地震正反演问题的基本工具之一.无单元法作为一种新兴的偏微分方程数值计算方法,已经在材料力学、热传导等领域取得了显著的成功.由于抛弃了单元的概念及采用滑动最小二乘的拟合方法,使得无单元法具有前处理简单、精度高、独立变量解高次连续等优点.本文首先介绍无单元法求解波动方程的原理,指出影响其精度的主要因素.在算例的基础上详细讨论了无单元法用于实际波动问题的效果,并进一步尝试利用无单元法进行地震波数值模拟和反演成像的研究.模型计算的结果表明,无单元法能够较好的处理地震模拟和成像问题,精度和稳定性是令人满意的.     

方位各向异性粘弹性介质波场数值模拟

当地震信号通过复杂地球介质时,地层除了表现为各向异性,还表现为内在的粘弹性特征.因此,为准确描述地震波在地球介质中的传播特征,理想的地球介质模型应该能够模拟岩石的各向异性特征和衰减特征.本文给出了各向异性粘弹性介质模型的波动方程及其差分格式,并利用有限差分法实现了地震波波场数值模拟.结果表明了该介质模型中地震波场特征与各向异性主轴方位和介质的粘滞性参数之间的关系.     

在SeisDps实时系统中用零点、极点计算ML的过程

在SeisDps地震处理软件的实时系统中,不是用系统的幅频特性计算震级ML,而是利用台站的零点、极点和增益,将记录的数字地震波形仿真成伍德-安德森地震仪记录的波形,计算具有3 大优点的ML震级,本文主要介绍用零点、极点和增益计算ML的过程。     

分析云南地震资料给出的近震震级ML和面波震级Ms

云南省区域数字地震遥测台网自1999年建成以来,记录了大量地震波形资料,从2000年1月至2006年2月,云南省及邻区共发生ML≥4.0地震168次,对这些地震的震级重新进行了测定,结果表明,对于云南省的区域地震,近震的体波震级ML和面波震级MS之间不满足昆明地震台网多年使用的MS=1.13ML-1.08这一换算公式,且发现ML也不一定大于MS。     

非一致激励条件下工程场地地震动相干函数的数值模拟--Ⅰ分析原理和方法

本文将随机振动的虚拟激励原理与工程波动理论散射问题的求解方法相结合,建立了开放系统中非一致激励条件下工程场地地震动随机场的数值模拟方法。该方法将随机输入下的波动分析问题转换为多个虚拟激励下的确定性波动分析组合问题,从而可以方便地获得场地波动观测量之间的谱密度矩阵,进而计算给出工程场地的地震动相干函数。本文还用数值模拟的办法对所提出方法的合理性和稳定性进行了探讨。     

波浪作用下成层海床孔隙水压力响应的简化分析

采用不排水条件下孔隙水压力发展模式,作为Terzagh i一维固结方程中考虑波浪循环作用所引起的孔隙水压力源项,对于成层海床建立了推广的一维动力固结方程,运用数理方程中的分离变量法与G reen函数求解了成层海床在波浪作用下残余孔隙水压力的发展规律,进而对成层海床的液化势进行了评判。对比计算与分析表明,海床表层土的渗透性及其厚度对于海床的整体抗液化性能具有显著的影响,低渗透性的表层导致海床孔隙水压力的显著积累,此时表层置换法是防治液化的有效途径。     

3D viscous-spring artificial boundary in time domain

After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.