学位论文近场强震动加速度记录的校正处理方法.doc

《学位论文近场强震动加速度记录的校正处理方法.doc》由会员分享，可在线阅读，更多相关《学位论文近场强震动加速度记录的校正处理方法.doc（72页珍藏版）》请在三一办公上搜索。

1、图书分类号：P315 工学硕士学位论文近场强震动加速度记录的校正处理方法江汶乡中国地震局工程力学研究所中国哈尔滨二零一零年五月 Masters Degree ThesisCorrection method for digital strong motion acceleration records in near-fieldMaster Candidate:Jiang WenxiangTutor:Prof. Yu HaiyingSubject :Civil EngineeringSpecialty:Disaster Prevention and Reduction Engineering an

2、d Protective EngineeringResearch Field:Processing of Strong Motion RecordsInstitute of Engineering Mechanics，China Earthquake AdministrationHarbinChinaMay, 2010图书分类号：P315 单位编码：85406 学 号：2007235工学硕士学位论文近场强震动加速度记录的校正处理方法硕 士 研 究 生： 江汶乡导 师： 于海英研究员申请 学位 级别： 工学硕士学位学 科： 土木工程专 业： 防灾减灾工程及防护工程研 究 方 向： 强震记录处理中

3、国地震局工程力学研究所中国哈尔滨二零一零年五月独 创 性 声 明本人郑重声明：所呈交的学位论文，是本人在导师的指导下，独立进行研究所取得的成果。除文中已经注明引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写过的科研成果，也不包含为获得中国地震局工程力学研究所或其它教育机构的学位或证书而使用过的材料。对本文的研究作出重要贡献的个人和集体，均已在文中以明确方式标明并表示谢意。本人完全意识到本声明的法律责任由本人承担。论文作者签名： 签字日期： 学位论文版权使用授权书本人完全了解中国地震局工程力学研究所有关保留、使用学位论文的规定，同意学校保留或向国家有关部门或机构送交论文的复印件和电子版

4、，允许论文被查阅和借阅；本人授权中国地震局工程力学研究所可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或其他复制手段保存论文和汇编本学位论文，允许被查阅和借阅。(保密的学位论文在解密后适用本授权书)论文作者签名： 导师签名： 签字日期： 签字日期： 摘 要近一个世纪以来，破坏性大地震给人们的生产和生活带来了巨大的损失。但与此同时人类也积累了大量宝贵的强震记录。随着人们对地震研究的深入，近场强震动加速度记录也显得越来越重要。这也对加速度强震记录的精准性提出了更高的要求。由于近断层场地的震动最为激烈，往往会造成仪器墩不同程度的倾斜，而且倾斜的角度和方向是时变的。地震时近

5、断层场地还会有不同程度的旋转，在得到的强震记录中也会耦合旋转分量。这都给强震记录的校正造成了相当大的困难。国内外学者对近场记录的校正方法进行了大量的研究，得出了很多有用的结论。但是在汶川地震中得到的近断层强震记录体现出了更为复杂的误差引入因素。因此，本文在前人的研究基础上，以汶川地震和集集地震的加速度记录为主要研究对象，对近场强震动记录中由仪器倾斜引起的误差校正方法进行了研究。并且在MATLAB环境下实现了校正程序。具体内容如下：1. 本文对数字强震记录系统进行了误差分析；对近场强震动加速度记录的校正处理方法进行了系统的研究；并且对加速度记录的校正准则进行了论述。给出了适合近场强震动加速度记录

6、的校正方法和准则。2. 对5种型号数字强震仪系统进行振动台对比实验，并且对结果进行了分析。给出了加速度记录积分后速度和位移时程零线漂移的原因。并对振动台实验中强震仪记录到的加速度经过校正后再两次积分得出位移时程与实验时记录到的绝对位移进行比较。3. 利用MATLAB语言编写了加速度记录零线校正软件，对2008年5.12汶川8.0级大地震现场加速度记录两次积分后得出永久位移。对汶川强震记录校正后的位移与GPS同震位移进行了比较分析。对台湾集集地震的近场记录用本文方法进行校正，校正结果符合提出的校正准则。结果表明，该方法对大地震近场仪器墩会发生倾斜或产生永久位移的加速度记录的零线有明显的校正效果。

