基于綠色供應鏈視角的企業(yè)生態(tài)創(chuàng)新擴散演化博弈研究

楊國忠，周午陽

（中南大學商學院，湖南長沙 410083）

如今人類社會經(jīng)濟蓬勃發(fā)展，與之形成鮮明對比的即是生態(tài)環(huán)境的日益惡化，過度捕撈、亂砍濫伐時有發(fā)生，空氣污染、垃圾成山、物種滅絕的現(xiàn)象愈演愈烈，環(huán)境問題已然成為全球性的經(jīng)濟發(fā)展議題［1］。在低碳背景下，我國積極推動全球氣候治理多邊進程，堅持創(chuàng)新、協(xié)調(diào)、綠色、開放、共享的發(fā)展理念，實現(xiàn)社會經(jīng)濟可持續(xù)發(fā)展［2］。國務院在2021 年10 月24 日發(fā)布的《關于完整準確全面貫徹新發(fā)展理念做好碳達峰碳中和工作的意見》中明確指出要構建綠色低碳循環(huán)發(fā)展經(jīng)濟體系，提升能源利用效率。隨著我國經(jīng)濟發(fā)展進入新常態(tài)，生態(tài)創(chuàng)新已初見成效，但仍然面臨營商環(huán)境和政策機制不完善、研發(fā)成本高等障礙［3-4］，商業(yè)化擴散存在較大阻力。因此，研究企業(yè)生態(tài)創(chuàng)新擴散，既是實現(xiàn)我國生態(tài)文明建設的需要，也是實現(xiàn)企業(yè)可持續(xù)成長的需要。

1 文獻綜述

關于技術擴散的研究最早可追溯至20 世紀，Rogers［5］認為技術擴散是技術隨著時間的推移通過各種渠道傳播到社會系統(tǒng)的過程，擴散過程分為了解、興趣、評估、試驗和采納五個階段。Bass［6］提出了經(jīng)典的“S”形增長曲線的Bass 模型。這些傳統(tǒng)的技術擴散理論為生態(tài)創(chuàng)新擴散的研究提供了理論基礎。目前國內(nèi)外學者對生態(tài)創(chuàng)新擴散的研究頗有成果，針對生態(tài)創(chuàng)新擴散的定義，諸多學者各抒己見，本文更認同田紅娜等［7］的觀點，生態(tài)創(chuàng)新擴散是指在一定環(huán)境條件下，創(chuàng)新擴散源提供的生態(tài)創(chuàng)新成果通過中介渠道最終被創(chuàng)新接受體接收采納的過程。

現(xiàn)階段，國內(nèi)外學者對生態(tài)創(chuàng)新擴散的研究大都集中在兩方面，一是對擴散過程模型的研究，例如襲希等［8］構建低碳技術擴散的逾滲模型并做了仿真實驗。李真等［9］基于消費者互動網(wǎng)絡采用計算實驗方法構建模型，模擬相關參數(shù)變化對綠色產(chǎn)品擴散的影響。Kok 等［10］追蹤了美國48 個大都市地區(qū)綠色建筑的擴散情況，探討綠色建筑在空間和時間上擴散的速率以及決定因素。Wang 等［11］基于Bass模型構建出新能源產(chǎn)業(yè)的領先技術擴散模型。二是對擴散影響因素的研究，學者重點聚焦于政府管制對生態(tài)技術擴散的影響，例如，張倩等［12］分析了3種環(huán)境規(guī)制工具（排污稅、排污許可證和排放標準）和不同環(huán)境規(guī)制強度對企業(yè)綠色技術擴散的影響，為政府制定環(huán)境規(guī)制政策提供了新思路。宋英杰等［13］從橫向縱向兩方面分析了環(huán)境分權對環(huán)保技術擴散的影響。張躍勝［14］從地區(qū)差異角度，發(fā)現(xiàn)技術存量、技術吸收能力、碳減排技術補貼和碳排放稅均會對碳減排技術擴散產(chǎn)生影響。Cantonoa 等［15］結(jié)合消費者的傳染特征和主體特征建立了一個新技術擴散的網(wǎng)絡模型，分析了補貼政策和學習曲線效應對新生態(tài)技術擴散的影響。Foxon 等［16］定性分析了環(huán)境政策對國家可持續(xù)創(chuàng)新擴散的影響，研究提出應當結(jié)合創(chuàng)新現(xiàn)狀制定政策，并應該運用系統(tǒng)思維。Halkos 等［17］從國家層面分析外商投資、國民教育水平、體制效率和管理質(zhì)量對環(huán)保技術擴散的影響。

