论地震波速度在地震勘探中的应用

论地震波速度在地震勘探中的应用
论地震波速度在地震勘探中的应用

论地震波速度在地震勘探中的应用

论文提要

地震勘探是地球物理勘探中的一种重要的方法,始于19世纪中叶,近五十年来经历了光点记录、模拟磁带记录、数字磁带记录三个阶段,广泛的应用于石油和天然气资源勘查、煤田勘查、工程地质勘查、及某些金属矿的勘查等方面。地震勘探利用地震波在不同的岩石中的传播速度不同,研究地下的地质构造,判断油气藏等的可能储藏位置。在这个过程中,地震波速度是一个重要的参数,根据不同的实际情况,可将地震波速度分为均一速度、层速度、等效速度、叠加速度、均方根速度等具有不同意义的速度,可谓是贯穿整个地震资料处理的过程。正确的应用各种速度可以使所得资料更准确,更接近地下实际形态。

正文

一、地震勘探发展史

地震勘探是地球物理勘探中的一种最重要的方法。它的原理是由人工制造的强烈的震动(一般是在地下不深处的爆炸)所引起的弹性波在岩石中的传播时,常遇到岩层的分界面,便产生反射波或折射波,在它返回地面时用高度灵敏的仪器记录下来,根据波的传播路线和时间,确定发生反射波或折射波的岩层介面的埋藏深度和形状,认识地下地质构造,以寻找油气圈闭。

地震勘探始于19世纪中叶。1845年,R.马利特曾用人工激发的地震波来测量弹性波在地壳中的传播速度。这可以说是地震勘探方法的萌芽。在第一次世界大战期间,交战双方都曾利用重炮后坐力产生的地震波来确定对方的炮位。

反射法地震勘探最早起源于1913年前后R.费森登的工作,但当时的技术尚未达到能够实际应用的水平。1921年,J.C.卡彻将反射法勘探投入实际应用,在美国俄克拉何马州首次记录到人工地震产生的清晰的反射波。1930年,通过反射法地震勘探工作,在该地区发现了3个油田。从此,反射法进入了工业应用的阶段。

50-60年代,反射法的光点照相记录方式被模拟磁带记录方式所代替,从而可选用不同因素进行多次回放,提高了记录质量。70年代,模拟磁带记录又为数字磁带记录所取代,形成了以高速数字计算机为基础的数字记录、多次覆盖技术、地震数据处理技术相互结合的完整技术系统,大大提高了记录精度和解决地质的能力。

从70年代初期开始,采用地震勘探方法研究岩性和岩石孔隙所含流体成分。根据地震时间剖面振幅异常来判定气藏的“亮点”分析,以及根据地震反射波振幅与炮检距关系来预测油气藏的AVO(振幅随偏移距的变化,该技术利用CDP道集资料,分析反射波振幅随炮检距变化规律,估算界面的弹性参数柏松比,进一步推断地层的岩性和含油气情况)分析,已有很多成功的例子。从地震反射波推算地层波阻抗和层速度的地震拟测井技术,在条件有利时,可以取得有地质解释意义的实际效果。现代的地震勘探正由以构造勘探为主的阶段向着岩性勘探的方向发展。

中国于1951年开始进行地震勘探,并将其应用于石油和天然气资源勘查、煤田勘查、工程地质勘查、及某些金属矿的勘查。

二、与地震波速度有关的因素

地震发生时,震源区的介质发生急速的破裂和运动,这种扰动构成了一个波源。由于地球介质的连续性,这种波动就向地球内部及表层各处传播开去,形成了连续介质中的弹性波,所以地震波是指从震源产生向四外辐射的弹性波。地球内部存在着地震波速度突变的基干界面、莫霍面和古登堡面,将地球内部分为地壳、地幔和地核三个圈层(表1-1)。

表1-1 地球内部层圈波速表

地震勘探是以研究地震波在岩层中的传播规律为基础的,岩石的弹性性质(主要表现为地震波的传播速度)不同,地震波在其中的传播情况也就不同,地震勘探正式利用这种关系来研究地下地层的地质构造的。应用中的基本公式是H=1/2Vt,H是界面的深度,V是地震波的传播速度,t是地震波从地面垂直向下到界面在返回地面的旅行时,从这一基本关系式可以看到速度参数V的重要性。

而地震波在岩层中的传播速度和岩层的性质,如弹性常数或岩石的成分、密度、埋藏深度、地质年代,孔隙率等因素有关。一般说来,火成岩的地震波速度的变化范围比变质岩和沉积岩小,火成岩地震波速度的平均值比其他类型岩石的要高。大多数变质岩的地震波速度变化范围比较大。地震波速度在几种常见岩石中的变化范围可见表1-2。

表1-2 各类岩石的波速

从大量的资料(表1-3)可得出结论:地震波速度与岩石密度成正比,有时可把速度与密度表示成一种近似的线性关系,如某些石灰岩和砂页岩的关系可近似表示为:

V=6ρ-11

而地震波速度在沉积岩中的分布存在垂直梯度、水平梯度及成层分布等特点。

表1-3 各种岩石的波速和波阻抗

三、常用的几种地震波速度

地震波在地层中的传播速度是一个十分重要的参数,但又很难精确测定它的数值。根据生产工作的需要和地震勘探方法技术所能达到的水平,对极其复杂的实际情况作种种简化,建立各种简化介质模型,并引进了各种速度的概念。而且每种速度概念都有它的意义、引入的原因、计算或测定的方法以及适用范围等,不能混为一谈。

当地下介质均匀,界面水平时,地震波的传播速度不随深度变化,不随传播路径改变而发生变化,在界面上所有的各点处速度一样,故称做均一速度;界面倾斜时,根据实际速度与倾角的关系确定的替代速度称作等效速度。

