程 蕾,郝增超,曾維華(北京師范大學環(huán)境學院,北京 100875)

基于生態(tài)網(wǎng)絡分析法的煤化工水網(wǎng)絡結構解析

程 蕾,郝增超,曾維華*(北京師范大學環(huán)境學院,北京 100875)

由于煤化工水網(wǎng)絡結構的復雜性,水網(wǎng)絡的用水/排水去向以及結構的整體性分析已成為系統(tǒng)集成優(yōu)化的基礎工作.選取某典型煤制烯烴項目的水網(wǎng)絡體系為研究案例,采用生態(tài)網(wǎng)絡分析法深度解析煤化工水網(wǎng)絡體系的結構和用水效率.采用生態(tài)分析方法中的通量分析和效用分析對選取的煤化工水網(wǎng)絡結構進行研究.通量分析的結果顯示,在該煤化工水網(wǎng)絡結構中,通量占比最大的節(jié)點分別為第一冷卻水站(35.70%)和第二冷卻水站(26.35%).效用分析中用網(wǎng)絡共生指數(shù)(M)來衡量系統(tǒng)結構的穩(wěn)定性,計算得到該案例中M值為0.97 ,接近于1.從生態(tài)學的視角來看,網(wǎng)絡中各節(jié)點之間呈現(xiàn)競爭與共生的強度相當,網(wǎng)絡協(xié)作性較好.

生態(tài)網(wǎng)絡分析；煤化工；水網(wǎng)絡；結構解析；共生指數(shù)

煤化工行業(yè)的發(fā)展一直因為其高水耗、高污染的特點受到水資源的制約問題日益嚴重.我國煤炭資源和水資源總體呈逆向分布,由于產(chǎn)業(yè)布局受煤炭資源主導,使得產(chǎn)業(yè)發(fā)展中水資源配置的問題尤為凸顯.目前我國水資源和水環(huán)境問題成為煤化工產(chǎn)業(yè)發(fā)展最關鍵的瓶頸,應運而生的“廢水零排放”的目標又因為煤化工行業(yè)存在的一系列技術經(jīng)濟等各方面的問題而備受爭議[1-2].

目前關于煤化工行業(yè)的節(jié)水減排工作主要是關注煤化工工藝的優(yōu)化、廢水深度處理技術的開發(fā)從而來提高廢水的回用率[3],缺乏整體全局性的分析.從整體上關于煤化工水網(wǎng)絡體系優(yōu)化的研究主要是水網(wǎng)絡集成優(yōu)化理論[4-5],這方面的理論研究已經(jīng)比較成熟,然而實際應用相對較少.楊友麒等[6-7]提出的基于過程系統(tǒng)工程的節(jié)水三步曲是目前這方面的顯著研究成果之一,但是該方法體系全套應用成本相對龐大,還需在煤化工企業(yè)中進一步推廣.實施整體性的水網(wǎng)絡體系的集成優(yōu)化,需要對水網(wǎng)絡結構進行深入的剖析和解讀,從系統(tǒng)的角度對煤化工項目的整個水網(wǎng)絡結構進行整體性分析,考慮工程中所有不同水質的水流動過程,以圖的形式將整個流動過程表示出來,從而形成一個復雜的拓撲結構譜圖.同時需要采用一些結構分析的方法如結構解析模型、生態(tài)網(wǎng)絡分析等,結合投入產(chǎn)出以及生態(tài)學方面的理論,對煤化工水網(wǎng)絡結構進行逐步剖析.其中,生態(tài)網(wǎng)絡分析最早是起源于生態(tài)系統(tǒng)的結構分析方法[8],目前廣泛應用于城市能量流動以及社會經(jīng)濟[9-10]等方面的社會經(jīng)濟網(wǎng)絡結構分析.在工業(yè)領域,生態(tài)網(wǎng)絡分析方法多應用于生態(tài)工業(yè)園區(qū)的物質流動分析和穩(wěn)定性分析[11-13].

