綠色環(huán)保當(dāng)量稅額的協(xié)同效應(yīng)研究

李小梅 喻宣瑞 常 進(jìn)

內(nèi)容提要：基于全國31 個?。▍^(qū)、市）的綠色環(huán)保當(dāng)量稅額、污染氣體排放量、地區(qū)生產(chǎn)總值的面板數(shù)據(jù)，運用Copula 聯(lián)合概率預(yù)測模型，研究了三者之間的協(xié)同效應(yīng)與互動機制。研究結(jié)果表明：綠色環(huán)保當(dāng)量稅額對經(jīng)濟和環(huán)境的改善作用，存在著地區(qū)差異，與各地區(qū)產(chǎn)業(yè)結(jié)構(gòu)相關(guān)；綠色環(huán)保當(dāng)量稅額的最優(yōu)波動區(qū)間為2.1 元至3.3 元。

一、引 言

改革開放以來，我國經(jīng)濟快速發(fā)展，人民的生活水平得到了大幅度提高，公共基礎(chǔ)設(shè)施和公共服務(wù)也有所改善，但資源的不合理利用，資源浪費而引發(fā)的環(huán)境問題也在一定程度上阻礙了經(jīng)濟的可持續(xù)發(fā)展。黨的十九屆六中全會明確指出，要立足新發(fā)展階段、貫徹新發(fā)展理念、構(gòu)建新發(fā)展格局、推動高質(zhì)量發(fā)展。這就要求社會經(jīng)濟既要實現(xiàn)“量”的增長，也要推進(jìn)“質(zhì)”的提升，才能著實提高人民的幸福指數(shù)。

實現(xiàn)社會經(jīng)濟質(zhì)與量的共同提升離不開國家政策的引導(dǎo)。綠色稅收作為政府調(diào)節(jié)社會經(jīng)濟發(fā)展和生態(tài)環(huán)境的一種強制性手段，主要用于平衡社會經(jīng)濟和生態(tài)環(huán)境的發(fā)展，是實現(xiàn)二者可持續(xù)發(fā)展的關(guān)鍵點（李彪，2020）。綠色稅收，這一概念最早由瑞典等發(fā)達(dá)國家提出（何平林等，2019）。其最初源于英國經(jīng)濟學(xué)家庇古所提出的負(fù)外部性理論，即根據(jù)污染所造成的危害程度對排污者進(jìn)行征稅，用于彌補排污者生產(chǎn)的私人成本和社會成本之間的差距①［英］庇古著，朱泱等譯：《福利經(jīng)濟學(xué)（上卷）》，商務(wù)印書館2006 年版，第207-208 頁。。為此我國吸收借鑒了這一方式，于2018 年正式實施了《中華人民共和國環(huán)境保護稅法》，通過對污染排放企業(yè)進(jìn)行征稅，擬實現(xiàn)對社會經(jīng)濟和生態(tài)環(huán)境的靈活調(diào)控。但我國綠色稅收體系構(gòu)建相對較晚，尚不健全，尤其對綠色環(huán)保當(dāng)量稅額的設(shè)置，還缺乏一套行之有效的理論。綠色環(huán)保當(dāng)量稅額過低，對環(huán)境污染無法起到應(yīng)有的規(guī)制作用，而過高的綠色環(huán)保當(dāng)量稅額會加重傳統(tǒng)企業(yè)的稅收負(fù)擔(dān)，不利于經(jīng)濟的發(fā)展。如何實現(xiàn)綠色環(huán)保當(dāng)量稅額對環(huán)境和經(jīng)濟的合理調(diào)控，是近年來學(xué)術(shù)界關(guān)注的熱點與難點問題。

確定合理的綠色環(huán)保當(dāng)量稅額，實現(xiàn)對社會經(jīng)濟和生態(tài)環(huán)境的靈活調(diào)控，是實現(xiàn)綠色發(fā)展的關(guān)鍵環(huán)節(jié)?；诖?，本文從綠色環(huán)保當(dāng)量稅額、社會經(jīng)濟發(fā)展、生態(tài)環(huán)境這三個維度出發(fā)，分析了三者的互動過程，揭示了三者的協(xié)同效應(yīng)。

二、文獻(xiàn)綜述

為探究綠色環(huán)保當(dāng)量稅額對經(jīng)濟與環(huán)境的影響機制，相關(guān)學(xué)者對其進(jìn)行了深入的研究。概括起來，主要包括以下兩個方面：一方面，綠色環(huán)保當(dāng)量稅額對經(jīng)濟與環(huán)境的單向作用。關(guān)于綠色環(huán)保當(dāng)量稅額對經(jīng)濟增長的作用機制，已有的研究認(rèn)為，綠色環(huán)保當(dāng)量稅額過高將會加重傳統(tǒng)行業(yè)的稅負(fù)，進(jìn)而影響經(jīng)濟的發(fā)展（Fang et al.，2013），但適當(dāng)提高綠色環(huán)保當(dāng)量稅額可促進(jìn)傳統(tǒng)企業(yè)的轉(zhuǎn)型升級，亦可推動地區(qū)經(jīng)濟的發(fā)展（王燦等，2005）。此外，綠色環(huán)保當(dāng)量稅額對經(jīng)濟的調(diào)控還存在地域性差別，與各地區(qū)產(chǎn)業(yè)結(jié)構(gòu)形式相關(guān)，對于山西等重工業(yè)地區(qū)，其經(jīng)濟發(fā)展主要依賴于傳統(tǒng)能源行業(yè)，綠色環(huán)保當(dāng)量稅額的大幅度提高，將直接影響該地區(qū)經(jīng)濟的發(fā)展，而對于廣州等發(fā)達(dá)地區(qū)，高新技術(shù)產(chǎn)業(yè)占據(jù)了城市經(jīng)濟發(fā)展的主干，綠色環(huán)保當(dāng)量稅額的提高將會進(jìn)一步優(yōu)化其產(chǎn)業(yè)結(jié)構(gòu)（呂敏等，2018）。不僅如此，綠色環(huán)保當(dāng)量稅額對經(jīng)濟的調(diào)節(jié)作用還時常會受到地方政策、市場競爭等因素的影響（金艷鳴等，2007），而地方政策的調(diào)控和市場的導(dǎo)向性又與各地區(qū)的人文環(huán)境、資源配比等隨機因素相關(guān)。

