抗震設(shè)計(jì)中多遇地震的問(wèn)題和討論＊

石樹(shù)中 沈建文

1）中國(guó)杭州310013浙江省地震局

2）中國(guó)上海200062上海市地震局

抗震設(shè)計(jì)中多遇地震的問(wèn)題和討論＊

石樹(shù)中1）沈建文2），

1）中國(guó)杭州310013浙江省地震局

2）中國(guó)上海200062上海市地震局

指出作為截面抗震驗(yàn)算標(biāo)準(zhǔn)的“小震”混淆了兩種明顯不同的概念：其一是對(duì)應(yīng)于50年超越概率63%的地震動(dòng)參數(shù)（以下簡(jiǎn)稱(chēng)“多遇小震”）；其二是50年超越概率10%的地震動(dòng)參數(shù)經(jīng)過(guò)折算得到的參數(shù)（以下簡(jiǎn)稱(chēng)“折算小震”）.上述混淆導(dǎo)致抗震設(shè)計(jì)中實(shí)際上執(zhí)行著不同的標(biāo)準(zhǔn).文中討論了中、小地震動(dòng)參數(shù)之間關(guān)系的復(fù)雜性，包括從地震危險(xiǎn)性分析得到的基巖面中、小地震峰值加速度的離散性，一致概率譜的差異，以及中、小地震由于土層反應(yīng)非線(xiàn)性影響導(dǎo)致的場(chǎng)地相關(guān)譜的系統(tǒng)性差異.指出重大工程的抗震設(shè)計(jì)以“多遇小震”取代“折算小震”進(jìn)行截面驗(yàn)算可能導(dǎo)致反應(yīng)譜特征周期Tg明顯低估的問(wèn)題，并提出了相應(yīng)建議.

多遇地震 折算小震 設(shè)計(jì)地震動(dòng)參數(shù) 峰值加速度 反應(yīng)譜特征周期

引言

按照《中華人民共和國(guó)防震減災(zāi)法》，我國(guó)建設(shè)工程抗震設(shè)防要求的確定目前采用兩種做法：重大建設(shè)工程和可能發(fā)生嚴(yán)重次生災(zāi)害的建設(shè)工程，根據(jù)地震安全性評(píng)價(jià)（以下簡(jiǎn)稱(chēng)“安評(píng)”）的結(jié)果，確定抗震設(shè)防要求，進(jìn)行抗震設(shè)防；對(duì)于一般工業(yè)和民用建筑，則按照地震區(qū)劃圖規(guī)定的抗震設(shè)防要求，進(jìn)行抗震設(shè)防.

隨著安評(píng)工作的廣泛開(kāi)展，安評(píng)工作的水平和質(zhì)量對(duì)于我國(guó)的經(jīng)濟(jì)建設(shè)和防震減災(zāi)的影響也越來(lái)越大.為此，加強(qiáng)對(duì)于安評(píng)工作的管理，規(guī)范統(tǒng)一安評(píng)的技術(shù)要求應(yīng)是我國(guó)工程地震界的重要任務(wù).

目前，我國(guó)重大工程可能直接用安評(píng)得到的50年超越概率63%的“多遇小震”地震動(dòng)作為設(shè)計(jì)地震動(dòng)參數(shù).一般工業(yè)與民用建筑按《中國(guó)地震動(dòng)參數(shù)區(qū)劃圖》（GB18306—2001）給定的地震動(dòng)參數(shù)作為抗震設(shè)防要求.該區(qū)劃圖給出的地震動(dòng)參數(shù)相應(yīng)的超越概率為50年10%（以下簡(jiǎn)稱(chēng)“中震”）（國(guó)家質(zhì)量技術(shù)監(jiān)督局，2001）.實(shí)際上，抗震驗(yàn)算時(shí)使用的設(shè)計(jì)地震動(dòng)參數(shù)是經(jīng)過(guò)折算的“折算小震”.兩者的意義有所不同，地震動(dòng)參數(shù)的大小與反應(yīng)譜形態(tài)實(shí)際上可能有很大的差別.

