謝皓宇, 牛松山, 鄭萬山, 唐光武, 仉文崗

(1. 招商局重慶交通科研設計院有限公司 橋梁工程安全與韌性全國重點實驗室,重慶 400067;2. 重慶大學 土木工程學院,重慶 400045)

0 引 言

交通系統(tǒng)是城市居民日常生活與經(jīng)濟活動的重要的生命線。近年來,由于交通技術發(fā)展以及城市區(qū)域的不斷擴張,道路系統(tǒng)與橋梁網(wǎng)絡正變得愈加復雜,而橋梁結構通常被認為是交通系統(tǒng)基礎設施運營維護安全中最脆弱的一環(huán)。20世紀以來全球發(fā)生的諸多地震災害顯示,由于橋梁工程遭到嚴重破壞而切斷了震區(qū)交通生命線,造成救災工作的巨大困難并使次生災害加重,導致巨大的經(jīng)濟損失[1]。因此,評估橋梁網(wǎng)絡的震災防御性能,并制定相應加固優(yōu)先級策略顯得十分關鍵。

我國早期的橋梁抗震設計規(guī)范由于要求偏低或規(guī)定不夠細致等原因,1990年以前修建的橋梁在抗震設計方面可能存在不同程度的問題,現(xiàn)行的CJJ 166—2011《城市橋梁抗震設計規(guī)范》和JTG/T 231-01—20《公路橋梁抗震設計規(guī)范》,同《AASHTO LFRD橋梁抗震設計細則》[2]等國外標準類似,對橋梁抗震設計的要求提高。21世紀以來,學者對既有橋梁網(wǎng)絡地震風險的決策評估有持續(xù)的研究,M.LIU等[3]提出了一個基于概率的綜合數(shù)學模型,用于評估橋梁網(wǎng)絡關鍵節(jié)點連通性和橋梁抗震性能的概率分析;Y.C.YUE等[4]對HAZUS系統(tǒng)[5]進行了二次開發(fā),建立了意大利特倫蒂諾省的橋梁網(wǎng)絡地震風險的健康曲線;M.M.MESSORE等[6]針對既有橋梁網(wǎng)絡地震風險,提出了一個基于全生命周期成本的概率框架,用以評估不同空間分布的在役橋梁的地震風險,研究中地震風險根據(jù)基于成本目標閾值的年度超越概率制定;J.WANG等[7]使用改進的Order-Ⅱ-Dijkstra算法,為大型的國家級公路橋梁網(wǎng)絡建立了一個貝葉斯網(wǎng)絡地震風險評估系統(tǒng);邢云斐[8]制定了一個城市橋梁抗震性能評估方法,針對合肥地區(qū)城市橋梁的地震危險性級別做出了研究;劉振亮等[9]考慮災后救援決策,基于圖論和結構易損性分析提出了一個橋梁網(wǎng)絡的抗震加固優(yōu)先級決策方法;董優(yōu)等[10]基于抗震韌性理論給出了城市橋梁網(wǎng)絡抗震脆弱程度的分析方法。

綜上,既有橋梁網(wǎng)絡的地震風險評估與分級是近年橋梁工程及防災工程領域的一個熱點,但現(xiàn)有的研究中存在缺陷,包括橋梁網(wǎng)絡中社會屬性重要程度在綜合評估體系中的體現(xiàn)不夠明確,過于依賴于經(jīng)驗判斷以及橋梁結構易損性分析改進優(yōu)化等問題。筆者提出一個多尺度的既有橋梁網(wǎng)絡地震風險綜合評估方法,并改進橋梁結構抗震性能分析手段,最后基于CRITIC-熵的客觀賦權法解決不同的評估標準的加權統(tǒng)一的問題。

1 多尺度評估橋梁網(wǎng)絡地震風險

單一橋梁節(jié)點或路段在區(qū)域橋梁網(wǎng)絡中的地震風險評估需要考慮的重要性因素至少包括3個方面:①橋梁在社會屬性層面上的重要程度;②橋梁節(jié)點在橋梁網(wǎng)絡中幾何屬性的重要程度;③橋梁結構自身的抗震性能的強弱程度。橋梁幾何屬性重要度可以根據(jù)基于圖論的網(wǎng)絡模型和路段拓撲特性來表征,橋梁結構的抗震性能則通常使用易損性分析結果來表征,而橋梁社會屬性重要度與前兩者相比更抽象、難以量化。有鑒于此,筆者提出多尺度的橋梁網(wǎng)絡地震風險評估,將定性描述的橋梁社會屬性重要度作為第一尺度,將橋梁拓撲重要度和易損性分析作為第二尺度。

