The Background Variation of Natural Source ELF and Its EM Abnormal Phenomena in Yunnan Earthquakes

In order to add earthquake monitoring methods and develop new method research, the ELF Network for Earthquake Monitoring selected 30 stations in the Capital Circle and the Sichuan-Yunnan region. Finding electromagnetic field background variation is the basis of distinguishing the seismic electromagnetic anomalies. This paper introduces the data acquisition and selection of the Shexian, Anqiu, Lijiang and Dali stations which have recorded for longer time with better data and are located on the similar latitude. Then we use the natural source electromagnetic field''s auto-power spectrum to express the intensity of the electromagnetic field. By using power spectral data of many frequencies in the observation frequency band, after the data pre-processing and sliding average noising, the background variation of extremely low frequency stations and the range ability were acquired. Taking the Baoshan M5.1 earthquake on October 30, 2015 and Dali M5.0 earthquake on May 18, 2016 as examples, the authors analyzed the earthquake electromagnetic anomaly characteristic of ELF stations around the earthquakes.     

ELECTROMAGNETIC ANOMALOUS PULSES BEFORE AND AFTER EARTHQUAKES(M>6.0) IN CHINA CONTINENT IN 2008

In this paper,we provided the electromagnetic anomaly phenomena prior to five earthquakes with magnitude M_W>6.0 occurring in China continent in 2008.The electromagnetic data in frequency bands 800~0.1Hz are recorded at two stations located in Urumqi,Xinjiang Uygur Autonomous Region (URQ) and Tonghai,Yunnan Province (TH).The time series of four horizontal components of the electromagnetic field at the stations are converted to the spectra using Fast Fourier Transform.The relationship between the electromagnetic anomalous phenomena and the earthquake events is studied through analyzing the temporal variation of electromagnetic spectra and comparing them with earthquake events.This study provides the new examples for electromagnetic anomaly phenomena before the earthquakes.The following features can be found.①The obvious anomalous power spectrum density (PSD) of electromagnetic fields in frequency band of 128~0.5Hz appeared before earthquakes.The anomalous PSD is related to the earthquake magnitude and the epicenter distance.The anomaly size of magnetic PSD is about 1-3 orders of magnitude bigger than the background field before Wenchuan M_S8.0 earthquake.② The changes of parameters are also related with frequency of the signal;The ELF signal strength observed at seismic stations is significantly related with distance between the stations or the measured magnetic and electrical field components and the source of earthquake.We also found that the field strength attenuated faster in area closer to the source than in far area.It is clear that the magnetic signals have higher signal noise ratio than the electric data.③ The amplitudes are also related with the azimuth of the earthquake to station and with the frequency of the data.④ The anomalous pulses showed cluster phenomenon in time.     

地磁日变与地震活动关系的主成分分析——以日本岩手县北部6.1级地震为例

本文提出了将主成分分析方法和地磁日变分析方法相结合以从较强干扰背景中提取相对较弱地震地磁信息的一种新思路.具体以1998年日本岩手县北部6.1级地震为例,选用三个地磁观测台资料进行研究.利用调和分析方法得到了各台基于多次谐波拟合的地磁日变曲线,在此基础上对各台的地磁日变形态进行了研究,发现距震中最近地磁台的日变形态在地震前大约两周出现了明显的异常,而其他两个距离震中相对较远的观测台日变形态基本正常.最后,将主成分分析方法应用于上述拟合地磁日变结果,得到了各主成分及其所占能量比的时间变化,结果表明在地震前两周左右第二主成分所占能量比显著增加,而且距离震中最近的台站的上述变化明显高于另两个较远的台站.上述结果表明地震前两周左右检测到的地磁日变异常可能与震中附近地下电阻率的变化或孕震过程中产生的电磁信号存在一定关系,相关研究结果有助于加深对地震电磁现象的认识和理解.     

基于主成分方法分析2013年甘肃岷县漳县6.6级地震前地磁日变化异常