沈建文和石樹(shù)中（2004）通過(guò)危險(xiǎn)性分析的實(shí)際例子等，討論了50年超越概率2%—3%、10%和63%的地震動(dòng)，即“大震”、“中震”和“小震”的相對(duì)大小關(guān)系，指出了上述3者比例關(guān)系的不確定性.并提出建議：“由業(yè)主根據(jù)工程的性質(zhì)確定結(jié)構(gòu)系數(shù)，將中震進(jìn)行折算（在業(yè)主無(wú)法確定結(jié)構(gòu)系數(shù)時(shí)，暫取1／3作為折算系數(shù)），并與小震地震動(dòng)相比取保守值作為截面驗(yàn)算的標(biāo)準(zhǔn).”

從目前安評(píng)工作的實(shí)際情況來(lái)看，仍有許多研究者簡(jiǎn)單地建議以多遇小震作為截面驗(yàn)算的標(biāo)準(zhǔn).

本文在討論多遇小震與折算小震差異的基礎(chǔ)上，指出以多遇小震作為截面驗(yàn)算標(biāo)準(zhǔn)存在的問(wèn)題，提出安評(píng)工作中作為截面驗(yàn)算標(biāo)準(zhǔn)的地震動(dòng)參數(shù)的確定方法.

1 以多遇地震作為設(shè)計(jì)基準(zhǔn)的由來(lái)和現(xiàn)狀

1.1 三水準(zhǔn)二階段設(shè)計(jì)思想

以“三個(gè)水準(zhǔn)“作為抗震設(shè)防目標(biāo)，即“小震不壞，中震可修，大震不倒”，是由《建筑抗震設(shè)計(jì)規(guī)范（GBJ11—89）》（中華人民共和國(guó)城鄉(xiāng)建設(shè)環(huán)境保護(hù)部，1989）正式提出，其含義包括以下兩個(gè)方面：

其一為定義大、中、小震分別與危險(xiǎn)性分析超越概率2%—3%、10%和63%的烈度對(duì)應(yīng)，根據(jù)當(dāng)時(shí)中國(guó)建筑科學(xué)院抗震研究所對(duì)全國(guó)60多個(gè)城市地震危險(xiǎn)性分析結(jié)果的統(tǒng)計(jì)，大、中震烈度之間的差值平均為1度（峰值加速度之比為4∶3—2∶1），中、小震烈度之間的差值大體為1.55度（峰值加速度之比為3∶1）.

其二為三水準(zhǔn)二階段設(shè)計(jì)思想，即要求“小震不壞，中震可修，大震不倒”.三水準(zhǔn)的設(shè)防目標(biāo)通過(guò)“二階段設(shè)計(jì)”來(lái)實(shí)現(xiàn).第一階段設(shè)計(jì)是承載力驗(yàn)算，即采用第一水準(zhǔn)的地震動(dòng)參數(shù)，計(jì)算結(jié)構(gòu)的彈性地震作用標(biāo)準(zhǔn)值和相應(yīng)的地震作用效應(yīng)，進(jìn)行結(jié)構(gòu)構(gòu)件的截面承載力驗(yàn)算，以期滿(mǎn)足在第一水準(zhǔn)下具有必要承載力可靠度，即小震不壞；又滿(mǎn)足第二水準(zhǔn)的損壞可修，即“中震可修”的目標(biāo).對(duì)于有特殊要求的建筑，地震時(shí)易于倒塌的結(jié)構(gòu)等，除進(jìn)行第一階段設(shè)計(jì)外，還需進(jìn)行結(jié)構(gòu)薄弱部位的彈塑性變形驗(yàn)算，要求變形不超過(guò)容許的限值以防倒塌.上述思想在《建筑抗震設(shè)計(jì)規(guī)范（GB50011—2001）》（中華人民共和國(guó)建設(shè)部，2001）中得以延續(xù).

顯而易見(jiàn)，地震危險(xiǎn)性概率分析得到的多遇小震與建筑抗震設(shè)計(jì)規(guī)范中使用的折算小震，其含義與計(jì)算方法有明顯的不同，只是存在統(tǒng)計(jì)上的對(duì)應(yīng)關(guān)系.一般工業(yè)與民用建筑抗震驗(yàn)算使用的小震并不是50年超越概率63%的多遇小震，而是中震（50年超越概率10%的基本烈度或設(shè)計(jì)基本地震加速度）按0.35倍折算得到的折算小震.統(tǒng)計(jì)結(jié)果針對(duì)的是多遇小震，實(shí)際使用的是折算小震.

