川滇地區(qū)強(qiáng)震誘發(fā)斷層錯動量值的概率危險性分析

馬亞麗娜 史世波 舒恒 盛謙 崔臻

摘要：在中國川滇強(qiáng)震區(qū)興建的大型工程中，隧道常不可避免地橫跨多條區(qū)域性活動斷裂。受活斷裂錯動的影響，隧洞結(jié)構(gòu)面臨著嚴(yán)重的錯斷威脅。在跨活斷層隧洞抗斷設(shè)計(jì)中，斷層錯動量是重要的輸入?yún)?shù)。然而，傳統(tǒng)的確定性方法提供的斷層錯動量很難反映地震發(fā)生并引發(fā)斷層錯動的隨機(jī)性影響，因此常高估了錯動量大小。為此，提出了強(qiáng)震誘發(fā)斷層錯動的超越概率理論錯動量分析方法，引入了目標(biāo)區(qū)域地震活動性參數(shù)及地震引發(fā)斷層錯動的概率，以工程服役期限內(nèi)，可能發(fā)生的位錯量達(dá)到給定值的概率來評價活動斷層發(fā)生錯動的危險性，從而為穿越活斷裂隧洞的抗錯斷設(shè)計(jì)提供一定參考。為了說明該概率分析方法，以滇中引水工程香爐山隧洞為例，對隧洞跨越的龍蟠-喬后斷裂展開錯動量的估算，并進(jìn)行了參數(shù)分析。

關(guān) 鍵 詞：活動斷層; 錯動量值; 地震活動參數(shù); 概率分析方法; 滇中引水工程; 川滇地區(qū)

中圖法分類號： X43 ? 文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2022.01.028

0 引 言

巖石工程選址時遇到活動斷層通常是遵循避讓的原則。但受選線限制，在大量長距離輸水隧洞的建設(shè)中，特別是中國西南高地震烈度地區(qū)的輸水隧洞建設(shè)中，不可避免地面臨輸水隧洞鄰近或穿越活動發(fā)震斷裂的問題。地震作用會對隧洞結(jié)構(gòu)造成由振動引起的破壞和斷層錯動引起的破壞。

實(shí)踐表明，地震發(fā)生時沿斷層跡線的破壞最為嚴(yán)重，隧洞結(jié)構(gòu)的破壞多發(fā)生在圍巖質(zhì)量差和地層條件有較大變化的斷層破碎帶，國內(nèi)外典型跨斷層隧洞震害案例如表1所列[1-6]。這些錯斷破壞歷史案例，凸顯了開展隧洞抗錯斷研究的重要性，而斷層錯動量是結(jié)構(gòu)抗斷設(shè)計(jì)的輸入依據(jù)。但是，目前對于隧洞安全運(yùn)營考慮斷層錯動輸入的相關(guān)研究較少，缺乏合理的抗斷設(shè)計(jì)依據(jù)?？缁顒訑鄬铀矶纯箶嘣O(shè)防量的確定應(yīng)建立在既安全又經(jīng)濟(jì)的基礎(chǔ)上，隧洞工程設(shè)防量定值過高，會導(dǎo)致工程抗斷處理相對保守，費(fèi)用過高;設(shè)防量過低，又會給工程帶來安全隱患。因此，以考慮斷層錯動破壞作用的危險性為基礎(chǔ)，確定隧洞工程設(shè)防量值，是解決抗斷設(shè)計(jì)問題的關(guān)鍵，且需要通過斷層錯動分析獲得。