煤化工水網(wǎng)絡體系是一個龐大的系統(tǒng),各個組分之間有著復雜的直接或間接的聯(lián)系,呈現(xiàn)出的網(wǎng)絡結構可類比于環(huán)境生態(tài)網(wǎng)絡,在網(wǎng)絡內部進行著新陳代謝.水在代謝過程中出現(xiàn)各種變化,一旦某個組分的不利作用增強(即不利于水量的補充),整個水代謝系統(tǒng)就會出現(xiàn)紊亂;同時,與其相關的其它組分的有利作用如果可以抵消不利作用,則整個系統(tǒng)可以正常代謝實現(xiàn)水量的補充.本文基于此,從生態(tài)學的角度運用生態(tài)網(wǎng)絡分析法研究煤化工水網(wǎng)絡體系的網(wǎng)絡結構及用水效率,系統(tǒng)甄別現(xiàn)有煤化工企業(yè)內部結構存在的問題,從而為煤化工水網(wǎng)絡體系優(yōu)化研究開拓新的視野,為煤化工行業(yè)的節(jié)水減排工作提供理論支持和技術借鑒.

1 數(shù)據(jù)來源與研究方法

1.1 數(shù)據(jù)來源

本文主要從某個煤制烯烴項目的企業(yè)內部水平衡入手,建立基于水量變化的水網(wǎng)絡結構圖,分析網(wǎng)絡結構中各個節(jié)點之間的間接的潛在關系,以及各節(jié)點消耗水量的占比以及整體網(wǎng)絡體系穩(wěn)定性和效率.本文數(shù)據(jù)主要來自對該項目的實地調研,采用水平衡核算的方法最終得到具體各股水流的通量.

1.2 生態(tài)網(wǎng)絡分析法

網(wǎng)絡分析的基礎是鄰接矩陣[14]An×n=(aij)和流量平衡矩陣Fn×n=(fij),其中 n為網(wǎng)絡節(jié)點數(shù),若aij=1意味著節(jié)點i到j之間是有邊連接的,此時fij代表它們之間的單位時間內流量的大小.在網(wǎng)絡中,還存在著輸入、輸出和節(jié)點的內部存儲,分別以矢量 z,y,x來表示.那么,網(wǎng)絡的平衡約束可以表示為[15]：

式中：Tj即為通過節(jié)點j的凈輸入(輸出)量,任意節(jié)點的總輸入量和總輸出量是相等的.根據(jù)鄰接矩陣或流量平衡矩陣的不同迭代方式,網(wǎng)絡分析分為結構分析、通量分析、儲存分析和利用分析4個子類別.

1.2.1 通量分析 通量分析的基礎是流量密度矩陣G.通過G的迭代可以獲得兩點之間通過一定長度路徑的間接流量密度和綜合流量密度[16].

式中：Tj為通過節(jié)點j的凈通量.G0表示單位矩陣,它意味著流從某一節(jié)點出又回到相同的節(jié)點中;矩陣G1意味著網(wǎng)絡中的任何節(jié)點之間都是直接流;Gm(m≥2)意味著節(jié)點之間存在 m長度的間接流,N即為兩點之間通過長度為m的綜合流量密度矩陣[15].

無量綱綜合流矩陣 N通過右乘對角化通量矢量(T)以轉換N的縱向量成為有量綱的矩陣Y.通過計算矩陣Y的每個縱向量總數(shù)和,矩陣Y的列向量即得出[17].

W可以反映出節(jié)點j對整個系統(tǒng)的貢獻情況,從而得到每個節(jié)點的貢獻并以此確定每個節(jié)點的比重.

1.2.2 效用分析 網(wǎng)絡效用分析是用于表征網(wǎng)絡系統(tǒng)中不同組分間的定性定量關系[18],主要分析節(jié)點i和j之間的凈通量,即相互之間的利用程度.公式如式(6)：

對其進行迭代,則得到總利用密度矩陣U：

式中：D為直接效用矩陣,U為無量綱綜合效用矩陣.采用這些分析可用來研究生態(tài)網(wǎng)絡的一些主要性質,如：間接效應,放大效應、網(wǎng)絡均質性,以及網(wǎng)絡協(xié)作性.

指定矩陣 U的任意兩個節(jié)點之間的生態(tài)關系表示為 suij,則一對節(jié)點之間的直接關系可通過一對直接利用矩陣給出.如：(su21,su12) = (+,?)表示節(jié)點2利用節(jié)點1,如果(su21,su12) = (?,+),則是節(jié)點 2被節(jié)點 1利用.即 (0,0)表示中立狀態(tài), (?,?)表示競爭關系, (+,+)表示共生關系[19].效用分析的最終目的是研究系統(tǒng)的共生狀態(tài),引入共生指數(shù)M以定量系統(tǒng)的整個共生情況,J (U)代表積極關系數(shù)量與消極關系數(shù)量的比率.