我國(guó)64、74、78抗震規(guī)范中都引入了“結(jié)構(gòu)影響系數(shù)”的概念，采用對(duì)基本烈度按彈性反應(yīng)譜計(jì)算的地震作用進(jìn)行折減來(lái)確定設(shè)計(jì)地震作用的設(shè)計(jì)方法.結(jié)構(gòu)影響系數(shù)C值因結(jié)構(gòu)不同而異，延性較大、抗倒塌能力較強(qiáng)的結(jié)構(gòu)，結(jié)構(gòu)影響系數(shù)較小.我國(guó)一般認(rèn)為結(jié)構(gòu)的C值在0.25—0.5之間，常粗略認(rèn)為接近1／3.抗震設(shè)計(jì)規(guī)范GBJ11—89與GB50011—2001認(rèn)為，用折算小震取代中震作截面驗(yàn)算與用C值折減效果相當(dāng)，實(shí)際上相當(dāng)于對(duì)不同結(jié)構(gòu)統(tǒng)一粗略地取結(jié)構(gòu)影響系數(shù)0.35.由以上分析可見(jiàn)，對(duì)于不同類(lèi)型的工程結(jié)構(gòu)，結(jié)構(gòu)影響系數(shù)取值不同，對(duì)應(yīng)于多遇小震的超越概率水準(zhǔn)也會(huì)有所不同.

對(duì)于不同結(jié)構(gòu)取統(tǒng)一結(jié)構(gòu)影響系數(shù)的做法，不少學(xué)者持有疑義（翟長(zhǎng)海，謝禮立，2006）.我國(guó)某些行業(yè)抗震設(shè)計(jì)規(guī)范仍以中震為基準(zhǔn)，輔以折減系數(shù)，如水工建筑物抗震設(shè)計(jì)規(guī)范（DL5073—2000）（中華人民共和國(guó)經(jīng)濟(jì)貿(mào)易委員會(huì)，2001）.

1.2 重大工程抗震設(shè)防要求關(guān)于多遇地震的確定

從我國(guó)現(xiàn)行安評(píng)工作看，對(duì)于大部分重大工程，危險(xiǎn)性分析的超越概率通常包括3種超越概率，即50年63%、10%和2%—3%（吳健，高孟潭，2004）.我國(guó)安評(píng)規(guī)范中并未具體規(guī)定概率標(biāo)準(zhǔn).僅《工程場(chǎng)地地震安全性評(píng)價(jià)》宣貫教材13.2.3中提及“對(duì)一般建設(shè)工程的地震小區(qū)劃工作中，設(shè)防地震的概率水平通常取50年超越概率63%、10%和2%—3%”（盧壽德，2006）.既然安評(píng)給出3種超越概率的地震動(dòng)參數(shù)，多遇小震可能理所當(dāng)然地被推薦為截面驗(yàn)算標(biāo)準(zhǔn).

1.3 國(guó)際上關(guān)于抗震設(shè)防的概率標(biāo)準(zhǔn)

各國(guó)區(qū)劃和抗震設(shè)計(jì)并不采用50年超越概率63%多遇小震作為基準(zhǔn)，如美國(guó)UBC規(guī)范以50年超越概率10%的地震動(dòng)作為設(shè)計(jì)基準(zhǔn)，NEHRP2003規(guī)范以50年超越概率2%地震動(dòng)的2／3作為設(shè)計(jì)基準(zhǔn)，歐洲Eurocode 8以50年超越概率10%的地震設(shè)計(jì)為基準(zhǔn).

據(jù)作者掌握的資料，國(guó)外重大工程安評(píng)一般也不采用50年超越概率63%的地震動(dòng)參數(shù).

2 以多遇小震作為設(shè)計(jì)基準(zhǔn)產(chǎn)生的問(wèn)題

3種超越概率水平的地震動(dòng)參數(shù)之間的關(guān)系取決于多方面的因素.沈建文和石樹(shù)中（2004）在“大中小震與設(shè)防標(biāo)準(zhǔn)”一文中僅考慮危險(xiǎn)性分析的影響，且主要是峰值加速度.實(shí)際上，對(duì)于土層場(chǎng)地，三水準(zhǔn)地表地震動(dòng)參數(shù)的差異還會(huì)受到土層地震效應(yīng)的影響；地震動(dòng)參數(shù)也不僅是峰值加速度，還有反應(yīng)譜等，情況要復(fù)雜得多.