7、关键词： 数字强震仪；近场加速度记录；零线校正；永久位移；MATLABAbstractIn the last century, devastating earthquakes have brought about great lose to peoples lives and property. But at the same time, human have collected huge amount of Strong-Motion records. As research goes on, near-field permanent displacement became more and

8、 more noticed all around the earthquake engineering circles and the records with high accuracy are required. But it is the near field shakes most intensively, and this will force instrument platform to tilt. Moreover, the tilting degree and direction are time-varying. Even worse, near fields are lik

9、ely to spin in great earthquakes. And, somehow, the according records may include spin components. All of these bring a lot of difficult to records correction. Scholars have managed to find correction methods, and have drawn many useful conclusions. But the records acquired form Wenchuan earthquake

10、show more complicate records errors. Therefore, this paper research the correction method based on previous study, and compiled a program in MATLAB environment. Research specifics are as follow:1. Analyzed the error in digital strong motion seismograph, and studied the correction method of near fiel

11、d strong motion acceleration records.2. Revealed the reasons why the baselines usually drift in velocity and displacement tracks by conducted a shake-table experiment and compared the displacement integrated form acceleration and displacement measured directly.3. Compiled a correction software using

12、 MATALB, and acquired permanent displacement by integrated the corrected acceleration records from Wenchuan earthquake. And compared the displacement with GPS coseismic displacement. The result denoted that this method is effective against errors caused by instrument tilt. Key words: digital acceler

13、ograph, near-field acceleration records, base-line correction, permanent displacements, MATLAB目 录 摘 要IAbstractII目 录III第一章 绪 论11.1 引言11.2 国内外研究历程及现状11.3 论文立题背景和本文工作31.4 各章主要研究内容及安排4第二章 加速度记录的误差分析52.1 引言52.2 基线初始值52.3 低频误差52.4 随机噪声72.5 仪器倾斜92.6 两种基线误差的表现形式102.7 小结11第三章 现有的处理方法133.1引言133.2 加速度记录的零线调整（R

14、AP）133.3 USGS使用的基本加速度记录处理方法153.4 Iwan等的零线校正方法（Iwan方法）163.5 Graizer的多项式拟合法173.6 IEM使用的基本加速度记录处理方法173.7 小结19第四章 本文采用的校正方法214.1 引言214.2 零线校正准则214.3 近断层零线校正方法224.4 小结25第五章 振动台对比实验275.1 实验目的275.2 实验用数字强震仪275.3 振动台输入地震波295.4 标准位移计和记录系统305.5 实验步骤305.6 实验数据回收30第六章 程序实现356.1 数值积分方法356.2 主要代码366.2.1 计算斜率和倾斜时刻

15、T的函数366.2.2 校正函数366.2.3 计算傅里叶谱376.3 应用程序界面37第七章 校正结果分析437.1 振动台结果分析437.2 台湾集集地震记录校正效果487.3 汶川地震记录的校正效果517.4 小结60第八章 总结与展望61参考文献62致 谢63作者简介64第一章 绪 论1.1 引言强震记录对地震工作者有相当重要的意义。大量有价值的强震动记录的获取，可为研究震源机制、地震动衰减规律、场地和活断层对地震动的影响、土与结构的相互作用、典型结构的地震反应特性等提供可靠的基础资料。为我国的地震动参数区划图和各行业抗震设计规范的编制和修订提供依据，并为各种重要工程结构的地震反应分析

16、提供典型的输入地震动时程，从而使建设工程的抗震设防要求和抗震设计更为科学、合理。这将有效地减小未来地震发生时各类结构的破坏程度，减轻地震造成的生命财产损失。在强震记录中，近场强震加速度记录记录的是离断层较近的场地的地面运动加速度，所以它对研究震源机制和地震动衰减规律等有着更为突出的意义。集集地震以前，全世界获得的近断层强震记录很少。随着近断层记录的不断丰富，这些记录给人们提供了重要的途径去研究断层破裂过程，地震动的衰减特性和近断层地面永久位移。本文研究的正是这种近场强震动加速度记录，简称强震记录。但是，由于各种原因，记录到的大地震近断层记录的零线严重漂移。导致零线漂移的原因不仅很多而且很复杂，