此外，鄧琨（2019）、王濤等（2018）基于中國東部、中部和西部地區(qū)的稅收結(jié)構(gòu)和稅額，探討了綠色稅收對生態(tài)環(huán)境的影響，研究發(fā)現(xiàn)，不同稅種稅收競爭對生態(tài)環(huán)境的影響效應(yīng)不同，不同地區(qū)稅收競爭對生態(tài)環(huán)境的影響效應(yīng)存在差異。Qin et al.（2012）通過CGE 模型分析了水資源稅收標(biāo)準(zhǔn)對水資源利用效率的影響，發(fā)現(xiàn)水資源稅收標(biāo)準(zhǔn)的提高將有效地降低水資源的用量。Zhang et al.（2017）通過建立CGE 模型剖析了碳稅的變動對能源消耗量的影響，研究發(fā)現(xiàn)，環(huán)保當(dāng)量稅額的適當(dāng)提高將對節(jié)能減排具有十分顯著的效果，但當(dāng)量稅額超過一定區(qū)間后，其調(diào)節(jié)作用將不再顯著。此外，環(huán)保當(dāng)量稅額對環(huán)境的調(diào)節(jié)也將受到市場供需、要素需求、政策決策等因素的共同作用（翁智雄等，2021），而這些因素又與整個社會的環(huán)境相關(guān)，具有鮮明的地域性、隨機性、時變性。

另一方面，綠色環(huán)保當(dāng)量稅額對經(jīng)濟與環(huán)境的雙向調(diào)控。針對如何調(diào)節(jié)環(huán)保當(dāng)量稅額，實現(xiàn)綠色稅收雙重紅利的問題，相關(guān)學(xué)者展開了積極探索。如Tullock（1967）、van Heerden et al.（2006）對綠色稅收的雙重紅利做了較為詳細(xì)的闡述，將雙重紅利劃分為綠色紅利和藍(lán)色紅利，所謂綠色紅利是指在保護生態(tài)環(huán)境上所產(chǎn)生的額外收益，而藍(lán)色紅利是指在推動經(jīng)濟增長方面所產(chǎn)生的額外收益。唐明和明海蓉（2018）運用環(huán)境治理成本分?jǐn)偡?，測算了我國31 個?。▍^(qū)、市）環(huán)境保護稅大氣污染物的最優(yōu)稅率，發(fā)現(xiàn)各地當(dāng)前制定的稅額標(biāo)準(zhǔn)既無法有效彌補環(huán)境治理成本，也不利于推動綠色經(jīng)濟的發(fā)展。王佳鄧等（2021）以江蘇省為例，建立了CGE 模型，分析了環(huán)境保護稅對經(jīng)濟和碳排放量的影響。但以上研究，都僅從綠色環(huán)保當(dāng)量稅額對社會經(jīng)濟和生態(tài)環(huán)境的單向作用或雙重作用進(jìn)行分析，并未揭示三者之間的互動機制。從上述分析不難看出，綠色環(huán)保當(dāng)量稅額對社會經(jīng)濟與生態(tài)環(huán)境的調(diào)控易受多個隨機因素的影響，存在著較為復(fù)雜的非線性、非對稱關(guān)系（呂敏，2015），難以通過一次性轉(zhuǎn)換和傳統(tǒng)的線性回歸模型來描述三者間的互動機制。

Copula 函數(shù)作為直觀揭示各變量之間相依機制的工具，幾乎包含了隨機變量所有的相依信息，彌補了傳統(tǒng)線性模型的缺陷，能夠較好地展示各隨機變量的變化趨勢，反映變量之間的內(nèi)在關(guān)系。在經(jīng)濟宏觀調(diào)控（李述山，2010；吳建華等，2014），糧食安全風(fēng)險評估（李紅芳等，2021），水文水資源評估以及政策調(diào)控，地域差異對碳排放量的影響等方面都有運用（龔誼承，2012；萬埠磊，2014）。這些問題都存在一個共同的特點，變量之間的關(guān)聯(lián)性較為復(fù)雜，隨機性較強，存在非對稱響應(yīng)。而綠色環(huán)保當(dāng)量稅額對環(huán)境與經(jīng)濟的調(diào)控也是一個十分復(fù)雜的過程，與經(jīng)濟宏觀調(diào)控、糧食安全風(fēng)險評估等所面臨的困難相似，即如何揭示各變量間的復(fù)雜非線性關(guān)系是明確各變量調(diào)控機制的難點所在。Copula 函數(shù)可通過概率統(tǒng)計的方法，采用擬合優(yōu)度的檢驗形式，尋求綠色環(huán)保當(dāng)量稅額、社會經(jīng)濟和生態(tài)環(huán)境三個隨機變量的邊緣分布函數(shù)，通過引入邊緣概率分布函數(shù)，以概率的形式體現(xiàn)綠色環(huán)保當(dāng)量稅額、社會經(jīng)濟和生態(tài)環(huán)境三個變量間非線性、非對稱的相關(guān)關(guān)系，尤其容易捕捉到三者的分布尾部的相關(guān)關(guān)系（Patton，2002；Ioan 和Radu，2002；Werner，2003）。此外，由于不限制邊緣分布的選擇，可更好地構(gòu)造三者的多元分布模型，直觀揭示綠色環(huán)保當(dāng)量稅額、社會經(jīng)濟和生態(tài)環(huán)境三者之間的互動過程。為此，本文基于概率統(tǒng)計的思想，運用三維Copula 函數(shù)得到三者的互動規(guī)律，提出了綠色環(huán)保當(dāng)量稅額的最佳波動區(qū)間，這為進(jìn)一步探索實現(xiàn)綠色稅收的雙重紅利提供了新的視角和方法，具有一定的理論意義和實踐價值。

