水库诱发地震预测方法及在街面水库中的应用

收稿日期3

水库诱发地震预测方法及在街面水库中的应用

王章森

(福建省水利水电勘测设计研究院,福建福州　350001)

摘要:该文归纳了水库诱发地震的成因类型和判别标志以及水库诱发地震强度预测方法中的综合影响参数法、类比法、地震学法、统计分析法,并结合在街面水库中的应用加深对预测方法的理解。

关键词:水库诱发地震;预测方法;最大地震震级

中图分类号:P315175 文献标识码:B 文章编号:1002-3011(2008)01-0039-03

前言

水库诱发地震是指因水库蓄水而诱使坝区、水库库盆或近岸范围内发生地震。这是一种特殊的地震活动,其震中通常紧邻重要水工设施,特别是中强水库诱发地震多发生在库坝附近的深水库区及其周边,震源浅,震中烈度高,破坏性大,对大坝、发电设施和下游安全威胁大。

水库诱发地震预测是在查明库坝区地质条件的基础上,预测诱发地震的可能性和强度(最大震级),是前期地质勘察工作的一项重要内容。可能性预测一般采用判别标志法即可,强度预测方法则多种多样。笔者归纳、总结了水库诱发地震的成因类型和判别标志,以及在街面水库诱发地震预测研究中所采用的几种较适用的预测方法,供参考。

1　水库诱发地震的成因类型和判别标志111　水库诱发地震的成因类型

水库地震诱发机制相当复杂,有多方面的因素。一般认为库坝区的地质条件是决定一个水库是否发生诱发地震及强度大小的基本条件,库水的重力荷载作用和孔隙压力作用是诱震因素。

水库诱发地震具有不同的成因类型,按工程地质条件分类,主要有构造破裂型、岩溶塌陷型和地壳表层卸荷型三种类型。构造型水库地震有可能达到中等以上强度,破坏性水库地震绝大部分属于构造型水库地震;岩溶塌陷型水库地震只出现在碳酸盐岩分布的库段,与岩溶洞穴和地下管道系统的发育有关,多为弱震或中强震;在断裂发育、坚硬脆性的岩体中,具备一定的卸荷应力和水动力条件时即可发生地壳表层卸荷型水库地震,但其震级小,不会对大坝和周围环境造成危害,因此一般不作过多的研究。

112　水库诱发地震的判别标志

(1)构造破裂型水库诱发地震的判别标志

①区域性断裂或地区性断裂通过库坝区。②新构造断裂通过库坝区,其力学性质为张性或扭性;虽未发现新构造断裂,但深部存在重力梯级异常或磁异常带;断层有晚更新世以来活动的直接地质证据。

③库区或断层沿线有可靠的历史地震记载或仪器记录的地震活动。

④断裂带和破碎带有一定规模的导水能力,有可能成为通往地质体深处的水文地质结构面;地湿梯度较高或地温异常(如温泉出露),可促使库水的渗流活动加速进行。

⑤断裂带或通过旁侧断层、横断层与库水保持一定的水力联系。

⑥断层发育、有利于应力集中(如断裂交接、尖灭部位)和具较强透水性的坚硬性脆岩体(如花岗岩、致密火山岩)。

(2)岩溶塌陷型水库诱发地震的判别标志

①库区有大面积碳酸盐岩分布。

②现代喀斯特作用强烈,可见明显喀斯特管道系统,水库蓄水前已存在天然喀斯特塌陷或喀斯特地震记载。

③合适的喀斯特水文地质结构条件。

当库坝区地质条件符合以上条件越齐备,越典型,水库存在诱发地震的可能性就越大。据资料统计,坝高大于

100m 、库容大于5亿m 3的水库存在诱发地震的可能性大。2　水库诱发地震强度预测方法211　综合影响参数法

该方法把最大库深D 、库容V 和水库面积S 作为影响水库诱发地震的3个因素,它们之间的比例常数E 称为综合影响参数,E =DS/V 。通过统计21个M ≥415级的水库诱发地震震例,得出统计回归线性方程:

M s =-41725+11196E +112421nD (1)水库诱发地震最大震级:

M Max =M s ±σM

(2)

式中σM 为回归标准偏差值,σM =15。

该方法是根据水库规模和大量震例得出的,易操作、简

3:2007-12-0:0119

便,但也将水库诱发地震抽象化了。对于构造破裂型水库地震,该方法预测的可信度较高。

212　类比法

包括地震地质条件类比与震例区工程地质类型类比,前者是指水库所处的构造单元和地震构造特征与临近其它地区的构造单元或地震构造特征相类似时,且有较丰富的地震资料,可取历史上最大地震或预测今后可能的最大地震作为类比区未来可能发生的最大地震;后者是已有的水库地震震例工程地质特征与类比区进行比较,当特征相近时,以震例中的最大震级作为类比区未来可能发生的最大震级。

