Opportunities and Challenges of Strong Motion Observation Instruments
-
摘要: 本文介绍了国内外强震仪、力平衡加速度计、烈度计、MEMS加速度仪等强震动观测仪器的发展历史与现状，概述了强震动记录中典型异常波形产生的原因和力平衡加速度计的仪器响应误差及校正，建议了解决强震仪器缺陷的措施，分析了强震仪在强震动观测发展中面临的机遇和挑战。研究结果如下：① 进行强震动记录异常波形的研究是有针对性改进现有强震仪存在问题的很好途径；② 早期数字强震仪（力平衡加速度计频带范围0至30Hz）获取的强震动记录应进行仪器校正；③ 应跟踪强震动观测新方法和相关领域新技术，发展MEMS加速度仪和光纤加速度仪等新型仪器。Abstract: The development history of strong motion instruments, force-balance accelerometers, seismic intensity meter and MEMS accelerometers is described, and typical abnormal waveforms in strong motion records and instrument response errors correction for force-balance accelerometers are analyzed in brief, then challenges and opportunities of strong motion instruments are studied. The measures to solve the defects of strong motion instruments are suggested. The results are showed as following: (1) it is a good way to promote the development of strong motion instruments by studying abnormal waveforms in strong motion records; (2) for force-balance accelerometers with the band range from 0 to 30Hz, strong motion records obtained by early digital strong motion instruments should be corrected by instrument response; (3) new methods and technology in strong motion observation should be tracked, and new instruments (MEMS accelerometers and Opticalfibre accelerometers) should be developed to meet the demands of our country.
-
引言
桥梁是为保证道路的连续性而专门建造的人工构筑物，也是城市生命线工程中极其重要的一部分。随着我国国民经济的不断发展，桥梁已成为城市正常运行的不可或缺的建筑物。我国是多地震国家，地震的突发性和强破坏性，往往给我国城市桥梁带来巨大的损害。桥梁是城市生命线工程和城市交通系统的重要组成部分，在交通系统防震减灾工作中处于核心地位 (王东升等，2001)，如果桥梁在地震中受到严重的破坏，将会严重阻碍地震应急和灾区的救援，从而严重影响人们的生产生活和灾后的恢复重建工作。有效的桥梁震害预测可以提高城市综合抗震救灾能力，为城市抗震防灾规划提供科学依据，因此具有重要的现实意义。
国内外对桥梁震害预测的方法主要有4类 (王东升等，2001；2003；姜淑珍等，2004)：经验统计法 (包括久保庆三郎方法、日本土木工程学会方法、朱美珍方法、Buckle方法、回归统计法等)，规范校核法，Pushover (推覆分析) 法和大跨度桥梁定性与定量相结合的震害预测方法。庄丽等 (2007)基于VB语言和经验统计公式对青岛市桥梁做出了震害预测，并将该方法运用到青岛市地震应急指挥系统；柳春光等 (2008)将遗传算法和BP神经网络相结合，克服了BP神经网络容易陷入局部最优、运算速度慢等弱点，建立了遗传优化BP神经网络的桥梁震害预测方法；谭潇等 (2013)用支持向量机 (SVM) 原理建立了桥梁震害评估的SVM模型，同时也证明了SVM模型对桥梁震害的预测优于基于人工神经网的预测模型。本文充分利用粒子群算法和支持向量机的优点，借助Matlab软件和Libsvm工具箱，结合桥梁的震害影响因素，建立基于PSO-SVM的城市桥梁震害模型，使桥梁的震害预测方法更加完善和多元。
1. PSO的基本知识
1.1 粒子群算法的基本原理
粒子群算法 (Particle Swarm Optimization，PSO) 是一种结构简单、易于实现、通用的智能算法。通常粒子群 (田雨波，2104) 的数学描述为：假设在一个n维的搜索空 间中由m个粒子组成的种群，其中第i个粒子位置${x_i} = {({x_{i, 1}}, {x_{i, 2}}, \ldots, {x_{i, m}})^{\rm{T}}} $，其粒子速度${v_i} = {({v_{i, 1}}, {v_{i, 2}}, \ldots, {v_{i, n}})^{\rm{T}}}$。它的一个极值为${P_i} = {({p_{i, 1}}, {p_{i, 2}}, \ldots, {p_{i, n}})^{\rm{T}}}$，种群的全局极值为${P_{\rm{g}}} = {({p_{{\rm{g}}, 1}}, {p_{{\rm{g}}, 2}}, \ldots, {p_{{\rm{g}}, n}})^{\rm{T}}}$。粒子群算法在找到以上的2个极值后，根据公式 (1) 和 (2) 不断更新自己的速度和位置：
$$ v_{i, d}^{k + 1} = v_{i, d}^k + {c_1} \cdot \;{r_1} \cdot \;(p_{i, d}^k - x_{i, d}^k) + {c_2} \cdot {r_2} \cdot \;(p_{g, d}^k - x_{i, d}^k) $$ (1) $$ x_{i, d}^{k + 1} = x_{i, d}^k + v_{i, d}^{k + 1} $$ (2) 其中，c1和c2常被称为学习因子或加速常数；r1、r2为介于0和1之间的随机数；$v_{i, d}^k$和$x_{i, d}^k$是粒子在第k次迭代中第d维的速度和位置；$p_{i, d}^k$是粒子在第k次迭代中第d维的个体极值的位置；与个体极值相对应的$p_{g, d}^k$是群体在第k次迭代中第d维的全局极值的位置。
1.2 粒子群的优化及其SVM的参数设置
支持向量机 (Support Vector Machine，SVM) 在解决小样本和非线性高维模式识别中具有明显优势。一个非线性的SVM模型主要通过把原始空间的低维非线性数据映射到高维空间，将其转化为线性可分问题 (王建国等，2015)，这样在高维的空间中就对应其低维的非线性回归问题，其本质就是在训练样本中来寻找一个最优超平面，是一个二次规划问题，可以通过对偶问题来求解。
其具体的步骤是通过一定的映射关系φ(·)，在一定的条件下构成核函数$K({x_i}, {x_j}) = (\varphi ({x_i}) \cdot \varphi ({x_j}))$以避免高维空间的复杂计算。设样本集 (王雪刚，2014；王书舟，2009)：$\{ (xi, yi), i = 1, 2, \cdots, l, xi \in {R^n}, yi \in R\} $，R为欧式空间，n为样本输入维数，对于引入松弛因子${\xi _i}、{\zeta _i}$的二次优化问题：
$$ \min [\frac{1}{2}{\left\| \omega \right\|^2} + C\sum\limits_{i = 1}^l {({\xi _i} + {\zeta _i})} ] $$ (3) $$ {\rm{s}}{\rm{.t}}{\rm{.}}\left\{ \begin{array}{l} yi - [{\omega ^{\rm T}}\varphi (xi) + {\rm{b}}] \le \varepsilon + {\xi _i}\\ {\omega ^{\rm T}}\varphi (xi) + {\rm{b}} - yi \le \varepsilon + {\zeta _i}\\ \xi i, \zeta i \ge 0 \end{array} \right. $$ (4) 通过引入Lagrange乘子α, β, η，定义Lagrange函数的二次规划方程为：
$$ \begin{array}{l} \;\;L(\omega, b, \alpha, \beta) = \frac{1}{2}{\left\| \omega \right\|^2} + C\sum\limits_{i = 1}^l {({\xi _i} + {\zeta _i}) - \sum\limits_{i = 1}^l {{\alpha _i}({\xi _i} + \varepsilon - {y_i} + (\omega \cdot \varphi ({x_i})) + b)} } \\ - \sum\limits_{i = 1}^l {{\beta _i}} ({\zeta _i} + \varepsilon - {y_i} - (\omega \cdot \varphi ({x_i})) - b) - \sum\limits_{i = 1}^l {{\eta _i}({\xi _i} + {\zeta _i})} \end{array} $$ (5) 其中，$\alpha {\rm{, }}\beta {\rm{, }}\eta \ge 0$；xi，yi为模型的输入和输出；ω为权重向量，b为偏值；i=1，2，…，l，l为样本个数；C为惩罚参数，用于控制模型的复杂度和逼近误差；ε主要用于控制误差和模型的泛化能力。其函数的值条件为：
$$ \frac{{\partial L}}{{\partial \omega }} = 0 \to \omega = ({\beta _i} - {\alpha _i}) \cdot \varphi ({x_i}) $$ (6) $$ \frac{{\partial L}}{{\partial b}} = 0 \to \sum\limits_{i = 1}^l {({\beta _i} - {\alpha _i})} = 0 $$ (7) $$ \frac{{\partial L}}{{\partial {\xi _i}}} = 0 \to \gamma - {\alpha _i} - {\eta _i} = 0\;\;\;(i = 1, 2, \cdots, l) $$ (8) $$ \frac{{\partial L}}{{\partial {\zeta _i}}} = 0 \to \gamma - {\beta _i} - {\eta _i} = 0\;\;\;(i = 1, 2, \cdots, l) $$ (9) 解得：
$$ b = {y_i} - \sum\limits_{j = 1}^l {({\alpha _j} - {\beta _j})\; \cdot K({x_i}, {x_j}) - \varepsilon } $$ (10) 或：
$$ b = {y_i} - \sum\limits_{j = 1}^l {({\alpha _j} - {\beta _j}) \cdot K({x_i}, {x_j}) - \varepsilon } $$ (11) 根据这些条件可以得到引入Lagrange函数的二次规划的对偶问题：
$$ \mathop {{\rm{min}}}\limits_{{a_i}, {\beta _i}} \frac{1}{2}\sum\limits_{i = 1}^l {\sum\limits_{j = 1}^l {({\alpha _i} - {\beta _i})({\alpha _j} - {\beta _j})} } \cdot K({x_i}, {x_j}) - \sum\limits_{i = 1}^l {{y_i} \cdot ({\alpha _i} - {\beta _i}) + \varepsilon \sum\limits_{i = 1}^l {({\alpha _i} + {\beta _i})} } $$ (12) 令：
$$ \sum\limits_{i = 1}^l {({\alpha _i} - {\beta _i}) = 0} \;\;\left({{\alpha _i} \ge 0, C \ge {\beta _i}, i = 1, 2, \cdots l} \right) $$ (13) 在求得对偶问题的最优解，在最优解中对应的非零向量为α, β; $\alpha = (\alpha 1, \alpha 2, \cdot \cdot \cdot, \alpha_l)$, $\beta = (\beta 1, \beta 2, \cdot \cdot \cdot, \beta_l)$；进一步就得到该高维特征空间的最优超平面 (即非线性回归方程)：
$$ f(x) = \sum\limits_{i = 1}^l {({\alpha _i} - {\beta _i})} \cdot K({x_i}, {x_j}) + b $$ (14) 核函数为RBF (径向基函数)，采用其高斯形式 (白春华等，2013)：
$$ K(x, x') = {\rm{exp(}} - \frac{{{{\left\| {x - x'} \right\|}^2}}}{{2{\delta ^2}}}{\rm{)}} $$ (15) 在RBF为核函数的SVM模型中，SVM的性能在一定的方面容易受到惩罚参数C和核函数参数δ2的影响，在下文中核函数参数δ2用Γ表示，因此，结合粒子群算法的优点进行全局寻优来不断的优化惩罚参数C和核函数参数，使得预测结果更加准确。
2. 模型的建立
2.1 模型数据的输入
通过查阅资料，收集了唐山地震、丽江地震、汶川地震、海城地震、通海地震、台湾集集地震中的城市桥梁震害情况 (吴昊，2009；郭恩栋等，2014；刘恢先，1986)，根据经验统计法已有的研究成果，且城市桥梁建设广泛采用桩基础，所以模型中不考虑基础类型因素。模型选择影响桥梁震害等级的8个主要因素主要包括地震烈度、场地土类型、地基失效程度、上部结构、墩台高度、支座型式、桥梁跨数和桥梁的建成年代，用以上8个因素作为特征输入向量，并结合城市桥梁的相关特征。由于PSO-SVM模型只能对已量化的数字信息进行处理，但本文所选择的影响桥梁震害的8个因素都是非量化的，所以采用吴昊 (2009)对桥梁震害的影响因素进行量化的方法得到相应的值。城市桥梁主要震害因素和其量化值如表 1。
表 1 震害影响因素指标及量化值Table 1. Quantitative value of disaster factors影响因素 量化值 地震烈度 7度 1.0 8度 1.1 9度 1.2 10度 1.5 场地土类型 Ⅰ类 0.8 Ⅱ、Ⅲ类 1 Ⅳ类 1.2 场地失效程度 无 1.0 轻 1.5 重 1.8 上部结构 简支梁 1.0 连续梁 1.4 墩台高度 ≤5m 1.0 5-10m 1.1 ≥10m 1.2 支座型式 有防落梁措施 1.0 消能支座 1.2 一般支座 1.4 桥梁长度 跨数≤3 1 跨数＞3 1.2 桥梁的建成年代 1987年以前 1.1 1987-2006年 1.0 2006年以后 0.9 2.2 模型数据的输出
根据国内现有桥梁地震破坏等级的分类标准，通常把桥梁震害分为基本完好、轻微破坏、中等破坏、严重破坏和毁坏5个等级 (王再荣等，2010)，各个震害等级下破坏现象的描述如表 2。由于模型只能对已量化的数字信息进行处理，因此把桥梁震害等级的5种情况--基本完好、轻微破坏、中等破坏、严重破坏和毁坏分别对应模型的输出结果为1、2、3、4、5(如表 3)。输出结果的每一个标识对应着一种桥梁的破坏情况。
表 2 桥梁破坏等级划分Table 2. Classification criterion of bridge damage破坏等级 震害描述 基本完好 承重结构完好，桥面无明显变化，只有个别构件轻微损坏，其他构件无损，可正常通行 轻微破坏 桥台、桥面、桥墩、桥拱、桥塔、主梁等的混凝土部件表面出现细裂缝，非承重构件有破坏，桥面伸缩有变化，梁有轻微的移动，腹拱有细微裂缝，拱身有微裂，拱肋与拱波关系处松脱，墩台轻微变位，台背填土下沉等，可照常使用。 中等破坏 桥墩混凝土出现明显裂缝，梁端混凝土出现明显裂缝，支座与梁连接的螺栓部分剪断，梁移动，桥塔结构轻微变形，墩台轻微移动，出现明显裂缝梁，端混凝土脱落，腹拱错位，拱肋、拱波龟裂，墩台微移动或开裂，经适当修复可使用。 严重破坏 桥墩混凝土出现贯通裂缝、剥落，梁、拱出现贯通裂缝或破碎，桥塔结构变形，主要承重构件破坏、断裂，如梁、拱开裂或破碎，墩台滑移、断裂或严重倾斜，承载能力明显降低，必须大修。 毁坏 发生落梁、塌拱、墩台折断等破坏、倒塔，已不能使用，需重建。 表 3 桥梁震害等级标识Table 3. Classification criterion of bridge damage by earthquake桥梁破坏等级 基本完好 轻微破坏 中等破坏 严重破坏 倒毁 标识 1 2 3 4 5 2.3 粒子群优化算法
上文已经确定了影响支持向量机预测的基本因素是惩罚参数C和核函数的宽度，这两个参数的微妙变化都会影响预测结果的精度，本文主要比较用粒子群优化的前后模型预测准确度的变化，通过对结果的比较来研究更为准确的城市桥梁震害预测的模型。通过粒子群算法来寻找最佳的核函数参数和惩罚参数C，使模型的预测结果不断得到优化，相关的步骤的流程图如下：
图 1 PSO-SVM桥梁震害预测模型的构建Figure 1. Flowchart of PSO-SVM construction for seismic damage prediction model of urban bridge2.4 桥梁数据的归一化处理
在利用PSO-SVM模型进行桥梁震害预测时，需要对已量化的训练数据和测试数据进行归一化处理，这样就使得不同类型的数据之间有一定的可比性，同时在一定的程度上也消除了数据之间的量纲差异，使不同类型的桥梁震害影响因子的量化值处于同一个数量级，有利于对数据进行综合评价。
$$ {\overrightarrow x_{ik}} = \frac{{{X_{ik}} - X_k^{\min }}}{{X_k^{\max } - X_k^{\min }}} $$ (16) 公式 (16) 对数据进行线性转换，将桥梁震害数据的量化值转换到[0，1]，${\overrightarrow x _{ik}}$表示第i个样本、第k个因素归一化指标值，Xik表示第i个样本、第k个因素的量化值，$X_k^{\max }$和$X_k^{{\rm{min}}}$分别为k个因素列中的最大值和最小值。
2.5 PSO-SVM的迭代的收敛过程
利用Libsvm工具箱及Matlab软件设置PSO-SVM的初始参数。设粒子的种群s=40，最大迭代次数kmax=300，学习因子c1，c2取值分别为1.2和1.5。图 2可以看出粒子群的适应度在前期的变化的幅度比较大，这样可以保证粒子群算法具有较好的全局寻优能力，避免寻优结果达到局部最优。后期粒子寻优的最佳适应度的变化幅度较小，则表明粒子快速收敛已达到最优。同样粒子群优化的适应度是判定训练样本位置优化好坏的有力标准，优化后的惩罚参数C和Γ分别取值为53.88和0.01，从图 2中可以看出，经过参数优化的支持向量机预测模型具有较好的收敛效果。
3. 实验仿真结果与分析
PSO-SVM预测模型选择123个具有代表性的桥梁样本，其中103个作为训练样本，20个作为测试样本用于仿真检验。借助Matlab及相关软件对20个测试样本进行预测，将模型的预测结果与实际桥梁震害结果进行比较，更能直接地观测实验结果的准确度。相应的实验预测结果在图 3和图 4中进行了详细的描述。从表 4中得知，PSO-SVM模型预测结果只有序号为4、15、18的3座桥梁的震害与实际震害结果不一致，而SVM模型的预测结果有序号为4、6、9、15、18、20的6座桥梁与实际震害结果有差异。
表 4 模型预测结果差异项比较Table 4. Comparison of predicted results from various models预测模型 预测结果差异项序列号 差异项的预测震害指数 差异项的实际震害指数 模型预测精度 SVM 4 3 5 70% 6 3 4 9 2 3 15 2 3 18 2 4 20 3 4 PSO-SVM 4 4 5 85% 15 2 3 18 3 4 从实验的结果来看，PSO-SVM和SVM的预测结果都和桥梁震害的实际结果进行了比较，PSO-SVM的计算结果的精度高于SVM。为使结果更易于清晰可见，特制作表 4来对比2种模型的预测结果。
从表 4可以看出，PSO-SVM模型对桥梁的预测结果明显优于没有被PSO算法优化过的预测模型，同时模型没有出现过拟合现象，从而可以得知PSO-SVM预测模型用于城市桥梁震害预测和评估是切实可行的。
4. 结论
本文通过对比基于SVM和PSO-SVM的桥梁群体震害预测模型对测试桥梁样本的预测结果，结论显示，PSO-SVM模型不仅提高了城市群体桥梁震害预测的准确率，也给出了一种更为科学有效的桥梁震害评估模型，该模型充分利用粒子群算法和支持向量机算法的优点不断优化惩罚参数C和核函数参数δ2，使桥梁震害预测更准确。
但是，本文所提出的PSO-SVM桥梁群体震害评估模型还存在一定的不足之处：① PSO-SVM的桥梁震害模型的预测结果是相对准确的，但是其计算时间较长，不如SVM模型省时；② 模型训练样本和测试样本都基于小样本系统，没有经过较多样本的检验与测试。
-
图 1 第一台强震仪内部结构图（Trifunac，2009）
Figure 1. The internal structure of the first strong motion accelerograph (Trifunac, 2009)
-
段向胜, 周锡元, 2010.土木工程监测与健康诊断原理、方法及工程实例.北京:中国建筑工业出版社. 高光伊, 于海英, 李山有, 2001.中国大陆强震观测.世界地震工程, 17(4):13-18. http://www.cnki.com.cn/Article/CJFDTOTAL-SJDC200104002.htm 何先龙, 赵立珍, 2010.基于多重互相关函数分析剪切波速.岩土力学, 31(8):2541-2545. http://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201008032.htm 李沙白, 1987.强震加速度仪的发展近况.世界地震工程, 1:28-34. http://www.cnki.com.cn/Article/CJFDTOTAL-SJDC198701005.htm 李沙白, 1998.强震数字技术的最新进展及其在大震快速反应中的应用.自然灾害学报, 7(2):79-85. http://www.cnki.com.cn/Article/CJFDTOTAL-ZRZH802.011.htm 李小军, 温瑞智, 于海英, 2009.汶川8.0级地震余震固定台站观测未校正加速度记录.北京:地震出版社. 梁永朵, 姜金征, 李莹等, 2015.辽宁灯塔5.1级强震记录分析.地震地磁观测与研究, 36(5):25-29. http://www.cnki.com.cn/Article/CJFDTOTAL-DZGJ201505006.htm 曲明哲, 罗新恒, 温瑞智, 2014.MEMS加速度传感器噪声分析.地震工程与工程振动, 34(6):26-32. http://www.cnki.com.cn/Article/CJFDTOTAL-DGGC201406004.htm 王玉石, 李小军, 梅泽洪, 刘燕, 2013.几种仪器烈度算法在汶川地震与芦山地震中的可靠性比较.地震学报, 35(5):759-770. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXB201305014.htm 谢礼立, 彭克中, 1984. 强震动观测的数字纪元. 国际地震动态, (7): 4-8, 32-33. 王家行, 胡振荣, 1997.SLJ型宽频带大动态力平衡三分向加速度计的设计与研制.地震地磁观测与研究, 18(5):46-52. http://www.cnki.com.cn/Article/CJFDTOTAL-DZGJ705.007.htm 于海英, 2006.力平衡式加速度计仪器响应误差分析和校正方法.地震工程与工程振动, 26(6):200-203. http://www.cnki.com.cn/Article/CJFDTOTAL-DGGC200606032.htm 于海英, 韦良杰, 2007.强震动记录误差分析和校正方法综述.仪器仪表学报, 28(4 Supp1.1):175-177. http://www.cnki.com.cn/Article/CJFDTOTAL-DGGC200906000.htm 于海英, 杨永强, 王栋, 周宝峰, 解全才等, 2008. 汶川8. 0级地震强震动加速度记录初步分析处理. 汶川地震工程震害调查交流研讨会报告, 成都. 喻畑, 陆鸣, 2013.随机方法在地震烈度速报中的应用.地震地磁观测与研究, 34(5/6):271-276. http://www.cnki.com.cn/Article/CJFDTOTAL-DZGJ2013Z3050.htm 袁一凡, 1998.由地震动三要素确定地震动强度 (烈度) 的研究.哈尔滨:国家地震局工程力学研究所. 苑立波, 耿涛, 杨军, 2008. 双芯光纤集成式加速度计及测量方法: 中国, CN101368978A. 2009-02-18. 中国地震局震害防御司, 2008.汶川8.0级地震未校正加速度记录.北京:地震出版社. 周宝峰, 2012.强震动观测中的关键技术研究.哈尔滨:中国地震局工程力学研究所. 周雍年, 2011.强震动观测技术.北京:地震出版社. 朱建刚, 周朝晖, 赖敏, 黎大虎, 2006.有缺陷强震仪的使用及其数据处理.纪念松潘平武地震30周年交流文集:272-276. http://cdmd.cnki.com.cn/Article/CDMD-10486-1015305218.htm Trifunac M. D., 2009. 75th anniversary of strong motion observation-a historical review. Soil Dynamics and Earthquake Engineering, 29:591-606. doi: 10.1016/j.soildyn.2008.05.011 Kashima K., 2000. Strong Earthquake Motion Observation in Japan.//http://iisee.kenken.go.jp/staff/kashima/soa2000/soa.htm. -
下载:
下载:
点击查看大图
计量
- 文章访问数: 190
- HTML全文浏览量: 19
- PDF下载量: 26
- 被引次数: 0
