趙愛武，杜建國，關(guān)洪軍

1 江蘇大學(xué) 管理學(xué)院，江蘇 鎮(zhèn)江 212013 2 山東財經(jīng)大學(xué) 管理科學(xué)與工程學(xué)院，濟南 250014

環(huán)境稅情景下企業(yè)環(huán)境技術(shù)創(chuàng)新模擬

趙愛武1，杜建國1，關(guān)洪軍2

1 江蘇大學(xué) 管理學(xué)院，江蘇 鎮(zhèn)江 212013 2 山東財經(jīng)大學(xué) 管理科學(xué)與工程學(xué)院，濟南 250014

環(huán)境技術(shù)創(chuàng)新是解決經(jīng)濟發(fā)展和環(huán)境可持續(xù)兩難問題的重要途徑，為實現(xiàn)環(huán)境績效和經(jīng)濟績效的雙贏，環(huán)境稅應(yīng)以誘導(dǎo)企業(yè)環(huán)境技術(shù)創(chuàng)新為重要目的。

關(guān)于環(huán)境稅對企業(yè)環(huán)境技術(shù)創(chuàng)新的激勵作用，學(xué)者們基于不同的情景假設(shè)、研究樣本、分析方法和實驗?zāi)Ｐ?，取得了并不一致的結(jié)論。一方面，環(huán)境稅可以激勵企業(yè)實施環(huán)境技術(shù)創(chuàng)新，并通過創(chuàng)新補償機制，使企業(yè)實現(xiàn)經(jīng)濟績效和環(huán)境績效的雙贏；另一方面，環(huán)境稅會帶來企業(yè)生產(chǎn)成本的增加，削弱企業(yè)的創(chuàng)新能力?？梢?，環(huán)境稅對企業(yè)環(huán)境技術(shù)創(chuàng)新行為的影響是一個動態(tài)的、演化的過程，單純依靠實證或數(shù)理方法難以全面認識環(huán)境稅對企業(yè)環(huán)境技術(shù)創(chuàng)新行為的影響機理。

基于社會科學(xué)計算實驗方法，將環(huán)境稅情景下企業(yè)環(huán)境技術(shù)創(chuàng)新行為的研究置于復(fù)雜多變的動態(tài)市場環(huán)境，在識別影響企業(yè)環(huán)境技術(shù)創(chuàng)新行為關(guān)鍵因素及其傳導(dǎo)機制的基礎(chǔ)上，運用數(shù)理方法將微觀個體的特征、行為和規(guī)則描述為個體與環(huán)境的函數(shù)，進而構(gòu)建企業(yè)環(huán)境技術(shù)創(chuàng)新行為動態(tài)仿真模型，運用計算機系統(tǒng)模擬不同消費者選擇、產(chǎn)品競爭和環(huán)境稅情景下企業(yè)環(huán)境技術(shù)創(chuàng)新行為演化過程，探索環(huán)境稅征收時機及強度對企業(yè)環(huán)境績效和經(jīng)濟績效的影響。實驗結(jié)果表明，在市場機制、產(chǎn)品競爭和環(huán)境稅的協(xié)同作用下，當(dāng)可替代環(huán)境創(chuàng)新技術(shù)達到一定水平時，較高的環(huán)境稅能夠推動企業(yè)采納環(huán)境創(chuàng)新技術(shù)，并實現(xiàn)減排和增效的目標，促使產(chǎn)業(yè)向綠色經(jīng)濟方向發(fā)展；同時環(huán)境稅政策的實施時機對于產(chǎn)業(yè)向環(huán)境創(chuàng)新技術(shù)方向演化起關(guān)鍵作用，這與環(huán)境創(chuàng)新技術(shù)的發(fā)展水平密切相關(guān)。反之，較低的環(huán)境稅情景下，可能出現(xiàn)企業(yè)花錢買污染的狀況，與環(huán)境稅初衷背道而馳。

本模型同時可以作為環(huán)境稅設(shè)計的參考工具，通過改變模擬參數(shù)，觀察企業(yè)環(huán)境技術(shù)創(chuàng)新行為及其績效的演化。

環(huán)境技術(shù)創(chuàng)新；環(huán)境稅；計算實驗；經(jīng)濟績效；環(huán)境績效

1 引言

改革開放以來，中國經(jīng)濟發(fā)展取得了舉世矚目的成就，30年間年均GDP增長率高達9.8%。然而，在經(jīng)濟水平高速發(fā)展的背后，是粗放型傳統(tǒng)經(jīng)濟模式導(dǎo)致的資源匱乏和環(huán)境惡化，中國目前面臨的環(huán)境問題是發(fā)達國家在上百年工業(yè)化過程中分階段出現(xiàn)的[1]。這一現(xiàn)象必將威脅中國社會的和諧發(fā)展，成為經(jīng)濟可持續(xù)增長的瓶頸，并有損中國負責(zé)任的大國形象。發(fā)達國家尤其是經(jīng)濟合作與發(fā)展組織國家的實踐表明，環(huán)境稅能夠?qū)?jié)能減排、保護環(huán)境起到很好的促進作用。它在生產(chǎn)成本和市場價格的核算中考慮了環(huán)境污染和生態(tài)破壞導(dǎo)致的社會成本，并通過市場機制優(yōu)化環(huán)境資源配置，較之傳統(tǒng)的行政手段具有較高的經(jīng)濟效率，使企業(yè)能夠以最經(jīng)濟的方式對市場信號做出反應(yīng)[2]。中國仍然是發(fā)展中國家，城市化和工業(yè)化目標尚未實現(xiàn)，仍需大力發(fā)展經(jīng)濟以提高社會福利水平，解決大量貧困人口問題。因此，兌現(xiàn)全球減排義務(wù)的同時，必須考慮經(jīng)濟的增長問題，如何在發(fā)展經(jīng)濟的同時實現(xiàn)減排目標成為亟待解決的熱點問題。

2 相關(guān)研究評述

關(guān)于環(huán)境規(guī)制下企業(yè)環(huán)境績效和經(jīng)濟績效能否實現(xiàn)雙贏的問題，學(xué)者們基于不同的情景假設(shè)、研究樣本、分析方法和實驗?zāi)Ｐ?，取得了并不一致的結(jié)論。一方面，環(huán)境規(guī)制通過強制企業(yè)購置治污設(shè)備和技術(shù)或限制污染密集型產(chǎn)品的生產(chǎn)、減少采用特定投入要素組合等多種方式達到所要求的環(huán)境標準，會帶來企業(yè)生產(chǎn)成本的增加，從而削弱企業(yè)的創(chuàng)新能力和國際競爭力[3-5]。另一方面，環(huán)境規(guī)制和企業(yè)競爭力可以實現(xiàn)雙贏。PORTER et al.[6-7]認為，合理的環(huán)境規(guī)制能夠激發(fā)企業(yè)的創(chuàng)新意識，促使企業(yè)主動優(yōu)化資源配置效率，改造生產(chǎn)工藝或開發(fā)新產(chǎn)品，并通過創(chuàng)新補償效應(yīng)獲取更高收益，提高企業(yè)整體競爭力。

