考慮技術(shù)投資的供應(yīng)鏈網(wǎng)絡(luò)下可交易污染排放許可均衡問題

楊玉香，吳增源，黃祖慶

(中國計量大學(xué)經(jīng)濟(jì)與管理學(xué)院，浙江 杭州 310018)

?

考慮技術(shù)投資的供應(yīng)鏈網(wǎng)絡(luò)下可交易污染排放許可均衡問題

楊玉香，吳增源，黃祖慶

(中國計量大學(xué)經(jīng)濟(jì)與管理學(xué)院，浙江 杭州 310018)

構(gòu)建了包含多污染物、多污染接收點的供應(yīng)鏈網(wǎng)絡(luò)，考慮了生產(chǎn)技術(shù)和減排技術(shù)投資，對可交易的污染排放許可下的供應(yīng)鏈網(wǎng)絡(luò)進(jìn)行了均衡分析，給出了變分不等式框架，建立了均衡模型。模型明確地處理了空間差異，同時保證了通過初始污染排放許可分配就可達(dá)到指定的環(huán)境標(biāo)準(zhǔn)。提出了模型的求解算法，最后利用提出的模型及算法對算例進(jìn)行了計算，得到利潤最大化下供應(yīng)鏈網(wǎng)絡(luò)中各工廠的產(chǎn)品生產(chǎn)量、排放量、污染排放許可的均衡分配額和污染排放許可的交易價格、技術(shù)最優(yōu)投資額，并分析了成本函數(shù)參數(shù)以及工廠數(shù)量的變化對生產(chǎn)技術(shù)投資和減排技術(shù)投資的影響。

供應(yīng)鏈網(wǎng)絡(luò)；可交易污染排放許可；變分不等式；技術(shù)投資

1 引言

低碳時代，在供應(yīng)鏈間競爭的背景下，企業(yè)應(yīng)該從供應(yīng)鏈角度出發(fā)，優(yōu)化各節(jié)點企業(yè)及整個供應(yīng)鏈的低碳行為，通過鏈上企業(yè)相互合作來解決碳排放問題。然而，目前，我國企業(yè)并未足夠重視低碳化發(fā)展，在開發(fā)應(yīng)用各項節(jié)能減排的新技術(shù)上缺乏動力。為激勵企業(yè)使用低碳技術(shù)，推動節(jié)能減排，可利用市場機(jī)制的調(diào)節(jié)作用，即制定相關(guān)的經(jīng)濟(jì)激勵政策[1]。張保銀等[2]基于委托-代理模型分析了政府對企業(yè)的最優(yōu)激勵合約的制定問題。鄧峰[3]構(gòu)建政府是否選擇污染排放檢查與企業(yè)是否服從管制的靜態(tài)博弈模型并提出相關(guān)對策。常香云等[4]為引導(dǎo)企業(yè)選擇低碳制造/再制造技術(shù)，針對碳稅、補(bǔ)貼政策，建立系統(tǒng)動力學(xué)模型，分析了以上政策對企業(yè)決策的影響。張國興等[5]分析了政府與企業(yè)在節(jié)能減排補(bǔ)貼申請與發(fā)放過程中雙方策略的選擇與依存性變化，探討補(bǔ)貼額、違規(guī)成本的最優(yōu)邊界問題。然而，除政府補(bǔ)貼和稅收政策之外，污染排放許可交易機(jī)制也是一種有效的經(jīng)濟(jì)杠桿，能有效實現(xiàn)減少污染排放的目標(biāo)。