类比法的前提是两个水库要具有可比性,这种可比性主要表现在以下几个方面:水库规模相近;水库所处的构造单元和地震构造特征相当,或者被类比水库地震地质条件相对较差;被类比水库需蓄水运行多年,已发生过水库诱发地震,并且有完整的地震监测资料。

213　地震学法

以水库地区历史上已知最大地震作为未来水库诱发地震的震级上限,并将这个最大地震分配到最具条件的诱震区,估计最大震级。构造破裂型水库诱发地震的发震危险期在水库蓄水后10a～20a之内。按照构造应变能提前释放的概念,构造破裂型水库诱发地震的最大震级应小于或等于目前至蓄水后10a～20a的时段中天然地震的最大震级。

214　统计预测法

统计预测法有多种,国际上多采用5因子方案,我国有些工程采用了6因子方案和7因子方案,甚至还采用8因子方案等。下面仅对5因子方案做详细说明。

美国人佩克(D1R1Packer)和比切尔(G1B1Baec her)等人通过震例分析,认为水库地震的产生主要与库深、库容、构造应力环境、断层活动性及诱震区介质条件等5个因子有关。该方法的具体内容如下:

(1)根据库深(D)、库容(V)、应力状态(S)、断层活动性(F)以及介质条件(G)等5个因子进行表1的三种因子状态分类;根据世界上39座发震的和173座未发震的大型水库的资料进行统计,得出相应各因子在三种因子状态下的似然率比值如表2。

(2)在各因子相互独立的情况下,根据贝叶斯定理,发

表1　诱震因子及状态

诱震因子

因子状态

123

库深(D)

(m)d1:很深

d1>150

d2:深

150>d2>92

d3:浅

d3<92

库容(V) (×109m3)V1:很大

V1>10

V2:大

10>V2>112

V3:小

V3<112

应力状态(S)S1:逆断层S2:正断层S3:走滑断层断层活动性(F)F活动的F不活动的

介质条件(G)沉积岩变质岩3火成岩

表2　世界发震的与未发震的大型水库似然率比值

诱震因素

似然率比值

状态1状态2状态3 D213160177501921

V111111191601621

S019860155811707

F1150011000

G111450186011008

震和不发震的因素组合的条件概率与发震的和不发震的似然比值的关系式表示如下:

P(RIS|d.v.s,f,g)

P(RIS|d.v.s,f,g)

=

P(RIS)

P(RIS)

LR(d,v,s,f,g)(3) P(R IS|d,v,s,f,g)=1-P(R IS|d,v,s,f,g)(4)式中:P(RIS|d,v,s,f,g)为发生水库地震的条件概率;

P(RIS|d,v,s,f,g)为不发生水库地震的条件概率;

L R(d,v,s,f,g)为发震与未发震的似然率比;

P(RIS)为发生水库地震检验前概率,P(R IS)=01184;

P(R IS)为不发生水库地震检验前概率,P(R IS)=

01816。

(3)根据被预测水库的五个因子,按表1进行因子状态分类,查表2得似然率比值,利用(3)式和(4)式计算出P(RIS|d,v,s,f,g)值。当P(RIS|d,v,s,f,g)≥Pc(临界概率值取012)时,可以认为该水库可能会发生诱发地震。综合模糊评判,依据诱震因素状态的诱震概率有偏差和诱震边界是模糊的,把诱震地震强度评估分为3个等级,即M大(M s ≥5)、M小(M s<5〉和M0(不发震),模糊估计可能性最大的那一等级。

除以上几种方法外,水库诱发地震强度预测方法还有应变能法、有限元法、经验公式法、灰色预测法等。不同的预测方法各有其优缺点,采用多种方法预测水库诱发地震强度很有必要。当多种方法预测出水库诱发地震强度后,考虑地震基本烈度和水库规模,并对预测成果进行综合分析研究,赋予各方法以权重或重视专家意见,评估出最大诱发地震强度。

3　街面水库诱发地震预测

311　水库概况

街面水库面积S=36175km2,库容V=18124亿m3,坝高126m,最大库深D=120m。

在坝址近场25km范围内(含库区)发生过210级～219级地震11次,未记载到3级以上地震。库区地层主要为梨山组、漳平组、长林组的沉积岩。

地质构造以断裂为主,这些断裂晚更新世以来没有活动。政和～海丰断裂带的主干断裂未穿过库区。库区规模较大的次级断裂有两条,均以走滑断层为主。但断裂带和破碎带有一定规模的导水能力,有可能成为通往地质体深处的水文地质结构面。

1

:2:

g1:g2:g: 04