牛志廣，姜 巍，陸仁強，張宏偉

(天津大學環(huán)境科學與工程學院，300072天津)

城市配水系統(tǒng)是一個極其復雜的系統(tǒng)，很容易受到蓄意破壞、系統(tǒng)事故、自然災害等多種威脅.根據(jù)以往研究成果可知［1］，城市供水風險的大小由威脅發(fā)生的概率、嚴重程度及系統(tǒng)自身脆弱性3個因素決定.配水系統(tǒng)面臨威脅(本文以地震災害為例)的嚴重程度和發(fā)生概率是人為無法控制的，但是，供水部門可以改變結(jié)構(gòu)參數(shù)使配水系統(tǒng)的脆弱性降低，達到對供水風險的控制.另一方面，配水系統(tǒng)脆弱性是其在特定的威脅發(fā)生時，由于社會、經(jīng)濟、自然等多個因素的作用，所表現(xiàn)出的固有敏感性和人類應災能力，因此，配水系統(tǒng)的脆弱性是針對特定威脅的.以往在配水系統(tǒng)脆弱性評價方面的研究，多借助于數(shù)值分析模型，以機理性研究為主［2－3］，但是，配水系統(tǒng)脆弱性由管徑、管材、接口形式和周邊環(huán)境等多個控制變量決定，配水系統(tǒng)的復雜性使得這些機理性研究不太適合實際管理的需要.為此，根據(jù)地震發(fā)生時配水系統(tǒng)供水功能轉(zhuǎn)變的突變特征，結(jié)合以往的研究基礎，建立了基于突變理論的配水系統(tǒng)脆弱性綜合評價數(shù)學模型.

1 突變理論基礎

在社會或自然現(xiàn)象中，如地震的發(fā)生、水的沸騰等，這種由量變、漸變的過程發(fā)展到突變、質(zhì)變的過程隨處可見.為有效解決上述非連續(xù)的突變過程，法國數(shù)學家Rene Thom基于拓撲學、結(jié)構(gòu)穩(wěn)定性等數(shù)學基礎提出了突變理論.Rene Thom證明，當系統(tǒng)的控制變量ui的個數(shù)不多于4，狀態(tài)變量xi的個數(shù)不超過2時，最多有7種基本突變形式［4］.突變理論在研究某種過程從某一穩(wěn)定態(tài)躍遷到另一穩(wěn)定態(tài)時，首先假設系統(tǒng)的狀態(tài)可以用一組參數(shù)表征，當系統(tǒng)處于某一穩(wěn)定態(tài)時，對應這一穩(wěn)定態(tài)的函數(shù)便取唯一的極值，而當參數(shù)在一定范圍變化，其函數(shù)的極值大于1個時，該系統(tǒng)便處于不穩(wěn)定態(tài).當系統(tǒng)從一種穩(wěn)定態(tài)進入某個不穩(wěn)定態(tài)時，參數(shù)的微小變化便可使系統(tǒng)從不穩(wěn)定態(tài)轉(zhuǎn)入一種新的穩(wěn)定態(tài)，系統(tǒng)的狀態(tài)便在這一瞬間發(fā)生突變.突變理論的主要優(yōu)勢就是通過對臨界點附近非連續(xù)變化現(xiàn)象的分析來實現(xiàn)對突變過程的解釋.因此，突變理論在物理、化學及社會科學等多學科都得到廣泛應用［5-6］.

2 基于突變理論的配水系統(tǒng)脆弱性綜合評價模型

2.1 配水系統(tǒng)脆弱性評價指標體系

根據(jù)上述脆弱性的定義可知，配水系統(tǒng)的脆弱性指數(shù)(V)由災損敏感性(Vs)和應災能力(Vd)兩個評價指標共同決定.災損敏感性由配水系統(tǒng)的自身結(jié)構(gòu)特征決定，是其遭受某一特定威脅后的損失程度，而配水系統(tǒng)應災能力是供水部門為保障配水系統(tǒng)在特定威脅作用下的供水安全而采取的措施，反映了供水部門應對威脅發(fā)生時的主觀能動性.

根據(jù)大量歷史地震資料及劉遂慶等對2008年“5.12”特大地震現(xiàn)場調(diào)研資料分析可知［7－8］，配水系統(tǒng)的抗震能力主要由管徑、管材、連接方式、管齡及所在區(qū)域的土壤環(huán)境5個因素決定，因此，本文從這幾個方面對配水系統(tǒng)面臨地震威脅時的災損敏感性進行評價.而與配水系統(tǒng)相伴的人類應災能力從人力、財力、物力等3方面進行評價，考慮到這幾個因素均與配水系統(tǒng)所在城市的經(jīng)濟地位有直接關聯(lián)，通常配水系統(tǒng)所屬城市的經(jīng)濟實力越雄厚，用于防震減災及災后重建的資金就越有保障，配水系統(tǒng)恢復正常供水功能的時間就越短，建立如下配水系統(tǒng)脆弱性評價指標體系(見表1、2).

