周江存 孙和平 徐建桥 陈晓东

(中国武汉430077中国科学院测量与地球物理研究所,大地测量与地球动力学国家重点实验室)

海潮负荷对地球内部潮汐应力和应变的影响*

(中国武汉430077中国科学院测量与地球物理研究所,大地测量与地球动力学国家重点实验室)

详细介绍了海潮负荷影响的计算理论,基于PREM地球模型计算了地球内部的负荷勒夫数及负荷格林函数； 并以上海台和武汉台为例,计算了海潮负荷对不同深度处的应力和应变潮汐的影响． 结果表明： 深度是影响海潮负荷应力的一个重要因素,在靠近计算点的区域,应力负荷的影响随深度增大而减小； 而对于远离计算点的区域,应力负荷的影响却随深度增大而增大； 另外,深度会影响某些应力和应变潮汐分量时间变化的相位． 在沿海地区,海潮负荷对应力和应变的影响超过了应力和应变固体潮的影响,因此该影响在应力和应变测量中必须要加以考虑．

海潮负荷 应变 应力 深度

引言

随着观测精度的提高,海潮负荷效应作为大地测量观测中的一个改正项显得越来越重要． 目前的研究主要关注在地球表面,而对地球内部的研究甚少． 应变和应力测量作为研究地震的重要手段,研究海潮负荷对地球内部应力和应变的影响就显得非常重要． 海潮负荷是在日月引潮力作用下海水的涨落对地球的影响,通常包括海水的直接引力、 地球在引力作用下的变形,以及由于变形导致的地球内部质量的重新分布(Melchior,1978； 方俊,1985； 许厚泽,2010)等3部分效应； 因此海潮负荷可以认为是日月引潮力对地球的二次效应． 海潮负荷对于在地表或空间进行的大地测量学与地球物理学方面的观测都具有非常大的影响(周旭华,2001； 周江存,孙和平,2005,2007a,b),在沿海地区更是如此(周江存,孙和平,2005)．

对于应变或应力潮汐来说,它们只包含地球在潮水压力作用下地球的变形效应． 应变测量一直是研究地震及其相关科学的重要手段(周坤根等,1991； 蒋骏,张雁滨,1994； 蒋骏等,1995),在应变测量结果中具有非常明显的海潮负荷的信号,因此扣除海潮负荷的影响就显得非常重要． 而地球内部的应力则对于地震的孕育与发生都具有非常重要的作用,在大地震来临之前,地球内部断层附近处于一个临界状态,任何微小的外力作用都有可能诱发地震． 因此,早在19世纪末,就有关于潮汐触发地震的研究(Schuster,1897)． 目前也有更多的证据表明,潮汐应力与地震的发生具有一定的联系(Tsuruokaetal,1995； Tanakaetal,2002； Métivieretal,2009)． 目前关于海潮负荷的研究只限于地球的表面,而关于海潮负荷对地球内部影响的研究则比较少. 孙文科(1991)曾对地球内部的海潮负荷应力潮汐进行了研究,但仅仅给出了2阶负荷勒夫(Love)数的结果； 阎春恒等(2010)研究了海潮负荷对地震断层面上的应力潮汐的影响,但计算应变时采用了差分代替微分的近似算法. 另外,目前应力和应变测量有向地下深部发展的趋势(骆鸣津等,2008)． 因此,研究海潮负荷对地球内部应力和应变的影响具有较大的科学意义,希望通过本文的研究能够对认识海潮负荷如何影响地球内部的应力和应变状态提供参考．

1 计算理论

1.1 负荷勒夫数

对于海潮负荷的研究,通常基于球对称、 非自转、 弹性、 各向同性(spherically symmetric non-rotating elastic isotropic, 简写为SNREI)的地球模型． 通过计算格林函数和负荷质量的积分获得海潮负荷的影响． 其中,负荷勒夫数的计算是基于SNREI地球模型在球坐标系下求解如下微分方程组：

