考慮金屬蓄熱情況下超臨界循環(huán)流化床鍋爐水冷壁流動不穩(wěn)定性數(shù)值分析

謝貝貝，　龍　俊，　茆凱源，　王文毓，　楊　冬,　夏良偉，　殷亞寧

(1． 西安交通大學(xué) 能源與動力工程學(xué)院， 多相流國家重點實驗室， 西安 710049；2． 哈爾濱鍋爐廠有限責(zé)任公司 高效清潔燃煤電站鍋爐國家重點實驗室， 哈爾濱 150046)

?

考慮金屬蓄熱情況下超臨界循環(huán)流化床鍋爐水冷壁流動不穩(wěn)定性數(shù)值分析

謝貝貝1，龍俊1，茆凱源1，王文毓1，楊冬1,夏良偉2，殷亞寧2

(1． 西安交通大學(xué) 能源與動力工程學(xué)院， 多相流國家重點實驗室， 西安 710049；2． 哈爾濱鍋爐廠有限責(zé)任公司 高效清潔燃煤電站鍋爐國家重點實驗室， 哈爾濱 150046)

建立了考慮金屬蓄熱情況下超臨界循環(huán)流化床(CFB)環(huán)形爐膛鍋爐流動不穩(wěn)定性分析的數(shù)學(xué)模型，采用交錯網(wǎng)格對超臨界CFB環(huán)形爐膛鍋爐水冷壁流動不穩(wěn)定性進行模擬，計算分析了對回路施加熱負荷擾動后進、出口質(zhì)量流量及金屬壁溫隨時間的變化.結(jié)果表明：施加1.1倍、1.2倍和1.3倍熱負荷擾動后，進、出口質(zhì)量流量隨時間呈反相脈動，最終進、出口質(zhì)量流量相等且恢復(fù)到穩(wěn)態(tài)值，金屬壁溫隨時間波動后也維持在一常數(shù)值，表明該工況下水冷壁流動仍是穩(wěn)定的.

超臨界循環(huán)流化床鍋爐； 環(huán)形爐膛； 金屬蓄熱； 水冷壁； 流動不穩(wěn)定性； 數(shù)值分析

循環(huán)流化床(CFB)燃燒技術(shù)具有煤種適應(yīng)性強、排放低、燃燒效率高等優(yōu)點.為了實現(xiàn)高效低污染發(fā)電，CFB鍋爐發(fā)電技術(shù)正朝著高參數(shù)和大容量方向發(fā)展.與亞臨界CFB鍋爐相比，超臨界CFB鍋爐的運行參數(shù)高，工質(zhì)流動可能處于單相區(qū)，也可能處于汽水兩相區(qū)，鍋爐水動力不穩(wěn)定性分析的主要目的是防止發(fā)生傳熱惡化.目前，國內(nèi)外學(xué)者對兩相流動不穩(wěn)定性進行了廣泛而深入的研究[1-5]，但對考慮金屬蓄熱的鍋爐水冷壁流動不穩(wěn)定性的研究卻很少.由于金屬蓄熱會對鍋爐爐膛的熱流密度產(chǎn)生影響，進而影響爐膛的動態(tài)特性，因此對考慮金屬蓄熱的鍋爐水冷壁流動不穩(wěn)定性進行研究很有必要.筆者采用交錯網(wǎng)格、集總參數(shù)法建立了考慮金屬蓄熱的超臨界CFB環(huán)形爐膛鍋爐水冷壁流動不穩(wěn)定性分析與壁溫計算的數(shù)學(xué)模型，編寫了以Fortran語言為基礎(chǔ)的水冷壁流動不穩(wěn)定性分析和壁溫計算程序.針對超臨界CFB環(huán)形爐膛鍋爐，采用單通道計算模型，選取水冷壁中最易發(fā)生流動不穩(wěn)定性的回路作為研究對象，模擬了水冷壁中工質(zhì)的流動狀況，為鍋爐的安全運行提供指導(dǎo).

1　數(shù)學(xué)模型

通常采用的流動不穩(wěn)定性數(shù)值分析方法有時域法[6]和頻域法[7-8]2種，壁溫計算可采用集總參數(shù)法或分布參數(shù)法.本文數(shù)值分析方法采用時域法，壁溫計算采用集總參數(shù)法.

1.1基本假設(shè)