水平层状介质情况下,地震波在每一层中传播速度不同,但同一层中速度不变,地

震波在某一层中的速度称为层速度。它直接反映地层的岩性,能用来划分地层。由于此时地震波的传播时间和距离之间并不是简单的二次双曲线的关系,所以当炮检距较小或为零时,可将地下介质简单化,用地震波在地层中垂直传播的总时间除地层总厚度求出的平均速度代替层速度。若把反射波实际曲线近似地当作双曲线,则求出的波速称作这一水平层状介质的均方根速度。

而在对地震记录进行共反射点叠加时,能得到最佳叠加效果的速度称作叠加速度。叠加速度是射线参数P(t pi=z i/[v i(1-p2v i2)1/2] , 其中z是深度,p是射线参数,v是速度,t是时间)的函数,也是炮检距的函数,且随炮检距的增大而增大。

通常地震波速度的测量方法包括声波测井的直接测量法和地震勘探数据的间接测量法,根据这两方面的信息,可以得到以上这些不同物理意义和用途的地震速度。四、地震波速度在地震勘探中的应用

在资料处理和解释的过程中,速度在许多环节都是一个重要的参数。例如,在进行动校正时,要有叠加速度资料;进行偏移叠加时,要有偏移叠加速度。时深转换时,要有平均速度资料。通过速度谱分析,获得叠加速度,进而求取均方根速度、层速度。

由于在原始的野外地震记录中非零炮检距道上的波至时间与零炮检距道上的时间存在差值,所以在地震资料处理中需要校正时差,这个消除接收点偏离炮点而形成的时差的过程称为动校正。只要在动校正方程中采用了正确的速度,就能校正不同炮检距的双曲线的影响(如表1-4)。

表1-4 各种地层模型的NMO速度

在速度分析中,对一个CMP道集用一系列速度重复对道集作动校正,每校正一次得到一张图像,当所用速度高于介质速度,那么双曲线不能完全拉平,称之为欠校正;反之,所用速度低于介质速度,那么就得到过校正的结果;而使道集的反射波曲线拉得最平的速度就是叠前最佳动校正速度(如图1-1)。

图1-1 动校正速度的选择对动校正效果的影响

以上这种速度分析的方法称为速度扫描法,该方法精度高,但是工作量大;而对于速度谱不好或者地下地形复杂的地区可采用常速叠加法。

对于多个相邻的CMP 道集用一系列常速度进行校正叠加,叠加后同相轴振幅强且连续性好的所对应的速度即所求速度,这个过程称为常速叠加法。在选择该方法所用的速度时,除了要预计地下实际速度可能的范围以外,还应该考虑到叠加数据所需要的速度范围和试验这些叠加速度采用的间隔大小。在范围选择上,要考虑倾斜同相轴和有用的非平面反射,它们可能具有非常高的叠加速度。在选择等速度间隔时,应考虑到是按照不同炮检距上的动校正时差,不是按照速度来做速度估计,这样可防止对高速同相轴形成过密采样和对低速同相轴欠采样。

与常速叠加法相比,速度谱是根据CMP 道集中地震道的互相关而不是叠加同相轴的横向连续性,因此比常速叠加法更适应于处理具有多次反射问题的数据,但不太适合于解决复杂构造问题。

做速度谱的原理是:对共反射点时距曲线t 2

=t 02

+x 2

/v α2

固定t 0选取v αj 计算t i ,按t i 取各道记录振幅相加,当选定速度等于叠加速度 v α时,计算时间与实际记录时间的同相轴一致,由于同相叠加故振幅增大,如果选定速度不等于叠加速度时,计算时间与实际记录时间同相轴一致,不能同相叠加,振幅减小,一系列速度中,叠加振幅最大的即所求的叠加速度。

由于地幔中的某一个圈层的岩石强度减弱,导致该层地震波速度较之地幔的最外层减弱6%,一般为每秒数百米,该层称为低速层。若在低速带中激发,能量将被大量消耗,频率降低,低速带底面又是一个强反射面,可以形成多次反射。而且低速带使地震波到达时间观测点的时间延迟,当其厚度或速度变化较大时,观测时间的精度将大受影响,甚至使波形歪曲。同时,地表的不规则变化也会影响反射时间。这时,必须进行低速带的时间校正。这种消除地表高程影响及浅层低降速带影响的过程称为静校正。

在地震资料处理中,原始地震记录中由地下界面的构造复杂性引起的时差,绝大部分可通过野外静校正予以消除,而往往它们不包括地形突变、风化层基地和风化层速度变化的影响,这些影响必须通过剩余静校正加以消除。经野外静校正的CMP道集利用速度分析(作速度谱)中所提取的一组初始速度进行动校正,由于此时的动校正速度不准确以及存在地表、低速区域的变化影响,故应做剩余静校正。经过剩余静校正,旅行时间有偏差的CMP道集的反射排列得好多了,而那些不需要这种校正的道保持不变,同时,增强和未增强的叠加剖面上的反射连续性,都有了改进,并明显消除了假构造。紧跟着剩余静校正要重做速度分析,以更新速度提取,利用更新速度再次做动校正。最终动校正后的CMP道集才能进行叠加偏移,得到的各反射点的同相轴的连线即反映了地下界面的基本形态。

由上述处理过程可知,地震波速度在勘探中占有重要的地位,可以说贯穿整个地震处理的过程。

同时地震波速度又是最直接用于岩性解释的地震参数之一,如应用纵横波速度比划分岩性和检测油气、利用速度资料估算砂泥岩百分比等。处理后的地震信息经过地震解释可转变成地质成果,如做出构造解释、地层解释,岩性和烃类检测解释及综合解释,绘出有关的成果图件,对测区作出含油气评价,提出钻井位置等。