17、因而零线校正是一个很困难的问题，到目前为止还没有人提出比较完善的解决方法。人们对地震地面运动的强震记录研究还主要停留在峰值加速度和反应谱的阶段，基本上得不到真实的地面运动速度时程和位移时程。这在很大程度上阻碍了人们对地震近断层地面运动规律和极震区抗震设计的深入研究。尽管大多数学者认为由加速度强震记录得不到真实的地面运动速度和位移，但是我们仍可以期望通过对误差因素和校正方法的研究，可以使校正后的结果尽量接近于真实。1.2 国内外研究历程及现状从30年代取得世界上第一条强震动记录到80年代初，强震动记录的处理和分析方法主要是针对模拟强震动记录的。70年代初期， Trifunac提出了一套具有代表性

18、的强震动记录处理方法，比较合理地解决了长期未能解决的在一定频段内对加速度记录二次积分计算位移的问题。70年代中期，我国学者谢礼立等针对所采用的电流计记录式强震仪(主要仪器为RDZ1-12-66)和直接光记录式强震仪(SMA-1)等的特点，详细地分析了在处理模拟记录过程中的各项误差，给出了消除或抑制这些误差的方法，并形成了一套相应的强震动记录数据处理流程。该流程包括以下步骤：1. 加速度记录的调整，消除仪器记录过程和数字化过程中产生的系统误差。2. 校正加速度记录，即仪器校正和零线校正。进行仪器频率响应失真和记录零线长周期失真的校正，消除误差，恢复失真信号。3. 数值积分，取得速度和位移时程。4

19、. 计算反应谱和傅里叶谱。80年代初，随着数字观测技术的发展，数字磁带强震仪已在我国强震动台网和流动观测中使用，这类仪器主要有：美国Kinemetrics公司生产的PDR-1型、DSA-1型，中国地震局工程力学研究所研制的SCQ-1型。所获取的记录动态范围大，频带宽，有绝对时标，波形完整不丢头。和模拟记录不同，这些数字记录的原始记录是不可视的，必须使用专用的回放系统和计算机，专用的数字回放软件。通过对磁带上数据的识别、传送、校验和转换，才能得到可视的地震加速度图和可用的强震动数据。针对我国强震动台网中使用的数字磁带强震仪记录格式，建立了数据回放和数据处理的软件系统，包括数据校验和恢复、数据解码

20、和数据格式转换等。80年代后期，固态记录式强震仪问世。目前我国强震动台网中使用的固态记录式强震仪主要有：美国Kinemetrics公司生产的K2型和Etna型大动态固态数字强震仪；中国地震局工程力学研究所自行研制生产的固态地震动烈度记录仪GDQJ-I型和GDQJ-II型；港震公司生产的GSMA强震仪；瑞士的MR2002强震仪。数字强震仪具有许多模拟强震仪所无法比拟的优越性能。数字强震仪一般都采用力平衡式加速度计，仪器频带范围通常为0至50赫兹或80赫兹左右，高低频响应特性均优于普通模拟记录式强震仪。一般说来，不作仪器频率响应特性校正，就能满足地震工程感兴趣的频率范围要求。有的早期的数字强震仪频

21、带范围较窄（0至30赫兹左右），这种仪器的记录应该进行仪器校正，但仪器校正对低频误差没有什么影响。可以说，低频误差校正是数字强震仪记录处理中的主要问题。有不少学者对数字强震仪系统的低频误差和校正方法进行了广泛的研究。其中，美国学者Iwan对误差原因进行了初步研究21、10美国USGS的学者Boore18、我国的王国权5运用Iwan等的零线校正方法进行了广泛的讨论，得出了一些有价值的结论。台湾学者Hung-Chie Chiu20在文章中给出了一些集集地震近断层记录校正后的速度时程和位移时程。我国学者周雍年等15对常用力平衡式加速度计与各种型号的强震仪组成的系统进行了振动台实验分析，给出了消除数字