總之，學者們對生態(tài)創(chuàng)新擴散的研究較為翔實，但角度較為單一且基本是從宏觀視角來分析。然而，生態(tài)創(chuàng)新擴散涉及多個主體且擴散環(huán)境復雜多變，其內(nèi)在機制仍有待探索，于是，有學者引入博弈論來深入研究微觀層面的生態(tài)創(chuàng)新擴散機制。肖漢杰等［18］考慮到技術協(xié)同效應，引入環(huán)境友好技術互補系數(shù)和風險損益，構建低碳環(huán)境友好技術創(chuàng)新擴散的演化博弈模型。于麗靜等［19］和曹霞等［20］構建了物流企業(yè)、政府與消費者三方博弈模型，并用數(shù)值仿真分析不同參數(shù)變化對系統(tǒng)演化路徑的影響。羅明等［21］構建政府和企業(yè)的協(xié)同環(huán)境治理博弈模型，分析了環(huán)境稅率、系統(tǒng)突變率和主體學習能力對系統(tǒng)演化結(jié)果的影響。Mohsin 等［22］構建了制造商和零售商之間的完全信息微分博弈模型，比較集中決策和分散決策的均衡解并且設計出動態(tài)批發(fā)價格機制的新合作協(xié)調(diào)模型。Galasso 等［23］建立了一個綠色技術采納博弈模型，分析了在消費者補貼、企業(yè)生產(chǎn)成本不對稱等不同情形下的產(chǎn)品市場競爭博弈的結(jié)果。Carfi 等［24］將博弈論應用于低碳技術在國際合作和推廣方面，在雙贏的情況下尋找帕累托最優(yōu)解。Wang 等［25］構建了政府、創(chuàng)新供給企業(yè)和潛在需求企業(yè)三方博弈模型。Zhang 等［26］構建復雜網(wǎng)絡環(huán)境下企業(yè)綠色技術擴散的演化博弈模型，模擬分析補貼、環(huán)境稅和碳交易市場等市場導向政策對擴散的影響。

梳理文獻發(fā)現(xiàn)，國內(nèi)外學者對生態(tài)創(chuàng)新擴散的模型、影響因素和博弈相關方面研究較為充分。但現(xiàn)有研究過多關注生態(tài)創(chuàng)新擴散的宏觀層面，在微觀層面也僅僅關注了企業(yè)與利益相關方之間的聯(lián)系，忽視了綠色供應鏈層面企業(yè)的相互影響，尤其是供應鏈上下游企業(yè)間的演化博弈機理?；诖?，本文基于綠色供應鏈的視角，采用演化博弈方法構建供應商與生產(chǎn)商的博弈模型，分析綠色供應鏈上下游企業(yè)生態(tài)創(chuàng)新擴散的演化穩(wěn)定性，從而提出政策建議。

2 基本假設與模型構建

生態(tài)創(chuàng)新包括管理創(chuàng)新、產(chǎn)品創(chuàng)新、末端處理技術創(chuàng)新和清潔流程工藝創(chuàng)新［27］。管理創(chuàng)新是指重新定義生產(chǎn)和運營流程以提高內(nèi)部效率；產(chǎn)品創(chuàng)新是設計出更低碳環(huán)保的產(chǎn)品；末端處理技術創(chuàng)新注重回收處理在生產(chǎn)過程中產(chǎn)生的廢棄材料和副產(chǎn)品；清潔流程工藝創(chuàng)新是在生產(chǎn)流程的各個環(huán)節(jié)減少污染和浪費、提高能源利用效率。在供應鏈中加入“綠色”一詞意味著從與自然環(huán)境結(jié)合的角度來看待供應鏈，綠色供應鏈層面的生態(tài)創(chuàng)新擴散主要涉及管理創(chuàng)新和清潔流程工藝創(chuàng)新兩方面。同時，為了發(fā)展綠色經(jīng)濟，政府也將采取必要的宏觀調(diào)控措施來引導企業(yè)生態(tài)創(chuàng)新擴散。