五、地球物理勘探的发展方向

地球物理勘探需要引进现代电子计算器技术,进一步压制干扰,提高分辨能力,提取更多的有用信息,发展反演的理论和技术,提高各类地质问题的地球物理解释、推断效果并不断提高地球物理数据处理的工作效率和图像处理技术。地球物理勘探仪器要向轻便化、高精度、多功能、数字化、系列化和智能化的方向发展。现代地质学理论的发展,使深部地质问题的研究愈显重要。应用于这方面研究的人工地震反射剖面、大地电磁测深、重力、磁法、地热等地球物理勘探方法,已显示出其潜力和优越性。

参考文献

姚长纲主编.石油地质学基础[M].石油物探学校教学用书

陆基梦主编.地震勘探原理[M].石油物探大学出版,1993.

李振春、张军华编.地震数据处理方法[M].中国石油大学出版社,2004.

顾功叙编.地球物理勘探基础[M].地质出版社.1990.

吴功建、林清湲﹑高锐编.地球物理方法在地质和找矿中的应用.地质出版社.1988.

P.V.莎玛著.王恕铭等译.地质学研究中的地球物理方法.地质出版社.1983.

在 ansys 中如何 施加 地震波

三向输入简化后的单向输入 首先,将三个方向的地震加速度放到一个文本文件里,如accexyz.txt,在这个数据文件里共放三列数据,每列为一个方向的地震加速度值,这里仅给出数据文件中前几行的数据: -0.227109E-02 -0.209046E+00 0.467072E+01 -0.413893E-02 -0.168195E+00 0.261523E+01 -0.574753E-02 -0.157890E+00 0.809014E-01 -0.731227E-02 -0.152996E+00 0.119975E+01 -0.876865E-02 -0.138102E+00 0.130902E+01 -0.101067E-01 -0.131582E+00 0.143611E+00 ....................... 然后,再建一个文本文件用来存放三个方向的地震加速度时间点,如time.txt,在这个数据文件里仅一列数据,对应于加速度数据文件里每一行的时间点,这里给出数据文件中前几行数据: 0.100000E-01 0.200000E-01 0.300000E-01 0.400000E-01 0.500000E-01 0.600000E-01 ....................... 编写如下的命令流文件,并命名为acce.inp *dim,ACCEXYZ,TABLE,2000,3 !01行 *vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行(3e16.6) !03行 *vread,ACCEXYZ(1,0),time,txt !04行 (e16.6) !05行 ACCEXYZ(0,1)=1 !06行 ACCEXYZ(0,2)=2 !07行,同上 ACCEXYZ(0,3)=3 !08行,同上 finish /SOLU ANTYPE,trans btime=0.01 !定义计算起始时间 etime=15.00 !定义计算结束时间 dtime=0.01 !定义计算时间步长 *DO,itime,btime,etime,dtime time,itime AUTOTS,0 NSUBST,1, , ,1 KBC,1 acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度 SOLVE

论地震勘探中几种主要地震波

论地震勘探中的几种主要地震波 论文提要 地震勘探,就是通过人工方法激发地震波,研究地震波在地层中传播的情况,以查明地下地质构造,为寻找油气田或其它勘探目的服务的一种方法。也可以理解为就是利用地震子波从地下地层界面反射回地面时带回来的旅行时间和形状变化的信息,用以推断地下的底层构造和岩性。地震勘探在勘探已有的各种物探方法中,是最有效地方法。在地震勘探中用炸药激发时,一声炮响之后会产生各种各样的地震波。按波在传播过程中质点震动的方向来区分,可以纵波和横波;根据波动所能传播的空间范围而言,地震波又可以分为体波和面波;按照波在传播过程中的传播路径的特点,又可以把地震波分为直达波、反射波、透射波、折射波,等等。地震勘探在石油勘探中除了能产生来自地层界面有用的反射波外,还会产生各种各样的干扰波。因此,我们要更好的了解各种波的产生、特点、用途,等等。下面简单介绍几种地震勘探中产生的地震波。 正文 一、反射波 (一)反射波的形成 1、几何地震学的观点 当炸药在井中爆炸激发地震波时,在雷管引爆几百微妙之内爆炸便完成了,在接近爆炸点的压强是一个延续时间很短的尖脉冲,爆炸脉冲向外传播,压强逐渐减少,地层开始产生弹性形变,形成地震波。地震波继续传播,由于介质对高频的吸收,地震波信号减小。当波入射到两种介质的分界面时(当上层介质波阻抗与下层介质波阻抗不等时,弹性地震波才会发生反射;上层介质波阻抗与下层介质波阻抗差别越大,反射波越强——反射波条件),一部分波回到第一种介质中,这就是所谓的反射波。如图所示 2、物理地震学观点 地震波从震源出发以球面波的方式向下传播,到达反射界面S,S可以就看成有许多

时程分析中地震波输入位置的讨论