22、强震记录的低频误差校正方法，并应用于当时的国家强震动台网常规数据处理工作中。但这种校正方法应用在有永久位移的台站记录中很不合理，得不到地面运动的永久位移。于海英研究员对国家强震动台网入网的5种型号数字强震仪系统做了振动台对比实验8，分析了该系统加速度记录积分后速度和位移时程零线漂移的原因。提出了加速度记录零线漂移的校正方法和校正准则。为了印证校正方法的可靠性，对振动台实验中强震仪记录到的加速度用提出的校正方法进行校正，然后对校正后的加速度进行两次积分得出位移时程。将校正后的位移时程与实验时记录到的绝对位移时程进行比较。结果表明，校正后的位移时程和振动台绝对位移时程的一致性良好；对2008年5.

23、12汶川8.0级大地震近场校正后的加速度记录两次积分后得到的永久位移与GPS同震位移进行比较，两者相符合，校正效果明显；对1999年台湾9.21集集7.6级地震的近场加速度记录用提出的方法进行校正，校正后的位移末端趋于平稳，完全符合校正准则。所有的结果都表明，提出的校正方法对大地震近场仪器墩会发生倾斜的记录有明显的校正效果，基本上能得到符合实情的永久位移。可以给出加速度积分后的速度和位移并符合校正准则。1.3 论文立题背景和本文工作近断层强震记录是非常宝贵的地震研究基础资料，但是记录中往往会包含各种误差，使得强震加速度记录不能反映真实的震时地面运动。国内外学者对近场记录的校正方法进行了大量的研

24、究，得出了很多有用的结论。在2008年5.12汶川8.0级地震中，国家强震台网中心获得了大批近断层强震动记录，受到国际地震工程界的普遍重视。和其它的大震记录一样，这次地震近断层记录的零线严重漂移。导致零线漂移的原因不仅很多而且又很复杂，因而零线校正是一个很困难的问题，到目前为止还没有人能提出比较完善的解决方法。这在很大程度上阻碍了人们对这次地震近断层地面运动规律的研究。本文对以往零线校正方法进行了梳理和研究。在前人的研究基础上，以中国地震局公益性行业科研专项（200708010）“强震动记录误差分析方法与校正处理技术研究”为支撑，利用国家强震动台网入网的5种型号数字强震仪系统振动台对比实验结果

25、8，分析了该系统加速度记录积分后速度和位移时程零线漂移的原因。在于海英老师的指导下，以汶川地震和集集地震为主要对象，对近场强震动记录的校正方法进行了研究，提出了加速度记录的零线漂移校正方法和校正准则。为了印证本文校正方法的可靠性，利用Matlab语言编写了校正处理软件，将校正后的位移时程与振动台实验记录的绝对位移时程进行比较，结果表明，校正位移时程和振动台绝对位移时程一致性较好。对2008年5.12汶川8.0级大地震近场加速度记录两次积分后得到的永久位移与GPS同震位移进行比较，校正效果明显。对1999年台湾集集7.6级地震的近场加速度记录用提出的方法进行校正，校正后的位移末端趋于平稳，完全符

26、合校正准则。所有的结果都表明，提出的校正方法对大地震近场仪器墩会发生倾斜的记录有明显的校正效果，基本上能得到符合实情的永久位移。目前，本文编写的加速度记录校正处理软件已经应用于国家强震动台网中心日常处理工作中。1.4 各章主要研究内容及安排本文共分为八章，第一章是绪论，简单介绍强震记录和基线误差及研究历程。第二章分析加速度记录的误差。第三章介绍现有的处理方法。第四章介绍本文采用的方法。第五章介绍振动台实验。第六章为程序实现。第七章为校正结果分析。第八章是总结与展望。第二章 加速度记录的误差分析2.1 引言数字强震记录中同时存在高频误差和低频误差，但只有低频误差影响零线。低频误差的主要来源有初始