本文考慮的供應鏈模型是在政府參與模式下的由上游多個供應商和下游多個生產(chǎn)商構成的兩層供應鏈結(jié)構，在供應鏈視角下的生態(tài)創(chuàng)新擴散存在兩種機制，其一是供應商為創(chuàng)新提供方，生產(chǎn)商是創(chuàng)新采納方；其二是生產(chǎn)商為創(chuàng)新提供方，供應商是創(chuàng)新采納方。若是下游生產(chǎn)商研發(fā)出新的清潔工藝技術或者掌握著先進的管理創(chuàng)新知識技能，出于獲取垂直分工利益的考慮，下游生產(chǎn)商將自身的生態(tài)創(chuàng)新出售給上游供應商以保證綠色原料的供應；而上游生產(chǎn)商因為引入生態(tài)創(chuàng)新需付高額的費用，其本身并不愿意采用，在政府管制下，必須使其采用生態(tài)創(chuàng)新后的收益不小于采用之前方可實現(xiàn)擴散成功。當上游供應商研發(fā)出新的清潔工藝技術或者掌握著先進的管理創(chuàng)新知識技能，若將創(chuàng)新提供給下游企業(yè)，可使其具備應用上游綠色產(chǎn)品系統(tǒng)的能力從而進一步開拓市場，但下游企業(yè)獲得該技術提升自身能源利用效率后或許會減少對上游原料的需求量。

綜上所述，供應鏈上下游企業(yè)之間存在相互制約的影響。供應商和生產(chǎn)商各自無法獲得對方全部信息，只能對有限的信息作出有限理性決策，經(jīng)過多次重復且動態(tài)的博弈，最終達到穩(wěn)定的策略。根據(jù)李昆［28］和方竺乾等［29］的研究框架，在此基礎上運用演化博弈方法并加入供應鏈的特征加以改進，本文提出如下假設：

假設1：若供應商和生產(chǎn)商均不進行生態(tài)創(chuàng)新，則雙方的收益分別為若供應商和生產(chǎn)商進行生態(tài)創(chuàng)新，則雙方的研發(fā)成本分別為供應商進行創(chuàng)新擴散（采納）的概率為，生產(chǎn)商進行創(chuàng)新采納（擴散）的概率為，且為時間t的函數(shù)。

假設2：若是供應商研發(fā)出清潔工藝創(chuàng)新技術或掌握先進的管理創(chuàng)新知識技能，保障了原材料的綠色度，相應提高了供應鏈終端產(chǎn)品的綠色度，消費者需求也隨之增加，供應商增加的經(jīng)濟收益計為。生產(chǎn)商若想引入供應商的技術和知識需要支付專利費，生產(chǎn)商采納后增加的經(jīng)濟收益為，并且生產(chǎn)商通過采納生態(tài)創(chuàng)新提升了自身能源利用效率后，會相應減少上游原料的需求量，這使供應商遭受經(jīng)濟損失；而生產(chǎn)商若不采納新的技術和知識，也將分享供應商碳減排帶來的收益即“搭便車”行為，生產(chǎn)商的搭便車收益為，且。

假設3：若是生產(chǎn)商研發(fā)出清潔工藝創(chuàng)新技術或掌握先進的管理創(chuàng)新知識技能，減少了加工環(huán)節(jié)碳排放進而提高終端產(chǎn)品綠色度，消費者需求隨之增加，生產(chǎn)商增加的經(jīng)濟收益計為，同時生產(chǎn)商提升自身能源利用效率會相應減少上游原料的需求，這使供應商遭受經(jīng)濟損失。供應商若想引入新的技術和知識需要支付專利費，采納后供應商增加的經(jīng)濟收益為；供應商若不采納新的技術和知識，也可獲得生產(chǎn)商碳減排帶來的共享收益即“搭便車”行為，供應商的搭便車收益為，且。

