1．均方根振幅（RMS Amplitude ）
均方根振幅是将振幅平方的平均值再开平方。由于振幅值在平均前平方了，因此，它对特别大的振幅非常敏感。适合于地层的砂泥岩百分比含量分析，也用于地层岩性相变分析，计算薄砂层厚度，识别亮点、暗点，指示烃类显示，识别火成岩等特殊岩性。
2．平均绝对值振幅（Average Absolute Amplitude ）
平均绝对值振幅没有均方根振幅那样，对特别大的振幅敏感。
适于地层的岩性变化趋势分析，地震相分析，也可用于地层岩性相变分析，计算薄砂层厚度，识别亮点、暗点，指示烃类显示，识别火成岩等特殊岩性。
3．最大波峰振幅（Maximum Peak Amplitude ）
最大波峰振幅的求取方法是，对于每一道，PAL 在分析时窗里做一抛物线，恰好通过最大正的振幅值和它两边的两个采样点，沿着这曲线内插可得到最大波峰值振幅值。
√
PAL画一个使这三个采样点适合曲线并且
沿这一曲线确定出最大值。
最大波峰振幅= 125
最大波峰振幅是分析时窗内的最大正振幅，最适合绘制层序内或沿着特定的反射体上的振幅异常图；这些异常可能是由于气体和流体的聚集，不整合，或是调谐效应而引起的。
适于沿某一层面进行储层分析，也可用于地层岩性相变分析，计算薄砂层厚度，识别亮点、暗点，指示烃类显示，识别火成岩等特殊岩性。
4．平均波峰振幅 (Average Peak Amplitude)
平均峰值振幅是对每一道在分析时窗里的所有正振幅值相加，得到总数除以时窗里的正振幅值采样数得到的。
适合研究某一层的岩性变化，也可用于地层岩性相变分析，计算薄砂层厚度，识别亮点、暗点，指示烃类显示，识别火成岩等特殊岩性。
5．最大波谷振幅 (Maximum Trough Amplitude)
最大波谷振幅的求取方法是，对于每一道，PAL在分析时窗里做一抛物线，恰好通过最大负的振幅值和它两边的两个采样点，沿着这曲线内插可得到最大波谷振幅值。
PAL 画一个适合这三个采样点的曲线
并且沿着这一曲线确定出最大值。
最大波谷振幅= |-90| = 90
最大波谷振幅是指定分析时窗内的最大负振幅，所以最适合绘制层序内或沿着特定的反射体上的振幅异常图；这些异常可能是由于气体和流体的聚集，不整合，或是调谐效应而引起的。
它与最大峰值振幅极性相反，应用相同。
6．平均波谷振幅(Average Trough Amplitude)
平均波谷振幅是对每一道在分析时窗里的所有负振幅值相加，得到总数除以时窗里的负振幅值采样数得到的。
与平均波峰振幅极性相反，应用相同。
7．最大绝对值振幅 (Maximum Absolute Amplitude)
计算每道的最大绝对值振幅的求取方法是，首先在分析时窗内计算出波峰和波谷的值，得出最大的波峰或波谷值，然后，PAL画一抛物线，恰好通过最大波峰或波谷振幅值和它两边的两个采样点，沿着这曲线内插可得到最大绝对值振幅值。
适合岩性分析，砂岩百分比研究，也可用于地层岩性相变分析，计算薄砂层厚度，识别亮点、暗点，指示烃类显示，识别火成岩等特殊岩性。
PAL 画一个适合这三个采样点的曲线
并且沿着这一曲线确定出最大值。
最大绝对值振幅= 123.6
8．总绝对值振幅 (Total Absolute Amplitude)
总绝对值振幅是计算确定时窗内的所有道的绝对值振幅值。
总绝对值振幅 = 振幅绝对值之和
= 1045
适合大套地层变化趋势分析和某一岩性的含量分析，也可用于地层岩性相变分析，计算薄砂层厚度，识别亮点、暗点，指示烃类显示，识别火成岩等特殊岩性。
9．总振幅 (Total Amplitude)
每一道的总振幅是，在层内对采样点求取总的振幅值。
适合大套地层变化趋势分析，也可用于地层岩性相变分析，计算薄砂层厚度，识别亮点、暗点，指示烃类显示，识别火成岩等特殊岩性。
总振幅= 振幅之和
= 559
三向输入简化后的单向输入 首先，将三个方向的地震加速度放到一个文本文件里，如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可以就看成有许多
SIS 软件软件技术应用技术应用技术应用之一之一 斯伦贝谢伦贝谢科技服务科技服务科技服务（（北京北京））有限公司 2007年3月 GeoFrame 地震属性分析和应用地震属性分析和应用
1 地震属性分析和应用 应用地震属性开展储层横向预测是地震资料综合解释的重要研究内容。随着地球物理理论、数学理论的不断发展，通过各种计算方法能够提取和分析的地震属性越来越多，如何从众多的地震属性中选择能够反映客观地质现象的属性对目的层储层开展分析，这是地球物理人员在实际工作中面对的一个主要问题。 GeoFrame 综合地学平台为地球物理人员开展储层横向预测研究提供了一套完善的工具。SATK 、SeisClass 、LPM 以及GeoViz 的组合应用，可以帮助研究人员完成从属性提取、属性优化、定性分析到定量计算的储层预测全过程。本文重点阐述GeoFrame 储层预测的基本思路及地震属性的地质应用。 1、地震属性储层预测的基本思路 地震地层学原理假定，地震剖面上的反射波同相轴具有年代分界面的意义，要研究地层岩性和沉积相主要依据的是地震反射特征及其横向变化，也就是地震属性的变化,这是应用地震属性进行储层预测的基本理论依据。 应用地震属性进行储层横向预测要解决的主要问题是多解性问题，即：一种地震属性参数的变化受多种地质因素的影响，而一种地质现象的改变，也会造成多种地震属性的异常。 因此，在对地震属性分析预测过程中，如何从众多的地球物理参数中选取能反映地质特征变化的参数，是地震属性预测的主要问题。实际工作表明，必须做好以下两项工作： ① 正确认识地震属性 正确认识地震属性是做好属性预测的基础，不同的地震属性参数，它的地球物理含义、数学含义不一样，反映的地质规律也不一样。如：半时能量和总能量，尽管都是振幅类参数，但具体的展布规律却不一样（图1）。 图1 1 相同地区相同地区相同地区半时能量半时能量半时能量和和总能量总能量对比图对比图对比图 半时能量半时能量（（Energy half-time ） 总能量总能量（（Total Energy ）
时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系，能真实地反应结构地震相应随时间变化的全过程，是抗震分析的一种重要方法[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的乘积即为单自由度体系的地震影响系数。 因此，从原理上讲，时程分析是将地震波的加速度时程曲线作用到结构的基础约束处，得到上部结构的各种地震反应。但是在不同的软件中，其实现方
常用地震属性的意义 欧阳家百（2021.03.07） 地震反射波来自地下地层，地下地层特征的横向变化，将导致地震反射波特征的横向变化，进而影响地震属性的变化，因此，地震属性中携带有地下地层信息，这是利用地震属性预测油气储层参数的物理基础。随着地震属性处理及提取技术的大量涌现，属性种类多达几百种，实际应用人员应用起来遇到了很大困难，迫切需要按实用的角度，总结各地震属性参数与储层特征参数间的内在联系，为进一步研究建立地震信息与储层参数之间的关系提供可靠的前提条件，做到信息提取有方向、有目标。为了达到这一目的，首先按类别较全面总结了目前常用地震属性，从算法开始，分析了各属性所表达的在地震波波形上的意义，从正向上分析地震属性变化与油气储层特征变化的关系，进而探讨总结了它的潜在地质应用。 1、属性体、属性剖面 这类属性是按剖面（或体）处理的，是一个体文件（或剖面文件），属性值对应空间位置，即（x、y、t0、属性值），可以用于常规地震剖面的方式显示与使用，常用的属性有：相干体（方差体、相似体等）、波阻抗、道积分数据体，经希尔伯特变换得到的瞬时属性体、倾角、倾向数据体等，这些属性体可以直接应用于解释，也可以用解释层位提取出来转变为属性层，下表为常用属性体属性意义及潜
2、沿层地震属性 这种属性是用解释层位在地震数据体（剖面）中提取出来的属性，它的数值对应一个层位或一套地层，每个属性值对应一个x、y 坐标。提取方式有两类：沿一个解释层开一个常数时窗，在此时窗内
提取地震属性，提取方式有4种（图2-1a）。用两个解释层提取某一段地层对应的地震属性，提取方式也有4种（图2-1b）。 常用地震属性的计算方法总结如下： (1)、均方根振幅（RMS Amplitude） 均方根振幅是将振幅平方的平均值开平方。由于振幅值在平均前平方了，因此，它对特别大的振幅非常敏感。 (2)、平均绝对值振幅（Average Absolute Amplitude） 平均绝对值振幅没有均方根振幅那样，对特别大的振幅敏感。 (3)、最大波峰振幅（Maximum Peak Amplitude） 最大波峰振幅的求取方法是，对于每一道，PAL在分析时窗里做一抛物线，恰好通过最大正的振幅值和它两边的两个采样点，沿着这曲线内插可得到最大波峰值振幅值。 PAL画一个使这三个采样点适合曲线并且 沿这一曲线确定出最大值。 MaximumPeak Amplitude = 125 （4）、平均波峰振幅 (Average Peak Amplitude) 平均峰值振幅是对每一道在分析时窗里的所有正振幅值相加，得到总数除以时窗里的正振幅值采样数得到的。 （5）、最大波谷振幅 (Maximum Trough Amplitude) 最大波谷振幅的求取方法是，对于每一道，PAL在分析时窗里做一抛物线，恰好通过最大负的振幅值和它两边的两个采样点，沿着这曲线内插可得到最大波谷振幅值。 PAL 画一个适合这三个采样点的曲线 并且沿着这一曲线确定出最大值。
地震波使用说明 此目录下提供了四类场地土的地震波时程曲线和上海人工波。 按照场地土类型（1，2，3或4），选择时程曲线。在定义时程工况时，对于多遇或罕遇地震，按比例调整时程曲线的最大值。中国抗震规范规定，作为抗震计算中底部剪力法和振型分解反应谱法的补充方法，对于特别不规则，特别重要的和较高的结构应采用时程分析法进行多遇地震下的补充计算。 可取多条时程曲线的计算结果的平均值与振型分解反应谱法计算结果的较大值。 采用时程分析法时，应咱建筑场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。 其加速度时程最大值可按规范中对于多遇和罕遇地震在不同烈度下的值。 弹性时程分析时，每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计算结果的65％，多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80％ 。 可使用弹塑性时程分析法计算罕遇地震下结构的变形。 时程分析是一个承受随时间变化的指定荷载结构的逐步动态反应分析，可以是线性或非线性的。 此章对时程分析进行一般的描述，特别是线性时程分析。 定义时程函数 用户可使用“从文件中添加函数”，导入已定义的文本文件，即实测的时程曲线；也可使用程序内置的时程函数。
时程函数定义对话框 时程函数定义对话框中的条目解释如下： ?函数名 通过在编辑框中直接键入以指定或修改时程函数的名称。 ?函数文件 1.在函数文件域点击浏览按钮以调出一个对话框，在此可找出包含时程函数的 文本文件名。注意文件名显示在文件名框中 2.在 "要跳过的标题行" 编辑框中输入一个希望ETABS在文本文件中跳过的 行数。 3.在 "每行要跳过的前缀字符" 编辑框中输入一个希望ETABS在文本文件中 每行要跳过的字符数。 4.在 "每行的点数" 编辑框中输入一个数告诉ETABS文本文件每行的绘图点 数。
施加地震波： 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 ………………
【全文】地震属性分析技术综述 [摘要] 地震属性是从地震资料中提取的隐藏有用信息，因而地震属性分析技术近几年在油气勘探开发中得到了广泛的应用与研究。本文对地震属性分析技术的发展状况进行了归纳、总结，简单阐述了地震属性分析技术的在不同时期所用到的基本原理和方法。特别对新地震属性进行了具体介绍。最后对该技术进一步的研究工作进行了总结和展望。 摘要:在勘探和开发周期的各个阶段，地震资料在复杂油藏系统的解释过程中，扮演着至关重要的角色。然而，缺少一种有效地将地质知识应用于地震解释中的上具。随着一系列属性新技术的出现，对地震属性进行充分研究，就给地质家提供了快速地从三维地震数据中获得地质信息的能力。尤其在用常规解释手段难以识别日的储层的情况下，属性分析技术更是给地质上作人员指出了新的方向。 [关键词] 地震属性储层预测叠前数据叠后数据 关键词:储层;波形分析;地震属性 1.引言 地震属性是指叠前或叠后的地震数据经过数学变换而导出的有关地震波的几何形态、运动学特征、动力学特征和统计学特征的特殊度量值。地震属性的发展大致从20世纪60年代的直接烃类检测和亮点、暗点、平点技术开始，经历了70年代的瞬时属性(主要是振幅属性)和复数道分析，90年代的多维属性(特别是相干体属性)分析，21世纪的地震相分析等阶段[1一SJ。随着地震属性分析技术的发展与研究，该技术已广泛应用于储层预测、油气藏动态监测、油气藏特征描述等领域，并取得了很好的效果。总之，地震属性分析技术可以从地震资料中提取隐藏其中的多种有用信息，这为油气勘探与开发提供了丰富宝贵的资料，也为解决复杂地质体评价提供了实用的分析手段。因此，对该技术进行深人调查研究具有很强的现实意义。 地震属性是指从地震数据中导出的关于儿何学、运动学、动力学及统计特性的特殊度量值。它可包括时问属性、振幅属性、频率属性和吸收衰减属性，不同的属性可指示不同的地质现象。地震属性分析则是从地震资料中提取其中的有用信息，并结合钻井资料，从不同角度分析各种地震信息在纵向和横向上的变化，以揭示出原始地震剖面中不易被发现的地质异常现象及含油气情况。 地震属性分析技术的研究已由线、面信息扩展到三维体信息，从分类提取扰化发展为一项系统的应用技术。随着地震技术的日趋成熟，地震属性技术近儿年也发展迅速，其中有多属性联合解释技术、波形分析技术、吸收滤波技术等。应用地震属性分析技术去完善勘探生产中的油藏描述工作，已经成为油藏地球物理的核心内容。利用地震属性分析技术预测岩性和有利储集体，描述油藏特征及孔隙度变化，寻找难以发现的隐蔽油区，以至于监测流体运动和进行其它综合研究，一直是石油工作人员追求的目标。 1波形分析技术的研究与应用 通常的层段属性只是表示了某儿个地震信号的物理参数(振幅、相位、频率等)，但它们没有一个能够单独描述地震信号的异常，而地震信号的任何物理参数的变化总是对应着反映地震道形状的变化，所以，研究和分析地震资料中代表各种属性总体特征的地震道形状(波形)，应该能有非常不错的效果[，]。 1. 1波形分析技术的原理及处理过程
对于地震波的输入，可以把荷载记录作成文件，利用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】功率谱，反应谱和傅里叶谱，地震波选取，地震持续时间确定功率谱 功率谱是功率谱密度函数的简称。对于一般情况的随机振动，其时间历程具有明显的非周期性，具有连续的多种频率成分，每种频率有对应的功率或能量，用图像来表示这种关系，称为功率在频率域内的函数，简称功率谱密度。加速度功率谱是对地震动加速度时程进行快速傅里叶变换（FFT）得到的[1]。 对于非平稳随机过程，功率谱密度的单位是G的平方/频率。G指的是随机过程。 对于加速度功率谱，加速度的单位是m/s2，则功率谱密度的单位是（m/s2）2/Hz，Hz的单位是1/s，故加速度功率谱密度的单位为m2/s3。加速度功率谱密度函数曲线下方的面积代表随机加速度的总方差，即加速度功率谱可以理解为“随机加速度方差的密度分度”。 参考文献 [1] 庄表中. 随机振动入门.科学出版社，1981. 反应谱和傅里叶谱 反应谱（earthquake response spectrum），是单自由度弹性系统对于某个实际地震加速度的最大反应（可以是加速度、速度和位移）和体系的自振特征（自振周期或频率和阻尼比）之间的函数关系。 反应谱是地震工程中分析结构和设备在地震中的性能的非常有用的工具，因为许多主要表现为简单的振荡器(也称为单自由度系统)。因此，如果能找出结构的固有频率，那么建筑的峰值响应可以通过从地面响应谱中读取相应频率的值来估计。在地震区域的大多数建筑规范中，这个值构成了计算结构必须抵抗的力的基础(地震分析)。如前所述，地面响应谱是在地球自由表面所做的响应图。如果建筑物的响应与地面运动(共振)的组成部分“协调”，可能会发生重大的地震破坏，这些成分可以从响应谱中识别出来。 傅里叶谱，全称为傅里叶振幅谱。地震波是在时间上连续的随机过程，地震
中国地质大学（北京） 地震属性综述 报告名称: 地震属性综述 学生姓名：王丹 学号：2010120052 所在院（系）：地球物理与信息技术学院
地震属性分类及其地质意义 地震勘探是在地表激发人工震源，由震源所引起的震动以地震波的形式向地下传播，并在一定的条件下向上反射传回地表，然后由地表的仪器(检波器)记录反射回来的地震波，从而得到地震记录(也叫地震资料)；之后对地震资料进行相关的处理与解释便可以间接地反映和得到地下相关信息。由于地下介质是地震波传播的载体，所以地下介质的物理性质，如岩性、孔隙度、密度以及流体性质等都会对传播中的地震波产生影响，如地震波的能量、波形、振幅、频率、相位等将在传播过程中发生变化。而这种影响和变化又将在地震记录中保留相应的信息。所以，通过对地震记录(地震资料)的“深加工”或者特殊处理，将会从地震资料中获取更多的有用信息以为地质服务。在早期的油气资源勘探中，地震勘探的目标主要是寻找地下有利的大尺度的构造圈闭，所以只需利用有限的地震资料信息便可达到目的。但是，随着油气勘探与开发难度的加大，人们迫切地需要更多地了解地下地层的岩性、流体性质等信息。这就促使人们运用新的技术和思想去从地震资料中发掘出更多的有用信息。从而，也就推动了地震属性技术的出现与发展。地震属性技术延伸了人类的视觉,从而有助于人们发现更多的隐藏于地震资料中的信息，也有助于人们从多角度去获取和分析地下地质信息，从而实现对地下地质的充分与准确认识。 1地震属性的发展与分类 随着油气勘探、开发工作的深入，也为了充分、有效地利用获取不易的地震资料，现今的地震解释人员需要从地震数据中提取越来越多的信息，然后利用这些信息综合解释地下构造、地层和岩性特征以及流体性质，最终定义精确的油藏模型，用于钻井决策、估计地质储量和可采储量。由于生成地震属性是获取所需信急的一条重要捷径，因此，长期以来地震属性技术一直是地震特殊处理和解释的主要研究内容。 地震属性是叠前或者叠后地震数据，经数学变换而导出的有关地震波的几何形态、运动学特征、动力学特征和统计学特征。长期以来以来地震数据的使用仅仅局限于对地震波同相轴的拾取，以实现面对油气储集体的几何形态、构造特征的描述。但是地震数据中隐藏着更加丰富的有关岩性、物性及流体成分等相关信
地震波的定义 地震是地壳的一切颤动，是一种自然现象。其主要能源来自地球的内部，是由地球内部自然力冲击引起的。地壳或地幔中发生振动的地方称为震源。震源在地面上的垂直投影称为震中。震中到震源的距离称为震源深度。地震波是指从震源产生向四外辐射的弹性波。地球内部存在着地震波速度突变的基干界面、莫霍面和古登堡面，将地球内部分为地壳、地幔和地核三个圈层。 发生原理 英文seismic wave.由地震震源发出的在地球介质中传播的弹性波。地球内 地震波 部存在着地震波速度突变的基干界面、莫霍面和古登堡面，将地球内部分为地壳、地幔和地核三个圈层。地震震源发出的在地球介质中传播的弹性波。地震发生时，震源区的介质发生急速的破裂和运动，这种扰动构成一个波源。由于地球介质的连续性，这种波动就向地球内部及表层各处传播开去，形成了连续介质中的弹性波。 概念介绍 地震波是指从震源产生向四外辐射的弹性波。地球内部存在着地震波速度突变的基干界面、莫霍面和古登堡面，将地球内部分为地壳、地幔和地核三个圈层。传播方式 地震波按传播方式分为三种类型：纵波、横波和面波[1]。纵波是推进波，地壳中传播速度为5．5~7千米/秒，最先到达震中，又称P波，它使地面发生上下振动，破坏性较弱。横波是剪切波：在地壳中的传播速度为3.2～4.0千米/秒，第二个到达震中，又称S波，它使地面发生前后、左右抖动，破坏性较强。面波又称L波，是由纵波与横波在地表相遇后激发产生的混合波。其波长大、振幅强，只能沿地表面传播，是造成建筑物强烈破坏的主要因素。 纵波和横波 现象介绍 我们最熟悉的波动是观察到的水波。当向池塘里扔一块石头时水面被扰乱，
利用PEER网站选取地震波的方法 云南省土木建筑学会建筑结构专业委员会 PEER是Pacific Earthquake Engineering Research Center（太平洋地震工程研究中心）的简称，设立在美国的加州大学伯克利分校（University of California, Berkeley）。由其运作的PEER Ground Motion Database（PEER地震动数据库）提供了大量的世界各地的地震记录，可自由下载。该网站提供了丰富的查询手段，可按距离、场地、震源类型等条件选择地震记录，也提供了按目标反应谱选择的手段。 该网站是https://www.doczj.com/doc/736221842.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”。
Landmark主要地震属性及其地质意义利用地震进行储层预测时主要从振幅属性及其延伸属性出发，分析属性的变化特征，然后与钻井和地质进行标定，赋予属性地质意义。 为了将已知井上的岩性信息，在整个工区进行有效的外推，需要优选出在该区对岩性参数和含油气性反映敏感的属性，我们通过两个层次来完成这一个工作。第一个层次是选择对岩性变化相对敏感的地震属性，这部分工作在属性提取时已完成，其最基本的理论基础是：时间派生的属性有利于对构造的细节进行解释；振幅和频率派生的属性用于解决地层和储层特征； 一般认为振幅是最稳健和有价值的属性；频率属性更有利于揭示地层的细节； 混合属性包含振幅和频率的因素，因此更有利于地震特征的测量；同时在对所提取的地震属性的物理意义的理解也有助于对地震属性的提取第二个层次是使用数学和信息学的方法优选属性。“地震属性和井数据采样伪相关在独立的井数据较少或者参加考虑的独立的地震属性过多时产生的概率较大”（CYNTHIA T. KALKOMEY），由于对于该区已知的独立井信息多数情况下较少，勉强满足统计分析的样本要求，单纯使用相关分析方法产生伪相关的概率较大，因此我们在经过第一个层次的筛选之后，采用数据相关和信息优化组合方法进行属性优选。 目前属性种类很多，属性软件也非常多，这里转列landmark软件中的PAL 属性，供大家参考选择使用：Average Reflection Strength 平均反射强度：识别振幅异常，追踪三角洲、河道、含气砂岩等引起的地震振幅异常；指示主要的岩性变化、不整合、天然气或流体的聚集；该属性为预测砂岩厚度的常用属性； Slope Half Time 能量半衰时的斜率：突出砂岩/泥岩分布的突变点；预测砂岩厚度的常用属性； Number of Thoughs 波谷数：可以有效的识别薄层，为预测砂岩厚度的常用属性；Average Trough Amplitude 平均波谷振幅：用于识别岩性变化、含气砂岩或地层。可以有效的区分整合沉积物、丘状沉积物、杂乱的沉积物等；预测含油气性的常用属性； Average Instantaneous Phase 平均瞬时相位：由于相位的横向变化可能与地
地震波的选取方法 建筑抗震设计规范(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后的)，然后将一组人
地震波的频率和振幅 时间:2010-06-05 20:18来源:unknown 作者:wowglad 点击:7次 2008年12月19日 地震波的频率和振幅 1、地震波的频谱及其分析 频谱：谐和振动的振幅和初相位则随频率的改变而改变的关系,统称为地震波的频谱。 频谱分 2008年12月19日 地震波的频率和振幅 1、地震波的频谱及其分析 频谱：谐和振动的振幅和初相位则随频率的改变而改变的关系,统称为地震波的频谱。 频谱分为： 振幅谱：振幅随频率变化的关系称为振幅谱。 相位谱：初相位随频率的变化关系称为相位谱。 作用：频率分析,根据有效波和干扰波的频段差异 ①指导野外工作方法的选择 ②给数字滤波和资料等工作提供依据。 频谱分析的方法： 为了研究地震波的频谱特征,可用傅立叶变换把波形函数a（t）变换到频率域中,得到振幅随频率的变化函数A(f),这个变换过程称之为频谱分析方法。 假设波形函数a(t) ------------------(1.3.1)--
--傅氏正变换 --------------------（1.3.2）-- --傅氏反变换 这两式是等价的,即A(f)与a(t)是一一对应的。 ① δ脉冲函数Aδ(t) ② 函数： ③ 函数： 可以看出：不同时间函数具有不同的频谱。 图1.3.52、地震波的频率特征 地震波是人工激发的振动,具有连续的频谱,如图1.3.6所示。
图1.3.6主频f0：振幅谱曲线极大值所对应的频率。 频带的宽度：若|A（f）|最大值为1，则可找|A（f）|=0.707的两个频率f1和f2,两者之差△f=f2-f1为频带宽度。 大量的实际观测和分析,各种不同类型的地震波的能量主要分布频带是不同的。如图1.3.7所示。 图1.3.7 3、地震波的振幅及其衰减规律 影响地震波激发和接收时振幅和波形的因素： ① 激发条件。 ② 地震波在传播过程中受到影响。 ③ 接收条件的影响。 ④ 其它如地下岩层界面的形态和平滑状态。
地震波的激发和接收 论文提要 在地震勘探的野外工作中，第一步要用人工方法激发地震波，为了适应各种地表条件及具体工作特点，震源及激发方式是多种多样的。使用与地震勘探的震源基本上分为两大类。一类是炸药震源，另一类是非炸药震源。目前以炸药震源为主。地震波的接受问题就是使用专门的仪器设备，采用合适的工作方法，把地震波传播情况纪录下来。当我们激发地震波时，既产生有效波，也会产生干扰波。人们往往利用有效波和干扰波的差异，在野外条件下采用不同的仪器手段和观测方式，来压制干扰波、突出有效波的。 正文 一、地震波的激发 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