层状各向同性多孔介质中井间地震的数值计算

依据Ｂｉｏｔ理论,提出计算理论地震图的相应传播矩阵算法,并针对井间地震震源和检波器的不同配置进行数值计算,结果表明,本文算法高频可达１０００Ｈ２以上．文中算例说明利用震源和检波器相对于储油层的不同位置,可从反射震相记录波形上的差异,判定储油层存在与否．     

层状多孔介质的地震电磁异常数值模拟方法研究

本文针对电激发模型在各向同性分层介质中利用传播矩阵方法对电磁异常产生过程的数值模拟.层状模型中一部分为饱和多孔介质,另一部分为普通线弹性介质,数值模拟过程中将固体场作为主要场单独进行计算,而将其形变作为等效流体源计算耦合电磁场.模拟结果显示,地震波在经过饱和多孔介质边界的时候会产生电磁场,而其在流体层内部传播的过程中则未产生电磁异常.模拟结果可以用于解释一些伴随地震产生的电磁异常,包括地震发生时以及地震波到达时所产生的电磁异常现象,另一方面对于前震及其它一些无感地震所产生的电磁异常的解释也具有一定的参考价值.      

横向各向同性多孔介质中的地震波传播

基于各向异性多孔介质中的广义Ｂｉｏｔ理论,导出了横向各向同性多孔介质中波传播的特征方程．指出在多孔介质中有４种类型的频散和耗散波传播：准纵波ＱＰ１（快纵波）、准纵波ＱＰ２（慢纵波）、准横波ＱＳＶ和横渡ＳＨ．文中给出了４种波速度的解析表达式．数值计算频率曲线和衰减曲线与Ｓｃｈｍｉｔｔ（１９８９）用均值处理得到的结果类似．还给出了波传播过程中３种类型准体波之间的耦合系数（或称转换系数）．     

横向各向同性多孔介质中的地震波传播

基于各向异性多孔介质中的广义Ｂｉｏｔ理论，导出了横向各向同性多孔介质中波传播的特征方程．指出在多孔介质中有４种类型的频散和耗散波传播：准纵波ＱＰ１（快纵波）、准纵波ＱＰ２（慢纵波）、准横波ＱＳＶ和横渡ＳＨ．文中给出了４种波速度的解析表达式．数值计算频率曲线和衰减曲线与Ｓｃｈｍｉｔｔ（１９８９）用均值处理得到的结果类似．还给出了波传播过程中３种类型准体波之间的耦合系数（或称转换系数）．     

弹性多孔介质中应力—孔隙压耦联理论在新丰江水库地震研究中的应用

本文用二维有限元法,将弹性多孔介质中应力-孔隙压耦联理论应用于新丰江水库诱发地震机制的研究。计算中,不仅考虑了均匀、非均匀模型的耦联,还考虑了断层的性质、产状对耦联作用的影响及库水位随时间的变化。 结果表明:附加应力与孔隙压耦联作用是诱发水库地震的重要因素。新丰江水库地震不仅与库区存在的断裂、高密度基底有关,还与库水的加、卸载方式有关。     

横观各向同性层状半空间中的弹性位错

本文在柱向量函数系下,利用传播矩阵法求解了层状横观各向同性半空间由内部点源位锚引起的变形;对六个基本点源位错,以等价体力法推出了横观各向同性情形下的点源函数,并且给出了内部任意剪切位错源引起的地表位移的积分表达式。为研究地球的层状结构,特别是其上部的横观各向同性对地表的地震位移、应变以及倾斜场的影响提供了计算公式。     

横观各向同性层状半空间中的弹性位错

本文在柱向量函数系下,利用传播矩阵法求解了层状横观各向同性半空间由内部点源位锚引起的变形;对六个基本点源位错,以等价体力法推出了横观各向同性情形下的点源函数,并且给出了内部任意剪切位错源引起的地表位移的积分表达式。为研究地球的层状结构,特别是其上部的横观各向同性对地表的地震位移、应变以及倾斜场的影响提供了计算公式。     

