碎裂礦段開采與空區處理協同研究

陳慶發，周科平，胡建華，張世超

(1. 廣西大學 資源與冶金學院，廣西 南寧，530004；2. 中南大學 資源與安全工程學院，湖南 長沙，410083)

對于松軟破碎復雜難采的高品位礦體，常規開采方法主要有分層充填法和分段空場嗣后充填法[1-3]。前者存在充填接頂質量差、后期礦柱回收安全性差、生產效率低等缺點；后者存在回采過程中容易出現垮方冒落等安全事故、難以適應礦體厚度多變的不足。針對這 2類方法的缺陷，周旭等[4]通過優化采場結構參數，采取單體水壓支柱護頂和光面爆破控頂方式對水平分層充填采礦方法進行改進；李興尚等[5]通過調整充填材料及配合比，提高了膠結充填體強度的方式，改進了水平分層充填采礦方法；吳賢振等[6]通過調整礦塊結構參數的方式，優化了分段鑿巖階段礦房嗣后充填采礦法。這些改進和優化措施并沒有從根本上實現采礦方法的變革。與此同時，賦存在其中的復雜空區也對這些隱患資源的安全高效開采提出了更高的要求。袁志安[7]提議采用階段崩落法對采空區群下的部分富礦體進行連續采礦的方案；任鳳玉等[8]提出了用斜面蹬碴落礦采礦方法回收活動空區邊緣的部分冒落殘礦；周科平等[9-11]針對大范圍隱患區下低品位厚大礦體開采條件，提出了頂板誘導崩落連續采礦新技術，并將其應用于工程實踐，獲得了成功。但是，對于賦存有復雜采空區群的高價值碎裂礦段的安全高效開采，國內外尚無現成的成套技術和經驗可供借鑒。

1 碎裂礦段開采與空區處理協同的提出

自20世紀60年代Ansoff提出協同理念以來，該理念便一直成為企業界、工程界研究很多問題的指導原則[12-14]。所謂協同指的是事物與事物之間的一種關系，一種相互之間的和諧與正向配合的關系。系統協同指的是通過某種方法來組織和調控所研究的系統，尋求解決矛盾或沖突的方案，使系統從無序轉換到有序，達到協同或和諧的狀態。系統協同的目的就是減少系統的負效應，提高系統的整體輸出功能和整體效應。

碎裂礦段開采與空區處理協同是指二者在礦山整體回采過程中的合作、協調和同步。過去人們對二者的協同關系研究并不多，基本認為二者是一種不協同或低協同的相互關系。事實上，從系統的觀點來看，二者并非完全孤立，包含著相互影響、相互制約、相互促進的協同關系，具有較強的耦合性，在礦山開采期間存在著較高程度的協同效應和協同空間。如果礦山設計的采礦方法能夠適應各種規模和形狀的采空區處理的情況，也就是說，存在的采空區能夠被納入整個礦山的開采布局之中，作為開采系統中的部分切割工程、自由爆破空間或采場，那么，礦山將取得較好的協同效果，輸出較高的協同效應。

2 碎裂礦段及其空區賦存概況

高峰礦是廣西華錫集團控股的大型礦山基地之一，其第100號礦體埋藏在+690～-79 m之間，目前已基本開采完畢。第105號礦體在-79 m標高以下，為第100號礦體的延伸礦體，該礦體在-114～-145 m之間為碎裂礦段，碎塊未膠結，易脫落。第105號礦體碎裂礦段受民采干擾，形成許多不規則的空區，而原民采對該礦段的開采主要為淺采，空區體積不均勻，有少量重疊貫通，大多不連續。采空區結構尺寸統計結果如圖1所示。經現場地質調查[15]，測得碎裂礦段主要有4組節理，其分布形式如圖2所示。傾角為25?的連續節理1組，節理間距為0.5 m；傾角為65?的節理有3組，屬于斷續節理，平均間距約為1.5 m。這2種傾角的節理是碎裂礦段介質不連續的主要表現形式。

圖1 采空區結構尺寸統計Fig.1 Statistical results of goafs structure sizes

圖2 碎裂礦段節理分布形式Fig.2 Joint distribution form in cataclastic ore section

3 碎裂礦段采礦方法設計