可見，環(huán)境稅開征是否有助于完成環(huán)境保護和經(jīng)濟發(fā)展的雙重任務(wù)，關(guān)鍵在于創(chuàng)新補償效應(yīng)的出現(xiàn)[8]。也就是說，環(huán)境稅設(shè)計應(yīng)以誘導(dǎo)企業(yè)環(huán)境創(chuàng)新發(fā)展為重要目的[9]，而不只是稅收所得。梁偉等[10]認為，簡單意義上的環(huán)境稅征收很難實現(xiàn)節(jié)能減排和經(jīng)濟增長的“雙重紅利”，環(huán)境稅會給宏觀經(jīng)濟帶來一定的負面影響；OUESLATI[11]和ABDULLAH et al.[12]的研究證實，長期看，環(huán)境稅可以增加財政收入，促進經(jīng)濟增長，但短期內(nèi)其總的社會福利效應(yīng)是負的。因此，從促進經(jīng)濟可持續(xù)發(fā)展的視角出發(fā)，ASHFORD et al.[13]和HAHN[14]強調(diào)環(huán)境規(guī)制制定者必須根據(jù)不同行業(yè)的可替代技術(shù)發(fā)展條件，認真考慮規(guī)制設(shè)計的嚴格性、靈活性和實施的時機；許士春等[15]認為，企業(yè)減排的主要影響因素是綠色技術(shù)創(chuàng)新能力和環(huán)境規(guī)制的嚴厲程度；而BERGEK et al.[16]研究發(fā)現(xiàn)，環(huán)境規(guī)制在環(huán)境創(chuàng)新技術(shù)擴散中發(fā)揮著關(guān)鍵性的作用，經(jīng)濟政策有利于激勵漸進式創(chuàng)新，而只有嚴格的環(huán)境規(guī)制可以推動新技術(shù)研發(fā)。同時，環(huán)境政策必須有一定的連續(xù)性，VAN VLIET et al.[17]研究發(fā)現(xiàn)，持續(xù)的政策工具在可替代能源采納中發(fā)揮了關(guān)鍵性的作用。由此可見，行業(yè)特點、技術(shù)條件、開征時機及稅率設(shè)置等多重因素影響環(huán)境稅在綠色經(jīng)濟發(fā)展中扮演的角色。

在環(huán)境稅的定量研究中，通常采用可計算的一般均衡方法(computable general equilibrium，CGE)作為政策模擬工具，研究環(huán)境稅對環(huán)境保護和宏觀經(jīng)濟的影響。然而，不同學(xué)者通過CGE模型得出不同的研究結(jié)論。李洪心等[18]的研究結(jié)果表明環(huán)境稅能夠產(chǎn)生雙贏效應(yīng)，而劉曄等[19]認為環(huán)境稅的實施與環(huán)境質(zhì)量的改善并沒有必然的聯(lián)系。污染大戶中的“富人”有可能用更多的錢換取對環(huán)境污染的超容，環(huán)境稅還有可能被轉(zhuǎn)嫁給消費者。何建武等[20]研究發(fā)現(xiàn)，以工業(yè)污染物排放全部達標為標準征收環(huán)境稅，將影響經(jīng)濟增長1%左右，并使制造業(yè)就業(yè)量減少1.8%，出口量下降1.7%；而LEE et al.[21]認為，環(huán)境稅的征收涉及到各方利益，是一個復(fù)雜的社會問題，環(huán)境稅開征后能夠?qū)崿F(xiàn)的減排效果以及對社會和經(jīng)濟產(chǎn)生的影響，在不同的國家、地區(qū)和行業(yè)是有差異的。

已有研究從宏觀層面探討環(huán)境稅對環(huán)境保護和宏觀經(jīng)濟的影響效果，但從微觀層面揭示環(huán)境稅征收對企業(yè)環(huán)境技術(shù)創(chuàng)新行為的影響機制、分析環(huán)境稅征收時機及強度等對創(chuàng)新補償效應(yīng)影響的研究尚不多見。本研究探尋不同環(huán)境稅情景下，企業(yè)選擇環(huán)境技術(shù)創(chuàng)新行為的動因，揭示在環(huán)境稅、市場機制及產(chǎn)品競爭的協(xié)同作用下，環(huán)境技術(shù)創(chuàng)新影響企業(yè)環(huán)境和經(jīng)濟績效的機制，打開企業(yè)的“黑箱”，探討產(chǎn)業(yè)向綠色可持續(xù)方向發(fā)展的動力源泉。該類研究難以通過面板數(shù)據(jù)得到準確答案，社會科學(xué)計算實驗方法[22-23]可以在調(diào)查研究取得實證數(shù)據(jù)的基礎(chǔ)上，通過構(gòu)建具有自適應(yīng)能力的多智能體模型，模擬外部環(huán)境變化對系統(tǒng)內(nèi)各主體行為及決策的影響，進而觀察整個系統(tǒng)的宏觀“涌現(xiàn)”。作為一種自下而上的研究工具，計算實驗已應(yīng)用于社會科學(xué)的多個研究領(lǐng)域[24-27]。本研究通過構(gòu)建企業(yè)環(huán)境技術(shù)創(chuàng)新行為動態(tài)仿真模型，對企業(yè)環(huán)境技術(shù)創(chuàng)新行為與其內(nèi)生要素和外部驅(qū)動之間的相互作用機理進行動態(tài)仿真，探索環(huán)境技術(shù)創(chuàng)新行為的動因，為制定政策提供參考依據(jù)，進而通過相應(yīng)政策的實施，推動企業(yè)自覺采納環(huán)境創(chuàng)新技術(shù)，步入綠色可持續(xù)發(fā)展之路，并實現(xiàn)環(huán)境績效和經(jīng)濟績效的雙贏。

3 模型構(gòu)建

3.1 情景描述

本研究以對企業(yè)生產(chǎn)過程中排放的CO2征收環(huán)境稅為背景，研究征稅前后企業(yè)環(huán)境技術(shù)創(chuàng)新行為的演變，并分析環(huán)境稅征收時機和強度對可替代環(huán)境創(chuàng)新技術(shù)發(fā)展和擴散的影響。系統(tǒng)初始階段該行業(yè)內(nèi)所有企業(yè)均采用傳統(tǒng)技術(shù)，生產(chǎn)工藝較為成熟，生產(chǎn)成本較低，產(chǎn)品質(zhì)量較好，但CO2排放水平較高。雖然企業(yè)可以通過末端治理等手段減排，但受限于技術(shù)手段，減排效果有限。與之對應(yīng)的環(huán)境創(chuàng)新技術(shù)CO2排放較低，但可替代原材料及創(chuàng)新技術(shù)的采用需要一定的設(shè)備投入和生產(chǎn)流程改造，由于生產(chǎn)工藝尚不成熟等原因，生產(chǎn)成本較高，且產(chǎn)品質(zhì)量不如傳統(tǒng)技術(shù)產(chǎn)品。但可以通過技術(shù)研發(fā)和改造等手段不斷改良生產(chǎn)工藝，從而降低生產(chǎn)成本，提高產(chǎn)品質(zhì)量。

環(huán)境政策設(shè)計的重要目的在于能夠誘導(dǎo)企業(yè)自覺選擇環(huán)境創(chuàng)新技術(shù)，孫曉華等[28]研究發(fā)現(xiàn)，消費者異質(zhì)性偏好為新產(chǎn)品提供了利基市場，在產(chǎn)業(yè)演化中發(fā)揮重要作用；趙愛武等[29]進一步從消費者價格敏感度、從眾心理、環(huán)保態(tài)度和消費習(xí)慣等要素出發(fā)，探索消費者綠色購買行為中各要素的作用機制，分析企業(yè)選擇環(huán)境創(chuàng)新技術(shù)的市場基礎(chǔ)，發(fā)現(xiàn)對不同政策設(shè)計效率的考察需要結(jié)合企業(yè)間競爭和消費者市場選擇等多種復(fù)雜的自相關(guān)演化機制。由于實證方法難以找到不同環(huán)境政策背景的對照樣本，因此，本研究采用社會科學(xué)計算實驗方法，參考AFAROUI et al.[30]和劉小峰等[31]構(gòu)建多智能體模型的部分設(shè)計思路，根據(jù)現(xiàn)實原型構(gòu)建計算實驗?zāi)Ｐ?，模擬企業(yè)在環(huán)境稅、市場機制和產(chǎn)品競爭復(fù)雜環(huán)境下的環(huán)境技術(shù)創(chuàng)新行為動態(tài)過程。

