0 引言

受印度板块与欧亚板块的碰撞作用，新疆塔里木盆地西南边缘的西昆仑地区构造运动十分强烈(潘家伟，2011)。位于该地区的西昆仑地震带地震活动频繁，1985—2011年共发生MS≥4.0地震近500次，其中4.0≤MS≤5.0地震425次，MS>5.0地震70次。该地区的地震活动及其发震机理研究已取得了不少的进展(朱传庆等，2008；万永革等，2010；程佳等，2014；姜迪迪等，2014)。吴小平等(2009)结合构造应力特点与潮汐应力触发问题，研究了固体潮与该区地震活动的关联性，得出该区的潮汐应力触发效应属于中等水平；李金等(2014)提出新疆地区地震活动存在一定的优势潮汐相位。

当某区域或断层上的构造应力积累接近临界状态时，若有恰当的潮汐应力与该构造应力叠加，总的应力叠加效果超过了触发阈值，就有可能触发地震。已有很多学者对潮汐触发问题进行了研究。丁中一等(1983)选取了1957—1983年中国或中国邻近地区发生的70次较大的地震，研究地震与潮汐应力的关系，发现有61%的地震在一定程度上受到了潮汐应力的触发。孙长青等(2014)把青藏高原东部及其相邻地区划分成5个构造应力场分区，研究潮汐应力触发效应，结果表明该区域潮汐触发效应的强弱与震级大小、震源深度、构造应力场等有关。张晶等(2007)选取了1970—2007年多个7级以上强震序列，分析了潮汐应力的触发效应，发现共有10个强震序列与引潮力相关性比较好，且在某些特殊区域其水平分量的方位角在多个方向上存在优势。潮汐触发效应与地震发展机理研究仍是目前地震研究的热点问题。

本文选取发生在新疆西昆仑地区的2014年于田MS7.3、2015年皮山MS6.5地震及其余震序列为研究对象，首先进行地震目录完备性分析，以确定最小完备震级MC；然后计算主震震源断层面上的潮汐库仑破裂应力值，以此分析2次主震与潮汐应力之间的触发关系。在此基础上，确定序列中每个地震的潮汐相位角，进而运用Schuster统计检验方法分析地震序列与潮汐应力的统计相关性，探讨潮汐应力对于田和皮山这2个余震序列的触发效应，进一步研究新疆地震活动发震机理。

1 研究区构造背景

新疆是地震多发区，地震主要集中在昆仑山、塔里木盆地等构造活动活跃的地区，且呈现条带状分布(陈建波，2008)。新生代以来，西昆仑地区由于受到欧亚板块与印度板块相互作用所造成的挤压应力，该地区的构造线呈现出北西-南东或接近东西方向(潘家伟，2011)。由于这种挤压作用的存在，产生了康西瓦、普鲁、西昆仑山前等多条较大规模的断裂。邹长桥等(2012)认为，位于西昆仑山南边的西昆仑地震活动区，由于地壳与岩石圈上地幔的解耦使得区域构造背景更为复杂。姜迪迪等(2014)在研究中国西部地区的地壳结构特征与强震活动的关系中，发现西昆仑地区地震的震源深度达到了莫霍面深度。该地区特殊的构造背景导致大震、强震较多。2014年2月12日，于田县发生MS7.3地震，震中位于阿什库勒断裂北段(程佳等，2014)；2015年7月3日，皮山县发生MS6.5地震，主震位于西昆仑山前断裂北缘。图1给出了2次地震序列的空间分布、主震的震源机制解、研究区域内的主要断裂分布等情况，其中高程数据来源于USGS公布的全球数字高程模型GTOPO30，其它数据的来源见表1。

图1 研究区构造背景及2014年于田MS7.3地震 序列、2015年皮山MS6.5地震序列的分布图 Fig.1 Tectonic structure in the study area and distribution of 2014 Yutian MS7.3 and 2015 Pishan MS6.5 earthquake sequences