时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系,能真实地反应结构地震相应随时间变化的全过程,是抗震分析的一种重要方法[1]。目前有限元软件可以实现结构的时程分析,但是在不同的软件中,其实现方式不同,主要区别在地震波的输入位置不同。本文通过有限元软件ABAQUS采用不同的地震波输入位置对同一结构进行时程分析分析,对比结构相同位置的时程位移曲线,结果表明结构在采用不同地震波输入位置的时程分析中,结构的地震响应基本一致。 关键词:时程分析、有限元软件、钢筋混凝土剪力墙 Abstract: The time history analysis method to analyze the available structure through direct power to the relationship between the corresponding changes over time, truly reflect the structure of earthquake corresponding to the whole process of change over time, is an important method of seismic analysis [1]. Finite element software can be time-history analysis of the structure, but in different software in different ways, the main difference between the different positions in the seismic wave input. In this paper the finite element software ABAQUS using different seismic wave input location on the same structure, process analysis analysis, contrast structure the same location of when the process displacement curve, the results show that the structure using different seismic waves enter the position time history analysis, the seismic response basically the same. Keywords: time history analysis, finite element software, reinforced concrete shear walls 一、引言 在时程分析等动力学问题中,地震力以加速度形式从基础固定处输入。由于结构的刚度不是无限大,在结构上的加速度反应与基础输入的加速度并不相同。在很多时候,结构的加速度比基础输入的加速度更大,即对输入的加速度有一个动力放大效应。在单自由度弹性体系中,体系最大绝对加速度与地面运动最大加速度的比值,即称为动力系数[2] (1) 动力系数与结构的动力学特性和输入的地震波的频率特性有关。它与地震系数k的乘积即为单自由度体系的地震影响系数。 因此,从原理上讲,时程分析是将地震波的加速度时程曲线作用到结构的基础约束处,得到上部结构的各种地震反应。但是在不同的软件中,其实现方

地震波使用说明

地震波使用说明 此目录下提供了四类场地土的地震波时程曲线和上海人工波。 按照场地土类型(1,2,3或4),选择时程曲线。在定义时程工况时,对于多遇或罕遇地震,按比例调整时程曲线的最大值。中国抗震规范规定,作为抗震计算中底部剪力法和振型分解反应谱法的补充方法,对于特别不规则,特别重要的和较高的结构应采用时程分析法进行多遇地震下的补充计算。 可取多条时程曲线的计算结果的平均值与振型分解反应谱法计算结果的较大值。 采用时程分析法时,应咱建筑场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线,其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。 其加速度时程最大值可按规范中对于多遇和罕遇地震在不同烈度下的值。 弹性时程分析时,每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计算结果的65%,多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80% 。 可使用弹塑性时程分析法计算罕遇地震下结构的变形。 时程分析是一个承受随时间变化的指定荷载结构的逐步动态反应分析,可以是线性或非线性的。 此章对时程分析进行一般的描述,特别是线性时程分析。 定义时程函数 用户可使用“从文件中添加函数”,导入已定义的文本文件,即实测的时程曲线;也可使用程序内置的时程函数。

时程函数定义对话框 时程函数定义对话框中的条目解释如下: ?函数名 通过在编辑框中直接键入以指定或修改时程函数的名称。 ?函数文件 1.在函数文件域点击浏览按钮以调出一个对话框,在此可找出包含时程函数的 文本文件名。注意文件名显示在文件名框中 2.在 "要跳过的标题行" 编辑框中输入一个希望ETABS在文本文件中跳过的 行数。 3.在 "每行要跳过的前缀字符" 编辑框中输入一个希望ETABS在文本文件中 每行要跳过的字符数。 4.在 "每行的点数" 编辑框中输入一个数告诉ETABS文本文件每行的绘图点 数。

地震波运动学理论

第二章地震波运动学理论 一、名词解释 1. 地震波运动学:研究在地震波传播过程中的地震波波前的空间位置与其传播时间的关系,即研究波的传播规律,以及这种时空关系与地下地质构造的关系。 2. 地震波动力学:研究地震波在传播过程中波形、振幅、频率、相位等特征的及其变化规律,以及这些变化规律与地下的地层结构,岩石性质及流体性质之间存在的联系。 3. 地震波:是一种在岩层中传播的,频率较低(与天然地震的频率相近)的波,弹性波在 岩层中传播的一种通俗说法。地震波由一个震源激发。 4. 地震子波:爆炸产生的是一个延续时间很短的尖脉冲,这一尖脉冲造成破坏圈、塑性带,最后使离震源较远的介质产生弹性形变,形成地震波,地震波向外传播一定距离后,波形逐渐稳定,成为一个具有2-3个相位(极值)、延续时间60-100毫秒的地震波,称为地震子波。地震子波看作组成一道地震记录的基本元素。 5.波前:振动刚开始与静止时的分界面,即刚要开始振动的那一时刻。 6.射线:是用来描述波的传播路线的一种表示。在一定条件下,认为波及其能量是沿着一条“路径”从波源传到所观测的一点P。这是一条假想的路径,也叫波线。射线总是与波阵面垂直,波动经过每一点都可以设想有这么一条波线。 7. 振动图和波剖面:某点振动随时间的变化的曲线称为振动曲线,也称振动图。地震勘探中,沿测线画出的波形曲线,也称波剖面。 8. 折射波:当入射波大于临界角时,出现滑行波和全反射。在分界面上的滑行波有另一种特性,即会影响第一界面,并激发新的波。在地震勘探中,由滑行波引起的波叫折射波,也叫做首波。入射波以临界角或大于临界角入射高速介质所产生的波 9.滑行波:由透射定律可知,如果V2>V1 ,即sinθ2 > sinθ1 ,θ2 > θ1。当θ1还没到90o时,θ2 到达90o,此时透射波在第二种介质中沿界面滑行,产生的波为滑行波。 10.同相轴和等相位面:同向轴是一组地震道上整齐排列的相位,表示一个新的地震波的到达,由地震记录上系统的相位或振幅变化表示。 11.地震视速度:当波的传播方向与观测方向不一致(夹角θ)时,观测到的速度并不是波前的真速度V,而是视速度Va。即波沿测线方向传播速度。 12 波阻抗:指的是介质(地层)的密度和波的速度的乘积(Zi=ρiVi,i为地层),在声学中称为声阻抗,在地震学中称波阻抗。波的反射和透射与分界面两边介质的波阻抗有关。只有在Z1≠Z2的条件下,地震波才会发生反射,差别越大,反射也越强。 13.纵波:质点振动方向与波的传播方向一致,传播速度最快。又称压缩波、膨胀波、纵波或P-波。 14.横波:质点振动方向与波的传播方向垂直,速度比纵波慢,也称剪切波、旋转波、横波或S-波,速度小于纵波约0.7倍。横波分为SV和SH波两种形式。 15.体波:波在无穷大均匀介质(固体)中传播时有两种类型的波(纵波和横波),它们在介质的整个立体空间中传播,合称体波。 16共炮点反射道集:在同一炮点激发,不同接收点上接收的反射波记录,称为共炮点道集。在野外的数据采集原始记录中,常以这种记录形式。可分单边放炮和中间放炮。 17.面波:波在自由表面或岩体分界面上传播的一种类型的波。 18.纵测线和非纵测线:激发点与接收点在同一条直线上,这样的测线称为纵测线。用纵测线进行观测得到的时距曲线称为纵时距曲线。激发点不在测线上,用非纵测线进行观测得到的时距曲线称为非纵时距曲线。

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 ………………