本模型中主要包含消費者和生產(chǎn)企業(yè)兩大類主體，消費者根據(jù)自身偏好選擇不同生產(chǎn)企業(yè)的產(chǎn)品，而生產(chǎn)企業(yè)間自由競爭，并根據(jù)各自的決策規(guī)則選擇產(chǎn)品技術(shù)路線。為更好地觀察不同環(huán)境政策情景下生產(chǎn)企業(yè)的環(huán)境創(chuàng)新軌跡，本模型盡可能考慮現(xiàn)實系統(tǒng)的復(fù)雜性，并對其進行抽象和簡化。

本研究所建模型旨在作為研究工具模擬不同環(huán)境稅情景下、在行業(yè)競爭和市場選擇的共同作用下環(huán)境創(chuàng)新技術(shù)的擴散過程，模型僅提供對復(fù)雜環(huán)境各影響因素協(xié)同作用下企業(yè)采納環(huán)境創(chuàng)新技術(shù)的微觀動力機制和環(huán)境創(chuàng)新技術(shù)宏觀動態(tài)擴散路徑的模擬，因此，對主體個性的設(shè)置盡可能簡化。如系統(tǒng)初始階段僅考慮自由競爭市場結(jié)構(gòu)，并假設(shè)企業(yè)主體規(guī)模和資金狀況相同。同時，本研究討論的環(huán)境技術(shù)創(chuàng)新屬于生產(chǎn)工藝創(chuàng)新，考慮到信息不對稱情景下消費者不了解企業(yè)內(nèi)部生產(chǎn)工藝細節(jié)，因此，沒有考慮消費者對環(huán)保產(chǎn)品的偏好屬性。另外，對產(chǎn)品不變報酬型的設(shè)置參考劉小峰等[31]構(gòu)建的生產(chǎn)企業(yè)模型，企業(yè)技術(shù)創(chuàng)新過程的設(shè)計及部分參數(shù)的設(shè)置參考AFAROUI et al.[30]構(gòu)建的化工企業(yè)技術(shù)改造模型。

3.2 基本假設(shè)

根據(jù)3.1中的情景描述，系統(tǒng)基本假設(shè)如下。

(1)G1為傳統(tǒng)技術(shù)，G2為環(huán)境創(chuàng)新技術(shù)，兩種技術(shù)生產(chǎn)產(chǎn)品分別為G1技術(shù)產(chǎn)品和G2技術(shù)產(chǎn)品，兩類產(chǎn)品的價格和質(zhì)量不同，生產(chǎn)過程中單位產(chǎn)品的CO2排放也不同，且通過技術(shù)改造能夠達到的最低生產(chǎn)成本、最高產(chǎn)品質(zhì)量和最低CO2排放也不同。

(2)系統(tǒng)中共有m家生產(chǎn)企業(yè)，分別以企業(yè)主體編號i表示，i=1, 2, …,m，企業(yè)初始階段規(guī)模相同，全部采用G1技術(shù)，但由于G2技術(shù)代表了未來發(fā)展的方向，根據(jù)i企業(yè)對G2技術(shù)的關(guān)注程度，將一定比例的R&D投入用于G2技術(shù)前期研發(fā)和學(xué)習(xí)。當(dāng)滿足一定條件時，企業(yè)開始采納G2技術(shù)正式生產(chǎn)，但G1技術(shù)可以同時并存，達到一定條件時企業(yè)放棄G1技術(shù)，企業(yè)對G2技術(shù)的采納閾值A(chǔ)dmin和G1技術(shù)的淘汰閾值A(chǔ)bandmin不同。為保證系統(tǒng)中生產(chǎn)企業(yè)總數(shù)穩(wěn)定，維持自由競爭格局，假設(shè)當(dāng)生產(chǎn)企業(yè)虧損達到一定程度時退出市場，同時有新的加入者進入市場，不考慮企業(yè)獲取貸款等外部資金的情況。

(3)系統(tǒng)中共有n個消費者，分別以消費者主體編號j表示，j=1, 2, …,n。假設(shè)產(chǎn)品為易耗必需品，消費者每個周期均需購買一件產(chǎn)品。由于消費者不了解生產(chǎn)工藝細節(jié)，根據(jù)ZEPPINI et al.[32]的模型，創(chuàng)新產(chǎn)品擴散中消費者主要考慮產(chǎn)品價格、性能、消費習(xí)慣和產(chǎn)品口碑，并具有一定的社會模仿能力。因此，假定消費者根據(jù)產(chǎn)品價格和質(zhì)量選擇不同企業(yè)、不同技術(shù)的產(chǎn)品，而不同消費者對產(chǎn)品價格和質(zhì)量的偏好不同，且購買決策受其他消費者影響。同時，考慮到消費者受消費習(xí)慣的影響，具有一定的路徑依賴屬性，當(dāng)前期選擇的生產(chǎn)企業(yè)產(chǎn)品價格和質(zhì)量在可容忍范圍內(nèi)時(以消費者容忍度Tol表示)，消費者本期仍選擇原企業(yè)，消費者對產(chǎn)品價格和質(zhì)量的容忍度不同。

(4)產(chǎn)品為不變報酬型，即生產(chǎn)效率不會隨生產(chǎn)規(guī)模的擴大而提高，只有通過技術(shù)改造才能提高生產(chǎn)效率，從而降低生產(chǎn)成本。因此，假定產(chǎn)品定價基于生產(chǎn)成本，即p=c(1+μ)，p為產(chǎn)品價格，c為生產(chǎn)成本，μ為生產(chǎn)企業(yè)的滿意利潤水平(考慮生產(chǎn)者為有限理性)。假定所有生產(chǎn)企業(yè)滿意利潤水平一致，可以通過生產(chǎn)企業(yè)的生產(chǎn)效率水平反映其產(chǎn)品價格水平。

3.3 主體行為規(guī)則

3.3.1 消費者主體行為規(guī)則

消費者選擇產(chǎn)品時同時考慮產(chǎn)品價格和質(zhì)量兩方面因素，且只有當(dāng)兩者達到消費者對價格的最高承受能力和對質(zhì)量的最低要求時，消費者才會考慮購買。同時，不同消費者對價格和質(zhì)量的偏好不同，且受其他消費者選擇影響表現(xiàn)出從眾傾向。因此，消費者選擇產(chǎn)品的效應(yīng)函數(shù)為

{(Xk,i,t-QT)(Msi,t-1+u(0,0.10)]inf}ξXj

(1)

圖1 消費者主體產(chǎn)品選擇流程圖Figure 1 Product Choice Flow Chart of Consumer Agent

首次選擇產(chǎn)品時，消費者根據(jù)對每種產(chǎn)品的效應(yīng)函數(shù)，確定其選擇概率，隨機選擇產(chǎn)品。后續(xù)周期，消費者依據(jù)路徑依賴原則，首先觀察當(dāng)期產(chǎn)品最低價格和最高質(zhì)量，當(dāng)原來選擇的產(chǎn)品性價比在消費者可容忍范圍內(nèi)時，消費者選擇原生產(chǎn)企業(yè)的產(chǎn)品；否則，根據(jù)各產(chǎn)品效應(yīng)函數(shù)重新確定產(chǎn)品選擇概率，隨機選擇產(chǎn)品。