2 地震数据选取

基于西昆仑地区复杂的构造背景和地震空间分布特点，本文采用新疆地震局提供的2014年于田MS7.3和2015年皮山MS6.5地震序列精定位后的震源参数进行分析。这2个序列的主震震级均在6.0级以上，余震序列的个数均在2 000个以上，适合统计分析。表1给出2个地震序列的主要参数和数据来源。

根据上文所说，针对当前的市场营销战略管理工作中存在的几个缺陷和问题进行了探究分析，加强管理方面的工作，可以保证企业获得更好的发展前景。而且，也要对市场营销战略管理工作进行优化升级根据自身的企业发展情况，制定出更加科学合理的市场营销战略方案，并且从不同的视角和途径来保证企业市场营销管理工作的顺利完成，以此提高企业发展能力，这对于企业自身来说也具有积极的现实意义。

张晶，郗钦文，杨林章，等.2007.引潮力与潮汐应力对强震触发的研究[J].地球物理学报，50(2)：448-454.

图2 2015年皮山MS6.5(a)和2014年于田 MS7.3地震序列(b)的目录完备性分析 Fig.2 Analysis of completeness of catalogues of 2015 Pishan MS6.5(a)and 2014 Yutian MS7.3(b) earthquake sequences 表1 2014年于田MS7.3和2015年皮山MS6.5地震序列主震及余震序列震源参数 Tab.1 Source parameters of 2014 Yutian MS7.3 and 2015 Pishan MS6.5 mainshocks and their aftershock sequences

地震事件主震的震源参数余震序列的参数发震时间(北京时间)φN/(°)λE/(°)震级震中位置节面Ⅰ/(°)节面Ⅱ/(°)走向倾角滑动角走向倾角滑动角余震个数φN/(°)λE/(°)时间范围震级范围震源参数来源于田MS7.3地震2014-02-1217:19:5736.1382.517.3于田24190-2233168-180773035.64～37.3080.69～83.942014-02-12～2014-03-290.1～5.8(李金,王琼,2015)皮山MS6.5地震2015-07-0309:07:4537.5378.176.5皮山11523723146897295436.93～38.0377.07～79.552015-07-03～2016-01-030.1～5.0皮山6.5级地震专题①

注：①http：//www.cea.gov.cn/publish/dizhenj/468/553/101527/101538/20150703221340444810340/index.html.

3 研究原理和方法

3.1 潮汐应力计算方法

如果N个地震事件是随机发生的，其潮汐相位角也是随机分布的。当N足够大(N>10)时，和向量的长度大于或等于|D|的概率(Tsuruoka et al，1995)为：

采用上述潮汐应力计算方法，首先计算球坐标系中于田MS7.3和皮山MS6.5主震的理论潮汐应变张量，通过胡克定律转换为应力张量，然后结合表1中列出的主震震源参数，计算出震源断层面上前后一天的潮汐剪应力τ、潮汐正应力σn和潮汐库仑破裂应力Tcfs随时间变化的关系。根据相关研究(房立华等，2015)，2014年于田MS7.3主震的震源断层面为北东向的节面，因此采用表1中节面1的参数进行计算；同理，根据相关研究(张广伟等，2016)，2015年皮山MS6.5主震的震源断层面采用表1中节面1的参数进行计算。通过计算可以得到于田MS7.3和皮山MS6.5主震发震时，震源断层面上的潮汐应力，结果如图3和表2所示。

(1)

式中：Vn为n阶起潮力位；g为地球的平均重力值；φ为P点的地理纬度；r为P点到地心的距离；Hn(r)和Ln(r)为Love数。表述应力与应变关系的胡克定律为：

σij=2μeij+λθδij

当时拖顶到奔子栏还不通公路，仅有的一条山路还特别难走，年少的陈莲曲珠跟着村里的其他尼姑，磨破双脚在寺院和家之间走着。她们需要徒步两三天才能到达。

