謝江麗 阿布都瓦里斯?阿布都瓦衣提 張文秀

摘要：利用1920—2020年地震資料，對(duì)柴達(dá)木—阿爾金地震帶采用多方案（不同時(shí)間段、震級(jí)段）統(tǒng)計(jì)計(jì)算并綜合分析，確定b值為0.78，較《中國(guó)地震動(dòng)參數(shù)區(qū)劃圖》中確定的b值0.84有所降低，對(duì)地震帶重大建設(shè)工程地震安全性評(píng)價(jià)與地震小區(qū)劃的工程設(shè)防有重要意義。本文確定b值作為《中國(guó)地震動(dòng)參數(shù)區(qū)劃圖》的補(bǔ)充，可為科學(xué)建立重大建設(shè)工程的地震危險(xiǎn)性計(jì)算模型提供基礎(chǔ)依據(jù)。

關(guān)鍵詞：柴達(dá)木—阿爾金地震帶；b值；多方案；地震動(dòng)參數(shù)

doi：10.16256/j.issn.1001-8956.2023.03.010

b值的統(tǒng)計(jì)結(jié)果受時(shí)空范圍、震級(jí)上限和震級(jí)的不確定性等多種因素影響^{［1，2］}，選擇不同的計(jì)算方法結(jié)果也不同^［3］。常用的方法為最大似然法和最小二乘法，其中最大似然法對(duì)樣本數(shù)量要求比較高^［4］，利用最小二乘法合理估計(jì)b值已經(jīng)相當(dāng)成熟，且計(jì)算結(jié)果也比較合理^［5］。本文研究區(qū)域歷史記錄時(shí)間較短，地震記錄較少，綜合實(shí)際資料情況，采用最小二乘法擬合b值比較準(zhǔn)確。利用柴達(dá)木—阿爾金地震帶1920—2020年地震資料，根據(jù)不同的統(tǒng)計(jì)時(shí)間段和震級(jí)段進(jìn)行計(jì)算，并用地震年發(fā)生率約束計(jì)算結(jié)果，綜合分析確定b值^［6］。此方法充分考慮了地震資料不完備性和認(rèn)識(shí)不確定性，同時(shí)利用歷史地震資料和現(xiàn)代儀器記錄地震資料，充分發(fā)掘有價(jià)值的信息，克服樣本資料不完備可能導(dǎo)致的認(rèn)識(shí)偏差，更好的體現(xiàn)柴達(dá)木—阿爾金地震帶的地震活動(dòng)水平。通過補(bǔ)充近年來地震資料重新統(tǒng)計(jì)計(jì)算b值，在下一代區(qū)劃圖編制實(shí)施前，可為科學(xué)建立重大建設(shè)工程地震危險(xiǎn)性計(jì)算模型提供基礎(chǔ)依據(jù)，對(duì)該地震帶重大建設(shè)工程地震安全性評(píng)價(jià)與地震小區(qū)劃工程均有重要意義。

1 柴達(dá)木—阿爾金地震帶地震資料統(tǒng)計(jì)分析

柴達(dá)木—阿爾金地震帶西南從瓊木孜塔格一帶起，東北至甘肅與新疆交界，由柴達(dá)木盆地和阿爾金兩個(gè)地震構(gòu)造區(qū)組成^［7］，總體走向NWW。該地震帶中部地區(qū)地震活動(dòng)水平較低，大震主要集中在地震帶的東段和西段。

通過對(duì)柴達(dá)木—阿爾金地震帶歷史地震資料收集整理，結(jié)合其地震活動(dòng)時(shí)空分布特點(diǎn)^［8］，使用多方案的統(tǒng)計(jì)方法對(duì)柴達(dá)木—阿爾金地震帶的b值進(jìn)行統(tǒng)計(jì)計(jì)算，汪素云等人對(duì)震級(jí)標(biāo)度轉(zhuǎn)換進(jìn)行了專門的研究，認(rèn)為M_L≈M_S^［9］，因此，本文在編制地震目錄時(shí)不再進(jìn)行面波震級(jí)與地方震級(jí)的轉(zhuǎn)換^［10］，且統(tǒng)一表示成M。

本文中采用的地震資料來源如下：《中國(guó)歷史強(qiáng)震目錄（公元前23世紀(jì)—公元1911年）》、《中國(guó)近代地震目錄（公元1912—1990年，M_S≥4.7）》、《中國(guó)地震臺(tái)網(wǎng)觀測(cè)報(bào)告》、《新疆維吾爾自治區(qū)地震目錄（1970—1984年）》和新疆地震臺(tái)網(wǎng)。

