基于耗散結(jié)構(gòu)的礦山生態(tài)環(huán)境保護(hù)Busselator模型研究

劉烈武，宋煥斌，楊 偉

（1.昆明理工大學(xué)國(guó)土資源工程學(xué)院，云南 昆明 650093；2.紅河學(xué)院商學(xué)院，云南 蒙自 661100）

礦山生態(tài)環(huán)境的污染破壞已成為影響我國(guó)社會(huì)經(jīng)濟(jì)發(fā)展的重要因素。我國(guó)現(xiàn)有國(guó)有礦山企業(yè)8000多個(gè)，個(gè)體礦山企業(yè)12萬多個(gè)，這樣大規(guī)模的礦產(chǎn)資源開發(fā)對(duì)礦山周圍的生態(tài)環(huán)境造成的破壞是相當(dāng)驚人的。因此，加強(qiáng)礦山生態(tài)環(huán)境的保護(hù)具有重要的現(xiàn)實(shí)意義。學(xué)者們對(duì)礦山生態(tài)環(huán)境保護(hù)與治理的機(jī)制及對(duì)策進(jìn)行了大量研究［1－3］，對(duì)礦山生態(tài)環(huán)境保護(hù)做出了積極的貢獻(xiàn)。但現(xiàn)有研究主要集中在對(duì)策方面的定性研究，缺乏對(duì)礦山生態(tài)環(huán)境保護(hù)的基礎(chǔ)理論方面的研究。由此，礦山生態(tài)環(huán)境保護(hù)缺乏理論方面的指導(dǎo)。

耗散結(jié)構(gòu)理論作為一個(gè)用于分析處理遠(yuǎn)離平衡的開放性系統(tǒng)的有力工具，已經(jīng)廣泛地應(yīng)用于諸多領(lǐng)域的研究［4－9］。耗散結(jié)構(gòu)理論中的Brusselator模型，由于能模擬廣泛的系統(tǒng)自組織行為和系統(tǒng)演化動(dòng)力行為，在資源利用、經(jīng)濟(jì)管理等方面得到了廣泛應(yīng)用。如張鐵男等利用耗散結(jié)構(gòu)的Brusselator模型研究企業(yè)的生存能力［10］；魏遙、雷良海利用Brusselator模型研究產(chǎn)融集團(tuán)演化［11－12］。耗散結(jié)構(gòu)理論及Brusselator模型為解決經(jīng)濟(jì)、管理及資源利用等方面的問題提供了新的思路和理論依據(jù)［13］。基于此，本文利用耗散結(jié)構(gòu)理論分析礦山生態(tài)環(huán)境保護(hù)的理論機(jī)制，通過對(duì)Brusselator模型的變量進(jìn)行轉(zhuǎn)義來定量研究礦山生態(tài)環(huán)境保護(hù)的耗散結(jié)構(gòu)的形成，從理論上和定量分析上研究礦山生態(tài)環(huán)境保護(hù)。

1 礦山生態(tài)環(huán)境保護(hù)的耗散結(jié)構(gòu)演化機(jī)理

耗散結(jié)構(gòu)理論是由比利時(shí)科學(xué)家伊里亞·普里戈津（I.Prigogine）20世紀(jì)70年代提出的。他把非平衡結(jié)構(gòu)持續(xù)不斷地從外界交換物質(zhì)和能量，通過相變轉(zhuǎn)變?yōu)椤坝行颉苯Y(jié)構(gòu)稱為耗散結(jié)構(gòu)。耗散結(jié)構(gòu)是一種“活”結(jié)構(gòu)，需不斷地與外界交換物質(zhì)、能量來維持其有序結(jié)構(gòu)。耗散結(jié)構(gòu)的形成需具備四個(gè)條件：①系統(tǒng)是開放的；②外界對(duì)系統(tǒng)的作用必須達(dá)到一定的閾值；③系統(tǒng)遠(yuǎn)離平衡狀態(tài)；④遠(yuǎn)離平衡態(tài)的漲落導(dǎo)致系統(tǒng)向耗散結(jié)構(gòu)轉(zhuǎn)換。自然中，廣泛地存在著這種有序的耗散結(jié)構(gòu)。如，天上的云在氣流或溫差等外力作用下在藍(lán)天上形成藍(lán)白相間條紋的“天街”；液體持續(xù)加熱后，形成微小對(duì)流形式的六角形花紋結(jié)構(gòu)等。