(2)

其中：λ，μ为拉梅常量；θ为体应变；σij为应力分量；eij为应变分量。

将(1)式代入(2)式后，便可将P点的潮汐应变分量由胡克定律转换为P点的应力张量。在此基础上，根据主震的震源断层面参数，将潮汐应力张量通过矩阵旋转投影到震源断层面上，得到潮汐剪应力τ、潮汐正应力σn；然后根据库仑破裂准则，计算得出主震断层面上的潮汐库仑破裂应力：

(3)

其中：Tcfs为断层面上因太阳和月亮引力影响所产生的库仑破裂应力，本文中μ′取值0.5，代表断层的有效摩擦系数(Harris，1998；解朝娣等，2015)。

3.2 统计检验方法

传统的潮汐相位分析方法是计算诸如Mf、M2等特定周期的潮汐波(解朝娣等，2015)。Schuster(1897)将Rayleigh在1880年提出的时间序列统计检验原理应用到地震活动与日月相关性的统计检验当中。Schuster统计检验(Tanaka et al，2006，2010)建立了潮汐现象随时间变化的函数，并且这个函数代表了作用在地震断层上的潮汐波的总和。该检验方法通过计算地震事件所对应的潮汐相位角，统计整个地震目录中地震事件相位角的分布情况，分析某一时间段或某一地区是否存在优势相位角，进而判断研究区域内的地震活动是否受到了潮汐的触发(Tanaka et al，2006，2010；孙长青等，2014)。具体而言，首先将发震时间转换为国际时，然后采用Gtoic2(Matsumoto et al，2001)按照前文所述的方法计算潮汐应变，并确定该地震事件所对应的潮汐相位角，再将地震序列中每个地震事件对应的潮汐相位角一一标出，最后统计分析潮汐相位角是否在特定角度附近集中。

在Schuster检验中，规定第j个地震事件为单位向量ej(方向由该事件对应的相位角确定，向量的模为1)。然后将N个地震事件的向量相加得到向量D，表示为：

(4)

式中：N为地震事件的总数；θj为第j个地震事件的潮汐相位角。

假设地球内部的物理性质具有球对称性，在前人研究的基础上(蒋骏，张雁滨，1995；Agnew，2007；李智蓉，张晓东，2011；熊子瑶等，2015)，地球内部或者表面某一点P在球坐标下由潮汐作用下产生的应变分量eθθ、eλλ、eθλ为：

(5)

式中：D2即为向量D模的平方；P值为定量描述潮汐效应触发地震活动的指标，P值范围为0～1，P值越小，潮汐对地震的触发效应越明显(李金等，2014；解朝娣等，2015)。前人通常取P<5%作为判断潮汐能否触发地震的标准(Tanaka，2010；解朝娣等，2015)。

3.3 ETAS模型分析方法

由大森公式发展得到的ETAS模型(Ogata，2005，2011；Lei et al，2008，2011)可以分析序列中活动性的异常变化，这对于研究地震活动性和外部触发因素有帮助。ETAS模型假设每个地震都有触发余震的能力，此能力具有震级依赖性(Lei et al，2013)；并给出总的地震发生率λ(t)可以由先前地震触发的次级余震发生率和表示“背景”地震活动的发生率λ0(t)之和来表述：

腌制过程：将新鲜猪肉拌盐腌好，经过3～4 d，待肉充分吸收盐分后，将其串起晾干，再用甘蔗渣烧火熏烤，经过一段时间后，将其放置在干燥阴凉的地方，即可制成腊肉。经过甘蔗渣熏制而成的腊肉，皮黄肉亮、味道鲜美，是别具一格的地方风味食品。

(6)