消費者主體各周期的產(chǎn)品選擇流程見圖1。

3.3.2 生產(chǎn)企業(yè)主體行為規(guī)則

根據(jù)假設(shè)，初始階段系統(tǒng)中所有生產(chǎn)企業(yè)規(guī)模和資金狀況相同，每個模擬周期，企業(yè)從生產(chǎn)銷售活動中獲取利潤，并為提高產(chǎn)品競爭力開展技術(shù)改造。本模型不考慮貸款等其他途徑獲取資金的情況，則生產(chǎn)企業(yè)的各周期可支配資金總額為

Bi,t=Bi,t-1+Πi,t-1-RDi,t-1

(2)

其中，Bi,t為i企業(yè)在t周期的可支配資金總額，Πi,t-1為i企業(yè)在(t-1)周期的當(dāng)期收益，RDi,t-1為i企業(yè)在(t-1)周期的R&D支出。

對于新的G2技術(shù)采納者，需要支付額外的技術(shù)轉(zhuǎn)換成本(如設(shè)備投入、人員培訓(xùn)費用等)，因此

Bi,t=Bi,t-1+Πi,t-1-RDi,t-1-SCi,t

(3)

其中，SCi,t為i企業(yè)在t周期采納G2技術(shù)時的相關(guān)轉(zhuǎn)換成本。

生產(chǎn)企業(yè)主體各周期利潤為

Πi,t=μ·Ci,t·Qi,t-FC

(4)

其中，Ci,t為i企業(yè)在t周期的產(chǎn)品變動成本，Qi,t為i企業(yè)在t周期的產(chǎn)品銷量，F(xiàn)C為產(chǎn)品固定成本。

(1)生產(chǎn)企業(yè)進入/退出規(guī)則

當(dāng)生產(chǎn)企業(yè)可支配資金小于一定水平時，企業(yè)宣告破產(chǎn)并從市場退出，同時一個新的生產(chǎn)企業(yè)進入市場。新進入的生產(chǎn)企業(yè)所采用的技術(shù)路線模仿系統(tǒng)中現(xiàn)存的某個企業(yè)，模仿目標的選擇概率基于各企業(yè)的市場份額。新企業(yè)模仿目標企業(yè)的技術(shù)路線，其學(xué)習(xí)吸收能力為介于[0.80,1.20]的隨機數(shù)[30]，產(chǎn)品價格、質(zhì)量和CO2排放等各項指標在被模仿企業(yè)產(chǎn)品指標基礎(chǔ)上乘或除(正指標乘，逆指標除)以隨機數(shù)，因此可以低于或超越被模仿企業(yè)。新企業(yè)的初始可支配資金B(yǎng)和固定成本FC與其他企業(yè)初始時類似，知識存量K和G2技術(shù)的轉(zhuǎn)換成本SC取行業(yè)平均值。

(2)生產(chǎn)企業(yè)技術(shù)路線選擇規(guī)則

每周期企業(yè)計算其對G2技術(shù)成熟度的感知，即

(5)

采納G2技術(shù)并不意味著一定放棄G1技術(shù)，根據(jù)假設(shè)，在同一企業(yè)兩者可以并存。是否淘汰G1技術(shù)，取決于G2技術(shù)的產(chǎn)品收入在企業(yè)總收入中所占的比重，即

(6)

(3)企業(yè)R&D活動規(guī)則

每周期各企業(yè)通過R&D活動提高產(chǎn)品各項性能，R&D投資額為

RDi,t=δ·Bi,t

(7)

其中，RDi,t為i企業(yè)在t周期的R&D支出，δ為R&D投資比例，前提是企業(yè)當(dāng)期可支配資金B(yǎng)i,t>0。

企業(yè)的R&D投資按比例分別用于G1和G2技術(shù)研發(fā)，即

(8)

(9)

1-exp(-αw·RDwatchi,t)≥u(0,1)

(10)

其中，αw為模式參數(shù)，反映前期研發(fā)和學(xué)習(xí)取得實質(zhì)性進展的難易系數(shù)。u均勻隨機分布在[0,1]之間，反映了現(xiàn)實世界中創(chuàng)新活動的不確定性。如果滿足條件，表明研發(fā)活動獲得階段性成果，G2技術(shù)知識存量增加，且轉(zhuǎn)換為G2技術(shù)的成本降低。知識存量和G2技術(shù)轉(zhuǎn)換成本為

Ki,t=Ki,t-1+αk·u(0,1)·(Kmax-Ki,t-1)

(11)

SCi,t=SCi,t-1-αSC·u(0,1)·(SCi-1-SCmin)

(12)

其中，αk為模式參數(shù)，反映知識存量增加速度；αSC為模式參數(shù)，反映轉(zhuǎn)換成本降低速度；Kmax為知識存量的最大值，SCmin為轉(zhuǎn)換成本的最小值。

生產(chǎn)活動中的技術(shù)改造過程與G2技術(shù)前期研發(fā)和學(xué)習(xí)過程類似，而技術(shù)改造能否成功取決于是否滿足

(13)

(14)

(15)

(16)

生產(chǎn)企業(yè)主體各周期的工作流程見圖2。

3.4 參數(shù)設(shè)置

設(shè)置系統(tǒng)模擬需要的公共參數(shù)及各主體的個性化參數(shù)等盡量考慮其客觀現(xiàn)實，對于現(xiàn)實中難以具體量化的參數(shù)(如消費者對產(chǎn)品價格和質(zhì)量偏好的量化)，采用基本模型(無政策干預(yù)場景)進行虛實聯(lián)動，通過調(diào)整參數(shù)反復(fù)模擬，觀察中間結(jié)果和最終結(jié)果，使各項結(jié)果與現(xiàn)實情況吻合。確定模型基本參數(shù)后，再引入政策參數(shù)，觀察政策工具對模擬結(jié)果的影響。

3.5 參數(shù)討論

具體參數(shù)值的設(shè)置影響系統(tǒng)的具體模擬結(jié)果，但在一定范圍內(nèi)并不會改變系統(tǒng)的演化趨勢，每組參數(shù)均有一定的適用范圍。

(1)G1和G2技術(shù)產(chǎn)品質(zhì)量和成本屬性

G1和G2技術(shù)產(chǎn)品的質(zhì)量和成本(決定產(chǎn)品價格)影響消費者選擇。根據(jù)假設(shè)，G1和G2技術(shù)產(chǎn)品能夠達到的質(zhì)量上限和成本下限相同，因此，兩類技術(shù)產(chǎn)品的質(zhì)量和成本極限值相同。但兩者初始值不同，G1技術(shù)產(chǎn)品的技術(shù)成熟、質(zhì)量和成本均接近極限值，而G2技術(shù)產(chǎn)品質(zhì)量和成本的初始值需要設(shè)置。分別以G1與G2技術(shù)產(chǎn)品質(zhì)量或成本的比值反映兩類技術(shù)產(chǎn)品質(zhì)量或成本的差異，設(shè)G1與G2技術(shù)產(chǎn)品質(zhì)量初始值比值為φ1，G1與G2技術(shù)產(chǎn)品成本初始值比值為φ2，則當(dāng)其他參數(shù)不變時，φ1和φ2對產(chǎn)品最終市場份額的影響見圖3(50次模擬的平均值)。

圖2 生產(chǎn)企業(yè)主體工作流程圖Figure 2 Work Flow Chart of Enterprise Agent