為實現復雜條件下軟破礦體安全高效回采，古德生院士繼提出“無間柱連續采礦[16-18]”之后，創造性地提出了“開采環境再造[19-21]”這一新的采礦科學命題，其基本內涵是“突破傳統的采礦方法設計思想的限制，應用新的理論、方法和技術，營造一個良好的礦巖開采環境，最終實現礦石資源的高效回采”。

對第105號礦體碎裂礦段的開采，若仍沿用礦山原采第100號礦體所采用的上向水平分層充填法和分段空場事后充填法[3]，則對于占礦段中大部分的走向長度短、厚度為中等的礦體，采用上向水平分層充填法，采場生產能力將嚴重制約著礦山的生產規模；而采用分段空場事后充填法因要考慮相鄰采場的回采，必須留礦柱和在采場底部構筑人工進路，采場全高度方向上需進行高強度充填，從而使充填成本較高。

針對碎裂礦段的特點以及原采礦方法生產的難點，本著提高生產效率、降低生產成本、保證作業安全、確保較低的資源貧化損失率、與空區處理協同等指導思想，基于采礦環境再造和連續開采理論，綜合提出了一種新的采礦方法，即采礦環境再造無間柱分段分條連續采礦法。該方法按礦體厚度不同(≤15 m，15～30 m，≥30 m)具有不同的表現形式。礦體厚度≤15 m時的采礦工藝如圖3所示。

當階段高度30～50 m時，將階段劃分為2～3個分段：沿礦體長度方向布置礦房和間柱，礦房長度為40～60 m，間柱長度為12～15 m，不留頂柱、底柱。整體上采用兩步驟回采，先回采間柱，進行較高配比的水泥砂漿膠結充填；然后，回采礦房，進行低配比的水泥砂漿充填和小部分高配比膠結充填。間柱和礦房的回采工藝大致相同，都是上、下分段同時爆破，在下分段出礦，全部回采完畢后再進行充填，不同的是間柱的拉槽布置在一端，礦房的拉槽布置在中央，間柱是一端往另外一端后退式回采，礦房是中間往兩端回采。出礦均經下盤的出礦進路從下分段的鑿巖道或平底采用鏟運機出礦(采空區比較大時采用遙控鏟運機出礦)。間柱和礦房的充填均通過上階段的運輸平巷，經出礦進路、充填通風聯絡道下放充填管路進行。間柱充填采用砂漿配比為1∶4的水泥。礦房底部4～6 m采用砂漿配比為 1∶4的水泥砂漿進行充填，其余部分采用砂漿配比為1∶10的水泥進行充填。鑿巖時，新鮮風流自下階段運輸平巷經出礦進路、鑿巖道、切割通風天井、充填通風聯絡道、上階段出礦進路后，排至上階段運輸巷道。回采時，新鮮風流自下階段運輸平巷經出礦進路、鑿巖道、空區、充填通風聯絡道、上階段出礦進路后，排至上階段運輸巷道。

圖3 采礦環境再造無間柱分段分條連續采礦法示意圖(礦體厚度≤15 m)Fig.3 Sketch map of mining environment reconstructing non-pillar sublevel-strip continuous mining method(orebody thickness is less than 15 m)

與上向水平分層充填法和分段空場事后充填法相比，新采礦法具有以下優點：(1)繼承了上向水平分層充填法和分段空場充填法的優點，不留頂、底柱，采場生產能力大，資源貧化損失低，能夠適應不同厚度的傾斜、急傾斜高品位礦體的開采；(2)巧妙地利用了上階段的出礦川作充填、通風聯絡道，使工藝的通風和充填可靠性較好；(3)人員不進入空區作業，且回采間柱和礦房時，周邊均有較好的應力環境，安全性較好；(4)工藝簡單、可靠。

4 采空區穩定性離散元分析

4.1 離散單元法簡介

離散元法是一種基于牛頓第二定律的巖石力學計算分析方法，特別適合于節理比較發育的工程巖體，在采礦工程、隧道工程、地下工程等領域得到了廣泛應用[22]。UDEC 是針對非連續介質開發的二維離散元程序，在數學求解方式上采用了有限差分法，力學上則增加了對接觸面的非連續力學的模擬，被普遍用于研究非連續面占主導地位的工程問題[23-24]。

4.2 計算模型