选用2013—2014年兰州台、 天水台、 寺滩台3个台站地磁垂直分量分钟值数据, 利用调和分析剔除1~4次地磁日变化谐波成分, 在此基础上对3个台站地磁Z分量数据进行主成分分析, 计算结果显示: 2013年7月22日岷县漳县6.6级地震前, 最小的第三主成分在2013年4月20日至6月10日出现明显的增强变化, 超出阈值, 异常持续50 d左右; 同时, 在信号的子空间中特征矢量的方向也出现同步的异常扰动, 震后恢复到原来的水平, 其他时段第三主成分未出现类似变化。 据此认为, 这些异常信号与2013年岷县漳县6.6级地震相关。 研究结果表明, PCA方法对地震电磁异常的提取是有效的。     

Establishment of the ELF Network in Yunnan and Electromagnetic Precursory Monitoring Results of the Yangbi MS5.1 Earthquake on March 27, 2017

This paper summarizes the layout of the Yunnan seismic ELF electromagnetic observation network, site selection, ELF electromagnetic instrument system, data processing and other construction. The principle and method of using the ELF electromagnetic wave technique to monitor and predict earthquakes are expounded. The long term monitoring of ELF electromagnetic fields is carried out in the Yunnan earthquake prone area, and at the same time, the changes in electrical parameters and spatial electromagnetic fields of the regional crustal medium structure are monitored. The functions such as automatic, quasi real time, remote monitoring, network monitoring, data processing specialization, data service, data sharing and industrialization of the ELF electromagnetic observation data have been realized. In order to capture the deep electromagnetic precursory information of the earthquakes, service for earthquake prediction research, which has broad application prospects and development potential. Through the research of the seismicity of Yunnan in the trial run period of the project, the preliminary results of the extreme low frequency electromagnetic observation of the Yunnan earthquake in recent years are given. The electromagnetic precursors and the electromagnetic effects of the Yangbi earthquake are recorded. In the 3-month period before the earthquake, the power spectrum of the electric and magnetic fields, the apparent resistivity and the impedance phase in the observed signals are all abnormal, and gradually increased with time. The maximum value is reached 20 days before the earthquake, and an earthquake occurs when the change is restored to normal.     

东北地区中强地震与地磁总场强度异常变化关系研究

应用东北三省及内蒙古东部地区的12个地磁台站总场强度整点值资料,利用空间相关法研究了总场强度数据与黑龙江及邻近地区中强地震的关系,给出了该方法在黑龙江地区的异常参数,统计了黑龙江省及邻区地区2008年以来发生的级地震前的地磁空间相关异常的次数,对震前异常进行了系统的统计和分析,发现震前空间相关法都有较好的映震效能。     

几次大震前的地面和空间电磁场变化

本文应用非线性谱估计方法(MEM方法)、电/磁场脉冲能量统计方法,处理了四川汶川MS8.0 (2008)、新疆于田MS7.3 (2008)和青海玉树MS7.1 (2010) 地震及强余震附近的地面电/磁场和Demeter卫星电离层磁场观测数据,研究了地震电/磁场变化,得到:在震中附近及周围出现了震前地电/地磁场极低频成...     

小波极大值方法及其在电磁异常信号提取中的应用

小波极大值法作为一种数据挖掘技术已经受到关注,并在地震科学研究领域开始进行应用研究。它能在长时间的观测数据中快速提取出不同时间、不同频率的有关异常信息。文中介绍了如何应用小波极大值法提取电磁场异常信息,即利用连续小波变换计算小波系数,计算各尺度小波系数的极大值,结合表征信号奇异性的Lipschitz指数分析异常的真实性,进而研究极大值在时间分布上的特点或规律,捕捉有关的电磁异常现象。利用该方法对2013年4月20日四川芦山7.0级地震后35d的连续观测的电磁场数据进行了分析,发现存在3次电磁场异常,探讨了这些异常可能与地震活动性之间存在的关系。尽管用这种方法尚不能确定发现的电磁场异常现象与地震有某种特定的关系,但是证明了小波极大值方法在连续观测数据中提取电磁异常现象的有效性。     

