横向各向同性介质中地震波走时模拟

横向各向同性介质是地球内部广泛分布的一种各向异性介质.针对这种介质,我们对各向同性介质的最小走时树走时模拟方法进行了推广,推广后的方法可适用于非均匀、对称轴任意倾斜的横向各向同性介质模型.为保证计算效率,最小走时树的构建采用了一种子波传播区域随地震波传播动态变化的改进算法.对于弱各向异性介质,我们使用了一种新的地震波群速度近似表示方法,该方法基于用射线角近似表示相角的思想,对3种地震波（qP, qSV和qSH）均有较好的精度.应用本文地震波走时模拟方法对均匀介质、横向非均匀介质模型进行了计算,并将后者结果与弹性波方程有限元方法的模拟结果进行了对比,结果表明两者符合得很好.本文方法可用于横向各向同性介质的深度偏移及地震层析成像的深入研究.     

计算横向非均匀介质体波理论地震图的Maslov渐近法

本文运用Maslov渐近理论编写的二维横向非均匀介质中的理论地震图程序,与其它类型算法作了精确对比,结果表明,在层状介质模型中,本程序的结果无论振幅还是波形对比都与反射率法基本相同。对于二维横向非均匀模型,在射线理论的非奇点处,本算法与射线方法基本一致,在射线理论的奇异点处,Maslov方法消除了射线理论所固有奇点,提高了计算精度。     

计算横向非均匀介质体波理论地震图的Maslov渐近法

本文运用Maslov渐近理论编写的二维横向非均匀介质中的理论地震图程序,与其它类型算法作了精确对比,结果表明,在层状介质模型中,本程序的结果无论振幅还是波形对比都与反射率法基本相同。对于二维横向非均匀模型,在射线理论的非奇点处,本算法与射线方法基本一致,在射线理论的奇异点处,Maslov方法消除了射线理论所固有奇点,提高了计算精度。     

层状介质中有限震源引起的地面运动的计算对点源情形的拓展

从分层均匀介质中地震波传播的基本理论出发, 参考已有的利用广义反射透射系数矩阵方法计算分层均匀介质中点源引起的地面运动的方法和程序, 考虑了任意几何特征的有限震源及有限震源的震源机制随时间和空间发生变化的可能性, 并考虑了并行计算的发展方向, 本文对点源情况下的计算流程进行了改进, 重新设计编写了计算程序GRTMatSyn, 使其不但适应于点源的情形, 也适应于任意几何形状的、 有限的、 震源机制随时间和空间变化的有限震源的情形． 为了验证该程序的可靠性和计算效率, 设计了必要的有限震源模型和观测点分布, 计算了地面运动, 并与已有的被广泛认可的计算有限平面断层引起的地面运动的程序CompSyn的计算结果进行了对比分析． 结果表明, GRTMatSyn的计算结果可靠、 计算效率较高．      

运用改进的地震定位交切法重定位2001—2010年云南地震

传统地震定位方法利用震源轨迹确定震源位置,但基于均匀或横向均匀介质模型必然导致定位误差。为此对传统方法进行改进,发展适用于三维复杂地壳速度模型的地震定位交切法。利用最小走时树射线追踪技术,以离散方式准确计算三维复杂地壳速度模型中的震源轨迹,将震源定位于震源轨迹交汇的密集点。将该方法应用于云南地区地震重定位,得到较高定位精度。     

复杂介质结构中折射界面的哈格多恩原理波前成像

在城市活断层探测中 ,浅层结构常常表现为强烈的非均匀性 ,界面横向强烈起伏 ,层内速度变化较大 ,传统的基于平界面均匀层模型的折射资料处理方法不能适用。研究开发能应用于复杂介质结构中折射资料处理的方法就显得十分必要。文中基于惠更斯原理 ,用波前扩张法对波场作正演计算 ,根据哈格多恩折射波前成像原理 ,在lecomte算法和Hole有限差分计算程序的基础上 ,开发出 1种复杂介质结构中折射资料的处理方法与软件 ,并用此方法处理了福州城市活断层折射探测试验中在义序完成的 2条折射剖面资料。结果表明 :探测区浅层为 3层结构 ,分别为盖层、强风化层和基岩。基岩顶界面的埋深约为 5 8～ 5 2m ,盖层P波速度变化较大     

地震定位的广义反演方法及计算程序

本文介绍一种利用近震P、S波到时,在任意垂直非均匀介质中利用广义反演理论确定地震事件空间位置的概率分布的方法及计算程序。做为一个例子,给出了用这个方法对滦县近场地震观测资料的数值计算结果,并与程序HYPO81所得的结果作了比较。     