表1 配水系統(tǒng)脆弱性評價指標體系1

表2 配水系統(tǒng)脆弱性評價指標體系2

2.2 配水系統(tǒng)控制變量的選擇及尖點突變模擬

根據(jù)突變理論可知，配水系統(tǒng)不連續(xù)構(gòu)造的個數(shù)只由控制變量ui的個數(shù)確定，而與狀態(tài)變量xi無關.本文將配水系統(tǒng)的整體脆弱性指數(shù)V作為狀態(tài)變量，將災損敏感性指數(shù)Vs、應災能力指數(shù)Vd作為配水系統(tǒng)的控制變量.根據(jù)歷史資料分析可知，配水系統(tǒng)脆弱性主要由Vs決定，而Vd的貢獻次之，因此，將Vs作為主控制變量u1，Vd作為次控制變量u2.

根據(jù)突變理論，對于配水系統(tǒng)尖點突變現(xiàn)象，描述該配水系統(tǒng)的勢函數(shù)為

因此，配水系統(tǒng)的供水功能狀態(tài)可以用3維空間中的相點(V，Vs，Vd)表示，為避免Vs、Vd出現(xiàn)負值，可得到如下配水系統(tǒng)供水功能的尖點突變模型［9］:

勢函數(shù):F(V)=(V-V0)4+(Vs-Vs0)×(V-V0)2+(Vd-Vd0)×(V-V0)，

突變流形M:4(V-V0)3+2×(Vs-Vs0)×(V-V0)+(Vd-Vd0)=0，

分叉集 B:8×(Vs-Vs0)3+27×(Vd-Vd0)2=0.式中V0、Vs0、Vd0為配水系統(tǒng)固有的特征參數(shù).

根據(jù)一元三次方程求根公式和盛金判別法可知，上述配水系統(tǒng)的分叉集方程便是配水系統(tǒng)尖點突變模型的穩(wěn)定判據(jù)，方程4(V-V0)3+2×(Vs-Vs0)×(V-V0)+(Vd-Vd0)=0可以有1個或者3個實根，實根的數(shù)目可由判別式Δ=8×(Vs-Vs0)3+27×(Vd-Vd0)2的值來確定:(1)當Δ＜0時，配水系統(tǒng)位于不穩(wěn)定的狀態(tài)區(qū)域，有可能發(fā)生突變;(2)當Δ＞0時，配水系統(tǒng)控制參數(shù)的穩(wěn)定變化使得狀態(tài)參數(shù)V也隨之穩(wěn)定變化，配水系統(tǒng)功能是穩(wěn)定的，不會發(fā)生突變;(3)當Δ=0時，配水系統(tǒng)處于臨界平穩(wěn)狀態(tài).

根據(jù)上述突變模型的分叉集方程還不能對配水系統(tǒng)的脆弱性進行評價，需借助歸一公式對配水系統(tǒng)脆弱性的主次關系進行演算［10］，目的是將控制變量及狀態(tài)變量的取值范圍限制在［0，1］內(nèi)，尖點突變模型的的歸一公式為

其中:VVs，VVd分別表示災損敏感性指數(shù)Vs、應災能力指數(shù)Vd對配水系統(tǒng)整體脆弱性指數(shù)的貢獻度.

2.3 配水系統(tǒng)脆弱性的突變評價準則

在采用突變理論對城市配水系統(tǒng)脆弱性進行綜合評價時，需考慮“非互補”準則和“互補”準則.當主次控制變量ui對于狀態(tài)變量x的貢獻不可相互替代時，按照上述歸一公式(2)確定狀態(tài)變量x的值，要依據(jù)“非互補”的取值原則，即將突變級數(shù)值xui中的最小值作為配水系統(tǒng)的x值.當主次控制變量ui對于狀態(tài)變量x的貢獻可以互相彌補不足時，則依據(jù)“互補”的取值原則，即取各xui的平均值作為配水系統(tǒng)的x值.而根據(jù)上述配水系統(tǒng)突變模型的主、次控制變量之間關系可知，表1、2中的脆弱性評價指標之間遵循“互補”原則.

2.4 配水系統(tǒng)脆弱性分級標準

根據(jù)上述分析可知，配水系統(tǒng)整體脆弱性的大小由災損敏感性指數(shù)Vs和人類應災能力指數(shù)Vd共同決定，由于Vs和Vd對于V值的貢獻具有相反的作用，在計算V值時，首先對Vs和Vd作同向歸一處理，V值便介于［0，1］內(nèi)，V值越大表示配水系統(tǒng)的脆弱性越高，即越容易被破壞.同時，基于以往的研究結(jié)果，根據(jù)V值的大小可確定如表3所示的配水系統(tǒng)脆弱性分級標準.依據(jù)上述計算的脆弱性指數(shù)V值及表3，可對城市配水系統(tǒng)的脆弱性等級進行評價.