华北地区地磁加卸载响应比方法应用研究

运用地磁加卸载响应比方法(LURR)对华北地区23个地磁台站2008—2014年的数字化地磁相对记录数据进行计算分析,结果显示华北地区7年内发生了38次共计10组加卸载响应比异常,其中7组有震例对应,在选定的24次震例中16次有震前异常。经统计发现各次异常发生时其异常面积与最大异常强度成正比关系,各组异常内累计异常台站数量与组内地震辐射能量有一定相关性。对有震例的7组异常分析后认为多数情况下组内各次异常之间有一定联系,主要表现为位置的连续性及继承性,并且发现震中位置和异常位置分布有一定关系。     

门源MS6.9地震前地电场优势方位角与加卸载响应比异常变化分析

2022年1月8日在青海门源县发生M_S6.9地震,基于震中300 km范围内地电场近5年观测资料,综合分析选取9个观测站,根据大地电场岩体裂隙水(电荷)渗流(移动)模型计算其优势方位角,并尝试以地电场为响应量,通过库仑应力触发模型的加卸载响应比(LURR)计算方法,计算地电场LURR值。结果显示：(1)两种不同方法计算的地电场异常站在空间分布上具有一致性,其中古丰、黄羊川、寺滩和兰州站地电场优势方位角、LURR在震前皆出现异常变化,而山丹等其他站方位角、LURR均看不出明显的异常变化。(2)兰州和寺滩站两种计算方法的结果时序变化较为吻合,表现出准同步性。进一步结合震源机制解对异常观测站优势方位与区域主压应力P走向的关联进行分析,结果表明基本符合岩石物理学理论,这在一定程度上可增强地电场优势方位角方法在分析地震前兆异常中的可信度。地电场优势方位角以及LURR值两种计算方法在机理上具有关联性,综合分析其异常演化特征可能有助于进一步认知地震孕育的物理过程。     

2022年1月8日青海门源MS6.9地震前地磁垂直强度极化异常变化特征

基于甘青川地区的14个地磁台站秒采样资料,对2022年1月8日青海门源M_S6.9地震前地磁垂直强度极化异常的时空变化特征进行分析研究。计算结果显示,在2021年10月底出现了地磁极化超阈值高值异常变化,异常台站主要分布在青海、甘肃和四川地区,异常出现后73天发生了门源地震,震中位于极化高值异常阈值线附近。研究还发现,此次地磁极化异常具有一定的时空变化特征,时间上,不同台站极化异常具有较好的时间同步性;空间上,极化异常高值区呈现出沿震中附近出现后不断扩展最终再向震中收缩的特点。此外,各地磁极化异常台站的归一置零极化值、异常持续时间与震中距存在较好的负相关性,异常台站距离地震震中越近,其归一置零极化值越高,异常持续时间也越长,这一特征符合地震电磁扰动信号的衰减特征。根据极化异常和地震的时空关系分析认为,此次地磁极化高值异常对应了之后在异常高值区边缘发生的门源M_S6.9地震。     

2022年青海门源MS6.9地震前地电场频谱特征研究

应用最大熵谱法研究2022年1月8日青海门源M_S6.9地震震中周围松山、山丹、高台、瓜州地电场台站的观测数据，分析地震发生前后地电场功率谱密度变化。结果表明：(1)门源地震发生前，松山等4个台站的地电场低频成分功率谱密度值出现增大的现象；(2)4个台站在门源地震前谱值变化的时间基本同步；(3)结合微裂隙机制及“多点场”观点，对比分析认为震源孕育激发的电磁辐射是造成震前电磁异常现象的主要原因。     

SHIFTING SELF-CORRELATION METHOD INITIALY USED IN THE APPARENT RESISTIVITY OBSERVATION DATA