式中：Mi和MC分别表述第i个地震事件的震级和地震序列的最小完备震级。一般情况下，“背景”地震活动的发生率λ0为常数，不是时间t的函数；但如果外部触发因素存在时，λ0(t)就包括构造背景活动和外部触发地震活动，分析λ0(t)就有助于分析外部触发地震活动。 K0，c，α和P这些模型参数都是通过最小AIC准则检验给出，合适的平滑化参数ne也可以通过最小AIC准则来确定。

4 结果分析

4.1 主震断层面上的潮汐库仑破裂应力计算结果

从表1中看出，当小正方形的边长从1 cm变化到3 cm时，盒子的容积逐渐增大，在小正方形的边长等于3 cm时盒子的容积达到最大，其后随着小正方形边长的增加容积逐渐减小，当边长为10 cm时容积最小为0.而从图1中发现，当小正方形的边长在3～4 cm之间时盒子的容积达到最大.为了进一步计算最大的容积，在小正方形的边长3～4 cm间，可依次取0.5 cm、0.25 cm、…、0.1 cm的间隔值进行计算.下表2与图2是在3～4 cm间取0.1 cm的间隔值时所求得的无盖长方体形盒子的容积V随小正方形边长x的变化情况.

图3 于田MS7.3(a)和皮山MS6.5(b)主震 震源断层面上的潮汐应力随时间变化曲线 Fig.3 The curve of tidal stress on the fault plane of 2014 Yutian MS7.3(a)and 2015 Pishan MS6.5 (b)mainshocks varied with time 表2 于田MS7.3、皮山MS6.5主震断层面上的理论潮汐应力值 Tab.2 The theoretical tidal stresses on the fault planes of the Yutian MS7.3 and Pishan MS6.5 mainshocks

地震事件发震时间(世界时)潮汐应力剪应力/Pa正应力/Pa库仑破裂应力/Pa于田MS7.3地震2014-02-12 09:19:57-380.78-1429.98-1095.77皮山MS6.5地震2015-07-03 01:07:451111.25-803.00709.75

结合表2和图3可见，在于田主震发震时刻，于田主震震源断层面上的潮汐剪应力处于波谷值与波峰值之间的上升阶段，为-380.78 Pa，而潮汐正应力和潮汐库仑破裂应力位于波峰值到波谷值的下降阶段，分别为-1 429.98 Pa、-1 095.77 Pa。在皮山主震发震时刻，皮山主震震源断层面上的潮汐剪应力和潮汐库仑破裂应力处于波峰值到波谷值的上升阶段，分别为1 111.25 Pa、709.75 Pa，潮汐正应力处于波峰值到波谷值的下降阶段，为-803.00 Pa。

4.2 余震序列与潮汐应力相关性统计分析结果

使用表1中列出的震源参数，对2个地震序列中震级大于MC的所有余震事件进行潮汐相位分析，确定各事件的潮汐相位角，计算得到潮汐相位角的频率分布图，于田和皮山主震的震源断层面均采用表1中节面1进行分析，结果如图4所示。图中显示，统计计算得到的P值几乎都小于5%。再结合图3和表2中潮汐库仑破裂应力分析结果，可以发现潮汐应力对皮山和于田地震的余震序列均有触发作用，但是潮汐应力对于田主震的触发效果不明显，因为于田主震断层面上受到的潮汐库仑破裂应力为负值。

VRIO模型最早由杰恩·巴尼提出[12]。在《从内部寻求竞争优势》一文中，巴尼概括了该模型的核心思想：可持续竞争优势不能通过简单地评估环境机会和威胁，然后仅在高机会、低威胁的环境中通过经营业务来创造，它还依赖于独特的资源和能力，企业可把这些资源和能力应用于环境竞争中。他认为，判断企业特定的资源和能力是优势还是劣势需要回答以下四个问题：(1)企业的资源和能力通过开发机会和抵御威胁能否增加价值？(2)有多少竞争企业已经获得了这些有价值的资源和能力？(3)与已经获得资源和能力的企业相比，不具有某些资源和能力的企业是否面临获取它的成本劣势？(4)企业是否被组织起来开发利用其资源和能力？

