新疆地區(qū)電磁預(yù)測指標(biāo)體系在2023年1月30日沙雅M_S6.1地震中的應(yīng)用

艾薩?伊斯馬伊力 翟世龍 張桉赫 馮麗麗 賀曼秋 馮志生

摘要：2023年1月30日新疆沙雅縣（82.29°E，40.01°N）發(fā)生M_S6.1地震，震源深度50 km。地震發(fā)生前，中國局地震短臨預(yù)報(bào)跟蹤電磁學(xué)科組和新疆地震局電磁學(xué)科組依據(jù)中國地震局《地震電磁分析預(yù)測技術(shù)方法工作手冊》、《新疆地區(qū)地震分析預(yù)測技術(shù)方法工作手冊》和《地震電磁分析預(yù)報(bào)方法清單》對新疆地區(qū)電磁臺站的數(shù)據(jù)進(jìn)行異常梳理，發(fā)現(xiàn)新疆地區(qū)電磁臺站存在流動地磁、地磁垂直強(qiáng)度極化、地磁逐日比、地磁低點(diǎn)位移等多項(xiàng)異常。震前異常綜合預(yù)測意見與實(shí)發(fā)地震基本一致，震前中短期漸進(jìn)式異常演化特征較突出，2017—2020年新疆地區(qū)地震電磁學(xué)科預(yù)測指標(biāo)體系建設(shè)工作已初見成效。

關(guān)鍵詞：沙雅M_S6.1地震；電磁異常；異常指標(biāo)；地震預(yù)報(bào)

doi：10.16256/j.issn.1001-8956.2023.03.008

2017―2020年，中國地震局地震短臨跟蹤電磁學(xué)科組聯(lián)合新疆電磁學(xué)科組對新疆及周邊地區(qū)電磁資料及主要提取和識別異常的分析方法進(jìn)行了回溯性震例分析研究。此項(xiàng)研究梳理出多項(xiàng)地震預(yù)報(bào)效能和眏震效能較好的方法和手段：地震地磁場日變化異常分析方法（地磁垂直分量（Z）日變幅加卸載響應(yīng)比法（以下簡稱地磁加卸載響應(yīng)比）、地磁垂直分量（Z）日變化幅度逐日比法（以下簡稱地磁逐日比）、地磁低點(diǎn)位移法）、地震地磁擾動異常分析方法（地磁垂直強(qiáng)度極化法）和地磁測深視電阻率（地磁諧波振幅比）異步法、地磁垂直分量（Z）歸零法（以下簡稱地磁歸零法）、直流視電阻率形態(tài)法、地電場優(yōu)勢方位角法等^［1-15］。國內(nèi)外許多學(xué)者對電磁分析方法和手段的異常機(jī)理進(jìn)行了分析研究^{［1，3，5，12］}。理論上這些異常現(xiàn)象在地震孕育過程中地下介質(zhì)電磁參數(shù)突變或出現(xiàn)短時感應(yīng)電流，會導(dǎo)致地磁低點(diǎn)位移、地磁逐日比及地磁加卸載響應(yīng)比的突變分界線或閾值線，在空間上呈現(xiàn)高曲率形態(tài)并在較大范圍呈現(xiàn)２個或2個以上獨(dú)立的異常區(qū)。地磁日變化類異常原始數(shù)據(jù)分析發(fā)現(xiàn)，原始數(shù)據(jù)形態(tài)上均存在不同程度的地磁垂直分量“反相位”現(xiàn)象。因此地磁日變化低點(diǎn)位移法的異常突變分界線、地磁逐日比和地磁加卸載響應(yīng)比閾值線高曲率段往往與孕震區(qū)高導(dǎo)體的走向基本一致^［1-3］。在此基礎(chǔ)上，建立新疆及周邊區(qū)域的地磁日變化異常方法和綜合地震預(yù)測指標(biāo)體系，對新疆地區(qū)不同構(gòu)造區(qū)域各地磁方法識別的日變化異常判定指標(biāo)和地震預(yù)測規(guī)則進(jìn)行詳細(xì)論證和總結(jié)，同時編寫《地震電磁分析預(yù)測技術(shù)方法工作手冊》、《新疆地區(qū)地震分析預(yù)測技術(shù)方法工作手冊》和《地震電磁分析預(yù)報(bào)方法清單》等資料，目前以上資料已應(yīng)用到地震電磁異常日常震情跟蹤分析工作^{［12，16-19］}。中國地震局地震短臨跟蹤電磁學(xué)科組和新疆電磁學(xué)科組分別依據(jù)《地震電磁分析預(yù)測技術(shù)方法工作手冊》和《新疆地區(qū)地震分析預(yù)測技術(shù)方法工作手冊》新疆維吾爾自治區(qū)地震局：新疆地震分析預(yù)測技術(shù)方法手冊.2019.，對新疆地區(qū)2023年度上半年的電磁場資料進(jìn)行跟蹤分析，并據(jù)此對2023年上半年新疆地區(qū)震情形勢作出綜合預(yù)測。2023年1月30日新疆沙雅縣發(fā)生M_S6.1地震，中國地震臺網(wǎng)中心測定的微觀震中為82.29°E，40.01°N，震源深度為50 km，此次地震與上述預(yù)測結(jié)果相吻合。