表3 配水系統(tǒng)脆弱性分級標準

3 實例分析

3.1 研究對象

以華北某區(qū)配水系統(tǒng)為例，該區(qū)配水管網(wǎng)共56.975 km，最大管徑為 DN1 000，最小管徑為DN50，共包括108個節(jié)點和125根管段，管段分別用1～125表示.根據(jù)該區(qū)供水部門提供的數(shù)據(jù)統(tǒng)計分析可知，該區(qū)配水管道的材質(zhì)種類較單一，主要是灰口鑄鐵管和球墨鑄鐵管.球墨鑄鐵管采用承插式楔型接頭和柔性接口，灰口鑄鐵管采用石棉水泥剛性接頭和承插式接口.該區(qū)配水管網(wǎng)按照不同管材、管齡、管徑的累計管長及所占比例見表4～6.根據(jù)勘察資料可知，該區(qū)配水管網(wǎng)所在地區(qū)的地層場地屬于III類土［11］，為抗震不利地段.

表4 華北某配水管網(wǎng)不同管材管長的統(tǒng)計結(jié)果

表5 華北某配水管網(wǎng)不同管齡管長的統(tǒng)計結(jié)果

表6 華北某配水管網(wǎng)不同管徑管長的統(tǒng)計結(jié)果

3.2 配水系統(tǒng)脆弱性評價結(jié)果

根據(jù)華北某區(qū)供水管理地理信息系統(tǒng)統(tǒng)計資料可得每一根管段的管材、管徑、管齡、管道連接方式，該區(qū)配水系統(tǒng)所處地區(qū)地層場地為III類土，同時該區(qū)為省會級城市.根據(jù)表1和2確定上述災損敏感性評價指標及應災能力評價指標的值，利用尖點突變模型分別計算災損敏感性指數(shù)Vs和人類應災能力指數(shù)Vd后，便可以根據(jù)式(2)得出該區(qū)配水系統(tǒng)的整體脆弱性指數(shù)V為0.1637，根據(jù)表3可知，該區(qū)配水系統(tǒng)的脆弱性等級為B類.

4 結(jié)語

針對配水系統(tǒng)面臨地震威脅時其供水功能轉(zhuǎn)變的突變性特征，提出一種新的城市配水系統(tǒng)脆弱性評價方法.首先，根據(jù)歷史震害資料及配水系統(tǒng)抗震機理分析結(jié)果，提出配水系統(tǒng)脆弱性評價指標體系及分級標準;通過對配水系統(tǒng)狀態(tài)變量及主次控制變量的確定，建立配水系統(tǒng)尖點突變模擬模型，并通過歸一公式的導出及取值原則的確定，提出配水系統(tǒng)整體脆弱性指數(shù)的計算方法.應用結(jié)果表明，提出的城市配水系統(tǒng)脆弱性評價數(shù)學模型具有實際可行性，可作為供水管理部門進行配水系統(tǒng)的安全設計或改造時的決策依據(jù).

［1］陸仁強，牛志廣，張宏偉.城市供水系統(tǒng)風險評價研究進展［J］.給水排水，2010，36(增刊):4－8.

［2］LEE M，MCBEAN E A，GHAZALI M，et al.Fuzzy-logic modeling of risk assessment for a small drinking-water supply system ［J］.Journal of Water Resources Planning and Management，2009 135(6):547 －552.

［3］OGAWA Y，KOIKE T.Structural design of buried pipelines for severe earthquakes［J］.Soil Dynamics and Earthquake Engineering，2001，21(3):199－209.

［4］FORBES G J，HALL F L.The application of catastrophe theory in modelling freeway traffic operations［J］.Transportation Research Part A:General，1990，24(5):335－344.

［5］KNEMEYER A M，ZINN W，EROGLU C.Proactive planning for catastrophic events in supply chains ［J］.Journal of Operations Management，2009，27(2):141-153.

［6］YANG Kun，WANG Tongxu，MA Zhitao.Application of cusp catastrophe theory to reliability analysis of slopes in open－pit mines［J］.Mining Science and Technology，2010，20(1):71－75.

［7］王紹偉，王永，錢正華，等.“5.12”地震災區(qū)供水系統(tǒng)受損及恢復調(diào)研與分析［J］.中國給水排水，2009，25(20):1－5.

［8］熊易華，羅萬申，陳洵.汶川“5.12”地震給排水設施震損情況與啟示［J］.給水排水，2009，35(10):21－24.

［9］康宇虹，徐照宇.股票價格突變模型的建立與分析［J］.哈爾濱理工大學學報，2007，12(1):100-104.

［10］BARUNIK J，VOSVRDA M.Can a stochastic cusp catastrophe model explain stock market crashes［J］.Journal of Economic Dynamics and Control，2009，33(10):1824-1836.

［11］彭艷菊，唐榮余，呂悅軍，等.天津濱海場地土類別特征及其對地震動的影響［J］.地震工程與工程振動，2004，24(2):46－52.