假設4：若生態(tài)技術得到擴散即生產(chǎn)商和供應商都具有新的生態(tài)技術，這對整條供應鏈的環(huán)保建設有所貢獻，供應鏈的綠色度大幅提升、整體競爭力顯著增強，創(chuàng)造出“1+1 ＞2”的協(xié)同效應［30］，供應商和生產(chǎn)商獲得的協(xié)同效益分別為

假設5：政府為了發(fā)展低碳經(jīng)濟鼓勵企業(yè)進行生態(tài)創(chuàng)新擴散，對積極擴散創(chuàng)新成果的企業(yè)按其研發(fā)投入給予各種優(yōu)惠補貼，補貼系數(shù)（i=1 表示供應商，i=2 表示生產(chǎn)商；，對不采納生態(tài)技術的企業(yè)給予一定的懲罰（i=1 表示供應商，i=2 表示生產(chǎn)商），同時對采取創(chuàng)新保護主義的企業(yè)推行一定的限制政策使其受到損失（i=1 表示供應商，i=2 表示生產(chǎn)商）。

根據(jù)以上假設構建的支付矩陣如表1和表2所示。

表1 供應商是技術擴散方而生產(chǎn)商是采納決策方

表2 生產(chǎn)商是技術擴散方而供應商是采納決策方

3 模型求解及分析

3.1 供應商為擴散方的博弈模型

供應商采取生態(tài)技術保護的收益為：

供應商的平均收益為：

生產(chǎn)商采納生態(tài)技術的收益為：

生產(chǎn)商不采納生態(tài)技術的收益為：

生產(chǎn)商的平均收益為：

根據(jù)Malthusian 動態(tài)方程原理，得到供應商和生產(chǎn)商生態(tài)技術擴散的復制者動態(tài)方程分別為：

應用雅克比矩陣的局部穩(wěn)定性來進行系統(tǒng)均衡點的穩(wěn)定性分析。該博弈系統(tǒng)的雅克比矩陣如下所示：

當平衡點使得detJ（J行列式的值）＞0 且trJ（J行列式的跡）＜0 時，平衡點處于局部穩(wěn)定狀態(tài)。以此為判定依據(jù)得出穩(wěn)定點及其所對應系統(tǒng)演化狀態(tài)的推論，具體分析見表3～表6 所示。

表3 推論1 情況下的穩(wěn)定性分析

表4 推論2 情況下的穩(wěn)定性分析

表5 推論3 情況下的穩(wěn)定性分析

表6 推論4 情況下的穩(wěn)定性分析

圖1 供應商為擴散方的系統(tǒng)演化相位圖

3.2 生產(chǎn)商為擴散方的博弈模型

供應商采納生態(tài)技術的收益為：

供應商不采納生態(tài)技術的收益為：

供應商的平均收益為：

生產(chǎn)商采取生態(tài)技術擴散的收益：

生產(chǎn)商采取生態(tài)技術保護的收益：

生產(chǎn)商采取生態(tài)技術保護的收益：

根據(jù)Malthusian 動態(tài)方程原理，得到供應商和生產(chǎn)商生態(tài)技術擴散的復制者動態(tài)方程分別為：

當平衡點使得detJ（J行列式的值）＞0 且trJ（J行列式的跡）＜0 時，平衡點處于局部穩(wěn)定狀態(tài)。以此為判定依據(jù)得出穩(wěn)定點及其所對應系統(tǒng)演化狀態(tài)的推論，具體分析見表7 和表8 所示。

表7 推論5 情況下的穩(wěn)定性分析

表8 推論6 情況下的穩(wěn)定性分析

圖2 生產(chǎn)商為擴散方的系統(tǒng)演化相位圖

4 數(shù)值分析

下面運用MATLAB2017b 軟件進行仿真，結(jié)合現(xiàn)實狀況分別對供應商為擴散方、生產(chǎn)商為采納方的系統(tǒng)a 和生產(chǎn)商為擴散方、供應商為采納方的系統(tǒng)b 進行數(shù)值分析，探求各參數(shù)變化對各系統(tǒng)演化的影響。假設系統(tǒng)演化的初始點