目前的研究大多基于地震條件確定情況下的斷層破裂長度，通過室內(nèi)試驗(yàn)[7-9]和數(shù)值分析[10-11]來評價斷層錯動的概率及進(jìn)行危險性分析，或是根據(jù)統(tǒng)計(jì)方法[12-14]，即以地震震害資料為參考，通過統(tǒng)計(jì)回歸建立震級-地表破裂長度、震級-地表位錯量的關(guān)系。這些方法對斷層未來一定時段內(nèi)發(fā)生位錯的可能性進(jìn)行分析，評估其最大潛在位錯量，本質(zhì)上屬于確定性方法。按照這種方法，所有工程只按照最大潛在地表位錯設(shè)防，難以反映地震發(fā)生并引發(fā)斷層錯動的隨機(jī)性影響。此外，強(qiáng)震誘發(fā)斷層錯動量值的確定性分析方法，與地面震動概率危險性研究在形式上難以協(xié)調(diào)。以地震發(fā)生的隨機(jī)性與地震活動參數(shù)不確定性為基礎(chǔ)，提出強(qiáng)震誘發(fā)斷層錯動量值的概率危險性評價方法，是實(shí)現(xiàn)跨活動斷層隧洞抗斷設(shè)計(jì)科學(xué)合理化的重要前提。

本研究以斷層錯動危險性分析方法為基礎(chǔ)，考慮目標(biāo)區(qū)域地震活動性參數(shù)和引入地震引發(fā)錯動的概率，提出強(qiáng)震誘發(fā)斷層錯動量值的概率危險性分析方法。該方法將地震活動性與斷層錯動聯(lián)系起來，可以反映地震在發(fā)生和傳播過程中的隨機(jī)性，并以工程安全服役期內(nèi)，斷層錯動引發(fā)的工程位錯量超過工程結(jié)構(gòu)允許值的概率來評價活動斷層潛在錯動量發(fā)生的危險性，采用“錯動量超越概率”這一術(shù)語表達(dá)計(jì)算結(jié)果。最后，依托滇中引水工程展開計(jì)算，并對關(guān)鍵參數(shù)進(jìn)行敏感性分析，為跨活動斷層隧洞的抗斷設(shè)計(jì)提供一定參考。

1 研究方法

1.1 地震引發(fā)錯動概率的確定

地震引發(fā)錯動概率的確定，其實(shí)是計(jì)算在工程安全服役期內(nèi)，斷層錯動引發(fā)工程位移量超過工程允許值的概率。

由于地震年平均發(fā)生率不出現(xiàn)在1 a超越概率與錯動量關(guān)系的計(jì)算中，僅僅對多年超越概率與錯動量關(guān)系產(chǎn)生一定影響，所以本文以50 a超越概率為分析對象。從圖3可以看出，不同地震年平均發(fā)生率條件下，50 a超越概率與錯動量關(guān)系曲線變化趨勢一致;當(dāng)錯動量值在0.01～6.5 m的范圍內(nèi)時，50 a超越概率隨地震年平均發(fā)生率的增大而小幅增大。

3.2 震級密度分布系數(shù)β

β為震級密度分布系數(shù)，是指某潛在震源區(qū)內(nèi)所有m0≤m≤m1地震所遵循的分布中的系數(shù)，β值可由古登堡-里克特公式lgN=a-bM推得，即β=b×ln10，其中b為震級-頻度關(guān)系斜率的負(fù)值[32]。β的大小對震級M的概率密度函數(shù)f（M）、地震引發(fā)錯動的概率PM和震級M引發(fā)錯動時斷層錯動量u的概率密度函數(shù)f（u）都具有一定影響。本文結(jié)合中國大陸[33]、川滇地區(qū)[26]、華北地區(qū)[34]的震級-頻度關(guān)系表達(dá)式，獲得不同研究區(qū)域的β值（見表6），對震級密度分布系數(shù)β進(jìn)行參數(shù)敏感性分析，所得結(jié)果如圖4所示。

以1 a超越概率和50 a超越概率為分析對象。從圖4可以看出：不同震級密度分布系數(shù)條件下，1 a超越概率曲線變化趨勢相同;同一斷層位錯量下，1 a超越概率隨著研究區(qū)域β值的減小而降低;當(dāng)位錯量小于分米量級時，1 a超越概率因β值不同而產(chǎn)生較大差距。另外，在不同β值條件下，50 a超越概率與錯動量關(guān)系曲線變化趨勢一致;當(dāng)錯動量值在厘米至米的量級范圍內(nèi)時，50 a超越概率隨震級密度分布系數(shù)的減小而降低。