礦山生態(tài)環(huán)境保護(hù)實(shí)際上也是一個(gè)系統(tǒng)。礦產(chǎn)資源開采促進(jìn)社會(huì)經(jīng)濟(jì)的發(fā)展，但受社會(huì)經(jīng)濟(jì)規(guī)律的制約和生態(tài)環(huán)境和自然條件的影響。礦山周圍的生態(tài)環(huán)境服從自然規(guī)律，但受到人類社會(huì)活動(dòng)的嚴(yán)重影響。礦山、地方社會(huì)經(jīng)濟(jì)和周圍的生態(tài)環(huán)境三者之間相互聯(lián)系、相互影響，相互制約，構(gòu)成了一個(gè)復(fù)雜的巨系統(tǒng)。這個(gè)復(fù)雜的系統(tǒng)是個(gè)開發(fā)的系統(tǒng)。礦山與其他二者之間有序、共贏、協(xié)調(diào)發(fā)展，我們就可以說礦山生態(tài)環(huán)境得到很好的保護(hù)。因此，礦山生態(tài)環(huán)境保護(hù)，實(shí)質(zhì)上就是礦山與社會(huì)經(jīng)濟(jì)、生態(tài)環(huán)境所構(gòu)成的系統(tǒng)動(dòng)態(tài)演化為耗散結(jié)構(gòu)的過程。根據(jù)耗散結(jié)構(gòu)理論，系統(tǒng)內(nèi)部的矛盾因素相互作用，總是產(chǎn)生“正熵”，“正熵”的增加使系統(tǒng)更加無序混亂；外部環(huán)境對(duì)系統(tǒng)輸入“負(fù)熵”，“負(fù)熵”可以抵消系統(tǒng)內(nèi)部的“正熵”，從而能使系統(tǒng)通過漲落作用遠(yuǎn)離平衡狀態(tài)，進(jìn)而演化為耗散結(jié)構(gòu)。礦山生態(tài)環(huán)境保護(hù)的耗散結(jié)構(gòu)演化機(jī)理可看成如圖1所示的過程。

圖1 礦山生態(tài)環(huán)境保護(hù)耗散結(jié)構(gòu)形成過程

系統(tǒng)內(nèi)部的正熵H＋是因素集ui的相互作用產(chǎn)生，外部系統(tǒng)產(chǎn)生的負(fù)熵H－是因素集Vi相互作用，系統(tǒng)的總熵變dH是負(fù)熵抵銷正熵的程度?？傡刈凃?qū)動(dòng)著礦山生態(tài)系統(tǒng)從無序到低層次或從低層次有序向高層次有序變遷。系統(tǒng)的變遷能否突變?yōu)楦叨扔行虻暮纳⒔Y(jié)構(gòu)，需總熵變小到某個(gè)閾值。尋求礦山生態(tài)系統(tǒng)耗散結(jié)構(gòu)突變的閾值，分析內(nèi)外部要素對(duì)系統(tǒng)的熵變大小，才能針對(duì)性地采取對(duì)策，引導(dǎo)礦山生態(tài)環(huán)境保護(hù)取得成效。

2 礦山生態(tài)環(huán)境系統(tǒng)正負(fù)熵的指標(biāo)體系