ANSYS地震波的输入

对于地震波的输入,可以把荷载记录作成文件,利用apdl的读取功能读入数据库中。下面的例子是自己编的一个小文件。修改一下可以更简洁。 Fini /config,nres,1000 *dim,aceX,TABLE,3000,1 *dim,aceY,TABLE,3000,1 *dim,aceZ,TABLE,3000,1 *creat,ff *vread,aceX(1,1),acex,txt,,1 (e16.6) *vread,aceX(1,0),acexTT,txt,,,1 (e16.6) ACEX(0,1)=1 *end /input,ff *creat,ff *vread,aceY(1,1),txt,,1 (e16.6) *vread,aceY(1,1),ACETT,,,1 (e17.6) ACEY(0,1)=1 *end /input,ff *creat,ff *vread,aceZ(1,1),txt,,1 (e16.6) *vread,aceZ(1,0),ACETT,,,1 (e17.6) ACEZ(0,1)=1 *end /input,ff !地震波时程记录分成了3个文件,每个文件是一列。分别记录x,y,z方向的加速度。Accett是时间记录。 这样就可以把加速度记录读取到ansys数据库中作为数组。 也可以把加速度记录作成一个文件,这样程序就简单多了。 下面是计算部分语句: /SOLU ANTYPE,trans !求解其自己选了 TM_START=0.01 TM_END=15.00 TM_INCR=0.01 *DO,TM,TM_START,TM_END,TM_INCR TIME,tm

第六章练习题

第六章练习题 一、名词解释 1.平均速度 2.叠加速度 3.均方根速度 4.等效速度 5.层速度 二、 填空题 1. 地震波在石油中传播速度为________m/s至________m/s;在石灰岩中传播速度为_________m/s至___________m/s. 2. 地震波的速度与孔隙度成__________;同种性质的岩石,孔隙度越大地震波速度越____________反之则越__________. 3.描述地震波速度与岩石孔隙度经验公式是_________平均方程.公式为1/V=(1-Ф)/Vm+Ф/Vl.式中V是__________Vl是孔隙中__________Ф是岩石_________. 4.地震波在岩石中传播速度与岩石的孔隙度成______比例;与岩石的密度成_____________. 5. 岩石孔隙中充满水的时的速度_______充满油时的速度,充满油时的速度 ________充满气时的速度. 6.地震波速度,一般随地层深度的_______而增大, 随地层压力的增大而_______. 7. 岩石年代越老, 其速度越___________,反之则_________. 8. 在速度谱上拾取的速度是___________在时一深转换尺上读取的速度是_______________. 9. 分析叠加速度谱拾取________速度, 主要的是便于________和水平叠加. 10. 用VSP测井能得到的速度资料包括____________和______________资料. 11. 一般进行时深转换采用的速度为_________________.研究地层物性参数变化需采用__________________. 12.用于计算动校正量的速度称为______________速度,它经过倾角校正后即得到________________.

利用PEER网站选取地震波的方法

利用PEER网站选取地震波的方法 云南省土木建筑学会建筑结构专业委员会 PEER是Pacific Earthquake Engineering Research Center(太平洋地震工程研究中心)的简称,设立在美国的加州大学伯克利分校(University of California, Berkeley)。由其运作的PEER Ground Motion Database(PEER地震动数据库)提供了大量的世界各地的地震记录,可自由下载。该网站提供了丰富的查询手段,可按距离、场地、震源类型等条件选择地震记录,也提供了按目标反应谱选择的手段。 该网站是https://www.360docs.net/doc/d114660383.html,/peer_ground_motion_database,首页如下:

如果要按目标反应谱(例如我国GB50011-2010的地震影响系数曲线)选取地震波,需要事先准备好目标反应谱的数据文件。云南省土木建筑学会建筑结构专业委员会和昆明理工大学提供了Microsoft Excel 2003格式的文件Spectrum-2010.xls,可用于生成我国建筑、公路、水工等抗震设计规范规定的反应谱对应的数据文件。用法简述如下: 一、生成中国规范的目标谱 1、打开该文件后,在“图形”工作表(左下角选择)上进行最大地震影响系数(αmax)、地震分组、场地类别的选择(均为下拉菜单)。 2、在“表格”工作表(左下角选择)上即可得到所需要的各种数据。例如,PEER地震动数据库的目标谱格式为“周期(s)-谱加速度(g)”,复制拷贝其中的A列和D列即可。