2.1 三水準(zhǔn)地震動(dòng)峰值加速度比值的不確定性

1992年，高孟潭等采用編制1990區(qū)劃圖的大量基礎(chǔ)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，研究表明，50年超越概率10%與50年超越概率63%之間，50年超越概率2%—3%與50年超越概率10%之間，烈度差值較為離散，中、小震的烈度差為0.5—2.8度，大、中震為0.3—1.2度，其分布與區(qū)域有關(guān)（高孟潭，韓煒，1992）.

美國(guó)加利福尼亞海岸，50年10%與2%概率水平地震動(dòng)的差別，明顯比美國(guó)東部和中部地震較不活躍地區(qū)的相應(yīng)地震動(dòng)的差別?。∟ordenson，Bell，2000）.

大、中、小地震差值離散的情況在國(guó)外也存在，但是由于我國(guó)采用地震危險(xiǎn)性分析的綜合概率法，地震帶內(nèi)各潛在震源的地震活動(dòng)性參數(shù)可以根據(jù)地震預(yù)報(bào)的經(jīng)驗(yàn)作人為二次調(diào)整，故上述離散可能更顯著.以上海市地震動(dòng)參數(shù)區(qū)劃為例，根據(jù)地震危險(xiǎn)性概率分析計(jì)算結(jié)果，對(duì)比上海市行政區(qū)較小范圍內(nèi)選擇的4個(gè)場(chǎng)地的基巖峰值加速度，小震峰值加速度僅略有差異（0.190—0.213m／s2），中震的差別卻增大了1倍（0.635—1.339m／s2），而大震的差別則已經(jīng)增加了2倍（1.125—3.453m／s2）（上海市地震局，同濟(jì)大學(xué)，2004）.由此可見(jiàn)，大、中、小震之間的關(guān)系并不統(tǒng)一.在這種情況下若以多遇小震作為設(shè)計(jì)基準(zhǔn)，則顯然對(duì)4個(gè)場(chǎng)地滿(mǎn)足“中震可修”和”大震不倒”的抗震設(shè)防要求是不同的.這種人為不當(dāng)?shù)奶幚砜赡茉斐煽拐鹪O(shè)計(jì)標(biāo)準(zhǔn)的重大失誤.

在絕大多數(shù)情況下，我國(guó)安評(píng)得到的中震的地震動(dòng)參數(shù)顯然會(huì)受中國(guó)地震動(dòng)參數(shù)區(qū)劃圖的影響.但在安評(píng)給出的50年超越概率10%的基本烈度或設(shè)計(jì)基本地震加速度與國(guó)家地震區(qū)劃圖完全相同的條件下，即使不考慮土層反應(yīng)的影響，僅由于地震危險(xiǎn)性分析導(dǎo)致的中震和小震的離散性，就可以使不同工程多遇小震的烈度和峰值加速度產(chǎn)生重大差異.

2.2 三水準(zhǔn)地震動(dòng)反應(yīng)譜的不確定性

上面僅討論了烈度和基巖面的峰值加速度，實(shí)際上，多遇小震與折算小震在基巖面危險(xiǎn)性分析的環(huán)節(jié)，對(duì)基巖面反應(yīng)譜也可以造成很大的差異.基巖面一致概率譜取決于潛在震源的分布、參數(shù)以及衰減規(guī)律，超越概率不同，一致概率譜肯定也會(huì)有差異.

吳健和高孟潭（2004）統(tǒng)計(jì)了經(jīng)過(guò)國(guó)家地震安全性評(píng)價(jià)委員會(huì)評(píng)審的56份工程場(chǎng)地地震安全性評(píng)價(jià)報(bào)告和地震小區(qū)劃報(bào)告，共涉及全國(guó)19個(gè)省、市、自治區(qū)的130個(gè)場(chǎng)點(diǎn)，得到了Ⅱ類(lèi)場(chǎng)地3個(gè)不同超越概率水平的統(tǒng)計(jì)結(jié)果，分別對(duì)應(yīng)抗震設(shè)計(jì)規(guī)范中的大、中、小地震，如表1所示.

表1 Ⅱ類(lèi)場(chǎng)地上場(chǎng)地相關(guān)設(shè)計(jì)反應(yīng)譜特征周期Tg統(tǒng)計(jì)特征（單位：s）Table 1 Statistical features of site－related design spectrum characteristic period Tg（unit：s）for classⅡsite