基本假設(shè)如下：(1) 工質(zhì)沿管道一維軸向流動；(2) 同一橫截面上溫度場和速度場分布均勻；(3) 只考慮徑向傳熱；(4) 在能量守恒方程中忽略黏性耗散、動能和勢能的影響.

1.2控制方程

質(zhì)量守恒方程為

(1)

動量守恒方程為

高寒沙區(qū)降水量少，風(fēng)沙大，這樣的惡劣環(huán)境下在沙丘(地)造林，成活率和成林率是技術(shù)上的難題。已往的造林常采用裸根苗或機械固沙這兩種措施。裸根苗造林技術(shù)育苗周期長、沙丘上不易挖穴、栽植難度大、受造林季節(jié)限制，不適宜長期在高寒沙區(qū)推廣使用。機械固沙措施采用麥草、粘土、草繩、石礫、尼龍網(wǎng)和PVC編織袋等材料做沙障，短期內(nèi)可以達到固定流動沙丘的目的，然而其有效作用期限短，結(jié)合生物固沙措施才能發(fā)揮最佳治沙效果。

在招標(biāo)文件的編制方面，甲方要根據(jù)項目特點和功能需要進行編制，應(yīng)包括項目建設(shè)的所有實質(zhì)性要求和條件?；ú块T還應(yīng)做好招標(biāo)文件內(nèi)容的審查工作，一是避免出現(xiàn)傾向性、地方保護性條款等違反規(guī)定的；二是招標(biāo)文件中的合同條款的擬定應(yīng)參照住建部門發(fā)布的示范文本，條款應(yīng)具體明確，避免簽訂合同時發(fā)生爭議或糾紛；三是對招標(biāo)圖紙與工程量清單進行嚴(yán)格的控制與審查，盡可能減少遺漏項，避免施工過程中的變更和增項過多；四是對招標(biāo)單位的資質(zhì)進行審核，選擇整體素質(zhì)高、社會效益好、業(yè)內(nèi)信譽好、施工團隊強的施工方，避免工程的再次轉(zhuǎn)包。

Kinδd(z)+Koutδd(z-L)+

(2)

能量守恒方程為

(3)

狀態(tài)方程為

(4)

金屬蓄熱方程為

(5)

工質(zhì)對流側(cè)傳熱方程為

(6)

在全面提升服務(wù)的同時，黑龍江銷售還緊盯油品質(zhì)量不松懈?！坝推烦鰩鞎r有鉛封、編號，保證無其他油品混入。經(jīng)化驗的油品到達油庫后，將進行再次檢測并且留樣。在油品配送過程中，各分公司實行統(tǒng)一配送、全程監(jiān)控、規(guī)定路線、規(guī)定時間的方式，確保秋收農(nóng)業(yè)用油品質(zhì)。”張文興告訴記者。

2　數(shù)值計算方法

2.1離散方法

的確，作為一個專精于打造全地形SUV的品牌，路虎一直做得很出色。但如果僅僅依靠簡單可靠的衛(wèi)士系列和身價有些高不可攀的攬勝家族，這個品牌也絕對發(fā)展不到今天的規(guī)模。在上世紀(jì)80年代末，一個在野外和公路，實用與舒適之間平衡得恰到好處的全新車型的面世，讓路虎品牌有了長足的發(fā)展和進步，也讓更多消費者體驗到了SUV這個在當(dāng)時還算新穎的概念的獨到之處。

圖1　網(wǎng)格劃分

將方程離散化，對主控制體以及速度方向控制體進行研究，可得到以下方程.

質(zhì)量守恒方程為

(7)

能量守恒方程為

(8)

對于沿流動方向的U控制體，則有

動量守恒方程為

檢查患者的循環(huán)系統(tǒng)以及呼吸系統(tǒng)的相關(guān)功能，以保障氧氣的足夠供給，是患者成功進行手術(shù)的必要條件和基礎(chǔ)。如果患者存在心肺功能的相關(guān)問題，會在一定程度上阻礙氧氣的吸入，導(dǎo)致患者出現(xiàn)缺氧。

Kinδd(z)+Koutδd(z-L)+Kjbδd(z-zjb)]

(9)

金屬蓄熱方程為