3、利用具有“列编辑模式”功能的文本编辑软件(例如UltraEdit 等。也可直接使用Excel,注意粘贴时采用“选择性粘贴”-“数值”), 得到如下形式的文本: 4、将其保存为“.csv”后缀的文本文件。 二、利用PEER地震动数据库获得地震记录 1、上述准备工作完成后,进入PEER地震动数据库首页。点击 “Scaled”。

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.强度调整。将地震记录的加速度值按适当的比例放大或缩小,使其峰值加速度等于事先所确定的设计地震加速度峰值。即令 其中a(为记录的加速度值为调整后的加速度值;A众为设计地震加速度峰值;。为记录的加速度峰值。这种调整只是针对原记录的强度进行的,基本上保留了实际地震记录的特征。也就是所说的(强度修正。将地震波的加速度峰值及所有的离散点都按比例放大或缩小以满足场地的烈度要求)

2.频率调整考虑到场地条件对地震地面运动的影响,原则上所选择的实际地震记录的富氏谱或功率谱的卓越周期乃至形状,应尽量与场地土相应的谱的特性一致。如果不一致,可以调整实际地震记录的时间步长,即将记录的时间轴“拉长”或“缩短”,以改变其卓越周期而加速度值不变也可以用数字滤波的方法滤去某些频率成分,改变谱的形状。另外,为了在计算中得到结构的最大反应,也可以根据建筑结构基本自振周期,调整实际地震记录的卓越周期,使二者接近。这种调整的结果,改变了实际地震记录的频率结构,从物理意义上分析是不合理的。 另外,在测定场地土和建筑结构的卓越周期时,运用不同的测试仪器和测试技术,往往得到不同的结果。即使是对同一个测试结果,在频谱上确定卓越周期时,不同的分析方法也会导致不同的结果。有的选取谱的第一个峰值所对应的周期作为卓越周期,有的选最大峰值时的,也有的取某一段周期等,很不一致。对如何确定地震加速度记录的卓越周期,也是各行其是,有的利用加速度反应谱,有的用伪速度谱,有的用富氏谱,结果当然是不一样的。上述各种作法在工程中引起了一些混乱。 王亚勇认为,用脉动测试方法测定场地土和结构的卓越周期及自振周期时,应采用速度摆型或加速度摆型的地震仪测定地运动和结构振动,然后计算其富氏谱或功率谱,以谱的最大峰值所对应的周期作为卓越周期和自振周期比较合适。反应而相应地根据记录的位移谱或速度谱。 这也就是所谓的滤波修正。可按要求设计滤波器,对地震波进行时域或频域的滤波修正。这样修正的地震资料不仅卓越周期满足要求,功率谱的形状和面积也可控制。卓越周期修正。将地震波的离散步长按人为比例改变,

地震波的激发和接收

地震波的激发和接收 论文提要 在地震勘探的野外工作中,第一步要用人工方法激发地震波,为了适应各种地表条件及具体工作特点,震源及激发方式是多种多样的。使用与地震勘探的震源基本上分为两大类。一类是炸药震源,另一类是非炸药震源。目前以炸药震源为主。地震波的接受问题就是使用专门的仪器设备,采用合适的工作方法,把地震波传播情况纪录下来。当我们激发地震波时,既产生有效波,也会产生干扰波。人们往往利用有效波和干扰波的差异,在野外条件下采用不同的仪器手段和观测方式,来压制干扰波、突出有效波的。 正文 一、地震波的激发 1. 对激发的要求 激发的有效地震波要有足够强的能量,良好的频谱特性和较高分辨能力,这样才能查明地下几千米深度范围的一整套地层的构造形态。此外还须指出,在激发出有效波的同时还会产生各种各样的波,如干扰波,异常波等。应使地震有效波具有较强能量、显著的频谱特性和较高的分辨能力。以利于纪录有效波。地震勘探的震源基本上分为两大类型,一类是炸药震源,另一类是非炸药震源。目前陆上主要以炸药震源、可控震源、气动震源为主,海上用电火花震源、空气枪震源、无气泡蒸汽枪震源等,其中炸药震源 是最常用的。因此,以炸药震源为例,介绍地震波激发试验。 2. 陆上用炸药震源 从20年代开始到现在,地震勘探方法一直采用炸药为主要震源。 炸药震源激发的效果主要取决于井深,药量、激发岩性因素的选择与使用,因此,激发岩性试验阶段不但要进行干扰波的调查,观测因素的选择,还要进行激发因素的试验。 1) 激发岩性 爆炸时所产生的波的频率谱很大程度上决定于激发岩石的物理性质。若在松软的干燥岩层(如砂层)或松散的岩层(如淤泥)中爆炸,频率很低,爆炸能量大部分被松散岩层所吸收,会产生极高频率,这种高频的振动很快被吸收掉,而且在爆炸点周围产生很大破碎带,转换成弹性能量不多,因此,激发岩性 应选取潮湿的可塑性岩层(如胶泥、粘土、湿 砂)。对大庆黑鱼泡地区岩性试验,见图 2.4.1。激发点选在灰胶泥层。 2) 激发深度 关于激发深度,以反射波来说,要选在 潜水面以下,最好是潜水面以下3~5米的粘 土层,或泥岩中爆炸,这样可使激发的频谱适 中,且由于激发点离上面的潜水面不远,潜水 面又是一个强反射界面,激发的能量由于潜水 面强烈反射作用大部分往下传播,从而增强有 效波的能量,减少了干扰波的能量,通常采用 民用水井、捞取钻井岩样或通过电测井、也可用浅层折射波法求低速带深度来进行潜水面深度、激发岩性的调查。也可在试验阶段分别选择固定药量,井深来进行井深试验获

新输入地震波