表1 主要變量及初始賦值規(guī)則Table 1 Main Parameters and Initial Value

圖3 產(chǎn)品質(zhì)量和成本初始設(shè)置敏感測試結(jié)果Figure 3 Results of Sensitivity Test for Initial Setting of Product Quality and Cost

(2)G1和G2技術(shù)產(chǎn)品CO2排放量屬性

(3)其他參數(shù)的敏感度測試

企業(yè)初始資金與固定成本的比值可以反映企業(yè)的相對財務(wù)狀況，敏感度測試發(fā)現(xiàn)，當(dāng)兩者的比值為[4,30]時，可以維持無環(huán)境稅情景下1%～5%的淘汰率，這與實證結(jié)果一致，且該比值范圍內(nèi)模擬實驗結(jié)論一致。本研究設(shè)B0=12，F(xiàn)C=2，即企業(yè)初始資金與固定成本的比值為6，在[4,30]范圍內(nèi)。另外，新技術(shù)轉(zhuǎn)換成本SC∈[12,40]，相應(yīng)的SCmax∈[6,30]時，測試結(jié)果與本研究設(shè)置的SC=20、SCmin=10一致。

其他參數(shù)的設(shè)置基于上述參數(shù)的具體設(shè)置值，并考慮其合理性，經(jīng)過反復(fù)實驗獲取。如αw，當(dāng)取值為[0.10,0.40]時，可以使新產(chǎn)品轉(zhuǎn)換成本的降低速度維持在合理水平。

4 模擬實驗和結(jié)果分析

本模型以無政策工具干預(yù)的基本情景為基礎(chǔ)，主要分析環(huán)境稅征收對企業(yè)環(huán)境技術(shù)創(chuàng)新行為和環(huán)境創(chuàng)新技術(shù)采納的影響。根據(jù)環(huán)境稅開征時機和征收強度設(shè)計如下場景。

O為無環(huán)境稅征收的基本情景。

H0為系統(tǒng)初始階段開始征收環(huán)境稅，此時可替代的環(huán)境創(chuàng)新技術(shù)(G2技術(shù))轉(zhuǎn)換成本較高，且由于創(chuàng)新技術(shù)發(fā)展時間短，G2技術(shù)產(chǎn)品成本較高且質(zhì)量較差；根據(jù)環(huán)境稅征收強度，設(shè)置H0_0.20、H0_0.40、H0_0.60、H0_0.80、H0_1.00和H0_1.50等子場景，分別對應(yīng)環(huán)境稅征收強度(或稅率)為0.20、0.40、0.60、0.80、1.00和1.50的情形。

H100為在環(huán)境稅征收預(yù)期(未正式開征)下，系統(tǒng)運行100周期，并從101周期開始正式征收環(huán)境稅，此時可替代的環(huán)境創(chuàng)新技術(shù)(G2技術(shù))轉(zhuǎn)換成本較低，且創(chuàng)新技術(shù)已得到一定發(fā)展，G2技術(shù)產(chǎn)品成熟度較高；根據(jù)環(huán)境稅征收強度設(shè)置H100_0.20、H100_0.40、H100_0.60、H100_0.80、H100_1.00和H100_1.50等子場景。

4.1 不同情景下企業(yè)技術(shù)路線及排放結(jié)果分析

環(huán)境稅征收情景下，當(dāng)模擬周期到達征收時機設(shè)定值時，開始對生產(chǎn)過程中的CO2排放征收環(huán)境稅。因此，當(dāng)期利潤(4)式調(diào)整為

(a)H0情景下企業(yè)技術(shù)路線 (b)H100情景下企業(yè)技術(shù)路線

(c)企業(yè)總排放量 (d)300周期時企業(yè)當(dāng)期排放量

Πi,t=μ·Ci,t·Qi,t-FC-CO2i,t·Qi,t·f

(17)

其中，f為單位產(chǎn)品單位CO2排放征收的環(huán)境稅，代表環(huán)境稅的征收強度(或稅率)。

當(dāng)f分別為0、0.20、0.40、0.60、0.80、1.00和1.50，環(huán)境稅實施時機分別為系統(tǒng)開始階段(H0)和某中間階段(H100)時，其他參數(shù)設(shè)置與基本情景相同(此時f=0，即基本情景O)，系統(tǒng)運行300周期，模擬結(jié)果見圖4。

由圖4(a)和圖4(b)可知，采納環(huán)境創(chuàng)新技術(shù)的企業(yè)數(shù)量隨環(huán)境稅征收強度增加，而由圖4(c)和圖4(d)可知，隨著征收強度的增大，企業(yè)總排放量和300周期的當(dāng)期排放量呈下降趨勢，并當(dāng)征收強度大于0.80后，排放量急劇下降；由圖4(a)和圖4(b)可知，當(dāng)征收強度達到1.00時，所有企業(yè)均采納環(huán)境創(chuàng)新技術(shù)；由圖4(c)可知，當(dāng)征收強度小于或等于0.60時，CO2排放量雖小于未征收環(huán)境稅(圖中征稅強度為0)時的情景，但在0.20～0.80時排放量隨征收強度變化不大；對比兩種環(huán)境稅征收時機，從圖4(a)、圖4(b)和圖4(d)可知，征收強度大于0.80時，兩類情景下的企業(yè)技術(shù)路線和300周期當(dāng)期排放量均差別不大，但圖4(c)顯示，總排放量H0情景小于H100情景；而征收強度小于1.00時，雖然圖4(d)顯示兩類情景下300周期當(dāng)期企業(yè)排放量差別不大，但H0情景下有更多企業(yè)最終選擇環(huán)境創(chuàng)新技術(shù)；不過意外的是，圖4(c)中總排放量在征收強度小于0.60時H0>H100，即早期采用環(huán)境稅情景下的總排放量反而高于100周期后才開始施行環(huán)境稅情景下的總排放量。為進一步分析原因，選取征收強度為0.40時，對比H0和H100下各周期G2技術(shù)產(chǎn)品所占市場份額，見圖5。

圖5 環(huán)境稅0.40下G2技術(shù)產(chǎn)品市場占有率Figure 5 Market Share of G2under Environmental Tax of 0.40

(a)300周期時累計淘汰企業(yè) (b)市場集中度

(c)H0情景下產(chǎn)業(yè)結(jié)構(gòu) (d)H100情景下產(chǎn)業(yè)結(jié)構(gòu)

由圖5可知，H0和H100情景下G2技術(shù)產(chǎn)品市場占有率差別不大，均小于40%，也就是說，更多的企業(yè)選擇傳統(tǒng)技術(shù)進行生產(chǎn)。然而，在H0情景下，由于環(huán)境稅征收降低了企業(yè)技術(shù)改造的投入，在不采用環(huán)境創(chuàng)新技術(shù)的情況下，企業(yè)實際上是在花錢買污染，因而圖4(c)中早期開始征收環(huán)境稅的H0情景總排放量反而大于100周期后才開始征收環(huán)境稅的H100情景?？梢?，只有當(dāng)環(huán)境稅強度達到一定水平，才能促使企業(yè)積極選擇環(huán)境創(chuàng)新技術(shù)，真正實現(xiàn)減排的目標。

4.2 環(huán)境稅影響下的產(chǎn)業(yè)演化分析

通過4.1的分析可知，只有高強度的環(huán)境稅才能真正推動企業(yè)采納環(huán)境創(chuàng)新技術(shù)，但過高的環(huán)境稅會加重企業(yè)負擔(dān)，并導(dǎo)致企業(yè)面臨生存危機。不同情景下企業(yè)生存狀況及環(huán)境創(chuàng)新技術(shù)產(chǎn)業(yè)演化結(jié)果見圖6。