1 震前異常及提取技術(shù)概述

地磁日變化是由太陽風(fēng)在電離層引起的環(huán)狀電流及其在地殼內(nèi)的感應(yīng)電流合并而成，地磁內(nèi)源場是由外源場在上地幔和地殼內(nèi)的臨時感應(yīng)電流產(chǎn)生的感應(yīng)場。因此，地磁垂直分量（Z）日變化畸變是由地震孕育過程中地殼電性結(jié)構(gòu)變化引起的^{［1，5，10］}。地磁低點(diǎn)位移法、地磁逐日比法是提取地磁日變化畸變的有效方法，在不同程度上可以圈定地殼和上地幔中的高導(dǎo)帶的走向。地磁垂直強(qiáng)度極化法是利用地磁場三分量秒采樣數(shù)據(jù)，計(jì)算地磁垂直分量與水平分量譜幅度之比的高值。異常機(jī)理為一次源在地殼上、周期在1 s左右，垂直強(qiáng)度大于水平強(qiáng)度的磁場信號，利用二者的比值既可突出巖石圈異常信號，也可抑制外源場的電磁信號。由于孕震區(qū)的電磁輻射為來自地殼的一次源電磁信號，因此該方法可以提取震前孕震區(qū)的電磁輻射信號，異常特征為極化高值超閾值^{［12，20-21］}。

新疆沙雅縣發(fā)生M_S6.1地震前，運(yùn)用上述方法對新疆及周邊地區(qū)電磁觀測數(shù)據(jù)（部分地磁方法采用全國地磁數(shù)據(jù)），按照手冊中的異常判據(jù)，進(jìn)行時空掃描處理分析^{［8，12］}，梳理并認(rèn)定新疆及其周邊地區(qū)存在7項(xiàng)地震前兆異常，其中地磁垂直強(qiáng)度極化3項(xiàng)、流動地磁1項(xiàng)、地磁逐日比１項(xiàng)、地磁低點(diǎn)位移2項(xiàng)（表1）。

2 異常分析與地震預(yù)測

2.1 地磁垂直強(qiáng)度極化異常

按照地震電磁分析預(yù)測技術(shù)方法工作手冊、《中國地震局電磁學(xué)科推薦分析方法清單（2020）》和異常判定規(guī)則，地磁垂直強(qiáng)度極化類異常屬性為區(qū)域電磁異常，異常項(xiàng)目為地震地磁場擾動異常。其異常判定依據(jù)為多個地磁臺站地磁垂直強(qiáng)度極化值出現(xiàn)超2倍均方差，并同時滿足以下5個條件時則視為1次異常：① 異常持續(xù)3 d以上；② 全國1/5以上地磁臺站出現(xiàn)準(zhǔn)同步高值異常；③ 歸一置零極化值高于0.2的面積大于6×10⁴km²；④ 2個月內(nèi)出現(xiàn)的異常區(qū)有重疊的視為同一組異常；⑤ 強(qiáng)震后40 d內(nèi)出現(xiàn)且該強(qiáng)震震中或附近有異常的異常為震后效應(yīng)。統(tǒng)計(jì)中國110°E以西區(qū)域發(fā)生的M_S≥6.0地震，地磁垂直強(qiáng)度極化異常預(yù)報(bào)效能R=0.462，R₀=0.317，通過置信度97.5%的效能檢驗(yàn)。新疆及周邊地區(qū)2015年1月1日以來，共有16組異常，其中3組屬于震后效應(yīng)，因此共有13組異常。異常發(fā)生后在預(yù)測區(qū)內(nèi)或附近地區(qū)發(fā)生M_S≥6.0地震7組，對應(yīng)率64%。根據(jù)以往震例及地震電磁分析方法預(yù)測技術(shù)手冊，極化異常拐點(diǎn)、重疊和反復(fù)出現(xiàn)異常的區(qū)域邊緣發(fā)震概率相對較高^{［12，21-26］}。

