Spatial domain division based on superposition pattern of seismic waves and spatial variation of motion:Case of SV wave incidence
-
摘要: 针对地震动空间差异问题,以半无限空间内平面SV波入射为例,采用波动理论和叠加原理相结合的方法研究地震动的空间特征。首先,分析入射波和反射波在空间内形成质点运动的叠加模式差异,并依此对空间域进行划分;其次,针对不同泊松比和SV波入射角情形研究叠加区与分离区的分界线控制情况;最后,对比分析不同空间域内的质点运动在峰值、持时等方面的特征。研究结果显示:当SV波入射时,叠加区与分离区的分界线通常由z3 (反射P波与反射SV波的分离线)控制;同时存在两种特殊情况,当SV波垂直入射时,分界线由z1 (入射SV波与反射SV波的分离线)控制,当反射SV波幅值为零时,分界线由z2 (入射SV波与反射P波的分离线)控制。在入射波和反射波的影响下,质点运动时程的形状具有水平不变性。三波贡献时段只出现在深度小于z1的质点的运动时程中,且持时随着深度的增加线性减少;双波贡献时段出现在位于叠加区内的质点的运动时程中,持时沿深度先增加后减少;单波贡献时段随着深度的增加而逐渐加长,在分离区达到最大值。质点运动总持时随深度逐渐增加,在分离线z2和z3深度处存在两个拐点。在质点运动峰值方面,靠近自由面的叠加区质点运动峰值变化较大,深度较大的叠加区和分离区的质点运动峰值一般不变。Abstract: For the investigation of spatial variation of seismic motion, the wave theory and superposition principle are used to illustrate the spatial variation of motions caused by seismic waves. And plane SV wave’s incidence in semi-infinite space is used as a case in this paper. Firstly the difference in superposition pattern of waves in different position is investigated and the space region is divided into different domains according to superposition pattern. Then the border line between superposition domain and separation domain with different Poisson’s ratio and SV wave’s incident angle is discussed. Finally the difference in peak value and duration time of particle motion in different domains is investigated. The result shows that, the border line is usually controlled by z3 (separation line generated by reflected P wave and reflected SV wave). Particularly, there exist two special cases, one is the vertical incidence of SV wave in which the border line is controlled by z1 (separation line generated by incident SV wave and reflected SV wave), the other is the zero amplitude of reflected SV wave in which the border line is controlled by z2 (separation line generated by incident SV wave and reflected P wave). The time history of different particle motion in one depth has the same shape with the influence of incident and reflected waves. In time history of particle motion, the period contributed by three waves only exists in depth lower than z1 and the duration time decreases with depth. The period contributed by two waves exists in all superposition domain and the duration time increases firstly and then decreases with depth. The period contributed by one single wave increases with depth and reaches its maximum in separation domain. The increase of particle motion duration has two slope change in depths of z2 and z3. The peak value of particle motion changes dramatically in superposition domain near free surface and remains constant in deep superposition domain and separation domain.
-
引言
为了尽可能地减少地震灾害造成的损失,多年来地震学家一直专注于地震前兆信息提取和地震预测研究(屈春燕等,2006;马晓静等,2009;张建国等, 2012, 2013;Zhang et al,2017)。Gorny等(1988)分析得出,1984年加兹利MS7.1地震震前曾出现显著的卫星热红外亮温异常,此后,与地震相关的卫星热红外亮温异常研究发展迅速,逐步产生了多种卫星热红外数据分析方法。例如卫星技术稳健估计方法(robust satellite technique,简写为RST)(Tronin et al,2002;Filizzola et al,2004;Ouzounov,Freund,2004;邓志辉等, 2005, 2014;单新建等,2005)、内外温差法(陈顺云等, 2006, 2013;Pulinets et al,2006;陈梅花等,2007;Saraf et al,2008)和相对功率谱方法(张元生等,2010,2011;Blackett et al,2011;Xie et al,2013;解滔等,2015a,b;张治广等,2017a)等。Filizzola等(2004)认为在地震孕育过程中,震中的地表温度高于其它区域,因此采用一种基于统计的RST方法,该方法将研究区的长期背景温度值从研究时段内研究区的亮温数据中减去,从而得到了研究时段的显著温度异常变化,并于1999年雅典MS5.9地震前提取到了显著的异常。为了进一步研究地震孕育过程中断裂带内外温度差异,陈梅花等(2007)以断裂带为界,分析了断裂带内外的温度差异及其与地震的关系,并分析得出2000年姚安MS6.5地震前,红河断裂内外温差呈稳定上升趋势,震后迅速回落至正常水平;张元生等(2010,2011)利用相对功率谱方法提取到了2011年日本MS9.0和2008年汶川MS8.0地震前的热红外亮温异常;Xie等(2013)沿用该相对功率谱方法,采用具有可变时频窗的连续小波变换计算了亮温数据的能谱,不同频段热红外相对小波能谱时空演化特征表明,该方法能够同时兼顾能谱分析的频率和时间分辨率,得到的相对小波能谱变化时频域更为平滑,并能够很好的抑制边界效应。
据中国地震台网中心测定,2018年5月28日1时50分,吉林省松原市宁江区发生MS5.7地震,震中位置为(124.71°E,45.27°N),震源深度为13 km。震中位于松原油气开发区,且处于第二松花江断裂前郭部分与扶余—肇东断裂的交会区(图 1)。地震孕育过程中,断层活动必然会导致CO2,CH4等“温室气体”释放增强,而此类气体极易吸收太阳的直接辐射和地表再发射的红外辐射,故使得地表温度逐渐升高,因此若能够检测到红外异常,一方面可以支持前人关于震前热红外亮温高值异常的研究结果(张元生等,2010;解滔等,2015a,b;张治广等,2017a),另一方面可以更有力地说明红外异常的温室效应机制。此外,由于此次地震震中位于两条断裂带的交会区,不同断裂对本次地震的贡献可通过相应地区热红外演化特征来进行揭示。
基于此,本文拟利用相对小波能谱分析方法,对松原MS5.7地震前(119°E—134°E,40°N—50°N)范围内的静止气象卫星FY-2G,FY-2E在2016年5月29日至2018年5月28日的连续观测亮温数据予以分析,以获取与此次地震相关的地震热红外异常演化特征,以期为进一步探索卫星热红外亮温异常与地震的关系积累研究震例。此外,由于本次地震的余震空间分布并未呈现出与扶余—肇东断裂与第二松花江断裂前郭部分相近的优势分布,因此此次地震的发震构造一直无法确定。本文研究分析得出,热红外相对小波能谱沿着两条断裂衰减,其中沿扶余—肇东断裂衰减较慢,因此可初步认为,扶余—肇东断裂对本次地震的贡献较为显著,可能为进一步揭示发震构造提供佐证。
1. 数据选取和处理方法
1.1 数据选取
本文选取中国静止气象卫星FY-2E和FY-2G (FY-2G为FY-2E的替代卫星)观测的热红外遥感亮温产品数据进行分析(中国气象局国家卫星气象中心,2008)。FY-2E和FY-2G两颗卫星分别于2008年6月15日和2014年12月31日发射,定点于105°E赤道上空,每1小时或30分钟对全球约1/3面积观测1次,红外波段分别为10.3—11.3 μm和11.5—12.5 μm (张璇等,2016;张治广等,2017a,b )。为了避免白天太阳照射导致的地表大幅升温造成的干扰,选取地方时夜间即北京时间00:00—04:00 (GMT 16:00—20:00)时段的5次亮温观测资料进行分析。
1.2 数据处理
小波变换作为一种线性时频分析方法,在频率域和时间域均具有很好的局部化性质,能够较好地解决频率分辨率与时间分辨率的矛盾,因此被广泛应用于地球物理勘探、地震信号处理及地震前兆分析等各项研究中(解滔等,2015a,b;王黎宏,2016)。对于一组有限的时间序列f (t),定义其小波变换为
$ {W_\varPsi }f(a{\text{,}}b) {\text{=}} \int\nolimits_{ {\text{-}} \infty }^{ {\text{+}} \infty } {f(t)\varPsi ^*_{a{\text{,}}b}} (t)\;{\rm{d}}t{\text{,}} $
(1) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$ \varPsi {\rm{(}}\omega {\rm{)}} {\text{=}}\frac {1} {\sqrt[{4}]{{\text{π}}}}{\rm{exp}}{\rm{[}}{ {\text{-}} {\frac {{{(\omega - {\omega _0})}^2}} {2}}}{\rm{]}} {\text{,}} $
(2) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$ {R_W}(a{\text{,}}b) {\text{=}}{\frac {|{W_\varPsi }f(a{\text{,}}b){|^2}} {{\overline W ^2}(a{\text{,}}b)}}{\text{,}} $
(3) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
${\overline W ^2}(a{\text{,}}b) {\text{=}} \displaystyle\frac{1}{N}\sum\limits_{l {\text{=}} 0}^{N {\text{-}} 1} {|{W_l}f(a{\text{,}}b){|^2}} {\text{,}} $
(4) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
Torrence和Compo (1998)对连续小波变换的边界效应的研究结果表明:小波能谱在接近信号边界处呈明显衰减状态,且衰减时间随着频率降低而变长;而在边界处的实际能谱值相较于利用小波变换获取的相对能谱值高很多,即如果信号在边界表现出能谱高值异常,实际的能谱异常幅值将大于计算结果,因此利用连续小波变换进行数据处理时并不会由于边界效应导致虚假异常的产生(解滔等,2015a,b)。连续小波变换为非正交变换,不仅在时域和频域有混叠,频域内相邻频段信息中也含有重叠的成分,因此本文取频率间隔对尺度因子a=0.5进行离散化处理(Torrence,Compo,1998;解滔等,2015b);然后对所选观测资料的每个像元进行计算,数据时间长度取N=730 (闰年时N=731);之后应用Xie等(2013)改进的连续小波变换对FY-2E和FY-2G卫星亮温观测数据进行处理,计算得到信号相对小波能谱时频信息后,对研究区域(119°E—134°E,40°N—50°N)进行时空扫描,分析异常区域的时空演化特征及其与本次松原MS5.7地震的关系。
当研究区域的像元上空有云层覆盖时,静止卫星FY-2E和FY-2G观测到的数据为云顶的温度,而云层下地表温度要远高于云顶的温度,为了去除云层影响,采用补窗法进行处理(张元生等,2011;解滔等,2015b),即以补窗后整点值的平均值作为像元的日均值。对所有像元进行相同处理,计算得到像元的日均值序列,采用低通滤波器提取温度日均值数据中周期T≥1年的趋势成分,计算二者的均方差,并以1.5倍均方差作为数据下限阈值,去除云层的影响(Xie et al,2013;解滔等,2015b)。本文按照上述数据处理方法,计算了松原MS 5.7地震前两年间所有像元的相对小波能谱值,并对计算结果进行了时频扫描,得到地震前后热红外亮温小波相对能谱的时空演化结果。
2. 结果分析
2.1 2018年5月28日松原MS5.7地震热红外异常分析
为了尽可能地排除离散噪声异常像元的影响,分析可能与地震等构造活动有关的热红外亮温高值异常现象,本文按照以下3条准则进行异常信息判别(解滔等,2015a,b):① 异常像元须聚集成一定范围的异常区域,而非零星离散的分布;② 显著异常持续时间须在15天以上;③ 异常区域须沿构造断裂带分布,特别是活动断裂带。
对所选数据的小波相对能谱的时空变化特征(图 2)进行分析可以看出:2018年5月28日松原MS5.7地震前,震中附近的断裂带出现了明显的热红外亮温异常现象。从空间上看,2018年2月初在震中附近NW向的第二松花江断裂前郭部分出现异常;2018年3月,异常区开始沿NE向的扶余—肇东断裂展布,3月中旬异常空间分布范围达到峰值,随后异常区域逐渐缩小;至2018年4月仅NE向的扶余—肇东断裂区域存在条带状分布的异常区域,直到2018年4月中旬异常消失。异常空间分布呈现出NW向与NE向的条带交会状,即“V”型分布特征,未来震中位于两条带状热红外异常区域在迁移过程中形成的交会带附近。从相对小波能谱异常幅度上看,2018年2月初在震中附近NW向的第二松花江断裂前郭部分出现异常,随后异常幅度逐渐增大,3月中旬异常幅度达到最大,能谱相对变化达到均值的7倍以上,之后异常幅值逐渐衰减,直至2018年4月中旬异常消失,异常持续时间达3个月。本次地震热异常呈现出与以往震例类似的演化过程,即由初始小范围异常发展至大区域高能量异常,随后开始急速减弱的过程(张璇等,2016)。因此,2018年2月至4月于松原地区出现的热红外亮温异常可能与此次松原MS5.7地震有关。
2.2 研究区历史震例分析
尽管地震前震中及附近区域表现出的地表增温异常现象已被大量的观测结果所证实,但由于地震孕育机理的复杂性和地震热红外异常的多样性,地震热红外成因机理的研究仍然是一个难点(解滔等,2015a)。郭增建和秦保燕(1979)认为,大震震源区锁住的断层会发生预滑移,沿断层带的温度会升高,因而产生震前热红外亮温异常现象。Cicerone等(2009)讨论了地震前卫星热红外亮温异常现象,结果显示,在富含地热和油气的区域地震前出现热红外亮温异常现象的可能性更大,因此解滔等(2015b)分析得出,地震孕育过程中应力逐渐增大,当应力接近岩石临界破裂状态时,会出现大量贯通地表的裂隙,地下逸出气体也随之增多,导致地下岩层与地表的热对流加强。大量裂隙的出现及热对流的增强导致了地震前热红外异常的产生,且异常区沿震中附近的断裂带集中分布(Tronin,1996,2000;马瑾等,2006;强祖基等,2009,2010;解滔等,2015b)。
需要特别说明,热红外异常只是表明在异常时段该区域地表温度较正常背景出现了异常变化,而温度变化并非仅由地震所致,因此,地震热异常分析过程中不可避免地会出现虚假地震异常现象。本文将2013年以来研究区域内热红外亮温异常与地震对应关系列于表 1。在2013年1月至2018年5月期间,依据本文的异常判别准则,该研究区内共出现了7组热红外亮温相对能谱异常,其中3组异常出现后无地震对应,另外4组分别对应于2013年4月22日内蒙古与辽宁交界处的MS5.3地震、2013年6月20日内蒙古与黑龙江交界处MS5.0地震、2013年10月31至11月23日吉林松原市前郭MS5.8震群和本次2018年5月28日松原MS5.7地震,这4组地震前均在震中附近出现了亮温热红外相对能谱高值异常,且地震均发生在异常区边缘。在2013年1月至2018年5月期间,研究区域内共发生MS≥5.0地震5组,除上述4组地震外,还有2016年1月2日发生在黑龙江省牡丹江市林口县的MS6.4地震,该次地震的震源深度为580 km,为深源地震,但风云卫星没有记录到该地震产生的亮温异常。为便于对比分析,选取异常幅值最大的日期作为异常极值日,有震异常和无震异常的极值空间分布如图 3所示。可以看出,有震异常空间分布较为聚集,异常幅值较大,而无震异常空间分布较有震异常更零散,异常幅值相对较小。
表 1 2013—2018年研究区热红外亮温异常与地震对应情况Table 1. The thermal infrared anomaly and corresponding earthquakes in the studied area异常出现时间 异常空间范围 对应地震 2013年1—3月 121°E—123°E,42°N—44°N 2013年4月22日内蒙古与辽宁交界MS5.3地震 2013年2—3月 124°E—126°E,46°N—50°N 2013年6月20日内蒙古与黑龙江交界MS5.0地震 2013年10—12月 121°E—124°E,42°N—46°N 2013年10月31日—11月23日吉林省松原市前郭MS5.8震群 2014年4—7月 129°E—131°E,44°N—46°N 无 2015年3—5月 121°E—122°E,47°N—49°N 无 2016年6—7月 121°E—122°E,40°N—41°N 无 2018年2—4月 125°E—127°E,44°N—46°N 2018年5月28日吉林松原MS5.7地震 图 3 历史地震及热红外相对能谱极值的空间分布(a) 有震异常相对能谱极值空间分布;(b) 无震异常相对能谱极值空间分布Figure 3. Spatial distribution of thermal infrared extreme value anomaly and history earthquakes(a) Spatial distribution of thermal infrared extreme value anomaly and corresponding earthquakes;(b) Spatial distribution of thermal infrared extreme value anomaly without corresponding earthquakes距离松原MS5.7地震最近的为2013年10月31日至11月23日松原MS5.8震群,两次地震震中相距仅为88 km。本文以该震群为例,对比分析了两组地震前后热红外亮温异常特征及其差异,如图 4所示。从异常及发震时间上看,松原MS5.7地震发生于热红外亮温异常消失后约1个月,而2013年10月31日至11月23日松原MS5.8震群发生在异常出现至异常幅值逐渐增大的过程中,地震后异常幅值逐渐减小。从空间上看,两组地震均发生在异常区边缘,且邻近异常衰减最慢区域。不同的是,2018年5月28日松原MS5.7地震前出现的热红外亮温异常呈现出显著的沿发震构造分布的特点,即异常主要沿第二松花江断裂前郭地区和扶余—肇东断裂展布,而2013年松原震群前异常首先出现在震中距较远的嫩江断裂,其后才出现在北东向的扶余—肇东断裂处。这种大范围远距离的热红外亮温异常可能与地下逸出气体增多以及地下岩层与地表的热对流增强有关(徐秀登等,1991;强祖基等,1992;Tronin,1996,2000)。
3. 讨论与结论
本文针对2018年5月28日吉林松原MS5.7地震,利用连续小波相对能谱方法对发震前两年期间中国静止气象卫星FY-2E (FY-2G)连续观测的亮温资料进行分析,分析结果显示,在震中附近的第二松花江断裂前郭段及扶余—肇东断裂区域均在震前出现显著的相对小波能谱高值异常。震前两个月左右异常幅值和规模开始衰减并于4月中旬消失。地震震中位于异常区边缘的活动断层上,发震时间为异常峰值后的两个月内,与前人的研究结果基本一致。因而认为本文讨论的2018年2月至4月松原地区热红外异常可能与松原MS5.7地震有关。通过对2013年以来研究区历史地震及其与热红外亮温相对小波能谱高值异常的统计可初步发现,有震异常的相对小波能谱高值异常在空间上通常较为集中,呈聚集分布,且异常幅值较高;而无震异常的相对小波能谱高值异常空间上分布较为零星,异常幅值也较低。但由于统计样本个数的限制,目前还无法定量给出两者之间的差别,更多异常特征和异同需要基于更大样本的统计。因此,下一步拟基于中国大陆地区历史地震开展有震异常和无震异常的定量分析。
-
表 1 不同深度质点地震波贡献时段
Table 1 Contribution time interval of waves on particle motion at different depth
质点 深度/km 单波贡献时段/s 双波贡献时段/s 三波贡献时段/s Si Pr Sr (Si,Pr) (Pr,Sr) (Si,Pr,Sr) D1 0 0—2.00 D2 1.0 0—0.89 2.89—3.17 0.89—1.17 2.00—2.89 1.17—2.00 D3 1.8 0—1.61 2.00—2.11 3.61—4.11 1.61—2.00 2.11—3.61 D4 3.0 0—2.00 2.68—3.52 4.68—5.52 3.52—4.68 D5 9.0 0—2.00 8.05—10.05 10.57—12.57 -
陈厚群,徐泽平,李敏. 2008. 汶川大地震和大坝抗震安全[J]. 水利学报,39(10):1158–1167. doi: 10.3321/j.issn:0559-9350.2008.10.002 Chen H Q,Xu Z P,Li M. 2008. Wenchuan earthquake and seismic safety of large dams[J]. Journal of Hydraulic Engineering,39(10):1158–1167 (in Chinese).
贺春晖,王进廷,张楚汉. 2017. 基于震源—河谷波场数值模拟的坝址地震动参数确定方法[J]. 地球物理学报,60(2):585–592. He C H,Wang J T,Zhang C H. 2017. Determination of seismic parameters for dam sites by numerical simulation of the rupture-canyon wave field[J]. Chinese Journal of Geophysics,60(2):585–592 (in Chinese).
何卫平,周宜红,赵春菊. 2019. 平面P波入射半无限介质空间域划分及运动偏振特征[J]. 岩石力学与工程学报,38(2):321–331. He W P,Zhou Y H,Zhao C J. 2019. Spatial division of semi-infinite space and polarization characteristics of particle motion for incident plane primary wave[J]. Chinese Journal of Rock Mechanics and Engineering,38(2):321–331 (in Chinese).
金星. 1993. 地震动场研究状况简介[J]. 世界地震工程,(3):19–22,52. Jin X. 1993. Introduction of research on earthquake motion field[J]. World Earthquake Engineering,(3):19–22,52 (in Chinese).
金星,廖振鹏. 1994. 地震动随机场的物理模拟[J]. 地震工程与工程振动,14(3):11–19. Jin X,Liao Z P. 1994. Physical simulation on the random field of seismic motion[J]. Earthquake Engineering and Engineering Vibration,14(3):11–19 (in Chinese).
林皋. 2017. 地下结构地震响应的计算模型[J]. 力学学报,49(3):528–542. Lin G. 2017. A computational model for seismic response analysis of underground structures[J]. Chinese Journal of Theoretical and Applied Mechanics,49(3):528–542 (in Chinese).
缪惠全,李杰. 2018. 基于物理机制的工程场地地震动相干函数模型[J]. 中国科学:技术科学,48(2):209–216. Miao H Q,Li J. 2018. The coherence function model of ground motion in engineering sites based on physical mechanism[J]. Scientia Sinica:Technologica,48(2):209–216 (in Chinese).
张楚汉,金峰,王进廷,徐艳杰,潘坚文. 2016. 高混凝土坝抗震安全评价的关键问题与研究进展[J]. 水利学报,47(3):253–264. Zhang C H,Jin F,Wang J T,Xu Y J,Pan J W. 2016. Key issues and developments on seismic safety evaluation of high concrete dams[J]. Journal of Hydraulic Engineering,47(3):253–264 (in Chinese).
张翠然,陈厚群,李敏. 2011. 采用随机有限断层法生成最大可信地震[J]. 水利学报,42(6):721–728. Zhang C R,Chen H Q,Li M. 2011. Generation of the maximum credible earthquake by using the stochastic finite fault method[J]. Journal of Hydraulic Engineering,42(6):721–728 (in Chinese).
钟万勰,林家浩,吴志刚,孙东科,张亚辉. 2000. 大跨度桥梁分析方法的一些进展[J]. 大连理工大学学报,40(2):127–135. doi: 10.3321/j.issn:1000-8608.2000.02.001 Zhong W X,Lin J H,Wu Z G,Sun D K,Zhang Y H. 2000. Development of analytical method for long-span bridges[J]. Journal of Dalian University of Technology,40(2):127–135 (in Chinese).
Hao H,Oliveira C S,Penzien J. 1989. Multiple-station ground motion processing and simulation based on SMART-1 array data[J]. Nucl Eng Des,111(3):293–310. doi: 10.1016/0029-5493(89)90241-0
Loh C H,Yeh Y T. 1988. Spatial variation and stochastic modelling of seismic differential ground movement[J]. Earthq Eng Struct Dyn,16(4):583–596. doi: 10.1002/eqe.4290160409
Luco J E,Wong H L. 1986. Response of a rigid foundation to a spatially random ground motion[J]. Earthq Eng Struct Dyn,14(6):891–908. doi: 10.1002/eqe.4290140606
Shearer P M. 2009. Introduction to Seismology[M]. 2nd ed. Cambridge: Cambridge University Press: 215−240.
Zerva A. 2009. Spatial Variation of Seismic Ground Motions: Modeling and Engineering Applications[M]. Boca Raton: CRC Press: 1–64.
Zerva A,Harada T. 1997. Effect of surface layer stochasticity on seismic ground motion coherence and strain estimates[J]. Soil Dyn Earthq Eng,16(7/8):445–457.
-
期刊类型引用(3)
1. 李继业,胡澜缤,李营,马龙辰,王强,张思萌,李冬妮. 松原M_S5.1地震前断层土壤气H_2、Hg地球化学特征与热红外异常响应研究. 地震工程学报. 2023(04): 933-945 . 百度学术
2. 李利波,张翔,李智蓉. 2021年云南省3次中强地震前后热红外异常特征分析. 地震地磁观测与研究. 2022(05): 16-22 . 百度学术
3. 米日阿依·买土地,玉素甫江·如素力,喀迪阿依·阿力木,艾尔肯·图尔孙. 基于STL和GESD算法的于田两次Ms7.3地震热异常时空分布特征识别研究. 河南科学. 2021(03): 421-428 . 百度学术
其他类型引用(2)