由圖6(a)可知，在H0情景下，當(dāng)環(huán)境稅強度大于0.80時，累計被淘汰企業(yè)大量增加，表明環(huán)境稅的強度已超出企業(yè)的承受能力；在H100情景下，環(huán)境稅強度的增加對累計淘汰企業(yè)數(shù)量影響不大。這是由于H0情景下，征收高額環(huán)境稅雖然能夠促使企業(yè)采納環(huán)境創(chuàng)新技術(shù)，但由于環(huán)境創(chuàng)新技術(shù)不成熟，市場競爭力差，因此大量先期采納環(huán)境創(chuàng)新技術(shù)的企業(yè)被淘汰。現(xiàn)實場景中，這勢必損害企業(yè)環(huán)境技術(shù)創(chuàng)新的積極性，且增加產(chǎn)業(yè)風(fēng)險，難以出現(xiàn)模擬中大量后續(xù)企業(yè)進入的情況，也難以出現(xiàn)最終環(huán)境創(chuàng)新技術(shù)完全擴散的情景。由圖6(b)可知，環(huán)境稅強度在0.40～0.80之間時，H0情景下市場集中度(以赫芬達爾-赫希曼指數(shù)的倒數(shù)衡量)提高，出現(xiàn)少量企業(yè)壟斷市場的局面，且由圖6(c)可知，壟斷企業(yè)為采用組合技術(shù)(G1/G2)的企業(yè)，而完全采納環(huán)境創(chuàng)新技術(shù)的企業(yè)數(shù)量盡管已超過半數(shù)(圖4(a))，但所占市場份額較低，企業(yè)生存狀況欠佳，面臨被“擠出”的風(fēng)險。而由圖6(b)可知，H100情景下，環(huán)境稅強度大于1.00時，才出現(xiàn)壟斷局面，而由圖6(d)可知，壟斷企業(yè)全部采用環(huán)境創(chuàng)新技術(shù)。也就是說，環(huán)境創(chuàng)新技術(shù)獲得成功擴散。

(a)O情景下各周期企業(yè)累計利潤變化情況

(b)H100_1.00情景下各周期企業(yè)累計利潤變化情況

可見，環(huán)境稅政策的實施時機對于產(chǎn)業(yè)向環(huán)境創(chuàng)新技術(shù)方向演化起著關(guān)鍵作用。當(dāng)以促進企業(yè)環(huán)境創(chuàng)新技術(shù)采納為前提時，應(yīng)該選擇環(huán)境創(chuàng)新技術(shù)較成熟時推出環(huán)境稅，且環(huán)境稅需要達到一定強度才能真正發(fā)揮作用。綜合分析對不同情景的模擬，由圖4(b)可知，H100_1.00情景下征稅強度為1.00時全部企業(yè)采用G2技術(shù)；由圖4(d)可知，該情景下排放量最低；由圖6(d)可知，該情景下產(chǎn)業(yè)結(jié)構(gòu)為綠色；由圖6(a)可知，此時淘汰率也不高。所以，H100_1.00是最優(yōu)選擇，能夠有效促進環(huán)境創(chuàng)新技術(shù)擴散，降低CO2排放。

4.3 最優(yōu)情景下企業(yè)經(jīng)濟效益分析

根據(jù)上述分析，環(huán)境創(chuàng)新技術(shù)發(fā)展到一定階段后，較強的環(huán)境稅情景H100_1.0在產(chǎn)業(yè)轉(zhuǎn)型升級和減排方面最優(yōu)。對比無環(huán)境稅的O情景和H100_1.00情景下各周期企業(yè)累計利潤變化情況，見圖7。

由圖7可知，O情景下企業(yè)各周期平均累計收益基本在10左右，而H100_1.00情景下，盡管企業(yè)間收益差距拉大，但平均累計收益大于20。也就是說，環(huán)境稅的征收使企業(yè)累計收益大幅提高，這也在一定程度上驗證了波特假說，即企業(yè)通過創(chuàng)新補償，獲得了環(huán)境績效和經(jīng)濟績效的雙贏。

5 結(jié)論

本研究運用社會科學(xué)計算實驗方法，模擬環(huán)境稅征收不同情景下企業(yè)環(huán)境創(chuàng)新行為的動態(tài)選擇過程，并分析不同情景下企業(yè)排放水平和經(jīng)濟效益的變化。模擬實驗從征收時機和稅率設(shè)置兩個維度探討環(huán)境稅對企業(yè)環(huán)境技術(shù)創(chuàng)新行為選擇以及對企業(yè)環(huán)境效益和經(jīng)濟效益的短期和長期影響，并分析環(huán)境稅與產(chǎn)業(yè)演化的關(guān)系。研究結(jié)果表明，較高的環(huán)境稅率對企業(yè)環(huán)境技術(shù)創(chuàng)新行為具有促進作用，并能夠促使企業(yè)向環(huán)境創(chuàng)新技術(shù)轉(zhuǎn)型。但在環(huán)境創(chuàng)新技術(shù)尚不成熟時，盲目的技術(shù)轉(zhuǎn)換會增加企業(yè)經(jīng)營風(fēng)險，并損傷企業(yè)環(huán)境技術(shù)創(chuàng)新的積極性，反而阻礙環(huán)境創(chuàng)新技術(shù)的發(fā)展和擴散，不利于減排和經(jīng)濟發(fā)展。實驗結(jié)果為環(huán)境技術(shù)創(chuàng)新企業(yè)和政府管理部門決策者提供了積極的管理啟示。

(1)環(huán)境創(chuàng)新技術(shù)的采納應(yīng)考慮技術(shù)成熟度。當(dāng)環(huán)境創(chuàng)新技術(shù)尚不成熟時，高強度的環(huán)境稅雖然能促使企業(yè)采納環(huán)境創(chuàng)新技術(shù)，但早期的少數(shù)環(huán)境技術(shù)創(chuàng)新企業(yè)在與傳統(tǒng)技術(shù)企業(yè)的競爭中處于劣勢，甚至慘遭淘汰(H0下征收強度大于0.80的情景)，這必將影響企業(yè)環(huán)境技術(shù)創(chuàng)新的積極性。因此，制定環(huán)境稅政策應(yīng)充分考慮征收企業(yè)的行業(yè)特點，結(jié)合可替代環(huán)境創(chuàng)新技術(shù)的發(fā)展水平，從加大環(huán)境創(chuàng)新技術(shù)的引進和研發(fā)力度入手，積極扶持環(huán)境技術(shù)創(chuàng)新企業(yè)，促使環(huán)境創(chuàng)新技術(shù)盡快走向成熟。

(2)當(dāng)可替代環(huán)境創(chuàng)新技術(shù)成熟時，可選準時機出臺較高強度的環(huán)境稅，通過市場機制和企業(yè)間競爭等的協(xié)同作用，促使企業(yè)主動選擇環(huán)境創(chuàng)新技術(shù)，并實現(xiàn)環(huán)境績效和經(jīng)濟績效的雙贏。環(huán)境稅征收的強度要能夠達到促使企業(yè)主動放棄傳統(tǒng)技術(shù)，否則可能出現(xiàn)企業(yè)花錢買污染的狀況，一定程度上還會降低企業(yè)技術(shù)改造的能力，增加總排放量(H0下征收強度小于0.60的情景)，與環(huán)境稅的初衷背道而馳。適當(dāng)?shù)沫h(huán)境稅征收時機和稅率設(shè)計可以在產(chǎn)業(yè)演化中發(fā)揮積極作用，通過環(huán)境創(chuàng)新技術(shù)替代傳統(tǒng)技術(shù)，可促使產(chǎn)業(yè)向綠色經(jīng)濟方向發(fā)展。