图4 于田MS7.3(a)和皮山MS6.5(b) 地震序列潮汐相位的频率分布图 Fig.4 Frequency distribution diagram of tidal phase of Yutian MS7.3(a)and Pishan MS6.5(b) earthquake sequences

4.3 ETAS模型计算结果

从图5可以看出，通过最小AIC准则得到的2个序列的ne均为 14。“背景”地震活动率(f.s.)和模型的其他参数如图5所示。结果表明2个余震序列以受迫成分为主(f.s.>80%)，说明于田、皮山余震序列中80%以上的事件是来自于外部因素的影响，诸如潮汐应力和其他影响因素，其余的是大森型自触发地震。于田和皮山地震序列的模型参数α分别为1.069和0.918，小的α值表明了序列对地震震级的弱依赖性，可以从背景地震活动中区分开(Ogata，1999；Lei et al，2008)。

图5 于田MS7.3(a)和皮山MS6.5(b)地震 序列目录的ETAS分析结果 Fig.5 ETAS analysis results of catologs of Yutian MS7.3(a)and Pishan MS6.5(b) earthquake sequences

5 结论与讨论

本文选取新疆西昆仑地区的2014年于田MS7.3、2015年皮山MS6.5两次中强地震及其余震序列进行地震潮汐触发研究。

(1)首先通过计算主震震源断层面上的潮汐剪应力、潮汐正应力，得到了主震震源断层面上的潮汐库仑破裂应力。结果显示：于田主震断层面上的潮汐库仑破裂应力值为负值，而皮山主震断层面上的潮汐库仑破裂应力值为正值。依据地震触发的相关研究结果(Harris，1998)，正的库仑破裂应力变化超过触发阈值时，可能起到地震触发作用；而负的库仑破裂应力变化没有触发作用，甚至可能起到抑制后续地震发生的作用。但潮汐应力的地震触发作用除了与潮汐库仑破裂应力值的正负以及大小(是否超过触发阈值)有关，还与潮汐应力的时间震荡性有关(解朝娣等，2015)；潮汐触发作用的研究不能单一的以潮汐库仑破裂应力的正负作为唯一判断标准。因此，本文结果表明潮汐应力对皮山主震有明显的触发作用，而对于田主震的触发作用不明显。

(2)通过Schuster统计检验方法，计算2个余震序列潮汐库仑破裂应力的相位角分布发现，潮汐应力对皮山地震和于田地震的余震序列均有触发作用。最后结合ETAS模型分析，也发现西昆仑地区的地震活动受到来自于外部因素的触发作用，诸如潮汐应力和其他影响因素，与大森型自触发的余震活动不一致。

房立华，吴建平，王未来，等.2015.2014年新疆于田MS7.3级地震序列重定位[J].地球物理学报，58(3)：802-808.

新疆西昆仑地区由于受到欧亚板块与印度板块相互作用所造成的挤压应力，产生了多条较大规模的断裂，以及一系列的活动断层，而断层交汇处，是流体活动的通道，所以该地区的地热流体活动很丰富的。而流体对潮汐的拉伸作用敏感，可以凸显潮汐触发作用。这可能是该地区潮汐触发作用较为明显的主要原因。

参考文献：

陈建波.2008.新疆地震构造特征研究[D].兰州：中国地震局兰州地震研究所.

程佳，杨文，刘杰，等.2014.2014年2月12日新疆于田MS7.3地震序列及其构造背景研究[J].地震学报，36(3)：350-361.

丁中一，贾晋康，王仁.1983.潮汐应力对地震的触发作用[J].地震学报，5(2)：46-58.

我的伯父当年也是一位英雄，骑马、打枪绝对的好。后来虽然已经五十岁了，但是风采犹存。我们都爱伯父的，伯父从小也就爱我们。诗、词、文章，都是伯父教我们的。翠姨住在我们家里，伯父也很喜欢翠姨。今天早饭已经开好了。催了翠姨几次，翠姨总是不出来。