三、機理分析

2018 年我國正式實施《中華人民共和國環(huán)境保護稅法》，逐步完善了我國的綠色稅收體系。而綠色環(huán)保當(dāng)量稅額是以綠色稅收體系為基礎(chǔ)，通過不斷優(yōu)化升級來平衡社會經(jīng)濟發(fā)展與生態(tài)環(huán)境保護的重要手段，理應(yīng)得到重點關(guān)注。為探討綠色環(huán)保當(dāng)量稅額對經(jīng)濟和生態(tài)環(huán)境發(fā)展的作用機理，首先要明確綠色環(huán)保當(dāng)量稅額與經(jīng)濟和生態(tài)環(huán)境的互動機制。

第一，綠色環(huán)保當(dāng)量稅額的設(shè)置不僅可減輕經(jīng)濟發(fā)展中環(huán)境資源的過度開發(fā)、環(huán)境污染的不斷惡化等問題，使企業(yè)、廠商加強對資源的循環(huán)利用（馮俏彬，2018），而且可運用稅收杠桿，進(jìn)一步調(diào)節(jié)社會收入分配，鼓勵特定行業(yè)和區(qū)域發(fā)展，提高綠色經(jīng)濟發(fā)展的均衡性和可持續(xù)性。在綠色發(fā)展理念的引導(dǎo)下，為綠色經(jīng)濟增長注入新鮮動力（武蕊君、張馨，2020）。第二，綠色環(huán)保當(dāng)量稅額的設(shè)置，一方面可鼓勵企業(yè)運用新技術(shù)和新生產(chǎn)模式，助力產(chǎn)業(yè)綠色化發(fā)展。通過實施稅收優(yōu)惠政策，鼓勵企業(yè)主動向綠色發(fā)展轉(zhuǎn)型，為建立健全綠色循環(huán)發(fā)展的經(jīng)濟體系提供政策引導(dǎo)。另一方面通過抑制生產(chǎn)消費負(fù)面行為，推動綠色經(jīng)濟發(fā)展。通過適當(dāng)提高綠色環(huán)保當(dāng)量稅額，實現(xiàn)稅收引導(dǎo)生產(chǎn)消費行為不斷向高效化、綠色化方向轉(zhuǎn)變，也為綠色行業(yè)和環(huán)保行業(yè)的發(fā)展提供了平臺和契機，引領(lǐng)了市場發(fā)展的方向（張彩云，2017）。

由上述分析可知，綠色環(huán)保當(dāng)量稅額的設(shè)置，一方面強有力地融合創(chuàng)新戰(zhàn)略，促使產(chǎn)業(yè)轉(zhuǎn)型升級和綠色發(fā)展，起到了節(jié)能減排的作用；另一方面不斷優(yōu)化生產(chǎn)、調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，也推動了社會經(jīng)濟的高質(zhì)量發(fā)展?？梢姡G色環(huán)保當(dāng)量稅額對社會經(jīng)濟與環(huán)境保護的共同發(fā)展起到了舉足輕重的作用，但其中面臨的一個難題就是如何調(diào)控綠色環(huán)保當(dāng)量稅額以實現(xiàn)社會經(jīng)濟與生態(tài)環(huán)境的協(xié)調(diào)發(fā)展。

四、研究設(shè)計

在綠色發(fā)展理念的指引下，本文基于概率統(tǒng)計的方法，借助Copula 函數(shù)構(gòu)建了不同于現(xiàn)有文獻(xiàn)的綠色環(huán)保當(dāng)量稅額模型，充分考慮綠色環(huán)保當(dāng)量稅額、社會經(jīng)濟、生態(tài)環(huán)境三個變量之間的隨機性、聯(lián)動性以及地域差異性。通過擬合優(yōu)度的檢驗方法，得到了三者的聯(lián)動規(guī)律，就如何確定最佳綠色環(huán)保當(dāng)量稅額問題展開了深入探討。（1）本文的目的在于探索綠色環(huán)保當(dāng)量稅額、社會經(jīng)濟、生態(tài)環(huán)境三者的聯(lián)動效應(yīng)，得到綠色環(huán)保當(dāng)量稅額的合理調(diào)控區(qū)間。（2）通過揭示綠色環(huán)保當(dāng)量稅額、社會經(jīng)濟、生態(tài)環(huán)境三者的聯(lián)動規(guī)律，實現(xiàn)綠色環(huán)保當(dāng)量稅額對經(jīng)濟和環(huán)境發(fā)展的合理控制，深入挖掘復(fù)雜過程背后的潛在規(guī)律，探尋問題背后的理論邏輯，為實現(xiàn)綠色稅收的雙重紅利奠定基礎(chǔ)。（3）本文以我國31 個?。▍^(qū)、市）為研究對象，充分考慮了各地綠色環(huán)保當(dāng)量稅額、經(jīng)濟狀況以及污染排放量的差異，解釋了對象的復(fù)雜性和特殊性，有助于對比和識別不同情景下綠色環(huán)保當(dāng)量稅額調(diào)控的異同，從而全景式地呈現(xiàn)綠色環(huán)保當(dāng)量稅額的調(diào)整過程，以適應(yīng)新時代綠色發(fā)展的現(xiàn)實要求。

（一）指標(biāo)選取

綠色環(huán)保當(dāng)量稅額、社會經(jīng)濟、生態(tài)環(huán)境的協(xié)調(diào)發(fā)展是一個多系統(tǒng)、多變量相互作用、相互影響復(fù)雜迭代的過程，易受多個因素的影響。但選擇多個指標(biāo)評價三者之間的聯(lián)動規(guī)律，將會忽視指標(biāo)的重要性和實用性。如何選擇合理的指標(biāo)來反映三者的發(fā)展?fàn)顩r，對模型的搭建十分重要。