(3)消費者對環(huán)境創(chuàng)新技術(shù)的選擇會直接影響系統(tǒng)的演化結(jié)果，在征收環(huán)境稅的同時，可通過對環(huán)境技術(shù)創(chuàng)新產(chǎn)品補貼等方式，降低產(chǎn)品成本，提高消費者對環(huán)境技術(shù)創(chuàng)新產(chǎn)品的選擇概率。同時，隨著消費者環(huán)保意識水平的提高，可以通過增加信息披露力度的方式，降低消費者與企業(yè)的信息不對稱程度，引導(dǎo)消費者主動選擇環(huán)境友好產(chǎn)品。

本研究尚存在不足，需要在后續(xù)研究中進一步深入。①本實驗假設(shè)系統(tǒng)初始階段企業(yè)規(guī)模和資金狀況相同，沒有考慮不同市場結(jié)構(gòu)的情況；②模型設(shè)計中，企業(yè)進入機制設(shè)計時沒有考慮產(chǎn)業(yè)風(fēng)險對企業(yè)進入意愿的影響；③消費者產(chǎn)品購買決策僅考慮了產(chǎn)品價格、質(zhì)量、消費習(xí)慣和從眾效應(yīng)，沒有考慮其對產(chǎn)品環(huán)境屬性的偏好。

[1]張紅鳳,周峰,楊慧,等.環(huán)境保護與經(jīng)濟發(fā)展雙贏的規(guī)制績效實證分析.經(jīng)濟研究,2009,44(3):14-26.

ZHANG Hongfeng,ZHOU Feng,YANG Hui,et al.Regulation performance of the win-win of environmental protection and economic development.EconomicResearchJournal,2009,44(3):14-26.(in Chinese)

[2]陳詩一.邊際減排成本與中國環(huán)境稅改革.中國社會科學(xué),2011(3):85-100.

CHEN Shiyi.Marginal abatement cost and environmental tax reform in China.SocialSciencesinChina,2011(3):85-100.(in Chinese)

[3]DENISON E F.Accounting for slower economic growth:the United States in the 1970s.SouthernEconomicJournal,1981,47(4):1191-1193.

[4]GOLLOP F M,ROBERTS M J.Environmental regulations and productivity growth:the case of Fossil-fueled electric power generation.JournalofPoliticalEconomy,1983,91(4):654-674.

[5]GRAY W B.The cost of regulation:OSHA,EPA and the productivity slowdown.TheAmericanEconomicReview,1987,77(5):998-1006.

[6]PORTER M E.America′s green strategy.ScientificAmercian,1991,264(4):168.

[7]PORTER M E,VAN DER LINDE C.Green and competitive:ending the stalemate.HarvardBusinessReview,1995,73(5):120-134.

[8]張成,陸旸,郭路,等.環(huán)境規(guī)制強度與生產(chǎn)技術(shù)進步.經(jīng)濟研究,2011,46(2):113-124.

ZHANG Cheng,LU Yang,GUO Lu,et al.The intensity of environmental regulation and technological progress of production.EconomicsResearchJournal,2011,46(2):113-124.(in Chinese)

[10] 梁偉,朱孔來,姜巍.環(huán)境稅的區(qū)域節(jié)能減排效果及經(jīng)濟影響分析.財經(jīng)研究,2014,40(1):40-49.

LIANG Wei,ZHU Konglai,JIANG Wei.The regional energy conservation and emission reduction effect and economic impact of environmental tax.JournalofFinanceandEconomics,2014,40(1):40-49.(in Chinese)

[11] OUESLATI W.Environmental tax reform:short-term versus long-term macroeconomic effects.JournalofMacroeconomics,2014,40:190-201.

[12] ABDULLAH S,MORLEY B.Environmental taxes and economic growth:evidence from panel causality tests.EnergyEconomics,2014,42:27-33.

[13] ASHFORD N A,AYERS C,STONE R F.Using regulation to change the market for innovation.HarvardEnvironmentalLawReview,1985,9(2):419-466.

[14] HAHN R W.Economic prescriptions for environmental problems:how the patient followed the doctor′s orders.JournalofEconomicPerspectives,1989,3(2):95-114.

[15] 許士春,何正霞,龍如銀.環(huán)境規(guī)制對企業(yè)綠色技術(shù)創(chuàng)新的影響.科研管理,2012,33(6):67-74.

XU Shichun,HE Zhengxia,LONG Ruyin.The effects of environmental regulations on enterprise green technology innovation.ScienceResearchManagement,2012,33(6):67-74.(in Chinese)

[16] BERGEK A,BERGGREN C.The impact of environmental policy instruments on innovation:a review of energy and automotive industry studies.EcologicalEconomics,2014,106:112-123.

[17] VAN VLIET O,DE VRIES B,FAAIJ A,et al.Multi-agent simulation of adoption of alternative fuels.TransportationResearchPartD:TransportandEnvironment,2010,15(6):326-342.

[18] 李洪心,付伯穎.對環(huán)境稅的一般均衡分析與應(yīng)用模式探討.中國人口·資源與環(huán)境,2004,14(3):19-22.

LI Hongxin,FU Boying.General equilibrium analysis and application pattern inquiry to environmental taxation.ChinaPopulation,ResourcesandEnvironment,2004,14(3):19-22.(in Chinese)

[19] 劉曄,周志波.完全信息條件下寡占產(chǎn)品市場中的環(huán)境稅效應(yīng)研究.中國工業(yè)經(jīng)濟,2011(8):5-14.

LIU Ye,ZHOU Zhibo.Studies on effects of environmental taxes under oligopolistic markets with complete information.ChinaIndustrialEconomics,2011(8):5-14.(in Chinese)

[20] 何建武,李善同.節(jié)能減排的環(huán)境稅收政策影響分析.數(shù)量經(jīng)濟技術(shù)經(jīng)濟研究,2009,26(1):31-44.

HE Jianwu,LI Shantong.The impacts of energy tax and environmental tax on China′s economy.TheJournalofQuantitative&TechnicalEconomics,2009,26(1):31-44.(in Chinese)

[21] LEE C F,LIN S J,LEWIS C,et al.Effects of carbon taxes on different industries by fuzzy goal programming:a case study of the petrochemical related industries,Taiwan.EnergyPolicy,2007,35(8):4051-4058.

[22] LAZER D,PENTLAND A,ADAMIC L,et al.Computational social science.Science,2009,323(5915):721-723.

[23] 盛昭瀚,張維.管理科學(xué)研究中的計算實驗方法.管理科學(xué)學(xué)報,2011,14(5):1-10.

SHENG Zhaohan,ZHANG Wei.Computational experiments in management science and research.JournalofManagementSciencesinChina,2011,14(5):1-10.(in Chinese)

[24] 李悅雷,張維,熊熊.最小報價單位對市場流動性影響的計算實驗研究.管理科學(xué),2012,25(1):92-98.

LI Yuelei,ZHANG Wei,XIONG Xiong.Impact of tick size on market liquidity by agent-based modeling approach.JournalofManagementScience,2012,25(1):92-98.

[25] RIXEN M,WEIGAND J.Agent-based simulation of policy induced diffusion of smart meters.TechnologicalForecasting&SocialChange,2014,85:153-167.

[26] 吳術(shù),李心丹,張兵.基于計算實驗的賣空交易對股票市場的影響研究.管理科學(xué),2013,26(4):70-78.