根據采空區規模、圍巖狀況和原巖應力場等條件的差異，選擇Ⅰ號、Ⅱ號、Ⅲ號采空區分別代表較小、中等、較大規模的3種空區(Ⅰ號空區面積為300 m2，高度為6 m；Ⅱ號空區面積為900 m2，高度為14 m；Ⅲ號空區面積為3 600 m2，高度為22 m)。所有模型兩側面限制水平移動，底面限制豎向移動，頂部模擬上覆巖重，建立3種空區的離散元模型進行計算。數值模擬中礦巖體和節理的物理力學參數分別如表1和表2所示，3種空區的模型尺寸如表3所示，離散元計算模型如圖4所示。

4.3 計算結果分析

碎裂礦段中空區圍巖的破壞情況主要以塑性區分布的形式表現出來。3種空區圍巖塑性區分布如圖 5所示?？梢姡?/p>

(1)Ⅰ號采空區(圖5(a))頂板礁灰巖中間偏右1 m處有一狹長深度約為11 m的塑性區，其主要原因是空區底部碎裂礦體變形較大，頂板巖體受拉破壞造成。頂板條帶狀塑性區偏離中心，主要是節理組方位和地應力共同影響的結果；底板塑性區范圍較廣，但分布相對零散，深度在10 m以內。采空區兩側礦柱塑性區呈零散分布，且范圍較小，在4 m以內；由于規模小，Ⅰ號空區整體較為穩定。

(2)Ⅱ號采空區(圖5(b))頂、底板塑性區深度分別為7 m和6 m，分別位于左側頂板和右側底板。礦柱沒有出現塑性區，說明礦柱是穩定的。頂板多條塑性帶水平寬度大于1/2頂板長度，由于塑性帶尚未連通，不可能發生大規模冒路事故，但有可能會出現松動巖塊下落現象。隨著規模擴大，塑性區具有連通的趨勢，頂板穩定性變差。

(3)Ⅲ號采空區(圖5(c))頂、底板塑性區深度分別為20 m和15 m。頂板塑性區水平寬度大于頂板長度的1/2，由于塑性區已連通，大規模冒落事故的發生是必然的，礦柱未出現塑性區。

表1 礦(巖)體物理力學參數Table 1 Physical mechanical parameters of ore (rock)body

表2 各種節理力學參數Table 2 Physical mechanical parameters of joints

表3 計算模型尺寸Table 3 Calculation model sizes

圖4 空區穩定性離散元分析模型Fig.4 Discrete element analysis models of goafs stability

圖5 空區圍巖塑性區分布規律Fig.5 Distribution principle of surrounding plastic zone

總體來看，處于圍巖和碎裂礦段分界處的Ⅰ號采空區規模較小，頂板均質巖體中的塑性區呈線型，深度為10 m，空區總體上較穩定。完全分布在碎裂礦段內Ⅱ號和Ⅲ號采空區塑性區分布具有明顯的規律：塑性區分布均在空區頂板的左半部分，隨著空區規模增大，塑性區范圍增大，頂板穩定性變弱；礦柱均未出現塑性區；中等規模空區頂板塑性區深度為7 m，最大規?？諈^頂板塑性區深度為20 m。

5 采空區協同處理方案

采空區處理是控制和預防采空區災害的非常重要的部分。目前，采空區處理主要有充填、崩落、支撐、封閉等方法，其實質是緩和巖體應力集中的程度，轉移應力集中部位，或使應力達到新的相對平衡，從而控制地壓，促使安全生產[25-28]。

從碎裂礦段的賦存特點可知：賦存在其中的采空區如直接采用常規的空區處理方法，難以獲得系統最優。碎裂礦段采空區在未處理前，所在位置、賦存地質環境、空區規模和空區穩定性等因素是基本確定的，由采礦環境再造無間柱分段分條連續采礦法知：碎裂礦段內包括采場寬度和高度、階段高度和數目及部分切割工程及自由爆破空間等開采布局是可作適當調整的，這就為采空區協同處理奠定了良好的基礎。

5.1 較小或中等規?？諈^處理方案