從表1中可以看出，不同的超越概率水平下，特征周期Tg取值有顯著差異.隨著超越概率的減小，特征周期值明顯有增加的趨勢(shì).

事實(shí)上，對(duì)于Ⅱ類(lèi)場(chǎng)地，由于土層的非線(xiàn)性不強(qiáng)，上述特征周期Tg的取值差異還不顯著.對(duì)Ⅳ類(lèi)場(chǎng)地，由于土層較強(qiáng)的非線(xiàn)性影響，3個(gè)超越概率的地表場(chǎng)地相關(guān)譜將有更大的差異，且規(guī)律性更明顯.例如，上海地區(qū)大部分為Ⅳ類(lèi)場(chǎng)地，上海市地震動(dòng)參數(shù)區(qū)劃收集了上海市的安評(píng)資料，包括上海行政區(qū)內(nèi)不同位置，不同超越概率的基巖地震動(dòng)，以及不同土層柱狀圖土層反應(yīng)的組合結(jié)果共328個(gè).地表水平加速度反應(yīng)譜的特征周期Tg與基巖地震動(dòng)水平峰值加速度X的回歸關(guān)系式（1），及散點(diǎn)圖與回歸直線(xiàn)均示于圖1（上海市地震局，同濟(jì)大學(xué)，2004）.

圖1 地表加速度反應(yīng)譜的特征周期Tg與基巖地震動(dòng)峰值加速度X的關(guān)系Fig.1 Relationship between characteristic periods（Tg）of surface acceleration response spectra and bedrock peak acceleration（X）

從圖1可以看出，地表水平加速度反應(yīng)譜的特征周期Tg與基巖水平峰值加速度X有較好的線(xiàn)性關(guān)系.基巖水平峰值加速度越小，特征周期Tg越??；隨著基巖地震動(dòng)加大，特征周期Tg則明顯加大.按照回歸結(jié)果，當(dāng)X取圖中最小值0.252m／s2時(shí)，特征周期Tg為0.53s；而當(dāng)X取圖中最大值2.15m／s2時(shí)，特征周期Tg為1.32s.實(shí)際上，隨著基巖輸入的增大，由于土層非線(xiàn)性的影響和土層對(duì)于輸入地震動(dòng)某頻段的共振等效應(yīng)，在地震波中高頻成分縮減的同時(shí)，某頻段的反應(yīng)譜得到了放大，地表反應(yīng)譜峰值向較長(zhǎng)周期部分移動(dòng)，上述現(xiàn)象也不難理解.

3 討論與結(jié)論

關(guān)于大、中、小地震動(dòng)峰值加速度比值的離散性早已提出，但至今未引起重視.“小震不壞”也可以有兩種不同的理解：一是確實(shí)有多遇地震不壞的要求，如特殊結(jié)構(gòu)或其中的設(shè)備確需保證遭遇多遇地震時(shí)不壞；二是僅試圖通過(guò)“小震不壞”的設(shè)計(jì)達(dá)到“中震可修”和“大震不倒”的目的.對(duì)絕大多數(shù)工程而言，“小震不壞”取第二種理解.小震有多遇小震和折算小震之分.此時(shí)，合理的做法顯然應(yīng)該是：對(duì)一般工程，采用折算小震參數(shù)作為設(shè)計(jì)基準(zhǔn)；對(duì)于需要確?！岸嘤龅卣鸩粔摹钡闹卮蠊こ蹋瑧?yīng)取多遇小震與折算小震的較大地震動(dòng)峰值加速度和較長(zhǎng)的反應(yīng)譜特征周期作為設(shè)計(jì)基準(zhǔn).

高孟潭，韓煒.1992.抗震設(shè)計(jì)中大、中、小地震的確定［M］∥地震工程研究文集.北京：地震出版社：14－20.

國(guó)家質(zhì)量技術(shù)監(jiān)督局.2001.GB18306—2001中國(guó)地震動(dòng)參數(shù)區(qū)劃圖［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社：1－2.

盧壽德主編.2006.GB17741—2005《工程場(chǎng)地地震安全性評(píng)價(jià)》宣貫教材［M］.北京：中國(guó)標(biāo)準(zhǔn)出版社：111－112.