影響礦山生態(tài)環(huán)境系統(tǒng)內(nèi)部的正熵產(chǎn)生的因素具有層次性和結(jié)構(gòu)復(fù)雜性，外部環(huán)境對(duì)生態(tài)環(huán)境保護(hù)的投入也是多方面的。為能全面、準(zhǔn)確地描述礦山生態(tài)環(huán)境系統(tǒng)的耗散結(jié)構(gòu)演變的作用要素，我們參考了大量的礦山生態(tài)環(huán)境保護(hù)的文獻(xiàn)，同時(shí)咨詢相關(guān)專家，遵循指標(biāo)的相關(guān)性、綜合性、可操作性原則，采用分類的方法，經(jīng)過層層篩選，系統(tǒng)綜合，最后建立了影響礦山生態(tài)環(huán)境系統(tǒng)的正負(fù)熵變的定性指標(biāo)體系，如表1所示。

3 礦山生態(tài)環(huán)境系統(tǒng)正負(fù)熵變計(jì)算

i＝1，2，…，m；為指標(biāo)的個(gè)數(shù)；j＝1，2，…，5

各指標(biāo)的權(quán)重Wi的確定可以采用專家打分法、層次分析法以及熵權(quán)法等。專家打分與層次分析方法都具有很強(qiáng)的主觀性。本文已經(jīng)計(jì)算出每個(gè)指標(biāo)的熵值，因此采用熵權(quán)法能較快地得到每個(gè)指標(biāo)的權(quán)重，計(jì)算方法如下。

指標(biāo)的差異系數(shù)

則指標(biāo)的權(quán)重為

由此，可得到指標(biāo)的正熵值分別為

同樣的方法可確定負(fù)熵變

表1 影響礦山生態(tài)系統(tǒng)正熵變的指標(biāo)體系

4 礦山生態(tài)環(huán)境保護(hù)耗散結(jié)構(gòu)形成的判據(jù)

礦山生態(tài)環(huán)境保護(hù)耗散結(jié)構(gòu)形成的定量分析的關(guān)鍵問題，是建立數(shù)學(xué)模型并找到礦山系統(tǒng)成為耗散結(jié)構(gòu)的熵變判據(jù)。普里高津提出的Brusselator模型為我們提供了方法論和耗散結(jié)構(gòu)形成的熵變閾值的判據(jù)［14］。

4.1 Brusselator模型的基本形式

普里高津提出的Brusselator三分子模型如式（5）所示。

式中：A、B為初始反應(yīng)物；D、E為反應(yīng)產(chǎn)物，他們保持不變；X、Y為中間組分，它們可以有隨時(shí)問變化的濃度。Brusselator動(dòng)力學(xué)模型的反應(yīng)——擴(kuò)散動(dòng)力學(xué)方程可寫作

該方程存在唯一的均勻的定態(tài)解

則成為耗散結(jié)構(gòu)的動(dòng)力學(xué)臨界值條件為

4.2 對(duì)Brusselator模型的變量進(jìn)行轉(zhuǎn)義

現(xiàn)將“布魯塞爾器”模型轉(zhuǎn)義，將其中A、B、D、E、X、Y轉(zhuǎn)變?yōu)榈V山生態(tài)保護(hù)的耗散結(jié)構(gòu)研究的概念。

設(shè)：A、B為熵變，A為系統(tǒng)內(nèi)部的正熵變，B為外界的負(fù)熵變；D、E正、負(fù)熵作用下的狀態(tài)，D為非耗散結(jié)構(gòu)，E為耗散結(jié)構(gòu)；X、Y為礦山生態(tài)系統(tǒng)的中間狀態(tài)，X為狀態(tài)Ⅰ，Y為中間狀態(tài)Ⅱ。

經(jīng)過上面的概念轉(zhuǎn)義，可建立礦山生態(tài)系統(tǒng)的Brusselator模型，具體的礦山生態(tài)系統(tǒng)耗散結(jié)構(gòu)的形成用“布魯塞爾器”描述如下。