(1)

式中,Y=[y1, y2, y3, y4, y5, y6]T是由径向应力、 径向位移、 水平应力、 水平位移、 重力位及重力组成的向量,A是由地球模型参数组成的矩阵,r是到球心的距离(Longman,1962,1963； Farrell,1972)． 许多研究人员对此进行了研究(汪汉胜等,1996； Sun,Sjöberg,1999； 徐建桥,孙和平,2003),即采用龙格-库塔方法求解并给出了地表的负荷勒夫数． 由于我们需要研究地球内部的应力和应变特征,因此需要获得地球内部的负荷勒夫数,实际上在求解地表负荷勒夫数的过程中其内部的解也是一并求出的． 我们采用Sun和Sjöberg (1999)关于内部负荷勒夫数的定义

(2)

(3)

式中,f是引力常数,ρ是密度． 由此,其地表的边界条件变得更加简洁,为

(4)

在计算过程中,为了防止计算溢出,采用汪汉胜等(1996)的变换方法,并同时对数值积分的每一步进行归一化处理． 另外, 对0阶和1阶计算的特殊处理可参考Sun和Sjöberg (1999)文章．

1.2 应力和应变负荷格林函数

(5)

式中,μ和λ为拉梅(Lamé)常数,由地球模型给出,是r的函数；y1对r的导数可以通过式(1)求得,它是y1,y2和y3以及拉梅常数的函数．Pn是勒让德(Legendre)多项式. 对于应变的负荷格林函数可以通过应力-应变的关系求得,也可通过垂直和水平位移求得,即

(6)

图1 球面坐标示意图Fig.1 Sketch map of spherical coordinate

同样,y3对r的导数也可以通过式(1)求得,它是y1,y3和y4以及拉梅常数的函数．

1.3 海潮负荷计算

上述负荷格林函数是基于位于北极点的负荷的,而实际上海水负荷质量是全球分布的,因此在实际计算时需要将每个单元负荷质量对应力的影响进行坐标变换并进行叠加,即将以负荷质量所处位置为北极的坐标系转换到地球坐标系(图1)． 图1中N表示地球的北极,A表示负荷质量所处位置,B是计算点,则以负荷点为北极的坐标系下的应力张量为

(7)

根据应力的坐标变换公式即可求出在地球坐标系下的应力张量的各个分量,其转换矩阵为

(8)

即以B点处的r轴为旋转轴,逆时针旋转α角度,可得在地球坐标系中6个独立的应力格林函数为

(9)

由式(9)可以看出,除了rr分量的格林函数外,其它分量的格林函数都与方位角有关． 因此,由此计算的海潮负荷应力分量也不像rr分量那样,只与负荷所在位置和计算点之间的角距离有关．

有了负荷格林函数,就可以计算海潮负荷的影响. 其计算公式(Agnew,1997)为

(10)

式中,L是负荷应力,ρw是海水密度,H是海潮的瞬时潮高,G是式(9)表示的格林函数． 计算时,我们采用NAO99b全球海潮模型(Matsumotoetal,2000),并采用积分格林函数方法,其定义(Goad,1980； Agnew,1997)为

(11)

图2 地球地表及内部负荷勒夫数随球谐阶数的变化Fig.2 Changes of surface and internal load Love numbers with harmonic degree

通过求出的y1,y2,y3和y4,以及勒让德函数及其导数,就可以算出积分格林函数(Zhouetal,2013)． 由于格林函数在近区的变化比较剧烈,因此近区的积分负荷格林函数的计算需要非常小的积分步长．

表1 计算积分负荷格林函数的积分步长

2 数值结果

我们首先采用PREM地球模型(Dziewonski, Anderson,1981； 汪汉胜等,1996)进行计算,获得了地球地表及内部的负荷勒夫数(图2)． 由图2可知： 当阶数很大时,地球内部的负荷勒夫数都趋近于0,这与地表的负荷勒夫数有很大的不同； 当深度变大时,负荷勒夫数随着阶数的增大更快地趋近于0,这是因为它们都含有(r/R)n的因子．

