确定地壳形变观测曲线趋势转折点的定量方法与实例.pdfVIP

第27卷 第3期 内 陆 地 震 V01．27 No．3 2013年 9月 INLAND EARTHQUAKE Sep． 20l3 文章编号：1001-8956(2013)03-0251-06 中图分类号：P315．72 文献标识码：A 确定地壳形变观测 曲线趋势转折点 的定量方法与实例 刘冠中，刘远征 (1．中国地震局地壳应力研究所 ，北京 100085；2．华北水利水电大学资源与环境学院，河南 郑州 450008) 摘要 ：趋势转折是地壳形变变化的一个重要信息，为定量 、精确地提取地壳形变观测资料的趋势转 折点 ，基于改进的有序聚类分析法 ，采用时窗滑动的方法，提取时间序列的趋势转折点，并编制了相 应的程序。以新源地震台地倾斜南北向1992—2006年的观测数据为例 ，以4年窗长、2年步长滑动 计算其趋势转折点。研究结果表明：(1)采用改进的有序聚类分析法提取时间序列的趋势转折点具 有明确的统计意义，可减少人为判断的主观和任意性。(2)该方法算法简单 ，易于实现，可用于 日常 的地壳形变观测资料的分析预报工作。(3)新源地震台地倾斜南北向观测数据计算结果表明，在 2001年之前 ，趋势转折比较频繁，变化也 比较明显，而在2001年之后，趋势转折较少，变化更接近线 性 关键词 ：地壳形变 ；趋势转折 ；地倾斜；新源地震台 地壳形变观测结果提供了具有明确物理意义的地壳运动的直接证据和定量特征，在防震 减灾工作中是一项不可缺少的重要手段 ¨J。在地壳形变资料的分析 中，观测曲线的趋势转折 (上升、下降或转平)是一个重要的信息，它表示观测的物理量偏离了原有的长期背景，可能说 明观测点的地壳形变发生了较大的变化或异常，但是，在转折时间的精确定位方面，目前主要 依靠形态分析和经验判定 2 ，比较缺少数学上的、定量的方法。对于时间序列的特征点，如极 值点、拐点的位置，在局部范围内可以通过求导获得，但如果在某一时段存在多个局部的特征 点，就需要利用最优化的方法寻找最佳分割点，最大程度地反映转折点两侧不同的变化趋势。 针对这一问题，为减少主观判断带来的偏差，提取客观实际的转折点，本文将在陈远中等 关 于趋势转折点的研究基础上，研究地壳形变观测曲线趋势转折时间的定量判别方法，并编制相 应的可 以滑动计算的程序。 新疆新源地倾斜仪 自1992年观测至今，积累了 10多年精度较高的连续资料，趋势变化 比 较显著，适合作提取趋势转折点的实例，并且，在观测期间台站周边发生多次5级以上中强震， 分析其趋势转折的特征对震例研究具有参考意义。如果该方法的可靠性和有效性得以验证 ， 可为地壳形变异常信息提取、异常时间判定提供支持。 收稿 日期：2013-01—10；修 回日期 ：2013-05-31． 课题项 目：中国地震局地壳应力研究所基本科研业务专项 “鲜水河断裂带断层蠕变与复杂几何结构的关系及非稳态蠕 滑事件的动力学含义”(ZDJ2Ol1—18)项 目资助． 作者简介：刘冠中(1980～)，男，助理研究员，硕士 ，2006年毕业于中国地震局地壳应力研究所 ，现主要从事地壳形变观 测与地震研究 ．E—mail：bjyanjiu@163．Corn 252 内 陆 地 震 27卷 1 方法 陈远 中等对有序聚类分析法进行了改进 ，使其更加实用于序列转折点或突变点的提 取 。其基本思想简述如下，对于存在趋势变化的时间序列 y~(N为有限的数据长度)，取一 点 Yi(3≤i≤Ⅳ一2)将数据分为前后两段 ，采用最小二乘法分别对两段数据进行线性拟合，之 后求取前后两部分总的离差平方和，当总离差平方和最小时，对应的转折点为最优二分割点 (图1为示意图，纵坐标数值为新源南北 向倾斜量)，该方法主要利用 同类之间的离差平方和 最小，而类与类之间的离差平方和最大的分割思想。陈远 中等 已从数学上证 明该方法求得的 距离平方和更小，详细推导参见文献 [5]，由于本方法提取出的转折点具有统计意义，从