删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于时空影响域的加权地震网络拓扑特性分析

本站小编 Free考研考试/2020-03-23

徐艳杰, 任涛, 齐义
东北大学 软件学院, 辽宁 沈阳 110169
收稿日期:2017-08-21
基金项目:国家自然科学基金资助项目(61473073, 61104074);中央高校基本科研业务费专项资金资助项目(N161702001);辽宁省高校优秀人才计划项目(LJQ2014028)。
作者简介:徐艳杰(1992-), 女, 河南周口人, 东北大学博士研究生;
任涛(1980-), 男, 辽宁沈阳人, 东北大学教授, 博士生导师。

摘要:针对加利福尼亚地震网络, 为了研究网络中节点间的相互影响关系, 首先基于时空影响域, 以平均震级比对边权值进行定义, 从而生成加权地震网络.其次, 选取零模型为参照物, 分析了加权地震网络的拓扑特性.同时对边权值和节点权值的分布进行分析.结果发现:加权地震网络具有无标度和小世界特性, 且节点和边权值都具有幂律分布特性; 节点的权值与其最大震级值存在正相关.本文基于地震数据所构建的加权地震网络更符合实际情况.
关键词:地震网络复杂网络时空影响域无向/有向网络零模型
Topological Characteristic Analysis of Weighted Seismic Network Based on Space-Time Influence Domain
XU Yan-jie, REN Tao, QI Yi
School of Software, Northeastern University, Shenyang 110169, China
Corresponding author: XU Yan-jie, E-mail: Nancy_Xuyanjie@163.com
Abstract: The seismic network of the south California is taken as the research object. In order to study the influence relationship among nodes in seismic network, the weighted seismic network is generated based on the space-time influence domain and the mean magnitude ratio. Then, based on the reference effect of null model, the topological characteristic of weighted seismic network is analyzed. At the same time, the weight distribution of the edge and the node are analyzed. It is found that the weighted seismic network has scale-free and small-world characteristics, and both the node and the edge weight have power-law distribution. In addition, there is a positive correlation between the weight of the node and its maximum magnitude. Based on the seismic data, the weighted seismic network can reflect the actual situation well.
Key words: seismic networkcomplex networkspace-time influence domainundirected/directed networknull model
地震的研究是各国关注的焦点.地球每时每刻都有地方发生震动, 即便是一个小的地震多发带, 每年的地震次数也有几千甚至上万次.地震是影响人类生命和财产的最重要的自然现象之一, 因此, 对地震的研究有着悠久的历史.其中最著名的是描述地震震级和频数关系的Gutenburg-Richter定律[1]和描述余震发生频率比与时间间隔关系的Omori定律[2]这2个地震定律.近年来, 物理学家试图将地震动力学描述为尺度不变过程.Baiesi等[3]将不同的震颤事件定义为网络的节点, 其中一对节点如果它们之间的相关性超过特定阈值则被连接.而另一方面, Abe等[4-11]则认为每对连续的地震事件都是有联系的.
当前, 复杂网络在处理数据方面已经得到各界的认可, 并且已经应用到各个方面.而在分析网络拓扑特性研究中, 零模型的参照物作用不容忽视.从统计学的角度来看, 为了验证“具有H特征的G网络还具有I的特征”这个零假设是否成立, 而生成的与初始网络具有相同规模相同M特征的随机化网络模型就是零模型[12].在日常应用中, 人们往往习惯了将随机网络模型看作参照系, 在分析网络的小世界特征和聚类效应时, 都参照了随机化网络模型.
为构建地震网络, Abe等[5]在2004年提出了一种对地震区域进行网络建模的方法:先将要研究的地震区域划分为一个个的单元格, 如果单元格有地震发生, 则将这个单元格定义为节点, 若在一次地震事件中两个节点都发生了地震, 则定义这两个节点之间有边, 在一次地震事件中, 如果同一个节点有两次震动发生则直接给这个节点添加自环; 在研究加利福尼亚和日本的地震数据后, 通过统计网络节点度分布特性发现这两个区域的地震网络均具有无标度特性.何璇则在Gardner等[13]研究的基础上提出了基于时空影响域[14-15]构建地震网络.Lin等[16]基于加利福尼亚的震级时间序列建立了地震复发网络, 利用网络熵的特点去测定主震.通过上述方法构建的网络是无权网络, 而在实际网络中, 不同的连边所扮演的角色具有不同的重要性, 所以在研究地震网络时, 需要引进与连边相关联的“权重”这个新属性, 进而构建加权地震网络.
本文针对加利福尼亚地震区域的时空相关性, 根据实际统计数据中地震震级比对地震网络的边权值进行了定义, 生成了加权地震网络.然后, 分析了加权地震网络的边权值和节点权值的无标度特征.通过零模型的参照物作用研究了地震网络的无标度特征和小世界特征.最后分析得到节点权值与节点最大震级存在线性关系.
1 加权地震网络加利福尼亚位于环太平洋地震带南部, 发生在此地的地震多为浅源地震, 所以在划分地震区域的时候, 可以只按照经纬度划分, 不考虑深度.本文主震节点默认为震级至少是2.5级, 且文中所有节点的震级都为1级以上.
本文划分节点时采用经纬度划分节点的方法, 将加利福尼亚地震数据所包含的地理区域(114.0 W至122.0 W, 32.0 N至37.0 N)划分为若干个小格子, 除去不属于加利福尼亚州的小格子, 然后将这些地理区域进行编号, 根据南加利福尼亚地震数据中心统计的1992~2014年地震数据(数据来源于http://service.scedc.caltech.edu/ftp/catalogs/SCEC_DC/), 除去不属于加利福尼亚的地震数据和震级为1级以下的小区域, 从而得到相应的网络节点.
对于网络节点i, 它的震级与影响直径和影响持续时间不是线性关系, 它们之间的关系可由式(1), 式(2)[13-15]表示:
(1)
(2)
其中:Mi表示发生在节点i的地震震级; Ti表示节点i发生Mi级地震后的影响时间; Li表示节点i发生Mi地震后的最大影响距离, 即影响直径; a1, a2, b1, b2是根据统计所得的常数.
一次地震事件中, 节点i和节点j之间的时间和空间的联系由式(3)和式(4)来定义:
(3)
(4)
其中:dij为节点i, j之间的距离; ti, tj表示节点i, j的地震时刻; ilat, jlat表示节点i, j发生地震的经度; ilon, jlon表示节点i, j发生地震的维度.默认震源的震级大于其所影响的区域震级.
网络中节点之间连边的定义:在一次地震事件中, 如果节点i, j先后发生地震且节点j在节点i的影响范围和影响时间之内(ΔtTidijLi), 则两个节点之间存在连边(有向网络中i有一条边指向j); 如果节点i, j, v都依次发生了地震且节点v在节点ij的影响时间和影响范围之内, 则节点v与节点ij都有连边(有向网络中ij都有边指向v), 排除重边,同时排除节点受自身影响发生的地震, 即自环.
以往的加权地震网络在定义边权值时只考虑共同地震次数或者地震影响次数, 而忽略了地震区域间震级的影响, 因此, 本文以平均震级比定义边权值, 第m次地震时的边权值sijm定义为
(5)
其中:Mim表示第m次地震时震源节点i的震级; Mjm表示第m次地震时节点j的震级.
1992~2014年的加利福尼亚南部的震级不小于1的地震记录共有228 393次.为了便于研究, 在地震网络中, 选取众多边权值的期望表示加权地震网络的边权值, 即
(6)
其中, ni表示节点i影响节点j所发生的地震次数, 因为是针对具体数据所进行的计算, 所以不会出现数据无穷大的情况.
而在无向网络中, 式(6)可转换为
(7)
2 基于零模型地震网络拓扑特性分析置乱算法是最近经常使用的生成零模型的方法.置乱算法不仅能够进行断开连边重连改变网络的节点连接关系, 还能够在保持原先连接关系的基础之上随机化网络的某些因素.因此, 使用置乱算法生成对应的零模型, 通过各种统计量相对值的对比可以进一步深入理解网络.在各阶零模型中, 因为0阶的零模型与原网络具有相同的节点数目和连边数目, 也就是说只要维持节点和连边的数目一致, 网络其他特性可以随意改变.因此, 0阶零模型一般都用于对比分析网络节点度分布和聚类系数分布特性.
基于随机置乱算法[17]生成相应的0阶零模型, 从而对本文的地震网络分为无向和有向进行拓扑特性研究.因为地震网络边权值相等的很少, 因此在生成0阶零模型时, 选择无权地震网络作为原始网络.
图 1所示, 原始无向地震网络的度分布符合幂律分布, 而其0阶零模型的度分布不具有幂律分布特征, 则可得无向地震网络具有无标度特性, 且计算可得其幂指数为1.217 1.
图 1(Fig. 1)
图 1 双对数坐标下无向地震网络及其0阶零模型节点的概率度分布Fig.1 Nodes probability degree distribution of the undirected seismic network and its 0 order null model in lg-lg coordinate

通过计算可得, 无向地震网络的聚类系数为0.368 7, 它的0阶零模型的聚类系数为0.038 0, 与其节点平均度相同的完全随机网络的聚类系数为0.008 6.因为原始无向地震网络的聚类系数远远大于与它对应的完全随机网络和它的0阶零模型, 从而可以得出无向地震网络具有小世界特性.
对有向地震网络进行分析.如图 2所示, 原始有向地震网络的度分布符合幂律分布, 且出度的幂指数为1.112 9, 入度的幂指数为1.253 6.而与之相比, 其0阶零模型的出度和入度分布都不具有幂律分布特性.从而可得有向地震网络具有无标度特性.
图 2(Fig. 2)
图 2 双对数坐标下原始有向网络及其0阶零模型节点的度分布概率Fig.2 Nodes probability degree distribution of the directed seismic network and its 0 order null model in lg-lg coordinate (a)—入度分布;(b)—出度分布.

通过计算可得, 有向地震网络的聚类系数为0.368 7, 它的0阶零模型的聚类系数为0.075 7, 与其节点平均度相同的完全随机网络的聚类系数为0.008 7.因为原始有向地震网络的聚类系数远远大于与它对应的完全随机网络和它的0阶零模型, 从而可以得出有向地震网络具有小世界特性.
综上所述, 本文以经纬度划分节点, 基于时空影响域构造的地震网络, 无论是无向还是有向网络, 都具有无标度和小世界特性.
3 加权地震网络节点和边权值分布特性分析表 1为基于平均震级比的加权地震网络边权值概率.由表 1可得,无向网络中不存在权值大于等于0.8的连边, 权值大于0.5的连边占总体5%(表 1中边权值为0.1时的概率代表网络中边权值大于0.1小于等于0.2的概率); 而在有向网络中, 不存在权值大于等于1的连边, 权值大于0.5的连边占14%.由此可得, 基于平均震级比的加权地震网络中, 震源节点影响其他节点发生与它震级相近地震的可能性很小.
表 1(Table 1)
表 1 基于平均震级比的加权地震网络边权值概率Table 1 Seismic network edge weight probability based on average magnitude ratio
网络 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
无向网络 0.67 0.15 0.07 0.04 0.02 0.02 0.01 0.01 0.01 0
有向网络 0.51 0.15 0.09 0.06 0.05 0.04 0.03 0.02 0.03 0.02


表 1 基于平均震级比的加权地震网络边权值概率 Table 1 Seismic network edge weight probability based on average magnitude ratio

通过对边权值和其所占概率取双对数, 进而研究边权值是否符合幂律分布特性.如图 3所示, 无论是在无向网络还是有向网络中, 边权值都服从幂律分布, 对于无向网络, 边权值的幂指数为1.387 0.而对于有向网络, 它的幂指数为0.934 1.这说明网络中边权值小的连边居多, 边权值大的连边只是极少数.
图 3(Fig. 3)
图 3 双对数坐标下加权地震网络边权值概率分布图Fig.3 Probability distribution of the weighted seismic networks edge weights in lg-lg coordinate (a)—无向网络;(b)—有向网络.

本文定义无向网络节点的权值为节点边权值的总和, 即.而有向网络节点的权值分为出边权值和入边权值.通过计算可得, 无向和有向地震网络节点权值分布都具有幂律分布特性, 也就是富人俱乐部现象.无向网络节点权值的幂指数为1.085 6, 有向网络根据出边权值和入边权值计算所得的节点的权值幂指数分别为0.928 9和0.980 9.且通过统计可得, 无向加权地震网络的节点最大震级与其权值呈正相关, 而有向加权地震网络的节点最大震级与其基于出边所计算的权值呈正相关, 根据数值仿真可得它们之间的比值分别拟合为(0.103, 3.126)和(0.016, 3.729).这说明节点权值越高, 节点的影响性越大, 即它发生大型地震的可能性越大.
4 结论1) 本文对美国南加利福尼亚州1992~2014年的地震数据进行建模分析.首先, 按照经纬度划分节点, 基于时空影响域构造节点连边, 并根据节点间的平均震级比定义边权值(现有的加权地震网络是以共同地震次数或者地震影响次数作为边权值), 构建了无向和有向加权地震网络.
2) 通过零模型的参照物作用对地震网络的度分布以及聚类系数进行了分析, 发现无向和有向地震网络都具有无标度和小世界特性.
3) 通过对边权值以及节点权值的分析, 发现无向和有向加权地震网络的边权值和节点权值都符合幂律分布特性.
4) 对加权地震网络的节点权值与节点发生过的最大震级进行分析, 发现二者之间存在正相关, 这与实际情况相符, 本文的工作将有助于对地震网络进行研究.
参考文献
[1]Gutenburg B, Richter C F. Frequency of earthquakes in California[J].Bulletin of the Seismological Society of America, 1944, 34(4): 185–188.
[2]Omori F. On the aftershocks of earthquakes[J].Journal of the College of Science, 1894, 7: 111–119.
[3]Baiesi M, Paczuski M. Scale-free networks of earthquakes and aftershocks[J].Physical Review E, 2004, 69(6): 066106–1.DOI:10.1103/PhysRevE.69.066106
[4]Abe S, Suzuki N. Scale-free network of earthquakes[J].Europhysics Letters, 2004, 65(4): 581–586.DOI:10.1209/epl/i2003-10108-1
[5]Abe S, Suzuki N. Small-world structure of earthquake network[J].Physica A:Statistical Mechanics & Its Applications, 2004, 337(1/2): 357–362.
[6]Abe S, Suzuki N. Complex earthquake networks:hierarchical organization and assortative mixing.[J].Physical Review E, 2006, 74(2): 26113–1.DOI:10.1103/PhysRevE.74.026113
[7]Abe S, Suzuki N. Dynamical evolution of clustering in complex network of earthquakes[J].European Physical Journal B, 2007, 59(1): 93–97.DOI:10.1140/epjb/e2007-00259-3
[8]Abe S, Suzuki N. Violation of the scaling relation and non-Markovian nature of earthquake aftershocks[J].Physica A, 2009, 388(9): 1917–1920.DOI:10.1016/j.physa.2009.01.031
[9]Abe S, Suzuki N. Scaling relation for earthquake networks[J].Physica A, 2009, 388(12): 2511–2514.DOI:10.1016/j.physa.2009.02.022
[10]Abe S, Suzuki N. Determination of the scale of coarse graining in earthquake network[J].Europhysics Letters, 2009, 87(4): 48008–48012.DOI:10.1209/0295-5075/87/48008
[11]Abe S, Pastén D, Suzuki N. Finite data-size scaling of clustering in earthquake networks[J].Physica A, 2011, 390(7): 1343–1349.DOI:10.1016/j.physa.2010.11.033
[12]汪小帆, 李翔, 陈关荣. 复杂网络理论及其应用[M]. 北京: 清华大学出版社有限公司, 2006.
( Wang Xiao-fan, Li Xiang, Chen Guan-rong. Complex network theory and its applications[M]. Beijing: Tsinghua University Press Co, Ltd, 2006.)
[13]Gardner J K, Knopoff L. Is the sequence of earthquakes in Southern California, with aftershocks removed, Poissonian?[J].Bulletin of the Seismological Society of America, 1974, 64(5): 1363–1367.
[14]何璇, 赵海, 蔡巍, 等. 基于时空影响域的地震网络构造方法[J].东北大学学报(自然科学版), 2014, 35(10): 1395–1399.
( He Xuan, Zhao Hai, Cai Wei, et al. Construction method based on time-space influence domain for earthquake network[J].Journal of Northeastern University(Natural Science), 2014, 35(10): 1395–1399.)
[15]赵海, 张娅, 何璇, 等. 基于时空影响域的地震网络动力学演化特征分析[J].东北大学学报(自然科学版), 2015, 36(9): 1232–1236.
( Zhao Hai, Zhang Ya, He Xuan, et al. Dynamic evolution analysis of the earthquake network based on the time-space influence domain[J].Journal of Northeastern University(Natural Science), 2015, 36(9): 1232–1236.DOI:10.3969/j.issn.1005-3026.2015.09.004)
[16]Lin M, Fan X X, Wang G, et al. Network structure entropy and its dynamical evolution for recurrence networks from earthquake magnitude time series[J].The European Physical Journal B, 2016, 89(5): 1–7.
[17]尚可可, 许小可. 基于置乱算法的复杂网络零模型构造及其应用-本期"复杂性科学"专栏评述[J].电子科技大学学报, 2014, 43(1): 7–20.
( Shang Ke-ke, Xu Xiao-ke. Construction and application of complex network zero model based on scrambling algorithm-this issue of "complexity science"[J].Journal of Electronic Science and Technology University, 2014, 43(1): 7–20.DOI:10.3969/j.issn.1001-0548.2014.01.002)