(5)对需要调节流量及其他参数的工艺设备，在集中控制中利用变频器进行适时调控，满足生产工艺调控需要，根据需要调整电机输出功率，也起到节电的作用。

通过负荷勒夫数,我们就可以计算获得积分负荷格林函数,在此过程中我们运用了增强级数收敛的方法(Farrell,1972),并且积分步长也采用随角距离而变化的值(表1),在近区采用较小的步长,在远区采用相对大一点的步长．

图3给出了不同深度处的4个独立积分负荷格林函数. 由于分成4种积分步长,函数值变化达几个数量级,因此为了看清各部分的细节,分成4个部分分别显示． 从图3可以看出： 深度越小,近区的积分负荷格林函数值越大,而远区的函数值则越小； 深度越大,积分负荷格林函数最大值出现的地方离负荷点越远,但是其幅度要比深度小的最大值小； 对于ψψ和αα分量,当角距离小于10°时,积分负荷格林函数值对深度比较敏感； 当角距离超过10°时,对于不同深度的积分负荷格林函数值基本保持一致; 而对于rr和rψ分量,无论是在远区还是在近区,函数值对深度都异常敏感．

图3 不同深度处的积分负荷格林函数

由于海潮负荷的影响与计算点到海洋的距离有非常大的关系,因此我们选择了近海的上海台和远离海洋的武汉台分别计算了海潮负荷对应力潮汐的影响(取2013年1月1日为中心的40小时时间段, 开始于2012年12月31日4时),结果如图4和图5所示．

由图4可知： 对于远离海洋的武汉台,海潮负荷对应力潮汐的影响量级在102—103Pa之间,呈现半日波的特征； 总体说来,海潮负荷对浅部的应力影响很小,而对深部的应力影响较大,这一点也可从图3中远区的格林函数看出； 除了λλ分量外,尽管在同一个经纬度,但是在不同深度处,海潮负荷对应力潮汐影响的相位是不一致的,有些深度甚至是反相的,这反映了深度对于海潮负荷应力潮汐计算的重要性．

图5给出了海潮负荷对上海台不同深度处应力潮汐变化的影响. 通过与图4对比可以看出, 除了表现为半日波特征外明显不同： 海潮负荷对应力潮汐的影响量级在103—104Pa之间,比在武汉台的影响大一个数量级；rr分量的应力潮汐振幅随深度增大而减小,这正好与武汉台相反； 3个正应力潮汐的相位随深度的变化而保持不变,而切应力潮汐的相位随深度变化有很大的不同,有些也出现反相的情形． 同样,图6和图7给出了上述两个台海潮负荷对应变潮汐的影响. 由于切应变与应力只相差一个比例系数2μ,因此其变化特征与应力相同,所以图中只给出了正应变的结果． 对比两图可知： 海潮负荷对武汉台的应变潮汐影响的振幅在10-9,而对上海台的影响振幅在10-8—10-7； 对于远离海洋的武汉台，

图4 海潮负荷对武汉台不同深度处的应力潮汐变化的影响

图5 海潮负荷对上海台不同深度处的应力潮汐变化的影响

图6 海潮负荷对武汉台不同深度 处的应变潮汐变化的影响

海潮负荷对其深部的应力影响较大,而对其浅部的影响较小,对于离海洋较近的上海台情况正好相反； 对于武汉台,海潮负荷对θθ分量的正应变潮汐影响的相位随深度变化较大,而对其它两个分量的正应变潮汐影响的相位基本不随深度变化； 对于上海台,海潮负荷对3个正应变潮汐影响的相位基本都不随深度变化,这是因为近区的影响占主导地位,它确定了正应变潮汐的相位,这也可从图3的积分负荷格林函数看出,近区的影响要比远区的影响大得多．