3.績效考核體系具有客觀公正性，因此可以很好的區(qū)分職工工作狀態(tài)是高效還是低效，根據(jù)不同職工工作水平?jīng)Q定他們不同的工資和福利待遇，合理的分配獎金和晉升機會，這樣才能可以對員工起到激勵作用。

(10)

傳熱方程為

(11)

狀態(tài)方程為

采用控制容積法沿流動方向進行離散化，沿管長方向?qū)⑹軣峁芫鶆騽澐殖蒼個固定網(wǎng)格，利用交錯網(wǎng)格劃分的回路見圖1.

編者按：非洲是全球竹資源集中分布區(qū)之一，竹林面積約占全球竹林面積的11%。非洲竹資源利用及產(chǎn)業(yè)發(fā)展落后，資源優(yōu)勢未能轉(zhuǎn)化為經(jīng)濟優(yōu)勢。為促進非洲竹資源合理開發(fā)，實現(xiàn)綠色可持續(xù)發(fā)展，2018年9月3日在中非合作論壇北京峰會開幕式上，習(xí)近平主席宣布在埃塞俄比亞首都成立“中非竹子中心”，幫助非洲開發(fā)竹藤產(chǎn)業(yè)。鑒于竹藤領(lǐng)域的這一重大事項，為使讀者更多地了解非洲竹資源及竹產(chǎn)業(yè)發(fā)展?fàn)顩r，本刊特邀請國際竹藤組織東非地區(qū)辦事處主任、駐埃塞俄比亞辦公室主任傅金和博士撰寫專稿，本刊也將持續(xù)報道非洲竹產(chǎn)業(yè)的發(fā)展?fàn)顩r，并跟蹤“中非竹子中心”的建設(shè)進展，以饗讀者。

式中：L為管長，m；d為內(nèi)徑，m；A為橫截面積，m2；ρ為流體密度，kg/m3；h為比焓，J/kg；qm為質(zhì)量流量，kg/s；p為壓力，Pa；θ為流動方向與水平面的角度；Kin為進口節(jié)流系數(shù)；Kout為出口節(jié)流系數(shù)；Kjb為局部阻力系數(shù)；f為摩擦阻力系數(shù)；δd為狄拉克函數(shù)；qL、q1、q2為線密度，W/m；α為工質(zhì)對流傳熱系數(shù),W/(m2·K)；Tj為壁溫，℃；Tf為工質(zhì)溫度，℃；cj為金屬比熱容，J/(kg·K)；mj為單位管長金屬質(zhì)量，kg/m；t為時間，s；F2為單位管長的受熱面積，m2；z為管長方向的軸向坐標(biāo),m；zJb為彎頭處管子的軸向坐標(biāo),m.

(12)

在保證計算精度下，傳熱方程可簡化為

(13)

2.2邊界條件和數(shù)值計算步驟

數(shù)值計算的邊界條件如下：(1)整個回路進、出口的壓降差為一定值；(2)進口壓力與進口焓值不隨時間變化.

對回路施加熱負荷擾動后，數(shù)值計算步驟如下：

(1) 求解穩(wěn)態(tài)時工質(zhì)的物性參數(shù),包括壓力、焓值、密度、質(zhì)量流量、溫度以及進、出口壓降.

1.3.2 葡萄糖標(biāo)準(zhǔn)曲線的繪制。精密稱取干燥恒重的標(biāo)準(zhǔn)葡萄糖10 mg，定容于100 mL容量瓶中，配制成濃度為0.1 mg/mL的葡萄糖標(biāo)準(zhǔn)溶液。精密移取0.2、0.4、0.6、0.8、1.0、1.2 mL的葡萄糖標(biāo)準(zhǔn)品溶液，分別置于具塞試管中，依次加水定容至2 mL，另取2.0 mL蒸餾水作空白對照。然后加入1.0 mL濃度5%的苯酚溶液，搖勻，迅速加入5 mL濃硫酸，搖勻后待試管冷卻至室溫即可。在波長490 nm處測定吸光值，以所測吸光值為縱坐標(biāo)，標(biāo)準(zhǔn)葡萄糖濃度為橫坐標(biāo)，繪制標(biāo)準(zhǔn)曲線，并求出線性回歸方程及相關(guān)系數(shù)。