1.綠色環(huán)保當(dāng)量稅額。根據(jù)呂敏（2015）等人的研究，采用綠色環(huán)保當(dāng)量稅額來度量綠色稅收政策對經(jīng)濟與環(huán)境的影響，綠色環(huán)保當(dāng)量稅額的高低，反映了國家征稅力度的大小。綠色環(huán)保當(dāng)量稅額越高，說明國家的征稅力度越大，反之則越小。因此，采用綠色環(huán)保當(dāng)量稅額作為核心變量，能夠較好反映綠色稅收對經(jīng)濟和環(huán)境的調(diào)控能力。

2.社會經(jīng)濟。根據(jù)周思齊等（2021）的研究，采用各地區(qū)的GDP 來反映地區(qū)經(jīng)濟的發(fā)展?fàn)顩r。GDP 是國民經(jīng)濟核算的核心指標(biāo)，也是衡量一個國家或地區(qū)經(jīng)濟狀況和發(fā)展水平的重要指標(biāo)。GDP增長越快，說明該地區(qū)的經(jīng)濟發(fā)展越好，反之則表示該地區(qū)經(jīng)濟發(fā)展還有待進(jìn)一步提高。故采用各地區(qū)GDP 這一核心指標(biāo)，能夠較為直觀地反映地區(qū)經(jīng)濟的變化與發(fā)展。

3.生態(tài)環(huán)境。借鑒劉文玉（2018）、劉潔和李文（2013）的研究，一個地區(qū)生態(tài)環(huán)境的優(yōu)劣程度與該地區(qū)的污染排放量緊密相關(guān)。污染排放量的多少將直接影響該地區(qū)的生態(tài)環(huán)境狀況。限制污染氣體的排放，是我國節(jié)能減排的重要目標(biāo)之一。因此，選取污染氣體排放量來評價一個地區(qū)的生態(tài)環(huán)境較為合理。

（二）綠色環(huán)保當(dāng)量稅額、污染氣體排放量、各地區(qū)生產(chǎn)總值的聯(lián)合概率模型

采用Copula 函數(shù)來建立綠色環(huán)保當(dāng)量稅額、各地區(qū)生產(chǎn)總值（GDP）以及污染氣體排放量之間的聯(lián)合分布模型。常見的Copula 函數(shù)模型大致分為三種：橢圓型、阿基米德族和二次型。其中阿基米德族Copula 函數(shù)適用性較廣、結(jié)構(gòu)形式簡單。阿基米德族Copula 函數(shù)是通過算子φ（又稱生成函數(shù)，一個完全單調(diào)的函數(shù)）構(gòu)造而成（吳建華等，2014）。本文一共有三個隨機變量，故阿基米德族Copula 函數(shù)模型如表1 所示。

表1 Copula 函數(shù)模型的分類

式中φ 是一個連續(xù)、嚴(yán)格遞減凸函數(shù)，被稱為算子式，常見的算子式一共有四種，分別對應(yīng)四種Copula 函數(shù)模型。u1表示綠色環(huán)保當(dāng)量稅額的概率密度邊際分布函數(shù)，u2表示各地區(qū)GDP 的概率密度邊際分布函數(shù)，u3表示各地區(qū)污染氣體的排放量（包括二氧化硫、氮氧化物等）的概率密度邊際分布函數(shù)。

（三）綠色環(huán)保當(dāng)量稅額、污染氣體排放量、各地區(qū)生產(chǎn)總值邊際函數(shù)分布的擬合優(yōu)度檢驗

通過不同的概率分布形式對三個隨機變量的分布形式進(jìn)行擬合，并采用擬合優(yōu)度的方法來獲取三者最優(yōu)的分布形式，確定三個隨機變量的邊緣概率密度函數(shù)，具體可通過概率邊際分布函數(shù)的D值，殘差值（RMSE值）進(jìn)行反映。如式（1）和式（2）所示。

式中ck為觀測樣本的模擬值，mk表示觀測樣本滿足x＜xk的樣本個數(shù)，n 表示樣本總量。殘差值RMSE是指實際樣本值與估計值之間的差值，其中MSE 為均方差，Pei表示經(jīng)驗頻率，Pi為理論計算頻率。邊際分布函數(shù)的精度通過D 值、殘差值等指標(biāo)來反映，D 值和殘差值越小，R2越接近1，則表示所選取的概率邊際分布函數(shù)模型的模擬精度越高。通過以上流程可得到綠色環(huán)保當(dāng)量稅額、污染氣體排放量、各地區(qū)生產(chǎn)總值（GDP）最優(yōu)的邊際分布函數(shù)。

（四）變量的關(guān)聯(lián)性分析及Copula 函數(shù)的合理選擇

運用斯皮爾曼相關(guān)性系數(shù)（Sperman）對三者之間的空間相關(guān)性進(jìn)行檢驗，根據(jù)各類Copula 函數(shù)的適用特征（Gumbel-Hougard Copula 函數(shù)和Clayton Copula 函數(shù)適用于具有正相關(guān)關(guān)系的隨機變量，F(xiàn)rank Copula 函數(shù)對具有正負(fù)相關(guān)關(guān)系的隨機變量都較為適用，Ali-Mikhail-Haq Copula 函數(shù)主要適用于弱相關(guān)關(guān)系的隨機變量），選取合適的Copula 函數(shù)來反映三者之間的聯(lián)動規(guī)律。斯皮爾曼相關(guān)性系數(shù)（Sperman）的求解方式如式（3）所示。

式中：ρs表示自相關(guān)系數(shù)，n 表示三個變量中數(shù)據(jù)的個數(shù)，di表示秩次的差值。

本文的分析流程如下所示：

圖1 模型研究流程

五、實證分析

（一）數(shù)據(jù)來源

