用块体应变能预测强震位置与震级研究.docVIP

用块体应变能预测强震位置与震级研究 　　摘要：综合地质构造、岩体力学性质、应力应变状态及其演化历史等各种因素将其统归为应变能积累特征对强震活动进行研究。对利用形变空问特征变化判断强震震源区位置的方法进行了讨论并确认其可行性，提出估算块体应变能的方法。由新构造划分块体边界、地壳测深确定底面进而确定块体体积、形变测量确定年平均应变速率、地壳速度构造确定弹性常数、历史强震确定能量积累起点（起始时问），由此估算块体积累总应变能，进而依据能量估算地震震级。最后以实例论述应变能积累特点对强震位置和震级预测的良好效果以及对地震安全一性评价的重要意义。 　　关键词：活动块体；应变能；巴颜喀拉块体；川滇菱形块体；强震危险区；汶川地震 　　中图分类号：P315. 32 文献标识码：A 文章编号：1000-0666（2015）Ol-0105-09 　　O 引言 　　地震安全性评价的一个重要内容是确定潜在震源区的位置及最大震级，这直接影响抗震设防要求的确定。设防过高，将使工程造价猛增；设防偏低，将会付出惨痛的代价。《工程场地地震安全性评价》（GB17741-2005）为工程场地地震安全性评价制定了一系列原则、标准和方法，其中规定了如何利用构造地质、大地形变测量、地球物理场以及历史地震等进行研究。以往研究多分别用地应力或应变（形变）进行地震预测研究，很少有人进行综合考虑。理论及实验研究显示，潜在震源区的位置和最大震级不仅与构造和历史地震有关，更与应力状态，包括方向、大小、增加速率、集中位置等，以及岩体自身性质，包括强度、杨氏模量、应变大小及增长速率、裂缝（断层）大小和数量等有关，同时还与作用方式有关。许昭永等（2010）综合各种因素，将其统归为一个物理量――应变能积累释放特征进行初步研究。陈祖安等（2009）以应变能密度的空间分布讨论2008年汶川8.0级地震孕震机理。本文将进一步以实例论述如何简便准确地计算活动块体应变能，并将其应用于潜在强震区及其最大震级的预测研究。 　　1 应变能密度及破裂判据 　　依据尹祥础（2011）固体力学的论述，块体应变能（积累） 　　1.1 应变能密度 　　1.1.1 完整块体应变能密度 　　对于线弹性体，应变能密度为 　　1.2 破裂判据 　　最简单的判据就是当K≥Kc时，裂纹开始失稳扩展（破裂）。Kc为断裂韧性，是材料的一个重要参数，表示裂纹开始失稳扩展时应力强度因子的临界值。 　　更常用的破裂判据是应变能密度因子和能量释放率，二者其实是一致的，只是表述方式不同。应变能密度因子准则认为，裂纹沿最小应变能密度因子Smin的方向扩展，即沿满足条件 的方向扩展；当Smin达到一定的临界值Sc时，裂纹开始失稳扩展。即 　　对于完整岩石的剪切破坏，一般用库仑准则，其也适用于塑性破坏，或用于判定结构体沿软弱面滑动破坏。 　　参照应变能密度因子判据和最大能量释放率准则，当所积累应变能密度达到或超过块体同一区域材料所能承受的极限（临界）值时，即发生局部破坏或破裂（地震）。 　　2 关于强震位置和震级的预测方法 　　2.1 潜在震源区位置预测 　　发震位置预测，有多种方法可行。如梅世蓉等（2009）、薛艳等（2009）和陈章立等（2009）就以中强地震的时空活动特征确定汶川地震的发震位置，李延兴等（2009）则认为汶川地震发生在剪切应变高梯度变化区。笔者更倾向以现代形变特征区确定强震位置。强震前会出现一系列地形变方面的特征变化，强震就发生在应变集中区或高应变梯度区（或其附近）。但往往有几个应变集中区或高应变梯度区，究竟哪个区先发生强震需另行研究（李延兴，郭良迁，2002）。这方面资料很多，且技术较为成熟（许昭永等，2010），本文不再赘述。 　　利用形变场的空间特征判定强震发生位置直观、简便，因此本文并不强调用应变能密度的空间分布来判定强震位置。 　　2.2 潜在强震震级预测 　　上述几种理论计算并非都能运用于实际观测中。为方便准确的利用（1）～（3）式估算块体应变能，这里再作几项假定： 　　（1）应力及应变速率都为常数，并且应力和应变速率也为线性关系： 　　（2）块体内应变能密度均匀分布。块体体积为V，于是（1）式简化为 　　（3）以单一形式分别计算各部分应变能，而后再迭加。 　　从式（2）、（3）可看出，应变能包括两部分：体积应变能3项，剪切应变能6项。理论上只要测得3个主应变方向的应变，其它6个（3对）剪应变都可由此导出。实际测量时经常给出两个方向的主应变（最小和最大主应变），因此仅能得到一个（一对）最大剪应变。平均两个方向的伸缩应变能，再乘以3得到体积变形总应变能密度。以一个最大剪应变能密度乘以3.5得到剪切总应变能密度。这是因为中间主应变居于最小和最大主应变之间：它若趋近最