其中：

如果 M＞1,則矩陣正號數(shù)目大于負號數(shù)目,說明系統(tǒng)的節(jié)點之間建立了更多積極的關系,各節(jié)點之間的共生互利關系較多,網(wǎng)絡協(xié)作性較好.相反,如果負號多于正號,則系統(tǒng)建立了更多消極的關系,節(jié)點之間競爭關系比較多,網(wǎng)絡協(xié)作性較差.

2 案例分析

2.1 網(wǎng)絡結構概況

本文選取的網(wǎng)絡結構圖主要依據(jù)是項目環(huán)評報告中的水平衡圖,并加以總結概括,從網(wǎng)絡分析的角度將水平衡圖中涉及到的所有工程環(huán)節(jié)概括為生產(chǎn)工藝環(huán)節(jié)和公用工程.整個水網(wǎng)絡系統(tǒng)(圖1)共包含20個網(wǎng)絡節(jié)點,主要分為生產(chǎn)工藝環(huán)節(jié)和公共工程,其中公用工程又包括循環(huán)水系統(tǒng)和供水環(huán)節(jié)以及其它公用裝置.在水網(wǎng)絡結構圖(圖2)中,主要供水環(huán)節(jié)中的生產(chǎn)水系統(tǒng)給生產(chǎn)工藝環(huán)節(jié)的節(jié)點提供新鮮補充用水,回用水系統(tǒng)主要是進行中水的深度處理,并為循環(huán)水系統(tǒng)以及除鹽水站提供補充用水,循環(huán)水系統(tǒng)主要為各生產(chǎn)環(huán)節(jié)提供循環(huán)冷卻水,各生產(chǎn)工藝環(huán)節(jié)則需要耗用一定的新鮮用水或者回用水,同時產(chǎn)生一定量的廢水,并且產(chǎn)生的中間廢水可以進入另一個生產(chǎn)環(huán)節(jié)進行重復利用,另外生產(chǎn)過程中每個工藝環(huán)節(jié)都會需要大量的循環(huán)冷卻水.因為這些不同水質和水量的水流在各個節(jié)點之間流入流出,各節(jié)點之間以直接或間接的方式互相關聯(lián).統(tǒng);

圖1 水網(wǎng)絡系統(tǒng)示意Fig.1 Thesystematic diagram of water network

圖2 水網(wǎng)絡結構示意Fig.2 The structural diagram of water network

2.2 通量分析

通過對該水網(wǎng)絡結構的直接流量矩陣迭代計算得到的比重分布如圖3所示.整體體系中占最大比重的節(jié)點分別為節(jié)點○14第一冷卻水站(,35.70%)、節(jié)點○15第二冷卻水站(26.35%),兩者之和比例占到了 62.05%,其次是節(jié)點⑩空分空壓環(huán)節(jié),占比為13.18%、節(jié)點⑨丙烯制冷+聚丙烯環(huán)節(jié)占比 10%.由網(wǎng)絡結構模型來看,第一冷卻水站提供5個主要工藝環(huán)節(jié)的循環(huán)冷卻水,提供的循環(huán)冷卻水量有65747t/h,占總循環(huán)冷卻水量41.4%.循環(huán)冷卻水是該水網(wǎng)絡體系中流量最大的一種水流.同時該環(huán)節(jié)處其它損失水量達到745t/h,需要回用水裝置補充水量278t/h,新鮮水補充量 1103t/h,并向清凈下水裝置排放污水338t/h.同樣第二冷卻水站為 5個主要工藝環(huán)節(jié)提供循環(huán)冷卻水 56456t/h,占總循環(huán)冷卻水量35.3%,該處其它損失量有 206t/h,排放廢水為285t/h,需要新鮮水補充水量 1137t/h.而節(jié)點⑨和節(jié)點⑩都是由2個子環(huán)節(jié)合并而成的,用到的循環(huán)冷卻水也就相對單個工藝環(huán)節(jié)較多,因為循環(huán)冷卻水的貢獻,從而最終占比也就較大.

圖3 水量變化各節(jié)點所占比重Fig.3 The proportion of each node during water variations

