魯祖亮 李林 曹龍舟 黃飛

摘?要?利用擬合有限體積法研究京津地區(qū)跨界空氣污染問題，提出碳排放交易條件下的隨機微分博弈模型，利用隨機最優(yōu)控制理論推導(dǎo)出合作博弈模型的哈密頓-雅可比-貝爾曼方程，并利用京津地區(qū)數(shù)據(jù)實證分析了算法的實用性和有效性.

關(guān)鍵詞?跨界空氣污染問題;京津地區(qū);哈密頓-雅可比-貝爾曼方程;擬合有限體積法

中圖分類號?O242.1?文獻(xiàn)標(biāo)識碼?A

Abstract?We study transboundary air pollution of Beijing-Tianjin Region by using the fitted finite volume method. We present a stochastic differential game to model the problems with carbon emission permits trading. By using the stochastic optimal control theory， we ?derive the system of Hamilton-Jacobi-Bellman equations for the cooperative games. Finally， we give an empirical study to illustrate the efficiency and the usefulness of this method based ?on the data of Beijing-Tianjin region.

Key words?transboundary air pollution; Beijing-Tianjin region; Hamilton-Jacobi-Bellman equation; fitting finite volume method

1?引?言

隨著社會經(jīng)濟(jì)的迅速發(fā)展，二氧化碳污染和碳排放問題成為我國重要污染問題之一.二氧化碳的污染是指人類生產(chǎn)生活向大氣排出污染物，其含量超過環(huán)境承受能力，使大氣質(zhì)量發(fā)生惡化，使人們的生活、工作、健康以及生態(tài)環(huán)境等遭受惡劣響.京津地區(qū)的大氣污染問題愈發(fā)嚴(yán)重，尤其是二氧化碳的污染問題，故對京津地區(qū)二氧化碳的跨界污染問題研究是十分有必要的.

在京津地區(qū)，北京大氣環(huán)境中大約23%的二氧化碳是來自北京周邊地區(qū)，如天津、秦皇島、唐山和廊坊一帶.根據(jù)中國科學(xué)院的研究結(jié)果：北京市二氧化碳排放大約30%～40%來自原始排放，30%～40%來自周邊區(qū)域輸送，北京280 m高空的二氧化碳濃度遠(yuǎn)遠(yuǎn)高于8 m的濃度.北京城區(qū)大量污染物通過低空氣團(tuán)自東南方向近距離輸送至北京，與當(dāng)?shù)氐奈廴練怏w混合后加重污染.京津地區(qū)的二氧化碳跨界污染造成的影響如此之大，如果無法區(qū)分跨界污染的主體責(zé)任，就可能造成管理責(zé)任不清，當(dāng)?shù)卣疅o法對本區(qū)域內(nèi)的大氣質(zhì)量負(fù)責(zé).北京市和天津市作為中國的兩個直轄市、超大城市和經(jīng)濟(jì)文化中心，每年的生產(chǎn)總值在全國范圍內(nèi)遙遙領(lǐng)先，一直領(lǐng)跑全國的經(jīng)濟(jì)水平，但京津地區(qū)所產(chǎn)生的二氧化碳污染也是全國最嚴(yán)峻的，霾污染事件的頻繁發(fā)生為我國環(huán)境危機拉響了警報，解決經(jīng)濟(jì)發(fā)展與大氣環(huán)境污染之間的矛盾勢在必行.跨界大氣污染造成的影響十分巨大，而一個城市的二氧化碳排放和經(jīng)濟(jì)的發(fā)展有著密不可分的關(guān)系，因此對京津地區(qū)跨界二氧化碳的污染問題進(jìn)行研究，通過數(shù)值計算提出有效建議，將有助于改善京津地區(qū)的環(huán)境，促進(jìn)京津區(qū)域的可持續(xù)發(fā)展.

2007年，Yeung率先推導(dǎo)出合作微分博弈中的時間均衡解，并且研究了隨機微分博弈框架下的污染管理問題[1].Youssef考慮了非合作且對稱的三階段兩公司博弈問題，結(jié)果表明投資研發(fā)的溢出效應(yīng)和公司之間的競爭可以促進(jìn)非合作國家更好地將跨界污染問題國際化[2].Benchekroun等人分析了非合作跨界污染博弈中，采用清潔技術(shù)所帶來的后果[3].Li在Yeung提出理論的基礎(chǔ)之上，考慮了排放權(quán)交易機制對各博弈者決策的影響[4].

3?實證分析

這一節(jié)，討論京津地區(qū)跨界二氧化碳污染問題的數(shù)值模擬.首先，給出北京市和天津市2006-2015年的經(jīng)濟(jì)數(shù)據(jù)，如表1所示：

其中GDP表示國內(nèi)生產(chǎn)總值，北京市和天津市的碳排放交易價格是20～50元，年碳排放總量是300-2000萬噸，再通過文獻(xiàn)[7-8]可以大致算出2006-2015年的北京市和天津市二氧化碳排放數(shù)據(jù)如表2：

在圖3中，給出了京津地區(qū)的總價值函數(shù)隨污染量和排放權(quán)價格的變化圖像，在圖4～5中，分別給出了北京市和天津市的二氧化碳排放量隨污染量和排放權(quán)價格的變化圖像.從這3個圖像可以看到，隨著碳排放價格增加產(chǎn)生的價值也會有所增加，當(dāng)價格增長至一定時刻后會有一定的下降趨勢，京津地區(qū)的排放權(quán)交易量會隨著碳排放價格的增加而減小.由此可見，圖5～6驗證了模型的有效性.

為了更好驗證價值函數(shù)的有效性，抽取了t=0時刻（碳排放量為400萬噸）與t=5時刻（碳排放量為400萬噸）的價值函數(shù)數(shù)據(jù)，如圖6所示，從圖像可知，當(dāng)排放權(quán)價格大于5時t=0時刻的價值函數(shù)會明顯小于t=5時刻的價值函數(shù)，這和理論分析是完全吻合的.

參考文獻(xiàn)

[1]?YEUNG D. Dynamically consistent cooperative solution in a differential game of transboundary industrial pollution[J]. Journal of Optimization Theory and Applications， 2007， 134（1）： 143-160.

[2] YOUSSEF S. Transboundary pollution， R&D spillovers and international trade[J]. The Annals of Regional Science， 2009， 43（1）： 235-250.

[3] BENCHEKROUN H ， CHAUDHURI A. Transboundary pollution and clean technologies[J]. Resource and Energy Economics， 2014， 36（2）：601-619.

[4] LI S. A differential game of transboundary industrial pollution with emission permits trading[J]. Journal of Optimization Theory and Applications，2014， 163（2）：642-659.

[5] LABRIET M， LOULOU R. In the framework of linear programming coupled climate destruction and greenhouse gas emission reduction costs[J]. Environmental modeling and evaluation， 2003， 8（2）： 261-274.

[6] HUANG C， WANG S， CHEN C， LI Z. A radial basis collocation method for Hamilton-Jacobi-Bellman [J]. 2006， 42（12）：2201-2207.

[7] 王卉彤， 慕淑茹. 北京市能源消費總量、結(jié)構(gòu)與碳排放的趨勢研究[J]. 城市發(fā)展研究， 2010， 17（9）：55-60.

[8] 王雅捷， 何永. 基于碳排放清單編制的低碳城市規(guī)劃技術(shù)方法研究[J]. 中國人口·資源與環(huán)境， 2015， 25（6）：72-80.