基于Pade逼近的Cole-Cole频散介质GPR有限元正演

针对Cole-Cole频散介质中的复介电常数是jω的分数次幂函数,传统的时域有限元法难以离散及计算时间域分数阶导数,本文采用Pade逼近算法将含有时间分数阶导数的Cole-Cole频散介质电磁波方程推导为一组整数阶辅助微分方程,提出了一种适用于Cole-Cole频散介质的GPR有限元正演模拟算法.在复数伸展坐标系下,通过在频率域Cole-Cole频散介质电磁波方程中引入2个中间变量,并将其变换到时间域,从而以变分形式将PML边界条件加载到Cole-Cole频散介质GPR有限元方程组中,并给出了详细的求解公式.在此基础上,编制了基于Pade逼近的Cole-Cole频散介质GPR有限元正演程序,利用该程序对均匀模型进行计算,并与解析解进行对比,验证了本文构建的GPR有限元正演算法的正确性和有效性.设计了一个复杂Cole-Cole频散介质GPR模型,利用本文构建的GPR有限元正演算法进行模拟并与非频散介质模型的模拟结果进行对比,分析了电磁波在Cole-Cole频散介质中传播衰减增强、子波延伸,分辨率降低等传播规律,有助于实测雷达资料更可靠、更准确的解释.模拟结果表明,基于Pade逼近的GPR有限元正演算法可用于复杂Cole-Cole频散介质结构模拟,且具有较高的计算精度.     

宽角反射地震波走时模拟的双重网格法

在研究地壳结构的人工源宽角反射地震资料解释中,常规宽角反射波走时和射线路径计算大都假定地壳模型为层状块状均匀介质.为了逼近实际地壳结构模型,要求模型尺度较大,为了提高地震资料解释的可靠性,须减小模型离散单元的尺寸,但同时计算量大大增加,使资料解释的效率较低.为此,本文尝试同时提高宽角反射地震资料解释效率和可靠性的方法,即使用双重网格计算宽角反射地震波走时和射线路径的最小走时树方法.双重网格法在均匀介质内部仅计算大网格节点,在速度变化点、震源点和检波点区域,同时计算小网格节点;在界面边界点使用比介质内部节点更大的子波传播区域.模型计算结果表明,对于大尺度的层状块状均匀介质模型,在保证精度的条件下,本文所提出的双重网格射线追踪方法的计算效率比单网格方法显著提高.     

用传输线模型法计算非均匀介质中的电场响应

介绍了一种非均介质中电场响应的数值模拟方法.传输线模型法是以惠更斯原理为基础,通过建立麦克斯韦方程中电场和磁场与传输线方程中电压和电流的等效关系,按时间步进的办法来模拟电磁场响应的一种电磁场数值计算方法.从一维介质中电场传播的物理过程出发,导出了点电流源情形下地质介质中电场响应迭代方程;编程计算了非均匀介质中的电场响应,其结果表明,用这种方法对非均匀介质中的电场响应进行数值模拟是行之有效的.     

余震震源机制变化的原因

大震以后,跟随着一系列的余震,其震源机制有的和主震一致,有的则大不一样,这表现在P波初动的反向,反应了震源应力场的变化。 余震震源机制变化的原因,既不是岩块的错动过头,也不是液体流出的下陷,而是发生大震的滑动岩块对其前后邻接岩块产生的转换应力场所引起的剪切破裂。     

徕远广场结构设计地震波的确定

选择合适的人造地震动时程是高层建筑采用时程分析法计算地震作用的基础。介绍了徕远广场结构设计地震波的确定方法及过程。鉴于各条地震波输入进行时程分析的结果的离散性,提出可以采用场地土层反应分析得到的地表加速度时程及按照《高层技术规程》设计规准谱拟合人造波作为结构设计地震波。     

西南天山柯坪塔格前缘断裂带东段晚第四纪活动特征

柯坪塔格前缘断裂东段是柯坪推覆构造系前缘的一条活动断裂,野外调查获得了其晚第四纪错断洪积扇、冲沟阶地面的证据,实测了变形地貌面上的断层陡坎,分析了断层的形变量,通过采样测年估算了断层的缩短速率。由7个观测点的断层陡坎剖面测量,计算了观测点处断层的水平缩短速率,结果表明,断裂弧顶部位的五道班、三间房以西及其大山口道班附近,断层错断了Ⅰ级和Ⅱ级洪积扇(阶地)。断层在这些地点最新活动强烈,水平缩短平均速率全新世以来为0·35~0·44mm/a,更新世晚期末以来为0·16~0·30mm/a,而在非弧顶部位的巴楚磷矿、三岔口以北及大山口北断层只错断了更新世晚期Ⅲ级洪积扇,且水平缩短速率较小,断层水平缩短平均速率更新世晚期以来为0·05~0·07mm/a     

