三向输入简化后的单向输入 首先，将三个方向的地震加速度放到一个文本文件里，如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，3或4），选择时程曲线。在定义时程工况时，对于多遇或罕遇地震，按比例调整时程曲线的最大值。中国抗震规范规定，作为抗震计算中底部剪力法和振型分解反应谱法的补充方法，对于特别不规则，特别重要的和较高的结构应采用时程分析法进行多遇地震下的补充计算。 可取多条时程曲线的计算结果的平均值与振型分解反应谱法计算结果的较大值。 采用时程分析法时，应咱建筑场地类别和设计地震分组选用不少于二组的实际强震记录和一组人工模拟的加速度时程曲线，其平均地震影响系数曲线应与振型分解反应谱法所采用的地震影响系数曲线在统计意义上相符。 其加速度时程最大值可按规范中对于多遇和罕遇地震在不同烈度下的值。 弹性时程分析时，每条时程曲线计算所得结构底部剪力不应小于振型分解反应谱法计算结果的65％，多条时程曲线计算所得结构底部剪力的平均值不应小于振型分解反应谱法计算结果的80％ 。 可使用弹塑性时程分析法计算罕遇地震下结构的变形。 时程分析是一个承受随时间变化的指定荷载结构的逐步动态反应分析，可以是线性或非线性的。 此章对时程分析进行一般的描述，特别是线性时程分析。 定义时程函数 用户可使用“从文件中添加函数”，导入已定义的文本文件，即实测的时程曲线；也可使用程序内置的时程函数。
时程函数定义对话框 时程函数定义对话框中的条目解释如下： ?函数名 通过在编辑框中直接键入以指定或修改时程函数的名称。 ?函数文件 1.在函数文件域点击浏览按钮以调出一个对话框，在此可找出包含时程函数的 文本文件名。注意文件名显示在文件名框中 2.在 "要跳过的标题行" 编辑框中输入一个希望ETABS在文本文件中跳过的 行数。 3.在 "每行要跳过的前缀字符" 编辑框中输入一个希望ETABS在文本文件中 每行要跳过的字符数。 4.在 "每行的点数" 编辑框中输入一个数告诉ETABS文本文件每行的绘图点 数。
论地震勘探中的几种主要地震波 论文提要 地震勘探，就是通过人工方法激发地震波，研究地震波在地层中传播的情况，以查明地下地质构造，为寻找油气田或其它勘探目的服务的一种方法。也可以理解为就是利用地震子波从地下地层界面反射回地面时带回来的旅行时间和形状变化的信息，用以推断地下的底层构造和岩性。地震勘探在勘探已有的各种物探方法中，是最有效地方法。在地震勘探中用炸药激发时，一声炮响之后会产生各种各样的地震波。按波在传播过程中质点震动的方向来区分，可以纵波和横波；根据波动所能传播的空间范围而言，地震波又可以分为体波和面波；按照波在传播过程中的传播路径的特点，又可以把地震波分为直达波、反射波、透射波、折射波，等等。地震勘探在石油勘探中除了能产生来自地层界面有用的反射波外，还会产生各种各样的干扰波。因此，我们要更好的了解各种波的产生、特点、用途，等等。下面简单介绍几种地震勘探中产生的地震波。 正文 一、反射波 （一）反射波的形成 1、几何地震学的观点 当炸药在井中爆炸激发地震波时，在雷管引爆几百微妙之内爆炸便完成了，在接近爆炸点的压强是一个延续时间很短的尖脉冲，爆炸脉冲向外传播，压强逐渐减少，地层开始产生弹性形变，形成地震波。地震波继续传播，由于介质对高频的吸收，地震波信号减小。当波入射到两种介质的分界面时（当上层介质波阻抗与下层介质波阻抗不等时，弹性地震波才会发生反射；上层介质波阻抗与下层介质波阻抗差别越大，反射波越强——反射波条件），一部分波回到第一种介质中，这就是所谓的反射波。如图所示 2、物理地震学观点 地震波从震源出发以球面波的方式向下传播，到达反射界面S，S可以就看成有许多
第二章地震波运动学理论 一、名词解释 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.纵测线和非纵测线：激发点与接收点在同一条直线上，这样的测线称为纵测线。用纵测线进行观测得到的时距曲线称为纵时距曲线。激发点不在测线上，用非纵测线进行观测得到的时距曲线称为非纵时距曲线。
第一节河道砂体形态研究 河道砂是河流相储层中最主要和最重要的油气富集场所，因此河道沉积的研究比较系统和深入，国内外相关的文献和研究成果十分丰富，针对研究目标和对象有现代河流沉积、露头剖面砂体、石油钻井三种不同研究体系。因研究的对象和目标的差异，对河道的认识有着不同的的侧重，因侧重面的不同导致对河道砂体几何形态的认识上有很大的差异。利用现代河流沉积和露头沉积研究成果去认识古代河道容易导致认识的偏差，人们对古代河道的认识停留在用现代河流臆测古代河流的程度上；对聚集油气的河道砂体的认识则受勘探阶段和钻井密度的限制，在不同的井网密度阶段下，对河道砂体的认识程度存在很大的不同，因此，我们利用与汪家屯气田河流沉积有关的文献和研究成果，分析河道砂体沉积环境、形成过程、保存条件，总结砂体的四维时空规律，充分认识汪家屯气田河道砂体在空间上可能的展布规律和存在的三维几何形态，通过地震属性成像技术和属性数字地质统计相关分析，从几何形态模型和数字地质方面，识别河流沉积保留下的砂体。 1.国外研究现状 河流分类Leopold(1957)和Wolman（1957）最初将河道体系划分为顺直河、辫状河和曲流河，已为沉积学者所通用。Schumm(1968)根据搬运方式又提出一种分类方案，将河道分为三种类型：即悬载河、混载河、底载河，对现代河流较适用，对古代河流意义不大。B.R.Rust（1978）对河型分类，分为辫状河、曲流河、网状河、顺直河，这一方案在石油业得到广泛认可。 早在50年代Leopold(1957)和Wolman（1957）就建立了河道宽度和曲流带长度的关系以及曲流半径和曲流长度的关系，Strokes(1961)测量了Mesaverd 组河道砂体大小，厚0.6-30.5m，宽1.5-61.0m，长4.5-12km。 Knutson(1971,1976)研究Colorado西部河道砂得出曲流河的宽厚比为14：1。Campbell研究新墨西哥Morrison组Westwater Canyon 砂岩段，总结低弯度辫状河，单个河道平均宽度183m，厚度4m，厚宽比46：1。 Cowan(1991)研究表明，如果砂岩厚度大于12m，就是由几个河道砂体复合而成，指出辫状河道体系平均宽度500m，厚度7m，厚宽比70：1。在曲流河点坝中一个侧积体的最大宽度大致是2/3河宽。 以上河道砂学者的研究现代沉积、露头，还是停留在对河道的静态平面二维或者是剖面二维认识认识上，虽然对于河道的认识有指导作用，但对储集油气的河道砂来说研究还很不够。 只有对砂体的三维空间几何形态认识的基础上，我们才能够识别和寻找河道砂体。 2．大庆油田钻井研究砂体平面形态 钻井研究河道几何形态需要有足够的井网密度，钻井揭示河道是一维的，多口井连的剖面是二维，无论研究的如何精细也不能给出河道在空间的展布几何形态，只有在满足研究需要的足够大的一定区域和密集井网的测井、录井、取心充足的资料背景下，通过特定的砂体对比和组合手段，才可能完整认识保留下的河流沉积砂体类型、微相、三维形态。我们有幸得到了大庆萨中开发区高密度的井网砂体资料。 大庆油田萨中开发区为萨尔图油田中部，面积116km2，萨、葡、高油层高密度的井网为认识砂体提供了很好的例证。 高、葡、萨油层是在青山口组-姚家组-嫩江组早期，松辽盆地由水退转水进
时程分析中地震波输入位置的讨论 摘要:时程分析法通过直接动力分析可得到结构相应随时间的变化关系，能真实地反应结构地震相应随时间变化的全过程，是抗震分析的一种重要方法[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 *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 ………………
V ol 35No.1 Feb.2015 噪 声与振动控制NOISE AND VIBRATION CONTROL 第35卷第1期2015年2月 文章编号：1006-1355(2015)01-0240-04 地震监测中异常次声波的识别方法 左明成，武云 (中国地质大学（武汉）计算机学院，武汉430074) 摘要：地震异常次声波的监测是地震监测中的重要手段和途径。但是，该次声波在接收过程中受到了众多噪声的干扰。为了找到一种地震异常次声波识别监测的有效方法，根据已经收集到的次声波数据进行了分析与研究，按照去噪、特征抽取、信号筛选、分类决策的过程鉴别异常的次声波。在实验中此方法和思路不仅仅得到了较好的识别效果，而且在监测过程中也可识别出矿山爆炸信号和巴东地震次声波信号。从而说明该方法是地震次声波自动识别与监测的一条有效途径。由此，不仅可以较大地减轻地震监测的工作量，实现异常次声波的自动监测，还可以应用到其他的地质灾害的监测和地震的震前监测程序之中。 关键词：声学；信号分析；地震监测；地震次声波中图分类号：TB132；TN911.6；TP18 文献标识码：A DOI 编码：：10.3969/j.issn.1006-1335.2015.01.049 Recognition Method of Abnormal Infrasound in Earthquake Monitoring ZUO Ming-cheng ,WU Yun (China University of Geosciences (Wuhan)Computer College,Wuhan 430074,China ) Abstract :Monitoring the earthquake abnormal infrasound is the important means for seismic monitoring,but the infra-sound is often disturbed by many noises in the receiving process.In order to find an effective method for abnormal sound recognition in seismic monitoring,the collected infrasound wave data was analyzed.With the process of de-noising,feature extraction,signal filtering,classification and determinasion,the abnormal infrasound was detected.Application of this idea and method can get a better recognition effect in the experiment.As an example,this method was applied to identify the mine explosion signal and Badong earthquake infrasound signal.It shows that this method is effective for automatic recogni-tion and monitoring of earthquake infrasound.This method can realize automatic monitoring of abnormal infrasound effec-tively and economically,and can be applied to the monitoring programming for earthquake monitoring and some other geo-logical disasters monitoring. Key words :acoustics ;signal analysis ;seismic monitoring ;seismic infrasound waves 地震发生时震源会向大气中辐射有明显特征的异常次声波[1–3]，这为异常次声波的识别分离提供了先决条件。通过对地震异常次声波的特征进行相关研究，就可以掌握异常次声波的大体形态特征。而 收稿日期：2014－08－01基金项目：国家级大学生创新创业训练项目（201310491060）作者简介：左明成(1992－)，男，山东莱阳人，本科生，空间信 息与数字技术专业，主要研究方向为地震监测、数据挖掘；数字图像处理、三维可视化。E-mail:1317085693@https://www.doczj.com/doc/e615761813.html, 通讯作者：武云，讲师，计算机应用系研究生。 E-mail:23753648@https://www.doczj.com/doc/e615761813.html, 地震在发生之前通常也会产生异常次声波，对于地 震前的监测和其它地质灾害的监测而言，异常次声波的特征研究就显得十分必要了。 1数据源 为了得到次声波数据，在湖北省境内安装次声波接收仪器的方式接收次声波，次声波接收仪使用的是中科院声学所研制的In SYS 2008型号的次声波传感器，仪器可以长时间稳定地接受信号。将16台次声波接收仪接收到的数据通过仪器转换存储在二进制文件中，可以更直观的对次声波数据进行研究。图1中是次声波传感器。
第31卷 第2期 O IL &GAS GEO LOGY 2010年4月 收稿日期:2009-03-04。 第一作者简介:黄军斌(1970 ),男,博士生,能源地质工程。 文章编号:0253-9985(2010)02-0225-07 辫状河流相薄砂体地震子波效应和识别方法 黄军斌1 ,高利军2 ,高 勇3 ,王德发 1 (1.中国地质大学,北京100083; 2.哈里伯顿公司,北京100004;3.中国石油天然气股份有限公司大港油田公司石油研究院,天津300280) 摘要:在地震资料频率一定的情况下,地震子波旁瓣效应会增加薄砂岩识别和解释的难度。为了研究子波旁瓣的影响程度,设计系列正演模型,利用不同频率的零相位雷克子波合成人工地震剖面,确定砂体位置、频率与地震子波旁瓣效应的关系,研究振幅类属性与薄砂体描述的关系。结果表明,综合利用多频带地震剖面和振幅类属性,可以可靠地解释和描述厚度大于1/8波长的辫状河流相薄砂岩储层。 关键词:雷克子波;振幅属性;子波旁瓣;旁瓣效应;薄砂体中图分类号:TE122.2 文献标识码:A W avelet effect and recognition m ethods of thi n sand bodies i n brai ded ri ver facies H uang Junbin 1 ,Gao L ij u n 2 ,Gao Y ong 3 ,and W ang Defa 1 (1.Chi na Universit y of Geoscie nces ,Beiji ng 100083,China ; 2.H alli burton D igit al an d Consultin g Soluti ons ,B eiji ng 100004,China ; 3.Exp loration and D evel opment Rese arch Instit ute ,PetroChi na D a gan g O il f i el d Co mpany,T ia nji n 300280,China) Abst ract :Under a fi x ed seis m ic frequency ,w avelet si d e lobesm ay i n crease t h e d ifficulty o f inter preting and recog n izi n g thin sand bod i e s .To i n vestigate the i m pacts ofw avelet si d e lobe effects ,w e design a series o f for w ard m ode ls and pr oduce synthetic se is m ic profiles using zero phase R icker w avelets o f d ifferent frequencies .The pur pose is to i d entify t h e re lationship of sand body l o cati o n and frequency w ith side lobe effects .This study a lso ai m s at explo ring t h e relati o nsh i p bet w een a m plitude attri b utes and t h e descri p tion of th i n sand bod ies .The results sho w tha,t an i n tegrated use of m ulti frequency profile and a m plitude attributes can re liab l y i n terpret and describe thin sand bod i e s of braided ri v er facies which are thicker t h an 1/8w aveleng t h . K ey w ords :R icker w avele,t a m plit u de a ttri b ute ,si d e lobe of w avele,t si d e lobe effec,t th i n sand body 随着勘探开发的深入,越来越多的难题摆在了勘探者面前,薄砂体描述就是其中之一。如何解决此类问题,一直是值得探讨的课题。很多勘探工作者对此进行了研究,提出接近或大于分辨率(1/4波长)的砂体,可以通过反射特征和属性 参数来解释和描述[1~3] ;而小于分辨率的砂体,集中讨论砂体自身顶底面调谐效应的影响,指出调谐振幅可以用来预测砂体厚度[3~5] 。这些成果适 用于砂体顶面可清晰识别的情况[6~13] ,而当砂体顶面与其他反射相互影响而不易识别时,又如何来进行薄砂体描述呢?毫无疑问,提高地震分辨 率,让薄砂体顶面反射变得清晰,是根本的解决之道。但是,由于技术、资金等条件限制,很难在短时期内从采集和处理方面入手而提高分辨率。而在解释阶段,又需要识别薄砂体,并且描述其分布范围和空间展布特征。在这种情况下,如何立足当前分辨率有限的地震资料,确定影响薄砂体描述的因素,找出一条可行之路呢? 近几年,在研究中国西部某盆地一些气田时,发现辫状河流相薄砂岩的识别和解释与地震子波关系甚大,其中子波旁瓣是影响砂体描述的关键因素,在很大程度上影响着储层的识别和解释。
对于地震波的输入，可以把荷载记录作成文件，利用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
菲涅尔体和透射波 摘要 在地震成像实验中，通常使用基于波动方程高频渐进解的几何射线理论，因此，通常假设地震波沿着空间中一条连接激发点和接受点的无限窄的线传播，称为射线。事实上，地震记录有非常多的频率成分。地震波频率的带限性就表明波的传播应该扩展到几何射线周围的有限空间。这一空间范围就成为菲涅尔体。在这片教案中，我们讲介绍关于菲涅尔体的物理理论，展示适用于带限地震波的波动方程的解。波动方程的有限频理论通过敏感核函数精确地描述了带限透射波和反射波的旅行时与振幅和地球介质中慢度扰动之间的线性关系。菲涅尔体和有限频敏感核函数可以通过地震波相长干涉的概念联系起来。波动方程的有限频理论引出了一个反直觉的结论-在三维几何射线上的点状速度扰动不会不会造成波长的相位扰动。因此，这说明在射线理论下的菲涅尔体理论是波动方程有限频理论在有限频下的一个特例。最后，我们还澄清了关于菲涅尔体宽度限制成像实验分辨率的误解。 引言 在地震成像技术中，射线理论通常在正演和反演中被用有构建正反演波长算子。射线理论之所以收到欢迎部分是由于计算机速度和内存的限制，因为射线理论具有较高的计算效率并且对于各种地震成像方法的应用也比较容易。而另一方面，地震成像实验清晰的表明，射线理论，由于他对波场传播的近似描述，对于散射效应严重的波场的成像是不完备的。Cerveny 给出了对于地震波射线理论的一个全面的理解。 在地震成像实验中，记录到的透射波和反射波信号都是由一个主要由低频信号组成的宽带震源激发产生的，因为地震波的高频信号在地层中很容易衰减。但是射线理论是基于高频近似的，这表明基于射线理论的成像技术和和测量波场这件之能会存在方法上的冲突。这个围绕射线且对带限地震波的传播起主要影响的空间范围就被叫做菲涅尔体。射线理论在地下构造尺度大于记录波场的第一菲涅尔带的介质中能够取得较好的效果。对于低频反射波（频率成分在10-70Hz之间）和透射波（频率成分在300-800Hz之间），第一菲涅尔体的宽度可以分别达到500m和50m的量级。这个宽度要大于我们在陆地和海洋的反射波地震勘探以及井间和垂直地震剖面中想要成像的地下地质特征。 在这篇教案中，我们将看到如何将地震分辨率扩展到识别体积小于第一菲涅尔带的不均匀体。我们将展示如把射线理论下的旅行时和振幅公式扩展到更精确的、可以应用与带限反射和透射地震信号波场近似理论。波动方程的有限频理论提出了反射和透射地震波的敏感核函数（也称作Frechet核函数）。这些有限频Frechet核函数将速度扰动和旅行时与振幅的扰动线性的联系起来。有限频波长近似被直接应用到各种地震成
利用PEER网站选取地震波的方法 云南省土木建筑学会建筑结构专业委员会 PEER是Pacific Earthquake Engineering Research Center（太平洋地震工程研究中心）的简称，设立在美国的加州大学伯克利分校（University of California, Berkeley）。由其运作的PEER Ground Motion Database（PEER地震动数据库）提供了大量的世界各地的地震记录，可自由下载。该网站提供了丰富的查询手段，可按距离、场地、震源类型等条件选择地震记录，也提供了按目标反应谱选择的手段。 该网站是https://www.doczj.com/doc/e615761813.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”。
地震波的频率和振幅 时间：2010-06-05 20:18 来源：unknown 作者：wowglad 点击：$ 7 次2008年12月19日地震波的频率和振幅 1、地震波的频谱及其分析 频谱：谐和振动的振幅和初相位则随频率的改变而改变的关系 ,统称为地震波的频 频谱分 2008年12月19日 地震波的频率和振幅 1、地震波的频谱及其分析 频谱：谐和振动的振幅和初相位则随频率的改变而改变的关系，统称为地震波的频谱。 频谱分为： 振幅谱：振幅随频率变化的关系称为振幅谱。 相位谱：初相位随频率的变化关系称为相位谱。 作用：频率分析，根据有效波和干扰波的频段差异 ①指导野外工作方法的选择 ②给数字滤波和资料等工作提供依据。 频谱分析的方法： 为了研究地震波的频谱特征，可用傅立叶变换把波形函数a(t )变换到频率域中，得到振幅随频率的变化函数A(f),这个变换过程称之为频谱分析方法。 假设波形函数a(t)
--傅氏正变换 ------------- ()-- --傅氏反变换 这两式是等价的,即A⑴与a(t)是一一对应的①S脉冲函数A S (t) ②函数： ③函数： 可以看出：不同时间函数具有不同的频谱。
图、地震波的频率特征 地震波是人工激发的振动，具有连续的频谱，如图所示 图主频fO :振幅谱曲线极大值所对应的频率。 频带的宽度：若|A（f）|最大值为1则可找|A（f）|=的两个频率fl和f2,两者之差△ f=f2 -f1为频带宽度。 大量的实际观测和分析，各种不同类型的地震波的能量主要分布频带是不同的。如图所示。
3、地震波的振幅及其衰减规律 影响地震波激发和接收时振幅和波形的因素: ①激发条件。 ②地震波在传播过程中受到影响。 ③接收条件的影响。 ④其它如地下岩层界面的形态和平滑状态。 图大地低通滤波器效应: 地震波在传播过程中随着距离（或深度）的增加，高频成分会很快地损失，而且波的振幅按指数规律衰减。实际地层对波的这种改造，称之为大地低通滤器效应。 ⑴波前扩散 球面扩散：在均匀介质中，点震源的波前为球面，随着传播距离的增大，球面逐渐扩展，但总能量仍保持不变，而使单位面积上的能量减小，振动的振幅将随之减小，这称之为球面扩散（或波前扩散）。 设某一时刻球面的波前面S,总能量为E,单位面积上的能量为e,则有
地震波的选取方法 建筑抗震设计规范(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后的)，然后将一组人