沙雅M_S6.1地震前新疆地區(qū)出現(xiàn)3次地磁垂直強(qiáng)度極化異常（圖1）：2022年5月8日新疆南部至甘青新交界地區(qū)出現(xiàn)1次地磁垂直強(qiáng)度極化異常，異常面積約為114.6×10⁴km²；2022年7月15日新疆東部至甘青新藏交界地區(qū)出現(xiàn)一次地磁垂直強(qiáng)度極化異常，異常面積約為222×10⁴km²；2022年11月8日新疆天山中段地區(qū)出現(xiàn)一次地磁垂直強(qiáng)度極化異常，異

常面積約為12×10⁴km²。以上3次地磁垂直強(qiáng)度極化異常，經(jīng)異?，F(xiàn)場核實(shí)認(rèn)為該異常信度較高，空間分布連續(xù)性較好，異常屬性為區(qū)域性地磁擾動類地震前兆異常，在震中周圍區(qū)域反復(fù)和重疊出現(xiàn)（圖1d）。沙雅M_S6.1地震發(fā)生在多次異常重疊區(qū)域附近，與預(yù)測目標(biāo)地震“時空強(qiáng)”基本吻合。

2.2 流動地震地磁異常

對2022年天山及周邊地區(qū)95個流動地震地磁觀測點(diǎn)數(shù)據(jù)采取日變通化、長期變數(shù)據(jù)改正、巖石圈磁場變化剝離等處理后，排除孤立測點(diǎn)異常和地磁單要素異常，綜合分析數(shù)據(jù)與構(gòu)造分區(qū)、構(gòu)造帶、斷層、地震活動性的關(guān)系，認(rèn)為2021—2022年阿克蘇―拜城地區(qū)出現(xiàn)了較顯著的巖石圈磁場局部異常變化，符合手冊中的異常Z矢量、H矢量、Ｋ圖、dD、dF、dI、dY、dZ、︱H︱、︱F︱的判定依據(jù)（圖2），根據(jù)預(yù)報(bào)規(guī)則震級大于6級屬于區(qū)域性中期異常^［12］。沙雅M_S6.1地震距離2022年度阿克蘇―拜城流動地磁異常區(qū)域120 km，已初步對應(yīng)本次地震。

2.3 地磁日變化類異常

2.3.1 地磁逐日比異常

地磁逐日比計(jì)算所采用數(shù)據(jù)為中國地磁臺網(wǎng)中心連續(xù)地磁觀測數(shù)據(jù)，其值為地磁垂直分量（Z）日變化幅度前后兩天之比的比值。新疆地區(qū)主要采用新疆及周邊地區(qū)地磁垂直分量（Z）數(shù)據(jù)計(jì)算地磁逐日比并開展日常震情跟蹤工作（圖3），屬于區(qū)域短期異常。異常判據(jù)：① 新疆地區(qū)地磁逐日比超過閾值2.6；② 剔除異常日前2天Dst小于（-70） nT或單臺K指數(shù)大于7；③ 兩次異常間隔時間2個月內(nèi)算為一組異常；④ 異常結(jié)束后有地震發(fā)生視為結(jié)束；⑤ 以最后一次異常日算為起點(diǎn)計(jì)算發(fā)震時間。2019年12月9日新疆天山中西段地區(qū)多臺站出現(xiàn)地磁日變幅逐日比值超閾值的現(xiàn)象，閾值線高曲率段位于庫車—阿克蘇（圖3），震例回溯統(tǒng)計(jì)顯示，2008年至今，已累計(jì)出現(xiàn)10次異常，其中有7次地震對應(yīng)，對應(yīng)率達(dá)到70%。優(yōu)勢發(fā)震期為前兩個月和異常后5個月。按照預(yù)測規(guī)則，一般情況下，地磁逐日比超過閾值后9個月內(nèi)，優(yōu)勢發(fā)震時間6個月內(nèi)，對應(yīng)率達(dá)到70%，閾值線高曲率（閾值線外凸、尖銳、內(nèi)凹）段存在發(fā)生M_S≥6.0地震的可能，震級與超閾值臺站的分布范圍呈正相關(guān)^［8］，預(yù)測效能R=0.35，