WU Shu,LI Xindan,ZHANG Bing.Impacts of short sale on stock market based on computational experiment.JournalofManagementScience,2013,26(4):70-78.

[27] LOPOLITO A,MORONE P,TAYLOR R.Emerging innovation niches:an agent based model.ResearchPolicy,2013,42(6/7):1225-1238.

[28] 孫曉華,郭少蓉.新技術(shù)、異質(zhì)性偏好與產(chǎn)業(yè)演化.系統(tǒng)工程學(xué)報,2014,29(3):334-342.

SUN Xiaohua,GUO Shaorong.Heterogeneity preferences,new technology and evolution of industry.JournalofSystemsEngineering,2014,29(3):334-342.(in Chinese)

[29] 趙愛武,杜建國,關(guān)洪軍.基于計算實驗的有限理性消費者綠色購買行為.系統(tǒng)工程理論與實踐,2015,35(1):95-102.

ZHAO Aiwu,DU Jianguo,GUAN Hongjun.Green purchase behavior of bounded rational consumers based on computational experiments.SystemsEngineering-Theory&Practice,2015,35(1):95-102.(in Chinese)

[30] AFAROUI N,BROUILLAT E,JEAN M S.Policy design and technological substitution:investigating the REACH regulation in an agent-based model.EcologicalEconomics,2014,107:347-365.

[31] 劉小峰,盛昭瀚,杜建國.產(chǎn)品競爭與顧客選擇下的清潔生產(chǎn)技術(shù)演化模型.管理科學(xué),2013,26(6):25-34.

LIU Xiaofeng,SHENG Zhaohan,DU Jianguo.Evolutionary model of cleaner production technologies under product competition and consumer choice.JournalofManagementScience,2013,26(6):25-34.(in Chinese)

[32] ZEPPINI P,FRENKEN K,KUPERS R.Thresholds models of technological transitions.EnvironmentalInnovationandSocietalTransitions,2014,11:54-70.

SimulationofEnterpriseEnvironmentalTechnologyInnovationGiventheContextofEnvironmentalTax

ZHAO Aiwu1，DU Jianguo1，GUAN Hongjun2

1 School of Management, Jiangsu University, Zhenjiang 212013, China 2 School of Management Science and Engineering, Shandong University of Finance and Economics, Ji′nan 250014, China

Enterprise environmental technology innovation is the key to solve the problem of sustainable economic development and environmental sustainability. In order to improve both environmental performance and economical performance, environmental tax should be designed to induce enterprise environmental technology innovation.

With regard to the incentive effect of environmental tax on enterprise environmental technology innovation, scholars have made different conclusions based on different scenarios hypothesis, research samples, analytic methods and experimental models. On the one hand, environmental tax can stimulate enterprises to actively carry out environmental technology innovation. By the innovation compensation mechanism, enterprises can improve both economic performance and environmental performance. On the other hand, environmental tax will increase the cost of enterprises and weaken their innovation abilities. The effects of environmental tax on enterprise environmental technology innovation is a dynamic, evolutionary process. It is difficult to fully understand the influence mechanism of environmental tax on enterprise environmental technology innovation behavior by simply relying on empirical or mathematical methods.

This research is based on computational experimental method in social science and takes into account the complex and changing market to study the environmental technology innovation behavior. First of all, key factors and their transmission mechanism are identified and verified. Then, individual characteristics, behaviors and rules are described as the “function” of individual and its circumstance using mathematical methods. Further, a dynamic enterprise environmental technology innovation adoption model is built and the corresponding computer system is employed to simulate the evolution of enterprise technology innovation given the context of different customer choice, product competition and environmental tax scenarios. The model is employed to explore the impacts of timing and intensity of environmental tax on enterprise environmental performance and economic benefit. The experimental results show that the synergistic effects of market mechanism, enterprise competition and environmental tax can promote enterprises to adopt environmental innovational technology when the tax rate is high enough and the alternative environmental technology innovation reaches a certain level. This would reduce emissions yet improve economic benefits, and promote the industry to become green and sustainable. Meanwhile, the timing of the implementation of environmental tax policy plays a key role in the evolution of an industry to environmental technology innovation. It is related to the development level of environmental innovational technology. On the contrary, when the environmental tax is quite low, enterprises may “buy” pollution, though this is not the original intention of environmental tax.

This model can be used as a reference tool to test environmental tax by changing relative simulation parameters, and to observe the evolution of enterprise environmental technology innovation behaviors and its performance.

environmental technology innovation；environmental tax；computational experiment；economical performance；environmental performance

Date:March 24th, 2015

DateNovember 18th, 2015

FundedProject:Supported by the National Natural Science Foundation of China(71171099,71471076,71201071), Joint Research of the NSFC-NRF Scientific Cooperation Program(71411140250), the Humanities and Social Science Planning Foundation of Education Ministry(14YJAZH025) and the Shandong Natural Science Foundation(ZR2013GM003)

Biography:ZHAO Aiwu, doctor in engineering, currently works as a lecturer in the School of Management at Jiangsu University. Her research interests include modeling of complex system and computational experiment in social science. Her research works have been published in theSystemEngineering-Theory&Practice,ChineseJournalofManagementScience, andJournalofIntelligent&FuzzySystems. E-mail：aiwuzh@126.com

DU Jianguo, doctor in management, is a professor in the School of Management at Jiangsu University. His research interests include theory and methodology of complex system. His research works have been published in theJournalofManagementScienceinChina,SystemEngineering-Theory&Practice, andJournalofSystem&Management. E-mail：jgdu2005@126.com

GUAN Hongjun, doctor in engineering, is a professor in the School of Management Science and Engineering at Shandong University of Finance and Economics. His research interests include theory and methodology of management science, e-commerce. His research works have been published in theSystemEngineering-Theory&Practice,ChineseJournalofManagementScience, andJournalofIntelligent&FuzzySystems. E-mail：jjxyghj@126.com

F273.1

A

10.3969/j.issn.1672-0334.2016.01.004

1672-0334(2016)01-0040-13

2015-03-24修返日期2015-11-18

國家自然科學(xué)基金(71171099,71471076,71201071)；國家自然科學(xué)基金中韓國際合作與交流項目(71411140250)；教育部人文社會科學(xué)研究規(guī)劃項目(14YJAZH025)；山東省自然科學(xué)基金(ZR2013GM003)

趙愛武，工學(xué)博士，江蘇大學(xué)管理學(xué)院講師，研究方向為復(fù)雜系統(tǒng)建模與社會科學(xué)計算實驗，曾在《系統(tǒng)工程理論與實踐》《中國管理科學(xué)》和《Journal of Intelligent & Fuzzy Systems》等雜志發(fā)表論文，E-mail：aiwuzh@126.com

杜建國，管理學(xué)博士，江蘇大學(xué)管理學(xué)院教授，研究方向為復(fù)雜系統(tǒng)理論與方法，曾在《管理科學(xué)學(xué)報》《系統(tǒng)工程理論與實踐》和《系統(tǒng)管理學(xué)報》等雜志發(fā)表論文，E-mail：djg@ujs. edu.cn

關(guān)洪軍，工學(xué)博士，山東財經(jīng)大學(xué)管理科學(xué)與工程學(xué)院教授，研究方向為管理科學(xué)理論與方法、電子商務(wù)等，曾在《系統(tǒng)工程理論與實踐》《中國管理科學(xué)》和《Journal of Intelligent & Fuzzy Systems》等雜志發(fā)表論文，E-mail: jjxyghj@126.com

□