目前，對于污染排放許可的分配和定價問題研究較多。首先，關(guān)于初始污染排放許可分配的研究。Montgomery[6]提出完全競爭的排放權(quán)市場，在總量一定條件下，排放權(quán)的最終配置與初始分配是互相獨(dú)立的，政府無需知道各個污染源的成本函數(shù)，只需根據(jù)環(huán)境容量確定排污總量，市場最終能實現(xiàn)均衡。而考慮不完全競爭市場、企業(yè)在產(chǎn)品市場的勢力以及排放權(quán)交易成本，則排污交易的效率將受到初始排放權(quán)分配的影響。徐玉高等[7]在不考慮和考慮各國歷史排放兩種情況下分別提出碳排放權(quán)分配機(jī)制。陳文穎和吳宗鑫[8]提出兩種綜合考慮公平和效率的碳權(quán)混合分配機(jī)制，在確定碳權(quán)交易價格的基礎(chǔ)上，對全球碳權(quán)交易情況進(jìn)行模擬，分析碳排放權(quán)交易對中國經(jīng)濟(jì)的影響以及不同碳權(quán)分配機(jī)制對全球碳權(quán)交易收益的影響。Jensen和Rasmussen[9]利用均衡模型研究排放權(quán)分配，發(fā)現(xiàn)利用排放權(quán)的盈利來降低現(xiàn)行稅會最大限度地降低福利成本，但也會減少能源密集部門的就業(yè)機(jī)會，根據(jù)市場份額分配排放權(quán)會減少部門調(diào)整程度，導(dǎo)致高的福利成本。Malik[10]考慮工廠具有市場支配力且存在一個或多個不服從工廠的排放權(quán)交易市場，結(jié)論：如果所有工廠都順從，則具有市場支配力的工廠可能選擇擁有超過其所需的排放權(quán)數(shù)量，并使剩余排放權(quán)從交易市場上退出。Sartzetakis[11]研究當(dāng)生產(chǎn)成本與削減成本不一致時，在不完全競爭的產(chǎn)品市場上，計劃分配優(yōu)于市場配置的條件。Park等[12]提出使用玻爾茲曼分布分配排放權(quán)，這是一種全新的分配方法。其次，關(guān)于污染排放許可的定價問題。排放權(quán)首先是一種環(huán)境容量資產(chǎn)，對其定價遵循了環(huán)境資源定價的方法。Rubin[13]使用均衡模型探討在連續(xù)且有限時間前景下最小化公司聯(lián)合成本，同時排放權(quán)可以借貸條件下的定價問題，表明減排邊際貼現(xiàn)成本應(yīng)不隨時間變動。由于均衡條件下，排放權(quán)應(yīng)該等于邊際成本，因此許可價格應(yīng)隨無風(fēng)險貼現(xiàn)率增長。Cronshaw等[14]擴(kuò)展了Rubin的工作，假定公司目標(biāo)是最小化期望總成本，分析了被管制者的排放權(quán)存儲行為，通過一個隨機(jī)的、連續(xù)且有限時間下的模型研究內(nèi)涵現(xiàn)貨價格，其結(jié)論是預(yù)期的排放權(quán)價格將以一個不高于貼現(xiàn)率的速度增加。Paolella等[15]和Benz和Trueck[16]研究了歐洲ETS中排放權(quán)的動態(tài)價格機(jī)制，表明處罰成本、交易的時間間隔、遞增的邊際減排成本等對價格造成了影響。Phelan等[17]建立碳排放與其產(chǎn)生損失之間的責(zé)任聯(lián)系，利用保險體系為碳排放權(quán)定價。此外，Jotzo和Pezzey[18]認(rèn)為未來經(jīng)濟(jì)發(fā)展路徑與排放量增長的不確定性對總量控制以及碳排放權(quán)交易效率是一個障礙，因而達(dá)成原定減排目標(biāo)的成本是不確定的，而確定最優(yōu)排放強(qiáng)度目標(biāo)可使碳排放權(quán)交易更有效。Li Shoude和Gu Mengdi[19]研究可交易的排放權(quán)對企業(yè)生產(chǎn)和庫存策略的影響，并將其引入Arrow-Karlin動態(tài)生產(chǎn)庫存模型，分析排放權(quán)交易對企業(yè)生產(chǎn)庫存決策的影響。Chin和Zhang Peng[20]應(yīng)用Cournot模型模擬航空業(yè)排放量并分析歐盟排放權(quán)交易系統(tǒng)和加強(qiáng)型歐盟排放權(quán)交易系統(tǒng)的效率，得到前者的排放量比后者要大同時效率更低。