在地震風險評估的第一尺度中,采用定性的方法,提出評估區(qū)域內(nèi)的特殊路段并將路段上的關鍵橋梁定義為特定等級的地震風險以及抗震加固優(yōu)先級。特殊路段定義為:經(jīng)過黨政軍部門、能源化工、醫(yī)院、學校等重要部門及重要基礎設施的主要道路。第一尺度具體的定性劃分原則參考表1,實際工程評估可根據(jù)對象的特定情況做出相應調(diào)整。經(jīng)過第一尺度的地震風險評估,可以得到橋梁網(wǎng)絡中不同社會屬性等級的橋梁分類,對于同等級橋梁則需要進行第二尺度的地震風險分析,進一步獲取目標區(qū)域橋梁網(wǎng)絡的地震風險等級、抗震加固的優(yōu)先順序。

表1 地震風險綜合評估方法第一尺度重要度劃分原則Table 1 Principles for the first scale importance classification of comprehensive seismic hazard assessment method

第二尺度的風險分析包括基于圖論得到的拓撲特性以及基于易損性分析得到的結構抗震性能2個部分。由于是不同的量綱,因此選擇合適的客觀賦權法將兩個因素加權統(tǒng)一。

筆者提出的多尺度分析方法步驟如圖1。

圖1 多尺度評估橋梁網(wǎng)絡地震風險流程Fig. 1 Flow chart of multi-scale assessment of bridge network’s seismic hazard

2 抗震性能分析及拓撲特征選取

2.1 既有橋梁抗震重要性系數(shù)衰減

對于既有橋梁結構抗震性能的評估,需要考慮既有橋梁的抗震設防標準。在橋梁網(wǎng)絡抗震風險評估與分級研究中,因為已服役一定年限的橋梁擁有更短的剩余使用壽命,其地震風險評估所使用的抗震設防標準應低于新建橋梁的標準。根據(jù)風險一致性原則,可以推導出不同剩余使用壽命的橋梁衰減過后的抗震重要性系數(shù)[11]。

根據(jù)重要性系數(shù)C及基本烈度Ie計算出對應的極值Ⅲ型概率分布pⅢ,a,如式(1):

(1)

式中:ω為最大地震烈度;ε為基本烈度值減去1.55。

根據(jù)極值Ⅲ型概率分布計算出設計基準期內(nèi)的超越概率p,如式(2):

(2)

式中:tL為設計基準期。

根據(jù)風險一致性原則,將式(2)計算得到的p用作一定剩余使用壽命t內(nèi)的超越概率計算,最終得到衰減后的抗震重要性系數(shù)Cr,如式(3):

Cr=10(Ib-Ie)·log102

(3)

式中:Ib為對應的地震烈度,Ib=ω-(ω-ε)·[-ln(1-pⅢ,b)]1/k;υb、pⅢ,b分別為在剩余使用壽命t下,超越概率為p所對應的年超越概率及極值Ⅲ型概率分布,νb=1-(1-p)1/t,pⅢ,b=1-(1-νb)50。

2.2 改進的易損性分析輸入地震動選取

通過目標結構有限元模型的大量的非線性動力時程分析,并選取特定的結構破壞狀態(tài)限值,計算結構的地震響應超出限值的條件概率,可以評估地震風險。

多自由度體系非線性結構的動力分析具有極高的不可預測性,包括剛度的本構模型、阻尼模型等結構參數(shù)的不確定性,設計地震動選取的不確定性及不同的微分方程組,計算求解方法帶來的結果不確定性。對于現(xiàn)有的橋梁結構易損性分析輸入地震動的選取,主要存在2個問題:①選取輸入地震波組過程中鮮有考慮多樣性程度,只有輸入能夠激勵目標模型不同動力響應特性的地震波組,才能降低易損性分析過程的不確定性;②需要大量的輸入地震波,極大提升了設計和研究成本。

地震動強度指標可用來直觀判斷地震波可能激勵結構的響應大小,如峰值加速度或結構特征周期對應的譜加速度等,但強度指標很難反映目標結構的非線性動力特性。因此,考慮簡化目標結構,利用具有相似非線性動力特性的簡化模型的損傷指標來體現(xiàn)輸入地震動的特性。R.HONDA等[12]提出一個以三折線本構的10自由度框架體系,模擬了10層的鋼筋混凝土框架結構,并選取了Park-Ang損傷指標評估地震動特性,但該簡化模型動力學參數(shù)的規(guī)定方法并沒有給出,通過10階模態(tài)參數(shù)推導10自由度彈簧體系的剛度取值也十分困難。筆者根據(jù)規(guī)則梁橋順橋向的動力特性,建立一個二折線非彈性本構的二自由度彈簧體系[13],通過求解模態(tài)方程得到非彈性本構的參數(shù),模擬出目標橋梁結構在順橋向的動力特性,該模型示意如圖2。