4.1 專利費對系統(tǒng)演化的影響

模型參數(shù)設為I1=18，I2=18，F(xiàn)1=25，F(xiàn)2=25，W1=20，W1=20，β1=5，β2=5，α1=0.15，α2=0.15，K1=2，K2=2，C1=2，C2=2，g1=15，g2=15。P1和P2分別取4、8、12、16、20，供應商為擴散方的系統(tǒng)a 演化情況和生產(chǎn)商為擴散方的系統(tǒng)b 演化情況分別如圖3（a）和圖3（b）所示。

圖3 專利費變化對系統(tǒng)演化的影響

觀察圖3（a）可以看出，當專利費為4 和8 時，系統(tǒng)a 演化至供應商保護創(chuàng)新而生產(chǎn)商采納創(chuàng)新的策略狀態(tài)，此時生態(tài)創(chuàng)新未能擴散；當專利費提高到12 和16 時，系統(tǒng)a 演化至供應商擴散創(chuàng)新且生產(chǎn)商采納創(chuàng)新的穩(wěn)定狀態(tài)，生態(tài)創(chuàng)新得以成功擴散；而當專利費提高到20 時，系統(tǒng)a 演化至供應商擴散創(chuàng)新而生產(chǎn)商不采納創(chuàng)新的穩(wěn)定狀態(tài)，生態(tài)創(chuàng)新擴散趨向失敗。觀察圖3（b）可以看出，當專利費為4、8、12、16 時，系統(tǒng)b 演化至生產(chǎn)商擴散創(chuàng)新且供應商采納創(chuàng)新的策略狀態(tài)，生態(tài)創(chuàng)新成功擴散；當專利費提高到20 時，系統(tǒng)b 演化至生產(chǎn)商擴散創(chuàng)新而供應商不采納創(chuàng)新的穩(wěn)定狀態(tài)，此時生態(tài)創(chuàng)新未能擴散。

觀察專利費對系統(tǒng)a 和系統(tǒng)b 中的采納方的決策行為影響可知，隨著費用的提高，采納方的策略逐漸從采納生態(tài)創(chuàng)新轉(zhuǎn)移至不采納生態(tài)創(chuàng)新。觀察擴散方?jīng)Q策行為影響可知，對于系統(tǒng)a 而言，專利費提高到一定閾值時擴散方才會選擇擴散生態(tài)創(chuàng)新；但是對于系統(tǒng)b 而言，不同的費用對擴散方的演化路徑影響不大，最終都會演化到擴散生態(tài)創(chuàng)新的策略狀態(tài)，這是因為供應商擴散生態(tài)創(chuàng)新會造成其一定程度的未來間接損失，所以需要一定的專利收入來彌補。

4.2 協(xié)同效益對系統(tǒng)演化的影響

模型參數(shù)設為P1=16，P2=12，I1=18，I2=18，F(xiàn)1=25，F(xiàn)2=25，W1=20，W2=20，α1=0.15，α2=0.15，K1=2，K2=2，C1=2，C2=2，g1=15，g2=15。β1和β2分別取0、1、3、5、8 時，供應商為擴散方的系統(tǒng)a演化情況和生產(chǎn)商為擴散方的系統(tǒng)b 演化情況分別如圖4（a）和圖4（b）所示。

圖4 協(xié)同效益變化對系統(tǒng)演化的影響

由圖4（a）可知，當協(xié)同效益為0、1、3 時，系統(tǒng)a 演化至供應商保護創(chuàng)新而生產(chǎn)商采納創(chuàng)新的策略狀態(tài)，生態(tài)創(chuàng)新未能擴散。當協(xié)同效益增長至5 和8 時，系統(tǒng)a 演化至供應商擴散創(chuàng)新且生產(chǎn)商采納創(chuàng)新的穩(wěn)定狀態(tài)，生態(tài)創(chuàng)新成功擴散。由圖4（b）可知，當協(xié)同效益為0、1、3 時，系統(tǒng)b 演化至生產(chǎn)商擴散創(chuàng)新而供應商不采納創(chuàng)新的策略狀態(tài)，生態(tài)創(chuàng)新未能擴散。當協(xié)同效益逐漸增加至5 和8 時，系統(tǒng)b 演化至生產(chǎn)商擴散創(chuàng)新且供應商生態(tài)創(chuàng)新的穩(wěn)定狀態(tài)，此時生態(tài)創(chuàng)新得以成功擴散。