In this paper,we propose a method of seismic prediction using the geo-electric resistivity shifting self-correlation (SSC),and a numerical test is carried out using random time series analysis to verify the validity of the method.The SSC method is applied to the actual observation data of three geo-electric resistivity stations,and results are obtained as follows:(1) SSC coefficient changes in Ganzi and Shandan stations have good correspondence to earthquake,which is represented mainly by the phased increase of correlation coefficient appearing six months to a year before the earthquake.At the same time,the correlation coefficient anomalies of the two stations also exhibit strong anisotropy.(2) Although Chengdu geo-electric resistivity station had suffered serious disturbance,the correlation coefficient anomaly also has a good correspondence with earthquake.In addition to the validity of the SSC method,it may also be attributed to the magnitude of the earthquake event,the smaller distance of epicenter,and the time of the earthquake.Anisotropy also exists in the anomaly at Chengdu station.(3) By comparing the characteristics of different magnitudes of earthquakes,the results are obtained that,when the magnitude of the selected characteristic earthquake is relatively small,the amplitude of the anomaly before earthquake is different,but when the magnitude is larger,for example M_S ≥ 5.0,the impact on the results of this study is very limited.In addition,we briefly discussed the anisotropy of seismic geoelectrical resistivity anomalies and the selection of the characteristic earthquake.     

Shifting Correlation Between Earthquakes and Electromagnetic Signals: A Case Study of the 2013 Minxian–Zhangxian ML 6.5 (MW 6.1) Earthquake in Gansu,China

The shifting correlation method (SCM) is proposed for statistical analysis of the correlation between earthquake sequences and electromagnetic signal sequences. In this method, the two different sequences were treated in units of 1 day. With the earthquake sequences fixed, the electromagnetic sequences were continuously shifted on the time axis, and the linear correlation coefficients between the two were calculated. In this way, the frequency and temporal distribution characteristics of potential seismic electromagnetic signals in the pre, co, and post-seismic stages were analyzed. In the work discussed in this paper, we first verified the effectiveness of the SCM and found it could accurately identify indistinct related signals by use of sufficient samples of synthetic data. Then, as a case study, the method was used for analysis of electromagnetic monitoring data from the Minxian–Zhangxian M_L 6.5 (M_W 6.1) earthquake. The results showed: (1) there seems to be a strong correlation between earthquakes and electromagnetic signals at different frequency in the pre, co, and post-seismic stages, with correlation coefficients in the range 0.4–0.7. The correlation was positive and negative before and after the earthquakes, respectively. (2) The electromagnetic signals related to the earthquakes might appear 23 days before and last for 10 days after the shocks. (3) To some extent, the occurrence time and frequency band of seismic electromagnetic signals are different at different stations. We inferred that the differences were related to resistivity, active tectonics, and seismogenic structure.     

2017年西藏米林MS6.9地震前地磁垂直强度极化异常特征

基于中国大陆65 个地磁台站秒采样观测资料,利用地磁垂直强度极化方法提取了2017年西藏米林M_S6.9 地震震中及附近区域超低频电磁异常信号并分析其时空演化特征。结果显示:米林M_S6.9 地震前在青藏高原出现了大范围的地磁垂直强度极化高值异常,异常过程共持续9天;高值异常开始于2017 年10 月30 日,高值持续3 天之后出现短时间下降,随即转折上升,再次出现高值异常并持续了4 天,在此过程中各高值台站的时序曲线呈现出单峰或双峰的形态。空间分布图显示高值异常在中国大陆西部区域反复出现,尤其是在青藏高原巴颜喀拉地块与羌塘地块的交界处。异常呈“出现—扩大—收缩—消失—扩大—消失”的演化过程,10月31日异常面积达到最大值。异常过程结束后10天发生了西藏米林M_S6.9地震,震中距离10 月31 日异常阈值线5 km。综合分析认为,此次异常与米林地震具有较强的时空相关性,是可靠的地震电磁前兆异常。      

Modeling spatially correlated spectral accelerations at multiple periods using principal component analysis and geostatistics