R₀=0.252，通過置信度97.5%的效能檢驗(yàn)，異常信度為A類^{［3，26］}。本次逐日比異常為震前提出，異常后有效期內(nèi)在預(yù)測區(qū)附近發(fā)生2023年1月30日沙雅M_S6.1地震（距閾值線約250 km）。分析認(rèn)為，該異?？梢詫?yīng)沙雅M_S6.1地震。

2.3.2 地磁低點(diǎn)位移異常

地磁低點(diǎn)位移法是使用中國地磁臺網(wǎng)地磁垂直分量（Z）（國際時）分鐘采樣資料計(jì)算日變化極小時間（低點(diǎn)時間），在地圖對應(yīng)臺站上標(biāo)注時間，然后判定異常分界線。如果存在兩個或兩個以上區(qū)域，在各自區(qū)域內(nèi)低點(diǎn)時間基本一致，而各區(qū)域之間低點(diǎn)時間存在2小時以上差距，形成明顯的低點(diǎn)時間突變分界線，則判斷出現(xiàn)低點(diǎn)位移異常，屬于地磁場日變化類區(qū)域異常^［1］。

2022年11月26日中國大陸出現(xiàn)低點(diǎn)位移異常，近東西向穿越新疆中部；2023年1月9日中國大陸出現(xiàn)低點(diǎn)位移異常，異常圖分界線突變穿過新疆和田至阿克蘇一帶，形成一個弧形異常區(qū)，為短期區(qū)域性異常，前后兩條低點(diǎn)位移分界線在天山中段重疊（圖4）。預(yù)測時段為異常后60天，優(yōu)勢預(yù)測時段為2022年2月5日和2月19日前后4天；

預(yù)測區(qū)域?yàn)楫惓Ｍ蛔兎纸缇€附近300 km范圍內(nèi)；預(yù)測震級6級左右，預(yù)測效能R=0.179，R₀=0.109，通過置信度97.5%的效能檢驗(yàn)，異常信度A類^［12］?？紤]到地磁逐日比和低點(diǎn)位移前后超過閾值的臺站分布范圍不大，判定閾值線高曲率段所處的新疆阿克蘇―庫車一帶在2023年上半年存在發(fā)生M_S6.0左右地震的可能。

本次低點(diǎn)位移異常為震前提出，異常后預(yù)測區(qū)內(nèi)60天內(nèi)發(fā)生了2023年1月30日沙雅M_S6.1地震（距閾值線約250 km）。分析認(rèn)為，該異?？梢詫?yīng)沙雅M_S6.1地震。

3 電磁異常綜合分析

3.1 異常時空分布特征分析

沙雅M_S6.1地震前震中300 km范圍內(nèi)定點(diǎn)電磁臺站柯坪地電阻率、烏什地電場震前均無異常，該區(qū)域?qū)儆谟^測站點(diǎn)稀疏，監(jiān)測能力較弱地區(qū)，但區(qū)域異常顯著，空間短臨異常較為突出，這是與以往震例明顯不同之處。震前發(fā)現(xiàn)7項(xiàng)電磁中短期異常（表1），其中地磁垂直強(qiáng)度極化和流動地磁異常等4項(xiàng)為中期異常，流動地磁和2022年5月8日、7月15日地磁垂直強(qiáng)度極化異常為年度異常，其他3項(xiàng)次均為短期異常。以上7項(xiàng)中短期異常疊加區(qū)域位于天山中西段，沙雅M_S6.1地震發(fā)生在異常疊加區(qū)附近；地磁低點(diǎn)位移異常2022年11月26日和2023年1月9日穿越天山中段，且呈現(xiàn)低點(diǎn)位移異常閾值線集中分布的現(xiàn)象，優(yōu)勢對應(yīng)期為2個月，沙雅M_S6.1地震發(fā)生在異常附近地區(qū)；2022年12月9日逐日比異常，優(yōu)勢對應(yīng)期為震后2個月和第5～6個月，異?？臻g范圍較大，但是庫車―輪臺―拜城地區(qū)（以下簡稱庫輪拜地區(qū)）異常閾值線呈現(xiàn)外突尖或內(nèi)凹，地震發(fā)生在異常突尖或內(nèi)凹凸區(qū)附近，判定時空強(qiáng)指示意義明顯。以上7項(xiàng)電磁異常綜合判定與突發(fā)地震對比情況見圖5。