模型在K1階段，礦山在持續(xù)生產(chǎn)、經(jīng)營(yíng)的同時(shí)，不斷地對(duì)周圍的生態(tài)環(huán)境進(jìn)行污染破壞，此時(shí)，礦山生態(tài)系統(tǒng)在內(nèi)部正熵作用下演化為中間狀態(tài)Ⅰ，在K2階段，外部環(huán)境對(duì)系統(tǒng)輸入物質(zhì)能量與信息作用等，系統(tǒng)生態(tài)環(huán)境得到一定的保護(hù)，系統(tǒng)演變?yōu)橹虚g狀態(tài)Ⅱ和非耗散結(jié)構(gòu)D；在K3階段，兩種中間狀態(tài)相互作用，導(dǎo)致中間狀態(tài)加劇演變。最終，在K4階段，中間狀態(tài)演變?yōu)榉€(wěn)定、有序的耗散結(jié)構(gòu)E。

4.3 礦山生態(tài)環(huán)境保護(hù)耗散結(jié)構(gòu)形成的判據(jù)

轉(zhuǎn)義后的Brusselator模型，很好地說明了礦山生態(tài)環(huán)境的動(dòng)態(tài)演化過程。因此，礦山耗散結(jié)構(gòu)的判據(jù)也應(yīng)該為B＞1＋A2。也就是說：外部環(huán)境的負(fù)熵B的輸入達(dá)到一定閾值時(shí)，礦山系統(tǒng)才突變?yōu)楹纳⒔Y(jié)構(gòu)；低于這個(gè)閾值，礦山系統(tǒng)處于一種臨界狀態(tài)，生態(tài)環(huán)境尚未得到有效的保護(hù)。因此，我們只要控制正負(fù)熵變，滿足上面的關(guān)系式就能夠使影響礦山生態(tài)系統(tǒng)的漲落起正向作用，向耗散結(jié)構(gòu)突變。

5 實(shí)證研究

本文選取云南紅河州某A礦山進(jìn)行實(shí)證研究。近年來，該礦制定和完善促進(jìn)生態(tài)和諧的綠化規(guī)劃，加大污水和尾礦的處理力度，提高生產(chǎn)裝備水平，加強(qiáng)礦容礦貌整治，計(jì)劃“十二五”期間，把礦山建設(shè)成綠色礦山，積極向國(guó)家級(jí)綠色礦山目標(biāo)邁進(jìn)。

以調(diào)查問卷的方式獲取基本數(shù)據(jù)。發(fā)出調(diào)查問卷50份，回收有效問卷42份。數(shù)據(jù)統(tǒng)計(jì)如表2所示。

表2 各正負(fù)熵變指標(biāo)調(diào)查統(tǒng)計(jì)表

根據(jù)式（1）～（4）計(jì)算礦山生態(tài)環(huán)境系統(tǒng)正負(fù)熵變值如表3所示。

由表3的計(jì)算，A＝H＋＝0.7206；B＝H－＝0.8134。

由此，B－（1＋A2）＝0.48134－（1＋0.72062）＝－0.7058

故B＜1＋A2。

由上計(jì)算可知，該礦山生態(tài)環(huán)境保護(hù)系統(tǒng)未達(dá)到耗散結(jié)構(gòu)狀態(tài)。究其原因，一方面是內(nèi)部要素產(chǎn)生的正熵變過大，正熵變偏高的要素主要有：礦產(chǎn)資源開發(fā)條件0.0519，交通條件0.0520，企業(yè)人員素質(zhì)0.0527，循環(huán)經(jīng)濟(jì)發(fā)展水平0.0622，監(jiān)督機(jī)制完善程度0.0582，企業(yè)環(huán)保機(jī)制和環(huán)保理念等。另一方面，是外部輸入系統(tǒng)的負(fù)熵偏低，礦山環(huán)境補(bǔ)償資金使用效益0.0337，完善礦山環(huán)境恢復(fù)保證金機(jī)制0.0499，產(chǎn)學(xué)研結(jié)合程度0.0395，生態(tài)環(huán)境保護(hù)資金投入0.0465。為此，該礦山為創(chuàng)建綠色礦山，必須從企業(yè)內(nèi)、外兩方面改善。內(nèi)部系統(tǒng)：提高礦產(chǎn)資源的采選率、改善資源的開發(fā)條件，提高企業(yè)人員的素質(zhì)，發(fā)展循環(huán)經(jīng)濟(jì)，提高廢水處理率，內(nèi)部管理機(jī)制更為健全，加大宣傳，提高環(huán)保理念；外部方面：政府要加強(qiáng)礦山環(huán)境補(bǔ)償資金使用效益，完善礦山環(huán)境恢復(fù)保證金機(jī)制，把產(chǎn)學(xué)研得到的成果盡快地轉(zhuǎn)化為生產(chǎn)力，同時(shí)加大生態(tài)環(huán)境保護(hù)資金的投入。