姜迪迪，江为为，胥颐，等.2014.中国西部地区地壳结构特征与强震活动相关性研究[J].地球物理学报，57(12)：4029-4040.

其中Pi表示并购方选择现金支付的概率。β0表示常数项，β1,β2,…,β8为回归系数。xip(p=1,2,…,8)为解释变量。表1描述了各解释变量的含义。

蒋骏，张雁滨.1995.发震断层面上潮汐有效剪切应力增量的计算与加卸载的确定[J].中国地震，(1)：72-83.

解朝娣，Lei X L，吴小平，等.2015.潮汐应力对2007年MS6.4宁洱地震震源断层成核失稳过程的影响[J].中国科学：地球科学，45(9)：1409-1420.

1.3 神经功能评价 (1)NIHSS：分别在入院时、溶栓24 h及入院7 d时进行评估；总分0～42分，小于1分为趋于正常或正常，1～4分为轻度神经功能缺损，5～15分为中度神经功能缺损，>15分为重度神经功能缺损。(2)改良Rankin量表评分(mRS)：在患者出院时及发病90 d时进行神经功能恢复状况评估，分为5级，其中0～2级为预后良好，3～5级为预后不良及死亡。(3)溶栓结果评价：出血事件、发病90 d死亡率、发病90 d mRS评分。

李金，蒋海昆，桂荣，等.2014.新疆喀什—乌恰交汇区潮汐触发地震活动的统计检验[J].中国地震，30(1)：64-73.

2.1.2 双主动脉弓患儿 在标准的三血管-气管切面叠加彩色多普勒均不能显示正常的“V”字形结构，从胎儿左侧开始向右侧分别为动脉导管、主动脉左弓、气管、主动脉右弓、上腔静脉，而主动脉左弓和主动脉右弓分别位于气管的左侧和右侧，两者包绕气管，气管位于左弓和右弓之间，形成典型“O”形血管环。

李金，王琼.2015.2014年2月12日于田7.3级地震序列震源机制特征分析[J].中国地震，31(1)：110-120.

李智蓉，张晓东.2011.潮汐应力与大震关系研究[J].地震，31(4)：48-57.

潘家伟.2011.西昆仑构造地貌与阿什库勒地区活动构造研究[D].北京：中国地质科学院.

孙长青，阎春恒，吴小平，等.2014.青藏高原东部及邻区地震断层面上的潮汐应力触发效应[J].地球物理学报，57(7)：2054-2064.

万永革，沈正康，盛书中，等.2010.2008年新疆于田7.3级地震对周围断层的影响及其正断层机制的区域构造解释[J].地球物理学报，53(2)：280-289.

示踪剂的加入是根据之前设置好的投料点坐标作为球心，将半径为15 mm的球体所包含区域浓度值定义为1，混合室内其他区域的浓度定义为0，利用Fluent中初始化(initialization)打补丁的功能(Patch)以补丁方式加入到计算过程中。

吴小平，冒蔚，黄雍，等.2009.基于不同构造分区中国地震的潮汐应力触发效应及相关天文特征[J].中国科学：物理学 力学 天文学，39(6)：901-912.

熊子瑶，解朝娣，雷兴林，等.2015.云南大姚、姚安及附近地区地震活动与潮汐应力的触发关系[J].地震研究，38(1)：66-74.

张广伟，张洪艳，孙长青.2016.2015年新疆皮山MS6.5地震震源机制及余震序列定位[J].地震地质，38(3)：711-720.

为使潮汐相关性的统计分析结果更准确，本文首先对所选地震序列目录的完备性进行分析，结果如图2所示。分析得到皮山地震序列的MC=2.0，序列中M≥2.0地震共有1 088个，b值约为0.53；于田地震序列的MC=2.0，M≥2.0地震共1 297个，b值约为0.50。