綜合分析認(rèn)為，在M_S6.1地震發(fā)生前，電磁異常具有中期逐漸向短（臨）期逼近的過程。此外本次地震發(fā)生在區(qū)域地球物理異常疊加交匯處，提供了后續(xù)強(qiáng)震的發(fā)生地點(diǎn)和時間的判定依據(jù)。

3.2 綜合判定結(jié)果與實(shí)發(fā)地震對比

震前電磁手段共出現(xiàn)了7項(xiàng)中短期異常，從年度學(xué)科會商的3項(xiàng)異常逐漸增加到7項(xiàng)，這種地球物理觀測資料群體性快速增加，震前均有一定程度的短期預(yù)測。綜合考慮新疆中部地區(qū)地震活動性和地震構(gòu)造背景，并結(jié)合每一項(xiàng)異常的預(yù)測結(jié)論，對新疆地區(qū)地震形勢進(jìn)行綜合預(yù)測，認(rèn)為2023年上半年天山中段尤其是庫輪拜地區(qū)可能發(fā)生M_S6.0左右地震的危險(xiǎn)性（圖6）。2023年1月30日新疆沙雅發(fā)生M_S6.1地震，其發(fā)震時間和強(qiáng)度與預(yù)測意見完全一致，震中位于預(yù)測區(qū)邊緣（圖5），與綜合預(yù)測意見稍有偏差，綜合預(yù)測時間進(jìn)程見圖6。

4 結(jié)論與討論

地球磁場中各種成分的變化，能靈敏地反映地球內(nèi)外部空間發(fā)生的與電磁過程有關(guān)的各種物理過程，因此中強(qiáng)地震前存在電磁異?，F(xiàn)象，并且已被國內(nèi)外大量震例所證實(shí)。2023年1月30日新疆沙雅M_S6.1地震發(fā)震前，電磁學(xué)科組共監(jiān)測到7項(xiàng)異常，并提出了明確的預(yù)測意見，為電磁學(xué)科日常震情跟蹤及地震監(jiān)測預(yù)報(bào)提供了較好的應(yīng)用案例。此次新疆沙雅M_S6.1強(qiáng)震是對電磁學(xué)科指標(biāo)體系建設(shè)成果的有效檢驗(yàn)，《地震電磁分析預(yù)測技術(shù)方法工作手冊》在此次地震預(yù)測和判定過程中發(fā)揮著重要作用。

縱觀此次沙雅M_S6.1地震電磁異常的時間進(jìn)程：2021—2022年阿克蘇―拜城地區(qū)出現(xiàn)了比較顯著的巖石圈磁場局部異常變化（圖2），預(yù)示著該區(qū)域存在強(qiáng)震背景；2022年5月8日和7月15日新疆地區(qū)出現(xiàn)較大范圍的地磁垂直強(qiáng)度極化異常并與庫輪拜地區(qū)兩次異常線交匯，預(yù)示著異常交匯區(qū)域有發(fā)生M_S≥6.1地震的可能，發(fā)震時間預(yù)測進(jìn)程上則進(jìn)入中期；2022年11月8日天山中段區(qū)域出現(xiàn)地磁垂直強(qiáng)度極化異常，極化異常區(qū)明顯縮小，異常區(qū)集中在天山中西段區(qū)域，預(yù)示著發(fā)震概率的進(jìn)一步增加；2022年12月9日天山中西段出現(xiàn)地磁逐日比異常，發(fā)震時間預(yù)測進(jìn)程進(jìn)入短期階段；2022年11月26日和2023年1月9日地磁低點(diǎn)位移出現(xiàn)異常且在天山中段呈現(xiàn)重疊區(qū)，預(yù)示著發(fā)震臨震階段。2023年1月30日發(fā)生沙雅M_S6.1地震，時空演化特征較明顯，即距離發(fā)震時間越近，異常數(shù)量明顯增加，預(yù)示著該地區(qū)的發(fā)震概率增大（圖6）。2017年精河M_S6.6、2018年精河M_S5.4和2020年伽師M_S6.4地震前也出現(xiàn)過類似地磁異常變化^［20-22］，有一定相似性，表明庫輪拜地區(qū)在此次沙雅M_S6.1地震震前7～8個月左右開始呈現(xiàn)出介質(zhì)電性結(jié)構(gòu)的改變^［10-11］，體現(xiàn)了多方法綜合預(yù)測的優(yōu)勢。眾所周知，電磁單方法異常在時間軸上均為一段時間的持續(xù)變化，一般較難判定地震是否進(jìn)入臨近時段或發(fā)震危險(xiǎn)性在時間軸上的變化，同時還需考慮發(fā)震地點(diǎn)和強(qiáng)度預(yù)測不確定等因素。而電磁綜合預(yù)測可有效縮小預(yù)測“時空強(qiáng)”誤差范圍。

