模拟地震仪器
abaqus如何施加地震波
施加地震波： 1 *amplitude,name=amp,input=seismicdata.dat 输入地震波 2 *boundary,type=acceleration,amplitude=amp施加荷载 方法：module选load，在tools-----amplitude-----creat默认的continue在Edit A mplitude里面输入时间和加速度，点OK。点creat boundary condition，涌现对 话框creat boundary condition，选择acceleration/angular acceleration,continu e---选择要施加的边界---done----涌现对话框edit bondary condition对话框，在 amplitude里选择你所定义的时间和加速度。点ok就完工了。 在网上查了些方法： module选load，在tools-----amplitude-----creat默认的continue在Edit Amplitude 里面输入时间和加速度，点OK。点creat boundary condition，出现对话框creat boundary condition，选择acceleration/angular acceleration,continue---选择要施加的边界---done----出现对话框edit bondary condition对话框，在amplitude里选择你所定义的时间和加速度。点ok就完工了。 这是在CAE里输入地震波的方式，我用的方法是直接在inp文件里加地震波的。 首先在CAE里建好模型，定义两个分析步。 第一个分析步是加自重，采用线性加载的方式。 (a) 加载方式：ABAQUS在施加Gravity时，默认为Instantaneous（瞬时加载），如果把结构自重以瞬间加载方式加到结构上，相当于对结构施加了一个脉冲荷载，会引起结构在竖向的振动，在不考虑结构阻尼的情况，这种振动会一直持续下去。如果是混凝土结构，这种竖向振动也会造成混凝土受拉损伤，所以这种加载方式不太合理。 (b)新建加载方式：创建一个新的Amplitude，Type=smooth tpye，0时刻Am=0，然后再选择一个0.5s~1s时刻，Am=1，在这个区间内线性插值，实现幅值从0到1。这种方式加载要优于上述瞬时加载，但是在起初的0.5s（或者1s，即smooth tpye中设置的终点时间）内计算结果是不准确的，所以要把这部分的计算结果剔除，剔除方法就是，创建2个step，第一个step主要分析自重作用，待自重稳定后开始第二个step地震时程反应分析。 第二个分析步就是加地震波。 输入地震波有两种方法： 1、在如下位置加入下面加黑的字体部分。格式如下：时间,地震波,时间,地震波,时间, 地震波,时间, 地震波…………每行8个数据（我下到的地震波文件是不带时间的，自己用C++处理了一下）。%%%%%%%%%%%%%%%%%%%%%% *End Assembly *Amplitude, name=Amp-1 0.005, -7.5e-08, 0.01, -3.55e-07, 0.015, -7.03e-07, 0.02, -4.53e-07 0.025, 1.82e-06, 0.03, 7.01e-06, 0.035, 1.5e-05, 0.04, 2.49e-05 0.045, 3.54e-05, 0.05, 4.5e-05, 0.055, 5.2e-05, 0.06, 5.5e-05 ………………
模拟地震波传播可视化
模拟地震波传播的可视化研究 摘要：实验中选取了与地壳平均波速相近的光学玻璃作样品，利用动态光弹的成像系统，来观测波在光学玻璃及波从光学玻璃透射到水中的传播过程，并记录0～50μs内的波的传播过程，以此来模拟地震波在地壳中反射、透射等传播行为。 abstract： the experiments selected optical glass similar with average velocity， and used imaging system of dynamic photoelasticity to observe the communication process of light in optical glass and light refraction from optical glass to water， and record the wave transmission during 0～50μs，for simulating reflection and transmission of seismic wave in crust. 关键词：地震波；动态光弹；反射；透射 key words： seismic waves；dynamic photoelastic；reflection；transmission 中图分类号：p315.3+1 文献标识码：a 文章编号：1006-4311（2013）04-0297-02 0 引言 在地震勘探中，通常是通过检波器来记录地下地震波带来的信息，根据相应的数学和物理模型进行复杂的计算机处理以获得地下的构造情况，虽然地震勘探的相关理论有很大的发展，但是由于理论结果难以获得，并且对于复杂形状的结构，解析方法变得相当繁
Midas地震波的选取方法
地震波的选取方法 建筑抗震设计规范(GB 50011-2001)的5.1.2条文说明中规定，正确选择输入的地震加速度时程曲线，要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。 频谱特性可用地震影响系数曲线表征，依据所处的场地类别和设计地震分组确定。这句话的含义是选择的实际地震波所处场地的设计分组(震中距离、震级大小)和场地类别(场地条件)应与要分析的结构物所处场地的相同，简单的说两者的特征周期Tg值应接近或相同。特征周期Tg值的计算方法见下面公式(1)、(2)、(3)。 加速度有效峰值按建筑抗震设计规范(GB 50011-2001)中的表5.1.2-2采用。地震波的加速度有效峰值的计算方法见下面公式(1)及下面说明。 持续时间的概念不是指地震波数据中总的时间长度。持时T d的定义可分为两大类，一类是以地震动幅值的绝对值来定义的绝对持时，即指地震地面加速度值大于某值的时间总和，即绝对值｜a(t)｜＞k*g的时间总和，k常取为0.05；另一类为以相对值定义的相对持时，即最先与最后一个k*a max之间的时段长度，k一般取0.3～0.5。不论实际的强震记录还是人工模拟波形，一般持续时间取结构基本周期的5～10倍。 说明: 有效峰值加速度 EPA＝Sa/2.5 （1） 有效峰值速度 EPV＝Sv/2.5 （2） 特征周期 Tg = 2π*EPV/EPA（3） 1978年美国ATC－3规范中将阻尼比为5％的加速度反应谱取周期为0.1-0.5秒之间的值平均为Sa，将阻尼比为5%的速度反应谱取周期为0.5-2秒之间的值平均为Sv(或取1s附近的平均速度反应谱)，上面公式中常数2.5为0.05组尼比加速度反应谱的平均放大系数。 上述方法使用的是将频段固定的方法来求EPA和EPV，1990年的《中国地震烈度区划图》采用了不固定频段的方法分析各条反应谱确定其相应的平台频段。具体做法是：在对数坐标系中同时做出绝对加速度反应谱和拟速度反应谱，找出加速度反应谱平台段的起始周期T0和结束周期T1，然后在拟速度反应谱上选定平台段，其起始周期为T1(即加速度反应谱平台段的结束周期T1)，结束周期为T2，将加速度反应谱在T0至T1之间的谱值求平均得Sa，拟速度反应谱在T1至T2之间的谱值求平均得Sv（注：生成谱的时候一定要用对数谱），加速度反应谱和拟速度反应谱在平台段的放大系数采用2.5，按公式(1)、（2）、（3）求得EPA、EPV、Tg。 在MIDAS程序中提供将地震波转换为绝对加速度反应谱和拟速度反应谱的功能(工具>地震波数据生成器，生成后保存为SGS文件)，用户可利用保存的SGS文件(文本格式文件)根据上面所述方法计算Sv、Sa、Tg＝Sv/Sa。通过Tg值可判断该地震波是否适合当地场地和地震设计分组，然后将抗震规范中表5.1.2-2中的EPA值与Sa相比求出调整系数（即放大系数），将其代入到地震波调整系数中。将地震波转换为绝对加速度反应谱和拟速度反应谱时注意周期范围要到6秒(建筑抗震规范规定)。 建筑抗震设计规范5.1.2条中规定，采用时程分析方法时，应按照场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。所谓“在统计意义上相符”指的是，其平均影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在各周期点上相差不大于20%。 在MIDAS程序中，可选取两组实际强震记录生成两个SGS文件(调整Sa后的)，然后将一组人
反应谱理论与人工模拟地震波技术简介
第３３卷第２６期?１０６?２００７年９月山西建筑 ＳＨＡＮＸＩＡＲＣＨＩＴＥＣＴＵＲＥ Ｖｄ３３Ｎｏ．２６ Ｓｅｐ．２００７ 文章编号：１００９—６８２５｛２００７）２６—０１０６—０３ 反应谱理论与人工模拟地震波技术简介 邱玉国王玉富 摘要：介绍了反应谱理论的发展历程和国内外研究现状，分析了研究问题的思路，指出了利用反应谱理论来解决实际工程时遇到的问题，并简单介绍了国外对人工模拟地震波技术的应用和研究，为抗震理论提供了参考依据。 关键词：反应谱理论，地震波，随机振动，非弹性地震波 中图分类号：ＴＵ３５２文献标识码：Ａ １概述 反应谱理论是建筑结构抗震设计的重要理论基础之一。从２０世纪５０年代开始，反应谱理论逐渐成为结构抗震设计的重要方法，经过５０多年的发展，目前这种方法已经为世界上大多数国家的设计规范所采用。但是，由于地震产生机理和作用效果的复杂性，采用反应谱理论进行分析和设计与工程实践还存在很多与实际不相符合之处。此外，对于反应地震重要特性的时间问题，反应谱法也无能为力。 人工模拟地震波技术是近年来才发展起来的一项新的结构抗震设计的技术手段，目前主要用于计算机模拟和特别重要结构模型的振动台试验。它能够通过模拟地震波的特性来用于对结构进行时程分析，是～种新兴的、具有革命性意义的试验手段。 图２数值模拟结果２．３计算结果分析 通过数值模拟和试验得到瓦斯管承载力等数值如表２所示。 表２数值模拟和试验结果 Ｉ研究方法承载力仆但ａ最大应变／％最大剪应力／ＳＰａＩ数值模拟７．１４Ｏ．０８４２１６０室内试验６．６２０．０９６４ ３结语 通过对丁集煤矿瓦斯管材质和整体抗外压的试验研究以及数值模拟分析，可以获得如下重要结论： １）通过对管材材质的试验研究表明：工作管材质采用Ｑ３４５，尺寸为柘３０ｒｆｌｌＴｌ×１４ｉｎｌｎ，能够满足强度和稳定性要求。 ２）瓦斯管整体抗外压试验结果表明：工作管抗外压承载力为６，６２ＭＰａ；通过大变形有限元数值计算，采用变形稳定性控制其承载力，结果为７．１４ＭＰａ，两者数值十分接近，说明用文中方法模拟大直径瓦斯管的承载力是可行的。 参考文献： ［１］李正来．瓦斯抽排钻孔定向技术的改进［Ｊ］．安徽科技，２００６（３）：４９—５０． ［２］汪东生．瓦斯抽排技术治理本煤层采空区瓦斯涌出的实践［Ｊ］．煤矿安全，２００６（１）：１３—１５． ［３］张敦伍，任胜杰．瓦斯抽排钻孔防偏斜实践［Ｊ］．矿业安全与环保，２００５（８）：６７—６８． ［４］刘克功，范再良，赵新华．采空区瓦斯抽排法治理综放面瓦斯超限［Ｊ］．煤，１９９８（２）：４８—５０． Ｓｔｕｄｙｉｎｇｏｎｒａｄｉａｌｓｔａｂｉｌｉｔｙｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｌａｒｇｅｐｉｐｅｉｎｍｉｎｅ ＴＯＮＧＷｅｎ－ｌｉｎ Ａｂｓｔｒａｃｔ：ＴｈｅｅｘｐｅｒｉｍｅｎｔａｌａｎｄｖａｌｕｅｓｉｍｕｌａｔｉｏｎｍｅｔｈｏｄｓｈａｖｅｓｔｕｄｉｅｄｔｈｅＤｉｎｇｉｉｃｏａｌｍｉｎｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅｕｎｄｅｒｍｅｃｈａｎｉｃｓｃｈａｒａｃｔｅｒ—ｉｓｔｉｅ．Ｒｅｓｕｌｔｉｎｄｉｃａｔｅｄ：ｔｈｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅｉｓｐｒｅｓｅｎｔｅｄｓｔａｂｉｌｉｔｙｆａｉｌｕｒｅｍｏｄｅｌｉｎｅｎｃｉｒｃｌｅｓｐｒｅｓｓｅｓｓｈａｐｅ，ｉｔｓｓａｆｅｔｙｆａｃｔｏｒｒｅａｃｈｅｓ３．０，ｉｔｉｓｄｅｓｉｇｎｔｈｅｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅａｎｄｔｈｅｃｏｎｓｔｒｕｃｔｐｍｖｉｄｅｓｔｈｅｒｅｆｅｒｅｎｃｅ． Ｋｅｙｗｏｒｄｓ：ｌａｒｇｅｄｉａｍｅｔｅｒｇａｓｔｕｂｅ，ｅｘｐｅｒｉｍｅｎｔａｌｉｎｌａｂ，ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ，ｓｔａｂｉｌｉｔｙｆａｉｌｕｒｅｍｏｄｅｌ 收稿日期：２００７．０４．０６ 作者简介：邱玉国（１９７３。），男，工程师，辽宁工程技术大学软件学院，辽宁阜新１２３０００ 王玉富（１９７０．），男，工程师，中铁十九局集团第三工程有限公司，辽宁辽阳１１１０００
碳酸盐岩储层地震波数值模拟影响因素分析
碳酸盐岩储层地震波数值模拟影响因素分析 通过对比分析已有井的钻测井资料，建立了基于单井的正演模型以及区域连井地质-地球物理模型，并且两者的储层正演响应特征规律性一致。分析讨论了模型建立过程中需考虑的影响因素：子波选择与旁瓣，围岩的尺度、位置、形状。揭示了发育不同厚度时的地震响应特征，进而正确认识了储层的地震相，在实际应用中取得了良好的效果。 标签：地震波数值模拟；有限差分法；碳酸盐岩储层；影响因素 1 概述 在地震勘探中，地震波数值模拟又称地震正演，可供正确认识储层的地震响应特征，为储层预测提供基础。通过分析不同厚度、岩性组合对地震响应的影响，建立储层和地震响应特征之间的联系，為应用地震资料进行储层预测提供一定的依据。地震波数值模拟方法主要分为射线追踪法和波动方程法两类，而其中波动方程法因其能够提供更丰富的波场信息而得到了更加广泛的应用。基于波动方程的数值模拟按照算法不同又分为有限差分法、伪谱法、有限元法及谱元法等，其中有限差分法是最为流行的方法之一[1]。文章采用地震波数值模拟的最常用的波动方程有限差分法正演模拟对下二叠统的储层特征进行了正演影响因素分析实验。 2 基本原理 3 储层正演影响因素分析 在研究区范围内，栖霞组以深灰色厚层状石灰岩为主，含泥质条带及薄层，具灰黑色生物碎屑灰岩、藻灰岩、藻团粒灰岩互层。栖霞组与下伏梁山组黑色含煤岩系及上覆茅口组浅灰色块状灰岩均为整合接触。结合区域地质认识、地震、钻井、测井资料及已有研究成果，建立如图1所示的正演模型。茅口组整体发育大套灰岩，在茅口组底部普遍性发育的一套泥灰岩，由于物性差异较大，对实验结果影响较大。模型仅在透镜体一侧设计了一定厚度的泥灰岩，从实验结果中可以得到效果对比。储层发育在栖霞组上部，储层厚度透镜状变化由中间70米向两侧逐渐减薄，直至储层不发育。在下伏地层中，梁山组黑色含煤系地层虽然很薄（十米左右），但地震波阻抗差异更大，同样不可忽视。 根据上述建立的地质-地球物理正演模型，选用接近实际地震资料的子波进行正演实验。实验选用了30Hz理论Puzirov子波和Riker子波两种不同子波，其中，Puzirov子波波形与Riker子波波形相似均为零相位子波，但旁瓣能量较弱并且能量延续时间较短，具有更高的分辨率。两种不同子波模型正演结果分别如图2所示，图2a是选用30Hz Puzirov子波的结果，图2b则是同一频率常用的Riker 子波的正演结果。总体而言，选用Puzirov子波的正演剖面中，波形信息更加丰富，具有更高的分辨率。在细节刻画方面，图2a中随着储层厚度增大，储层顶
地震波数值模拟方法研究综述.
地震波数值模拟方法研究综述 在地学领域，对于许多地球物理问题，人们已经得到了它应遵循的基本方程(常微分方程或偏微分方程)和相应的定解条件，但能用解析方法求得精确解的只是少数方程性质比较简单，且几何形状相当规则的问题。对于大多数问题，由于方程的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析解。这类问题的解决通常有两种途径。一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。但这种方法只是在有限的情况下是可行的，过多的简化可能导致很大的误差甚至错误的解答。因此人们多年来寻找和发展了另一种求解方法——数值模拟方法。 地震数值模拟(SeismicNumericalModeling)是地震勘探和地震学的基础，同时也是地震反演的基础。所谓地震数值模拟，就是在假定地下介质结构模型和相应的物理参数已知的情况下，模拟研究地震波在地下各种介质中的传播规律，并计算在地面或地下各观测点所观测到的数值地震记录的一种地震模拟方法。地震波场数值模拟是研究复杂地区地震资料采集、处理和解释的有效辅助手段，这种地震数值模拟方法已经在地震勘探和天然地震领域中得到广泛应用。 地震数值模拟的发展非常迅速，现在已经有各种各样的地震数值模拟方法在地震勘探和地震学中得到广泛而有效
的应用。这些地震波场数值模拟方法可以归纳为三大类，即几何射线法、积分方程法和波动方程法。波动方程数值模拟方法实质上是求解地震波动方程，因此模拟的地震波场包含了地震波传播的所有信息，但其计算速度相对于几何射线法要慢。几何射线法也就是射线追踪法，属于几何地震学方法，由于它将地震波波动理论简化为射线理论，主要考虑的是地震波传播的运动学特征，缺少地震波的动力学信息，因此该方法计算速度快。因为波动方程模拟包含了丰富的波动信息，为研究地震波的传播机理和复杂地层的解释提供了更多的佐证，所以波动方程数值模拟方法一直在地震模拟中占有重要地位。 1地震波数值模拟的理论基础 地震波数值模拟是在已知地下介质结构的情况下，研究地震波在地下各种介质中传播规律的一种地震模拟方法，其理论基础就是表征地震波在地下各种介质中传播的地震波传播理论。上述三类地震波数值模拟方法相应的地震波传播理论的数学物理表达方式不尽相同。射线追踪法是建立在以射线理论为基础的波动方程高频近似理论基础上的，其数学表形式为程函方程和传输方程。积分方程法是建立在以惠更斯原理为基础的波叠加原理基础上的，其数学表达形式为波动方程的格林函数域积分方程表达式和边界积分方程表达式。波
模拟地震波测试进行设备抗震性测量
地震是自然界最具破坏力的力量之一，有可能造成毁灭性的生命和财产损失。除了建筑物和结构，地震损坏的设备可能直接或间接地为人们或环境造成危害。为了保障员工和公众的安全，设计安全措施（例如关闭反应堆）或在一定时间内在一定水平的地震作用下的设备应通过抗震鉴定或抗震验证。 为了模拟受地震作用的结构中的设备，地震波的振动测试利用不同类型的运动来有效地模拟假定的地震环境。单频和多频是这这类运动的两大分类。每个类别包括不同情况下的多个波形类型。 单频试验模拟地面在主要频率下的运动振动。它包括以下的运动波形，它们都可以用EDM振动控制系统（VCS）来模拟。 ●连续正弦试验:驻留正弦控制 ●正弦拍频试验：瞬态时程控制 ●衰退正弦测试(Decaying-sine test):瞬态时程控制 ●正弦扫描试验: 扫频正弦控制
地震地面运动是宽频带振动。当不被建筑物或地面过滤时，产生的影响设备的地板运动,趋向于保持宽频带特性。在这样的条件下，多频测试使用复杂的波形，来模拟宽频带地板运动去测试设备。 在EDM软件中，以下波形可用作测试运动，在设备安装时模拟特定的地震波激励(excitation)。一些测试类型将优于其他模拟特定类型的设备应激反应。 杭州锐达数字技术有限公司是美国晶钻仪器公司中国总代理，负责产品销售、技术支持与产品维护，是机械状态监测、振动噪声测试、动态信号分析、动态数据采集、应力应变测试等领域的供应商，提供手持一体化动态信号分析系统、多通道动态数据采集系统、振动控制系统、多轴振动控制系统、三综合试验系统和远程状态监测系统等。更多详情请拨打联系电话或登录杭州锐达数字技术有限公司咨询。
正确选取地震波
地震波的选取方法（MIDAS (2009-05-16 22:51:32) 转载▼ 分类：结构专业 标签： 杂谈 建筑抗震设计规范(GB 50011-2001)的5.1.2条文说明中规定，正确选择输入的地震加速度时程曲线，要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。 频谱特性可用地震影响系数曲线表征，依据所处的场地类别和设计地震分组确定。这句话的含义是选择的实际地震波所处场地的设计分组(震中距离、震级大小)和场地类别(场地条件)应与要分析的结构物所处场地的相同，简单的说两者的特征周期Tg值应接近或相同。特征周期Tg值的计算方法见下面公式(1)、(2)、(3)。 加速度有效峰值按建筑抗震设计规范(GB 50011-2001)中的表5.1.2-2采用。地震波的加速度有效峰值的计算方法见下面公式(1)及下面说明。 持续时间的概念不是指地震波数据中总的时间长度。持时Td的定义可分为两大类，一类是以地震动幅值的绝对值来定义的绝对持时，即指地震地面加速度值大于某值的时间总和，即绝对值｜a(t)｜＞k*g的时间总和，k常取为0.05；另一类为以相对值定义的相对持时，即最先与最后一个k*amax之间的时段长度，k 一般取0.3～0.5。不论实际的强震记录还是人工模拟波形，一般持续时间取结构基本周期的5～10倍。 说明: 有效峰值加速度EPA＝Sa/2.5 （1） 有效峰值速度EPV＝Sv/2.5 （2） 特征周期Tg = 2π*EPV/EPA （3） 1978年美国ATC－3规范中将阻尼比为5％的加速度反应谱取周期为0.1-0.5秒之间的值平均为Sa，将阻尼比为5%的速度反应谱取周期为0.5-2秒之间的值平均为Sv(或取1s附近的平均速度反应谱)，上面公式中常数2.5为0.05组尼比加速度反应谱的平均放大系数。 上述方法使用的是将频段固定的方法来求EPA和EPV，1990年的《中国地震烈度区划图》采用了不固定频段的方法分析各条反应谱确定其相应的平台频段。具体做法是：在对数坐标系中同时做出绝对加速度反应谱和拟速度反应谱，找出加速度反应谱平台段的起始周期T0和结束周期T1，然后在拟速度反应谱上选定平台段，其起始周期为T1(即加速度反应谱平台段的结束周期T1)，结束周期为T2，将加速度反应谱在T0至T1之间的谱值求平均得Sa，拟速度反应谱在T1至T2之间的谱值求平均得Sv，加速度反应谱和拟速度反应谱在平台段的放大系数采用2.5，按公式(1)、（2）、（3）求得EPA、EPV、Tg。 在MIDAS程序中提供将地震波转换为绝对加速度反应谱和拟速度反应谱的功
midas地震波选取
地震波的选取方法 建筑抗震设计规范(GB 50011-2001)的5.1.2条文说明中规定，正确选择输入的地震加速度时程曲线，要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。 频谱特性可用地震影响系数曲线表征，依据所处的场地类别和设计地震分组确定。这句话的含义是选择的实际地震波所处场地的设计分组(震中距离、震级大小)和场地类别(场地条件)应与要分析的结构物所处场地的相同，简单的说两者的特征周期Tg值应接近或相同。特征周期Tg值的计算方法见下面公式(1)、(2)、(3)。 加速度有效峰值按建筑抗震设计规范(GB 50011-2001)中的表5.1.2-2采用。地震波的加速度有效峰值的计算方法见下面公式(1)及下面说明。 持续时间的概念不是指地震波数据中总的时间长度。持时T d 的定义可分为两大类，一类是以地震动幅值的绝对值来定义的绝对持时，即指地震地面加速度值大于某值的时间总和，即绝对值｜a(t)｜＞k*g的时间总和，k常取为0.05；另一类为以相对值定义的相对持时，即最先与最后一个k*a max 之间的时段长度，k一般取0.3～0.5。不论实际的强震记录还是人工模拟波形，一般 说明: 特征周期 Tg = 2π 用了不固定频段的方法分析各条反应谱确定其相应的平台频段。具体做法是：在对数坐标系中同时做出绝对加速度反应谱和拟速度反应谱，找出加速度反应谱平台段的起始周期T0和结束周期T1，然后在拟速度反应谱上选定平台段，其起始周期为T1(即加速度反应谱平台段的结束周期T1)，结束周期为T2，将加速度反应谱在T0至T1之间的谱值求平均得Sa，拟速度反应谱在T1至T2之间的谱值求平均得Sv，加速度反应谱和拟速度反应谱在平台段的放大系数采用2.5，按公式 (1)、（2）、（3）求得EPA、EPV、Tg。 在MIDAS程序中提供将地震波转换为绝对加速度反应谱和拟速度反应谱的功能(工具>地震波数据生成器，生成后保存为SGS文件)，用户可利用保存的SGS文件(文本格式文件)根据上面所述方法计算Sv、Sa、Tg。通过Tg值可判断该地震波是否适合当地场地和地震设计分组，然后将抗震规范中表5.1.2-2中的EPA值与Sa相比求出调整系数，将其代入到地震波调整系数中。将地震波转换为绝对加速度反应谱和拟速度反应谱时注意周期范围要到6秒(建筑抗震规范规定)。 建筑抗震设计规范5.1.2条中规定，采用时程分析方法时，应按照场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。所谓“在统计意义上相符”指的是，其平均影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比，在各周期点上相差不大于20%。 在MIDAS程序中，可选取两组实际强震记录生成两个SGS文件(调整Sa后的)，然后将一组人工模拟的加速度时程曲线也保存为SGS文件，将三个SGS文件的数值取平均后与振型分解反应谱
地震波应用
建筑抗震设计规范(GB 50011-2001)的5.1.2条文说明中规定，正确选择输入的地震加速度时程曲线，要满足地震动三要素的要求，即频谱特性、有效峰值和持续时间要符合规定。 频谱特性可用地震影响系数曲线表征，依据所处的场地类别和设计地震分组确定。这句话的含义是选择的实际地震波所处场地的设计分组(震中距离、震级大小)和场地类别(场地条件)应与要分析的结构物所处场地的相同，简单的说两者的特征周期Tg值应接近或相同。特征周期Tg值的计算方法见下面公式(1)、(2)、(3)。 加速度有效峰值按建筑抗震设计规范(GB 50011-2001)中的表5.1.2-2采用。地震波的加速度有效峰值的计算方法见下面公式(1)及下面说明。 持续时间的概念不是指地震波数据中总的时间长度。持时Td的定义可分为两大类，一类是以地震动幅值的绝对值来定义的绝对持时，即指地震地面加速度值大于某值的时间总和，即绝对值｜a(t)｜＞k*g的时间总和，k常取为0.05；另一类为以相对值定义的相对持时，即最先与最后一个k*amax之间的时段长度，k一般取0.3～0.5。不论实际的强震记录还是人工模拟波形，一般持续时间取结构基本周期的5～10倍。 说明: 有效峰值加速度 EPA＝Sa/2.5 （1） 有效峰值速度 EPV＝Sv/2.5 （2） 特征周期Tg = 2π*EPV/EPA （3） 1978年美国ATC－3规范中将阻尼比为5％的加速度反应谱取周期为0.1-0.5秒之间的值平均为Sa，将阻尼比为5%的速度反应谱取周期为0.5-2秒之间的值平均为Sv(或取1s附近的平均速度反应谱)，上面公式中常数2.5为0.05组尼比加速度反应谱的平均放大系数。 —1 —
地震波波动方程数值模拟方法
地震波波动方程数值模拟方法 地震波波动方程数值模拟方法主要包括克希霍夫积分法、傅里叶变换法、有限元法和有限差分法等。 克希霍夫积分法引入射线追踪过程，本质上是波动方程积分解的一个数值计算，在某种程度上相当于绕射叠加。该方法计算速度较快，但由于射线追踪中存在着诸如焦散、多重路径等问题，故其一般只能适合于较简单的模型，难以模拟复杂地层的波场信息。 傅里叶变换法是利用空间的全部信息对波场函数进行三角函数插值，能更加精确地模拟地震波的传播规律，同时，利用快速傅里叶变换(FFT)进行计算，还可以提高运算效率，其主要优点是精度高，占用内存小，但缺点是计算速度较慢，对模型的适用性差，尤其是不适应于速度横向变化剧烈的模型. 波动方程有限元法的做法是：将变分法用于单元分析，得到单元矩阵，然后将单元矩阵总体求和得到总体矩阵，最后求解总体矩阵得到波动方程的数值解；其主要优点是理论上可适宜于任意地质体形态的模型，保证复杂地层形态模拟的逼真性，达到很高的计算精度，但有限元法的主要问题是占用内存和运算量均较大，不适用于大规模模拟，因此该方法在地震波勘探中尚未得到广泛地应用。。 相对于上述几种方法，有限差分法是一种更为快速有效的方法。虽然其精度比不上有限元法，但因其具有计算速度快，占用内存较小的优点，在地震学界受到广泛的重视与应用。 声波方程的有限差分法数值模拟 对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波方程描述： )()(2222222t S z u x u v t u +??+??=?? （4-1） (,)v x z 是介质在点（x , z ）处的纵波速度，u 为描述速度位或者压力的波场，)(t s 为震源函数。 为求式（4-1）的数值解，必须将此式离散化，即用有限差分来逼近导数，用差商代替微商。为此，先把空间模型网格化（如图4-1所示）。 设x 、z 方向的网格间隔长度为h ?，t ?为时间采样步长，则有： z ?,i j 1,i j +2,i j +1,i j -