柴達(dá)木—阿爾金地震帶最早1次地震記錄為公元1832年8月昌馬5?級(jí)地震，歷史上最大震級(jí)為7.3，分別為2008年3月21日和2014年2月12日于田M7.3地震。由于柴達(dá)木—阿爾金地震帶歷史地震資料較少、時(shí)間記載較短，資料缺失嚴(yán)重，1920年以前僅記錄到4次地震。1920年以后M≥5.0地震記錄基本完整^{［11、12］}。截至2020年12月，該地震帶共發(fā)生901次M≥4.0地震，其中5次7級(jí)地震，24次6級(jí)地震，169次5級(jí)地震和703次4級(jí)地震。5次7級(jí)地震中4次發(fā)生在地震帶西段，包括時(shí)間最近的2014、2008年于田的兩次7.3級(jí)地震；1次7級(jí)地震發(fā)生在地震帶東段，地震帶西段強(qiáng)震發(fā)生率高于東段。柴達(dá)木—阿爾金地震帶M≥6.0 地震震源參數(shù)見表1。

2 地震帶地震目錄完整性及活動(dòng)趨勢(shì)分析

在研究地震活動(dòng)性時(shí)，地震目錄完整性尤其重要，其對(duì)地震活動(dòng)性參數(shù)計(jì)算結(jié)果產(chǎn)生較大影響。地震監(jiān)測(cè)能力在不同歷史時(shí)期的差異性較大，不少學(xué)者多次對(duì)新疆及相鄰地區(qū)地震目錄進(jìn)行完整性分析^［13］。為保證地震資料的可靠性、連續(xù)性和充分性^［14］，結(jié)合本次研究區(qū)域特點(diǎn)，使用統(tǒng)計(jì)基本資料為1920年以來柴達(dá)木—阿爾金地震帶內(nèi)M≥5.0歷史地震和1970年以來區(qū)域臺(tái)網(wǎng)記錄M≥4.0地震。

利用柴達(dá)木—阿爾金地震帶地震目錄繪制地震時(shí)序圖（圖1），隨著地震監(jiān)測(cè)手段逐步完善，能記錄到的地震最小完整性震級(jí)降低，地震分布越來越密集，間接說明歷史地震遺漏較多。另外地震活動(dòng)平靜與活躍交替出現(xiàn)，存在成叢性特點(diǎn)，在時(shí)間上分布很不均勻。從地震發(fā)生趨勢(shì)上看，本文認(rèn)為該地震帶正處于地震活躍期，在此基礎(chǔ)上進(jìn)行b值計(jì)算。

3 地震帶余震刪除

地震活動(dòng)性統(tǒng)計(jì)特征與余震關(guān)系密切，余震對(duì)地震獨(dú)立性有較大影響。假設(shè)震源之間具有相互獨(dú)立性，余震被認(rèn)為是和主震相關(guān)的事件，與主震之間不是獨(dú)立關(guān)系，如果余震去除的不完整，地震的不獨(dú)立性會(huì)使其不符合泊松分布。未刪除余震會(huì)造成地震活動(dòng)性統(tǒng)計(jì)分析不準(zhǔn)確^［15］，故需先刪除余震再對(duì)地震目錄進(jìn)行統(tǒng)計(jì)分析。目前刪除余震的方法已有許多種。利用M（主震震級(jí)）的經(jīng)驗(yàn)公式（M=3.3+2.1logR）與R（km）斷層長(zhǎng)度確定M（主震震級(jí)）與r（km）余震刪除空間半徑的關(guān)系為^［16］

logr=0.48M-1.57 .（1）

式中：r是余震刪除空間半徑，M是主震震級(jí)。根據(jù)地震實(shí)際定位精度，取r=5 km作為空間窗半徑的最小值。采用陳凌等人研究的不同震級(jí)余震時(shí)間范圍^［17］，確定4.0≤M＜4.5余震時(shí)間為42天，4.5≤M＜5.0為83天，5.0≤M＜5.5為155天，5.5≤M＜6.0為290天，6.0≤M＜6.5為510天，6.5≤M＜7.0為790天。

從本次統(tǒng)計(jì)基本資料選取2008年3月21日至2009年3月21日內(nèi)記錄的地震，剔除不在此區(qū)域范圍內(nèi)的地震，繪制于田7.3級(jí)地震余震分布圖^［18］（圖2），根據(jù)余震的空間半徑r與主震震級(jí)M的關(guān)系式得出本次地震余震空間半徑r=85 km，由圖2可以看出，余震基本都在此范圍內(nèi)，證明余震分布符合此關(guān)系式。

本文選取了柴達(dá)木—阿爾金地震帶1920—2020年12月M≥4.0地震目錄為基礎(chǔ)資料，余震空間窗和時(shí)間窗進(jìn)行掃描，手動(dòng)刪除余震，同時(shí)進(jìn)行刪除前后地震數(shù)目對(duì)比（表2），本次研究選取的數(shù)據(jù)資料為刪除余震后的地震目錄。

