古田地震水口水电站重力坝强震观测记录HHT研究.docVIP

古田地震水口水电站重力坝强震观测记录HHT研究摘要：利用HHT变换研究2008年3月6日福建古田ML46地震中，水口水电站重力坝强震反应台阵获取的强震反应观测资料的时频特性及重力坝结构动力特性，得到如下结论：（1）强震记录能量集中在0～15 Hz频段和10～15 s时段；（2）重力坝顺河向主频为37 Hz。分析结果对认识水库地震近场地震动特性和重力坝地震反应有一定的实际意义。 关键词：重力坝；边际谱；Hilbert能量谱；瞬时能量谱 中图分类号：P3159文献标识码：A文章编号：1000-0666（2012）02-0236-04 0引言 长期以来，地震工程领域对结构强震记录最常用的处理方法就是傅立叶分析和小波分析，但从根本上来说，这两种方法都建立在信号为稳态的基础上，处理结构强震记录这类非平稳性很强的信号可能会产生一些比较严重的问题，如能量泄漏等（公茂盛等，2007）。HHT（HilbertHuang Transform）变换是一种全新的数据处理方法，具有简单高效、自适应性强、高分辨率、多分辨率等优点（Huang et al，1998）。该方法提出后，在大气、海洋、医学和机械等各个工程领域得到迅速推广和应用。在地震工程领域中，王彬等（2005）、李洪泉等（2003）、李书进等（2004）、张郁山等（2005）对强震观测资料采用HHT变换识别结构动力特性进行了系统研究，并认为对于结构强震记录这种非平稳性很强的数据，可以用边际谱识别结构自振频率。鉴于此，本文首先对HHT变换原理作了介绍，然后用HHT方法研究2008年3月6日福建古田ML46地震水口水电站重力坝强震记录，得到了结构强震记录时频幅值三维分布特性、边际谱、Hilbert能量谱和瞬时能量谱；又根据坝基和坝顶强震记录Hilbert边际谱识别了水口水电站重力坝顺河向自振频率，将坝基和坝顶强震记录能量谱和瞬时能量谱进行对比并阐述了差异。 1HHT变换原理 HHT变换是由美国宇航局的Huang等（1998）提出的一种全新的信号处理方法。它的核心思想是依据数据本身的时间尺度特征将信号分解为有限个固有模态函数（intrinsic mode function，IMF），然后对各模态分量进行HHT变换构造解析信号，得到信号的瞬时频率和振幅，进而得到三维Hilbert时频谱h(ω)=∫0TH(ω,t)dt（7）将振幅的平方对时间积分，可定义Hilbert能量，它表达了每个频率在整个时间长度内所累积的能量ES(ω)=∫0TH2(ω,t)dt（8）将振幅的平方对频率积分，可定义瞬时能量IE（ω），它提供了信号能量随时间的变化情况IE(ω)=∫ωH2(ω,t)dω（9）地震研究35卷第2期方嘉治等：古田地震水口水电站重力坝强震观测记录HHT分析2强震记录及台阵简述 张力飞和邢国良（1996）就水口水电站重力坝强震反应台阵的台阵设计、监测仪器选型、强震设台、台阵验收等方面做了介绍，张锦福和蔡辉腾（2009）详细论述了此次地震的震源参数以及水口水电站重力坝强震反应台阵的强震记录情况，本文不再赘述。水口水电站重力坝高101 m，强震反应台阵由7个测点21个频道组成。每个测点由竖直向、顺河向（上游至下游）和横河向（左岸至右岸）三个频道组成，测点具体的布设方案如图1所示，其中测点C7布设在大坝左岸基岩自由场上。图1台阵布设方案示意图 Fig1Schematic diagram of the array installation图2C1测点顺河向记录及各IMF分量（单位：Gal） Fig2Riveralong records of strong motion and the IMF component decomposed by EMD for measuring point C1（unit：Gal） 3强震记录HHT分析 据2008年福建古田ML46地震强震记录可知，坝体中部C2测点和C3测点因仪器原因未能获取记录，而C5测点和C6测点位于坝体左侧，边界效应影响较大，因此笔者仅选取C1和C4测点的强震记录作为分析对象（张锦福，蔡辉腾，2009）。 31C1测点顺河向记录分析 C1测点的强震记录（顺河向）经EMD分解后，得到10个IMF分量和1个趋势项，如图2所示。EMD分解首先将记录中的高频分量分离出来，之后所分解出的IMF分量的频率依次降低，波长逐渐增加，残余项并不是常数，说明该记录有一个趋势。对EMD分解结果进行Hilbert分解后得到的振幅—频率—时间分布关系、边际谱、能量谱、瞬时能量谱分别如图3、图4所示。从图4a和图4b可以看出，大部分能量集中在0～15 Hz以内，主要能量突出体现在37 Hz，也即重力坝结构顺河向