當礦體內存在規模較小或中等規模采空區時，可根據空區存在的具體情況，采取適當的方法將空區調整為拉槽、拉底、自由爆破空間、分段鑿巖道的一部分加以利用。斜交采場長度方向較小規?？諈^處理方案如圖6所示?？梢姡号c采場長度方向呈角度a的較小規?？諈^調整為拉底空間加以利用。對于中等規?？諈^(如圖 7所示)，由于空區上部的塑性區分為 6～7 m，可將高度方向稍大的中等規?？諈^調整為自由爆可考慮將空區下部調整到下一階段，但深度一般不宜超過5 m，利用周邊隆口崩落部分礦體，崩落礦石堆積高度至出礦進路水平時，將剩余的上部空區調整為自由爆破空間加以利用。在選擇空區周邊的隆口進行崩礦作業時，根據碎裂礦段空區圍巖塑性區分布規律，應選擇主節理組下方的隆口更安全。在對空區周邊礦體進行打眼放炮等作業時，應根據空區圍巖塑性區分布規律選擇安全部位進行施工。

圖6 斜交采場長度方向較小規模空區處理方案Fig.6 Treatment scheme of small scale goaf intersecting with stope

圖7 高度方向稍大的中等規?？諈^處理方案Fig.7 Treatment scheme of middle scale goaf with a large size in height direction

5.2 較大規?？諈^處理方案

礦山大規模的采空區不多，Ⅲ號采空區是最大規模的采空區，目前，基本被崩塌的廢石及礦石充滿，所有能夠進入采空區的坑道都已被封堵。對該類型采空區，不可能在此基礎上再布置采場，只能根據采空區規模將其視為已經回采的采場或采區。在采空區區域較大或相當于采場規模，應回采空區周邊礦體時，需要在空區精密探測的基礎上，采用空區上部打充填井或鉆孔充填的辦法完成對空區的治理，必要時，還可以用加壓注漿的辦法增強充填體的穩定性。

6 結論

(1)為在生產高效和作業安全上取得系統整體最優，提出了將碎裂礦段開采與空區處理協同理念。

(2)針對廣西高峰 105號礦體碎裂礦段賦存特點和開采技術條件，基于采礦環境再造和連續開采理論，提出了一種適合碎裂礦段開采的新方法，即采礦環境再造無底柱分段分條連續采礦法。該方法巧妙地利用了將上階段的出礦川作充填、通風聯絡道，采礦和充填綜合生產能力高，貧化損失率低，作業安全性好，適應能力強，是一種開采類似礦體的較好采礦方法。

(3)對于完全嵌入在碎裂礦段種的空區，塑性區分布均在空區頂板的左半部分，礦柱上無塑性區；小規模空區整體上穩定，僅可能出現零星掉塊現象；中等規?？諈^頂板塑性區范圍為7 m，塑性帶跨度大于頂板長度的1/2，但未連通，不可能發生大規模冒落事故，但有可能會出現松動巖塊下落現象；大規模空區頂板塑性區為20 m；塑性區跨度大于頂板長度的1/2，塑性區已連通，大規模冒落事故不可避免。

(4)在采空區穩定性分析的基礎上，靈活調整開采布局，將采空區調整為部分切割工程或自由爆破空間，從而確定各種規?？諈^的協同處理方案。

(5)對碎裂礦段采礦方法與空區處理進行協同研究可為類似礦床的安全高效開采提供參考。

[1]趙勇, 張花, 張兆林, 等. 上向分層膠結充填采礦法在侯莊礦區的應用[J]. 礦業快報, 2006, 445(6): 38-40.ZHAO Yong, ZHANG Hua, ZHNAG Zhao-lin, et al. Application of upwards lamination stoping-cemented fill method at Houzhuang Mining Area[J]. Express Information of Mining Industry, 2006, 445(6): 38-40.

[2]谷新建, 胡磊, 陳澤呂, 等. 新龍礦業公司極破碎難采礦體采礦方法優化研究[J]. 礦業研究與開發, 2007, 27(4): 1-2, 89.GU Xin-jian, HU Lei, CHEN Ze-lü, et al. Optimization of mining method for extremely fractured difficult-to-mine orebody in Xinlong mineral Co. Ltd[J]. Mining Research and Development, 2007, 27(4): 1-2, 89.

[3]黃應盟. 廣西高峰礦業公司深部礦體采礦方法的改進[J]. 礦業研究與開發, 2002, 22(3): 12-14.HUAND Ying-meng. A modification on the mining method for the deep orebody in Guangxi Gaofeng mining Ltd[J]. Mining Research and Development, 2002, 22(3): 12-14.

