上地幔各向异性介质对固体潮及负荷潮的影响

讨论了上地幔各向异性介质中的潮汐运动方程,根据Deiewonski提供的地球模型参数,利用经典的Runge－Kutta数值积分方法,计算了固体潮勒夫数和负荷勒夫数．结果表明,考虑上地幔介质各向异性与否对固体潮勒夫数的影响较小（约为0．06％）,而对负荷勒夫数的影响较大（2．5％）,进一步说明了中低阶负荷勒大数对上地慢介质特性的敏感性．     

岩溶管道介质空隙率变化（视表征体元）及分形研究

本文参照Ｊ．Ｂｅａｒ对于多孔介质提出的表征体元,以杭州瑶琳洞为例给出了岩溶管道介质空隙率变化曲线,并应用分形理论研究空隙率变化规律。这种作为基础理论的研究,具重要的岩溶水文地质意义。     

利用量子力学势能数据预测温度、压力、盐度和沉积物孔径对甲烷水合物形成和分解的影响

介绍一个预测不同温度、压力、盐度和沉积物毛细管孔径条件下甲烷水合物 溶液 气体多相平衡模型。该模型以Van der Waals和 Platteeuw热力学模型、量子力学从头算粒子相互作用势能、DMW 92状态方程和Pitzer电解质理论为基础,能在很宽广温压范围内预测温度、压力、盐度和毛细管力对甲烷水合物形成和分解的影响。通过对比本模型的预测结果与实验数据,可知本模型能够准确地预测海水和多孔介质中甲烷水合物的相平衡条件。对于一定盐度下多孔介质中甲烷水合物的形成温压条件的在线计算可浏览: www.geochem model.org/models.htm。     

川南、滇北地震台网区域壳幔介质分层结构的初步研究

本文利用川南、滇北21个地震台短周期地震仪的深源远地震波形作为观测资料,通过调整地壳和上地幔介质模型并计算相应的综合地震图的方式,直至获得与观测波形可相比拟的综合地震图。据此,推断该地震台网区域的地壳中下部由一组高速和低速相间的薄层所组成,且大部分台站下的地壳厚度达到或超过54km。此外,文中还研究了该区域壳幔介质的吸收特点。     

一种适于近场理论地震图宽频域计算的矩阵分解算法

理论地震图方法在地震震源过程的研究中得到了广泛的应用。为了研究震源过程的细节,必须利用近场地震资料的高频信息,但是在应用Haskell矩阵法计算近场理论地震图时,格林函数的高频成分的数值不稳定性是计算宽频域理论地震图的一个基本困难。因此,目前仅有不超过10Hz的算例。本研究采用Haskell矩阵的一种新的分解组合形式,在数值计算时有效地避免了在计算过程中数值结果的溢出,实现了近场理论地震图的宽频域计算。同时,由于利用了矩阵运算中的解析关系,减少了运算的次数,从而提高了计算的速度。根据本算法建立的Fortran程序,在Univoc-1100计算机上进行了数值检验。结果表明,计算得到的格林函数至少在0—40Hz频率域范围内仍保持良好的数值稳定性。本文给出的这种Haskell矩阵分解组合形式的矩阵元素简单,且具有一定的对称性。对于这种分解所包含的物理意义,将有待进一步深入研究。     

关于非均匀介质中电导率按深度变化时的势及其视电阻率函数

本文共分六部分内容。在前言中,简述了电测深方面有关研究的情况,尤其叙述了E.S.Sampaio（1976）的结果。 因本文采用的方法不同,从而求解了适合要求的微分方程定解问题,获得了电导率以 （其中α、σ0、σ1均为正实数,β为实数）按深度变化时势的分布。 其次,当β=0时,对电导率σ=（σ0+σ1z）α按深度变化时势的分布,应用Bessel函数的积分表达时,明显地改进和推广了E.S.Sampaio（1976）结果中的重要部分。 最后,结合实际作出了第一类、第二类标准化视电阻率函数。这种视电阻率函数,推广了Wenner和Schlumberger分别作出的视电阻率函数的有关部分,对电阻率法和激发极化法,都有其理论价值和实际意义。     

介质的磁化率与包裹体浓度的关系

在第五卷第一期的地球物理学报上,我们曾讨论介质的磁化率与包裹体形状的关系,但是对于介质磁化率与包裹体浓度的关系没有详细讨论。在本文中,将根据巳知的材料,对这个问题加以讨论,作为“介质的磁化率与包裹体形状的关系”一文的补充。     

含气介质临界点、流体和骨架弹性参数的数值计算方法

目前临界孔隙度及孔隙介质的研究主要为测试分析的实验方法,如何运用数值计算方法求取孔隙流体介质临界点和组分弹性参数一直是备受关注的课题.本文提出了求取孔隙介质临界点、孔隙内流体和骨架弹性参数的数值计算公式及方法,并结合含气样品测试数据实现了这种数值计算.文章首先给出了3个包含孔隙度的线性方程,通过有机组合每个线性方程中的两个系数得到求取有关弹性参数的数值计算公式.然后,详细阐述了计算的方法和步骤以及需要注意的问题,把室内含气砂岩样品实测的综合介质密度、纵波速度和横波速度作为数值计算的输入数据,求得了临界点和流体及骨架等相关弹性参数的具体数值.通过比较分析本文数学方法求得的计算数据与实验方法测得的实测数据,表明了该数值计算公式的正确性和数值实现方法的有效性.     

饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法

本文基于Biot的饱和多孔介质本构方程,考察具有辐射阻尼的外行球面波,推导了饱和多孔介质三维黏弹性人工边界的法向和切向边界方程;在已有的饱和多孔介质二维显式有限元数值计算方法基础上,提出该理论的三维方法,并开发了实现该三维方法的有限元程序.算例表明饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法具有较好的精度和稳定性.     

横向各向同性介质弹性波多分量叠前逆时偏移

随着油气勘探程度的提高,隐蔽油气藏在增储上产方面起到了重要作用,因此发展基于各向异性介质的多分量偏移方法是非常必要的.本文基于横向各向同性（VTI）介质,从二维弹性波速度\|应力方程出发,通过在时间上的二阶差分和空间上的交错网格高阶差分对方程进行离散,得到弹性波交错网格高阶差分的多分量逆时偏移算子.在激发时间成像条件的应用过程中引入Poynting矢量进行成像并消除逆时偏移所引起的低频干扰,在此基础上实现了VTI介质中二维弹性波叠前多分量逆时深度偏移.理论模型的偏移处理表明,该方法能够对地层进行准确成像,并可以消除逆时偏移所引起的低频噪声.