一种具有永久位移的强震动位移时程仿真方法技术

本发明专利技术公开了一种具有永久位移的强震动位移时程仿真方法，首先加速度记录基线初始化，然后采用二阶Butterworth低通滤波对基线初始化的加速度记录进行滤波；再然后对低通滤波后的加速度记录积分得到速度时程；最后计算位移时程末尾10％的位移平坦度。能够获得地震动的合理永久位移，防止地震动永久位移信息的损失和相关时程的畸变，消除了主观性，避免引入人为误差，较为客观准确地识别出永久位移。较为客观准确地识别出永久位移。

【技术实现步骤摘要】
一种具有永久位移的强震动位移时程仿真方法


[0001]本专利技术涉及一种具有永久位移的强震动位移时程仿真方法。

技术介绍

[0002]强烈地震发生时，除了地面强烈振动引起的灾害外，其伴随地面的位错、地破坏和山崩等地变动，对建筑构造物也引起了强大的破坏力。特别是地面在近距离高震级的地震作用下，较长时间的剧烈振动将导致地面产生不同程度的永久位移。这种现象在我国1999年的台湾集集地震、2008年的汶川地震，2010年的智利和新西兰地震，2011年的东日本大地震中都得到了证实。强烈地震时，地变动由于相互挤压变形、位差或基地破坏而导致建筑构筑物损坏或倒塌，生命线系统的损坏或断裂。因此，研究地震灾害除了研究地面的强烈振动外，还需研究地震动的永久位移特性及其分布，地表位移与相关的地下浅层断层错动也有着密不可分的关系。在我国的防震减灾事业中，要探讨地震灾害，强震时地震动永久位移的识别、特性及其分布也是重要的研究课题之一，永久位移的识别是其中关键的环节。根据当时地震灾情调查，发现评估地震灾情除了使用峰值加速度等传统方法外，最大地表位移或同震位移对震害造成的影响也不容忽视。有必要通过合理的理论研究，并结合强震数据提取地表永久位移，为工程建设及震源机制的研究提供必要的参考。
[0003]理论上，强震加速度记录在经过二次积分后可以得到合理的最大地表位移及可能发生的永久位移，但是由于强震过程中地表因强烈振动而抬升或倾斜，强震仪也会由于基线的漂移而失真。因此，在强震数据的处理过程，如果难以消除因强烈振动而导致地表倾斜抬升所带来的基线漂移，则无法获得合理的地表永久位移。并且现在的基线校正程序主观性较强，永久位移的识别结果因人而异。
[0004]总体说来，造成加速度记录基线偏移的影响因素有三类，有强震仪的原因，有环境背景噪声的原因，也有地震动的原因。初始参考基线通过加速度记录的适合部分来确定，或者是P波到达前的时间区间内的数据，如果没有事前记录，则是整条记录的平均值。通过用所有记录数据减去这个平均值完成基线初始化。然后，并不对加速度记录进行滤波，而是对记录进行积分得到速度，并检查速度的长周期漂移，可能表明存在参考基线的偏移，大多情况下，基线的偏移在加速度时程中并不很明显，对于基线偏移量比较小的数据来说，进行简单的滤波就可以有效地去除信号中的噪声，对于具有永久位移信息的记录而言，会削弱其中的永久位移信息，但是并非不可以采用滤波的方法进行永久位移的识别。在基线偏移比较大的情况下，会严重影响信号的长周期部分，为了对数据进行校正，需要对速度进行线性的分段拟合或低阶多项式的拟合，然后从加速时间序列减去这些拟合的偏移量。大多情况下，对于强震数据的处理，只能是一定程度上削弱了原始误差，由于处理方法本身就具有一定的主观性，因此很难避免不引入新的人为误差。这些综合因素的影响给地震动永久位移的识别带来了一定的挑战。因此，在近断层强震动记录的处理过程中，如何消除强震记录中的基线漂移成为识别永久位移的核心研究内容。

技术实现思路