3.3 震級M引發(fā)錯動的條件概率密度函數(shù)f（m）

f（m）為震級M引發(fā)錯動的條件概率密度函數(shù)，常根據(jù)歷史資料統(tǒng)計(jì)回歸獲得，并對地震引發(fā)錯動的概率PM的大小產(chǎn)生影響。前文中已提到參考馬東輝等[15]總結(jié)的西部地區(qū)地面斷裂的發(fā)生頻率與震級M的關(guān)系曲線，統(tǒng)計(jì)回歸得到在震級為6.5～8.0級條件下以震級大小為變量的表達(dá)式?，F(xiàn)針對中國大陸及西部地區(qū)統(tǒng)計(jì)資料，獲得地面斷裂發(fā)生頻率與震級M的關(guān)系曲線，統(tǒng)計(jì)回歸得到震級M引發(fā)錯動的條件概率密度函數(shù)（見表7），進(jìn)行參數(shù)敏感性分析，結(jié)果如圖5所示。

以1 a超越概率和50 a超越概率為分析對象。從圖5可以看出：不同條件概率密度函數(shù)條件下，不同研究區(qū)域的超越概率曲線變化趨勢相同;同一斷層位錯量下，中國大陸與西部超越概率相近，大于華北地區(qū)的超越概率。

4 結(jié) 論

本文以斷層錯動危險性分析方法為基礎(chǔ)，考慮目標(biāo)區(qū)域地震活動性參數(shù)和引入地震引發(fā)錯動的概率，提出基于區(qū)域地震活動性參數(shù)的超越概率理論錯動量分析方法。針對依托滇中引水工程展開計(jì)算，并進(jìn)行關(guān)鍵參數(shù)的敏感性分析，得到以下結(jié)論。

（1） 在跨活動斷層隧洞抗錯斷設(shè)計(jì)中，斷層錯動量的大小是重要的輸入?yún)?shù)。斷層錯動一般由強(qiáng)震誘發(fā)產(chǎn)生，但強(qiáng)震不一定伴隨斷層錯動的發(fā)生。強(qiáng)震引發(fā)錯動概率的確定是超越概率理論錯動量分析的前提。

（2） 強(qiáng)震誘發(fā)斷層錯動的超越概率錯動量估算，是一項(xiàng)復(fù)雜的工作，涉及多方面因素。將區(qū)域地震活動性參數(shù)引入計(jì)算過程，可以使地震活動性與斷層錯動聯(lián)系起來，反映目標(biāo)區(qū)域內(nèi)地震在發(fā)生和傳播過程中的不確定性，使結(jié)果更具合理性和針對性。

（3） 針對川滇地區(qū)穿越龍蟠-喬后斷層的香爐山輸水隧洞進(jìn)行超越概率理論錯動量分析，得到當(dāng)錯動量值在厘米至米的量級范圍內(nèi)時，超越概率隨斷層錯動量的增大而降低;50 a超越概率63%對應(yīng)于地震眾值烈度，斷層錯動量約0.1 m;50 a超越概率10%對應(yīng)于基本烈度，斷層錯動量約為1.2 m;50 a超越概率2%～3%對應(yīng)于罕遇地震，錯動量為2.7～3.3 m;根據(jù)隧洞設(shè)計(jì)規(guī)范，跨斷層隧洞抗斷設(shè)計(jì)按照基本烈度水平設(shè)防，對應(yīng)的斷層錯動量約為1.2 m。

（4） 相同斷層錯動量條件下，不同的地震年平均發(fā)生率v、震級密度分布系數(shù)β、震級M引發(fā)錯動的條件概率密度函數(shù)f（m）均會對超越概率大小產(chǎn)生顯著影響。

（5） 在不同的地震年平均發(fā)生率v、震級密度分布系數(shù)β、震級M引發(fā)錯動的條件概率密度函數(shù)f（m）條件下，超越概率與錯動量關(guān)系曲線變化趨勢一致;當(dāng)錯動量值在厘米至米的量級范圍內(nèi)時，50 a超越概率隨地震年平均發(fā)生率v的減小而降低，超越概率隨震級密度分布系數(shù)β的減小而降低，中國大陸與西部超越概率相近，大于華北地區(qū)的超越概率。