27、值、仪器噪声、背景噪声以及处理误差，与低频噪声有关的误差可通过高通滤波的方法在一定程度上得到抑制。但是在近场强震记录中存在永久位移信息，滤波以后会丢失这个位移信息。所以滤波的方法不适合于近场强震记录。对汶川8.0级地震而言，靠近断层的区域发生了强烈的地表变形(地表破裂、竖直抬升、水平位移等)，可能导致强震仪发生倾斜，强震仪微小的倾斜将导致零线产生很大的漂移。因此将导致零线漂移的原因归为两类:低频误差和强震仪的倾斜。2.2 基线初始值加速度记录中的事前部分记录的是地震波到达强震仪所在位置之前的地面的运动状态，事前部分的地面运动加速度应当接近零。但实际上由于各种原因使得记录事前部分的加速度值并不为

28、零，而是在一个水平线附近震动。这个水平线的值就是基线初始值。在所有的基线校正方法中，都应首先从原始记录全时程中减去事前部分的平均值。这样就使得事前部分的地面运动加速度接近零。这一操作只会使得加速度整体发生上下平移，没有改变加速度的相对大小，却使得加速度值在地震波到达之前停留在零值附近。因此这种方法也称作加速度时程的基线初始化，下文简称RAP(Removing Average of Pre-event)，也叫零线调整。这是对任何一条原始记录首先应该进行的一个处理。2.3 低频误差数字强震记录中的低频误差包含了仪器噪声和场地背景噪声两类。仪器噪声主要有电子噪声、可能因采样率和分辨率不足引起的误差、

29、传感器材料和电路的微小磁滞效应以及其它不明原因引起的噪声。现代数字强震仪的低频响应特性良好，一般不需要对仪器的低频响应进行校正。数字强震仪的采样率通常为200sps，对低频信号来说，这样高的采样率不会引起多大的低频误差。分辨率不足引起的低频误差也是很有限的。电子噪声的水平与仪器本身有关外，也受环境因素的影响。电子噪声会引起仪器记录的零基线产生微小的变动。场地背景噪声来自各种环境振动，包括海洋波浪、风、各种人为活动等。场地背景噪声的来源和频率范围很广，而且每个噪音来源对总体背景噪音的贡献很小，可以认为是服从高斯分布的白噪声。场地所处环境对背景噪声的幅值有很大的影响。邱宏智研究员4对数字化强震仪的

30、记录和模拟式强震仪的记录进行过对比研究。指出：尽管相对于模拟式记录，数字化强震记录的质量有明显的提高，但绝大多数导致模拟式记录产生低频误差的因素在数字化记录中仍然会出现。零线漂移源于仪器的记录过程中和其后对资料的处理过程中，前者称为系统误差，后者称为偶然误差。系统误差包括源于纸张或胶片的弯曲或变形、背景噪声、仪器本身产生的噪声、初始零线的漂移以及零线的不稳定等；偶然误差主要是指数字化过程中产生的误差。对于记录汶川8.0级地震的数字化强震仪来讲，则不存在由于纸张变形等产生的误差；由于数字强震仪能够完整地记录长达20s事前部分，初始零线可以说已经确定。因此，由于初始零线不确定而造成的误差也很小。对

31、于数字化强震记录来讲，在数据处理过程中产生的处理误差也几乎不存在。因此，数字化记录零线漂移主要源于记录过程中产生的系统误差。由于记录过程中产生的误差主要影响低频信息，因此，也称作低频误差。仪器误差主要来源于5个方面:仪器的不准确响应(imperfect instrument response)、仪器分辨率有限、采样率不足、电磁噪声、传感器的物质疲劳(material fatigue of the sensor).仪器的不准确响应主要发生在高频段，这种差对高频信号影响当大的，但对于低频信号影响很小。而使零线漂移的主要是低频信号，因此传感器的不准确反应对零线的影响很小。产生仪器误差的另一个原因是仪