結合對水平衡網(wǎng)絡中循環(huán)冷卻水分布的初步分析可以發(fā)現(xiàn),各節(jié)點總體通量所占比重與循環(huán)冷卻水消耗比重排序具有一定的同步性,因為相對其它新鮮用水以及廢水排放量,循環(huán)冷卻水量遠遠大于新鮮用水或回用水量,那么循環(huán)水量對整體水代謝系統(tǒng)的通量影響也是最大的.也就是節(jié)水的關鍵首先應是減少循環(huán)冷卻水的用量,實施工程節(jié)水改造的首要應該是從第一冷卻水站、第二冷卻水站開始,進而再從整體上對水網(wǎng)絡進行優(yōu)化配置.

2.3 網(wǎng)絡協(xié)作性分析

通過對水網(wǎng)絡結構的直接流量矩陣的迭代計算,得到無量綱綜合效用矩陣U,從而可通過該效用矩陣的符號矩陣表 1來表征各個節(jié)點之間的生態(tài)關系.

在該煤化工水網(wǎng)絡系統(tǒng)中,如果兩個環(huán)節(jié)需要消耗同一個供水環(huán)節(jié)的來水,那么這兩個節(jié)點之間便存在著一定程度的競爭關系,同時用水環(huán)節(jié)與供水環(huán)節(jié)也呈現(xiàn)著利用關系,如果該用水環(huán)節(jié)產(chǎn)生的廢水經(jīng)過回收處理之后再被其它供水環(huán)節(jié)回收利用,那么此時用水環(huán)節(jié)與供水環(huán)節(jié)之間便呈現(xiàn)了共生關系.

在表1中,空白處的值均為0,表征的是對應的兩個節(jié)點之間不存在生態(tài)關系中的任何一種關系((suij,suji) = (0,0)),也即③甲醇單元、⑨丙烯制冷+聚丙烯、⑩空分+空壓、○12酸堿站4個節(jié)點都與其他節(jié)點不存在生態(tài)關系.由網(wǎng)絡結構圖2也可以看出節(jié)點③甲醇單元、節(jié)點⑩空分+空壓均只與節(jié)點○15第二冷卻水站交換循環(huán)冷卻水,與其它節(jié)點之間即沒有連接,同樣節(jié)點⑨丙烯制冷+聚丙烯只與節(jié)點○14第一冷卻水站交換循環(huán)冷卻水,節(jié)點○12酸堿站只與○16第三冷卻水站交換循環(huán)冷卻水.其它16個節(jié)點相互形成120對生態(tài)關系,其中處于共生關系(淺灰色填充標出, (suij,suji) = (+,+))的有30對,處于競爭關系(深灰色填充, (suij,suji) = (-,-))的有32對,其他的則都為利用關系((suij,suji) = (+,-)或(suij,suji) = (-,+)),共58對.由此可以得出在該煤化工水網(wǎng)絡體系中主要的生態(tài)關系是利用關系,實際中煤化工廢水回用體系是一個復雜的系統(tǒng),各個生產(chǎn)環(huán)節(jié)之間聯(lián)系緊密,在水量與水質的變化過程中互相影響.通過效用分析法得出的大量節(jié)點之間的生態(tài)關系為利用關系,從側面說明從側面說明該煤化工企業(yè)內部的水回用環(huán)節(jié)比較多.同時系統(tǒng)中競爭關系和共生關系的比例相當,均占總生態(tài)關系的25%左右,一方面在該水網(wǎng)絡體系中的網(wǎng)絡節(jié)點主要都是用水環(huán)節(jié),用水環(huán)節(jié)之間必然因為新鮮補充水、循環(huán)水之間的分配表現(xiàn)出競爭關系.而因為生產(chǎn)環(huán)節(jié)之間的廢水的直接或間接的回用,相關的網(wǎng)絡節(jié)點呈現(xiàn)生態(tài)學意義上的共生關系.由此計算出共生指數(shù),M=J(U)=S+/S =128/132= 0.97＜1,該數(shù)值十分接近于 1,從生態(tài)學的視角來看網(wǎng)絡中呈現(xiàn)競爭與共生的強度相當,網(wǎng)絡的協(xié)作性比較好.但是共生指數(shù)值小于 1,說明系統(tǒng)中各節(jié)點之間呈現(xiàn)的消極關系強于積極關系[23],而這也符合了煤化工企業(yè)耗水量大,耗水環(huán)節(jié)多的特征.

表1 水代謝系統(tǒng)中各節(jié)點關系Table 1 The relationship among various nodes in the water metabolism system