綜上所述，以上文獻(xiàn)主要集中于從宏觀領(lǐng)域研究污染排放許可交易問題，關(guān)于污染排放交易政策對不同減排主體的影響問題，僅限于研究此激勵政策對單個企業(yè)的影響，值得關(guān)注的是，鮮有文獻(xiàn)基于供應(yīng)鏈背景研究此問題。因此，本文將分析污染排放交易系統(tǒng)對供應(yīng)鏈網(wǎng)絡(luò)中各決策者決策行為的影響，在供應(yīng)鏈網(wǎng)絡(luò)背景下構(gòu)建企業(yè)污染排放許可均衡模型，在企業(yè)服從規(guī)定的情況下研究政府如何對污染排放許可進(jìn)行初始分配以達(dá)到指定的環(huán)境標(biāo)準(zhǔn)。

2 模型的建立

2.1 模型描述

考慮包括M個工廠和N個需求市場的供應(yīng)鏈網(wǎng)絡(luò)，每個特定的工廠記為m，每個特定的需求市場記為n；排出T種污染物，每種特定的污染物記為t；有J個污染接收點，每個特定的接收點記為j。假設(shè)每個工廠生產(chǎn)的產(chǎn)品均為同質(zhì)產(chǎn)品，不同工廠生產(chǎn)的產(chǎn)品沒有區(qū)別。工廠從長遠(yuǎn)角度為了保持競爭力需考慮投入資源提高生產(chǎn)技術(shù)和減排技術(shù)，以提高生產(chǎn)能力同時減少污染排放。工廠在每個特定的污染接收點只有持有污染排放許可，才可以排放污染，因此，對于每個工廠，都必須在所有的接收點持有一定份額的污染許可。同時，假設(shè)污染許可市場是完全競爭的。

為方便建模，給出如下參量和決策變量：

2.2 工廠的優(yōu)化條件

供應(yīng)鏈網(wǎng)絡(luò)內(nèi)的每個工廠都必須考慮各種成本，比如生產(chǎn)成本、減排成本、污染許可購買成本、生產(chǎn)和減排技術(shù)采用成本。

(1)

s.t.

(2)

(3)

假設(shè)此供應(yīng)鏈網(wǎng)絡(luò)中的M個寡頭壟斷工廠在產(chǎn)品市場中為非合作競爭關(guān)系，此網(wǎng)絡(luò)的均衡為Nash-Cournot均衡，先給出如下定義。

(4)

2.3 需求市場的優(yōu)化條件

各工廠生產(chǎn)的產(chǎn)品在需求市場n的單位需求價格為p3n，N個市場的價格構(gòu)成向量p3，市場的需求記為dn=dn(p3)。則市場n的優(yōu)化條件即空間價格均衡條件可表示為：

(5)

(6)

(7)

2.4 污染許可的市場出清條件

控制污染許可市場出清的經(jīng)濟(jì)系統(tǒng)條件可表示為：對每一污染物t和每一污染接收點j:

(8)

式(8)描述的是：如果在污染接收點j處關(guān)于污染物t的污染許可價格是正的，那么在此處的污染許可市場必須出清；如果在某一接收點j某一特定污染物t的污染許可供應(yīng)過剩，那么污染許可的價格必須是零。

(9)

2.5 網(wǎng)絡(luò)均衡模型

(10)

證明：變分不等式(4)適用所有的工廠，變分不等式(7)適用于所有的市場，變分不等式(9)適用于所有污染物和所有污染接收點，將上述變分不等式相加即得(10)。

反過來，下面證明(10)的解也滿足變分不等式(4)、(7)和(8)。

(11)

(12)

式(12)即意味著式(4)和(7)的成立。

下面證明如果對于每種污染物和每個污染接收點的污染許可是固定的，那么供應(yīng)鏈網(wǎng)絡(luò)均衡模式獨(dú)立于污染許可的初始分配，同時為了保證達(dá)到指定的環(huán)境質(zhì)量標(biāo)準(zhǔn)，污染許可的總量應(yīng)該等于指定的標(biāo)準(zhǔn)。