（6） 在跨活動斷裂隧洞抗錯斷設(shè)計(jì)中斷層錯動量輸入條件方面，本研究的結(jié)果可以為其提供參考，但仍認(rèn)識到該研究工作還存在很多不足，如P（U影響位置/Eu，e）為給定地震E和隨機(jī)誤差e及錯動量U條件下，影響到某工程場地并導(dǎo)致錯動的概率，其大小與斷層上場地位置相關(guān)，同時也與斷層的邊界條件相關(guān)。本文只針對目標(biāo)工程，采用常用的邊界條件，即地震可發(fā)生在斷層的任何位置，但破裂不能超過斷層的端點(diǎn)。除了此條常見邊界條件外，還有其他類型邊界條件，將在以后工作中進(jìn)一步總結(jié)和討論。

參考文獻(xiàn)：

[1] 徐錫偉，郭婷婷，劉少卓，等.活動斷層避讓相關(guān)問題的討論[J].地震地質(zhì)，2016，38（3）：477-502.

[2] KENNER S J，SEGALL P.Postseismic deformation following the 1906 San Francisco earthquake[J].Journal of Geophysical Research，2000，105（B6）：13195.

[3] HUO H，BOBET A，F(xiàn)ERNANDEZ G，et al.Load Transfer Mechanisms between Underground Structure and Surrounding Ground：Evaluation of the Failure of the Daikai Station[J].Journal of Geotechnical and Geoenvironmental Engineering，2005，131（12）：1522-1533.

[4] WANG W L，WANG T T，SU J J，et al.Assessment of damage in mountain tunnels due to the Taiwan Chi-Chi Earthquake[J].Tunnelling & Underground Space Technology Incorporating Trenchless Technology Research，2001，16（3）：133-150.

[5] HWANG J H，LU C C.Seismic capacity assessment of old Sanyi railway tunnels[J].Tunnelling and Underground Space Technology，2007，22（4）：433-449.

[6] ROY N，SARKAR R.A Review of Seismic Damage of Mountain Tunnels and Probable Failure Mechanisms[J].Geotechnical and Geological Engineering，2017，35（1）：1-28.

[7] COLE D A，LADE P V.Influence Zones in Alluvium Over Dip-Slip Faults[J].Journal of Geotechnical Engineering，1984，110（5）：599-615.

[8] 蘇生瑞，李慧榮.發(fā)震斷裂上覆土層破壞機(jī)制及破裂擴(kuò)展規(guī)律的試驗(yàn)研究[J].地質(zhì)災(zāi)害與環(huán)境保護(hù)，2002，13（1）：48-53.

[9] 郭恩棟，馮啟民，薄景山，等.覆蓋土層場地地震斷裂實(shí)驗(yàn)[J].地震工程與工程振動，2001，21（3）：145-149.

[10] 郭恩棟，邵廣彪，薄景山，等.覆蓋土層場地地震斷裂反應(yīng)分析方法[J].地震工程與工程振動，2002，22（5）：122-126.

[11] SCOTT R F，SCHOUSTRA J J.Nuclear Power Plant Sitting on Deep Alluvium[J].Journal of Geotechnical Engineering Division，1974，4（100）：449-459.