3 讨论与结论

本文介绍了海潮负荷对应力和应变潮汐的影响,并计算了地球内部的负荷勒夫数及负荷格林函数． 通过积分格林函数方法计算了海潮负荷对武汉台和上海台不同深度处应力和应变潮汐的影响,并分析了其变化特征．

地球内部的负荷勒夫数与地表负荷勒夫数具有不同的变化特征,随着深度的增加,高阶负荷影响所起的作用越来越小． 除了与海洋的距离之外,深度也是影响海潮负荷大小的重要因素． 本文给出的负荷格林函数的结果表明： 在计算点附近的区域,应力负荷的影响随深度增大而减小； 而对于远离计算点的区域,应力负荷的影响则随深度增大而增大．

海潮负荷对应力和应变的影响是负荷格林函数与海潮潮高的全球积分,沿海地区的海潮负荷影响的量级甚至超过了固体潮影响的量级． 与固体潮不同,由于海潮空间分布的不规则性,导致了计算的海潮负荷的应力和应变分量在沿海和内陆台站各分量的复杂性,与计算点所处的位置有很大的关系． 对于应力或应变的某些分量来说,深度也会影响其时间变化序列的相位．

应力和应变测量是地震研究的重要手段,分析地震断层面上的潮汐应力对潮汐触发地震的研究以及地震监测预测研究也具有非常重要的意义,而海潮负荷的影响又是潮汐应力中的重要组成部分,因此海潮负荷的影响就显得非常重要,本文结果为计算海潮负荷对应力和应变的影响以及影响量级和特征提供了有益的参考．

审稿专家为本文提出了具有建设性的意见, 在此表示感谢．

方俊. 1985. 固体潮[M]. 北京： 科学出版社: 30--79.

Fang J. 1985.TheSolidEarthTides[M]. Beijing: Science Press: 30--79 (in Chinese).

蒋骏,张雁滨. 1994. 潮汐线应变组合观测的物理力学实质及其调和分析[J]. 地球物理学报, 37(增刊): 204--212.

Jiang J, Zhang Y B. 1994. The physics-mechanical significance of tidal strain combination and its harmonic analysis[J].ChineseJournalofGeophysics, 37(Suppl.): 204--212 (in Chinese).

蒋骏, 张雁滨,周翠屏. 1995. 我国应变固体潮观测及其长时间序列的调和分析结果[J]． 地壳形变与地震, 15(1)： 81--85.

Jiang J, Zhang Y B, Zhou C P. 1995. Strain tide observation and its harmonic analysis results of long series in China[J].CrustalDeformationandEarthquake, 15(1): 81--85 (in Chinese).

骆鸣津, 池顺良, 马鸿钧. 2008. 地层绝对应力测量与钻孔应变测量[M]. 北京: 地震出版社: 80--191.

Luo M J, Chi S L, Ma H J. 2008.AbsoluteStressandBoreholeStrainSurveysintheEarth’sInterior[M]. Beijing: Seismological Press: 80--191 (in Chinese).

孙文科. 1991. Love函数与地球内部的潮汐形变[J]. 地球物理学报, 34(3): 318--328.

Sun W K. 1991. Love functions and their derivatives and tidal deformation in the interior of the Earth[J].ChineseJournalofGeophysics, 34(3): 318--328 (in Chinese).

汪汉胜, 许厚泽, 李国营. 1996. SNREI地球模型负荷勒夫数数值计算的新进展[J]. 地球物理学报, 39(增刊): 182--189.

Wang H S, Hsu H T, Li G Y. 1996. Improvement of computations of load Love numbers of SNREI Earth model[J].ChineseJournalofGeophysics, 39(Suppl.): 182--189 (in Chinese).

徐建桥, 孙和平. 2003. SNREI地球对表面负荷和引潮力的形变响应[J]. 地球物理学报, 46(3): 328--334.