地磁逐日比、低點(diǎn)位移、地磁垂直強(qiáng)度極化、流動地磁等地磁類異常主要反映了上地幔和地殼內(nèi)高導(dǎo)帶附近的地震孕震體信息，即異常曲線或閾值線跟上地幔和地殼中高導(dǎo)帶走向基本一致^［7］，反映了上地幔和地殼中高導(dǎo)帶的穩(wěn)定狀態(tài)。2022年12月9日新疆天山地區(qū)多個臺站地磁逐日比超過閾值，閾值線高曲率段正好在震中周邊；2022年11月26日、2023年1月9日新疆地區(qū)出現(xiàn)低點(diǎn)位移異常閾值線集中分布的現(xiàn)象，2022年5月8日、7月15日和11月8日地磁垂直強(qiáng)度極化異常閾值線在流動地磁異常區(qū)重疊或交匯，其中一處異常集中線位于新疆庫輪拜地區(qū)。此次沙雅M_S6.1地震震源體正好位于殼內(nèi)高導(dǎo)層附近，震前日變化類異常的集中出現(xiàn)很大程度上反映出庫輪拜老震區(qū)震前7～8個月介質(zhì)電性結(jié)構(gòu)變化呈現(xiàn)加速態(tài)勢。

雖然電磁場異常是地震孕育和發(fā)生過程中的一種伴生現(xiàn)象，但并不是所有中強(qiáng)地震孕震和發(fā)生過程中都會出現(xiàn)電磁異?，F(xiàn)象^{［3，12，27］}。馮志生等^［28］研究發(fā)現(xiàn)地磁場日變化類異常其原始數(shù)據(jù)在日變化形態(tài)上會出現(xiàn)不同程度的“反相位”或“畸變”變化，基于電磁場理論推測“反相位”或“畸變”變化期間，有感應(yīng)電流集中分布于分界線下方，感應(yīng)電流集中分布會導(dǎo)致地磁垂直分量（Z）日變化幅度變大或變?。坏卮糯怪睆?qiáng)度極化類異常是由震前斷層蠕動或破裂產(chǎn)生的高頻電磁輻射所致。但短時感應(yīng)磁場的復(fù)雜性也會導(dǎo)致異常時段地磁場日變相位和幅度異常的表現(xiàn)形式多樣化。比如地磁低點(diǎn)位移法異常，在相位變化上可能更為敏感^［1，3］，而地磁逐日比法異常則在幅度變化上更為敏感^［12］，地磁垂直強(qiáng)度極化則更為全面^［24］，這也可以解釋不同顯著地震事件前捕捉到不同的電磁異常。

變化磁場地震異常的源電流定位是震前發(fā)生活動異常構(gòu)造的位置^［28］，因此理論上，依據(jù)電磁異常來預(yù)測地震時，可以將電磁各方法異常簡單的綜合疊加估計(jì)震中位置。但是，實(shí)際工作中震前異常識別受諸多因素影響，還需不斷增強(qiáng)開展震前異常產(chǎn)生的物理機(jī)制研究，進(jìn)一步提高震中判定精確度。