[结构分析] 地震波输入的问题 三向输入简化后的单向输入 首先,将三个方向的地震加速度放到一个文本文件里,如accexyz.txt,在这个数据文件里共放三列数据,每列为一个方向的地震加速度值,这里仅给出数据文件中前几行的数据:-0.227109E-02 -0.209046E+00 0.467072E+01 -0.413893E-02 -0.168195E+00 0.261523E+01 -0.574753E-02 -0.157890E+00 0.809014E-01 -0.731227E-02 -0.152996E+00 0.119975E+01 -0.876865E-02 -0.138102E+00 0.130902E+01 -0.101067E-01 -0.131582E+00 0.143611E+00 ....................... 然后,再建一个文本文件用来存放三个方向的地震加速度时间点,如time.txt,在这个数据文件里仅一列数据,对应于加速度数据文件里每一行的时间点,这里给出数据文件中前几行数据: 0.100000E-01 0.200000E-01 0.300000E-01 0.400000E-01 0.500000E-01 0.600000E-01 ....................... 编写如下的命令流文件,并命名为acce.inp *dim,ACCEXYZ,TABLE,2000,3 !01行 *vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行 (3e16.6) !03行 *vread,ACCEXYZ(1,0),time,txt !04行 (e16.6) !05行 ACCEXYZ(0,1)=1 !06行 ACCEXYZ(0,2)=2 !07行,同上 ACCEXYZ(0,3)=3 !08行,同上 finish /SOLU ANTYPE,trans btime=0.01 !定义计算起始时间 etime=15.00 !定义计算结束时间 dtime=0.01 !定义计算时间步长 *DO,itime,btime,etime,dtime time,itime AUTOTS,0 NSUBST,1, , ,1 KBC,1

时程分析时地震波的选取及地震波的反应谱化