表3 正負(fù)熵變值計(jì)算表

6 結(jié)論

礦山生態(tài)環(huán)境保護(hù)的耗散結(jié)構(gòu)的形成機(jī)制，有助于我們從新的視角認(rèn)識(shí)礦山生態(tài)環(huán)境保護(hù)。利用轉(zhuǎn)義后的Brusselator模型來判據(jù)礦山生態(tài)環(huán)境系統(tǒng)耗散結(jié)構(gòu)的形成是一種新的方法嘗試。通過計(jì)算影響礦山生態(tài)系統(tǒng)正負(fù)熵變，可以定量地知道在礦山生態(tài)環(huán)境保護(hù)上存在的問題和不足，從而采取針對(duì)性的對(duì)策來促進(jìn)礦山生態(tài)環(huán)境保護(hù)取得好的成效。總之，礦山生態(tài)環(huán)境保護(hù)的耗散結(jié)構(gòu)Brusselator模型，對(duì)礦山生態(tài)環(huán)境保護(hù)具有重要的指導(dǎo)意義。

［1］魏風(fēng)華，梁振戀，張樑.解決礦山生態(tài)環(huán)境保護(hù)問題的對(duì)策建議［J］.中國(guó)國(guó)土資源經(jīng)濟(jì)，2005（2）：25-35.

［2］顏永聰.石漠化礦區(qū)治理與礦山生態(tài)環(huán)境建設(shè)模式初探［J］.輕金屬，2009（4）：8-10.

［3］黃德林，汪琳.我國(guó)礦山生態(tài)環(huán)境保護(hù)法律制度完善問題研究［J］.中國(guó)人口·資源與環(huán)境，2008，18（5）：216-220.

［4］郁濱，等.系統(tǒng)工程理論［M］.合肥：中國(guó)科學(xué)技術(shù)大學(xué)出版社，2009.

［5］Takeo Maruyama，etc，Entropy-based assessment clustering of potential water resourses availability［J］.Journal of Hydrogy，2005（309）：104-113.

［6］Frank C，Krysiak.Entropy limits to growth，and the prospects for weak sustainability［J］.Ecological Economics，2006（58）：182-191.

［7］Rechberger H，Graedel T E.The contemporary European copper cycle：statistical entropy analysis［J］.Ecological Economics，2002（42）：59-72.

［8］Antoniou I，etc.On the efficient resources distribution in economics based on entropy［J］.Physicca A 2004（336）：549-562.

［9］Ernie Jowsey.A new basis for assessing the sustainability of natural resources［J］.Energy，2007（32）：906-911.

［10］張鐵男，程寶元，張亞娟.基于耗散結(jié)構(gòu)的企業(yè)管理熵Brusselator模型研究［J］.管理工程學(xué)報(bào)，2010，24（3）：103-108.

［11］魏遙.基于Brusselator模型的產(chǎn)融集團(tuán)生成機(jī)制研究［J］.管理評(píng)論，2010（8）：39-44.

［12］魏遙，雷良海.基于Brusselator模型的產(chǎn)融集團(tuán)演化機(jī)制研究［J］.科技進(jìn)步與對(duì)策，2009，26（18）：60-64.

［13］賀建勛.系統(tǒng)建模與數(shù)學(xué)模型［M］.福州：福建科學(xué)技術(shù)出版社，1995.