[4]周旭, 陽雨平, 王貽明, 等. 上向水平分層充填法的優化研究[J]. 礦業快報, 2008, 472(8): 63-66.ZHOU Xu, YANG Yu-ping, WANG Yi-ming, et al. Research on optimization of flat-back cut-and-fill stoping method[J]. Express Information of Mining Industry, 2008, 472(8): 63-66.

[5]李興尚, 吳法春, 許家林. 上向水平分層充填采礦法的優化研究[J]. 金屬礦山, 2006, 358(4): 1-3, 6.LI Xing-shang, WU Fa-chun, XU Jia-lin. Research on optimization of flat-back cut-and-fill stoping method[J]. Metal Mine, 2006, 358(4): 1-3, 6.

[6]吳賢振, 饒運章, 熊正明. 銅山銅礦分段鑿巖階段礦房嗣后充填采礦法工藝優化研究[J]. 南方冶金學院學報, 2003, 24(4):1-4.WU Xian-zhen, RAO Yun-zhang, XIONG Zheng-ming. The study on the improvement of the technology of sublevel stoping and filling mining method[J]. Journal of Southern Institute of Metallurgy, 2003, 24(4): 1-4.

[7]袁志安. 大空區下連續開采方案研究[J]. 有色礦冶, 2006,22(3): 13-16.YUAN Zhi-an. Continuously exploitation project research of large empty area[J]. Non-ferrous Mining and Metallurgy, 2006,22(3): 13-16.

[8]任鳳玉, 李楠. 團城鐵礦多空區礦體開采技術研究[J]. 金屬礦山, 2008, 381(3): 32-34, 109.REN Feng-yu, LI Nan. Research of mining techniques for orebodies with multi-m ined-out areas in Tuancheng iron deposit[J]. Metal Mine, 2008, 381(3): 32-34, 109.