在我們的觀察中有3例不成功，在鼻內(nèi)鏡下吻合口可見瘢痕過度生長，考慮為創(chuàng)面的過度增殖所致，我們總結(jié)在年輕患者及術(shù)后復(fù)查吻合口過度肉芽組織生長的患者中可在吻合口填塞的海綿中使用5-FU以減輕瘢痕形成。

2.3程序驗證

卓別林有一句很幽默的名言：“人們?yōu)槲覛g呼，是因為他們對我了如指掌；人們?yōu)閻垡蛩固箽g呼，是因為沒人弄得懂他?！钡拇_，電影不同于高深的物理理論，卓別林用最直接的方式向不同層次的觀眾們傳達著自己想要表達的觀點，這是人們喜歡他的原因。

The authors are grateful to CINVESTAV del I.P.N. for providing the resources to carry out this work.

隨著城鎮(zhèn)經(jīng)濟社會發(fā)展，京津冀、長三角和珠三角城市群的居民生活用電量逐漸上升，用電量差距越來越大。京津冀、長三角和珠三角城市群的居民生活用電量及變化率見圖1。

采用國外商業(yè)軟件Dynastab計算所得600 MW超臨界W火焰鍋爐的數(shù)據(jù)，來驗證本文所建立的數(shù)學(xué)模型與計算方法在模擬超臨界鍋爐水冷壁流動不穩(wěn)定性時的可靠性.Dynastab軟件由西門子公司開發(fā)，已轉(zhuǎn)讓給國內(nèi)相關(guān)鍋爐制造企業(yè)，但只有可執(zhí)行文件，其應(yīng)用受到很大限制.在相同的幾何結(jié)構(gòu)、熱負荷及熱力參數(shù)條件下，對600 MW超臨界W火焰鍋爐的水冷壁流動不穩(wěn)定性進行了模擬.圖2～圖4給出了本文程序計算結(jié)果與Dynastab軟件計算結(jié)果的對比.

由圖2可知，施加1.2倍熱負荷擾動后，本文程序計算所得進口質(zhì)量流量隨時間的變化曲線與Dynastab軟件計算所得進口質(zhì)量流量隨時間的變化曲線在整體趨勢上一致，說明利用本文程序來模擬超臨界W焰鍋爐水冷壁流動不穩(wěn)定性是可靠的.由圖3和圖4可知，在兩相區(qū)與單相區(qū)本文程序計算出的工質(zhì)溫度與Dynastab軟件計算所得工質(zhì)溫度的相對誤差均為0.04%，而壁溫的相對誤差均為0.2%，均在誤差允許范圍內(nèi)，因此利用本文程序?qū)べ|(zhì)溫度與壁溫進行模擬計算是合理的.

圖2　進口質(zhì)量流量隨時間的變化

圖3　兩相區(qū)工質(zhì)溫度與壁溫隨時間的變化

圖4　單相區(qū)工質(zhì)溫度與壁溫隨時間的變化

3　超臨界CFB鍋爐水冷壁流動不穩(wěn)定性分析

3.1環(huán)形爐膛鍋爐水冷壁結(jié)構(gòu)

受熱面布置困難是超臨界CFB鍋爐在大型化過程中的一個難題.為解決此問題，可采用環(huán)形爐膛即內(nèi)圈與外圈并聯(lián)的布置方案，這種結(jié)構(gòu)不僅可以布置更多的水冷壁來獲得足夠的換熱面積，而且可以降低爐膛高度，降低循環(huán)泵的功耗.所研究超臨界CFB環(huán)形爐膛鍋爐水冷壁外圈前墻6回路的熱負荷Q沿爐膛高度的分布情況與管路結(jié)構(gòu)見圖5.從進口到出口的總長度為55.6 m，光管規(guī)格為直徑32 mm、壁厚6.5 mm.根據(jù)熱負荷分布將水冷壁分為33個管組，每個管組內(nèi)的熱負荷相等，取空間網(wǎng)格長度為0.1 m左右，可將回路劃分為560個管段.

圖5　水冷壁結(jié)構(gòu)示意圖(6回路)

3.2水冷壁動態(tài)特性

爐膛動態(tài)不穩(wěn)定時，水冷壁中工質(zhì)流量波動在某些情況下會導(dǎo)致壁溫頻繁劇烈變化，可能引起管壁金屬材料疲勞失效，應(yīng)避免這種情況出現(xiàn).

