利用非均匀细胞自动机模拟震级与应力降关系.pdfVIP

第31卷第2期 地 震 V01．31。No．2 2011年4月 EARTHQUAKE Apr．，2011 利用非均匀细胞自动机模拟 震级与应力降关系。 李 锰，杨 峰，胡伟华 (新疆维吾尔自治区地震局，新疆乌鲁木齐830011) 摘要：为解释实际观测资料震源破裂过程与自相似(分形)模型的不一致性，本文基于观测结 果，构建了由81×81个细胞单元组成二维非均匀断层模型，并通过设计的细胞自动机模拟程 序进行了模拟试验。研究结果表明：断层结构非均匀性是影响孕震过程特征的重要因素，而且 地震强度分布并非简单的自相似。随着断层非匀质性增加，破裂过程出现由相对的脆性破坏向 塑性破坏特征变化的趋势。利用细胞自动机不仅能较好地解释震级一频度关系中的大、小震级 段低头现象，而且也可解释大震级事件具有相对恒定的应力降，得到了与实际观测研究相一致 的结果。 关键词：细胞自动机；非均匀性；震级一频度关系；应力降 中图分类号：P315 文献标识码：A 文章编号：1000—3274(2011)02—0129—08 引言 常被解释为地震断层尺度具有自相似行为的分形模型。该模型认为地震震级和破裂面积之 间存在一个标度律，震级是随破裂面积的增加而增加，同时在一个很宽的震级范围内，地 震数量与震级之间存在良好的线性关系，并且不同震级的地震具有相对恒定的应力降。然 而近年来，随着地震资料的积累和研究的深入，人们逐渐发现了G—R关系式存在的非线性 和分段性问题。李全林等[1]指出，低震级端的频度“低头”现象不一定是地震仪器的漏记， 孙文福等[2]通过对12个独立过程地震序列的分析，认为G-R关系只在分段的震级范围内 成立，归一化后的震级分别近似服从对数正态分布。此外，Seholz[33通过实验研究认为应 力水平是影响b值的主要因素，而Mogi[‘3的实验研究则认为b值主要取决于介质的非均匀 性。蒋海昆等[s]通过岩石破裂过程的声发射规律，得出b值的变化既取决于岩石的匀质度， 也与岩石所处的应力状态和破裂方式有关。 *收稿日期：2010—08—30；修改回日期：2010—09—28 基金项目：国家自然科学基金项目资助 作者简介：李锰(1965-)，男，河南淮滨人，副研究员，2002年获博士学位，主要从事地质、工程地震及非线性理 论在地学中的应用等研究。 130 地 震 31卷 由于地震环境及其孕育、发生过程的复杂性，造成了不同震例之间的可比性较差，同 时无论是天然地震的统计还是实验室岩石破裂实验，由于地震仪器和测试设备的响应能力 的限制，都不可能保证数据资料的完整性和独立性。然而利用数值模拟可以克服上述不利 因素。本文试图通过细胞自动机这一数值模拟方法来研究介质非均匀性对地震应力降和震 级以及G—R关系影响，以便对实际观测研究结果给于合理解释。 1 观测结果 图1分别给出了利用近震源(5,--,20kin)宽频带地震图，在充分考虑了局部地壳效应基 动研究结果[6]。由图1可以看出，对于较大的地震事件M≥3．5，震级增加是由于断层面积 或长度的增加，应力降保持相对恒定，表现出一种自相似过程；而对于小地震事件 (M3．5)，破裂断层以平均500m半径为特征尺度，应力降与地震矩释放(震级)呈显著的 正相关，但断层尺度大体保持不变，这些小地震事件的地震矩释放可达到3个数量级，但 是其震源尺度变化却在3倍以内，这与地震具有自相似过程表现出不一致性。因此，一条 暑井 图1 let本Matsushiro地区震级与断层尺度(半径)和应力降关系(据文献[6]改编) 断层既可以产生大地震也可以产生小地震，这 些行为的差异并不是由断层上岩石的性质简单 控制的