圖2 二自由度簡化模型示意Fig. 2 Sketch of the simplified two-degree-of-freedom model

假設質(zhì)量1與質(zhì)量2均取m,可以建立關于2個非線性彈簧的初始剛度k1及k2關于順橋向前2階頻率ω1和ω2的模態(tài)方程:

(4)

求解式(4),可得:

(5)

根據(jù)式(5),通過目標橋梁結構的有限元模型模態(tài)分析,得到順橋向前2階振型頻率,即可求解簡化模型的非線性彈簧初始剛度。而基于減隔震梁式橋的動力特性,一般認為減隔震梁式橋順橋向的前兩階振型為主梁振動以及下部結構振動,因此在目標結構為一般減隔震梁式橋的情況下,筆者提出的二自由度簡化模型的參數(shù),可直接根據(jù)橋梁設計參數(shù)計算得到。

對簡化模型進行動力時程分析,將質(zhì)量2的最大響應位移作為損傷指標,可在較短時間內(nèi)計算出作為地震波組中每條地震波所對應的強度指標。針對特定地震波組對應的離散數(shù)據(jù),將極值間等比例劃分區(qū)間,從而得到損傷指標的概率密度函數(shù)(PDF),使用信息熵作為工具以衡量PDF中多樣性的尺度。

作為信息論中,量化度量信源不確定度的常用工具,通過使用信息熵對輸入地震波組的PDF進行評估,選取信息熵最高的候選地震波組作為目標橋梁結構的易損性分析輸入地震波組,使輸入地震波激勵下的目標結構的響應盡量發(fā)散,以概括實際可能發(fā)生的全部地震響應狀況,算法如式(6):

(6)

式中:H(P)為概率密度函數(shù)P對應的信息熵;pn為PDF中第n個損傷指標區(qū)間的概率。

2.3 橋梁網(wǎng)絡拓撲特征選取

城市交通網(wǎng)絡是典型的復雜網(wǎng)絡,忽略網(wǎng)絡中節(jié)點的大小、位置、形狀、功能及節(jié)點間作用關系后,可抽象為由點和線組成的空間幾何圖形,稱為網(wǎng)絡拓撲結構圖。選取橋梁網(wǎng)絡拓撲結構圖中待評估關鍵節(jié)點或邊的度中心性、接近中心性、中介中心性作為評價指標,3項指標均衡量節(jié)點及邊的網(wǎng)絡局部特性,既為獨立指標又互相關聯(lián)?；谶@個特點,利用灰色關聯(lián)度理論[14]描述3項指標間的相對變換和關聯(lián)性,作為評估橋梁的拓撲特性的綜合評價指標。方法如下:

根據(jù)圖論中的定義,求出拓撲網(wǎng)絡中待評價的關鍵節(jié)點或邊的度中心性、接近中心性、中介中心性3項指標。確定節(jié)點的比較數(shù)列以及參考數(shù)列:

Y(j)={Xi(j)}max

(7)

式中:Xi(j)為第i個節(jié)點對應的第j個評價指標;Y(j)為第j個評價指標在所有節(jié)點中的最優(yōu)值。

采用均值法對3項指標進行無量綱化處理,再求出無量綱化后,比較數(shù)列與參考數(shù)列差值絕對值的最大值和最小值:

(8)

式中:Pi(j)為無量綱化后比較數(shù)列與參考數(shù)列差值。

令Pi(j)最大值為Δmax,最小值為Δmin。進一步計算出關聯(lián)系數(shù)ξi(j):

(9)

式中:ρ為分辨系數(shù),筆者研究中取0.5。

最后,將每個目標節(jié)點的3項指標對應的關聯(lián)系數(shù)求和,即得到路網(wǎng)中每座橋梁的關聯(lián)度γi:

(10)

3 基于CRITIC-熵權法的地震風險綜合評價

針對組合測算模型,熵權法及CRITIC法均為目前學術界常用的客觀賦權方法,其中CRITIC法可兼顧考慮各指標間的對比強度和沖突性,并由此綜合確定指標的客觀權數(shù)。然而單一的CRITIC法未能慮及指標間的離散性,熵權法則可有效彌補這一不足[15]。筆者使用CRITIC-熵權法客觀賦權,綜合評價橋梁網(wǎng)絡中關鍵節(jié)點基于拓撲特性以及結構抗震性能的地震風險分級。