推論1顯示污染許可的任一初始分配依賴于每一污染物和每一接收點指定的標(biāo)準(zhǔn)，這不影響均衡模式。

下面提出一種機(jī)制決定初始分配的總和可以滿足環(huán)境質(zhì)量標(biāo)準(zhǔn)。

證明：根據(jù)約束(2)，有：

(13)

另外，由均衡條件(9)，得出：

(14)

得證。

3 模型求解

本文利用MPM(Modified Projection Method，修正投影法)求解污染許可的市場均衡模型即變分不等式(10)，算法實現(xiàn)如下：

第二步：計算。

(15)

第三步：修正。

(16)

第四步：收斂性驗證。若

則停止，其中ε>0，為一指定的收斂標(biāo)準(zhǔn)；否則τ=τ+1，返回第二步。

4 算例

本部分通過算例來進(jìn)一步解釋第2部分中提出的模型，同時驗證第3部分提出的算法的性能。

4.1 算例1.

研究的供應(yīng)鏈網(wǎng)絡(luò)包括2個工廠、2個需求市場，每個工廠排放2種污染物，影響2個污染接收點。算例中將使用如下二次生產(chǎn)成本函數(shù)和排放成本函數(shù)。

每個工廠的生產(chǎn)成本函數(shù)為：

工廠m將產(chǎn)品運(yùn)往市場n支付的運(yùn)輸成本為kmn(xmn)=2(xmn)2。

每個工廠排放污染物支付的成本為：

表1 不同下迭代次數(shù)的變化情況

需求市場的需求函數(shù)為：d1(p3)=-2p31-1.5p32+500,d2(p3)=-p32-0.5p31+300。

根據(jù)第3部分提出的MPM算法，編寫基于Matlab的應(yīng)用程序，且由MPM可知，迭代次數(shù)τ依賴于收斂標(biāo)準(zhǔn)ε，表1給出了在參數(shù)α=0.02時不同收斂標(biāo)準(zhǔn)ε下迭代次數(shù)τ的變化情況。

設(shè)置收斂標(biāo)準(zhǔn)ε=10-4，參數(shù)α=0.02，運(yùn)行程序得到如下均衡結(jié)果：

工廠的污染許可均衡值分別為：

工廠的均衡邊際減排成本分別為：

影響2個接收點的污染許可均衡價格分別為：

接下來，進(jìn)一步調(diào)整工廠的生產(chǎn)成本函數(shù)和排放污染成本函數(shù)，分別為：

每個工廠的生產(chǎn)成本函數(shù)為：

每個工廠排放污染物支付的成本為：

4.2 算例2

供應(yīng)鏈網(wǎng)絡(luò)中的工廠增加為3個，仍然有2個需求市場，每個工廠排放2種污染物，影響2個污染接收點。

工廠3的生產(chǎn)成本函數(shù)為：

工廠3排放污染物支付的成本為：

網(wǎng)絡(luò)中的成員數(shù)量增加后，同樣根據(jù)第3部分提出的MPM算法，編寫基于Matlab的應(yīng)用程序，表2給出了在參數(shù)α=0.02時不同收斂標(biāo)準(zhǔn)ε下迭代次數(shù)τ的變化情況。

設(shè)置收斂標(biāo)準(zhǔn)ε=10-4，參數(shù)α=0.02，運(yùn)行程序得到如下均衡結(jié)果：

工廠的污染許可均衡值分別為：

工廠的均衡邊際減排成本分別為：

影響2個接收點的污染許可均衡價格分別為：

表2 不同下迭代次數(shù)的變化情況

可見，由于污染許可市場上參與者增加了，導(dǎo)致為提高生產(chǎn)技術(shù)和減排技術(shù)的總投資也提高了。

接下來，進(jìn)一步調(diào)整工廠的生產(chǎn)成本函數(shù)和排放污染成本函數(shù)，工廠1和2的成本函數(shù)與算例1相同，工廠3的相關(guān)成本函數(shù)調(diào)整為：

調(diào)整后的工廠3的生產(chǎn)成本函數(shù)為：