在考慮金屬蓄熱情況下對超臨界CFB環(huán)形爐膛鍋爐施加熱負荷擾動，并進行數(shù)值模擬.圖6給出了施加1.2倍熱負荷擾動后進、出口質(zhì)量流量隨時間的波動趨勢.由圖6可以看出，當(dāng)進口質(zhì)量流量增大時，出口質(zhì)量流量減小，當(dāng)進口質(zhì)量流量減小時，出口質(zhì)量流量增大，即進、出口質(zhì)量流量隨時間呈反相脈動；進、出口質(zhì)量流量脈動的振幅隨時間延長逐漸減小直至消失，最終進、出口質(zhì)量流量相等，恢復(fù)到穩(wěn)定狀態(tài)，這表明在1.2倍熱負荷擾動下水動力流動是穩(wěn)定的.圖7給出了不同熱負荷擾動下進口質(zhì)量流量隨時間的波動趨勢.由圖7可以看出，隨著熱負荷擾動倍數(shù)的增加，進口質(zhì)量流量脈動的振幅逐漸增大，但隨時間變化的趨勢是一致的.

圖6　考慮金屬蓄熱時進、出口質(zhì)量流量隨時間的變化

圖7　不同熱負荷擾動下進口質(zhì)量流量隨時間的變化

3.3考慮金屬蓄熱時的環(huán)形爐膛鍋爐水冷壁壁溫分析

在未施加熱負荷擾動時，回路初始穩(wěn)態(tài)進口壓力為10.16 MPa，進口焓值為1 343.097 kJ/kg，進口質(zhì)量流量為0.065 257 kg/s.穩(wěn)態(tài)時各管組的平均工質(zhì)溫度與壁溫以及施加熱負荷擾動后4.1 s(熱負荷最大，進口質(zhì)量流量最小)時和13.1 s(進口質(zhì)量流量最大，熱負荷恢復(fù)到初值)時各管組的平均壁溫見圖8.圖9和圖10分別給出了回路中第275管段(兩相區(qū))和第560管段(出口管段，單相區(qū))的溫度變化.圖8中，各曲線拐點處即為相變點，如穩(wěn)態(tài)時2個拐點分別在第8管組與第26管組處，說明在第1～第7管組中工質(zhì)均為單相液體，在第8～第26管組中工質(zhì)處于兩相區(qū)，此后工質(zhì)為過熱蒸汽；在單相液體區(qū)與兩相區(qū)，壁溫只比工質(zhì)溫度高1～3 K.由圖9和圖10可知，工質(zhì)溫度和壁溫先隨時間波動，最終維持在一常數(shù)值而不再隨時間變化，最高壁溫為434.88 ℃，處于允許范圍內(nèi)，因此該工況下鍋爐運行是安全可靠的.

圖8　各管組溫度分布

圖9　第275管段溫度隨時間的變化

圖10　第560管段溫度隨時間的變化

4　結(jié)　論

(1) 建立了適用于考慮金屬蓄熱的超臨界CFB環(huán)形爐膛鍋爐水冷壁流動不穩(wěn)定性分析的通用數(shù)學(xué)模型，并編寫了數(shù)值計算程序.在相同的參數(shù)條件下，與實際工況的數(shù)據(jù)進行了對比，計算結(jié)果較吻合，表明本文數(shù)學(xué)模型與計算方法用于模擬超臨界CFB環(huán)形爐膛鍋爐水冷壁流動不穩(wěn)定性和壁溫計

算是可靠的.

(2) 在考慮金屬蓄熱情況下對回路施加熱負荷擾動后，最終進、出口質(zhì)量流量相等且恢復(fù)到穩(wěn)態(tài)值，金屬壁溫在允許范圍內(nèi)，且隨時間波動后也維持在一常數(shù)值，表明該工況下水冷壁流動仍是穩(wěn)定的.

文中共有50組有效PPG序列參與回歸模型的建立，每組序列時長為5 s，含有2 000個樣本點，且各組序列對應(yīng)一組真實血壓值。實際上每個脈搏波即可對應(yīng)一組血壓值，考慮到算法的實用性與準(zhǔn)確性，本文針對同一種特征值在一個測試序列中，先求出每個脈搏波的特征值，再取平均值作為最終特征值。假設(shè)一個測試序列有m個脈搏波，對于任意特征值v，本文在第i個脈搏波中取出特征值，i=1,2,...,m，并通過式（8）求出該組序列的特征值v的最終數(shù)值kv。