觀察協(xié)同效益對系統(tǒng)a 和系統(tǒng)b 中的擴散方和采納方的決策行為影響可知，不論是系統(tǒng)a 還是系統(tǒng)b，隨著協(xié)同效益的提高，上游供應商的策略逐漸從保護（不采納）生態(tài)創(chuàng)新轉(zhuǎn)移至擴散（采納）生態(tài)創(chuàng)新；而不同的協(xié)同效益對下游生產(chǎn)商的演化路徑影響不大，最終都會演化到采納（擴散）生態(tài)創(chuàng)新的策略狀態(tài)。

4.3 “搭便車”收益對系統(tǒng)演化的影響

模型參數(shù)設為P1=16，P2=12，I1=18，I2=18，F(xiàn)1=25，F(xiàn)2=25，W1=20，W2=20，α1=0.15，α2=0.15，β1=5，β2=5，K1=2，K2=2，C1=2，C2=2。g1和g1分別取7、11、15、19、23 時，供應商為擴散方的系統(tǒng)a 演化情況和生產(chǎn)商為擴散方的系統(tǒng)b 演化情況分別如圖5（a）和圖5（b）所示。

圖5 “搭便車”收益變化對系統(tǒng)演化的影響

由圖5（a）可知，當“搭便車”收益低于19 時，系統(tǒng)a 演化至供應商擴散創(chuàng)新且生產(chǎn)商采納創(chuàng)新的策略狀態(tài)，此時生態(tài)創(chuàng)新成功擴散；當“搭便車”收益增長至23 時，系統(tǒng)a 演化至供應商擴散創(chuàng)新而生產(chǎn)商不采納創(chuàng)新的狀態(tài)，生態(tài)創(chuàng)新未能擴散。由圖5（b）可知，當協(xié)同效益為7、11、15 時，系統(tǒng)b 演化至生產(chǎn)商擴散創(chuàng)新且供應商采納創(chuàng)新的策略狀態(tài)，生態(tài)創(chuàng)新成功擴散；當“搭便車”收益逐漸增加至19 和23 時，系統(tǒng)b 演化至生產(chǎn)商擴散創(chuàng)新而供應商不采納創(chuàng)新的狀態(tài)，此時生態(tài)創(chuàng)新擴散失敗。

觀察“搭便車”收益對系統(tǒng)a 和系統(tǒng)b 中的采納方的決策行為影響可知，“搭便車”收益的增長會降低采納方采納生態(tài)創(chuàng)新的積極性。

4.4 政府補貼系數(shù)對系統(tǒng)演化的影響

模型 參數(shù)設為P1=16，P2=12，I1=18，I2=18，F(xiàn)1=25，F(xiàn)2=25，W1=20，W2=20，β1=5，β2=5，K1=2，K2=2，C1=2，C2=2，g1=15，g2=15。α1分別取0、0.05、0.1、0.15、0.3，供應商為擴散方的系統(tǒng)a 演化情況如圖6（a）所示，α2分別取0、0.15、0.3、0.5、0.7，生產(chǎn)商為擴散方的系統(tǒng)b 演化情況如圖6（b）所示。

圖6 政府補貼變化對系統(tǒng)演化的影響

由圖6（a）可知，當α1取0 和0.05 時，系統(tǒng)a演化至供應商保護創(chuàng)新而生產(chǎn)商采納創(chuàng)新的策略，生態(tài)創(chuàng)新擴散沒有成功。當α1取0.1 及以上時，系統(tǒng)a 演化至供應商擴散創(chuàng)新且生產(chǎn)商采納創(chuàng)新使得生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)。由圖6（b）可知，不管α2取何值，系統(tǒng)b 均演化至生產(chǎn)商擴散創(chuàng)新且供應商采納創(chuàng)新使得生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)，但隨著補貼系數(shù)的增加，生產(chǎn)商向擴散生態(tài)創(chuàng)新策略的方向演化速度略有加快。