[0005]为解决以上技术上的不足，本专利技术提供了一种具有永久位移的强震动位移时程仿真方法，能够较为客观准确地识别出永久位移。
[0006]本专利技术是通过以下措施实现的：
[0007]一种具有永久位移的强震动位移时程仿真方法，包括以下步骤：
[0008]步骤1，加速度记录基线初始化，即记录的每个采样点数值减去事前时间段T记录的平均值；
[0009]步骤2，采用二阶Butterworth低通滤波对基线初始化的加速度记录进行滤波；
[0010]步骤3，对低通滤波后的加速度记录积分得到速度时程；
[0011]步骤4，对于速度时程从事前记录时间的结束点作为第一个分界点t1，第二个分界点t2则搜索从第一个分界点的下一个采样点开始选取直到事后时间记录的开始点，两个分界点把记录分成了三部分，首尾两部分采用直线拟合，中间部分利用两点三次Hermite插值分别拟合；
[0012]步骤5，计算每个分界点t2所对应的积分位移时程；
[0013]步骤6，计算位移时程末尾10％的位移平坦度，对于t2时间点的确定，需要建立一个判断参数，作为地面推移至永久位移时间点的判断依据，也就是通过方程式自动迭代计算，搜索t1到记录结束这段时间内的时间点，使得t3时间点后的位移时程尽可能达到幅值不变的平坦直线状态；对于每一个t2值，都需要计算平坦度系数，用于显示校正位移在t3与记录结束之间的平坦程度；把t3与记录结束之间位移的平均值作为最终的永久位移，并计算其标准方差σ；利用最小二乘法对于t3与记录结束之间的位移作线性拟合，位移的线性校正系数为r，并计算其坡度，坡度系数为b，平坦度系数φ计算公式如下：
[0014][0015]其中，
[0016]x
i
代表校正后的位移时程末尾部分数据点，代表校正后的位移时程末尾部分数据点的均值，y
i
代表校正后的位移时程末尾部分数据点直线拟合后的各数据值，代表校正后的位移时程末尾部分数据点直线拟合后的各数据值的均值，n代表校正后的位移时程末尾部分数据点的个数；
[0017]步骤7，排除非正常位移时程，然后选取最大平坦度对应的位移时程即为校正后的位移时程，从而获得永久位移，与距离台站5km范围内同震位移台站获得的同震位移进行比对。
[0018]在步骤7中，排除非正常位移时程的依据如下：(1)加速度与速度时程的初末速度为零或接近零，并且不应出现明显的基线漂移现象；(2)位移时程在经过峰值位移PGD的对应时刻后，其位移时程趋势会趋于靠近时间轴，最终一段时间内保持不变，即平行于时间轴；(3)永久位移与峰值位移PGD的符号相一致，即同为正或同为负；(4)永久位移的绝对值
不超过10m；(5)峰值位移PGD出现的时刻位于PGA或PGV对应时刻附近，不出现在位移时程的末端。
[0019]在步骤6中，线性校正系数r用以判断数据的离散程度，坡度系数b为对于产生永久位移段的时程部分进行最小二乘法线性拟合所得直线的斜率，用以确定拟合直线与时间轴的位置关系；σ反映了t3时间点至记录结束这段时间内位移波形的变异系数，用于表征地震动位移接近稳定状态的程度；
[0020]在步骤1中，事前时间段T为10s或20s。
[0021]本专利技术的有益效果是：本专利技术基于低通滤波、Hermit插值、平坦度判定的方法用于强震动记录残余位移的识别，能够获得地震动的合理永久位移，防止地震动永久位移信息的损失和相关时程的畸变，消除了主观性，避免引入人为误差，较为客观准确地识别出永久位移。
具体实施方式
[0022]本专利技术公开了一种具有永久位移的强震动位移时程仿真方法，基于低通滤波、Hermit插值、平坦度判定的方法用于强震动记录残余位移的识别，要优于其它方法。具体包括以下步骤：
[0023]步骤1，加速度记录基线初始化，即记录的每个采样点数值减去事前时间段T记录的平均值；事前时间段T为10s或20s；
[0024]步骤2本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种具有永久位移的强震动位移时程仿真方法，其特征在于，包括以下步骤：步骤1，加速度记录基线初始化，即记录的每个采样点数值减去事前时间段T记录的平均值；步骤2，采用二阶Butterworth低通滤波对基线初始化的加速度记录进行滤波；步骤3，对低通滤波后的加速度记录积分得到速度时程；步骤4，对于速度时程从事前记录时间的结束点作为第一个分界点t1，第二个分界点t2则搜索从第一个分界点的下一个采样点开始选取直到事后时间记录的开始点，两个分界点把记录分成了三部分，首尾两部分采用直线拟合，中间部分利用两点三次Hermite插值分别拟合；步骤5，计算每个分界点t2所对应的积分位移时程；步骤6，计算位移时程末尾10％的位移平坦度，对于t2时间点的确定，需要建立一个判断参数，作为地面推移至永久位移时间点的判断依据，也就是通过方程式自动迭代计算，搜索t1到记录结束这段时间内的时间点，使得t3时间点后的位移时程尽可能达到幅值不变的平坦直线状态；对于每一个t2值，都需要计算平坦度系数，用于显示校正位移在t3与记录结束之间的平坦程度；把t3与记录结束之间位移的平均值作为最终的永久位移，并计算其标准方差σ；利用最小二乘法对于t3与记录结束之间的位移作线性拟合，位移的线性校正系数为r，并计算其坡度，坡度系数为b，平坦度系数φ计算公式如下：其中，x
i
代表校正后的位移时程末尾部分数据点，代表校正后的位移时程末尾部分数据点的均值，y
i
代表校...

【专利技术属性】
技术研发人员：周宝峰，
申请(专利权)人：中国地震局工程力学研究所，
类型：发明
国别省市：