根據(jù)調研獲取的生產(chǎn)數(shù)據(jù),該煤化工企業(yè)生產(chǎn)過程中的噸產(chǎn)品水耗在50t/t左右,中水回用率73%,廢水回用率為 39.3%,相比于煤化工行業(yè)的噸產(chǎn)品水耗 24～74t/t[21],該煤制烯烴項目的耗水量很大,相比于煤化工行業(yè)的平均中水回用率60%,該項目的中水回用率比較高,但是廢水回用率比較低.綜合得出,該煤化工項目耗水量大,新鮮用水效率不夠高,雖然中水回用率已經(jīng)比較高,但是整體的廢水回用率不高,該煤化工企業(yè)沒有充分實現(xiàn)廢水的循環(huán)多級利用,做到低水低用、高水高用.

煤化工行業(yè)的節(jié)水減排是一項復雜的系統(tǒng)工程,進行水平衡分析摸清工業(yè)的用水分布是實施節(jié)水工作的基礎,具體企業(yè)的節(jié)水應該是在項目的設計階段就開始實施節(jié)水規(guī)劃.結合本案例分析結果,該煤化工企業(yè)的主要耗水都是循環(huán)冷卻水站的水分蒸發(fā)量,如果采用空冷技術的將大大減低循環(huán)水蒸發(fā)損失量.廢水回用率的提高不僅需要完善用水工藝,更需要廢水的適當處理,針對性的處理特定限制因素污染物,避免回用水的過度處理.

3 結論

3.1 生態(tài)網(wǎng)絡分析方法作為一種投入產(chǎn)出分析法,適用于煤化工企業(yè)復雜系統(tǒng)的水網(wǎng)絡結構解析,既可以科學而詳細的分析出企業(yè)內部的主要耗水環(huán)節(jié),又可以整體上評價企業(yè)的用水效率和穩(wěn)定性.

3.2 通量分析的結果顯示,在該煤化工水網(wǎng)絡結構中,通量占比最大的節(jié)點分別為節(jié)點○14第一冷卻水站(35.70%)、節(jié)點○15第二冷卻水站(26.35%).說明該煤化工企業(yè)節(jié)水的關鍵首先應是減少循環(huán)冷卻水的用量,實施工程節(jié)水改造的首要應該是從第一冷卻水站、第二冷卻水站開始,以及從整體上對其的優(yōu)化配置.

3.3 效用分析結果顯示該煤化工水網(wǎng)絡結構體系的共生指數(shù)為 0.97,數(shù)值略微小于 1,說明該體系的網(wǎng)絡協(xié)作性很好,整個水網(wǎng)絡體系比較穩(wěn)定.水網(wǎng)絡體系中呈現(xiàn)的主要生態(tài)關系是利用關系,從側面說明該煤化工企業(yè)內部的水回用環(huán)節(jié)比較多,結合實際工程狀況,可知這主要是循環(huán)水的重復利用導致的.

3.4 該煤化工企業(yè)的主要耗水都是循環(huán)冷卻水站的水分蒸發(fā)量,建議可采用空冷技術,減低循環(huán)水蒸發(fā)損失量.煤化工企業(yè)要減少新鮮水耗,提高廢水回用率,不僅需要進行用水工藝的優(yōu)化,更重要的是從全局角度進行水系統(tǒng)集成優(yōu)化管理.

[1]曲風臣.煤化工廢水“零排放”技術要點及存在問題 [J]. 化學工業(yè), 2013,31(2/3)：18-24.

[2]童莉.煤化工廢水零排放的制約性問題 [J]. 化工環(huán)保, 2010, 30(5)：371-375.

[3]任慶偉,商登峰,劉青海,等.新型煤化工企業(yè)水處理工藝方案設計優(yōu)化 [J]. 工業(yè)水處理, 2010,30(7)：87-90.

[4]徐宏勇.煤化工行業(yè)水系統(tǒng)集成優(yōu)化及循環(huán)水軟化技術研究[D]. 華東理工大學, 2014.

[5]楊友麒,莊芹仙.節(jié)水減排的過程系統(tǒng)工程方法 [J]. 現(xiàn)代化工, 2008,28(1)：8-11.

[6]楊友麒.節(jié)能減排的全局過程集成技術的研究與應用進展 [J].化工進展, 2009,28(4)：541-548.