對熵權法以及CRITIC法2項指標標準化之后,使用CRITIC法計算客觀權重Rn:

(11)

(12)

(13)

使用熵權法計算客觀權重En:

(14)

(15)

(16)

式中,η′i,n為第i座橋的第n個評價指標;η′n為第n個指標序列;ηi,n為η′i,n對應的臨界標準化值;en為第n個評價評價指標的信息熵。

Rn與En的算數(shù)平均就是橋梁網(wǎng)絡地震風險評估指標的綜合權重。

4 算 例

選取某市主要城區(qū)的道路橋梁網(wǎng)絡作為研究對象,對14座跨江橋梁進行地震風險多尺度評估分析,形成了橋梁道路網(wǎng)絡拓撲特征提取以及綜合賦權法地震風險評估的算例。圖3為選定區(qū)域主要道路的平面網(wǎng)絡模型,空心圓圈代表主要道路上的交叉節(jié)點,實心圓圈代表14座需要地震風險評估和分級的關鍵節(jié)點橋梁。

基于多尺度的風險評估方法,第一尺度中以特殊路段的標準將14座橋梁分為3組,第1組包含2座橋梁(6號、11號),依據(jù)為該路段通過市政府;第2組包含3座橋梁(4號、9號、10號),依據(jù)為相應路段有通過次級的重要基礎設施;剩下的9座橋梁為一組,具有相對最低的優(yōu)先級。

計算出道路網(wǎng)絡中重要橋梁所在的邊的度中心性、接近中心性、中介中心性,并對14座橋梁做抗震易損性分析,易損性分析中輸入地震波的選擇遵循了2.1節(jié)中基于橋梁服役剩余壽命的重要性系數(shù)衰減,以及2.2節(jié)中考慮了相對較高信息熵的輸入地震波組。選取抗震易損性分析中8度對應的地震動峰值加速度(PGA)對應的中等損傷超越概率,根據(jù)CJJ 166—2011《城市橋梁抗震設計規(guī)范》有:

A=C·a

(17)

式中:a為設防烈度的加速度,8度對應取0.3g;C為抗震重要性系數(shù),根據(jù)《城市橋梁抗震設計規(guī)范》中0.3g抗震設防分類對應取值1.7,并根據(jù)每座橋梁的剩余使用壽命對抗震重要性系數(shù)按式(3)進行折減;A為易損性分析所取的時程峰值加速度,按式(17)計算得到為0.51g。

相應累計對數(shù)概率密度分布的易損性曲線如圖4,拓撲特性及損傷概率計算結果如表2。

表2 橋梁中等破壞易損性概率及拓撲特性Table 2 Vulnerability probability and topological characteristics of moderate damage in bridges

將表2中的數(shù)據(jù)代入式(12)～式(14),得到地震易損性和拓撲特性對應的CRITIC權重R1、R2分別為0.380、0.620;將數(shù)據(jù)代入式(15)～式(17),得到易損性和拓撲特性對應的熵權法權重E1、E2分別為0.420、0.580。計算2種權重的算數(shù)平均值,得到綜合權重S1、S2分別取0.400、0.600。將表2數(shù)據(jù)歸一化后再代入綜合權重則能計算出14座橋梁的綜合地震風險,結果見表3。

表3 14座橋梁抗震風險重要度Table 3 Importance of anti-seismic hazard for 14 bridges

將表3中的計算結果結合第一尺度的評估結果,即可得到14座橋梁的抗震風險等級,由高到低排列為:11號>6號>10號>4號>9號>8號>12號>5號>2號>13號>1號>7號>14號>3號?？勺鳛槌鞘袠蛄壕W(wǎng)絡抗震加固優(yōu)先級的決策依據(jù)。

5 結 論

1)針對城市道路網(wǎng)絡中橋梁地震風險評估的問題提出了一個多尺度的評估方法,考慮了既有橋梁抗震設防標準的修正,并改進了傳統(tǒng)的地震易損性分析,對輸入地震波組的選擇提出了以信息熵恒量地震波組多樣性的方法;提出了一個基于灰色關聯(lián)度理論的橋梁網(wǎng)絡拓撲特性的計算方法;綜合考慮了橋梁的結構抗震性能、拓撲特性以及對于城市社會屬性的影響,多尺度地將橋梁的地震風險進行分級,并基于CRITIC-熵權法客觀賦權計算出橋梁的綜合重要度。

2)從算例中可以看到,筆者方法實現(xiàn)了決策過程的量化,可以作為評估道路網(wǎng)絡中橋梁地震風險、判斷城市橋梁抗震加固優(yōu)先級的可靠的依據(jù)。