致謝：感謝福建局倪曉寅高級工程師、中國地震臺網(wǎng)中心姚麗高級工程師、王亞麗高級工程師給予本文的幫助和指導(dǎo)。本文地磁數(shù)據(jù)來自于中國地震局地球物理研究所國家地磁臺網(wǎng)中心，地震目錄來自新疆地震臺網(wǎng)中心，在此一并致謝。

參考文獻(xiàn)：［1］ 丁鑒海，車時，余素榮，等.地磁日變地震預(yù)報(bào)方法及其震例研究［M］.北京：地震出版社，2009：1-3.

［2］ 黃頌，李軍輝，姚麗，等.基于日變化相關(guān)的地磁低點(diǎn)位移異?？煽啃苑治龇椒ǎ跩］.中國地震，2021，37（4）：780-788.

［3］ 艾薩·伊斯馬伊力，毛志強(qiáng)，馮麗麗.2018年精河5.4級地震前地磁日變化畸變異常分析［J］.地震工程學(xué)報(bào)，2020，42（2）：368-375.

［4］ 艾薩·伊斯馬伊力，馮麗麗.新疆南天山西段地區(qū)地磁加卸載響應(yīng)比分析［J］.國際地震動態(tài)，2019，（8）：62-63.

［5］ 戴苗，馮志生，劉堅(jiān)，等.2020-07唐山M_S5.1地震前電磁資料綜合分析［J］.大地測量與地球動力學(xué)，2021，41（11）：1 146-1 151.

［6］ 艾薩·伊斯馬伊力，高麗娟.新疆塔什庫爾干M_S5.5地震前兆異常分析［J］.地震地磁觀測與研究，2018，39（5）：81-86.

［7］ 艾薩·伊斯馬伊力，黃恩賢，高麗娟.2017年精河M_S6.6地震前地磁異常特征分析［J］.中國地震，2017，33（4）：764-770.

［8］ 艾薩·伊斯馬伊力.喀什臺地磁異常與地震［J］.華南地震，2017，37（4）：54-60.

［9］ 艾薩·伊斯馬伊力，張治廣，趙彬彬，等.新疆阿克陶M_S6.7地震前喀什地磁資料變化分析與探討［J］.內(nèi)陸地震，2017，31（2）：213-221.

［10］馮志生，李琪，李鴻宇，等.地磁低點(diǎn)位移線兩側(cè)異常變化的反相位現(xiàn)象及其解釋［J］.中國地震，2009，25（2）：206-213.

［11］艾薩·伊斯馬伊力，馮志生，毛志強(qiáng)，等.新疆伽師M_S6.4地震前地磁日變化異常分析［J］.內(nèi)陸地震，2020，34（1）：79-86.

［12］中國地震局監(jiān)測預(yù)報(bào)司.地震電磁分析預(yù)測技術(shù)方法工作手冊［M］.北京：地震出版社，2020：117-143.

［13］李鴻宇，袁桂平，楊云，等.基于速率累加算法的地磁諧波振幅比異常分析［J］.防災(zāi)減災(zāi)工程學(xué)報(bào)，2019，39（6）：1 010-1 015.

［14］艾薩·伊斯馬伊力，馮志生，陳界宏，等.2020年新疆于田M_S6.4地震前地磁諧波振幅比變化特征分析［J］.地震，2021，41（2）：190-202.

［15］艾薩·伊斯馬伊力，陳界宏，次旦，等.和田臺地電場巖體裂隙優(yōu)勢方位角變化特征分析［J］.內(nèi)陸地震，2021，35（3）：245-252.

［16］馮志生，倪曉寅，戴苗，等.地震地磁日變化預(yù)測指標(biāo)體系建設(shè)進(jìn)展［J］.國際地震動態(tài)，2019，（8）：11-12.

［17］馮志生.地磁日變化畸變地震預(yù)報(bào)指標(biāo)研究進(jìn)展［J］.國際地震動態(tài)，2018，（8）：17-18.

［18］馮志生，王俊菲，沈紅會，等.地震短期綜合預(yù)報(bào)的一種思路［J］.國際地震動態(tài)，2017，（8）：102-103.

［19］馮志生，李琪，李鴻宇，等.地磁低點(diǎn)位移線兩側(cè)異常變化的反相位現(xiàn)象及其解釋［J］.中國地震，2009，25（2）：206-213.