4 地震帶b值計(jì)算方法與結(jié)果

在概率地震危險(xiǎn)性分析方法中，地震帶活動(dòng)性參數(shù)包括地震帶震級(jí)上限M_uz、起算震級(jí)M₀、地震年發(fā)生率V₀及震級(jí)頻度關(guān)系式中的b值^{［19，20］}，其統(tǒng)計(jì)表達(dá)式為

lgN=a-bM .（2）

式中：N為不同震級(jí)檔地震發(fā)生的次數(shù)；M為地震震級(jí)；a、b為常數(shù)。

V₀是指在一定范圍內(nèi)（如地震帶）平均每年發(fā)生震級(jí)等于或大于M₀（通常用M=4）的地震數(shù)，代表了統(tǒng)計(jì)區(qū)內(nèi)的地震活動(dòng)水平。V₀在本文中作為確定b值的一個(gè)約束參數(shù)。

為考慮地震資料的完備性，根據(jù)起算震級(jí)和地震完整性可信時(shí)段的不同選取多種統(tǒng)計(jì)方案，同時(shí)考慮每個(gè)起算震級(jí)的地震完整可信時(shí)段樣本量是最充分的，也為避免忽略復(fù)發(fā)周期長(zhǎng)的大地震對(duì)b值的影響，選取了1950—2020年M≥5.5、M≥5.0地震和1970—2020年M≥4.0、 M≥4.5等代表整個(gè)地震帶可信時(shí)間段的統(tǒng)計(jì)樣本。利用b值統(tǒng)計(jì)關(guān)系式和最小二乘法得出b值擬合曲線，圖3為部分樣本選取示意圖，圖4為部分b值擬合曲線圖。

根據(jù)上述選取樣本計(jì)算b值和其對(duì)應(yīng)的地震年發(fā)生率，忽略因地震不完備性造成的b值過高或過低的樣本，可信樣本為33個(gè)（表3），同時(shí)對(duì)b值和地震年發(fā)生率結(jié)果進(jìn)行散點(diǎn)繪圖（圖5）。

如圖5（a）所示，b值結(jié)果處于0.47～0.81范圍內(nèi)，基本符合震級(jí)越高b值越低的規(guī)律。如圖5（b）所示，不同震級(jí)的地震年發(fā)生率均在一定范圍內(nèi)變化，震級(jí)越高，其地震年發(fā)生率越低。75%～85%范圍內(nèi)的V₄所對(duì)應(yīng)b值都在0.78附近，綜合考慮b值結(jié)果和V₄，b值取0.78與之最接近V₄所對(duì)應(yīng)b值都在0.78附近，故b值取0.78。該結(jié)果低于《中國(guó)地震動(dòng)參數(shù)區(qū)劃圖》，可為判定該地區(qū)地震危險(xiǎn)性提供參考。

5 結(jié)束語

利用現(xiàn)有地震資料，對(duì)柴達(dá)木—阿爾金地震帶b值進(jìn)行統(tǒng)計(jì)研究，采用多方案統(tǒng)計(jì)計(jì)算得到b值0.78，低于《中國(guó)地震動(dòng)參數(shù)區(qū)劃圖》b值0.84，b值的降低，影響該地震帶重大建設(shè)工程地震安全性評(píng)價(jià)與地震小區(qū)劃工程設(shè)防參數(shù)的調(diào)整。針對(duì)地震記錄時(shí)間短，歷史地震資料中、小地震不全的地震帶，本文所用多方案統(tǒng)計(jì)方法克服了樣本資料不完備性可能導(dǎo)致的認(rèn)識(shí)偏差，b值取值相對(duì)合理，歷史地震資料情況類似的阿爾泰山地震帶等可以用本方案計(jì)算b值。

參考文獻(xiàn)：

［1］ 許洪泰，閔偉.利用地震矩和實(shí)時(shí)概率模型對(duì)中長(zhǎng)期地震危險(xiǎn)性進(jìn)行評(píng)價(jià)——以東昆侖和阿爾金斷裂為例［J］.地震研究，2007，30（2）：120-126.

［2］ 謝卓娟，呂悅軍，蘭景巖，等.b值和V₄的統(tǒng)計(jì)分析及其不確定性對(duì)地震危險(xiǎn)性分析結(jié)果的影響研究［J］.地震研究，2013，36（1）：86-92.

［3］ 鄭確，劉財(cái)，田有，等.地震活動(dòng)性中震級(jí)-頻度關(guān)系研究進(jìn)展與再認(rèn)識(shí)［J］.地球物理學(xué)進(jìn)展，2018，33（5）：1 879-1 889.