本文以我國31 個?。▍^(qū)、市）為基礎(chǔ)研究單元，依據(jù)《中國統(tǒng)計年鑒》《中國能源統(tǒng)計年鑒》《中國環(huán)境年鑒》《中國稅務(wù)年鑒》、國家稅務(wù)總局稅收統(tǒng)計數(shù)據(jù)和相關(guān)專業(yè)數(shù)據(jù)庫。采用綠色環(huán)保當(dāng)量稅額（主要指大氣環(huán)保當(dāng)量稅額）來反映綠色稅收的狀況。通過污染氣體排放量（包括二氧化硫、氮氧化合物等）來評價節(jié)能減排的效果，進(jìn)而量化我國各地區(qū)生態(tài)環(huán)境的變化狀況?；诟鞯貐^(qū)生產(chǎn)總值來反映我國各省市的經(jīng)濟發(fā)展?fàn)顩r?；跀?shù)據(jù)可得性，本文三個變量的具體數(shù)據(jù)如圖2、圖3、圖4 所示。

圖2 各地區(qū)GDP 狀況

圖3 各地區(qū)污染氣體排放狀況

圖4 各地區(qū)綠色環(huán)保當(dāng)量稅額

圖2 展示了我國各地區(qū)GDP 近4 年來的變化狀況，從圖中可看出，近幾年來各地區(qū)的經(jīng)濟都呈現(xiàn)出一定的增長趨勢，尤其是在江浙地區(qū)較為明顯。2018 年以后福建、江蘇、浙江等地的GDP 增長更為明顯。不僅如此，我國各地區(qū)的GDP 存在十分明顯的差異，江浙、廣東等沿海地區(qū)經(jīng)濟發(fā)展較為迅速，而云南、西藏、青海、甘肅、山西等地區(qū)的經(jīng)濟發(fā)展卻較為緩慢，與沿海地區(qū)的經(jīng)濟增長狀況存在顯著差異。但各地區(qū)的污染排放量與各地區(qū)經(jīng)濟的增長狀況不同，由圖3 可知，我國污染氣體排放量最為嚴(yán)重的城市主要集中在山西附近，這與各地區(qū)產(chǎn)業(yè)鏈的差異密切相關(guān)。山西、河北等地區(qū)主要以能源產(chǎn)業(yè)、礦產(chǎn)業(yè)為主，形成了以原材料為基礎(chǔ)的經(jīng)濟圈，造成了資源的不合理分配，既不利于生態(tài)環(huán)境的保護，也不利于地區(qū)經(jīng)濟的發(fā)展。

為此，我國從2018 年開始擬通過征收綠色環(huán)保當(dāng)量稅來限制污染氣體的排放，同時推動社會經(jīng)濟的綠色發(fā)展。從圖3 中可知，自2018 年后山西等重污染地區(qū)污染氣體排放量的增加尺度逐年減小，但每年的排放總量仍居高不下，山西等重污染地區(qū)所設(shè)置的綠色環(huán)保當(dāng)量稅額遠(yuǎn)低于我國平均水平，對污染氣體的排放未能達(dá)到有效的限制性作用。故在設(shè)置綠色環(huán)保當(dāng)量稅額時，不僅要考慮對經(jīng)濟發(fā)展的影響，其生態(tài)環(huán)境效應(yīng)也尤為重要。但總體來講，我國目前對綠色環(huán)保當(dāng)量稅額的設(shè)置還缺乏一套系統(tǒng)的定論，來平衡環(huán)境與經(jīng)濟的發(fā)展?？梢?，應(yīng)將生態(tài)環(huán)境保護、社會經(jīng)濟發(fā)展、綠色環(huán)保當(dāng)量稅額放在同一個維度進(jìn)行考慮，明確三者的互動規(guī)律，才能推動我國綠色經(jīng)濟的發(fā)展。

（二）綠色環(huán)保當(dāng)量稅額、地區(qū)生產(chǎn)總值、污染氣體排放量的概率分布規(guī)律

為探尋三者的互動機制，本文采用概率統(tǒng)計的方法對我國各地區(qū)的綠色環(huán)保當(dāng)量稅額、地區(qū)生產(chǎn)總值、污染氣體排放量的統(tǒng)計數(shù)據(jù)進(jìn)行分析（2018-2020 年），通過擬合優(yōu)度的方法得到各變量最優(yōu)的邊際分布函數(shù)，將擬合優(yōu)度檢驗的誤差D 值、殘差值RMSE以及擬合精度結(jié)果繪于表2。

表2 綠色環(huán)保當(dāng)量稅額、污染氣體排放量及地區(qū)生產(chǎn)總值分布形式的檢驗結(jié)果

續(xù)表

表2 中展示了各地區(qū)綠色環(huán)保當(dāng)量稅額、污染氣體排放量以及地區(qū)生產(chǎn)總值三者的獨立分布形式，從分析結(jié)果可直觀看出，采用對數(shù)正態(tài)分布和正態(tài)分布來模擬各地區(qū)污染氣體排放量、地區(qū)生產(chǎn)總值、綠色環(huán)保當(dāng)量稅額的誤差值和殘差值最小、精度最高。將三種分布形式展示如圖5。

圖5 綠色環(huán)保當(dāng)量稅額、污染氣體排放量以及地區(qū)生產(chǎn)總值的分布規(guī)律

由圖5 可知，我國的綠色環(huán)保當(dāng)量稅額普遍偏低。一方面是由于我國綠色稅收體系處于初步探索階段，與西方發(fā)達(dá)國家相比，我國的工業(yè)化進(jìn)程較晚，面臨的環(huán)境承載壓力日趨接近上限，這才使得我國開始構(gòu)建綠色稅收體系。另一方面是受我國現(xiàn)有產(chǎn)業(yè)結(jié)構(gòu)的影響。相較于西方發(fā)達(dá)國家，我國產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級效果不佳，高新技術(shù)產(chǎn)業(yè)主要集中于東部地區(qū)，中西部地區(qū)仍然普遍存在能源產(chǎn)業(yè)和原材料產(chǎn)業(yè)。雖然我國不斷深化供給側(cè)結(jié)構(gòu)性改革，推動經(jīng)濟發(fā)展在質(zhì)量、效率、動力上的變革，但改革是一項長期的發(fā)展過程，需要較長的時間。

（三）綠色環(huán)保當(dāng)量稅額、地區(qū)生產(chǎn)總值、污染氣體排放量的空間聯(lián)動關(guān)系