上海市地震局，同濟(jì)大學(xué).2004.上海市地震帶參數(shù)區(qū)劃［M］.北京：地震出版社：188－191.

沈建文，石樹(shù)中.2004.大中小震與抗震設(shè)防標(biāo)準(zhǔn)［J］.地震學(xué)報(bào)，26（5）：533－538.

吳健，高孟潭.2004.場(chǎng)地相關(guān)設(shè)計(jì)反應(yīng)譜特征周期的統(tǒng)計(jì)分析［J］.中國(guó)地震，20（3）：263－268.

翟長(zhǎng)海，謝禮立.2006.抗震規(guī)范應(yīng)用強(qiáng)度折減系數(shù)的現(xiàn)狀及分析［J］.地震工程與工程振動(dòng)，26（2）：1－7.

中華人民共和國(guó)城鄉(xiāng)建設(shè)環(huán)境保護(hù)部.1989.GBJ11—89建筑抗震設(shè)計(jì)規(guī)范［S］.北京：中國(guó)建筑工業(yè)出版社：18－24.

中華人民共和國(guó)建設(shè)部.2001.GB50011—2001建筑抗震設(shè)計(jì)規(guī)范［S］.北京：中國(guó)建筑工業(yè)出版社：10－11.

中華人民共和國(guó)經(jīng)濟(jì)貿(mào)易委員會(huì).2001.DL5073—2000水工建筑物抗震設(shè)計(jì)規(guī)范［S］.北京：中國(guó)電力出版社：274－277.

Nordenson G J P，Bell G R.2000.Seismic design requirements for regions of moderate seismicity［J］.Earthquake Spectra，16（1）：205－225.

石樹(shù)中 浙江省地震局高級(jí)工程師.2004年同濟(jì)大學(xué)防災(zāi)減災(zāi)工程及防護(hù)工程專(zhuān)業(yè)畢業(yè)，獲工學(xué)博士學(xué)位.一直從事地震工程方面的研究工作.現(xiàn)為中國(guó)地震學(xué)會(huì)會(huì)員.

注：孟令媛、李宇彤、鄭建常、蔣海昆、王秀英、宋美琴等的簡(jiǎn)介分別見(jiàn)本刊：Vol.33，No.4；Vol.30，No.4；Vol.29，No.4；Vol.17，No.4；Vol.31，No.6；Vol.33，No.5.

Discussion on the problem of frequent earthquakes in seismic design

Shi Shuzhong1）Shen Jianwen2），
1）Earthquake Administration of Zhejiang Province，Hangzhou310013，China
2）Earthquake Administration of Shanghai Municipality，Shanghai200062，China

“Small earthquakes”used as the standard of section seismic checking in seismic design confuses two significantly different concepts.One is the ground motion parameters corresponding to 63%exceeding probability in 50 years（hereinafter referred to as“frequent earthquakes”），and the other is the parameters

by reducing that corresponding to 10%exceeding probability in 50years（hereinafter referred to as“reduction small earthquakes”）.Above confusion leads to a result that two different standards are actually followed in seismic design.Complexity of the relation of ground motion parameters corresponding to medium and small earthquakes is discussed，such as dispersion of peak acceleration，differences of consistent probability spectra on the bedrock obtained from seismic hazard analysis，and systematic differences of site－related response spectra between medium and small earthquakes by soil non－linear effect.Substitution frequent earthquakes for reduction small earthquakes in section seismic checking for major engineering projects may result in underestima－tion of response spectrum characteristic periodTg.Some suggestions are put forward.

frequent earthquake；reduction small earthquake；design ground motion parameter；peak acceleration；characteristic period of response spectrum

10.3969／j.issn.0253－3782.2012.01.010

P315.9

A

石樹(shù)中，沈建文.2012.抗震設(shè)計(jì)中多遇地震的問(wèn)題和討論.地震學(xué)報(bào)，34（1）：105－110.

Shi Shuzhong，Shen Jianwen.2012.Discussion on the problem of frequent earthquakes in seismic design.Acta Seismologica Sinica，34（1）：105－110.

浙江省重點(diǎn)科研社會(huì)發(fā)展項(xiàng)目（2005C23075）資助.

2011－04－10收到初稿，2011－05－30決定采用修改稿.

e－mail：jianwenshen＠126.com