［20］Hayakawa M，Kawate R，Molchanov O A，et al.Results of ultra-low-frequency magnetic field measurements during the Guam Earthquake of 8 August 1993［J］.Geophysical Research Letters，1996，23（3）：241-244.

［21］馮志生，李琪，盧軍，等.基于磁通門秒數(shù)據(jù)的ULF磁場可靠信息提?。跩］.華南地震， 2010，30（2）：1-7.

［22］艾薩·伊斯馬伊力，馮志生，陳界宏，等.2020年于田M_S6.4地震前地磁極化異常特征分析［J］.內(nèi)陸地震，2020，34（3）：295-302.

［23］馮麗麗，管貽亮，樊文杰，等.2020年10月15日地磁垂直強(qiáng)度極化異常與瑪多M_S7.4地震［J］.高原地震，2021，33（4）：1-6.

［24］李霞，馮麗麗，趙玉紅，等.2017年西藏米林M_S6.9地震前地磁垂直強(qiáng)度極化異常特征［J］.地震學(xué)報(bào)，2021，43（5）：584-594+678.

［25］馮麗麗，馮志生，樊文杰，等.中國大陸西部強(qiáng)震前超低頻磁場時空變化特征［J］.地震學(xué)報(bào)，2021，43（3）：359-375.

［26］樊文杰，馮麗麗，李霞，等.中強(qiáng)地震前地磁極化異常特征及其與后續(xù)地震的關(guān)系［J］.地震，2021，41（2）：170-179.

［27］艾薩·伊斯馬伊力，馮志生，陳界宏，等.地磁逐日比與新疆及周邊地區(qū)強(qiáng)震關(guān)系［J］.空間科學(xué)學(xué)報(bào)，2021，41（4）：617-625.

［28］倪曉寅，洪旭瑜.華南地區(qū)地磁預(yù)測方法預(yù)測指標(biāo)體系的初步應(yīng)用［J］.華南地震，2022，42（4）：60-66.

［29］馮志生，姜楚峰，馮麗麗，等.短期重現(xiàn)性地磁日變化感應(yīng)電流集中分布與地震關(guān)系初步研究［J］.中國地震，2020，36（3）：502-516.

XINJIANG APPLICATION OF ELECTROMAGNETIC

PREDICTION INDEX SYSTEM IN SHAYA M_S6.1

EARTHQUAKE ON JANUARY 30， 2023

Aisa Yisimayili¹，? ZHAI Shi-long¹，?? ZHANG An-he¹，? FENG Li-li²，

HE Man-qiu³，? FENG Zhi-sheng⁴

（1.Earthquake Agency of Xinjiang Uygur Autonomous Region，Urumqi 830011，Xinjiang，China;

2. Qinghai Earthquake Agency， Xining 810001，Qinghai， China;

3.Chongqing Earthquake Agency，Chongqing 401147，China;4.Jiangsu Earthquake Agency， Nanjing 810001，Jiangsu， China;）

Abstract： On January 30， 2023， the M_S6.1 earthquake with focal depth of 50 km occurred in Xayar County （82.29°E， 40.01°N）， Xinjiang. Before the earthquake， the electromagnetic discipline group of short-term and imminent earthquake prediction tracking of China Earthquake Administration and the electromagnetic discipline group of Earthquake Agency of Xinjiang Uygur Autonomous Region sorted out the anomalies of the data of electromagnetic stations in Xinjiang according to Manual of Seismic Electromagnetic Analysis and Prediction Technical Methods， Manual of Seismic Analysis and Prediction Technical Methods in Xinjiang and List of Seismic Electromagnetic Analysis and Prediction Methods， and found that there were mobile geomagnetism multiple anomalies such as geomagnetic vertical intensity polarization， daily geomagnetic ratio， and geomagnetic low point displacement. The comprehensive prediction opinions of anomalies before earthquakes are basically consistent with the actual earthquakes. The evolution characteristics of medium and short-term progressive anomalies before earthquakes are prominent. The construction of the prediction index system of seismoelectromagnetic discipline in Xinjiang from 2017 to 2020 has achieved initial results.

Key words： Shaya M_S6.1 earthquake; Electromagnetic anomaly; Anomaly index; Earthquake prediction