[7]宏曉晶,劉雪鵬,吳盛文,等.過程系統(tǒng)工程方法在煤化工節(jié)水優(yōu)化中的應用 [J]. 工業(yè)水處理, 2015,35(8)：107-109.

[8]Fath B D, Patten B C. Review of the foundations of network environ analysis [J]. Ecosystems, 1999,2(2)：167-179.

[9]Lu W, Su M, Zhang Y, et al. Assessment of energy security in China based on ecological network analysis： A perspective from the security of crude oil supply [J]. Energy Policy, 2014,74：406-413.

[10]Zhang Y, Yang Z, Fath B D, et al. Ecological network analysis of an urban energy metabolic system： Model development, and a case study of four Chinese cities [J]. Ecological Modelling, 2010, 221(16)：1865-1879.

[11]Yu C, de Jong M, Dijkema G P J. Process analysis of eco-industrial park development-the case of Tianjin, China [J]. Journal of Cleaner Production, 2014,64：464-477.

[12]Lu Y, Su M, Liu G, et al. Ecological Network Analysis for a Low-Carbon and High-Tech Industrial Park [J]. The Scientific World Journal, 2012,2012：1-9.

[13]Zhang Y, Zheng H, Fath B D. Ecological network analysis of an industrial symbiosis system： A case study of the Shandong Lubei eco-industrial park [J]. Ecological Modelling, 2015,306：174-184.

[14]Zhang Y, Zheng H, Chen B, et al. Social network analysis and network connectedness analysis for industrial symbiotic systems：model development and case study [J]. Frontiers of Earth Science, 2013,7(2)：169-181.

[15]Fang D, Fath B D, Chen B, et al. Network environ analysis for socio-economic water system [J]. Ecological Indicators, 2014, 47：80-88.

[16]Fath B D, Patten B C. Quantifying resource homogenization using network flow analysis [J]. Ecological Modelling, 1999, 123(2/3)：193-205.

[17]Zhang Y, Zheng H, Fath B D, et al. Ecological network analysis of an urban metabolic system based on input-output tables： Model development and case study for Beijing [J]. Science of The Total Environment, 2014,468-469：642-653.

[18]Fath B D, Patten B C. Network synergism： Emergence of positive relations in ecological systems [J]. Ecological Modelling, 1998, 107(2/3)：127-143.

[19]Fath B D. Network mutualism： Positive community-level relations in ecosystems [J]. Ecological Modelling, 2007,208(1)：56-67.

[20]Zhang Y, Zheng H, Fath B D. Ecological network analysis of an industrial symbiosis system： A case study of the Shandong Lubei eco-industrial park [J]. Ecological Modelling, 2015,306(SI)：174-184.

[21]楊友麒.三步法可使煤化工水耗大幅度下降-國家煤化工網(wǎng)[EB/OL].http：//www.coalchem.org.cn/publication/html/80020811 /577.html.2017-04-17.

Analysis of the water network structure of coal chemical industry based on ecological network analysis.

CHENGLei, HAO Zeng-chao, ZENG Wei-hua*(School of Environment, Beijing Normal University, Beijing 100875, China). China Environmental Science, 2017,37(6)：2363～2368

Due to the complicity of the water network,the systematic analysis of water use/drainage and structure of the water network is fundamental for integrated system optimizations.In this paper, the water network system of a typical coal-to-olefins project was selected as a case study to analyse the structure and water use efficiency of the water network system based on the ecological network analysis method. Specifically, the network flux analysis and utility analysis were employed for the structure analysis of the selected water network. Results of flux analysis showed that the first cooling water station had the largest proportion of 35.70%, followed by the second cooling water station with a proportion of 26.35%. For the utility analysis, the networkmutualism index (M) was adapted as the measure of the network structural stability. For the case study, the M value was obtained as 0.97, which was very close to one. From the perspective of ecology, the strength of competition and symbiosis was approximately equivalent, indicating strong collaboration among nodes in this network.

ecological network analysis；coal chemical；water network；structure analysis；mutualism index

X703

A

1000-6923(2017)06-2363-06

程 蕾(1993-),女,湖北黃岡人,北京師范大學碩士研究生,主要從事環(huán)境規(guī)劃與管理研究.

2016-09-19

典型脆弱生態(tài)修復與保護研究重大專項(2016YFC050350403-1)

* 責任作者, 教授, zengwh@bnu.edu.cn