調(diào)整后的工廠3排放污染物支付的成本為：

綜上所述，通過上面的算例分析了成本函數(shù)參數(shù)的變化以及供應(yīng)鏈網(wǎng)絡(luò)中工廠數(shù)量對生產(chǎn)技術(shù)投資和減排技術(shù)投資的影響，可以發(fā)現(xiàn)，要降低生產(chǎn)成本和減排成本，勢必要增加生產(chǎn)技術(shù)和減排技術(shù)的投資。此外，生產(chǎn)技術(shù)和減排技術(shù)的投資與污染許可市場中參與的工廠數(shù)量密切相關(guān)。

5 結(jié)語

本文建立了多污染物、多污染接收點的供應(yīng)鏈污染排放許可變分不等式框架，給出了均衡模型，對污染許可和技術(shù)投資進(jìn)行分析及均衡計算，特別地，在模型中考慮了兩種技術(shù)：生產(chǎn)技術(shù)和減排技術(shù)。同時，提出了均衡模型的求解算法，最后通過算例驗證了模型及求解算法的有效性，并對參數(shù)的變化及工廠數(shù)量對生產(chǎn)技術(shù)投資和減排技術(shù)投資的影響進(jìn)行了分析。進(jìn)一步的研究可以將此模型擴(kuò)展到動態(tài)模型分析。

致謝：

浙江省哲社重點研究基地“產(chǎn)業(yè)發(fā)展政策研究中心”和浙江省人文社科基地“管理科學(xué)與工程”資助.

[1] 許士春, 何正霞, 龍如銀. 環(huán)境政策工具比較：基于企業(yè)減排的視角[J]. 系統(tǒng)工程理論與實踐, 2012, 32(11): 2351-2362.

[2] 張保銀, 汪波，吳煜. 基于循環(huán)經(jīng)濟(jì)模式的政府激勵與監(jiān)督問題[J]. 中國管理科學(xué), 2006, 14(1):136-141.

[3] 鄧峰. 基于不完全執(zhí)行污染排放管制的企業(yè)與政府博弈分析[J]. 預(yù)測, 2008, 27(1): 67-71.

[4] 常香云, 鐘永光, 王藝璇, 等. 促進(jìn)我國汽車零部件再制造的政府低碳引導(dǎo)政策研究——以汽車發(fā)動機(jī)再制造為例[J]. 系統(tǒng)工程理論與實踐, 2013, 33(11): 2811-2821.

[5] 張國興, 張緒濤, 汪應(yīng)洛, 等. 節(jié)能減排政府補(bǔ)貼的最優(yōu)邊界問題研究[J]. 管理科學(xué)學(xué)報, 2014, 17(11): 129-138.

[6] Montgomery D. Markets in licenses and efficient pollution control programs[J]. Journal of Economic Theory, 1972, 5(3): 395-418.

[7] 徐玉高, 郭元, 吳宗鑫. 碳權(quán)分配：全球碳排放權(quán)交易及參與激勵[J]. 數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究, 1997, (3): 72-77.

[8] 陳文穎, 吳宗鑫. 碳排放權(quán)分配與碳排放權(quán)交易[J]. 清華大學(xué)學(xué)報(自然科學(xué)版), 1998, 38(12): 15-18.

[9] Jensen J, Rasmussen TN. Allocation of CO2 emissions permits: a general equilibrium analysis of policy instruments[J]. Journal of Environmental Economics and Management, 2000, 40(2): 111-136.

[10] Malik A. Further results on permit markets with market power and cheating[J]. Journal of Environmental Economics and Management, 2002, 44(3): 371-390.

[11] Sartzetakis E. On the efficiency of competitive markets for emission permits[J]. Environmental and Resource Economics, 2004, 27(1): 1-19.

[12] Park J, Kin CU, Isard W. Permit allocation in emissions trading using the Boltzmann distribution[J]. Physica A, 2012, 391(20): 4883-4890.

[13] Rubin J. A model of intertemporal emission trading, banking, and borrowing[J]. Journal of Environmental Economics and Management, 1996, 31(3): 269-286.