相关话题/地震 网络

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 一种用于肺结节恶性度分类的生成对抗网络
    徐久强,洪丽萍,朱宏博,赵海东北大学计算机科学与工程学院,辽宁沈阳110169收稿日期:2017-07-31基金项目:国家科技支撑计划项目(2012BAH82F04)。作者简介:徐久强(1966-),男,辽宁北镇人,东北大学教授;赵海(1959-),男,辽宁沈阳人,东北大学教授,博士生导师。摘要:针 ...
    本站小编 Free考研考试 2020-03-23
  • 基于局域社团的人类脑功能网络生成模型
    司帅宗,刘晓,朱剑,赵海东北大学计算机科学与工程学院,辽宁沈阳110169收稿日期:2017-08-09基金项目:中央高校基本科研业务费重大科技创新项目(02190022117021)。作者简介:司帅宗(1987-),男,河南长垣人,东北大学博士研究生;赵海(1959-),男,辽宁沈阳人,东北大学教 ...
    本站小编 Free考研考试 2020-03-23
  • 基于全卷积网络迁移学习的左心室内膜分割
    齐林1,吕旭阳1,杨本强2,徐礼胜1,31.东北大学中荷生物医学与信息工程学院,辽宁沈阳110169;2.沈阳军区总医院放射科,辽宁沈阳110016;3.东北大学教育部医学影像计算重点实验室,辽宁沈阳110169收稿日期:2017-08-11基金项目:国家自然科学基金资助项目(61773110,61 ...
    本站小编 Free考研考试 2020-03-23
  • 基于全卷积网络的左心室射血分数自动检测
    徐礼胜,张书琪,牛潇,徐阳东北大学中荷生物医学与信息工程学院,辽宁沈阳110169收稿日期:2017-08-02基金项目:国家自然科学基金资助项目(61773110,61374015,61701099);中央高校基本科研业务费专项资金资助项目(N161904002)。作者简介:徐礼胜(1975-), ...
    本站小编 Free考研考试 2020-03-23
  • 遗传算法优化神经网络的热轧带钢弯辊力预报模型
    王振华,龚殿尧,李广焘,张殿华东北大学轧制技术及连轧自动化国家重点实验室,辽宁沈阳110819收稿日期:2017-08-21基金项目:国家重点研发计划项目(2017YFB0304100);国家自然科学基金资助项目(51634002,51774084,51704067)。作者简介:王振华(1990-) ...
    本站小编 Free考研考试 2020-03-23
  • 筒仓-贮料-地基相互作用系统地震响应研究
    李树蔚1,2,赵文1,国志雨31.东北大学资源与土木工程学院,辽宁沈阳110819;2.中国能源建设股份有限公司,北京100022;3.中国能源建设集团黑龙江省电力设计院有限公司,黑龙江哈尔滨150078收稿日期:2018-01-09基金项目:国家自然科学基金资助项目(51504123)。作者简介: ...
    本站小编 Free考研考试 2020-03-23
  • 基于复变自适应神经网络的电网相位估计方法
    李云路,王大志,宁一,回楠木东北大学信息科学与工程学院,辽宁沈阳110819收稿日期:2015-07-09基金项目:国家自然科学基金重点资助项目(61433004);国家自然科学基金资助项目(51467017)。作者简介:李云路(1986-),男,辽宁沈阳人,东北大学博士研究生;王大志(1963-) ...
    本站小编 Free考研考试 2020-03-23
  • 连续型公交专用道网络布局优化模型
    卢小林1,俞洁2,邹难1,刘好德31.山东大学控制科学与工程学院,山东济南250061;2.威斯康辛大学密尔沃基分校土木工程系,威斯康辛州密尔沃基53201;3.交通运输部科学研究院,北京100029收稿日期:2015-08-25基金项目:国家科技支撑计划项目(2014BAG03B04);山东大学基 ...
    本站小编 Free考研考试 2020-03-23
  • 一类求解优化问题的神经网络及其全局吸引性分析
    王占山1,康云云2,牛海莎11.东北大学信息科学与工程学院,辽宁沈阳110819;2.国网池州供电公司,安徽池州247100收稿日期:2015-10-18基金项目:国家自然科学基金资助项目(61473070,61433004);流程工业综合自动化国家重点实验室项目(2013ZCX01)。作者简介:王 ...
    本站小编 Free考研考试 2020-03-23
  • 基于深度置信网络的云应用负载预测方法
    马安香,张长胜,张斌,张晓红东北大学计算机科学与工程学院,辽宁沈阳110169收稿日期:2015-09-16基金项目:国家科技支撑计划项目(2014BAI17B00);国家自然科学基金资助项目(61572116,61572117,61502089)。作者简介:马安香(1979-),女,辽宁沈阳人,东 ...
    本站小编 Free考研考试 2020-03-23