朱传庆，杨书江，李同彬，等.2008.塔里木盆地主要断裂的地球物理特征及与天然地震的关系[J].地球物理学进展，23(6)：1685-1691.

邹长桥，贺日政，张智.2012.藏北高原地震活动性特征及其大地构造意义[J].地球物理学进展，27(2)：429-440.

女性在月经前反复出现一系列精神、行为及体质等方面的症状，月经来潮后症状消失的一类病症，称为“经前综合征”，本病在精神、情绪障碍方面更为突出，平素容易紧张，这一部分女性在分娩后发生抑郁的可能性也较大。

Agnew D C.2007.Earth Tides，in Treatise on Geophysics：Geodesy，T A Herring，ed[M].New York：Elsevier，163-195.

Harris R A.1998.Introduction to Special Section：Stress Triggers，Stress Shadows，and Implications for Seismic Hazard[J].J Geophys Res，103(B10)：24347-24358.

Lei X，Shengli M A，Jinrong S U，et al.2013.Inelastic triggering of the 2013 MW6.6 Lushan earthquake by the 2008 MW7.9 Wenchuan earthquake[J].Seismology And Geology，35(2)：411-422.

目前所处的时代还属于复合图书馆时代，纸质文献与电子文献能够形成互补关系。另外，从人力资源角度和读者意见来分析，人员素质和所处环境（或相应条件）决定了许多高级服务开展不起来或开展的效果并不理想，读者对图书馆的要求也都集中在传统服务上，所以我认为有必要在国内图书馆界对传统服务进行升级，其重要途径之一，就是引入新一代信息技术。

Lei X，Xie C，Fu B.2011.Remotely triggered seismicity in Yunnan，southwestern China，following the 2004 MW9.3 Sumatra earthquake[J].J Geophys Res，116(B8)：114-123.

Lei X，Yu G，Ma S，et al.2008.Earthquakes induced by water injection at ～3 km depth within the Rongchang gas field，Chongqing，China[J].J Geophys Res，113(B10：1-12)，doi：10.1029/2008JB 005604.

SSM（Spring+SpringMVC+Mybatis），是目前较为主流的企业级架构方案。标准的MVC设计模式，将整个系统划分为显示层、Controller层、Service层、Dao层四层，使用SpringMVC负责请求的转发和视图管理，Spring实现业务对象管理, MyBatis作为数据对象持久化引擎。

Matsumoto K，Takanezawa T，Ooe M.2001.GOTIC2：A program for computation of oceanic tidal loading effect[J].J Geod Soc Japan，47(1)，243-248.

Ogata Y.1999.Seismicity Analysis through Point-process Modeling：A Review[J].Pure and Applied Geophysics，155(2)：471-507.

Ogata Y.2005.Detection of anomalous seismicity as a stress change sensor[J].J Geophys Res，110(B5)：767-782.

Ogata Y.2011.Pre-seismic anomalies in seismicity and crustal deformation：case studies of the 2007 Noto Hanto earthquake of M6.9 and the 2007 Chuetsu-oki earthquake of M6.8 after the 2004 Chuetsu earthquake of M6.8[J].Geophys J Int，186(1)：331-348.

Schuster A.1897.On Lunar and Solar Periodicities of Earthquakes[J].Nature，61(1449)：321-321.

Tanaka S，Sato H，Matsumura S，et al.2006.Tidal triggering of earthquakes in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region，Japan[J].Tectonophysics，417(1-2)：69-80.

Tanaka S.2010.Tidal triggering of earthquakes precursory to the recent Sumatra megathrust earthquakes of 26 December 2004(MW9.0)，28 March 2005(MW8.6)，and 12 September 2007(MW8.5)[J].Geophys Res Lett，37(37)：195-205.

Tsuruoka H，Ohtake M，Sato H.1995.Statistical test of the tidal triggering of earthquakes：contribution of the ocean tide loading effect[J].Geophys J Int，122(1)：183-194.