Xu J Q, Sun H P. 2003. Deformation response of a SNREI Earth to surface loads and tidal forces[J].ChineseJournalofGeophysics, 46(3): 328--334 (in Chinese).

许厚泽. 2010. 固体地球潮汐[M]. 武汉： 湖北科学技术出版社: 1--26.

Hsu H T. 2010.TheSolidEarthTides[M]. Wuhan: Hubei Science and Technology Press: 1--26 (in Chinese).

阎春恒, 吴小平, 许亚吉, 付虹. 2010. 印尼地区地震断层面上的海潮负荷应力特征[J]. 地震学报, 32(4): 412--422.

Yan C H, Wu X P, Xu Y J,Fu H. 2010. Characteristics of ocean tide loading stresses on earthquake faults in Indonesia region[J].ActaSeismologicaSinica, 32(4): 412--422 (in Chinese).

周江存, 孙和平. 2005. 海潮负荷对GPS基线的影响[J]. 大地测量与地球动力学, 25(4): 27--32.

Zhou J C, Sun H P. 2005. Effect of ocean tide loading on GPS baseline measurement[J].JournalofGeodesyandGeodynamics, 25(4): 27--32 (in Chinese).

周江存, 孙和平. 2007a. 武汉台长周期重力固体潮观测的海潮负荷改正[J]. 大地测量与地球动力学, 27(6): 28--32.

Zhou J C, Sun H P. 2007a. Ocean tide loading correction on long period Earth tidal gravity observation at Wuhan station[J].JournalofGeodesyandGeodynamics, 27(6): 28--32 (in Chinese).

周江存, 孙和平. 2007b. 海潮对卫星重力场恢复的影响[J]. 地球物理学报, 50(1): 115--121.

现代学徒制是学校和企业共同培养技能、老师和师傅联合教授本领和知识，以学生的技能培养为主、知识培养为辅的现代化先进人才培养的新模式和新体系。学生具有学徒和学生的双重身份，学校和企业共同承担培养学生的责任。双证融通是指将酒店管理行业准入证书所需掌握的知识和技能融入常规的酒店管理学专业教学体系，将学历证明和职业准入证明进行有机结合，学校专业教育和职业专项培训相互融通。双向融合是指学校与企业、学习专业和企业职业之间的双相贯通和融合，将职业对岗位的要求和专业对教学的要求进行渗透融通。依托国家职业技能鉴定所，在酒店设立考务分中心，组织学徒制学生在酒店岗位上完成考核认定。

Zhou J C, Sun H P. 2007b. Effect of ocean tide on recovery of satellite gravity field[J].ChineseJournalofGeophysics, 50(1): 115--121 (in Chinese).

周坤根,Ducarme B,胡延昌. 1991. 黄石台ORBES-81型伸缩仪观测结果的调和分析及海潮改正[J]. 地球物理学报, 34(4): 474--481.

Zhou K G, Ducarme B, Hu Y C. 1991. Harmonic analysis and ocean loading correction of observation by ORBES-81 extensometer at Huangshi station[J].ChineseJournalofGeophysics, 34(4): 474--481 (in Chinese).

周旭华. 2001. 高精度大地测量中倾斜及应变观测的海潮改正[J]. 地壳形变与地震, 21(3): 21--25.

Zhou X H. 2001. Ocean tidal corrections for tilt and strain observation in high precision geodesy[J].CrustalDeformationandEarthquake, 21(3): 21--25 (in Chinese).

Agnew D C. 1997. A program for computing ocean tide loading[J].JGeophysRes, 102(B3): 5109--5110.

Dziewonski A M, Anderson D L. 1981. Preliminary reference Earth model[J].PhysEarthPlanetInt, 25(4): 297--356.

Farrell W E. 1972. Deformation of the Earth by surface loads[J].RevGeophys, 10(3): 761--797.