横向各向同性介质中的地震定位

地震定位是地震学中最基本的问题之一.地震定位的精度和参考模型的准确性密切相关.目前无论是区域地震定位还是全球地震定位,所使用的参考地质模型几乎都是各向同性的,但大量的观测结果表明地球内部介质广泛存在着各向异性.     

二维层状介质中衍射波的计算

本文讨论了衍射波的计算,它是文献[5]的继续和发展,是为了应用射线法解决二维地震人机联作工作的一部分。根据波的物理性质,引入了等效炮点、等效检波点和等效界面等概念,将Berryhill导出的衍射波计算公式推广到倾斜界面和楔形结构等情形。这一工作与文献[1-3]中的工作相结合,在实用上解决了二维地震人机联作解释中的正演计算方法问题。     

层状孔隙介质震电波场数值模拟高频不稳定性问题的解析处理方法

自从发现震电现象以来,众多学者进行了相关研究.其中,Pride提出了一套描述流体饱和孔隙介质中震电波场的耦合与传播的宏观控制方程组,该方程组后来被广泛地应用到相关的震电研究中.Chen发展了一套广义反透射系数方法并将其应用到层状介质合成地震图的研究当中,该方法数值计算效率高并且可以处理带弯曲界面层状介质这种复杂模型.基于Pride的震电波场控制方程组,我们将Chen的广义反透射系数方法推广应用到层状孔隙介质中震电波场的数值模拟研究中,但是在数值计算过程中发现,当含源层的厚度相对于地震波波长较大时（即高频情况）,会出现数值计算的不稳定,此即为高频不稳定性问题.针对高频不稳定性问题,一种自然的处理方法就是在原来的含源层中插入两个虚拟界面,构造出一个新的含源薄层,但是这会带来一些额外的计算量,此外,由于虚拟含源薄层的厚度是有限的,必须针对具体模型参数设定一个合适的厚度值.高频不稳定性问题同样存在于层状介质合成地震图的数值计算过程中,Chen提出了一种解析的处理方法,即在原含源层内引入一无限薄的虚拟含源薄层,通过解析的方法解决高频不稳定性问题,该方法不会降低计算效率且适用于任意参数模型.本文首先对层状孔隙介质中的震电波场数值计算公式进行分析,指出源项积分中的指数增长因子是导致高频不稳定性问题的根本原因;其次将Chen在合成地震图数值模拟研究中采用的解析处理方法推广到震电波场研究中,得到了适用于数值计算的公式;然后给出数值算例,并针对一个含源层过厚的模型,比较了自然处理方法和解析处理方法,两种方法得到的结果具有相当好的一致性,而解析处理方法计算效率更高,证实了本文给出的解析处理方法在解决层状孔隙介质震电波场数值模拟的高频不稳定性问题方面的有效性.     

二维层状介质中衍射波的计算

本文讨论了衍射波的计算,它是文献[5]的继续和发展,是为了应用射线法解决二维地震人机联作工作的一部分。根据波的物理性质,引入了等效炮点、等效检波点和等效界面等概念,将Berryhill导出的衍射波计算公式推广到倾斜界面和楔形结构等情形。这一工作与文献[1—3]中的工作相结合,在实用上解决了二维地震人机联作解释中的正演计算方法问题。     

井孔声场的本征模及其在分层介质井孔声场计算中的应用

研究无穷大介质中的井孔本征模问题．从有限半径的波导出发,分析半径无限增大的极限情况,证明无穷大介质中的井孔有离散和连续两类本征模,组成正交完备基．进而应用它们分析计算分层介质中井孔的声场,算得孔外分层介质在孔内产生的反射波．     

井孔声场的本征模及其在分层介质井孔声场计算中的应用

研究无穷大介质中的井孔本征模问题．从有限半径的波导出发,分析半径无限增大的极限情况,证明无穷大介质中的井孔有离散和连续两类本征模,组成正交完备基．进而应用它们分析计算分层介质中井孔的声场,算得孔外分层介质在孔内产生的反射波．     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。