［4］ 吳果，周慶，冉洪流.震級(jí)-頻度關(guān)系中b值的極大似然法估計(jì)及其影響因素分析［J］.地震地質(zhì)，2019，41（1）：21-43.

［5］ 常金龍.鶴崗及周邊地區(qū)b值計(jì)算和地震危險(xiǎn)性分析［J］.防災(zāi)減災(zāi)學(xué)報(bào)，2017，33（1）：32-35.

［6］ 潘華，高孟潭，謝富仁.新版地震區(qū)劃圖地震活動(dòng)性模型與參數(shù)確定［J］.震災(zāi)防御技術(shù)，2013，8（1）：11-23.

［7］ 姚遠(yuǎn)，伊力亞爾·阿不力孜，吳傳勇，等.2014年2月12日新疆于田M_S7.3地震發(fā)震構(gòu)造初步研究［J］.內(nèi)陸地震，2014，28（2）：121-126.

［8］ 謝江麗，阿布都瓦里斯·阿布都瓦衣提，黃帥堂.柴達(dá)木—阿爾金地震帶b值統(tǒng)計(jì)計(jì)算［J］.震災(zāi)防御技術(shù)，2022，17（1）：95-103.

［9］ 汪素云，俞言祥.震級(jí)轉(zhuǎn)換關(guān)系及其對(duì)地震活動(dòng)性參數(shù)的影響研究［J］.震災(zāi)防御技術(shù)，2009，4（2）：141-149.

［10］金花，閆新義，任林，等.新疆地區(qū)近震震級(jí)與面波震級(jí)經(jīng)驗(yàn)關(guān)系初步研究［J］.內(nèi)陸地震，2021，35（2）：150-156.

［11］單新建，韓京，許靜.新疆測(cè)震臺(tái)網(wǎng)歷史監(jiān)測(cè)能力及現(xiàn)狀［J］.內(nèi)陸地震，1996，10（1）：61-67.

［12］王海濤，李瑩甄，屠泓為.新疆歷史地震目錄完整性分析［J］.內(nèi)陸地震，2006，20（1）：10-17.

［13］徐偉進(jìn)，高孟潭.中國(guó)大陸及周緣地震目錄完整性統(tǒng)計(jì)分析［J］.地球物理學(xué)報(bào)，2014，57（9）：2 802-2 812.

［14］任雪梅，高孟潭，馮靜.地震目錄的完整性對(duì)b值計(jì)算的影響［J］.震災(zāi)防御技術(shù)，2011，6（3）：257-268.

［15］冉慧敏，上官文明，趙曉成.2020年新疆于田M_S6.4地震及余震序列定位研究［J］.內(nèi)陸地震，2020，34（3）：249-257.

［16］Keilis-Borok V I，Knopoff L，Rotvain I M.Bursts of aftershocks，long-term precursors of strong earthquakes［J］.Nature，1980，283：259-263.

［17］陳凌，陳颙，劉杰，等.地震活動(dòng)性的統(tǒng)計(jì)分析：由過去推測(cè)將來的可能性研究［J］.地球物理學(xué)報(bào)，1998，41（1）：61-70.

［18］楊成榮，唐明帥.新疆于田7.3級(jí)地震及部分強(qiáng)余震震源參數(shù)的重新測(cè)定［J］.內(nèi)陸地震，2010，24（1）：18-24.

［19］中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局，中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).GB 17741-2005 工程場(chǎng)地地震安全性評(píng)價(jià)［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2006.

［20］中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局，中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).GB 18306-2015 中國(guó)地震動(dòng)參數(shù)區(qū)劃圖［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2016.

ANALYSIS OF B VALUE IN QAIDAM ALTUN SEISMIC

BELT WITH MULTIPLE SCHEMES

XIE Jiang-li，? Abuduwalisi Abuduwayiti，? ZHANG Wen-xiu

（Earthquake Agency of Xinjiang Uygur Autonomous Region，Urumqi 830011，Xinjiang，China）

Abstract： Based on seismic data from 1920 to 2020， b value of the Qaidam Altun seismic belt was calculated and analyzed with multiple schemes （different time periods and magnitude periods）， the b value was determined to be 0.78.It is lower than the value of 0.84 in engineering earthquakes comes from “Seismic ground motion parameters zonation map of China”， which is of great significance for seismic safety evaluation of major construction projects in seismic zones and engineering fortification of seismic zoning.As a supplement to “Seismic ground motion parameters zonation map of China”， b value determined in this paper can provide a scientific basis for the establishment of seismic hazard calculation model of major infra-structure.

Key words： Qaidam Basin-altun seismic belt；b value；Multiple schemes；Statistical calculations