運用斯皮爾曼相關(guān)性系數(shù)（Sperman）線性檢驗方法，得到三者的空間線性關(guān)系，根據(jù)三者的空間線性關(guān)系并結(jié)合表1 中四種Copula 函數(shù)模型的特征，得到最為合適的Copula 函數(shù)預(yù)測模型，以便精準(zhǔn)地描述三者的互動規(guī)律，揭示三者的互動機制。通過計算得到污染氣體排放量與綠色環(huán)保當(dāng)量稅額，地區(qū)生產(chǎn)總值與綠色環(huán)保當(dāng)量稅額，污染氣體排放量與地區(qū)生產(chǎn)總值三者之間的空間關(guān)聯(lián)度，如表3 所示。

表3 各地區(qū)綠色環(huán)保當(dāng)量稅額、污染氣體排放量與地區(qū)生產(chǎn)總值的線性關(guān)系

從表3 的計算結(jié)果可知，采用Frank Copula 來模擬三者之間的互動機制最為適宜。此外，由三者間的相關(guān)系數(shù)可知，污染氣體排放量與綠色環(huán)保當(dāng)量稅額兩個變量呈線性負(fù)相關(guān)，說明綠色環(huán)保當(dāng)量稅額的制定對各地區(qū)的節(jié)能減排起到了一定作用，這與圖3 中所展示的結(jié)果較為吻合。由圖3 可知，自2018 年后我國大部分地區(qū)污染氣體排放量的增量出現(xiàn)了降低的趨勢。即便如此，污染氣體排放量與綠色環(huán)保當(dāng)量稅額二者的空間線性關(guān)系并不顯著，截至2020 年，其線性相關(guān)性僅為-0.403，絕對值小于0.7，說明目前我國綠色稅收政策對生態(tài)環(huán)境的推動效應(yīng)還不太顯著。此外，由污染氣體排放量與地區(qū)生產(chǎn)總值的空間線性相關(guān)性可知，二者的相關(guān)性雖出現(xiàn)了逐年下降趨勢，但二者的線性相關(guān)系數(shù)都大于0.7，進(jìn)一步說明我國大部分地區(qū)的經(jīng)濟發(fā)展與能源的消耗量密切相關(guān)。這與呂敏（2015）的研究成果較為吻合，也間接說明了該預(yù)測模型的準(zhǔn)確性。

（四）綠色環(huán)保當(dāng)量稅額、地區(qū)生產(chǎn)總值、污染氣體排放量的互動展示

根據(jù)線性相關(guān)的分析結(jié)果，結(jié)合Copula 函數(shù)的結(jié)構(gòu)特征，采用Frank Copula 函數(shù)預(yù)測模型來模擬綠色環(huán)保當(dāng)量稅額、污染氣體排放量、地區(qū)生產(chǎn)總值三者之間的變化關(guān)系較為合理。其參數(shù)估計值可采用IFM 方法（極大似然法）計算，本文運用matlab 軟件編程計算得到Frank Copula 函數(shù)參數(shù)估計值θ=2.217。通過Frank Copula 概率預(yù)測模型分別計算得到三者的聯(lián)合分布概率，分別通過模型1 和模型2 來展示三者的互動規(guī)律，如圖6 所示。

圖6 中通過模型1、模型2 分別描述了綠色環(huán)保當(dāng)量稅額、污染氣體排放量、地區(qū)生產(chǎn)總值三者之間的關(guān)系。模型1 中展示了綠色環(huán)保當(dāng)量稅額與污染氣體排放量的互動關(guān)系，通過圖上對應(yīng)轉(zhuǎn)折點可較為明顯看出，當(dāng)綠色環(huán)保當(dāng)量稅額低于2.1 元（第一個轉(zhuǎn)折點所對應(yīng)的橫坐標(biāo)）污染氣體的排放量大概率會出現(xiàn)陡增的變化趨勢（第三個轉(zhuǎn)折點前的變化趨勢）。當(dāng)綠色環(huán)保當(dāng)量稅額位于2.1 元至3.3 元（第二個轉(zhuǎn)折點對應(yīng)的橫坐標(biāo)）之間時，此時污染氣體排放量從轉(zhuǎn)折點3 逐漸向轉(zhuǎn)折點4 過渡，污染氣體排放量的變化趨勢平緩，說明綠色環(huán)保當(dāng)量稅額的適當(dāng)提高對各地區(qū)污染氣體的排放量起到了限制性作用。此外，模型2 中展示了綠色環(huán)保當(dāng)量稅額與地區(qū)經(jīng)濟發(fā)展之間的空間關(guān)系，結(jié)合模型1 的結(jié)論及分析，可知綠色環(huán)保當(dāng)量稅額低于2.1 元時，社會經(jīng)濟雖得到了暫時的發(fā)展，但不利于生態(tài)環(huán)境的保護。當(dāng)綠色環(huán)保當(dāng)量稅額高于3.3 元時，污染氣體排放量的增長趨勢雖得到了進(jìn)一步抑制，但社會經(jīng)濟的發(fā)展也受到了一定影響（如模型2 中轉(zhuǎn)折點4 過后的變化趨勢）。故當(dāng)綠色環(huán)保當(dāng)量稅額位于2.1 元至3.3 元之間時，既對污染氣體排放量具有較好的限制性作用，又能促進(jìn)社會經(jīng)濟的發(fā)展。但結(jié)合圖5（a）展示的結(jié)果，不難看出我國絕大部分地區(qū)的綠色環(huán)保當(dāng)量稅額都位于2 元以下，低于最佳環(huán)保當(dāng)量稅額區(qū)間。為了實現(xiàn)經(jīng)濟和環(huán)境的可持續(xù)發(fā)展，各地區(qū)需結(jié)合自身發(fā)展的實際情況，階梯式地提高環(huán)保當(dāng)量稅額。

圖6 綠色環(huán)保當(dāng)量稅額、污染氣體排放量及地區(qū)生產(chǎn)總值的聯(lián)合分布概率模型