蒙脱石的脱水作用对断层摩擦本构行为的影响

利用高温双轴摩擦装置，研究了含蒙脱石的断层带在不同温度下摩擦滑动的速度依赖性，以期了解脱水作用对摩擦行为的影响。结果表明，断层带摩擦强度随温度而升高，而速度依赖性较为复杂，以１．４ｕ／ｓ为界，室温和１００℃时，低滑动速率下表现为微弱的速度弱化，高滑动速率下则表现为速度强化；２００℃时均为速度强化；３００℃时高滑动速率下仍为速度强化，但低滑动速率下转变为速度弱化；４００℃以上，均为明显的速度弱化。摩擦行为的变化与脱水过程及相应的断层物质变形方式的变化密切相关     

前进中的青海省工程地震研究院

青海省工程地震研究院是青海省地震局下属具有独立法人资格的全民所有制科技开发企业,是地震工程、岩土工程和环境地质领域从事科研、生产、人才培养的经济实体,具有建设工程地震安全性评价、桩基检测甲级资质,岩土工程勘察、咨询乙级资质,地质灾害危险性评估丙级资质。经过十余年的发展,从队伍建设、能力建设作业绩等方面,都取得了长足的进步。     

建筑物基础隔震技术的进展

介绍了建筑物基础隔震技术的原理。基础隔震技术是通过设置在结构物底部与基础顶面之间的隔震消能装置,来增加结构的变形能力和滞变阻尼。目前结构基础隔震体系按隔震机理不同划分类型,主要有橡胶垫隔震体系、滑动摩擦隔震体系、复合基础隔震体系等类型。总结了基础隔震技术应用和理论研究的发展,指出隔震技术的未来发展趋势。     

SNREI地球模型负荷勒夫数数值计算的新进展

基于Ｌｏｎｇｍａｎ和Ｆａｒｒｅｌｌ的表面负荷形变理论，深入研究计算ＳＮＲＥＩ弹性地球模型表面负荷勒夫数的数值计算方法，取得新的进展：采用了一种新的代换方法─ｒ￣ｎ因子法，比传统的无量纲化法简单易行，且能达到很高的计算精度；确立了检验表面边界条件、微分方程积分线性性和计算结果的收敛性，用于评价数值计算结果精度的有效准则；给出了ＰＲＥＭ地球模型的一组高精度计算结果．     

水对断层摩擦滑动稳定性的影响

在不同含水条件下研究了光面、含脆性物质夹层和含延性物质夹层 3类断层的摩擦滑动行为 ,认为水对断层滑动稳定性和摩擦滑动的速度依赖性有重要影响。建立了含有岩性、围压 (深度 )和含水量 3因素的断层活动稳定性的模式 ,认为脆性断层更多地表现为速度弱化 ,少量水的存在易于出现黏滑 ,高孔隙压时会表现为稳滑 ;半脆性或半延性的断层在少量水或较低孔隙压的情况下会表现出速度弱化 ,也就可能出现黏滑 ,但孔隙压较高后 ,会表现出速度强化 ,滑动也就会是稳定的 ;延性断层多为速度强化 ,不出现黏滑。提出水的进入降低了系统的稳定性 ,水的存在增加了系统的稳定性。这是考虑水或一般流体对岩体或断层活动稳定性影响时需要考虑的两个方面     

储层砂岩声波速度预测

本文主要基于Ｇａｓｓｍａｎｎ方程和经验规律，提出了孔隙流体替代和孔隙度改变时对砂央地震波速度变化的估计，以及直接利用岩石矿物和孔隙流体的弹性性质计算砂岩地震波速度方法，利用已知的岩芯，测井或地震数据，运用这此方法，可合理地对储层砂岩地震波速度进行预测。     

Biot理论的唯象修正对S波特性的影响

将复模量引入Biot方程后，在一维条件下通过S波的波动方程研究了S波的传播特性，S波的数值分析显示在频率域或温度域上都能获得热弛豫衰减峰和Biot衰减峰. 在频率域上由于温度的变化引起两峰相向位移,在温度域上,因频率的变化也发生相对移动. 随着温度和频率的不断提高，两峰发生叠加，叠加后两峰互换位置. 低频或低温段的热弛豫峰移到了高频或高温段，高频或高温段的Biot峰移到了低频或低温段.由于两峰的衰减机制不同，导致S波波速随频率或温度变化规律的复杂性. 这些规律已部分被共振实验所证实，证实该理论模型具有实验基础.