时程分析时地震波的选取及地震波的反应谱化 发表时间:2017-12-29T15:40:37.810Z 来源:《防护工程》2017年第22期作者:金林飞 [导读] 目前我国规范要求结构计算中地震作用的计算方法一般为振型分解反应谱法。 常州市规划设计院江苏常州 213002 摘要:目前我国规范要求结构计算中地震作用的计算方法一般为振型分解反应谱法。时程分析法作为补充计算方法,在不规则、重要或较高建筑中采用。进行时程分析时,首先面临正确选择输入的地震加速度时程曲线的问题。时程曲线的选择是否满足规范的要求,则需要首先将时程曲线进行单自由度反应计算,得到其反应谱曲线,并按规范要求和规范反应谱进行对比和取舍。本文通过介绍常用的数值计算方法及计算步骤,实现将地震加速度时程曲线计算转化成反应谱曲线,从而为特定工程在时程分析时地震波的选取提供帮助。 关键词:时程分析,地震波,反应谱,动力计算 1 地震反应分析方法的发展过程 结构的地震反应取决于地震动和结构特性。因此,地震反应分析的水平也是随着人们对这两个方面认识的深入而提高的。结构地震反应分析的发展可以分为静力法、反应谱法、动力分析法这三个阶段。在动力分析法阶段中又可分为弹性和非弹性(或非线性)两个阶段。 [1] 目前,在我国和其他许多国家的抗震设计规范中,广泛采用反应谱法确定地震作用,其中以加速度反应谱应用得最多。反应谱是指:单自由度弹性体系在给定的地震作用下,某个最大反应量(如加速度、速度、位移等)与体系自振周期的关系曲线。反应谱理论是指:结构物可以简化为多自由度体系,多自由度体系的地震反应可以按振型分解为多个单自由度体系反应的组合,每个单自由度体系的最大反应可以从反应谱求得。其优点是物理概念清晰,计算方法较为简单,参数易于确定。 反应谱理论包括如下三个基本假定:1、结构物的地震反应是弹性的,可以采用叠加原理来进行振型组合;2、现有反应谱假定结构的所有支座处地震动完全相同;3、结构物最不利的地震反应为其最大地震反应,而与其他动力反应参数,如最大值附近的次数、概率、持时等无关。[1] 时程分析法是对结构物的运动微分方程直接进行逐步积分求解的一种动力分析方法。由于此法是对运动方程直接求解,又称直接动力分析法。可直接计算地震期间结构的位移、速度和加速度时程反应,从而描述结构在强地震作用下弹性和非弹性阶段的内力变化,以及结构构件逐步开裂、屈服、破坏甚至倒塌全过程。 根据我国《建筑抗震设计规范》(GB5011-2010)(以下简称《抗规》)第5.1.2-3条要求,特别不规则的建筑、甲类建筑和表5.1.2-1所列高度范围的高层建筑,应采用时程分析法进行多遇地震下的补充计算。此外《高层建筑混凝土结构技术规程》(JGJ3-2010) (以下简称《高规》)第4.3.4条也有相关要求。 2 时程分析时地震波的选取要求 在进行时程分析时,首先面临地震波选取的问题。所选的地震波需要符合场地条件、设防类别、震中距远近等因素。《抗规》对于地震波的选取主要有以下几点要求: 1、当取三组加速度时程曲线输入时,计算结果宜取时程法的包络值和振型分解反应谱法的较大值;当取七组及七组以上的时程曲线时,计算结果可取时程法的平均值和振型分解反应谱法的较大值(其中实际强震记录的数量不应少于总数的2/3)。 2、弹性时程分析时,每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计算结果的65%,多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80%。 3、多组时程曲线的平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。根据规范条文说明,所谓“统计意义上相符”指的是,多组时程波的平均地震影响系数曲线与振型分解反应谱法所用的地震影响系数曲线相比,在对应于结构主要振型的周期点上相差不大于20%。但计算结果也不能太大,每条地震波输入计算不大于135%,平均不大于120%。 4、时程曲线要满足地震动三要素的要求,即频谱特性、有效峰值和持续时间均要符合规定。其中频谱特性可用地震影响系数曲线表征,依据所处的场地类别和设计地震分组确定;加速度的有效峰值按《抗规》表5.1.2-2中所列地震加速度最大值采用;输入的地震加速度时程曲线的有效持续时间,一般从首次达到该时程曲线最大峰值的10%那一点算起,到最后一点达到最大峰值的10%为止。有效持续时间一般为结构的基本周期的(5~10)倍,即结构顶点的位移可按基本周期往复(5~10)次。 根据以上《抗规》对地震波选取的要求,在工程中选取地震波时,首先需将地震波进行单自由度动力求解,得出相应的反应谱曲线,并与规范反应谱所用的地震影响系数曲线进行对比。对比频谱特性,查看是否在统计意义上相符,若相符则可将此地震波用于结构计算进行底部剪力验证,否则则需要更换地震波重新分析。 3 求解地震波单自由度体系反应的数值计算方法 对一般动力荷载反应的求解方法可分为叠加法和逐步法两类。 叠加法包括时域分析和频域分析,总反应计算采用独立反应贡献的组合。时域方法中,荷载被考虑为短暂持续时间的脉冲序列,由每个脉冲自由振动反应的独立贡献得到后续时间的总反应;频遇方法中,假设荷载为周期的,并用Fourier变换为离散的谐振分量Pn。再由这些荷载分量乘以结构的频率反应系数Hn,得到与其相应的结构谐振反应分量Vn。最后,由组合谐振反应分量(Fourier逆变换)获得结构的总反应。由于结果都使用了叠加,因此不适用于非线性反应分析。[2] 逐步法有很多种,但所有方法都是将荷载和反应历程分成一系列时间间隔或“步”。在每步期间均以此步开始时存在的初始条件(位移、速度和加速度)和该步期间的荷载历程来计算反应。因此每步反应是一个独立的分析问题。[2] 积分法是逐步法中的一般性的方法之一。积分法对每一时间步,从初始最终条件应用积分向前进一步,速度的变化依赖于加速度历程的积分,而位移的变化依赖于相应的速度积分。 用积分法进行分析时,首先需要假设在时间步内加速度是如何变化的。根据假设的不同,积分法可分为“基于常平均加速度”法(如

地震波的选取方法

地震波的选取方法 2010-10-20 22:32:00| 分类:默认分类|举报|字号订阅 建筑抗震设计规范(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程序中提供将地震波转换为绝对加速度反应谱和拟速度反应谱的功能(工具>地震 波数据生成器,生成后保存为SGS文件),用户可利用保存的SGS文件(文本格式文件)根据上面所 述方法计算Sv、Sa、Tg。通过Tg值可判断该地震波是否适合当地场地和地震设计分组,然后将 抗震规范中表5.1.2-2中的EPA值与Sa相比求出调整系数,将其代入到地震波调整系数中。将地 震波转换为绝对加速度反应谱和拟速度反应谱时注意周期范围要到6秒(建筑抗震规范规定)。 建筑抗震设计规范5.1.2条中规定,采用时程分析方法时,应按照场地类别和设计地震分组 选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线,其平均地震影响系数曲

地震波的吸收

五、 地震波的吸收 1.对吸收的概括说明 弹性介质是实际介质的近似,实际介质对地震波有吸收。使波传播时能量很快衰减。地震波的能量不可逆地转化成了热能散发掉叫做吸收。 目前公认岩石颗粒之间的内摩擦力是吸收的主要原因,这种内摩擦力也叫粘滞力。 不同的介质中,内摩擦力所遵守的规律不同,其中沃尹特(Voigt )的假设比较有代表性。 2.沃尹特(Voigt )假设下的波动方程 (1)沃尹特假设 前面讨论的完全弹性介质中的波动方程是在应力与应变成正比(虎克定律)的条件下得到的。 沃尹特假设:应力与应变的关系包括两部分,一部分是应力与应变成正比的满足虎克定律的应变,另一部分是应力与应变的时间变化率成比例的粘滞效应。 (2)沃尹特假设的粘弹介质中的波动方程 t u t u t u ???+???+?+?+=??222231)(ηθημθμλρ (6.1-45) k z j y i x grad ??+??+??==? η是粘滞系数。 对(6.1-45)求散度,可得纵波的波动方程为: θημλθρ222]34)2[(???++=??t t (6.1-46) 对(6.1-45)求旋度,可得横波的波动方程为: u rot t u rot t 222)(???+=??ημρ (6.1-47) 3.沃尹特粘弹介质中一维谐波方程的求解

(1)求解 分析平面谐波沿X 方向传播,以纵波为例。 假设纵波的位移位是: )(0)(0),(Kx t j v x t j e e t x --==ωωφφφ (6.1-48) x u z w y v x u ??=??+??+??=θ (注:波沿x 方向传播,y v ??=0,z w ??=0) 又x u ??= φ (注:w v u x grad u ++=++??==00φφ) φφθ222K x x u -=??=??=∴ (6.1-49) 将(6.1-49)代入(6.1-46),得 22234)2(K j K ηωμλρω++= ωημλρωηωμλρω'++=++=∴j j K )2(3 )2(2 2 2 (ηη34'=) (6.1-50) 将上式有理化,得 }])2[(])2[(2{])2[(2122221222212222 2ωημλωηωημλμλωημλρω'++'-'+++'++=j K (6.1-51) 222)2(ωημλ'++ ωη' β μλ2+ βωημλμ λcos ])2[(221222='+++ βωημλω ηsin ])2[(21222='++' βμ λωηtan 2=+'

相关文档
最新文档