[14] Cronshaw B, Kruse J, Schennach S. The economics of pollution permit banking in the context of title IV of the 1990 clean air act amendments[J]. Journal of Environmental Economics and Management, 2000, 40(3): 189-210.

[15] Paolella M, Taschini L. An econometric analysis of emission trading allowances[J]. Journal of Banking and Finance, 2008, 32(10): 2022-2032.

[16] Benz E, Trueck S. Modeling the price dynamics of CO2 emission allowances[J]. Energy Economics, 2009, 31(1):4-15.

[17] Phelan L, Henderson-Sellers A, Taplin R. Climate change, carbon prices and insurance systems[J]. International Journal of Sustainable Development and World Ecology, 2010, 17(2): 95-108.

[18] Jotzo F, Pezzey JCV. Optimal intensity targets for greenhouse gas emissions trading under uncertainty[J]. Environmental and Resource Economics, 2007, 38(2): 259-284.

[19] Li Shoude, Gu Mengdi. The effect of emission permit trading with banking on firm’s production-inventory strategies[J]. International Journal of Production Economics, 2012, 137(2): 304-308.

[20] Chin ATH, Zhang Peng. Carbon emission allocation methods for the aviation sector[J]. Journal of Air Transport Management, 2013,28:70-76.

Marketable Pollution Permits Equilibrium Problem in Supply Chain Network with Technological Investment

YANG Yu-xiang, WU Zeng-yuan, Huang Zu-qing

(China Jiliang University, College of Economics and Management, Hangzhou 310018, China)

Low carbon development for supply chain is an important way to seek the development mode of low carbon economy. In this context, firms should optimize low carbon behavior of every node firm and the whole supply chain from the view of supply chain, to solve carbon emission problem through cooperation of firms in the supply chain, and to implement low carbon supply chain management. Currently, the theoretical research and legislative practice on low carbon economy are still in their initial stages in China. Many firms did not still pay enough attention to low carbon development, who are lack of power to apply the new technology of the energy saving and emission reduction. In order to encourage enterprises to use low carbon technology and promote the energy saving and emission reduction, the government can formulate economy prompting policies. Among these policies, the marketable pollution permits is an effective economic lever, which can realize to decrease pollution emission. In this paper, the effect of the marketable pollution permits on the decision behavior of each decision-maker is analyzed in the supply chain network, and an equilibrium model is proposed in supply chain context, to study how government should allocate initially the pollution permits to obtain assigned environmental goal. Hence, a supply chain network with multi-pollutant and multi receptor points is developed. The investment in production technology and emission-abasement technology is included. The equilibrium patterns for the supply chain network with marketable pollution permits are analyzed. The variational inequality framework is given. The equilibrium model is developed. The model deals explicitly with spatial differentiation and also guarantees that the imposed environmental standards are met through the initial allocation of licenses. A solution algorithm for the model is proposed. Finally, numerical examples are solved by using the proposed model and algorithm. The profit-maximized quantities of oligopolistic firms’ products and the quantities of emission, along with the equilibrium allocation of licenses and their prices, as well as the optimal investments in the technologies for the supply chain network are computed. Moreover, the effect of changes in investment cost parameters, the number of firms on the investment measures for production and emission-abatement technology are illustrated. The study can guide firms in the supply chain network to join the ranks of low carbon and implement production mode of low carbon economy, and provide decision-making basis for relevant departments of government to design the pollution permits system.

supply chain network; marketable pollution permits; variational inequality; technology investment

1003-207(2016)04-0074-09

10.16381/j.cnki.issn1003-207x.2016.04.009

2014-05-30；

2015-04-07

國家自然科學(xué)基金青年資助項目(71402173);教育部人文社會科學(xué)基金青年資助項目(14YJC630168)

楊玉香(1979-)，女(漢族)，吉林松原人，中國計量大學(xué)經(jīng)濟(jì)與管理學(xué)院博士，副教授，研究方向：供應(yīng)鏈管理，E-mail：yyx_bj2005@126.com.

F224

A