第３０卷第３期 ２ ０ １ ０年６月 大地测量与地球动力学 ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＤＥＳＹ ＡＮＤ ＧＥＯＤＹＮＡＭＩＣＳ Ｖｏｌ＿３Ｏ Ｎｏ．３ Ｊｕｎｅ，２０１０ 文章编号：１６７１－５９４２【２０１０）０３４３０４７４３４ 强震动加速度记录基线校正问题探讨 郑水明 周宝峰 ’ 温瑞智。’ 王 岚 ， ’ ４３００７１ 、 ／，１）中国地震局地震研究所，武汉ｌ ２）武汉地震工程研究院，武汉４３００７１ ＼３）中国地震局工程力学研究所，哈尔滨 ｌ１５００９０ 摘 要 针对汶川Ｍｓ８．０地震中３个典型的加速度基线漂移特性，分别采用基线初始化法、１ｗａｎ法及低频滤波法 进行了基线校正。通过校正，基本消除了加速度记录中的低频噪音干扰，使基线大致回到了零线。通过实例分析 了３种校正方法存在的问题。 关键词 强震动加速度；低频；ａｎ速度时程；基线漂移；基线校正 中图分类号：Ｐ３１５．６３ 文献标识码：Ａ ＤＩＳＣＵＳＳＩＯＮ ＯＮ ＢＡＳＥＬＩＮＥ ＣＯＲＲＥＣＴＩｏＮ ｏＦ ＳＴＲｏＮＧ ＭｏＴＩｏＮ ＤＡＴＡ Ｚｈｅｎｇ Ｓｈｕｉｍｉｎｇ ＇ ，Ｚｈｏｕ Ｂａｏｆｅｎｇ ，Ｗｅｎ Ｒｕｉｚｈｉ ）ａｎｄ Ｗａｎｇ Ｌａｎ ， ） ｒ１）Ｉｎｓｔｉｔｕｔｅ ｏｆＳｅｉｓｍｏｌｏｇｙ，ＣＥＡ，Ｗｕｈａｎ ４３００７１ ３）Ｉｎｓｔｉｔｕｔｅ ｏｆＥｎｇｉｎｅｅｒｉｎｇ Ｍｅｃｈａｎｉｃｓ，ＣＥＡ，Ｈａｒｂｉｎ １５００８０ 、 ｆ ２）Ｗｕｈａｎ Ｉｎｓｔｉｔｕｔｅ ｏｆＥａｒｔｈｑｕａｋｅ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｗｕｈａｎ ４３００７１ ｌ Ａｂｓｔｒａｃｔ Ｔｈｅ ｓｏｕｒｃｅｓ ｏｆ ｓｏｍｅ ｐｒｉｍａｒｙ ｌｏｗ ｆｒｅｑｕｅｎｃｙ ｎｏｉｓｅ ｅｆｆｅｃｔｉｎｇ ｏｎ ｓｔｒｏｎｇ ｍｏｔｉｏｎ ａｃｃｅｌｅｒａｔｉｏｎ ａｒｅ ａｎａ— ｌｙｚｅｄ．Ｔｏ ｒｅｍｏｖｅ ｔｈｅ ｎｏｉｓｅ，ｔｈｒｅｅ ｓｃｈｅｍｅｓ ｏｆ ｂａｓｅｌｉｎｅ ｃｏｒｒｅｃｔｉｏｎ ａｒｅ ｉｎｔｒｏｄｕｃｅｄ，ａｎｄ ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｔｈｅ ｓｃｈｅｍｅｓ ｏｎ ｔｈｅ ｖｅｌｏｃｉｔｙ ａｎｄ ｄｉｓｐｌａｃｅｍｅｎｔ ｔｉｍｅ ｓｅｒｉｅｓ ｉｎｔｅｇｒａｔｅｄ ｔｈｅ ｐｒｏｃｅｓｓｅｄ ａｃｃｅｌｅｒａｔｉｏｎ ｉｓ ｐｏｉｎｔｅｄ ｏｕｔ．Ａｐｐｌｙ ｔｈｅ ｓｃｈｅｍｅｓ， ｔｈｒｅｅ ｔｙｐｉｃａｌ ａｃｃｅｌｅｒａｔｉｏｎｓ ｏｆ Ｗｅｎｃｈｕａｎ Ｍｓ８．０ ｅａｒｔｈｑｕａｋｅ ａｒｅ ｃｏｒｒｅｃｔｅｄ，ａｎｄ ｓｏｍｅ ｑｕｅｓｔｉｏｎｓ ａｒｅ ａｌｓｏ ｐｕｔ ｆｏｒｗａｒｄ ｆ０ｒ ｔｈｅ ｃｏｎｓｉｄｅｒａｔｉｏｎ ｏｆ ｂａｓｅｌｉｎｅ ｃｏｒｒｅｃｔｉｏｎ． Ｋｅｙ ｗｏｒｄｓ：ｓｔｒｏｎｇ ｍｏｔｉｏｎ ａｃｃｅｌｅｒａｔｉｏｎ；ｌｏｗ ｆｒｅｑｕｅｎｃｙ；ａｃｃｅｌｅｒａｔｉｏｎ ｔｉｍｅ ｓｅｒｉｅｓ；ｂａｓｅｌｉｎｅ ｓｈｉｆｔ；ｂａｓｅｌｉｎｅ ｃｏｒｒｅｃｔｉｏｎ １前言 强震仪记录的不仅是地震时纯粹的地震地面运 动信息，还包含有复杂的噪音，其中的低频噪音会导 致加速度时程基线出现漂移…。基线漂移对加速 度时程本身可能影响很小（一般不超过峰值加速度 值的２％），但通过积分求速度、位移时程时，漂移被 逐步放大，从而对速度、位移时程产生很大的影 基线校正。为此，我们将用基线初始化法、１ｗａｎ法 及低频滤波法对汶川地震３个典型的加速度基线漂 移进行改正，并探讨其改正效果。 ２加速度记录基线漂移的根源及其影 响 对于数字强震仪记录而言，导致加速度记录基 响 。因此，在使用加速度记录时，一般需要进行 收稿日期：２００９．１２－３１ 基金项目：中国地震局地震研究所所长基金（ＩＳ２００７５６０５０） 作者简介：郑水明，男，１９７８年生，硕士，工程师，主要从事地震工程开发工作．Ｅ—ｍａｉｌ：ｚｈｅｎｇｓｈｕｉｍｉｎｇ２０＠ｓｏｈｕ．ｃｏｎ 大地测量与地球动力学 ３０卷 线漂移的根源主要为仪器记录过程中产生的误差， 如传感器的磁滞现象和背景噪声，以及传感器的倾 斜等。对于模拟强震仪记录，还存在仪器分辨率有 限等产生的干扰 。 传感器的磁滞效应主要是由于传感器的物质疲 劳所产生的。１ｗａｎ对美国凯尼公司生产的ＦＢＡ一 １３传感器进行测试时发现，当加速度小于５０ｃｍ／ｓ ， 传感器的输出基线变化很小，但超过这个值后基线 会突然变动，致使积分后的位移时程发生明显的漂 移 。 传感器倾斜主要发生在地震近场区强震观测 台。由于地震时地面的强烈运动，观测场地可能发 生不均匀变形、滑移，致使传感器倾斜，倾斜可能导 致加速度基线严重漂移。相对于垂直向记录而言， 倾斜对水平记录的基线的影响更大＿３＇７＿。 背景噪音的主要特征是频率丰富的随机波形， 与记录场地条件密切相关。背景噪音导致加速度记 录的初始值不为零，从而对加速度基线产生影 响 ， 。 强震仪分辨率不足对基线漂移可能会产生不小 的影响，也会导致波形在一定程度上的变形，其上限 一般不超过分辨率的一半 Ｊ。分辨率不足主要存 在于模拟强震仪中，典型的模拟记录最大分辨率是 ２ ｃｍ／ｓ ，而４－２ｇ量程、１２位数字强震仪的分辨率是 １ ｃｍ／ｓ ，２４位的分辨率则为０．０００２ ｃｍ／ｓ 。 ３基线漂移校正方法 目前，加速度时程基线校正大致有３种方法。 １）基线初始化。从加速度记录中减去震前记 录平均值，对于没有震前记录的加速度记录则减去 整个记录平均值。对于复杂的基线漂移进行处理， 其结果可作为进一步校正处理的基础。 ２）对加速度时程用方程或积分求得的速度时 程进行拟合，然后对应地从加速度时程中减去拟合 方程或拟合方程的微商。方程形式一般为分段线性 方程或二次方程，常用的为１ｗａｎ法及ｖ０法 。 １ｗａｎ法基本思路是将整个加速度时程分为３ 段，中间一段代表由磁滞效应引起的基线偏移发生 在强震部分，并且认为在同一段内基线的偏移是一 个常量，用速度时程基线的斜率来表示加速度基线 的偏移量。依分段方法的不同，１ｗａｎ法分为两种 （图１）。 （１）１ｗａｎ法① １）计算加速度时程震前部分的平均值，整个原 始加速度时程均减去这个值后积分求得速度时程。 ２）选取ｔ 、ｔ 时刻为分界点，将整个速度时程分 为３段，ｔｓ为加速度时程结束的时刻，其中ｔ 、ｔ 分 别为加速度首次和最后一次超过或达到５０ ｃｍ／ｓ 的对应的时刻。然后用直线拟合ｔ：～ｆ，时间段的速 度时程，口，为该直线的斜率，并从相应时段的加速 度时程中减去０，。 ｖｓ（ｔ）＝ ＋ （１） ３）在加速度时程的中问段（ｔ ～ｔ ）减去相应速 度时程拟合直线的斜率ａ 。ａ 确定原则是使加速 度时程基线校正后积分得到的速度时程的末尾速度 为０。当ｔ：点的速度ａ （ｔ 一ｔ ）等于步骤２中拟合 直线在ｔ 点对应的速度值ｖｉ（ｔ：）即可满足要求，即 ｖｉ（ｔ２）＝ａｍ（￡２一ｔ１）。 （２）１ｗａｎ法② 方法①中ｔ 是从对基线偏移的物理机制考虑 而确定的，其实在加速度首次超过５０ ｃｍ／ｓ 的时刻 ｔ 和加速度时程的终了时刻ｔｚ之间，有很多点都能 够象ｔ 一样满足最终速度为零的条件，如果 ２在ｔｌ 和矿之间自由移动，则相应地得到的最终位移在一 个很大的范围内变化。为此１ｗａｎ等人提出了第二 种基线校正的方法，和方法①不同之处在于在ｔ。和 ｔｓ之间选取ｔ ，使积分求得的位移时程的最终位移 最小，满足该条件的ｔ 可通过下式求得，值得指出， 这种选择ｔ，的方法没有理论根据。 （２Ｖｏｔｆ＋ 一ｒｏｔ１）一２Ｄｕ￣ｏ ／，＇、 ＋ 一 ％＋ １ ， Ｄ ＝。， Ｄ一 一０．５ｖｓ（ｔ２）（ｔ２一ｔ１）＋ｖｚ（ｔ２）（ｔｚ—ｔ２）＋ ０．５ａｚ（ｔｅ—ｔ２） （３） 式中，Ｄ一，和Ｄ…分别为由未经过基线校正的加速 度时程二次积分求得的位移时程的最终位移和经过 基线校正的加速度时程对应的位移时程的最终位 移。 Ｂｏｏｒｅ等人 用１ｗａｎ法对集集地震加速度记录 进行基线处理后发现，ｔ 的选择对地面永久位移影 响较大（图２）。因此，Ｂｏｏｒｅ等人对１ｗａｎ法进行了 修改（图１），步骤如下ｉ １）计算记录从零时刻至地震波初达时刻ｔ 时 间段的加速度平均值，整个加速度记录减去平均值。 ２）积分第一步校正后的加速度记录，获得速度 时程。 ３）从ｔ 至记录结尾用二次方程对速度时程进 行拟合，方程为 ＝ｃ （ｔ—ｔ ）＋ｃ （ｔ—ｔ ） 。方程 如此选择是为了满足ｔ＝ｔ 时，初始速度值为０。 ４）第一步校正后的加速度记录从ｔ。至记录结 尾减去第三步拟合方程的微商，ｃ ＋２ｃ （ｔ—ｔｐ）。 对于没有事前触发的加速度记录，这种方法的 初始速度为０，条件无法满足。因此，对拟合方程进 第３期 郑水明等：强震动加速度记录基线校正问题探讨 ４９ 行了修改。对减去整个记录平均值的加速度时程进 行拟合，方程为 ＝ｃ。＋ｃ ｔ＋Ｃ２ｔ ，然后用减去整 时序中加入了足够长的零段。 滤波过程中，加速度信息会被引人新的噪音信 个记录平均值的加速度时程减去方程的微商。 对于记录时间较长的加速度，利用方程拟合进 行基线校正时，拟合方差可能会很大，导致部分的记 录信息严重失真。 图１ １ｗａｎ法（实线）及Ｂｏｏｒ等基线校正（虚线）示意图 Ｆｉｇ．１ Ｓｋｅｔｃｈ ｏｆ ｂａｓｅｌｉｎｅ ｃｏｒｒｅｃｔｉｏｎ ｕｎｄｅｒ １ｗａｎ ｓｃｈｅｍｅ （ｒｅａｌ ｌｉｎｅ）ａｎｄ Ｂｏｏｒ ｓｃｈｅｍｅ（ｂｒｏｋｅｎ ｌｉｎｅ） 蓉 ＼ 图２ １９９９年台湾集集地震ＴＣＵ１２９台ＥＷ向加速度经 不同ｔ：选择的１ｗａｎ法、基线初始化及截止频率为 ０．０２ Ｈｚ低频滤波分别进行基线校正后积分获得 的位移时程 Ｆｉｇ．２ Ｄｉｓｐｌａｃｅｍｅｎｔｓ ｏｂｔａｉｎｅｄ ｂｙ ｄｏｕｂｌｅ ｉｎｔｅｇｒａｔｉｏｎ ｏｆ ｔｈｅ ＥＷ ｃｏｍｐｏｎｅｎｔ ｏｆ ａｃｃｅｌｅｒａｔｉｏｎ ｒｅｃｏｒｄｅｄ ａｔ ＴＣＵ１２９ ａｎｄ ｍｏｄｉｆｉｅｄ ｕｓｉｎｇ ｂａｓｅｌｉｎｅ ｃｏｒｒｅｃｔｉｏｎｓ ｏｆ ｒｅｍｏｖｅ ｍｅａｎ ｏｎｌｙ，ｑｕａｄｒａｔｉｃ ｉｆｔ ｔｏ ｖｅｌｏｃｉｔｙ，ｌｏｗ—ｃｕｔ ｆｉｈ ｏｎ— ｌｙ ＝０．０２ Ｈｚ）ａｎｄ １ｗａｎ ｓｃｈｅｍｅ ｔ２＝３０ Ｓ，５０ Ｓ， ７０ Ｓ ３．３低频滤波【８，９ Ｊ 低频噪音是导致加速度基线漂移主要原因，利 用低频滤波处理，能有效地清除其对加速度基线影 响。Ｔｒｉｆｕｎａｃ 首次提议利用低切滤波对低频信号 导致的加速度时程基线漂移进行校正，但过程很复 杂。Ｃｏｎｖｅｒｓｅ ａｎｄ Ｂｒａｄｙ 对其进行了改进，在记录 号，并且部分长周期信息或地面永久位移信息也可 能会被滤掉 ’ｍ 。 实际工作中，为尽可能地消减记录中的低频噪 音，可以将这３种方法相互结合进行基线校正。 ３．４基线校正处理效果判断 基线校正效果可以用来大致衡量的方法有３ 种。 １）用速度时程最终速度是否为零来衡量校正 处理效果。地震过后，地面运动终止，因此最终速度 为零，有时也可以分析积分求得的位移时程的合理 性来进行衡量。 ２）地震过后，地面运动终止。因此，对于有地 震事后记录的加速度，其位移时程震后的部分应与 基线平行。 ３）参考台站附近ＧＰＳ测得的同震变形，衡量位 移时程中的永久位移值。使用这种方法时，要了解 ＧＰＳ测得的同震变形中是否含有震后场地的蠕变变 形。 ４汶川地震加速度记录基线漂移校正 汶川Ｍｓ８．０地震中，由于噪音干扰，某些台站 记录的基线出现较为明显的漂移。图３～５中未校 正的记录曲线为此次地震中３条较为典型的加速度 基线漂移记录。针对这３条记录各自基线漂移特 性，分别选用了上述３种不同的基线校正方法进行 校正。 图３为四川大邑台基线校正前、后的加速度记 录图。从图中可以看出，虽然原始加速度时程基线 偏离零线，但与零线基本平行。针对这类记录基线 的漂移校正，用基线初始化进行即可获得很好的效 果。 图４为安徽六安台基线校正前、后的加速度记 录图。原始加速度时程基线明显呈线性漂移，对于 这类记录基线的漂移，用线性方程拟合加速度时程 的方法比较合适。 图５为湖北恩施台基线校正前、后的加速度记 录。原始加速度时程中存在很明显的长周期干扰噪 音，这类记录基线的漂移，只能采用低频滤波进行校 正。 在汶川地震加速度记录中，许多加速度记录基 线漂移不具有上述记录的明显特征，其校正一般需 要采用３种方法进行组合使用，才能达到基线校正 的较好效果。 ５０ Ⅱ。＼ 大地测量与地球动力学 最鞭 ⅡＩ。＼ 曩 甚＼ 最鞍 ＼ ３０卷 ｌｏｏ ５０ Ｏ －５０ ．１００ ５０ ０ ＩｌｌＩｌＩＩ ＩＩＩ』ＩＬｌ‘Ｉ－．Ｉ …………一 町 ＿５Ｏ ０ ２０ ４０ ６０ ８ｏ １００１２０１４０ｌ６０ ｌ８０２００２２０２４０２６０２８ｏ ｆ／Ｓ 图３汶川地震大邑台ＥＷ向基线校正前、后的记录曲线 Ｆｉｇ．３ Ｕｎｃｏｒｒｅｃｔｅｄ ａｎｄ ｃｏｒｒｅｃｔｅｄ ＥＷ ｃｏｍｐｏｎｅｎｔ ｏｆ ａｃｃｅｌ— ｅｒａｔｉｏｎ ｏｆ Ｗｅｎｃｈｕａｎ Ｍｓ８．０ ｅａｒｔｈｑｕａｋｅ ａｔ Ｄａｙｉ ｓｔａ— ｔｉｏｎ ｌ ０．５ Ｏ －０．５ ．１ ．１５ ０ ４０ ８０ １２０ １６０ ２０ｏ ２４０ ２８０ ３２０ ｔ，ｓ 图４汶川地震六安台ＥＷ向基线校正前、后的记录曲线 Ｆｉｇ．４ Ｕｎｃｏｒｒｅｃｔｅｄ ａｎｄ ｃｏｒｒｅｃｔｅｄ ＥＷ ｃｏｍｐｏｎｅｎｔ ｏｆ ａｃｃｅｌ— ｅｒａｔｉｏｎ ｏｆ Ｗｅｎｃｈｕａｎ Ｍｓ８．０ ｅａｒｔｈｑｕａｋｅ ａｔ Ｌｉｕａｎ ｓｔａｔｉｏｎ ； ４ ６ ４Ｌ Ｌ ｒ ０ ２ ４ ６ ８ １０ １２ １４ １６ １８ ２０ ２２ ２４ ２６ ２５ ｔ／ｓ 图５汶川地震恩施台基线校正前、后的记录曲线 Ｆｉｇ．５ Ｕｎｃｏｒｒｅｃｔｅｄ ａｎｄ ｃｏｒｒｅｃｔｅｄ ＥＷ ｃｏｍｐｏｎｅｎｔ ｏｆ ａｃｃｅｌ— ｅｒａｔｉｏｎ ｏｆ Ｗｅｎｃｈｕａｎ Ｍｓ８．０ ｅａｒｔｈｑｕａｋｅ ａｔ Ｅｎｓｈｉ ｓｔａｔｉｏｎ ５结论 １）基线校正处理没有固定模式，需依据各记录 基线漂移的具体特性而定。 ２）用拟合方法进行基线校正时，存在一定的拟 合误差，会导致部分地震信息失真。 ３）用低频滤波方法进行基线校正时，会丢掉低 频滤波频段的地震信息，即部分长周期信息会被滤 掉。而长周期地震信息对大跨度桥梁、管线工程等 抗震计算十分重要，因此，在选择低频滤波进行基线 校正时，需考虑加速度时程的用途。 ４）加速度记录的低频造影干扰复杂且无明确 的规律，基线校正不可能完全清除这些噪声的干扰。 从图５可以看出，校正后的恩施台ＥＷ向加速度记 录基线基本上与零线重合，但时程曲线中仍然存在 幅值很小的高频干扰。 参考文献 ｌ谢礼立，于双久．强震观测与分析原理［Ｍ］．北京：地震出 版社，１９８２．（Ｘｉｅ Ｌｉｌｉ ａｎｄ Ｙｕ Ｓｈｕａｎ￣ｉｕ．Ｓｔｒｏｎｇ ｍｏｔｉｏｎ ｏｂ— ｓｅｒｖａｔｉｏｎ ａｎｄ ａｎａｌｙｓｉｓ ｔｈｅｏｒｙ［Ｍ］．ＢｅＯｉｎｇ：Ｓｅｉｓｍｏｌｏｇｉｃａｌ Ｐｒｅｓｓ，１９８２） ２ Ｄａｖｉｄ Ｍ Ｂｏｏｒｅ ａｎｄ Ｊｕｌｉａｎ Ｊ Ｂｏｍｍｅｒｂ．Ｐｒｏｃｅｓｓｉｎｇ ｏｆ ｓｔｒｏｎｇ． ｍｏｔｉｏｎ ａｃｃｅｌｅｒｏｇｒａｍｓ：ｎｅｅｄｓ，ｏｐｔｉｏｎｓ ａｎｄ ｃｏｎｓｅｑｕｅｎｃｅｓ［Ｊ］． Ｓｏｉｌ Ｄｙｎａｍｉｃｓ ａｎｄ Ｅａｒｔｈｑｕａｋｅ Ｅｎｇｉｎｅｅｒｉｎｇ，２００５，２５：９３— １１５． ３ Ｂｏｏｒｅ Ｄ Ｍ．Ｅｆｆｅｃｔ ｏｆ ｂａｓｅｌｉｎｅ ｃｏｒｒｅｃｔｉｏｎｓ ｏｎ ｒｅｓｐｏｎｓｅ ｓｐｅｃｔｒａ ｆｏｒ ｔｗｏ Ｒｅｃｏｒｄｉｎｇｓ ｏｆ ｔｈｅ １９９９ Ｃｈｉ—Ｃｈｉ，Ｔａｉｗａｎ，ｅａｒｔｈｑｕａｋｅ ［Ｒ］．ＵＳＡ Ｇｅｏｌｏｇｉｃａｌ Ｓｕｒｖｅｙ，１９９９． ４ Ｄａｖｉｄ Ｍ．ｅｔ ａ１．Ｃｏｍｍｅｎｔｓ ｏｎ ｂａｓｅｌｉｎｅ ｃｏｒｒｅｃｔｉｏｎ ｏｆ ｄｉｉｇｔａｌ ｓｔｒｏｎｇ—ｍｏｔｉｏｎ ｄａｔａ：ｅｘａｍｐｌｅｓ ｆｒｏｍ ｔｈｅ １９９９ Ｈｅｃｔｏｒ Ｍｉｎｅ， Ｃａｌｉｆｏｒｎｉａ，Ｅａｒｔｈｑｕａｋｅ［Ｊ］．Ｂｕｌｌｅｔｉｎ ｏｆ ｔｈｅ Ｓｅｉｓｍｏｌｏｇｉｃａｌ Ｓｏ— ｃｉｅｔｙ ｏｆ Ａｍｅｉｒｃａ，２００２，９２（４）：１ ５４３—１ ５６０． ５ Ｓｉｎａｎ Ａｋｋａｒ ａｎｄ Ｄａｖｉｄ Ｍ Ｂｏｏｒｅ．Ｏｎ ｂａｓｅｌｉｎｅ ｃｏｒｒｅｃｔｉｏｎｓ ａｎｄ ｕｎｃｅｒｔａｉｎｔｙ ｉｎ ｒｅｓｐｏｎｓｅ ｓｐｅｃｔｒａ ｆｏｒ ｂａｓｅｌｉｎｅ ｖａｒｉａｔｉｏｎｓ ｃｏｍｍｏｎｌｙ ｅｎｃｏｕｎｔｅｒｅｄ ｉｎ ｄｉｇｉｔａｌ ａｃｃｅｌｅｒｏｇｒａｐｈ ｒｅｃｏｒｄｓ［Ｊ］． Ｂｕｌｌｅｔｉｎ ｏｆ ｔｈｅ Ｓｅｉｓｍｏｌｏｇｉｃａｌ Ｓｏｃｉｅｔｙ ｏｆ Ａｍｅｒｉｃａ，２００９，９９ （３）：１ ６７１—１ ６９０． ６ １ｗａｎ Ｄ，ｅｔ ａ１．Ｓｔｒｏｎｇ ｍｏｔｉｏｎ ｅａｒｔｈｑｕａｋｅ ｍｅａｓｕｒｅｍｅｎｔ ｕｓｉｎｇ ａ ｄｉｇｉｔａｌ ａｃｃｅｌｅｏｒｇｒａｐｈ［Ｒ］．Ｔｈｅ Ｎａｔｉｏｎａｌ Ｓｃｉｅｎｃｅ Ｆｏｕｎｄａｔｉｏｎ， １９８４． ７ Ｖｌａｄｉｍｉｒ Ｍ ａｎｄ Ｇｒａｉｚｅｒ．Ｅｆｆｅｃｔ ｏｆ ｔｉｌｔ ｏｎ ｓｔｒｏｎｇ ｍｏｔｉｏｎ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ［Ｊ］．Ｓｏｉｌ Ｄｙｎａｍｉｃｓ ａｎｄ Ｅａｒｔｈｑｕａｋｅ Ｅｎｇｉｎｅｅｒｉｎｇ， ２００５，２５：１９７—２０４． ８ Ｃｏｎｖｅｒｓｅ Ａ Ｍ．ＢＡＰ：Ｂａｓｉｃ ｓｔｒｏｎｇ－ｍｏｔｉｏｎ ａｃｃｅｌｅｒｏｇｒａｍ ｐｒｏ— ｃｅｓｓｉｎｇ ｓｏｆｔｗａｒｅ，ｖｅｒｓｉｏｎ １．０，Ｕ．Ｓ．Ｇｅｏ１．Ｓｕｒｖ．Ｏｐｅｎ－Ｆｉｌｅ Ｒｅｐｔ．［Ｒ］．９２—２９６． ９ Ｄａｖｉｄ Ｍ Ｂｏｏｒｅ．Ｏｎ ｐａｄｓ ａｎｄ ｆｉｈｅｒｓ：ｐｒｏｃｅｓｓｉｎｇ ｓｔｏｒｎｇ—ｍｏ－ ｔｉｏｎ ｄａｔａ［Ｊ］．Ｂｕｌｌｅｔｉｎ ｏｆ ｔｈｅ Ｓｅｉｓｍｏｌｏｇｉｃａｌ Ｓｏｃｉｅｔｙ ｏｆ Ａｍｅｒｉ— ｃａ，２００５，９５（２）：７４５—７５０． １０ Ｇｕｅｒｇａｎａ Ｍｏｌｌｏｖａ．Ｅｆｆｅｃｔｓ ｏｆ ｄｉｇｉｔｌａ ｆｉｌｔｅｒｉｎｇ ｉｎ ｄａｔａ ｐｒｏ— ｃｅｓｓｉｎｇ ｏｆ ｓｅｉｓｍｉｃ ａｃｃｅｌｅｒａｔｉｏｎ ｒｅｃｏｒｄｓ［Ｒ］．Ｈｉｎｄａｗｉ Ｐｕｂ— ｌｉｓｈｉｎｇ Ｃｏｒｐｏｒａｔｉｏｎ ＥＵＲＡＳＩＰ Ｊｏｕｒｎａｌ ｏｎ Ａｄｖａｎｃｅｓ ｉｎ Ｓｉｇ— ｎａｌ Ｐｒｏｃｅｓｓｉｎｇ，２００７．