Regional seismic risk assessments and quantification of portfolio losses often require simulation of spatially distributed ground motions at multiple intensity measures. For a given earthquake, distributed ground motions are characterized by spatial correlation and correlation between different intensity measures, known as cross‐correlation. This study proposes a new spatial cross‐correlation model for within‐event spectral acceleration residuals that uses a combination of principal component analysis (PCA) and geostatistics. Records from 45 earthquakes are used to investigate earthquake‐to‐earthquake trends in application of PCA to spectral acceleration residuals. Based on the findings, PCA is used to determine coefficients that linearly transform cross‐correlated residuals to independent principal components. Nested semivariogram models are then fit to empirical semivariograms to quantify the spatial correlation of principal components. The resultant PCA spatial cross‐correlation model is shown to be accurate and computationally efficient. A step‐by‐step procedure and an example are presented to illustrate the use of the predictive model for rapid simulation of spatially cross‐correlated spectral accelerations at multiple periods.     

2017年九寨沟7级地震前电磁扰动异常分析

地震前的电磁扰动现象曾被多次观测到, 但需进一步研究。 在中国地震局的支持下, 可以监测电磁扰动数据的多分量地震监测系统AETA在中国境内已经布设了200余套。 2017年8月8日, 四川省九寨沟县发生7级地震。 震中周边500 km范围内当时已经布设22套AETA设备。 基于滑动四分位距法(时间窗为27 d)分析了这些台站震前30 d的数据。 11个台站出现了异常, 其中震中周边200 km内的6个台站有4个发生异常。 同时, 正异常明显多于负异常。 另外, 大部分异常发生在夜间(20时~6时)。     

2017年吉林宁江M 4.9地震前电磁变化

2017年7月23日吉林宁江M 4.9地震前,地震影响范围内电磁台站出现多种电磁异常,其中4项电磁手段异常明显,分别是：地磁垂直分量日变幅加卸载响应比异常、地磁场总场强度F₂₁每日一值空间相关异常、绥化地震台地电阻率EW向反向下降异常及林甸电阻率EW向整点值变化异常簇异常。本文对地磁异常进行统计、分析,为东北地区未来利用电磁方法分析预测地震提供经验。     

盈江5.8级和缅甸7.2级地震前电磁异常

本文分析了2010年10月-2011年3月云南通海ELF地震台站观测的电磁场资料,研究发现该资料变化与2011年发生在中国大陆"3·10"盈江5.8级及"3·24"缅甸7.2级地震存在以下关系,(1)470～1 Hz的磁场自功率密度谱(PSD)在2011年盈江5.8和缅甸7.2级地震前出现了脉冲丛集异常,异常幅度分别比正常月份变化约1～3个数量级.异常的幅度与观测点的电极方位和信号的频率有密切的关系.(2)2011年3月10日盈江地震的异常变化是先从地壳的深部向浅部发展,在震前25天2.6 km处的电阻率开始发生变化,震前10天400 m处的电阻率变化,出现了勺形的变化形态,阻抗相位同步变化.3月24日的缅甸7.2级地震,震前10天深部、浅部的电阻率同步发生了大幅度变化,阻抗相位发生转折,这可能是强震的变化特性.     

航空电磁数据主成分滤波重构的噪声去除方法

主成分分析方法利用低阶主成分重构航空电磁数据,解决了航空电磁探测中噪声与数据在频谱重叠情况下的噪声压制问题,但是参与重构的低阶主成分仍包含高频空间噪声,影响数据成像精度.本文提出的主成分滤波重构去噪方法,根据自适应窗宽平滑算法,设计了主成分低通滤波器组,对参与重构的低阶主成分进行测线滤波,再将滤波后的低阶主成分重构为电磁信号,不仅可以去除低阶主成分中的高频空间噪声,而且去除了高阶主成分包含的不相关噪声.仿真数据的去噪结果表明,主成分滤波重构获得较高的信噪比,较常规测线滤波与主成分重构分别提高了10.96dB和2.52dB;电导率深度成像结果证明了主成分滤波重构方法能够提高地下深部异常体的识别能力.最后通过实测数据的成像结果进一步验证了本文研究的主成分滤波重构去噪方法的有效性.