總之，政府補貼對系統(tǒng)向生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)演化有明顯的促進作用。但對于系統(tǒng)a 而言，只有政府補貼達到一定的閾值后才能推動生態(tài)創(chuàng)新的成功擴散。

4.5 政府對擴散方懲罰力度對系統(tǒng)演化的影響

模型參數(shù)設為P1=16，P2=12，I1=18，I2=18，F(xiàn)1=25，F(xiàn)2=25，W1=20，W2=20，α1=0.15，α2=0.15，β1=5，β1=5，K1=2，K2=2，g1=15，g2=15。C1分別取0、1、2、3、4 時供應商為擴散方的系統(tǒng)a 演化情況如圖7（a）所示，C2分別取0、2、4、6、8 時生產(chǎn)商為擴散方的系統(tǒng)b 演化情況如圖7（b）所示。

圖7 懲罰擴散方的力度變化對系統(tǒng)演化的影響

由圖7（a）可知，當C1取0 時，系統(tǒng)a 演化至供應商保護創(chuàng)新而生產(chǎn)商采納創(chuàng)新的策略狀態(tài)，此時生態(tài)創(chuàng)新未能擴散，而C1加大到1 以上時，供應商因忌憚高罰金逐漸朝著擴散創(chuàng)新演化，系統(tǒng)a 演化生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)。由圖7（b）可知，不管C2取何值，系統(tǒng)b 均演化至生產(chǎn)商擴散創(chuàng)新且供應商采納創(chuàng)新使得生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)，但隨著補貼系數(shù)的增加，生產(chǎn)商向擴散生態(tài)創(chuàng)新策略的方向演化速度略有加快。

以上說明對擴散方的懲罰機制對生態(tài)創(chuàng)新的成功擴散有一定的推動作用。隨著懲罰力度的提高，系統(tǒng)a 擴散方的策略逐漸從保護生態(tài)創(chuàng)新轉(zhuǎn)移至擴散生態(tài)創(chuàng)新，系統(tǒng)a 也逐漸演化到生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)。

4.6 政府對采納方懲罰力度對系統(tǒng)演化的影響

模型 參數(shù)設為P1=16，P2=16，I1=18，I2=18，F(xiàn)1=25，F(xiàn)2=25，W1=20，W2=20，α1=0.15，α2=0.15，β1=5，β2=5，C1=2，C2=2，g1=15，g2=15。K1和K2分別取0、1、2、3、4 時，供應商為擴散方的系統(tǒng)a演化情況和生產(chǎn)商為擴散方的系統(tǒng)b 演化情況分別如圖8（a）和圖8（b）所示。

圖8 懲罰采納方的力度變化對系統(tǒng)演化的影響

由圖8（a）可知，當K2取0 時，系統(tǒng)a 演化至供應商擴散生態(tài)創(chuàng)新而生產(chǎn)商不采納生態(tài)創(chuàng)新的策略狀態(tài)，此時生態(tài)創(chuàng)新未能擴散，而K2逐漸加大到2 時，生產(chǎn)商因忌憚高罰金逐漸朝著采納生態(tài)創(chuàng)新演化，系統(tǒng)a 演化生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)。由圖8（b）可知，當K1取0 時，系統(tǒng)b 演化至生產(chǎn)商擴散創(chuàng)新而供應商不采納創(chuàng)新的策略，生態(tài)創(chuàng)新擴散沒有成功，而K1逐漸增加到2 以上時，供應商因忌憚高罰金逐漸朝著采納生態(tài)創(chuàng)新演化，系統(tǒng)b演化至生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)。

總的來說，政府對采納方的懲罰機制對系統(tǒng)向生態(tài)創(chuàng)新成功擴散的穩(wěn)定狀態(tài)演化有顯著的推動作用。相較于前兩種規(guī)制手段，政府對采納方懲罰力度效果更為顯著。

5 結(jié)論與建議