32、器的分辨率不足。由于仪器在数字化的过程中，小于一个数字化单位的数字值被忽略掉，这种情况产生的加速度误差经过2次积分后有可能对位移产生影响。但是由于分辨率不足产生的零线漂移的上限是半个数字化单位，而且产生正负误差的概率相等，由于正负误差的抵消，实际上对积分后得到的速度和位移的影响很小。国家强震动台网使用的5种型号强震仪（Etna/ES-T、GDQJ-II/SLJ-100、Geosig/SLJ-100、GSMA-2400/SLJ-100、MR2002/SLJ-100型）的记录系统的模数转换是24位。数字化单位已经足够小，产生的误差不会造成零线的漂移。对于大震近断层强震记录来讲，真实的位移是相当大的

33、，有时能达到2米，因而仪器分辨率不足造成的误差就更显得微不足道。由于采样率的不足，强震记录会在相邻的两个采样点之间丢失一些高频信号，但对低频信号影响不大。对于上述数字化强震仪，采样率为200sps，个别为100sps。采样率不足对低频地震信号产生的误差是很小的。电磁噪声的大小决定于仪器本身和其周围的环境，其对环境温度的变化很敏感，从波形上来讲是随机的，其频率分布在一个比较宽的频带范围内。由电磁噪声和场地背景噪声产生的误差是很难分开的，对于它们可以通过滤波的办法在一定程度上来减小。传感器物质疲劳会导致零线微小的、持续的漂移，会对零线产生比较大的影响，Iwan21等人在零线校正中考虑的磁滞效应主要

34、就是由于传感器的物质疲劳所产生的。2.4 随机噪声在一般的场地下，背景噪声的来源很多，每个影响因素都是随机的，且每个因素对总噪音的贡献很小。所以背景噪声的主要特征是频率成分丰富的随机波形(random waveform )，其分布基本服从正态分布。其频率成分在各个频率都基本相等，是一种白噪音。汶川地震卧龙台强震记录事前10秒钟包含的背景噪音分布如图2.1。可见，背景噪音的确是服从0均值的高斯噪音。背景噪音对位移的影响可以作以下分析16：假定数字加速度记录的采样间隔为。记录中包含有标准方差为的随机噪音。那么由此可得随机噪音在 时的误差可以近似为：这个方程是典型的随机游走。积分得到近似位移：联合此

35、两个方程可以得到：我们假设代表位移在记录末尾的值。又上式可知，末端位移是加速度值的加权和。所以是和随机加速度线性相关的随机变量。末端位移的期望为0，因为加速度的平均值接近0。在实际中末端位移接近0。 上式的加权和可以用积分代替得：其中是记录的时长。 图2.1 背景噪音分布图典型的汶川地震卧龙台强震记录事前10秒钟包含的随机噪音如图2.2。计算得到的位移为：-0.0032cm，若在整个记录的160秒时间重复以上随机噪音那么对位移的影响为-0.0512cm.这个位移值小的完全可以忽略。所以，理论分析和实际验算都支持这种观念：随机误差不会对位移产生较大影响，对位移的影响完全可以忽略。图2.2 汶川地

36、震卧龙台事前10秒记录2.5 仪器倾斜由于近断层的地面有不均匀变形，使得近断层的强震观测仪很有可能发生倾斜，仪器的微小倾斜都将导致强震记录零线的严重漂移。震后通过观察传感器上的水准泡得知部分近断层台站的仪器确实发生了一定程度的倾斜8。即使仪器发生很微小的倾斜，也将在位移时程中产生严重的零线漂移18。典型的强震仪器力平衡式加速度水平向传感器原理如图2.3，其中质量块受到约束，只能作沿弹簧轴线方向的运动。图2.3 典型的强震仪器力平衡式加速度水平向传感器原理传感器发角度的倾斜后，机械部分的侧面如图2.4。这时重力会有个沿弹簧方向的分量，大小为M*g*sin()。其中M为质量块的质量，g为重力加速度