由上述分析可知，模型1 和模型2 直觀展示了三者的互動規(guī)律，當(dāng)綠色環(huán)保當(dāng)量稅額低于2.1元時，由于綠色環(huán)保當(dāng)量稅額過低，對環(huán)境的約束效應(yīng)不太顯著，不僅對環(huán)境起不到相應(yīng)的保護作用，也變相地助長了排污企業(yè)的污染行為，使得該地區(qū)的產(chǎn)業(yè)結(jié)構(gòu)單一化、低端化，不利于實現(xiàn)經(jīng)濟的高質(zhì)量發(fā)展。為此，各地政府應(yīng)提高綠色環(huán)保當(dāng)量稅額來推進(jìn)傳統(tǒng)企業(yè)的轉(zhuǎn)型升級，降低企業(yè)的排污效應(yīng)，調(diào)整地區(qū)的產(chǎn)業(yè)結(jié)構(gòu)。但過高的綠色環(huán)保當(dāng)量稅額，將在短時間內(nèi)加重傳統(tǒng)企業(yè)的稅負(fù)，使企業(yè)難以實現(xiàn)轉(zhuǎn)型升級，從而引起一系列社會效應(yīng)。

鑒于此，根據(jù)本文的分析結(jié)果，我國綠色環(huán)保當(dāng)量稅額的較優(yōu)區(qū)間位于2.1 元至3.3 元之間，考慮到我國東西部地區(qū)產(chǎn)業(yè)結(jié)構(gòu)和經(jīng)濟發(fā)展的差異，建議西部地區(qū)的環(huán)保當(dāng)量稅額取2.1 元至2.7元之間，而東部地區(qū)可取2.7 元至3.3 元這一較高的區(qū)段。

為檢驗綠色環(huán)保當(dāng)量稅額、污染氣體排放量以及地區(qū)生產(chǎn)總值的聯(lián)合分布概率模型的準(zhǔn)確性，將其計算結(jié)果與經(jīng)驗公式（4）的計算結(jié)果作對比，對比結(jié)果如圖7 所示。

圖7 預(yù)測模型與經(jīng)驗公式對比結(jié)果

圖7 中的直線表示Frank Copula 概率預(yù)測模型與經(jīng)驗公式（4）的等值線，當(dāng)模型精度為100%時，所有離散點應(yīng)全部落在該等值線上。但在實際情況中，由于系統(tǒng)誤差的影響，這一事件發(fā)生的概率較低。故只要當(dāng)離散點能夠分布在對角線周圍有限的區(qū)間內(nèi)，就認(rèn)為該模型具有較高的精度，通常情況下可通過R2值的大小來判斷該模型精度的高低，R2的取值位于（0-1）之間，當(dāng)R2＞0.7 時，就認(rèn)為該模型具有一定的精度（汪榮鑫，2014）。本文中所得到的預(yù)測模型其R2大于0.9，說明通過Frank Copula 函數(shù)預(yù)測模型來揭示三者的互動規(guī)律具有較高的精度。

（五）穩(wěn)健性檢驗

環(huán)保稅率和當(dāng)量稅額都是政府對環(huán)境規(guī)制的手段，環(huán)保稅率和當(dāng)量稅額越高說明政府的干預(yù)力度越顯著。為了保證上述實證結(jié)果的可靠性，本部分通過更換核心變量的方式做穩(wěn)健性檢驗。綠色環(huán)保稅率的邊際分布函數(shù)檢驗結(jié)果如表4 所示。

表4 綠色環(huán)保稅率、污染氣體排放量及地區(qū)生產(chǎn)總值分布形式的穩(wěn)健性檢驗結(jié)果

續(xù)表

由表4 結(jié)果可知，綠色環(huán)保稅率的邊際函數(shù)與環(huán)保當(dāng)量稅額的分布形式相同，都服從對數(shù)分布函數(shù)。運用斯皮爾曼相關(guān)性系數(shù)（Sperman）線性檢驗方法，得到綠色環(huán)保稅率與地區(qū)生產(chǎn)總值、污染氣體排放量三者之間的空間關(guān)聯(lián)度，如表5 所示。

表5 綠色環(huán)保稅率與地區(qū)生產(chǎn)總值、污染氣體排放量的空間關(guān)聯(lián)度

表5 展示了環(huán)保稅率與地區(qū)生產(chǎn)總值、污染氣體排放量的空間關(guān)聯(lián)度。從分析結(jié)果可看出，環(huán)保稅率與地區(qū)生產(chǎn)總值呈正相關(guān)關(guān)系，說明環(huán)保稅率的設(shè)置對地區(qū)生產(chǎn)總值的提高有正向激勵作用，而對污染氣體的排放量呈負(fù)相關(guān)關(guān)系。相關(guān)系數(shù)逐年升高，說明環(huán)保稅率對社會經(jīng)濟和生態(tài)環(huán)境的調(diào)控能力逐漸提高，但各年的相關(guān)系數(shù)仍低于0.7，說明三者之間的關(guān)聯(lián)度還需進(jìn)一步加強。根據(jù)三者的關(guān)聯(lián)度分析結(jié)果，結(jié)合各種Copula 函數(shù)的選擇條件，發(fā)現(xiàn)采用Frank Copula 函數(shù)來模擬三者之間的互動機制最為適宜。

通過更換核心變量，采用綠色環(huán)保稅率代替綠色環(huán)保當(dāng)量稅額，結(jié)果表明：綠色環(huán)保稅率與社會經(jīng)濟和生態(tài)環(huán)境的相關(guān)性關(guān)系與綠色環(huán)保當(dāng)量稅額和二者的關(guān)聯(lián)度結(jié)果基本一致，即對地區(qū)生產(chǎn)總值的提高有正向激勵作用，而對污染氣體的排放量有相應(yīng)的抑制效應(yīng)，但二者對社會經(jīng)濟與生態(tài)環(huán)境的調(diào)控作用都較為有限。

六、結(jié)論與建議

