云鹏电站水库诱发地震分析

云鹏电站水库诱发地震分析 3 云鹏水电站及其周围地震尾波 Q 值研究 c 赵小艳, 付 虹, 毛玉平, 邬成栋 ()云南省地震局 , 昆明 650224 摘要 : 利用云鹏水电站周边 5个数字地震台站记录到的 2000年 1月至 2008年 12月间的 108个地震的数字波形资 料 , 采用 A k i单次散射模型 , 计算了各台站周围 423条地震记录的尾波 Q值 。结果显示 : 云鹏水电站地区尾波0 0194 ( ) 的 Q 值与频率的关系为 Q f= 9012 f 。该地区的 Q值比云南其它构造活动强烈地区的要高 , 因此属于中等 c c 0 构造活动地区 。云鹏水电站蓄水后弥勒台记录到的 108个地震的尾波 Q值及分频 Q 值有所下降 , 但下降幅度不 0 c 大 , 明水库蓄水对库区及其附近区域地下介质性质有所影响 , 但影响的程度有限 。 关键词 : 云鹏水电站 ; 尾波 Q 值 ; A k i单次散射模型 c ( ) A 中图分类号 : P31513 文章编号 : 1000 - 0666 201001 - 0055 - 08 文献标识码 : 方法研究了非洲地区尾波 Q 值的大尺度横向变化 c ( ) 特征 ; 秦嘉政和阚荣举 1986 运用云南区域 台 0 引言 网短周期地震仪器 记录 到的 小 震尾 波资 料 , 估 算 地震尾波是由地下介质的不均匀性所激发的 ,出昆明周边地区频率约 110 H z的短周期波的地壳 它的存在可以认为是介质不均匀性的明显证据之一 。Q值的平均值为 191 , 证明昆明属于低 Q值地区 ;0 0 () ( ) A k i 1969在研究地震尾波时 , 提出了一种统计模 王伟君和刘杰 2004 在研究 1999 年辽宁岫岩地 型来解释尾波的形成机制 , 即由于地球介质的不均 震序列尾波 Q 值时 , 发现其在主震前后有较明显 c 匀性 , 地震波在传播过程中遇到不均匀介质时就会 的变化 ———震前增高 , 震后降低 。产生散射 , 台站记录到这些延迟了的散射波 , 形成 云南地区宽带数字化台网的建立 , 为我们深入 ( ) 尾波 。A k i和 Choue t 1975 在研究中 , 假设地震 研究区域尾波衰减特征提供了大量的数字化记录资 波在经过这些不均匀体时只发生一次散射 , 从而得 料 , 使我们有可能对该地区尾波的 Q值分布特征及 0 出了计算尾波功率谱的公式 。在上述尾波理论的基 其与区域构造的关系进行较为精细的研究 。如马宏 础上 , 尾波的研究工作得到了迅速发展 , 国内外地 ( ) 生等 2006用云南地区 22个区域数字地震台站记震学家从不同方面对尾波开展了大量研究 : 用尾波 录到的 5 668 个地震的数字波形资料 , 计算了各台 探讨小震震源频谱以及尾波形成的物理机制 , 运用 站周围地震记录的尾波 Q值 , 并从地质构造 、地震 c 尾波衰减的时空特征探讨地球内部介质非均匀性 ,活动以及大地热流分布等方面对 Q值分布特征予以 c 研究尾波 Q值平面分布与区域构造活动性的关系 ;( ) c 初步解释 。王伟君等 2006 用云南区域数字地震 大震前后区域性尾波 Q值随时间的变化特征在地震 台网记录到的 c 22个近震资料 , 研究了云南地区地震预测 中 的 应 用 等 , 并 取 得 了 一 些 研 究 成 果 。如 尾波在 1,20 Hz内 6个频段的衰减特征 , 并给出云 ( ) Su teau和 W h itcom b 1979 利用美国南加利福利亚 南地区尾波 Q值分布特征 , 即滇东的 Q值 >滇西 c c 地区 C IT台网记录到的 1972,1975 年的地震资料 , ( ) 的 Q值 >滇西南的 Q值 。秦嘉政等 2001、李白 c c 从理论上给出了用 M 计算尾波震级 M 的理论公 ( ) ( ) L C 基和 秦 嘉 政 2004 、钱 晓 东 等 2004 、李 琼 等 ( ) 式 ; X ie和 N u ttli 1988 根 据 尾波 随机 模 型及 单 ( ) 2007, 2008用云南区域台网及大震后围绕震中建 次各向同性散射模型 , 提出 了 用单 台记 录 叠加 频设的临时数字台网资料研究了不同区域或不同地震 谱比值法计算尾波 Q 值的方法 , 并采用层析成像c 序列的尾波衰减特征 。 3 收稿日期 : 2009 - 06 - 101 ( ) 基金项目 : 国家科技支撑计划 “水库诱发地震监测与分析处理技术研究 ” 2006BA C 01B 04 - 02 - 01 和云南省科技 ( ) 攻关项目 “复杂构造区水库地震监测预警研究 ” 2006SG 01联合资助 1 () 2007年 8,9月 , 在云南省文山州的丘北县与对方程式 1取对数 , 稍加整理 , 移项 , 可得 红河州的泸西县交界 地 区发 生了 一 群小 地震 。从 πf ) ( )( ) ( 2 ln [A f, t3 t ] = lnS f- t1震群所处的区域构造 位 臵来 看 , 北 东向 的 南盘 江 Q c 断裂带应为其发震断 裂 。震 群 所在 区域 被 云南 数( ) 可以看出方程式 2 的右边是时间 t的直线方程 , 条著名的断裂所包围 : 震群西边 85 km 即为 SN 向 我们 可 以 对 某 个 地 震 记 录 在 某 一 频 率 上 测 出的小江断裂带 , 震群西北 50 km 为 NN E 向的弥勒 ) ( A f, t和 t, 即可以求出直线的斜率K, 从而可以 —师宗断裂 ; 震群南边 50 km 为 NW 向的 文山 断 得到品质因子裂 。从更大的区 域 来 看 , 震 群 所 在 的 区 域 为 由 小 ( ) πQ f=f / k1 c ( )3 江断裂和红河断裂所包围的相对稳 定的 华 南块 体 () 苏有锦等 , 1999 。川滇菱块深部物质向东南 滑 ( ) 由方程式 3 , 可以求出不同频率的尾波 Q 值 , c η移 , 却遇到了较 为 坚 固 稳 定 的 华 南 断 块 , 从 而 形 ( )利用 Q f = Q f, 可以拟合求出 Q 与频率 f 的 0 c c () 成该区孕震的主要力源 徐彦 , 2002 。对区域应 关系 。具体计算步骤如下 : 力场的研究发现 , 该地 区应 力 场压 应力 优 势方 向 首先 , 对原始数据进行滤波 。共计算 9个中心 () 为 SSE向 赵小艳等 , 2007 。 由于该震群所处的滤波频 率 点 : 110 H z、 115 H z、 215 H z、 410 H z、 丘北 —泸西一带历来为云南 515 H z、 810 H z、 1010 H z、 1210 H z、 1810 H z。使 少震区 , 历史上从未有过如此密集且频繁的小震群用高斯滤波器滤波 , 其定义为 活动 , 因此震群发生后 , 云南省地震局分析预报中 ω ω2 - c心研究人员从各个方面对该震群进行了调查研究 。 (ω) ( )G = exp 1 4 α- ω c在调查过程中我们了解到 , 在震群所在地有 1 座云 () 鹏水电站 , 水 库 最低 水位 死 水 位 高程 877 m , αω πω 式中 , = 2f 为圆频率 , 是中心点圆频率 , c 3 正常蓄水位高程 902 m , 总库容 31663 亿 m。工程 为一常数 , 用来确定滤波器频带宽度 , 可在 1 ,99 于 2007年 3 月实现截流蓄水 , 至 6 月 22 日坝前水 α 间选择 , 值越大 , 频带宽度越窄 。本文中 = 40。 3 位高程达 89811 m , 库容 31293 亿 m。丘北震群最 第二步 , 去噪和滑移计算 。从尾波到时开始 , 早的地震发生在 2007年 3月 , 与水库蓄水时间十分 对滤波后的数据 , 取窗长为 2 s, 以 018 s间隔滑动 吻合 。因此 , 我们判断该震群为云鹏水电站所诱发 。 (计算 , 求出不同时间点 以窗的中心点为当前 时 为了研究云鹏水电站及其附近区域介质性质及其在 ) 间点 的平均振幅为 蓄水前后是否发生变化 , 我们利用震群周边的数字 2 2 1 / 2 () ( )A= A- A1 5 C T N 地震台记录到的 2000 ,2008年 108个小震的近场数 式中 , A是窗口内信号的均方根振幅 , A是噪声 T N 字化地震波形资料 , 采用尾波单次散射模型 , 研究 均方根振幅 。从记录开始处 , 选取 P 波到时前 2 s了该地区地震尾波在 9 个中心滤波频率点的垂直向 长的记录为噪声计算均方根振幅 。设窗口 T 内有 n 数据的衰减特征 。 ( ) xi = 1,, n , 则该窗口个离散数据 , 振幅为 i 的均方根振幅定义为1 计算方法 n 2 ?x i依据尾 波单 次 散射 模型 , 地震 波 尾 波 峰 值 振 i = 1 1 ( )A=6 T n ( ) 幅 A f, t与尾波流逝时间 t之间满足 π- 1 - f t /Q 第三步 , 选择 拟 合 开 始 点 和 窗 长 , 用 最 小 二 c ( ) ( ) ( ) A f, t= S fte 11 乘法计算各个频率点的 Q 值 。在拟合中 , 拟合开 c ( ) 其中 , A f, t为扣除了仪器响应和台站场地效应的( 始点和拟合窗长对 Q 值是有较大影响的 安艺敬 c ( ) ) 台站的尾波振幅观测谱 ; S f代表尾波波源因子 , 一 , 1991; 黄才中 , 葛焕称 , 1995 。在 2 T后 一S 它与震源和频率有关 ; f 为尾波频率 ; t为发震 时 定范围内选择开始点对 Q 影响相对要小些 。为了 c 间与 S波到时差 。根据经验 , 尾波开始于 2 倍的 S 让 Q 值有可比性 , 并且不受有突跳点数据记录的 c ( 波走时 , 即 t > 2R / v R 为 震 中 距 , v 为 S 波 速 ( ) ( ) 影响 , 编制程序时采用在 F t〃t - t分布图上 ,S ) 度 ; Q 为尾波衰减系数 。 c 手工选择拟合开始 点和 结束 点 , 同 时 还要 满足 以 第 1期赵小艳等 : 云鹏水电站及其周围地震尾波 Q 值研究 c 57 P、 S拟合窗长尽量保分布在震群周围且记录到该震群清晰的 下条件 :计算开始点在 2 T附近 , S ( ) 持一致 ,一般取 15,30 s。波到时的台站共有 5 个 , 分别是通海台 02 TH 、 ( ) ( ) 弥 勒 台 03ML 、马 龙 台 04MA 、文 山 台 ( ) ( ) 21W S、个旧台 23 GJ 。这 5 个 台 站的 仪器 的 2 地震资料及分析处理 采样率均为 50 H z, 在 011,20 H z具有速度平坦的 211 资料选取 响应 。表 1给出了这 5个台站的基本参数 。地震震 据云南地震台网测定 , 2007 年8 ,9 月在云南 中位臵和周围相应台站的分布如图 1 a 所示 。在上 丘北 —泸西一带发生了一群小地震() 图 1 a, 其中( ) 述 5个数字地震台中 , 弥勒地震台 03ML 距离 1 ,119 级 43 次 , 2 ,219 级 122 次 , 3 ,319 级 23 该震群最近 , 记录到的地震波形数据最清晰 。 ( 次 , 最大 为 2007 年 9 月 11 日 M 316 地 震 图L 表 1 震群周边数字地震台主要参数 ) 1 b。我们以震群的中心点为圆心 , 选取半径为 50 Tab11 M a in p a ram e te rs of the d igita l se ism ic km 范围 内 2000 ,2008 年 共 131 个 地 震 的 数 据 资 sta tion s a round the ea rthquake s 料 , 对这 些 资料 进行 筛 选 、整 理 , 从 中 挑 选 出 信台号 台站名称 台站代码 台基 地震计型号 噪比高 、尾波持 续 时 间 相 对 长 且 振 幅 不 饱 和 、干 02 TH 02 FB S23A 通海 灰岩 扰小且 S波清晰的 108个地震事件约 423 条地震记 弥勒 灰岩 FB S23A 03 03ML 录参与拟合 Q 值 ,震级为 215 ,317。 c 马龙 砂岩 04 04MA FB S23A 文山 灰岩 21 21W S FB S23A 个旧 灰岩 FB S23A 23 23 GJ 当尾波仅 限于 距台 站 很近 而且 震 级较 小的 地 震时 , 它所反应的介质 特 性不 再依 赖 于震 源和 观 测点位臵 , 而是 其 所 在 区 域 平 均 物 性 的 表 现 。因 此尾波 Q 值主要反映的是以震源和记录台站为焦 c 点的椭球范围内介 质的 品质 因 子特 征 , 地 震尾 波 的研究及其结果可用于分析震源和 地 球内 部结 构 分布特征 。 212 分析处理结果 以弥勒台记录到的 2007年 10月 17 日 16时 43 分 46秒丘北 M 310 地震为例 , 笔者给出了此次地 L 震垂直向的原始波形及该地震不同 中 心频 率的 滤 () 波结果 图 2 。图 3 a为使用 A k i模型以及上述计 ( ) 算步骤后 , 得到的 ln [Af, t3 t ]与 t的关系曲线。 c 用最小二乘法求出 各频 率点 的 拟合 结果 后 , 可 以 ( ) 由 9个频率点的 Q 值拟合得到 Q f。图 3 b给出 c c ( ) 的是弥勒 台 03ML 记录 到 的 56 号 地震 的 尾 波 Q 值拟合结果 。 把某个地震台记录到 的所 有地 震 c 1 2007年丘北 —泸西小震群震中 、台站 、图 的不 同中 心 ( ) 断裂 、云鹏水电站位臵分布图 a及 频率都叠加到一张图上 , 即可 以得 到 该台 站记 录 ( )本研究所用地震 M 2T图 b η到的所有地震的 平 均 Q 值 和 值 。图 4 为 弥 勒 c F ig11 D istribu tion of ep icen te rs、 sta tion s、 fau lts 台记录到的 108 个地震的不同中心频率的尾波 Q cand Yunp eng H yd ropowe r Sta tion 值 , 通过 拟 合 得 到 弥 勒 台 Q 和 f 值 的 关 系 为c ( ) aand the M 2T d iagram of the 0187 ( ) Q f= 8414 f 。 ( )se lec ted ea rthquake s b c 图 2 丘北 M 310地震原始波形及滤波后各中心频率的波形图 L F ig12 P rim a ry wavefo rm of Q iube i M 310 ea rthquake and the wavefo rm L a t d iffe ren t cen tra l frequency po ts afte r filte red ( ) ) 图 3 中心频率为 10 H z的 ln [Af, t3 t]与 c ( ) ( ) t的数据点及回归线 a和地震尾波 Q 值拟合结果 bc ( ) ) ( )F ig13 Po ts and trop ic be tween ln [Af, t3 t] and t when cen tre frequency is 10 H z a c ( )and the ske tch m ap of the re la tion sh ip be tween Q va lue and frequency b c 第 1期赵小艳等 : 云鹏水电站及其周围地震尾波 Q 值研究 c 59 所以我们所测算的 Q 值可能对震群及云鹏水电站 c 所处区域都有效 。 按照上 述 方法 和步 骤 , 我 们 分别 拟 合 得 到 了 (η) 其它 4 个台站的 Q和 值 表 2 , 以及 Q 值与 0 c () 频率 f的关系曲线 图 5 。图 5 为云鹏水电站周 边 5个台站 9个分频的 Q 值取平均值后的分布结 c 果 。从图中可以看出 , 5 个台站各个分频的 Q 值 c 无明显偏离点 , Q 值随 f 值增加而趋于分散 、随 f c 值减小而趋于集中分布 。随着频率的增加 , Q 值 c 差异越来越大 。除了弥勒台外 , 另外 4 个台站的拟 4 弥勒地震台记录的 108个地震图 合曲线几乎重叠在 一起 , 这种 高度 的 一致 性可 能 不同中心频率尾波 Q 值c 是多个地震取平均 值的 结果 , 但也 可 以清 晰地 看 F ig14 Q va lue of 66 ea rthquake s reco rded c 出 , 弥勒台的尾波 Q 值最低 , 弥勒台与震群的距 by M ile d igita l se ism ic sta tion a t d iffe ren t c () 离最近 图 1 , 由于选取的计算 Q 值的流逝时间 cen tre frequency po ts c () 从发震 时 刻 开 始 计 算 比 较 小 , 反 映 的 深 度 较 浅 , 因此 Q 值 很 低 。利 用 表 2 中 各 个 台 站 的 Qc 0 3 计算结果与分析 值 , 可以算出云鹏水电站及其附近区域 Q值的平 0 η ( ) 地震波流 逝时 间的 大小 反 映了 椭球 区 域的 大Q f = Qf 进 行 拟 合 , 得均结 果 。采 用 公 式 c 0 0194 ( ) 小和深度 。根据 Pu lli 1984 的 结 果 , 当 流逝 时 ( ) 到 Q f= 9012 f 。 c ( ) 间为 t秒时 , 所测定的尾波 Q f值反映的是以震 c 源和台站为焦点的椭 球 体介 质的 平 均性 质 。椭 球 体方程为 2 2 2 x y z + + = 112 2 2 2 2 ( ) ( ) ( ) ( ) ( ) vt / 2vt / 2- r/ 2vt / 2- r/ 2 ( )7 式中 r为震群到台站的大致距离 , v是 S波传播速 度 , S 波 在 云 南 地 区 的 传 播 速 度 平 均 值 为 315 km / s。 x、 y 为地表坐标 , z 为深度 坐标 , t为 延迟时间 , 本文中延迟时间为 45 s。针对本文所研 究的地震 , 以弥 勒 台 为 例 , 由 上 式 可 以 得 出 该 台 图 5 不同台站记录的地震尾波 Q 值 c 所记录的 地 震 的 尾 波 Q 值 反 映 的 最 大 面 积 约 为 及其与频率的关系c 2 F ig15 R e la tion be tween Q va lue and frequency 18 000 km, 最大深度为 65 km , 结果主要反映的 c reco rded by d iffe ren t d igita l se ism ic sta tion s 是震区 地 壳 和 上 地 幔 顶 部 的 介 质 品 质 因 子 特 征 , 表 2 云鹏水电站及其附近区域尾波 Q 值计算结果 c Tab12 R e su lts of va lue of Q coda wave a round Yunp eng H yd ropowe r Sta tion c 频率点 拟合结果 台站 地震数 Q η 110 H z 115 H z 215 H z 410 H z 515 H z 810 H z 1010 H z 1210 H z 1810 H z 0 0197 84 89 125 225 400 522 724 966 990 1 275 93 通海 0187 108 73 125 182 355 393 516 587 724 967 84 弥勒 61 81 121 194 360 490 806 951 994 1 197 84 1100 马龙 0196 94 86 129 236 348 562 713 853 1 038 1 229 92 文山 0193 个旧 76 92 160 212 352 479 757 905 1 121 1 203 99 根 据 已 有 的 研 究 成 果 , 姚 安 地 震 的 尾 波有一部分 Q值畸高 , 这可能是由于计算或波形读 0 0195 () 秦嘉政 ,取有误造成 的 , 为 消 除 这 些 畸 异 点 的 影 响 , 我 们武定地震2004 , Q = 49 f 李白基 , c 0195 (δ) ( 采用 1 倍均方差的方法 , 取出总数据 95 %臵宁洱余震 ) 的尾波 Q = 59 f 秦嘉政等 , 2001 , c 0188 信度内的数据 。表 3给出了弥勒台记录到的蓄水前 () Q = 53 f 李琼等 , 2007 。与这些 序列的尾波 c 后 108 个地震的 215 H z、 8 H z与 12 H z频率点的尾值相比 , 云鹏水电站及附近地区的 Q 值较高 , 表 c Q 值的变化情况 , 表中 Q 为所 有 地震 的尾 波波 明云鹏水电站所处的丘北地区的构 造活 动 水平 比 c c 姚安 、武定等地 要 低 。历 史 地 震 活 动 与 这 一 结 论 Q 均 值 , Q表 示 水 库 蓄 水 前 即 2000 年 1 月 至 c 1 类似 ———姚安 、武定 地 区 历 史 上 都 有 6 级 地 震 活 2007年 2 月的平均尾波 Q 值 , Q表示蓄水后即 c 2 ? 动 , 而丘北地区历史上 有地 震 记录 以来 最 大的 一2007年 3 月至 2008 年 12 月的平均尾波 Q 值 , Q c 次地震仅为 2000年 1月 27日丘北 M 515地震 。根 ? S ()表示尾波 Q 值的变化幅度 , 即 Q =Q- Q/Q。 c 2 1 ( ) ( ) 据 Singn 和 H e rrm ann 1983 、W oodgo ld 1990 、 ( ) ( ) H a rtze ll等 1996 和秦嘉政等 2001 给出的国 内外不同地区不同中心频率的地震尾波衰减公式 , 我们绘制了图 6。从图中可以看到 , 丘北震群地区 的尾波衰减与构造地震活动相对较 弱的 美 国中 东 部地区和我国河北昌黎地区的尾波衰减明显不同 , 与构造地震活动频繁 、强烈 的 美国 加州 斯 通峡 谷 地区和我国 云 南 武 定 、施 甸 、宁 洱 等 地 区 也 不 尽 相同 , 是介于两 者 之 间 的 一 种 活 动 状 况 , 应 该 是 一种中强构造活动的表现 。 图 7 弥勒台记录到的地震尾波 Q值 0 随时间变化趋势图 ( ) ( ) a蓄水前 ; b蓄水后 F ig17 Coda Qva lue va ria tion w ith tim e reco rded 0 by M ile d igita l se ism ic sta tion ( ) abefo re the wa te r sto rage; ( ) bafte r the wa te r sto rage 表 3 弥勒台地震尾波 Q 值蓄水前后变化特征 c Tab13 Change of Coda Q va lue reco rded by M ile d igita l c ?美国中东部地区 ; ?美国西部地区 ; ?日本关东地区 ; se ism ic sta tion befo re and afte r impound wa te r ?日本大石山地区 ; ?美国加州斯通峡谷地区 ; ?云南 武定地区 ; ?美国加州帝王谷地区 ; ?北京昌黎地区 ; 频率点 ?云南施甸地区 ; ?云南宁洱地区 ; 参数 215 H z 810 H z 1210 H z λϖ?云南丘北地区 169111 155 150 501647Q c- 1 171103 537125 684156 Q 1 图 6 不同地区的尾波 Q 随频率的变化特征 c - 1 167104 494128 636183 Q 2F ig16 Change cha rac te ristic of the Q va lue and c ? - 2136 % - 8157 % - 7137 % Q frequency among d iffe ren t a rea s 从图 7和表 3 可以看出 , 云鹏水电站蓄水后 ,由上述方 法得 到某 台站 记 录到 的某 次 地震 事 108 个地震的尾波 Q值及分频尾弥勒台记录到的 件的尾波 Q值 , 就 可 以 对 它 进 行 时 间 序 列 作 分 0 0 析 。图 7分别给出了云鹏水电站蓄水前和蓄水后地 波 Q 值有所下降 , 表明水库蓄水对水电站库区及 c 震尾波 Q值随时间变化趋势图 。从图中可以看出 0 其周边区域 的 地 下 介 质 有 影 响 。研 究 表 明 , 在 水 第 1期赵小艳等 : 云鹏水电站及其周围地震尾波 Q 值研究 c 61 黄才中 ,葛焕称 119951中国东部的尾波衰减特征和 Q 值估计 [ J ]. c 库蓄水初期 , 短 时 间 内 水 库 水 位 快 速 升 高 , 达 到 ( ) 地震学报 , 17 2 : 196 - 2021设计蓄水位 。由 于 蓄 水 快 速 增 长 , 使 得 库 区 及 附 李白基 ,秦嘉政 120041云南姚安地区的尾波衰减 [ J ]. 地震学报 , 26 近地区岩体荷载迅速 增 长 。同 时由 于水 的 不断 渗 ( ) 1 : 47 - 511 透 , 使得岩体孔 隙 压 力 不 断 增 大 , 导 致 库 区 及 附 李琼 ,秦嘉政 ,钱晓东 1200712007 年宁洱 614 级地震余震序列尾波 近地区岩体水文地质 条 件产 生变 化 。但 从 蓄水 前 ( ) Q 值研究 [ J ]. 地震研究 , 30 4 : 337 - 3431 c 李琼 ,钱晓东 ,秦嘉政 1200812005 年云南文山地震余震序列尾波 Q 后的尾波 Q值及分频尾波 Q 值的变化幅度以及 c0 c ( ) 值研究 [ J ]. 地震地磁观测与研究 , 29 3 : 2 2 - 281 Q 值所反映的韧性层应力积累状况来看 , 水库蓄 c 马宏生 ,刘杰 ,张国民 ,等 120061云南地区尾波 Q 值的分布特征及 c 水对区域介质影响的 幅 度较 小 , 作 用于 此 区域 的 ( ) 其初步解释 [ J ]. 地震 , 26 3 : 37 - 431应力没有显著调整的 迹 象 , 故 水库 蓄水 短 时期 内 钱晓东 ,李白基 ,秦嘉政 1200412001 年云南施甸 519 级地震余震序 不足以改变该区地下 介 质的 原始 状 态 。致 于水 库 ( ) 列尾波 Q 值研究 [ J ]. 地震地磁观测与研究 , 25 1 : 9 - 171 c 蓄水长时间段有何影响 , 还需进一步跟踪研究 。 秦嘉政 ,阚荣举 119861用近震尾波估算昆明及其周围地区的 Q 值 c ( ) 及地震矩 [ J ]. 地球物理学报 , 29 2 : 145 - 1551 秦嘉政 ,李白基 ,钱晓东 ,等 120011云南武定地震余震近场尾波衰减 4 结论与讨论 ( ) 的区域特征 [ J ]. 地震学报 , 23 1 : 87 - 971 苏有锦 ,刘祖荫 ,蔡民军 ,等 119991云南地区强震分布的深部地球介 ( ) 质背景 [ J ]. 地震学报 , 21 3 : 313 - 3221 本文采 用 A k i单 次 散 射 模 型 , 使 用 高 斯 滤 波 器 , 利用云鹏水电站周边 5 个数字地震台 记录到的王伟君 ,刘杰 1200411999 年岫岩地震序列尾波 Q 的变化过程 [ J ]. c ( ) 地震 , 24 4 , 37 - 441108次地震事件的 423 个地震记录尾波 Q 值进行 c 王伟君 ,刘杰 ,陈凌 120061云南地区尾波衰减特征研究 [ J ]. 地震学 了计算 。结果表明 , 云鹏水电站及其周边地区 Q c ( ) 报 , 28 1 : 8 - 191 0194 徐彦 120021澄江地震的后续丘北 515 级地震研究 [ J ]. 西北地震学 ( ) 值与频率的关系为 Q f= 9012 f 。尾波 Q 值与 c c ( ) 报 , 24 1 : 92 - 961 地震活动 性 呈 负 相 关 关 系 。在 地 震 活 动 密 度 大 、赵小艳 ,苏有锦 ,付虹 ,等 120071欧亚地震带现代构造应力场及其分 强度 高 的 地 区 , Q 值 都 较 低 , 反 之 , Q 值 就 偏 ( ) 区特征 [ J ]. 地震研究 , 30 2 : 146 - 1511 c c A k i K119691A na lysis of the se ism ic coda of loca l ea rthquake s a s sca tte r 高 。本文计算结果表明 云鹏 水 电站 所处 区 域为 中 wave s[ J ]. JGR , 74: 615 - 6311 等构造活动地区 。该区 历来 为 云南 强震 的 弱活 动 A k i K, Chou t B1O rigin of coda wave s: sou rce, a ttenua tion and sca tte ring 地区 , 1900 年以来仅发生了 4 次 5 级地震 , 但自 effec ts[ J ]. JGR , 80: 3322 - 33421 2000 年 以 后 , 该 区 强 震 活 动 有 所 加 强 , 发 生 了H a rtze ll S, L eed s A , F ranke l A , et a l119961Site respon se fo r u rban Lo s 2000年 1 月 27 日丘北 M 515、 2005 年 8 月 13 日 S A nge le s u sing afte rshock s of the No rth ridge earthquake [ J ]. B SSA , 文山 M 513 地震 , 因此本文结果对 研究 云 南强 震 86: 168 - 1921 S Pu lli J J119841A ttenua tion of Coda W ave s in N ew England [ J ]. B SSA , 74 的弱活动地区的介质状态有一定意义 。 ( ) 4 : 1149 - 11661 对云鹏水电站蓄水前后尾波 Q值及分频尾波 0 Singn S, H e rrm ann R B119831R egiona liza tion of cru sta l coda Q in the Q 值的变化幅度作对比研究发现 , 蓄水后尾波 Q c 0 con tinen ta l U n ited State s[ J ]. JGR , 881527 - 5381 值及分频尾波 Q 值有所下降 , 但下降的幅度不是 c Su teau A M , W h itcom b J H119791A loca l ea rthquake coda m agn itude and its re la tion to du ration, mom en t M and loca l m agn itude M [ J ]. 很大 , 表明水库蓄水对 水电 站 库区 及其 周 边区 域 0 L ( ) B SSA , 69 2 : 353 - 3681 的地下介质 有 所 影 响 , 但 影 响 程 度 有 限 。关 于 水 W oodgo ld C D R119901E stim ation of Q in ea ste rn Canada u sing coda 目库蓄水对地震尾波 Q 值的影响及其作用机制 , c wave s[ J ]. B SSA , 881411 - 4291 前国内外研究甚少 , 有待于进一步研究 。X ie J , N u ttli O W 119881 In te rp reta tion of h igh2frequency coda at la rge d istance s: stocha stic mode ling and m e thod of inve rsion [ J ]. Geo2 参考文献 : ( ) p hysica l J1 In t, 95 3 : 579 - 5951 安艺敬一 119911 构 造 活 动 性 的 地 球 物 理 标 志 ———尾 波 品 质 因 素( ) [ J ]. 魏淳 ,译 1世界地震译丛 , 3 : 68 - 871 Re sea rch of C oda Q 2Va lue of Y un pen g Hydropower c S ta t ion an d Its A d ja cen t Reg ion s ZHA O X iao 2yan, FU H ong, M A O Y u 2p ing, W U C heng 2dong ( )Ea r thqua ke A dm in is tra tion of Yunna n P rovince, Kunm ing 650224, Yunna n, C h ina A b stra c t The coda 2w ave 2a ttenua tion 2qua lity fac to r Q in Y unp eng H yd rop ow e r S ta tion a rea and its ad jacen t reg ions is 0 es tim a ted us ing A k iπs s ing le 2sca tte ring a ttenua tion m ode l based on 108 loca l se ism ic even ts reco rded by 5 reg iona l d ig ita l se ism ic s ta tions f rom J anua ry 2000 to D ecem be r 2008. The resu lt show s tha t the re la tionsh ip be tw een Q 0 0194 ( ) 2 C om p a red w ith the ones in o the r a reas of Y unnan w ith s trong tecand f can be desc ribed as Q f= 9012 f . c ton ic se ism ic ac tiv ity, the va lues a re h ighe r in th is a rea w ith m ode ra te tec ton ic ac tiv ity. The Q and Q va lues 0 c dec lined a little af te r the Y unp eng rese rvo ir im p ounded w a te r; th is m eans tha t the w a te r s to rage has inf luenced the m edea cha rac te r in th is a rea, bu t the effec t is lim ited. Key word s: Y unp eng H yd rop ow e r S ta tion; Q 2va lue of coda; A k iπs s ing le sca tte ring a ttenua tion m ode l c