本文基于綠色供應鏈視角，利用演化博弈方法構建了供應商和生產(chǎn)商的兩種博弈模型，分析生態(tài)創(chuàng)新擴散的穩(wěn)定性問題，并運用MATLAB2017b 軟件仿真研究不同參數(shù)變化對生態(tài)創(chuàng)新擴散演化軌跡的影響。結(jié)果表明：（1）綠色供應鏈的生態(tài)創(chuàng)新擴散是適合規(guī)?；a(chǎn)的合作模式，即研發(fā)出生態(tài)核心技術或者掌握先進的管理創(chuàng)新知識技能的企業(yè)可以采用外包戰(zhàn)略，與上下游企業(yè)分享技術或知識從而獲取更多市場份額及企業(yè)影響力，獲得戰(zhàn)略高回報；而處于創(chuàng)新采納方的企業(yè)往往選擇引入生態(tài)創(chuàng)新，提高資源利用效率的同時也獲取更高收益。（2）生產(chǎn)商研發(fā)出生態(tài)創(chuàng)新更有利于該創(chuàng)新的擴散；供應商擴散自身研發(fā)的生態(tài)創(chuàng)新會有未來的間接損失，導致其需要一定的專利費或政府補貼等補償來彌補這些損失。（3）協(xié)同效益和專利費對生態(tài)創(chuàng)新擴散有一定的促進作用，尤其是當上游供應商企業(yè)選擇提供或采納生態(tài)創(chuàng)新存在正凈收益時，生態(tài)創(chuàng)新才能成功擴散；此外，“搭便車”收益會降低采納方選擇采納生態(tài)創(chuàng)新的積極性。（4）政府的獎罰機制對生態(tài)創(chuàng)新擴散的影響不同。相較于政府對擴散方的獎懲機制，政府提高對采納方懲罰力度更能有效促進供應鏈層面的生態(tài)創(chuàng)新擴散；政府對擴散方的獎懲機制中，對供應商作為擴散方時的作用較為顯著，而對生產(chǎn)商作為擴散方的影響有限。

基于以上研究及結(jié)論，為構建更有效的綠色供應鏈企業(yè)生態(tài)創(chuàng)新擴散機制，本文給出以下建議：（1）目前我國生態(tài)創(chuàng)新發(fā)展處于起步階段，研發(fā)成本較高、創(chuàng)新動力不足，政府應以財稅補貼為主、輔之適度的懲罰措施，并搭建產(chǎn)學合作的創(chuàng)新研發(fā)應用平臺，鼓勵科研機構投入到生態(tài)創(chuàng)新的研究，采用生態(tài)創(chuàng)新研發(fā)到產(chǎn)業(yè)化的一體化模式，降低企業(yè)研發(fā)成本的同時也促進生態(tài)創(chuàng)新的商業(yè)化應用。（2）政府制定補貼監(jiān)管政策要有針對性，根據(jù)上文博弈分析，上游供應商擴散生態(tài)創(chuàng)新會給自身帶來一定損失，政府應及時調(diào)整財政補貼、稅收優(yōu)惠、產(chǎn)業(yè)扶持等政策，加大對供應商的扶持；同時對創(chuàng)新采納方輔以一定監(jiān)管力度，完善各項管理措施的協(xié)調(diào)配合，充分調(diào)動采納方的積極性。（3）政府應建立起綠色化和市場化有機結(jié)合的市場經(jīng)濟制度，并且完善生態(tài)知識產(chǎn)權的保護機制，利用網(wǎng)絡媒體宣傳綠色創(chuàng)新的效益，培養(yǎng)公眾綠色消費習慣，提高市場需求，使企業(yè)能夠從生態(tài)創(chuàng)新的研發(fā)中獲得更多收益，激發(fā)企業(yè)生態(tài)創(chuàng)新的內(nèi)在動力并激勵企業(yè)自愿擴散技術和知識。（4）供應鏈中的上下游企業(yè)應從整條供應鏈的綠色可持續(xù)發(fā)展角度出發(fā)，建設供應鏈誠信機制，共同打擊“搭便車”行為，以綠色、生態(tài)、環(huán)保為基本理念，優(yōu)化生產(chǎn)環(huán)節(jié)，最大限度利用資源，提升供應鏈綠色度和競爭力從而實現(xiàn)更高利潤。