Goad C C. 1980. Gravimetric tidal loading computed from integrated Green’s function[J].JGeophysRes, 85(B5): 2679--2683.

Longman I M. 1962. A Green’s function for determining the deformation of the Earth under surface mass loads: 1. Theory[J].JGeophysRes, 67(2): 845--850.

Longman I M. 1963. A Green’s function for determining the deformation of the Earth under surface mass loads: 2. Computations and numerical results[J].JGeophysRes, 68(2): 485--496.

Matsumoto K, Takanezawa T, Ooe M. 2000. Ocean tide models developed by assimilating TOPEX/POSEIDON altimeter data into hydrodynamical model: A global model and a regional model around Japan[J].JOceanogr, 56(5): 567--581.

Melchior P. 1978.TheTidesofthePlanetEarth[M]. California: Pergamon Press: 57--80.

Métivier L, de Viron O, Conrad C P, Renault S, Diament M, Patau G. 2009. Evidence of earthquake triggering by the solid earth tides[J].EarthPlanetSciLett, 278(3/4): 370--375.

Schuster A. 1897. On lunar and solar periodicities of earthquakes[J].ProcRSocLondon, 61: 455--465. doi:10.1098/rspl.1897.0060.

Sun W K, Sjöberg L E. 1999. Gravitational potential changes of a spherically symmetric earth model caused by a surface load[J].GeophysJInt, 137(2): 449--468.

Tanaka S, Ohtake M, Sato H. 2002. Evidence for tidal triggering of earthquakes as revealed from statistical analysis of global data[J].JGeophysRes, 107(B10): 2211. doi:10.1029/2001JB001577.

Tsuruoka H, Ohtake M, Sato H. 1995. Statistical test of the tidal triggering of earthquakes: Contribution of the ocean tide loading effect[J].GeophysJInt, 122(1): 183--194.

Zhou J C, Hwang C, Sun H P, Xu J Q, Zhang W M, Kao R, Cheng T. 2013. Precise determination of ocean tide loading gravity effect for absolute gravity stations in coastal area of China: Effects of land-sea boundary and station coordinate[J].JGeodyn, 68: 29--36.

The effect of ocean tide loading on stress and strain in the Earth’s interior

(StateKeyLaboratoryofGeodesyandEarth’sDynamics,InstituteofGeodesyandGeophysics,ChineseAcademyofSciences,Wuhan430077,China)

The theory of ocean tide loading computation was introduced in detail, and the load Love numbers and the Green’s functions in the Earth’s interior were computed for the Earth model PREM. Furthermore, the effect of ocean tide loading on strain and stress in the Earth’s interior with different depths were computed at the stations Shanghai and Wuhan as examples. The numerical results show that the depth is a critical factor in the computation. The ocean tide loading effect commonly decreases with the depth increasing for the stations in the near field, whereas the effect is contrary for the stations in the far field. Additionally, the depth affects indeed the phases of the stress and strain tensor components. And the effect over coastal area is larger than that from solid Earth tides. As a result, the ocean tide loading effect must be taken into account in the stress and strain measurements.

ocean tide loading； strain； stress； depth

10.3969/j.issn.0253-3782.2014.03.009.

国家自然科学基金(41004009, 41074053)资助.

2013-05-27收到初稿,2013-07-20决定采用修改稿.

e-mail: zjc@whigg.ac.cn

10.3969/j.issn.0253-3782.2014.03.009

P312.4

A

周江存, 孙和平, 徐建桥, 陈晓东. 2014. 海潮负荷对地球内部潮汐应力和应变的影响. 地震学报, 36(3): 433--442.

Zhou J C, Sun H P, Xu J Q, Chen X D. 2014. The effect of ocean tide loading on stress and strain in the Earth’s interior.ActaSeismologicaSinica, 36(3): 433--442. doi:10.3969/j.issn.0253-3782.2014.03.009.