本文對我國近年來綠色環(huán)保當(dāng)量稅額、污染氣體排放量以及地區(qū)生產(chǎn)總值做了相應(yīng)統(tǒng)計，分析了三者間的關(guān)聯(lián)度，并運用三維Copula 函數(shù)建立了三者的聯(lián)動模型，以期為我國綠色稅收的改革提供相應(yīng)的參考。

通過對我國近三年來的綠色環(huán)保當(dāng)量稅額、污染氣體排放量和地區(qū)生產(chǎn)總值的現(xiàn)狀進(jìn)行分析，發(fā)現(xiàn)廣東、江蘇等少數(shù)發(fā)達(dá)地區(qū)既保持了經(jīng)濟的快速增長，又較為有效地控制了污染氣體的排放量。綠色環(huán)保當(dāng)量稅額對經(jīng)濟和環(huán)境的協(xié)同作用并不顯著，如山西、內(nèi)蒙古等地區(qū)污染排放量居高不下，而其經(jīng)濟發(fā)展也較為遲緩。相比廣東、江蘇等發(fā)達(dá)地區(qū)，山西等地區(qū)的產(chǎn)業(yè)結(jié)構(gòu)較為傳統(tǒng)和單一，既不利于經(jīng)濟的高質(zhì)量發(fā)展，也不利于環(huán)境的保護。

通過空間關(guān)聯(lián)度的分析，分別得到污染氣體排放量、綠色環(huán)保當(dāng)量稅額和地區(qū)生產(chǎn)總值三者的空間線性關(guān)系。結(jié)果表明：污染氣體排放量與地區(qū)生產(chǎn)總值的空間線性關(guān)系較高，說明我國大部分地區(qū)的經(jīng)濟增長仍與能源消耗密切相關(guān)，主要還是依賴于傳統(tǒng)資源行業(yè)的發(fā)展。而污染氣體排放量與綠色環(huán)保當(dāng)量稅額，綠色環(huán)保當(dāng)量稅額和地區(qū)生產(chǎn)總值二者的線性相關(guān)關(guān)系較弱，都低于0.7，說明我國現(xiàn)有的綠色稅收體系還有待進(jìn)一步完善。

通過三維Copula 函數(shù)模型，揭示綠色環(huán)保當(dāng)量稅額、污染氣體排放量、地區(qū)生產(chǎn)總值三者的互動關(guān)系，得到了綠色環(huán)保當(dāng)量稅額的最佳取值區(qū)間為2.1 元至3.3 元?？紤]到我國東西部地區(qū)產(chǎn)業(yè)結(jié)構(gòu)和經(jīng)濟發(fā)展的差異，可根據(jù)不同區(qū)域制定差異化的綠色環(huán)保當(dāng)量稅額標(biāo)準(zhǔn)。如山西、吉林等地區(qū)主要以第二產(chǎn)業(yè)為主，若綠色環(huán)保當(dāng)量稅額過高，不利于傳統(tǒng)企業(yè)的轉(zhuǎn)型升級，從而影響地區(qū)經(jīng)濟的發(fā)展。而上海、廣東等較為發(fā)達(dá)地區(qū)以高新產(chǎn)業(yè)為主，因此適當(dāng)提高綠色環(huán)保當(dāng)量稅額，有利于推動綠色環(huán)保技術(shù)的發(fā)展，可適當(dāng)增加綠色經(jīng)濟發(fā)展的比重。

基于上述的研究結(jié)論，本文得出以下幾點政策啟示：

第一，根據(jù)不同區(qū)域制定差異化的綠色環(huán)保當(dāng)量稅額。由于各地區(qū)經(jīng)濟發(fā)展程度和主體功能的差異，因此需實現(xiàn)對不同地區(qū)綠色環(huán)保當(dāng)量稅額的動態(tài)調(diào)整與把控。北上廣等地區(qū)的經(jīng)濟水平和治理能力較高，民眾環(huán)保意識較強，可適度提升環(huán)保稅額以促進(jìn)綠色低碳環(huán)保技術(shù)的發(fā)展，增加綠色經(jīng)濟發(fā)展的比重，確定其主導(dǎo)地位。作為產(chǎn)業(yè)轉(zhuǎn)移承接主體的中西部地區(qū)，若征收過高的環(huán)保稅額，不利于傳統(tǒng)企業(yè)的轉(zhuǎn)型升級，影響地區(qū)經(jīng)濟的穩(wěn)步發(fā)展。

第二，構(gòu)建以綠色發(fā)展為導(dǎo)向的政績考核評價體系和機制，彰顯各地政府發(fā)展綠色經(jīng)濟的決心，推動“唯GDP”的短期政績考核向“綠色GDP”的中長期政績考核轉(zhuǎn)變。同時，在新常態(tài)背景下推動經(jīng)濟的高質(zhì)量發(fā)展，需要將保護環(huán)境和發(fā)展綠色經(jīng)濟作為政府政績考核的要件，有力落實生態(tài)環(huán)境污染的終身責(zé)任制。

第三，拓展環(huán)境保護稅的征收范圍。當(dāng)前我國明確了“碳達(dá)峰”“碳中和”的時間表和路線圖，全國碳市場上線交易已經(jīng)正式啟動，二氧化碳排放也成為評價環(huán)境優(yōu)劣的重要參考指標(biāo)。為此，應(yīng)著力將碳稅作為實現(xiàn)這一目標(biāo)的重要經(jīng)濟手段，將其納入環(huán)境保護稅當(dāng)中，推動我國生態(tài)文明建設(shè)和實現(xiàn)美麗中國的目標(biāo)。

第四，制定相應(yīng)的運用返還政策。與綠色稅收“雙管齊下”，通過補貼的手段既降低不同主體的經(jīng)濟損失，也轉(zhuǎn)變了企業(yè)的發(fā)展方式、提升了企業(yè)的自主創(chuàng)新能力，從而形成與征收綠色稅收相配套的長效機制，推動環(huán)境與經(jīng)濟的可持續(xù)發(fā)展。