[1]李志偉，孫獻斌，時正海．600 MW超臨界CFB鍋爐的水動力計算[J]．熱力發(fā)電，2006, 35(12):7-9．

LI Zhiwei，SUN Xianbin，SHI Zhenghai．Hydrodynamic calculation of a 600 MW supercritical CFB boiler [J]．Thermal Power Generation，2006,35(12):7-9．

[2]李燕，李文凱，吳玉新，等．帶隔墻的600 MW 超臨界循環(huán)流化床鍋爐水冷壁水動力特性[J]． 中國電機工程學(xué)報，2008, 28(29): 1-5．

LI Yan，LI Wenkai，WU Yuxin，etal．Hydrodynamic characteristics of a 600 MW supercritical CFB boiler with partition wall[J]．Proceedings of the CSEE，2008, 28(29): 1-5．

[3]PAN Jie, YANG Dong, CHEN Gongming,etal．Thermal-hydraulic analysis of a 600 MW supercritical CFB boiler with low mass flux[J]．Applied Thermal Engineering, 2012, 32: 41-48．

[4]PAN Jie, WU Gang, YANG Dong．Thermal-hydraulic calculation and analysis on water wall system of 600 MW supercritical CFB boiler[J]．Applied Thermal Engineering, 2015, 82: 225-236．

[5]沈艷榮，楊冬，王龍,等．600 MW超臨界W火焰鍋爐流動不穩(wěn)定性數(shù)值分析[J]．西安交通大學(xué)學(xué)報，2015，49(5)：68-72．

SHEN Yanrong, YANG Dong, WANG Long,etal．Numerical analysis of flow instability in 600 MW supercritical W flame boiler[J]．Journal of Xi'an Jiaotong University，2015，49(5)：68-72．

[6]XIONG T，YAN X，HUANG S F，etal．Modeling and analysis of supercritical flow instability in parallel channels[J]．International Journal of Heat and Mass Transfer，2013，57(2)：549-557．

[7]薛愛軍，程旭．超臨界水熱力系統(tǒng)的穩(wěn)定性的簡化模型分析[J]．核動力工程，2009，30(5)：35-39．

XUE Aijun，CHENG Xu．Stability analysis of a simplified model of supercritical water-cooled system[J]．Nuclear Power Engineering，2009，30(5)：35-39．

[8]ORTEGA G T，CLASS A，LAHEY R T，etal．Stability analysis of a uniformly heated channel with supercritical water[J]. Nuclear Engineering and Design，2008，238(8)：1930-1939．

Numerical Analysis on Flow Instability of Supercritical CFB Boiler Water Wall Considering Metallic Heat Storage

XIEBeibei1,LONGJun1,MAOKaiyuan1,WANGWenyu1,YANGDong1,XIALiangwei2,YINYaning2

(1. State Key Laboratory of Multiphase Flow, School of Energy and Power Engineering, Xi'anJiaotong University, Xi'an 710049, China; 2. State Key Laboratory of Utility Boilers for Clean & Efficient Utilization of Coal, Harbin Boiler Co., Ltd., Harbin 150046, China)

A mathematical model was established for analysis on flow instability in annular furnace of a supercritical CFB boiler considering metallic heat storage, with which the flow instability was simulated using staggered grids, so as to calculate and analyze the variation of inlet and outlet mass flow as well as tube wall temperature along with time under thermal load disturbance on the loop. Results show that once the loop was applied by 1.1, 1.2 and 1.3 times of thermal load disturbance, the inlet and outlet mass flow would fluctuate reversely along with time, both of which would finally get equivalent at the value of steady state; besides, the wall temperature would also fluctuate with time and finally get stable at a constant value, indicating that the flow condition is still stable.

supercritical CFB boiler; annular furnace; metallic heat storage; water wall; flow instability; numerical analysis

2015-09-22

“十二五”國家科技支撐計劃資助項目(2015BAA03B01-01)

謝貝貝(1991-)，女，河南鄭州人，碩士研究生，研究方向為鍋爐水動力流動不穩(wěn)定性數(shù)值分析.電話(Tel.)：029-82668393；

E-mail：bbeixie@163.com.

1674-7607(2016)08-0589-05

TK122

A學(xué)科分類號：470.20