[9]周科平, 高峰, 胡建華, 等. 頂板誘導崩落預裂鉆孔裂隙發育監測與分析[J]. 巖石力學與工程學報, 2007, 26(5): 1034-1040.ZHOU Ke-ping, GAO Feng, HU Jian-hua, et al. Monitoring and analysis of fracture development pre-spliting hole of cave inducement of roof[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(5): 1034-1040.

[10]HU Jian-hua, ZHOU Ke-ping, LI Xi-bing, et al. Numerical analysis of application for induction caving roof[J]. Journal of Central South University of Technology, 2005, 12(Suppl. 1):146-149.

[11]GAO Feng, ZHOU Ke-ping, DONG Wei-jun, et al. Similar material simulation of time series system for induced caving of roof in continuous mining under backfill[J]. Journal of Central South University of Technology, 2008, 15(3): 356-360.

[12]Haken H. Synergetics: From pattern formation to pattern analysis and pattern recognition[J]. International Journal of Bifurcation and Chaos, 1994, 4(5): 1069-1083.

[13]Daffershofer A, Haken H. Synergetic computers for pattern recognition-a new approach to recognition of deformed patterns[J]. Pattern Recognition, 1994, 27(12): 1697-1705.

[14]吳永紅, 蘇懷智, 高培偉. 混凝土大壩裂縫光纖監測關鍵性基本問題的協同研究[J]. 水力發電學報, 2007, 26(4): 120-123.WU Yong-hong, SU Huai-zhi, GAO Pei-we. Integrated study on key basic problems of fiber optic monitoring cracks in concrete large dams[J]. Journal of Hydroelectric Engineering, 2007, 26(4):120-123.

[15]鄧金燦. 高峰礦深部開采工程地質調查與穩定性分析[J]. 采礦技術, 2001, 1(3): 66-67.DENG Jin-can. Stability analysis and engineering geological investigation of deep mining ore body in Gaofeng mine[J].Mining Technology, 2001, 1(3): 66-67.

[16]吳愛祥, 韓斌, 古德生, 等. 國內外地下金屬礦山連續開采技術研究的發展[J]. 礦冶工程, 2002, 22(3): 7-10.WU Ai-xiang, HAN Bin, GU De-sheng, et al. Progress in research on continuous mining technology of domestic and foreign underground metal mines[J]. Mining and Metallurgical Engineering, 2002, 22(3): 7-10

[17]鄧建, 李夕兵, 古德生. 無間柱連續采礦法礦段回采的地壓規律與控制技術[J]. 中國有色金屬學報, 2001, 11(4): 666-670.DENG Jian, LI Xi-bing, GU De-sheng. Ground pressure and control techniques in non-pillar continuous mining method[J].The Chinese Journal of Nonferrous Metal, 2001, 11(4): 666-670.

[18]古德生, 鄧建, 李夕兵. 地下金屬礦山無間柱連續采礦可靠性分析與設計[J]. 中國工程科學, 2001, 3(1): 51-57.GU De-sheng, DENG Jian, LI Xi-bing. Reliability analysis and design of non-pillar continuous mining in underground metal mines[J]. Engineering Science, 2001, 3(1): 51-57.

[19]古德生, 李夕兵. 現代金屬礦床開采科學技術[M]. 北京:冶金工業出版社, 2006.GU De-sheng, LI Xi-bing. Modern mining science and technology for metal mineral resources[M]. Beijing:Metallurgical Industry Press, 2006.

[20]周科平, 高峰, 古德生. 采礦環境再造與礦業發展新思路[J].中國礦業, 2007, 16(4): 34-36.ZHOU Ke-ping, GAO Feng, GU De-sheng. Mining environment regenerating and new thoughts on the development of mining industry[J]. China Mining Magazine, 2007, 16(4): 34-36.

[21]周科平, 朱和玲, 肖雄, 等. 采礦環境再造連續開采地壓演化過程的控制與仿真[J]. 中南大學學報: 自然科學版, 2008,39(3): 417-422.ZHOU Ke-ping, ZHU He-ling, XIAO Xiong, et al. Control and simulation of ground pressure evolutional process based on reconstructed mining environment and continuous caving method[J]. Journal of Central South University: Science and Technology, 2008, 39(3): 417-422.

[22]王泳嘉, 邢紀波. 離散單元法及其在巖土工程中的應用[M].沈陽: 東北工學院出版社, 1991.WANG Yong-jia, XING Ji-bo. The discrete element method and its application in soil and rock mechanics[M]. Shenyang:North-Eastern Institute of Engineering press, 1991.

[23]Sung O C, So-Keul C. Stability analysis of jointed rock slopes using the Barton-Bandis constitutive model in UDEC[J].International Journal Rock Mechanics and Mining Sciences,2004, 41(Suppl. 1): 581-586.

[24]Rajinder B, Kaare H. Parametric study on a large cavern in jointed rock using a distinct element model(UDEC-BB)[J].International Journal Rock Mechanics and Mining Sciences,1998, 35(1): 17-29.

[25]宮鳳強, 李夕兵, 董隴軍, 等. 基于未確知測度理論的采空區危險性評價研究[J]. 巖石力學與工程學報, 2008, 27(2):323-330.GONG Feng-qiang, LI Xi-bing, DONG Long-jun, et al.Underground goaf risk evaluation based on uncertainty measurement theory[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(2): 323-330.

[26]李俊平, 彭作為, 周創兵, 等. 木架山采空區處理方案研究[J].巖石力學與工程學報, 2004, 23(22): 3884-3890.LI Jun-ping, PENG Zuo-wei, ZHOU Chuang-bing, et al. Study on schemes for disposing abandoned stope in Mujia hill[J].Chinese Journal of Rock Mechanics and Engineering, 2004,23(22): 3884-3890.

[27]閆長斌, 徐國元, 中國生. 復雜地下空區綜合探測技術研究及其應用[J]. 遼寧工程技術大學學報, 2005, 24(4): 481-484.YAN Chang-bin, XU Guo-yuan, ZHONG Guo-sheng. Research of composite prospecting and its application in complicated underground mined-out areas[J]. Journal of Liaoning Technical University, 2005, 24(4): 481-484.

[28]劉占魁, 劉寶許, 來興平. 內蒙古霍各氣銅礦采空區穩定性分析與處理[J]. 北京科技大學學報, 2003, 25(5): 391-393.LIU Zhan-kui, LIU Bao-xu, LAI Xing-ping. Stability analysis and disposal of Huogeqi copper mine gob in inner Mongolia[J].Journal of University of Science and Technology Beijing, 2003,25(5): 391-393.