[12] TOEHER D.Earthquake energy and ground breakage[[J].Bulletin of the Seismological Society of America，1958，48（2）：147-153.

[13] 龍鋒，聞學(xué)澤，徐錫偉.華北地區(qū)地震活斷層的震級-破裂長度、破裂面積的經(jīng)驗(yàn)關(guān)系[J].地震地質(zhì)，2006，28（4）：511-535.

[14] WELLS D L，COPPERSMITH K J.New empirical relationships among magnitude，rupture length，rupture width，rupture area，and surface displacement[J].Bulletin of the Seismological Society of America，1994，84（4）：974-1002.

[15] 馬東輝，李剛，孫平善.強(qiáng)震地面斷裂評估的概率分析研究[J].世界地震工程，2005，21（1）：7-11.

[16] 蘇經(jīng)宇，周錫元.斷層錯動的危險性分析方法[J].地震工程與工程振動，1993，13（4）：14-21.

[17] 王凡，沈正康，王敏，等.應(yīng)用統(tǒng)計(jì)物理方法分析川滇地區(qū)地震危險性[J].地震地質(zhì)，2013（1）：101-112.

[18] 孫堯，吳中海，安美建，等.川滇地區(qū)主要活動斷裂的活動特征及其近十年的地震活動性[J].地震工程學(xué)報，2014，36（2）：320-330.

[19] ZHOU H，LIU H L，KANAMORI H.Source Processes of Large Earthquakes along the Xianshuihe Fault in Southwestern China[J].Bulletin of the Seismological Society of America，1983，2（73）：537-551.

[20] WALLACE，R.The geology of earthquakes[J].Eos，Transactions American Geophysical Union，1998，79（9）：115-115.

[21] ZHOU Q，XU X，YU G，et al.Width distribution of surface ruptures associated with the Wenchuan earthquake：Implication for the setback zone of the seismogenic faults in postquake reconstruction[J].BSSA，2010，5B（100）：2666-2668.

[22] 鄧起東，于貴華，葉文華.地震地表破裂參數(shù)與震級關(guān)系的研究[M].北京：地震出版社，1992，247-264.

[23] 吳熙彥，徐錫偉，于貴華，等.國家川滇實(shí)驗(yàn)場地震地表破裂帶分布圖編制[J].地震地質(zhì)，2018，40（1）：27-41.

[24] 顏天佑，崔臻，張勇慧，等.跨活動斷裂隧洞工程賦存區(qū)域地應(yīng)力場分布特征研究[J].巖土力學(xué)，2018，39（增1）：387-395.

[25] 崔臻，侯靖，吳旭敏，等.脆性巖體破裂擴(kuò)展時間效應(yīng)對引水隧洞長期穩(wěn)定性影響研究[J].巖石力學(xué)與工程學(xué)報，2014，33（5）：983-995.

[26] 程佳.川滇地區(qū)地震危險性預(yù)測模型[D].北京：中國地震局地質(zhì)研究所，2017.

[27] 陳培善，白彤霞，李保昆.b值和地震復(fù)發(fā)周期[J].地球物理學(xué)報，2003，46（4）：510-519.

[28] 聞學(xué)澤，陳農(nóng).四川省地震重點(diǎn)監(jiān)視防御區(qū)未來地震潛勢的研究[J].四川地震，1992（1）：1-25.

[29] WEN X Z，JIA J K.Recent sliprates，Earthquake Recurrence intervals and strong seismic hazards on the northwestern segment of the Xianshuihe falut zone，Earthquake Research in China，1988，2（4）：431-451.

[30] 中華人民共和國行業(yè)標(biāo)準(zhǔn)編寫組.水工建筑物抗震設(shè)計(jì)標(biāo)準(zhǔn)：GB 51247-2018[S].北京：中國計(jì)劃出版社，2018.

[31] 王麗萍，張偉，羅文文，等.基于區(qū)域地震活動性參數(shù)的斷層錯動量概率分析方法[J].自然災(zāi)害學(xué)報，2018，27（4）：166-170.

[32] 趙紀(jì)生，周正華.發(fā)震斷層的永久位移概率評估方法[J].巖石力學(xué)與工程學(xué)報，2009，28（增2）：3349-3356.

[33] 聞學(xué)澤，徐錫偉，龍鋒，等.中國大陸東部中-弱活動斷層潛在地震最大震級評估的震級-頻度關(guān)系模型[J].地震地質(zhì)，2007，2（4）：236-253.

[34] 黃瑋瓊，李文香，曹學(xué)鋒.中國大陸地震資料完整性研究之一：以華北地區(qū)為例[J].地震學(xué)報，1994，3（16）：273-280.

（編輯：黃文晉）