摘要:綜合分析了燃?xì)夤芫W(wǎng)安全評(píng)價(jià)中多因素相互作用指標(biāo)體系,基于耦合理論探討了安全影響因素一級(jí)評(píng)價(jià)指標(biāo)間的協(xié)調(diào)程度,建立了燃?xì)夤芫W(wǎng)多因素相互耦合模型,以哈爾濱某區(qū)燃?xì)夤芫W(wǎng)為例,進(jìn)行了多因素相互耦合的定量分析。
關(guān)鍵詞:燃?xì)夤芫W(wǎng);安全評(píng)價(jià);安全影響因素;耦合;耦合度
Mutual Coupling Analysis of Safety Impact Factors of Gas Network
WANG Qing-shu,TAN Yu-fei,QI Hua-yun
Abstract:The index system of multi-factor interaction in safety assessment of gas network is comprehensively analyzed. The coordination degree among first-level assessment indexes for the safety impact factors 1s discussed based on the coupling theory,and a model for mutual coupling among multi-factors of gas network is established. Taking the gas network in a district of Harbin City for example,the quantitative analysis of mutual coupling among multi-factors is carried out.
Key words:gas network;safety assessment;safety impact factor;coupling;coupling degree
1 概述
作為城市生命線之一的燃?xì)夤芫W(wǎng)是復(fù)雜的系統(tǒng)性綜合設(shè)施,其安全影響因素眾多,地質(zhì)災(zāi)害或外部干擾、腐蝕、管材或施工質(zhì)量問題、運(yùn)營(yíng)管理不當(dāng)?shù)纫蛩乜蓪?dǎo)致長(zhǎng)期服役的燃?xì)夤芫W(wǎng)發(fā)生失效事故,燃?xì)夤芫W(wǎng)安全問題日益突出,且燃?xì)饩哂幸兹肌⒁妆?、有毒等特點(diǎn),燃?xì)夤芫W(wǎng)事故導(dǎo)致的經(jīng)濟(jì)損失巨大,因此,對(duì)供氣的安全可靠性要求也越來越高[1]。燃?xì)夤芫W(wǎng)系統(tǒng)各安全影響因素相互耦合、相互影響,進(jìn)行安全性評(píng)價(jià)時(shí),不能僅對(duì)各因素進(jìn)行單獨(dú)分析,還應(yīng)研究各因素的相互耦合作用,對(duì)燃?xì)夤芫W(wǎng)系統(tǒng)整體安全進(jìn)行分析。本文將影響燃?xì)夤芫W(wǎng)安全的因素分為4個(gè)一級(jí)評(píng)價(jià)指標(biāo),建立相應(yīng)的二級(jí)評(píng)價(jià)指標(biāo),對(duì)燃?xì)夤芫W(wǎng)安全影響因素進(jìn)行多因素相互耦合的分析。
2 市政燃?xì)夤芫W(wǎng)安全評(píng)價(jià)指標(biāo)的建立
天然氣行業(yè)標(biāo)準(zhǔn)SY/T 6621—2005《輸氣管道系統(tǒng)完整性管理》將燃?xì)夤芫W(wǎng)的安全影響因素分為9類,包括外腐蝕危險(xiǎn)、內(nèi)腐蝕危險(xiǎn)、應(yīng)力腐蝕開裂危險(xiǎn)、制造危險(xiǎn)(管子和焊縫)、施工危險(xiǎn)(管子環(huán)焊縫、制造焊縫、折皺彎曲、螺紋脫扣、管子或接頭破裂)、設(shè)備危險(xiǎn)(墊圈和O形圈、控制、泄壓、密封、填料失效)、第三方損壞危險(xiǎn)(第三方損壞立即失效、故意破壞、管子的舊損傷)、誤操作危險(xiǎn)、與天氣有關(guān)的危險(xiǎn)和外力危險(xiǎn)(土體位移、暴雨或洪水、寒冷天氣、雷擊)。
燃?xì)夤芫W(wǎng)安全影響因素復(fù)雜,各影響因素間的相關(guān)性較強(qiáng),部分因素對(duì)管網(wǎng)安全影響具有不確定性[2],建立燃?xì)夤芫W(wǎng)安全評(píng)價(jià)指標(biāo)體系尤其是與人有關(guān)的定性指標(biāo)時(shí),應(yīng)該充分考慮指標(biāo)的模糊性[3、4];相互間沒有影響或影響不大的因素劃為同一級(jí)評(píng)價(jià)指標(biāo)體系,利于說明一級(jí)評(píng)價(jià)指標(biāo)間的相關(guān)性。本文在對(duì)燃?xì)夤芫W(wǎng)綜合分析基礎(chǔ)上,構(gòu)建出燃?xì)夤芫W(wǎng)安全評(píng)價(jià)指標(biāo)體系,共分為4個(gè)一級(jí)評(píng)價(jià)指標(biāo),分別為:與人有關(guān)的安全影響因素指標(biāo)(簡(jiǎn)稱人因因素),與管道本體結(jié)構(gòu)有關(guān)的安全影響因素指標(biāo)(簡(jiǎn)稱管體因素),與管網(wǎng)所處的內(nèi)外環(huán)境有關(guān)的安全影響因素指標(biāo)(簡(jiǎn)稱環(huán)境因素),與管網(wǎng)維修管理有關(guān)的安全影響因素指標(biāo)[5](簡(jiǎn)稱管理因素)。每個(gè)一級(jí)評(píng)價(jià)指標(biāo)包含多個(gè)二級(jí)評(píng)價(jià)指標(biāo),具體如下。
① 人因因素指標(biāo)F1
a. 第三方野蠻施工破壞F1,1:主要指挖掘?qū)е滤查g破壞、打樁導(dǎo)致穿孔破壞、施工致使管道移位或劃痕導(dǎo)致滯后失效等。
b. 設(shè)計(jì)人員失誤F1,2:包括設(shè)計(jì)合理性、材料選擇、安全系數(shù)設(shè)計(jì)等因設(shè)計(jì)人員經(jīng)驗(yàn)、資質(zhì)造成的潛在失誤。
c. 管道員工誤操作F1,3:由于經(jīng)驗(yàn)不足、知識(shí)欠缺或工作壓力等造成的管道操作員工自身失誤。
d. 員工素質(zhì)與經(jīng)驗(yàn)F1,4:包括職工的文化水平、獲得技術(shù)職稱的人數(shù)、安全技術(shù)的掌握程度、安全意識(shí)和安全態(tài)度等。
e. 公眾風(fēng)險(xiǎn)教育和認(rèn)識(shí)水平F1,5:主要指對(duì)城市居民和施工人員的教育宣傳或培訓(xùn)、管網(wǎng)安全知識(shí)普及等。
② 管體因素指標(biāo)F2
a. 管道埋深F2,1。
b. 管道內(nèi)腐蝕F2,2:主要指管道內(nèi)壁與輸送燃?xì)馑s質(zhì)間相互作用造成的腐蝕。
c. 管材制造與焊接缺陷F2,3:主要指管材本身破損、小裂紋或折皺彎曲,焊接時(shí)焊縫存在的缺陷或強(qiáng)度不足。
d. 管道壁厚F2,4。
e. 管道施工影響F2,5:主要指管網(wǎng)建設(shè)與維修期間的施工中可能存在的環(huán)焊縫缺陷、接頭破壞、彎頭應(yīng)力集中、外保護(hù)層破壞等。
f. 使用年限F2,6。
g. 運(yùn)行壓力波動(dòng)F2,7。
h. 管網(wǎng)設(shè)備故障F2,8:主要指控制泄壓設(shè)備故障、密封填料失效、0型墊片失效、管網(wǎng)廠站設(shè)備故障等。
③ 環(huán)境因素指標(biāo)B
a. 管道外腐蝕F3,1:主要有土壤腐蝕、大氣腐蝕、管-地電位腐蝕、散雜電流影響等。
b. 管道應(yīng)力腐蝕F3,2。
c. 與天氣有關(guān)因素F3,3:主要指氣溫過低、雷擊、暴雨洪水等。
d. 土體移動(dòng)F3,4。
e. 城市地質(zhì)災(zāi)害F3,5:包括地震、地面沉降、砂土液化等。
④ 管理因素指標(biāo)F4
a. 管道維修管理F4,1。
b. 安全防護(hù)應(yīng)急措施F4,2:包括應(yīng)急預(yù)案的編制、應(yīng)急設(shè)備配置、事故管理等。
c. 管網(wǎng)事故報(bào)警系統(tǒng)F4,3。
d. 人員巡線頻率F4,4。
e. 政策法律保障F4,5:管網(wǎng)安全規(guī)章和責(zé)任制是否健全、安全責(zé)任制是否明確、責(zé)任落實(shí)及監(jiān)督情況等。
3 燃?xì)夤芫W(wǎng)安全影響多因素相互耦合模型
3.1 耦合及耦合度
耦合是指兩個(gè)(或兩個(gè)以上)系統(tǒng)或要素間通過相互作用而彼此影響的現(xiàn)象,系統(tǒng)或要素間存在相互依賴、相互協(xié)調(diào)、相互促進(jìn)的動(dòng)態(tài)關(guān)聯(lián)關(guān)系。耦合度是描述系統(tǒng)或要素間相互影響的程度,系統(tǒng)由無序趨向有序,其關(guān)鍵在于系統(tǒng)內(nèi)部各變量之間的協(xié)同作用,它們決定著系統(tǒng)變化的特征與規(guī)律,耦合度正是反映這種協(xié)同作用的度量[6]。
燃?xì)夤芫W(wǎng)系統(tǒng)安全影響因素4個(gè)一級(jí)評(píng)價(jià)指標(biāo)間既相互制約又相互對(duì)應(yīng)發(fā)展,它們產(chǎn)生的相互影響程度(即耦合度)反映一級(jí)評(píng)價(jià)指標(biāo)間的協(xié)調(diào)程度。人因因素的控制、管理因素的加強(qiáng)可以減少管體因素和環(huán)境因素的影響,同時(shí)又會(huì)增加額外的成本,反過來制約人因因素和管理因素的改善,通過這4個(gè)一級(jí)評(píng)價(jià)指標(biāo)相互作用使燃?xì)夤芫W(wǎng)的安全達(dá)到相對(duì)穩(wěn)定平衡。
3.2 功效函數(shù)
① 確定最優(yōu)指標(biāo)向量
本文對(duì)二級(jí)評(píng)價(jià)指標(biāo)采用專家打分法賦值,設(shè)二級(jí)評(píng)價(jià)指標(biāo)Fi,j的值為ei,j,可得到一級(jí)評(píng)價(jià)指標(biāo)的初值向量如下:
Ei=(ei,1,ei,2,…,ei,j,…,ei,m)T (1)
式中Ei——一級(jí)評(píng)價(jià)指標(biāo)的初值向量
ei,j——二級(jí)評(píng)價(jià)指標(biāo)Fi,j的值
m——一級(jí)評(píng)價(jià)指標(biāo)包含的二級(jí)評(píng)價(jià)指標(biāo)個(gè)數(shù)
設(shè)emax,i,j為二級(jí)評(píng)價(jià)指標(biāo)的最優(yōu)值,則:
Emax,i=(emax,i,1,emax,i,2,…,emax,i,j,…emax,i,m)T (2)
式中Emax,i——一級(jí)評(píng)價(jià)指標(biāo)的最優(yōu)指標(biāo)向量
emax,i,j——二級(jí)評(píng)價(jià)指標(biāo)Fi,j的最優(yōu)值
② 確定評(píng)判向量
以Emax,i為參考向量,Ei為被比較向量,則二級(jí)評(píng)價(jià)指標(biāo)與最優(yōu)指標(biāo)值間的關(guān)聯(lián)度計(jì)算式為:

式中εi,j——二級(jí)評(píng)價(jià)指標(biāo)Fi,j的值與其最優(yōu)指標(biāo)值間的關(guān)聯(lián)度
k——二級(jí)評(píng)價(jià)指標(biāo)序數(shù),k=1,2,…,m
ρ——分辨系數(shù),一般取0.5或小于0.5
min——取最小值函數(shù)
max——取最大值函數(shù)
一級(jí)評(píng)價(jià)指標(biāo)的評(píng)判向量為:
Pi=(εi,1,εi,2,…,εi,j,…,εi,m)T (4)
式中Pi——一級(jí)評(píng)價(jià)指標(biāo)E的評(píng)判向量
③ 確定二級(jí)評(píng)價(jià)指標(biāo)權(quán)重
二級(jí)評(píng)價(jià)指標(biāo)權(quán)重由相關(guān)專家進(jìn)行評(píng)估,對(duì)兩個(gè)二級(jí)評(píng)價(jià)指標(biāo)之間的相對(duì)重要性進(jìn)行比較,采用層次分析法計(jì)算,求得各二級(jí)評(píng)價(jià)指標(biāo)的權(quán)重值,其值的大小即為該指標(biāo)被接受的優(yōu)先級(jí)。該工作共分為2個(gè)步驟:
第1步,建立成對(duì)比較矩陣。對(duì)二級(jí)評(píng)價(jià)指標(biāo)以一級(jí)評(píng)價(jià)指標(biāo)為準(zhǔn)則進(jìn)行成對(duì)比較,采用Satty 1~9標(biāo)度法(見表1)對(duì)各要素進(jìn)行賦值[7]。
表1 Satty 1~9標(biāo)度法
標(biāo)度
|
含義
|
1
|
兩者重要性相同
|
3
|
一個(gè)因素比另一個(gè)因素稍重要
|
5
|
一個(gè)因素比另一個(gè)因素較重要
|
7
|
一個(gè)因素比另一個(gè)因素很重要
|
9
|
一個(gè)因素比另一個(gè)因素極端重要
|
2、4、6、8
|
介于相鄰重要程度之間
|
成對(duì)比較矩陣Ai為正倒值矩陣,矩陣中各要素不僅為正數(shù),且具有倒數(shù)特征。成對(duì)比較矩陣如下:

式中Ai——一級(jí)評(píng)價(jià)指標(biāo)E的成對(duì)比較矩陣
ai,j——第i個(gè)與第j個(gè)二級(jí)評(píng)價(jià)指標(biāo)Satty 1~9標(biāo)度值
wi——第i個(gè)二級(jí)評(píng)價(jià)指標(biāo)準(zhǔn)則權(quán)重
wj——第j個(gè)二級(jí)評(píng)價(jià)指標(biāo)準(zhǔn)則權(quán)重
第2步,計(jì)算準(zhǔn)則權(quán)重向量并進(jìn)行權(quán)重一致性檢驗(yàn)。在計(jì)算矩陣Ai的一致性時(shí),以矩陣Ai的最大特征值來取代,則求解準(zhǔn)則權(quán)重向量公式如下:
(Ai-λmax,i)Wi=0 (9)
Wi=(w1,w2,…,wi,…,wm)T (10)
式中λmax,i——成對(duì)比較矩陣Ai的最大特征值
Wi——一級(jí)評(píng)價(jià)指標(biāo)Fi的準(zhǔn)則權(quán)重向量
專家所作判斷的合理性程度需要進(jìn)行一致性檢驗(yàn),以一致性指標(biāo)(Consistency Index,CI)與一致性比率(Consistency Ratio,CR)來檢驗(yàn)一致性。當(dāng)兩者均小于0.1時(shí),可認(rèn)為專家所作判斷合理。
一級(jí)評(píng)價(jià)指標(biāo)的綜合指數(shù)計(jì)算式如下:
xi=Wi·Ai (11)
式中xi——一級(jí)評(píng)價(jià)指標(biāo)E的綜合指數(shù)
④ 構(gòu)造一級(jí)評(píng)價(jià)指標(biāo)功效函數(shù)
一級(jí)評(píng)價(jià)指標(biāo)Fi對(duì)系統(tǒng)有序度的功效計(jì)算式如下:

式中ui——一級(jí)評(píng)價(jià)指標(biāo)Fi對(duì)系統(tǒng)有序度的功效值
αi、βi——系統(tǒng)穩(wěn)定時(shí)一級(jí)評(píng)價(jià)指標(biāo)Fi綜合指數(shù)的上、下限值
在系統(tǒng)中ui反映了各指標(biāo)達(dá)到目標(biāo)的滿意程度,ui=0時(shí)為最不滿意,ui=1時(shí)為最滿意,ui的范圍為[0,1]。
3.3 構(gòu)建耦合度函數(shù)
根據(jù)物理學(xué)中的容量耦合概念及容量耦合系數(shù)模型確定兩個(gè)系統(tǒng)的耦合度函數(shù)為:

式中uj——一級(jí)評(píng)價(jià)指標(biāo)Fj對(duì)系統(tǒng)有序度的功效值
Ci,j——第i個(gè)與第j個(gè)一級(jí)評(píng)價(jià)指標(biāo)的耦合度
推廣得到多個(gè)系統(tǒng)整體的耦合度模型為:

式中C——系統(tǒng)整體的耦合度
系統(tǒng)耦合度值的范圍為[0,1)。當(dāng)0≤C<0.3時(shí),市政燃?xì)夤芫W(wǎng)4個(gè)一級(jí)評(píng)價(jià)指標(biāo)間處于較低的耦合階段,此時(shí)的燃?xì)夤芫W(wǎng)規(guī)模較小,人因因素、管體因素、環(huán)境因素、管理因素都不突出,彼此間的影響較?。划?dāng)0.3≤C<0.7時(shí),4個(gè)一級(jí)評(píng)價(jià)指標(biāo)間處于磨合期,管體因素、環(huán)境因素導(dǎo)致燃?xì)夤芫W(wǎng)的安全性降低,需要加強(qiáng)管理,投入人力物力來提高安全性,但這種投入是有限度的,需要進(jìn)行綜合考慮,使得燃?xì)夤芫W(wǎng)的安全性得到提高;當(dāng)0.7≤C<1時(shí),4個(gè)一級(jí)評(píng)價(jià)指標(biāo)間處于高度耦合,彼此相互促進(jìn),良性發(fā)展,有效的管理、適當(dāng)?shù)娜肆ξ锪ν度胧沟霉荏w因素、環(huán)境因素對(duì)燃?xì)夤芫W(wǎng)安全性造成的不良影響得以降低,達(dá)到可接受的程度,保證燃?xì)夤芫W(wǎng)安全性處于低風(fēng)險(xiǎn)區(qū)。
4 實(shí)例研究
本文以哈爾濱某區(qū)燃?xì)夤芫W(wǎng)為例,首先采用調(diào)查表的形式,由11名該領(lǐng)域?qū)<覍?duì)二級(jí)評(píng)價(jià)指標(biāo)打分(每個(gè)二級(jí)評(píng)價(jià)指標(biāo)滿分為10分),取其平均值作為各二級(jí)評(píng)價(jià)指標(biāo)的初值,得到初值向量Ei,分別為:
E1=(9,6,7,7,8)T
E2=(7,8,5,4,6,6,4,8)T
E3=(9,6,6,5,8)T
E4=(9,8,9,7,7)T
參考向量中的元素值全為10,得到最優(yōu)指標(biāo)值向量Emax,i;以該區(qū)燃?xì)夤居涗浀慕鼛啄旯芫W(wǎng)歷史事故數(shù)據(jù)作參考,依據(jù)式(3)、(4),綜合計(jì)算出燃?xì)夤芫W(wǎng)安全影響因素各一級(jí)評(píng)價(jià)指標(biāo)的評(píng)判向量為:
P1=(1,0.5,0.6,0.6,0.75)T
P2=(0.83,1,0.63,0.56,0.71,0.71,0.56,1)T
P3=(1,0.55,0.55,0.47,0.78)T
P4=(1,0.71,1,0.56,0.56)T
由專家對(duì)二級(jí)評(píng)價(jià)指標(biāo)相對(duì)重要性進(jìn)行比較,取其平均值構(gòu)成成對(duì)比較矩陣:


通過mat lab計(jì)算得到各成對(duì)比較矩陣的最大特征值分別為:
λmax,1=5.06
λmax,2=8.08
λmax,3=5.03
λmax,4=5.01
二級(jí)評(píng)價(jià)指標(biāo)構(gòu)成的權(quán)重向量由式(9)計(jì)算,分別為:
W1=(0.44,0.06,0.11,0.11,0.28)T
W2=(0.07,0.27,0.13,0.07,0.20,0.03,0.03,0.20)T
W3=(0.33,0.17,0.17,0.08,0.25)T
W4=(0.28,0.17,0.17,0.28,0.10)T
一致性檢驗(yàn)指標(biāo)均滿足要求,一級(jí)評(píng)價(jià)指標(biāo)的綜合指數(shù)由式(11)計(jì)算,分別為:
x1=0.812
x2=0.829
x3=0.750
x4=0.784
各專家根據(jù)歷史經(jīng)驗(yàn)確定一級(jí)評(píng)價(jià)指標(biāo)變量的上、下限值,由式(12)計(jì)算得到各一級(jí)評(píng)價(jià)指標(biāo)對(duì)管網(wǎng)安全的功效分別為:
u1=0.838
u2=0.516
u3=0.259
u4=0.782
最后,根據(jù)相互耦合度模型,由式(13)、(14)計(jì)算一級(jí)評(píng)價(jià)指標(biāo)間及管網(wǎng)系統(tǒng)整體的耦合度,計(jì)算結(jié)果見表2。
表2 哈爾濱市某區(qū)燃?xì)夤芫W(wǎng)一級(jí)評(píng)價(jià)指標(biāo)間及管網(wǎng)系統(tǒng)整體耦合度
C1,2
|
C1,3
|
C1,4
|
C2,3
|
C2,4
|
C3,4
|
C
|
0.486
|
0.425
|
0.500
|
0.472
|
0.489
|
0.432
|
0.343
|
由表2可知,燃?xì)夤芫W(wǎng)安全影響因素一級(jí)評(píng)價(jià)指標(biāo)間均存在著耦合作用,尤其是人因因素和管理因素這兩個(gè)一級(jí)評(píng)價(jià)指標(biāo)間的耦合作用最強(qiáng),人的行為會(huì)對(duì)管網(wǎng)的安全造成影響,設(shè)計(jì)施工時(shí)都有可能出現(xiàn)失誤,加強(qiáng)管理可以導(dǎo)致這種失誤減少,但是同時(shí)也增加成本,燃?xì)夤緫?yīng)該進(jìn)行綜合考慮。因此,在對(duì)燃?xì)夤芫W(wǎng)進(jìn)行安全性評(píng)價(jià)時(shí),特別是在構(gòu)建評(píng)價(jià)指標(biāo)體系及設(shè)計(jì)算法時(shí),應(yīng)該從系統(tǒng)的整體功能出發(fā),重點(diǎn)分析一級(jí)評(píng)價(jià)指標(biāo)間的相互作用、相互影響,在此基礎(chǔ)上對(duì)燃?xì)夤芫W(wǎng)系統(tǒng)進(jìn)行安全性評(píng)價(jià)。燃?xì)夤芫W(wǎng)系統(tǒng)整體的耦合度為0.343,這表明燃?xì)夤芫W(wǎng)4個(gè)一級(jí)評(píng)價(jià)指標(biāo)間還處于磨合期,不太協(xié)調(diào),需要進(jìn)一步改進(jìn)。因此,應(yīng)該加強(qiáng)燃?xì)夤芫W(wǎng)的安全管理,努力改善燃?xì)夤芫W(wǎng)所處的內(nèi)外環(huán)境,如加強(qiáng)對(duì)職工的安全教育、提升安全設(shè)施的管理狀況、加強(qiáng)安全檢查和監(jiān)察等。
參考文獻(xiàn):
[1] 陳秋雄,周衛(wèi),楊印臣.城市燃?xì)夤艿腊踩u(píng)估中的腐蝕評(píng)價(jià)[J].煤氣與熱力,2004,24(8):9-12.
[2] 王志華,周偉國(guó).城市燃?xì)夤艿朗┕ぐ踩芾淼哪:C合評(píng)價(jià)[J].煤氣與熱力,2006,26(8):21-23.
[3] 付士根,許開立.城市燃?xì)廨斉湎到y(tǒng)安全性模糊綜合評(píng)價(jià)[J].煤氣與熱力,2006,26(11):20-22.
[4] 劉楊,段常貴,聶廷哲,等.城市埋地燃?xì)夤芫W(wǎng)安全性模糊綜合評(píng)價(jià)[J].煤氣與熱力,2007,27(5):12-14.
[5] 何淑靜,周偉國(guó),嚴(yán)銘卿,等.城市燃?xì)廨斉湎到y(tǒng)事故統(tǒng)計(jì)分析與對(duì)策[J].煤氣與熱力,2003,23(12):49-51.
[6] 宋紅麗,薛惠鋒,張哲,等.經(jīng)濟(jì)-環(huán)境系統(tǒng)影響因子耦合度分析[J].河北工業(yè)大學(xué)學(xué)報(bào),2008,37(3):84-88.
[7] VANEGAS L V,LABIB A W. Application of new fuzzy-weighted average(NFWA)method to engineering design evaluation[J].International Journal of Production Research,2001,39(6):1147-1162.
(本文作者:王清樹 譚羽非 齊華云 哈爾濱工業(yè)大學(xué) 黑龍江哈爾濱 150090)
贊 賞 分享
您可以選擇一種方式贊助本站
支付寶轉(zhuǎn)賬贊助
微信轉(zhuǎn)賬贊助

- 注解:本資料由會(huì)員及群友提供僅供閱讀交流學(xué)習(xí),不得用于商業(yè)用途!