37、。此时质量块沿弹簧轴线的受力状态等效于传感器没有倾斜时在水平方向上受到大小为M*g*sin()的力。也就等效于传感器水平状态时，有个g*sin()大小的水平加速度，那么加速度强震仪就会记录到这个加速度。但实际上仪器只是倾斜了，场地并没有加速度，所以仪器的倾斜会在加速度记录中引起大小为g*sin()的误差。以上分析的是一个记录方向上的倾斜造成的误差，其它2个方向也是同样道理。只要是利用了质量块惯性的强震加速度传感器，无论是惯性加速度传感器，还是力平衡式加速度传感器，以上分析都适用。只是力平衡式传感器利用了反馈技术，使得即便场地有较大加速度时，质量块都保持在初始位置附近。这样，非线性大大降低了，灵

38、敏度大幅提高，并有效地改善由于温度等环境影响所造成的误差。图2.4 传感器倾斜后机械部分的侧面图若仪器的倾斜角度为，则该倾斜导致水平方向加速度记录零线的漂移为：ap=g*sin()竖直方向的零线漂移为：av=g*(1-cos()一般情况下仪器的倾斜角很小。根据震后调查，记录到的倾斜角度较大的也只有0.0611rad22。而且sin与为等价无穷小。因此，对于水平方向，零线的漂移与成正比：ap =g*同理，在竖直方向零线的漂移与0.52成正比：av=g*0.52其中，a为记录中的加速度基线的漂移量，g为重力加速度。 假设传感器的边长为10cm，场地的不均匀变形使得传感器的一侧比另一侧高了一张纸的厚

39、度。这样的不均匀变形在强震动阶段是很容易发生的。100张纸厚约1cm，那么，一张纸厚度为0.01cm。这样传感器的倾斜角为0.001rad。那么ap=g*=1gal。加速度零线的微小漂移经过积分后，由于累积效应，会对速度时程产生较大影响，对位移时程将产生严重的影响。假设在强震记录的30秒处开始出现上述不均匀变形，则这种情况持续到90秒记录结束时，在水平方向位移时程中产生的零线的漂移最终可达18m。然而，仪器的微小倾斜对竖直方向零线的影响比水平方向小的多。在上例中，假设同样的不均匀变形持续60秒，在竖直向位移中产生的漂移只有0.9cm。因此从理论上来讲，仪器的微小倾斜对竖直方向零线的影响比对水平

40、方向的影响小的多。实际的近场记录中也是如此，近场记录的未校正速度时程和位移时程在水平方法的漂移较大；但是在竖直方向的漂移较小，在速度时程中基本看不到漂移，只有在位移时程中会发现漂移。2.6 两种基线误差的表现形式加速度记录中的误差相对加速度来说很小，是直接观察不到的。但是通过积分得到速度和位移，加速度中的误差就会在速度和位移中有明显的表现。加速度中有可能出现脉冲和台阶，表现形式如图2.5和2.6。从图中可见，脉冲误差在速度时程中表现为台阶，在位移时程中表现为斜线。台阶式误差在速度时程中表现为斜线，在位移时程中表现为抛物线。其中台阶是由仪器倾斜造成的，是近场强震加速度记录中最常见的误差。图2.5

41、 加速度中脉冲误差表现图2.6 加速度中台阶误差表现2.7 小结综合本章分析可知：强震记录中包含多种误差，其中，基线初始值由事前记录确定；记录中的背景噪音基本服从高斯分布，幅值很小，对记录的影响不大；仪器误差包含仪器的不完全响应、分辨率有限、采样率不足、传感器疲劳和电磁噪声，这几种误差来源对加速度记录的影响较小；强震记录中包含处理过程中的误差，但是数字化记录中的处理误差很小；仪器倾斜对加速度记录的影响最大，本文研究的校正方法主要是针对这种误差。第三章 现有的处理方法3.1引言强震加速度记录的校正方法是一个重要的课题。从人类获得第一条强震加速度记录的时候，就开始了对其进行校正的研究。这些研究取得

42、了一些进展，得到了不少有用的结论，非常值得借鉴。本章将分别详细介绍美国地质调查局USGS常用的零线校正方法BAP、加州理工大学Iwan等人的方法以及Graizer的方法和国家强震动台网中心常用的处理方法。3.2 加速度记录的零线调整（RAP）加速度记录中的事前部分记录的是地震波到达强震仪所在位置之前的地面运动状态，所以事前部分的地面运动加速度应当接近零。但实际上由于仪器初始偏置未调零、环境温度变化和传感器簧片变形因素，使得记录事前部分的加速度值并不为零，而是在一个水平线附近震动。这个水平线的值就是基线初始值。例如：什邡八角台51SFB东西方向记录的事前部分平均值为0.2013cm/s/s，而绵

43、竹清平台51MQP东西方向记录的事前部分平均值有0.1943 cm/s/s。可想而知，对这种记录积分得到的速度时程从0时刻开始就会有漂移趋势。直接积分得到速度和位移时程如图3.1和3.2所示。位移时程中什邡八角台高达-4000cm，绵竹清平台高达400cm。而且，什邡八角台的积分速度和位移曲线零线漂移严重，十分夸张，运动末速度不为零，而且积分位移时程曲线末尾不稳定，显然是错误的结果。绵竹清平台积分速度曲线的零线也有轻微漂移，使得积分位移曲线零线严重漂移。这些结果显然与实际情况不符，需要进行校正。在所有强震记录的零线校正方法中，首先应该从整个原始记录中减去事前部分的平均值，这样就使得事前部分的值

44、在零附近。这一步骤称为基线初始，也叫RAP(Removing Average of Pre-event)。这样处理后的记录在国家强震动台网中心也叫未校正加速度(uncorrected Acceleration)。可以肯定：加速度时程的基线初始化是一个必不可少的步骤，但只作这一步所得到的速度和位移时程还是很不理想，不过事前部分已经不会发生漂移。这是对任何加速度零线进行进一步校正的前提。接下来的零线校正方法都是对初始化后的加速度时程进行校正。图3.1 什邡八角台竖直向原始加速度记录和直接积分得到的速度和位移时程图3.2 绵竹清平台东西向原始加速度记录和直接积分得到的速度和位移时程3.3 USGS使

45、用的基本加速度记录处理方法美国地质调查局(USGS)提出了BAP（Basic Strong-Motion Accelerogram Processing）方法19。BAP的基本思路是用一条最佳直线拟合加速度时程零线，再从加速度时程中减去这条直线，然后进行滤波处理。这种方法实际上是认为零线漂移从地震动一开始就产生，并且保持漂移的增长率不变。具体步骤如下：步骤一：从整个加速度时程中减去事前部分(0-20s)的平均值(RAP)。步骤二：用用最小二乘法对加速度时程拟合出一条直线，从加速度时程中减去该直线。步骤三：用Buterworth滤波器进行高通滤波，得到校正后的加速度时程曲线。用滤波频带为0.01

46、-80Hz的BAP方法处理后，汶川主震绵竹清平台东西方向加速度记录积分速度和位移如图3.3示。值得注意的是，大地震时近断层有可能产生永久位移，但经过高通滤波处理后将会丢失地面的永久位移信息。由图3.3可见，位移时程绝大部分地面永久位移信息已经被滤掉，而且位移末尾很不稳定，呈现简谐振荡。说明该方法不适合处理大地震近断层加速度记录。图3.3 汶川主震绵竹清平台东西方向校正加速度、积分速度和位移时程（USGS方法）3.4 Iwan等的零线校正方法（Iwan方法）Iwan等人曾对Kinemetrics公司生产的PDR-1和FBA-13型数字化强震仪的性能作了一些实验研究21。结果，发现在加速度超过一定的界限时其记录的零线会发生跳跃现象。尽管这种跳跃对加速度时程本身的影响十分微小，但对由加速度时程积分得到的速度时程有一定的影响，对经过二次积分得到的位移时程将产生很大的影响。Iwan等人认为，这种现象可能是由于传感器系统机械或电路的微小磁滞效应引起的。他们发现，对于实验的PDR-1和FBA-13型数字强震仪，这种磁滞效应在加速度大于50gal时开始出现，因此加速度的零线漂移在强震阶段和其后的延续阶段应当有所不同，因此他们提出了将加速度记录分两段进行零线校正的方法。I