孫菲 李友俊 鞏艷芬 彭民

摘要：基于工業(yè)綠色發(fā)展的內(nèi)涵，分別構(gòu)建了工業(yè)綠色發(fā)展水平和效率兩子系統(tǒng)的評(píng)價(jià)指標(biāo)體系，運(yùn)用熵權(quán)-TOPSIS模型和SBM模型測(cè)度了油氣城市的工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率，運(yùn)用耦合協(xié)調(diào)度模型研究了兩子系統(tǒng)的耦合協(xié)調(diào)效應(yīng)。結(jié)果表明：各地區(qū)工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)顯著，隨時(shí)間變化呈梯度遞增格局;2004—2018年各油氣城市兩子系統(tǒng)耦合協(xié)調(diào)度總體呈上升趨勢(shì)，耦合協(xié)調(diào)效應(yīng)較強(qiáng);各油氣城市耦合協(xié)調(diào)效應(yīng)隨時(shí)間推移差距有縮小趨勢(shì)，耦合協(xié)調(diào)度類型分級(jí)明顯，不同油氣城市的耦合協(xié)調(diào)效應(yīng)存在差距?；谏鲜鼋Y(jié)論，發(fā)現(xiàn)優(yōu)質(zhì)人才、先進(jìn)技術(shù)、區(qū)域差異性、環(huán)境規(guī)制等是影響兩系統(tǒng)耦合協(xié)調(diào)效應(yīng)的重要因素。最后，從完善人才機(jī)制、提高綠色技術(shù)創(chuàng)新投入、構(gòu)建協(xié)同發(fā)展機(jī)制等方面，提出了油氣城市工業(yè)綠色發(fā)展的提升策略。

關(guān)鍵詞：油氣城市;工業(yè)綠色發(fā)展水平;工業(yè)綠色發(fā)展效率;耦合協(xié)調(diào)度

中圖分類號(hào)：F124

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1673-5595（2022）02-0027-13

油氣城市，即以石油和天然氣開采加工為支柱產(chǎn)業(yè)的資源型城市，在其產(chǎn)業(yè)結(jié)構(gòu)中資源采掘及與之相關(guān)的產(chǎn)業(yè)占有較大比重。新常態(tài)下，隨著我國(guó)資源環(huán)境趨緊，在資源環(huán)境承載力的約束下，油氣城市的工業(yè)綠色發(fā)展已成為必然。《“十四五”工業(yè)綠色發(fā)展規(guī)劃》中明確提出：到2025年，工業(yè)產(chǎn)業(yè)結(jié)構(gòu)、生產(chǎn)方式綠色低碳轉(zhuǎn)型取得顯著成效，綠色低碳技術(shù)裝備廣泛應(yīng)用，能源資源利用效率大幅提高，綠色制造水平全面提升?？梢姡苿?dòng)油氣城市的工業(yè)綠色發(fā)展意味著將綠色發(fā)展理念貫穿于工業(yè)生產(chǎn)全過(guò)程、全領(lǐng)域，降低工業(yè)生產(chǎn)對(duì)油氣資源的過(guò)度消耗及對(duì)環(huán)境的破壞，積極推進(jìn)工業(yè)能效提升，減少污染排放，實(shí)現(xiàn)城市工業(yè)增長(zhǎng)與環(huán)境保護(hù)的“雙贏”。區(qū)別于已有文獻(xiàn)對(duì)工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率進(jìn)行單獨(dú)研究，

本文以油氣城市為研究對(duì)象，將工業(yè)綠色發(fā)展水平和效率相結(jié)合，研究二者的耦合協(xié)調(diào)效應(yīng)，并分析形成原因，由此提出提升油氣城市工業(yè)綠色發(fā)展的相關(guān)建議。

一、文獻(xiàn)綜述

（一）工業(yè)綠色發(fā)展水平的相關(guān)研究

已有文獻(xiàn)中關(guān)于工業(yè)綠色發(fā)展水平的研究主要涉及評(píng)價(jià)指標(biāo)體系和評(píng)價(jià)方法?，F(xiàn)有的評(píng)價(jià)指標(biāo)體系不僅考慮了經(jīng)濟(jì)、資源、環(huán)境等指標(biāo)，也考慮了科技、社會(huì)和政府支持等指標(biāo)。有的研究借鑒了國(guó)家發(fā)展改革委制定的《綠色發(fā)展指標(biāo)體系》（2016）中的部分指標(biāo)，構(gòu)建工業(yè)綠色發(fā)展水平評(píng)價(jià)指標(biāo)體系，該指標(biāo)體系涵蓋了工業(yè)綠色增長(zhǎng)、工業(yè)資源消耗與環(huán)境影響和政府工業(yè)綠色政策支持等方面。[1-2]有的研究從工業(yè)增長(zhǎng)水平、工業(yè)研發(fā)強(qiáng)度、政府支持力度、工業(yè)能耗水平、污染治理效率5個(gè)方面構(gòu)建工業(yè)綠色發(fā)展水平評(píng)價(jià)指標(biāo)體系。[3]還有研究結(jié)合現(xiàn)代工業(yè)發(fā)展要求，從綠色生產(chǎn)、綠色排放、綠色保護(hù)、綠色創(chuàng)新、綠色質(zhì)效5個(gè)方面選取各類相關(guān)指標(biāo)，構(gòu)建工業(yè)綠色發(fā)展水平評(píng)價(jià)指標(biāo)體系。[4-5]關(guān)于工業(yè)綠色發(fā)展水平評(píng)價(jià)方法，一般側(cè)重采用TOPSIS法[6]、熵值法[4]、動(dòng)態(tài)面板的系統(tǒng)GMM估計(jì)方法[7]、因子分析法[3]和k-means聚類算法[5]等方法。

（二）工業(yè)綠色發(fā)展效率的相關(guān)研究

關(guān)于工業(yè)綠色發(fā)展效率的研究最早起源于綠色技術(shù)效率，Pittman[8]最早利用指數(shù)測(cè)度了包含非期望產(chǎn)出的綠色技術(shù)效率。此后，德國(guó)學(xué)者Schaltegger等[9]提出了工業(yè)生態(tài)效率的概念，這一概念現(xiàn)已被廣泛應(yīng)用于工業(yè)領(lǐng)域。針對(duì)生產(chǎn)中的非期望產(chǎn)出與投入產(chǎn)出變量的松弛性問(wèn)題，Tone[10]提出了SBM模型，該模型是非徑向和非角度的，能有效解決非期望產(chǎn)出與投入產(chǎn)出變量的松弛性問(wèn)題，在測(cè)定環(huán)境和資源效率方面得到廣泛應(yīng)用。目前國(guó)內(nèi)外對(duì)工業(yè)綠色發(fā)展效率的研究一般是從投入產(chǎn)出的角度，選擇適當(dāng)?shù)哪Ｐ瓦M(jìn)行測(cè)度研究，重點(diǎn)分析現(xiàn)狀并對(duì)未來(lái)進(jìn)行預(yù)測(cè)。也有研究在進(jìn)行效率評(píng)估時(shí)，采用了在規(guī)模報(bào)酬可變條件下的DEA模型以及超效率SBM模型。以王建民等[11]、張仁杰等[12]、王喜平等[13]、黃磊等[14]、丁顯有等[15]為代表的學(xué)者，采用全局超效率SBM模型，選取資本、勞動(dòng)力、技術(shù)、資源作為投入指標(biāo)，選取地區(qū)生產(chǎn)總值作為期望產(chǎn)出指標(biāo)，選取工業(yè)廢水、工業(yè)廢氣、工業(yè)固體廢棄物排放量等作為非期望產(chǎn)出，測(cè)度了我國(guó)30個(gè)?。ㄗ灾螀^(qū)、直轄市）和長(zhǎng)江經(jīng)濟(jì)帶的工業(yè)綠色發(fā)展效率，研究發(fā)現(xiàn)中國(guó)省域工業(yè)綠色發(fā)展效率在時(shí)間維度上總體呈現(xiàn)出波動(dòng)上升的趨勢(shì)，且工業(yè)綠色發(fā)展效率的地區(qū)差異十分明顯，以省為單位測(cè)度的工業(yè)綠色發(fā)展效率處于波動(dòng)增長(zhǎng)的狀態(tài)。

綜上所述，現(xiàn)有研究在研究維度上，側(cè)重工業(yè)綠色發(fā)展水平評(píng)估、工業(yè)綠色發(fā)展效率測(cè)度;在研究對(duì)象上，側(cè)重以全國(guó)、省域、城市群等為對(duì)象，進(jìn)行工業(yè)綠色發(fā)展程度評(píng)估;在研究方法上，側(cè)重采用定量分析工具測(cè)度工業(yè)綠色發(fā)展效率，分析其影響因素。目前已有研究多集中在以全國(guó)、區(qū)域或城市為研究對(duì)象的工業(yè)綠色發(fā)展水平或工業(yè)綠色發(fā)展效率的單方面研究上，將工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率相結(jié)合的研究相對(duì)較少，特別是關(guān)于油氣城市工業(yè)綠色發(fā)展的全面研究更為罕見。鑒于此，本文將構(gòu)建工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率的綜合評(píng)價(jià)指標(biāo)體系，對(duì)油氣城市的工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率進(jìn)行綜合評(píng)價(jià)，并在此基礎(chǔ)上分析油氣城市工業(yè)綠色發(fā)展水平和效率的耦合協(xié)調(diào)效應(yīng)及其形成原因，提出油氣城市工業(yè)綠色發(fā)展的提升對(duì)策。

二、工業(yè)綠色發(fā)展水平與效率耦合協(xié)調(diào)的內(nèi)涵

工業(yè)綠色發(fā)展是指工業(yè)綠色低碳轉(zhuǎn)型發(fā)展。中國(guó)社會(huì)科學(xué)院工業(yè)經(jīng)濟(jì)研究所課題組將“工業(yè)綠色發(fā)展”定義為：在工業(yè)領(lǐng)域有效利用資源、能源，降低污染物排放，以獲得經(jīng)濟(jì)效益和環(huán)境效益為目標(biāo)，通過(guò)綠色技術(shù)創(chuàng)新，最終實(shí)現(xiàn)工業(yè)生產(chǎn)各環(huán)節(jié)的綠色化。[16-18]因此工業(yè)綠色發(fā)展的衡量標(biāo)準(zhǔn)為提高資源利用效率、環(huán)境效率、全要素生產(chǎn)率，即提高工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率。

（一）工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率

工業(yè)綠色發(fā)展水平反映的是工業(yè)綠色發(fā)展過(guò)程的成效，是對(duì)現(xiàn)有工業(yè)綠色發(fā)展程度的綜合體現(xiàn)。綜合已有研究，發(fā)現(xiàn)工業(yè)綠色發(fā)展水平主要通過(guò)工業(yè)經(jīng)濟(jì)發(fā)展、資源消耗強(qiáng)度、環(huán)境影響及治理、技術(shù)創(chuàng)新4個(gè)方面來(lái)體現(xiàn)。工業(yè)經(jīng)濟(jì)發(fā)展是保障當(dāng)前工業(yè)綠色發(fā)展的動(dòng)力，它體現(xiàn)了工業(yè)綠色發(fā)展的先導(dǎo)條件，可以分為經(jīng)濟(jì)效益與經(jīng)濟(jì)結(jié)構(gòu)兩個(gè)維度。資源消耗強(qiáng)度指工業(yè)生產(chǎn)中企業(yè)生產(chǎn)單位工業(yè)增加值所占用和耗費(fèi)的資源。[4]環(huán)境影響及治理體現(xiàn)了污染排放對(duì)工業(yè)綠色發(fā)展的影響以及政府的效力對(duì)工業(yè)綠色發(fā)展的支撐，分為工業(yè)污染排放和工業(yè)污染治理兩個(gè)維度。技術(shù)創(chuàng)新是提升油氣城市工業(yè)綠色發(fā)展水平的核心驅(qū)動(dòng)力[19-20]，衡量技術(shù)創(chuàng)新，主要看技術(shù)創(chuàng)新投入高低和技術(shù)創(chuàng)新體系是否完善[14]。

工業(yè)綠色發(fā)展效率是決定工業(yè)綠色發(fā)展的內(nèi)生動(dòng)力。從宏觀層面看，工業(yè)綠色發(fā)展效率包括全要素生產(chǎn)率、技術(shù)效率和規(guī)模效率等[21];從微觀層面看，工業(yè)綠色發(fā)展效率是指工業(yè)企業(yè)在一定生產(chǎn)條件下，“投入”與“產(chǎn)出”兩者之間的比率關(guān)系。

（二）工業(yè)綠色發(fā)展水平與效率的耦合協(xié)調(diào)效應(yīng)

工業(yè)綠色發(fā)展水平與效率的耦合協(xié)調(diào)效應(yīng)是指工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率的互補(bǔ)發(fā)展程度，反映工業(yè)綠色發(fā)展水平對(duì)工業(yè)綠色發(fā)展效率的支撐作用及工業(yè)綠色發(fā)展效率對(duì)工業(yè)綠色發(fā)展水平的引致效應(yīng)。一般采用耦合協(xié)調(diào)度模型測(cè)度工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)。油氣城市工業(yè)綠色發(fā)展進(jìn)程中，工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率兩大子系統(tǒng)互為前提，且相輔相成，二者與工業(yè)綠色發(fā)展的關(guān)系如圖1所示。

1.耦合協(xié)調(diào)

耦合概念來(lái)源于物理學(xué)，反映的是兩個(gè)或兩個(gè)以上的系統(tǒng)相互作用、彼此影響進(jìn)而達(dá)到融合的現(xiàn)象。[22]協(xié)調(diào)指的是系統(tǒng)與要素在發(fā)展過(guò)程中彼此和諧一致的程度，體現(xiàn)系統(tǒng)由無(wú)序走向有序的趨勢(shì)。[23]耦合度與協(xié)調(diào)度是對(duì)系統(tǒng)或要素之間耦合與協(xié)調(diào)狀態(tài)、程度的描述和度量，體現(xiàn)了整個(gè)系統(tǒng)如何走向融合發(fā)展和良性循環(huán)的過(guò)程。[24]

協(xié)調(diào)就是正確處理組織內(nèi)外各種關(guān)系，為組織正常運(yùn)轉(zhuǎn)創(chuàng)造良好的條件和環(huán)境，促進(jìn)組織目標(biāo)的實(shí)現(xiàn)。

2.耦合協(xié)調(diào)度

耦合度反映系統(tǒng)無(wú)序與有序狀態(tài)之間的轉(zhuǎn)換過(guò)程，體現(xiàn)系統(tǒng)內(nèi)部參量間的相互作用，度量系統(tǒng)之間在發(fā)展過(guò)程中的和諧程度;而耦合協(xié)調(diào)度則是從整體上看工業(yè)綠色發(fā)展水平子系統(tǒng)與工業(yè)綠色發(fā)展效率子系統(tǒng)耦合關(guān)系，其結(jié)果比耦合度的分析結(jié)果更全面更合理。[25]工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率是工業(yè)綠色發(fā)展系統(tǒng)的兩個(gè)方面，二者之間是互補(bǔ)發(fā)展關(guān)系，這兩個(gè)子系統(tǒng)之間應(yīng)該共促共進(jìn)、協(xié)調(diào)一致、共同發(fā)展。通過(guò)測(cè)度二者之間的耦合協(xié)調(diào)度，能夠判斷出二者是否和諧發(fā)展。[26-27]定量分析工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率的耦合協(xié)調(diào)程度，可以更細(xì)致、更清楚地發(fā)現(xiàn)工業(yè)綠色發(fā)展過(guò)程中所產(chǎn)生的問(wèn)題，從而找到解決問(wèn)題的辦法，促進(jìn)油氣城市經(jīng)濟(jì)與環(huán)境的可持續(xù)發(fā)展。

油氣城市要實(shí)現(xiàn)工業(yè)綠色發(fā)展，需要生態(tài)環(huán)境趨于良性循環(huán)。工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率耦合協(xié)調(diào)度達(dá)到良好水平，才能實(shí)現(xiàn)人與自然的和諧共生。

三、工業(yè)綠色發(fā)展水平與效率耦合協(xié)調(diào)效應(yīng)模型構(gòu)建

（一）評(píng)價(jià)指標(biāo)體系構(gòu)建

1.油氣城市工業(yè)綠色發(fā)展水平評(píng)價(jià)指標(biāo)體系

本文參考和借鑒了李琳等[1]、雷勛平等[28]、張永凱等[2]、刁莉娟等[29]、張仁杰等[12]、李曉鐘等[30]建立的工業(yè)綠色發(fā)展水平評(píng)價(jià)指標(biāo)體系，結(jié)合當(dāng)前的政策導(dǎo)向、指標(biāo)的代表性以及各指標(biāo)數(shù)據(jù)的可得性，遵循合理性、規(guī)范性和科學(xué)性原則，經(jīng)過(guò)仔細(xì)遴選，確定了油氣城市工業(yè)綠色發(fā)展水平評(píng)價(jià)指標(biāo)體系的構(gòu)成，具體見表1。

2.油氣城市工業(yè)綠色發(fā)展效率評(píng)價(jià)指標(biāo)體系

由于城市工業(yè)綠色發(fā)展不僅要實(shí)現(xiàn)城市經(jīng)濟(jì)發(fā)展，還需要進(jìn)行城市環(huán)境污染治理，改善和恢復(fù)城市生態(tài)。因此，油氣城市工業(yè)綠色發(fā)展效率應(yīng)綜合反映經(jīng)濟(jì)增長(zhǎng)、社會(huì)發(fā)展、環(huán)境保護(hù)和資源節(jié)約四者的協(xié)調(diào)發(fā)展關(guān)系。效率的測(cè)度一般選擇DEA方法，評(píng)價(jià)指標(biāo)通常劃分為投入和產(chǎn)出兩類。考慮投入產(chǎn)出指標(biāo)數(shù)與決策單元數(shù)（本文決策單元為14個(gè)地級(jí)油氣城市）的關(guān)系，最終選擇以資本存量、勞動(dòng)從業(yè)總?cè)藬?shù)、能源消耗總量

為投入指標(biāo);以經(jīng)濟(jì)效益為期望產(chǎn)出指標(biāo)、環(huán)境污染為非期望產(chǎn)出指標(biāo)，分別用城市工業(yè)增加值、環(huán)境污染綜合指數(shù)（根據(jù)工業(yè)廢水、工業(yè)二氧化硫和工業(yè)煙塵的排放量計(jì)算）來(lái)表示[31-32]，具體的投入產(chǎn)出指標(biāo)見表2。

（二）綜合評(píng)價(jià)模型

1.工業(yè)綠色發(fā)展水平測(cè)度模型——熵權(quán)-TOPSIS模型

工業(yè)綠色發(fā)展水平評(píng)價(jià)方法屬于多指標(biāo)評(píng)價(jià)綜合法。通常會(huì)采用熵值法、AHP法、神經(jīng)網(wǎng)絡(luò)及TOPSIS等方法來(lái)測(cè)算區(qū)域發(fā)展水平。但這些方法在應(yīng)用過(guò)程中存在一些問(wèn)題，如熵值法的權(quán)數(shù)依賴于樣本本身，沒有考慮各指標(biāo)間的橫向比較;AHP法主觀性較強(qiáng)，定性成分多;神經(jīng)網(wǎng)絡(luò)運(yùn)算所需參數(shù)數(shù)量多，需要海量數(shù)據(jù)進(jìn)行訓(xùn)練從而建立模型;TOPSIS 法有時(shí)會(huì)使評(píng)價(jià)對(duì)象與正、負(fù)理想解集距離都近，不能真實(shí)反映評(píng)價(jià)對(duì)象的優(yōu)劣程度。[29]工業(yè)綠色發(fā)展水平是關(guān)于資源消耗、經(jīng)濟(jì)效益與結(jié)構(gòu)、技術(shù)投入與創(chuàng)新等的綜合系統(tǒng)，對(duì)其評(píng)價(jià)需要遵循系統(tǒng)分析的思想。基于此，本文擬將TOPSIS模型與熵值法結(jié)合，采用熵值法確定權(quán)重，進(jìn)而評(píng)價(jià)油氣城市的工業(yè)綠色發(fā)展水平。

TOPSIS模型是指在目標(biāo)空間中定義一個(gè)測(cè)度，計(jì)算某一評(píng)價(jià)對(duì)象與正理想解和負(fù)理想解之間的加權(quán)距離，以測(cè)量目標(biāo)靠近正理想解和遠(yuǎn)離負(fù)理想解的程度，以此評(píng)估工業(yè)綠色發(fā)展水平。TOPSIS是一種距離綜合評(píng)價(jià)法。[33]同時(shí)，作為客觀賦權(quán)法的熵權(quán)法，在確定權(quán)數(shù)時(shí)依賴指標(biāo)的客觀數(shù)值，規(guī)避了主觀判斷，結(jié)果具有客觀性。因此，本文選擇熵權(quán)法計(jì)算出各個(gè)指標(biāo)的權(quán)重[34]，基本步驟如下。

設(shè)m為城市數(shù)，n為評(píng)價(jià)指標(biāo)數(shù)，工業(yè)綠色發(fā)展水平的原始評(píng)價(jià)矩陣為：

X=X11…X1n

Xm1…Xmnm×n（1）

（1）數(shù)據(jù)標(biāo)準(zhǔn)化。

各個(gè)指標(biāo)具有不同的量綱，為了能夠?qū)λ兄笜?biāo)進(jìn)行統(tǒng)一的運(yùn)算處理，選用極值法對(duì)相關(guān)數(shù)據(jù)進(jìn)行處理。

正向指標(biāo)（越大越優(yōu)型）為：

Vij=Xij-min（Xj）max（Xj）-min（Xj） （i=1， 2，…， m;j=1，2，…， n）?? （2）

負(fù)向指標(biāo)（越小越優(yōu)型）為：

Vij=max（Xj）-Xijmax（Xj）-min（Xj） （i=1， 2，…， m;j=1， 2，…， n）?? （3）

（2）計(jì)算各項(xiàng)指標(biāo)的比重。

第j項(xiàng)指標(biāo)下第i個(gè)樣本對(duì)象占該指標(biāo)的比重為：

Pij=Vij∑mi=1Vij（i=1， 2，…， m;j=1， 2，…， n）?? （4）

（3）計(jì)算第j項(xiàng)指標(biāo)的熵值。

信息熵值為ej=-K×∑mi=1Pijln（Pij），常數(shù)K可取值為K=1/lnm，（m為樣本數(shù)），0≤ej≤1。

信息效用值為：

gj=1-ej，（i=1， 2，…， m;j=1， 2，…， n）?? （5）

（4）計(jì)算熵權(quán)。

wj=gj∑nj=1gj，（j=1，2，…，n）（6）

（5）計(jì)算基于熵權(quán)的加權(quán)決策矩陣。

利用熵權(quán)法確定的指標(biāo)權(quán)重向量wj來(lái)計(jì)算加權(quán)規(guī)范化決策矩陣，即R=（Rij）m×n通過(guò)矩陣V的每一行與其相應(yīng)的權(quán)重wj相乘得到，即加權(quán)決策矩陣R=（Rij）m×n=wj×V。

R=R11…R1n

Rm1…Rmn=

V11w1…V1nwn

Vm1wm…Vmnwn? （7）

（6）確定正、負(fù)理想解。

正理想解為各指標(biāo)的最優(yōu)指標(biāo)值，即最偏好的方案;負(fù)理想解為各指標(biāo)的最劣指標(biāo)值，即最不偏好的方案。令R+表示正理想解，R-表示負(fù)理想解，則有：

R+={maxRi1，maxRi2，maxRi3…maxRin}?? （8）

R-={minRi1，minRi2，minRi3…minRin}? （i=1， 2， …， m）（9）

（7）計(jì)算數(shù)據(jù)與正負(fù)理想解之間的歐氏距離。

與正理想解的歐氏距離為：

D+i=∑mi=1（Rij-R+j）2（i=1， 2， …， m）?? （10）

與負(fù)理想解的歐氏距離為：

D-i=∑mi=1（Rij-R-j）2（i=1，2…，m）?? （11）

式中：Rij為第i個(gè)樣本第j項(xiàng)指標(biāo)的加權(quán)規(guī)范化值;Rj+、Rj-分別為第j項(xiàng)指標(biāo)在評(píng)價(jià)中的最偏好值和最不偏好值。

（8）計(jì)算相對(duì)貼近度。

ηi=D-iD+i+D-i? （12）

式中：ηi為指標(biāo)對(duì)最優(yōu)解的貼近度，是介于0和1之間的數(shù)值。貼近度的值越大說(shuō)明越接近最優(yōu)解，即該地區(qū)的工業(yè)綠色發(fā)展水平越高;相反，貼近度的值越小說(shuō)明越接近最劣解，該地區(qū)的工業(yè)綠色發(fā)展水平越低。當(dāng)ηi=1時(shí)，工業(yè)綠色發(fā)展水平最高，工業(yè)綠色發(fā)展達(dá)到最優(yōu)狀態(tài);當(dāng)ηi=0時(shí)，工業(yè)綠色發(fā)展處于高度無(wú)序混亂狀態(tài)。[35]

2.工業(yè)綠色發(fā)展效率測(cè)度模型——SBM模型

當(dāng)前學(xué)術(shù)界常用的效率評(píng)價(jià)方法主要有非參數(shù)分析法和參數(shù)分析法，其中，非參數(shù)分析法以數(shù)據(jù)包絡(luò)分析法（DEA）為主，參數(shù)分析法則以隨機(jī)前沿分析法（SFA）為主。由于SFA假定隨機(jī)誤差服從特定概率分布，并依托于多重約束下的生產(chǎn)函數(shù)，難以滿足現(xiàn)實(shí)數(shù)據(jù)及條件要求[36]，因此學(xué)術(shù)界較少選用SFA法進(jìn)行效率評(píng)價(jià)。而DEA法借助于投入產(chǎn)出數(shù)據(jù)，從線性規(guī)劃的思想出發(fā)，對(duì)各個(gè)決策單元偏離生產(chǎn)前沿面的程度進(jìn)行比較，測(cè)度評(píng)價(jià)其相對(duì)效率。DEA不受生產(chǎn)函數(shù)的理論約束，更具嚴(yán)謹(jǐn)性和科學(xué)性，因此被廣大學(xué)者所接納和改進(jìn)。[37]

DEA中CCR和BCC等傳統(tǒng)模型均是從徑向和角度兩方面測(cè)度效率，忽視了投入產(chǎn)出變量的松弛性，使得計(jì)算結(jié)果虛高、效率值不夠精確。為解決這些問(wèn)題，Tone于2001年提出了基于非徑向、非角度的SBM 模型（slacks-based measure），這種方法的優(yōu)勢(shì)在于其投入和產(chǎn)出是一個(gè)無(wú)量綱，對(duì)效率值不會(huì)產(chǎn)生影響，而且效率值會(huì)隨著投入和產(chǎn)出松弛程度的變化而嚴(yán)格單調(diào)遞減，能夠有效地處理傳統(tǒng)DEA模型存在的缺陷。[10]本文采用SBM模型，測(cè)度油氣城市的工業(yè)綠色發(fā)展效率，模型如下：

假設(shè)在生產(chǎn)過(guò)程中有n個(gè)決策單元（DMU），每個(gè)DMU均有一個(gè)投入向量X，兩個(gè)產(chǎn)出向量，即期望產(chǎn)出Yg與非期望產(chǎn)出Yb，三個(gè)向量可分別定義為：

X=（x1，x2，…，xn）∈Rm×n，xi∈Rm;

Yg=（yg1，yg2，…，ygn）∈Rs1×n，ygi ∈Rs1;

Yb=（yb1，yb2，…，ybn）∈Rs2×n，ygi ∈Rs2? （13）

式中：xi≥0，yig≥0，yib≥0。在規(guī)模報(bào)酬不變前提下，可定義生產(chǎn)可能集P為：

P={（xi，ygi，ybi）xi≥Xλ，ygi≤Ygλ，ybi≥Ybλ，λ≥0}（14）

根據(jù)Tone給出的處理辦法，DMU的效率值ρ′可表示為：

ρ′=min1-1m∑mi=1si-xi01+1s1+s2∑s1r=1sgrygr0+∑s2r=1sbrybr0

s.t.x0=Xλ+s-

yg0=Ygλ-sg

yb0=Ybλ+sb

λ≥0，s-≥0，sg≥0，sb≥0（15）

式中：λ為權(quán)重向量;s-為投入的松弛變量，表示投入過(guò)剩;sg、sb為產(chǎn)出的松弛變量，表示產(chǎn)出不足和非期望產(chǎn)出過(guò)多。目標(biāo)函數(shù)ρ′關(guān)于s-、sg、sb嚴(yán)格遞減，且0≤ρ′≤1。對(duì)于特定的DMU，當(dāng)ρ′=0時(shí)，說(shuō)明該DMU是無(wú)效的;當(dāng)ρ′=1，s-=sg=sb=0時(shí)，該DMU是有效的。

3.工業(yè)綠色發(fā)展水平和效率耦合協(xié)調(diào)效應(yīng)模型——耦合協(xié)調(diào)度模型

（1）耦合度模型

工業(yè)綠色發(fā)展水平與效率的耦合度模型可表示為：

C=L×E（L+E）/2? （16）

式中：L為工業(yè)綠色發(fā)展水平;E為工業(yè)綠色發(fā)展效率;C為二者的耦合度，C越大，表明工業(yè)綠色發(fā)展水平與效率的耦合度越大，兩子系統(tǒng)的相關(guān)性越強(qiáng)。

（2）耦合協(xié)調(diào)度模型

在某些情況下，耦合度雖能反映工業(yè)綠色發(fā)展水平和效率的作用強(qiáng)度和方向，但無(wú)法準(zhǔn)確反映二者的整體功效和耦合協(xié)調(diào)效應(yīng)[14]，即可能出現(xiàn)“假性耦合”。為了避免類似假性耦合現(xiàn)象的發(fā)生，克服耦合度模型的缺陷，本文在耦合度模型的基礎(chǔ)上引入耦合協(xié)調(diào)度，構(gòu)建工業(yè)綠色發(fā)展水平和效率的耦合協(xié)調(diào)度模型。

T=αL+βE? （17）

D=C×T? （18）

式中：α和β為各子系統(tǒng)的權(quán)重，且滿足α+β=1。由于本文研究工業(yè)綠色發(fā)展水平與效率這兩個(gè)子系統(tǒng)之間的耦合協(xié)調(diào)關(guān)系，故認(rèn)為各子系統(tǒng)對(duì)油氣城市的工業(yè)綠色發(fā)展同等重要，將權(quán)重α和β均取值為0.5。T為二者的綜合協(xié)調(diào)指數(shù)。D為工業(yè)綠色發(fā)展水平和效率的耦合協(xié)調(diào)度，反映二者的耦合協(xié)調(diào)效應(yīng)。耦合協(xié)調(diào)度D介于[0，1]之間，D的值越接近1，說(shuō)明兩子系統(tǒng)相互作用、相互影響的效果越顯著，兩子系統(tǒng)之間的耦合協(xié)調(diào)關(guān)系越強(qiáng);當(dāng)D=1時(shí)，兩子系統(tǒng)達(dá)到良性互動(dòng)的最優(yōu)協(xié)調(diào)發(fā)展階段;D=0，表明工業(yè)綠色發(fā)展水平和效率兩子系統(tǒng)之間沒有必然聯(lián)系，對(duì)彼此的現(xiàn)狀和未來(lái)發(fā)展沒有影響。[38]

（3）耦合度和耦合協(xié)調(diào)度的等級(jí)劃分

耦合度和耦合協(xié)調(diào)度能夠表征兩子系統(tǒng)的耦合與協(xié)調(diào)關(guān)系，依據(jù)中值分段原理，對(duì)耦合度的等級(jí)進(jìn)行劃分，[38-39]標(biāo)準(zhǔn)如表3所示。

表3 工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率耦合度類型劃分

耦合度0≤C≤0.40.4≤C<0.50.5≤C<0.60.6≤C<0.70.7≤C<0.80.8≤C≤1

發(fā)展階段低水平耦合頡頏階段初步磨合階段中級(jí)磨合階段深度磨合階段高水平耦合階段

本文關(guān)于耦合協(xié)調(diào)度與耦合協(xié)調(diào)類型的劃分主要參考吳傳清等[6]的分類方法，將耦合協(xié)調(diào)發(fā)展?fàn)顩r分為三大類十大亞類（見表4）。

（三）數(shù)據(jù)來(lái)源

參考已有文獻(xiàn)所界定的油氣城市，結(jié)合其特征，列出目前我國(guó)29個(gè)油氣城市，主要包括大慶、濮陽(yáng)、盤錦、東營(yíng)、玉門等，其中東部地區(qū)有8個(gè)、中部地區(qū)有5個(gè)、西部地區(qū)有16個(gè)，其中新疆最多，有7個(gè)，占全國(guó)的四分之一。

為了更好地分析油氣城市工業(yè)綠色發(fā)展水平與效率的耦合協(xié)調(diào)效應(yīng)，本文選取的時(shí)間為近15年，即選擇2004—2018年為研究時(shí)段。考慮各城市數(shù)據(jù)的可比性，本文以地級(jí)油氣城市作為實(shí)證評(píng)價(jià)對(duì)象，29個(gè)油氣城市中地級(jí)市有18個(gè)，如表5所示。由于哈密、慶陽(yáng)、昭通、撫順的研究指標(biāo)數(shù)據(jù)缺失嚴(yán)重，考慮到數(shù)據(jù)的可得性，本文搜集了東部、中部和西部地區(qū)14個(gè)具有代表性的油氣城市2004—2018年的指標(biāo)數(shù)據(jù)進(jìn)行研究。數(shù)據(jù)來(lái)源于EPS數(shù)據(jù)庫(kù)、《中國(guó)城市統(tǒng)計(jì)年鑒》、各省市統(tǒng)計(jì)年鑒及官方網(wǎng)站。在收集整理數(shù)據(jù)的過(guò)程中，發(fā)現(xiàn)存在一定數(shù)量的缺失數(shù)據(jù)，本文使用SPASS軟件中的回歸分析填補(bǔ)缺失數(shù)據(jù)。

四、工業(yè)綠色發(fā)展水平與效率耦合協(xié)調(diào)效應(yīng)實(shí)證分析

（一）各地區(qū)及全國(guó)油氣城市工業(yè)綠色發(fā)展水平與效率的耦合協(xié)調(diào)效應(yīng)分析

根據(jù)上述模型測(cè)算油氣城市的工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率，結(jié)果見表6和表7。在已測(cè)算的工業(yè)綠色發(fā)展水平貼進(jìn)度和工業(yè)綠色發(fā)展效率值的基礎(chǔ)上，測(cè)算各地區(qū)及全國(guó)油氣城市工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率的耦合度和耦合協(xié)調(diào)度，結(jié)果見表8和表9。通過(guò)表8可以看出各地區(qū)兩子系統(tǒng)15年間的耦合度均大于0.7，說(shuō)明各地區(qū)的工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率之間具有比較密切的關(guān)系。東部和西部地區(qū)的工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率15年間始終處于高水平耦合階段（但東部地區(qū)耦合度大部分時(shí)間高于西部地區(qū)），中部地區(qū)由開始的深度磨合階段逐漸上升到高水平耦合階段。可見，工業(yè)綠色發(fā)展水平與效率的耦合度中部地區(qū)相對(duì)較弱、東部地區(qū)最高、西部地區(qū)居中。從全國(guó)來(lái)看，15年間兩子系統(tǒng)關(guān)系密切，一直處于高水平耦合階段。

2004—2018年全國(guó)油氣城市工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率兩子系統(tǒng)耦合協(xié)調(diào)效應(yīng)保持平穩(wěn)增強(qiáng)態(tài)勢(shì)（見表9），由2004年的0.712持續(xù)增長(zhǎng)至2016年的0.789，后稍回落至2018年的0.771，年均增長(zhǎng)0.57%，耦合協(xié)調(diào)發(fā)展類型始終處于中級(jí)協(xié)調(diào)階段。各地區(qū)工業(yè)綠色發(fā)展水平與效率互動(dòng)效應(yīng)良好，二者互為支撐，協(xié)同推進(jìn)。各地區(qū)油氣城市工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)顯著，隨時(shí)間變化呈梯度遞增格局。2004—2010年，東部地區(qū)工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)最強(qiáng)，西部地區(qū)耦合協(xié)調(diào)效應(yīng)逐漸增強(qiáng)，中部地區(qū)耦合協(xié)調(diào)效應(yīng)最弱;2010—2014年，西部地區(qū)耦合協(xié)調(diào)效應(yīng)明顯高于東部和中部地區(qū)，東部和中部地區(qū)耦合協(xié)調(diào)效應(yīng)基本平穩(wěn)且較接近;2014—2018年，中部地區(qū)耦合協(xié)調(diào)效應(yīng)繼續(xù)增強(qiáng)，高于東部和西部地區(qū)，東部地區(qū)耦合協(xié)調(diào)效應(yīng)自2016年開始小幅回落，西部則仍保持平穩(wěn)增長(zhǎng)。

東部地區(qū)工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)呈“M型”增中有降態(tài)勢(shì)，由2004年的0.747上升至2009年的0.772，后經(jīng)多次反復(fù)升降，降至2018年的0.725，整體年均下降0.22%，始終處于中級(jí)協(xié)調(diào)階段。中部地區(qū)兩子系統(tǒng)的耦合協(xié)調(diào)效應(yīng)由2004年的0.634持續(xù)提升至2018年的0.785，年均增長(zhǎng)1.54%，呈明顯的線性增長(zhǎng)態(tài)勢(shì)，由初級(jí)協(xié)調(diào)階段上升至中級(jí)協(xié)調(diào)階段。其中2016和2017年，達(dá)到了良好協(xié)調(diào)。西部地區(qū)工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)呈“V型”緩慢增長(zhǎng)態(tài)勢(shì)，由2004年的0.747增長(zhǎng)至2018年的0.8，年均增長(zhǎng)0.85%，由中級(jí)協(xié)調(diào)發(fā)展階段穩(wěn)步上升到良好協(xié)調(diào)階段。

分析呈現(xiàn)上述空間分布特征的原因，從工業(yè)綠色發(fā)展水平方面看，2004—2018年間，西部地區(qū)油氣城市的工業(yè)綠色發(fā)展水平顯著提高，東部地區(qū)明顯下降，中部地區(qū)雖略有提高，但整體仍處于全國(guó)中間位置。作為東部地區(qū)的代表城市——大慶和東營(yíng)，其工業(yè)綠色發(fā)展水平在2004—2018年間明顯下降;西部地區(qū)的代表城市新疆克拉瑪依，其工業(yè)綠色發(fā)展水平則波動(dòng)上升。2015年，國(guó)際原油價(jià)格持續(xù)低位運(yùn)行，我國(guó)原油產(chǎn)量大幅度降低，中石油、中石化和中海油紛紛減少在油氣勘探開發(fā)上的投入。作為長(zhǎng)期以來(lái)我國(guó)產(chǎn)量較高的兩大油田——黑龍江大慶油田和山東勝利油田，自2015年開始，原油產(chǎn)量也開始出現(xiàn)明顯下滑趨勢(shì)。大慶作為黑龍江乃至東北三省的重要油氣城市，為保證油氣資源的可持續(xù)開采，原油已由原來(lái)的高產(chǎn)、逐漸減產(chǎn)轉(zhuǎn)向現(xiàn)在的穩(wěn)產(chǎn)。數(shù)據(jù)顯示，自2015年大慶油田原油持續(xù)減產(chǎn)，至2018年原油產(chǎn)量減至3 200萬(wàn)噸，年均減產(chǎn)幅度達(dá)130萬(wàn)噸。大慶油田減產(chǎn)使地方經(jīng)濟(jì)持續(xù)承壓，經(jīng)測(cè)算，僅2015年黑龍江省GDP因大慶油田原油減產(chǎn)而減少約200億元，財(cái)政減收約60億元。2015年以后，大慶市的單位工業(yè)用地增加值和單位工業(yè)電耗增加值不斷下降，單位工業(yè)增加值能耗不斷提高，導(dǎo)致城市工業(yè)綠色發(fā)展水平持續(xù)下降。東營(yíng)市是勝利油田的誕生地，也是山東省重要的油氣城市，受原油減產(chǎn)影響，人才流失、產(chǎn)業(yè)結(jié)構(gòu)單一、人口老齡化嚴(yán)重等資源型城市共性問(wèn)題再次凸顯，也間接影響了城市工業(yè)綠色發(fā)展水平的提高。克拉瑪依作為中國(guó)西部重要的油氣城市，其原油產(chǎn)量波動(dòng)較小。近些年，克拉瑪依尤其注重產(chǎn)業(yè)綠色發(fā)展，加快工業(yè)領(lǐng)域的科技創(chuàng)新，以此提升經(jīng)濟(jì)發(fā)展質(zhì)量。如新疆油田公司吉慶油田作業(yè)區(qū)從2016年起，通過(guò)科技攻關(guān)，建成了多井叢、大平臺(tái)、井站一體式無(wú)桿泵采油示范平臺(tái)，實(shí)現(xiàn)了油田的智能化、集約式管理。因此，克拉瑪依的工業(yè)綠色發(fā)展水平15年間呈波動(dòng)上升趨勢(shì)。

從工業(yè)綠色發(fā)展效率方面看，東部、中部和西部的油氣城市工業(yè)綠色發(fā)展效率存在顯著差異，西部地區(qū)效率值相對(duì)較高，是其他地區(qū)趕超的目標(biāo)，而東部和中部油氣城市工業(yè)綠色發(fā)展效率則持續(xù)處于相對(duì)較低水平。東部地區(qū)工業(yè)經(jīng)濟(jì)發(fā)達(dá)，城市平均工業(yè)增加值高于中部和西部地區(qū)，長(zhǎng)期以來(lái)既是我國(guó)工業(yè)生產(chǎn)主力，也是能源消耗主力，同時(shí)GDP增長(zhǎng)伴隨著大量污染物的排放，導(dǎo)致環(huán)境污染相對(duì)嚴(yán)重，其工業(yè)綠色發(fā)展效率偏低。西部地區(qū)效率值偏高，說(shuō)明西部地區(qū)達(dá)到有效發(fā)展的油氣城市較其他兩地區(qū)所占比重大，西部6個(gè)油氣城市中只有克拉瑪依的效率均值小于1，主要是因?yàn)槲鞑康貐^(qū)油氣城市平均能耗低于中部和東部地區(qū)，且西部地區(qū)油氣城市環(huán)境污染綜合指數(shù)也明顯低于其他兩個(gè)地區(qū)，因此西部地區(qū)油氣城市的工業(yè)綠色發(fā)展效率明顯高于中部和東部地區(qū)。中部地區(qū)效率值在2014年后基本高于其他兩地，原因是中部地區(qū)環(huán)境污染在2004—2007年間較嚴(yán)重，2007—2014年間有所改善，2015年后環(huán)境污染治理已見成效。

綜上，工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率耦合協(xié)調(diào)效應(yīng)表現(xiàn)為：東部地區(qū)工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)始終處于中級(jí)協(xié)調(diào)階段，中部和西部目前均到達(dá)了良好協(xié)調(diào)階段。因此，應(yīng)在大力支持東部地區(qū)發(fā)展的基礎(chǔ)上，加快帶動(dòng)全國(guó)油氣城市的工業(yè)綠色發(fā)展。

（二）各油氣城市工業(yè)綠色發(fā)展水平與效率耦合協(xié)調(diào)效應(yīng)分析

大慶、盤錦、東營(yíng)、濮陽(yáng)、克拉瑪依，在2004—2018年間始終處于高水平耦合階段;唐山、鄂爾多斯、南陽(yáng)、瀘州、延安、榆林，初期處于深度磨合階段，后來(lái)向高度耦合階段轉(zhuǎn)變;松原、南充、達(dá)州，則是前期深度磨合階段和高度耦合階段交替出現(xiàn)，后期穩(wěn)定在高度耦合階段。各油氣城市15年間的耦合度均大于0.7，這說(shuō)明各油氣城市的工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率之間關(guān)系較密切。具體數(shù)據(jù)見表10和表11。

由表11可知，首先，從耦合協(xié)調(diào)效應(yīng)絕對(duì)水平看，油氣城市工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率兩子系統(tǒng)耦合協(xié)調(diào)效應(yīng)較強(qiáng)，大多數(shù)油氣城市呈波動(dòng)上升趨勢(shì)。勉強(qiáng)協(xié)調(diào)共計(jì)出現(xiàn)5次，分別是2004年和2005年濮陽(yáng)、2004年松原、2005年盤錦、2007年達(dá)州;優(yōu)質(zhì)協(xié)調(diào)出現(xiàn)4次，分別是2011年、2012年和2013年克拉瑪依與2009年?yáng)|營(yíng)。其他油氣城市在不同年份里表現(xiàn)為初級(jí)、中級(jí)和良好，其中中級(jí)較多，說(shuō)明各油氣城市在15年間兩子系統(tǒng)的耦合協(xié)調(diào)效應(yīng)均較強(qiáng)。但有的城市工業(yè)綠色發(fā)展水平呈波動(dòng)下降趨勢(shì)，也有城市工業(yè)綠色發(fā)展處于無(wú)效率狀態(tài)，這樣的城市兩子系統(tǒng)的耦合協(xié)調(diào)度相對(duì)較低，如大慶和東營(yíng)。因此，從總體趨勢(shì)看，雖然部分油氣城市的耦合協(xié)調(diào)度在不同年份均出現(xiàn)下降，但除了大慶和東營(yíng)最終減少外，其他油氣城市均是增加的，工業(yè)綠色發(fā)展協(xié)同態(tài)勢(shì)呈現(xiàn)良好上升趨勢(shì)，工業(yè)綠色發(fā)展向前推進(jìn)。

其次，從耦合協(xié)調(diào)效應(yīng)發(fā)展趨勢(shì)看，各油氣城市在某一年差距不大，隨時(shí)間推移差距有縮小趨勢(shì)，且各油氣城市的耦合協(xié)調(diào)度類型分級(jí)較明顯。差距最大的是2007年?yáng)|營(yíng)（0.882 7）和達(dá)州（0.541 8），相差0.340 9;最小的是2014年大慶（0.856 1）和盤錦（0.681 1），相差0.175。從變化上看，變異系數(shù)由2004年的0.13下降到2018年的0.08，各油氣城市耦合協(xié)調(diào)效應(yīng)差距減小。另外，從分級(jí)看，除2004—2007年間盤錦、松原、濮陽(yáng)和達(dá)州出現(xiàn)勉強(qiáng)協(xié)調(diào)，2009年?yáng)|營(yíng)、2011—2013年克拉瑪依出現(xiàn)優(yōu)質(zhì)協(xié)調(diào)外，其他年份均出現(xiàn)初級(jí)、中級(jí)、良好三類分級(jí)，分級(jí)類型較多且明顯，說(shuō)明不同油氣城市的耦合協(xié)調(diào)效應(yīng)存在差距。其中，唐山、盤錦、松原、南陽(yáng)、濮陽(yáng)、瀘州、達(dá)州均是前五年（2004—2008年）出現(xiàn)初級(jí)協(xié)調(diào)或者勉強(qiáng)協(xié)調(diào);延安前五年為中級(jí)協(xié)調(diào)，后十年耦合協(xié)調(diào)效應(yīng)明顯增強(qiáng)，出現(xiàn)中級(jí)或者良好協(xié)調(diào)。需要注意的是，大慶在2004—2016年間一直處于良好協(xié)調(diào)狀態(tài)，2017年降為中級(jí)協(xié)調(diào)，2018年降為初級(jí)協(xié)調(diào);克拉瑪依15年間多數(shù)處于中級(jí)和良好協(xié)調(diào)狀態(tài)，在2011—2013年間出現(xiàn)優(yōu)質(zhì)協(xié)調(diào)。大慶作為我國(guó)第一大油田所在地，近年受原油減產(chǎn)、國(guó)際油價(jià)持續(xù)低迷等影響，城市工業(yè)經(jīng)濟(jì)發(fā)展放緩，工業(yè)綠色發(fā)展水平呈下降趨勢(shì)，但工業(yè)綠色發(fā)展效率尚保持穩(wěn)定，居全國(guó)中游水平，因此兩子系統(tǒng)協(xié)調(diào)程度降低。克拉瑪依作為西部地區(qū)的油氣城市，其工業(yè)綠色發(fā)展水平較高，但工業(yè)綠色發(fā)展效率偏低。近年，克拉瑪依在經(jīng)濟(jì)社會(huì)發(fā)展的全過(guò)程中積極貫徹綠色發(fā)展理念，如積極引導(dǎo)企業(yè)研發(fā)或引進(jìn)低消耗、低排污的先進(jìn)技術(shù)，推進(jìn)廢棄物循環(huán)利用，提升了城市工業(yè)綠色發(fā)展效率，使兩子系統(tǒng)一度達(dá)到優(yōu)質(zhì)協(xié)調(diào)。

最后，綜合以上分析可以看出，兩子系統(tǒng)的耦合協(xié)調(diào)度在2004—2018年間總體呈上升趨勢(shì)，即2004—2010年小幅上升，2010—2016年相對(duì)平穩(wěn)，2016—2018年出現(xiàn)小幅回落。兩子系統(tǒng)的協(xié)調(diào)性和均衡性不斷增強(qiáng)，各油氣城市兩子系統(tǒng)的耦合協(xié)調(diào)效應(yīng)差異隨時(shí)間推移在變小，反映了二者的協(xié)調(diào)度越來(lái)越好，表明兩子系統(tǒng)處于動(dòng)態(tài)耦合協(xié)調(diào)發(fā)展?fàn)顟B(tài)。雖然兩子系統(tǒng)之間的動(dòng)態(tài)協(xié)調(diào)度不斷提升，但需要注意的是2010年前后兩子系統(tǒng)發(fā)展趨勢(shì)出現(xiàn)差異。2004—2010年間工業(yè)綠色發(fā)展水平和效率均呈上升趨勢(shì)，兩子系統(tǒng)耦合協(xié)調(diào)度呈上升趨勢(shì);2010—2016年間工業(yè)綠色發(fā)展效率出現(xiàn)下降，工業(yè)綠色發(fā)展水平仍保持小幅上升;2016—2018年間工業(yè)綠色發(fā)展效率上升，工業(yè)綠色發(fā)展水平下降，也就是說(shuō)，2010年之后兩子系統(tǒng)耦合協(xié)調(diào)度穩(wěn)中有降?？梢?，未來(lái)油氣城市的工業(yè)綠色發(fā)展不僅應(yīng)關(guān)注其綠色發(fā)展水平，也應(yīng)關(guān)注其綠色發(fā)展效率，降低對(duì)環(huán)境資源的破壞，關(guān)注環(huán)境質(zhì)量建設(shè)，真正實(shí)現(xiàn)工業(yè)綠色協(xié)調(diào)發(fā)展，達(dá)到工業(yè)綠色發(fā)展水平和效率的和諧統(tǒng)一。

五、結(jié)論與建議

（一）結(jié)論

（1）從各地區(qū)耦合協(xié)調(diào)度看，各地區(qū)工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)顯著，隨時(shí)間變化呈梯度遞增格局。中部和西部地區(qū)兩子系統(tǒng)均達(dá)到了良好協(xié)調(diào)階段，東部地區(qū)兩子系統(tǒng)的協(xié)調(diào)發(fā)展始終處于中級(jí)協(xié)調(diào)階段。因此，應(yīng)加強(qiáng)對(duì)東部地區(qū)工業(yè)綠色發(fā)展的支持，全面帶動(dòng)加快全國(guó)油氣城市的工業(yè)綠色發(fā)展。

（2）從各油氣城市耦合協(xié)調(diào)度看，2004—2018年各油氣城市工業(yè)綠色發(fā)展水平與工業(yè)綠色發(fā)展效率兩子系統(tǒng)耦合協(xié)調(diào)度總體呈上升趨勢(shì)，耦合協(xié)調(diào)效應(yīng)較強(qiáng)。表明油氣城市兩子系統(tǒng)的協(xié)調(diào)度越來(lái)越好，兩子系統(tǒng)處于動(dòng)態(tài)耦合協(xié)調(diào)發(fā)展?fàn)顟B(tài)，工業(yè)綠色發(fā)展向前推進(jìn)。

（3）從各油氣城市兩子系統(tǒng)耦合協(xié)調(diào)效應(yīng)的發(fā)展趨勢(shì)看，其在某一年的差距不大，且隨時(shí)間推移差距有縮小趨勢(shì)，但不同油氣城市兩子系統(tǒng)的耦合協(xié)調(diào)效應(yīng)存在差距，且耦合協(xié)調(diào)度類型分級(jí)較明顯。

（二）建議

基于上述分析，發(fā)現(xiàn)優(yōu)質(zhì)人才、先進(jìn)技術(shù)、區(qū)域差異性、環(huán)境規(guī)制等是影響兩系統(tǒng)耦合協(xié)調(diào)效應(yīng)的重要因素。因此，本文提出以下對(duì)策建議。

（1）完善人才培育機(jī)制。以油氣城市人才需求為出發(fā)點(diǎn)，充分發(fā)揮企業(yè)、高校、科研院所的作用，依托重大專項(xiàng)、重要工程，構(gòu)建“產(chǎn)學(xué)研”聯(lián)合體，在有條件的油氣城市建立高新技術(shù)企業(yè)創(chuàng)新創(chuàng)業(yè)平臺(tái)，吸引高科技人才到油氣城市創(chuàng)業(yè)發(fā)展。

（2）繼續(xù)提高綠色技術(shù)創(chuàng)新投入。首先，通過(guò)各種稅費(fèi)優(yōu)惠手段激勵(lì)企業(yè)研發(fā)關(guān)鍵技術(shù)，尤其是能夠提高油氣資源利用效率的技術(shù)。其次，引導(dǎo)創(chuàng)新要素向新興產(chǎn)業(yè)、服務(wù)業(yè)等環(huán)境友好型產(chǎn)業(yè)流動(dòng)，使過(guò)度集中于油氣產(chǎn)業(yè)的研發(fā)資金和高技術(shù)人才等創(chuàng)新要素流動(dòng)起來(lái)。同時(shí)，通過(guò)加大財(cái)政投入落實(shí)獎(jiǎng)勵(lì)機(jī)制，鼓勵(lì)企業(yè)的創(chuàng)新積極性，加快項(xiàng)目成果應(yīng)用，密切聯(lián)系高等院校，關(guān)注科技項(xiàng)目合作。如東營(yíng)于2019年成立了“油地校協(xié)同創(chuàng)新聯(lián)盟”，借此推動(dòng)油、地、校三方深度融合，以科技及科技成果轉(zhuǎn)化驅(qū)動(dòng)城市綠色轉(zhuǎn)型，為油氣城市工業(yè)綠色發(fā)展探索出一條新的“出路”。

（3）構(gòu)建區(qū)域協(xié)同發(fā)展新機(jī)制。西部地區(qū)近年工業(yè)綠色發(fā)展效率明顯提升，但其經(jīng)濟(jì)發(fā)展程度明顯落后于東部地區(qū)，導(dǎo)致該地區(qū)培育出的優(yōu)質(zhì)要素由于“虹吸效應(yīng)”流向東部發(fā)達(dá)地區(qū)。為此，應(yīng)構(gòu)建區(qū)域協(xié)同發(fā)展新機(jī)制，利用區(qū)位優(yōu)勢(shì)，建設(shè)開放合作實(shí)驗(yàn)區(qū)，東部、中部和西部油氣城市合作引進(jìn)新材料、現(xiàn)代物流、文化娛樂(lè)等產(chǎn)業(yè)項(xiàng)目，發(fā)展服務(wù)業(yè)和旅游業(yè)，打破區(qū)域壁壘，提高地區(qū)自身吸納人才、科技、產(chǎn)品等要素的能力，完善制度體系，實(shí)現(xiàn)區(qū)域優(yōu)質(zhì)要素從東部向中、西部流動(dòng)，最終實(shí)現(xiàn)工業(yè)綠色發(fā)展要素在東部、中部和西部地區(qū)的均衡分布。

（4）推行示范區(qū)試點(diǎn)政策。在油氣城市推行示范區(qū)試點(diǎn)政策，選擇示范區(qū)試點(diǎn)城市，加大資金支持力度，提升城市創(chuàng)新能力，以此助推城市產(chǎn)業(yè)轉(zhuǎn)型升級(jí)，使其產(chǎn)業(yè)結(jié)構(gòu)更趨于合理。首先，以老油氣城市為核心，穩(wěn)妥推進(jìn)產(chǎn)業(yè)轉(zhuǎn)型升級(jí)示范區(qū)試點(diǎn)范圍擴(kuò)展，提升示范區(qū)政策的實(shí)施效果。其次，加強(qiáng)示范區(qū)試點(diǎn)城市間的分工協(xié)作，對(duì)薄弱環(huán)節(jié)、薄弱地區(qū)加強(qiáng)政策設(shè)計(jì)，實(shí)施示范區(qū)獎(jiǎng)勵(lì)機(jī)制和退出機(jī)制。最后，加強(qiáng)示范區(qū)建設(shè)的政策支持力度，加快傳統(tǒng)產(chǎn)業(yè)改造，突出創(chuàng)新理念，培育壯大新興產(chǎn)業(yè)，提高創(chuàng)新質(zhì)量。

參考文獻(xiàn)：

[1] 李琳，張佳.長(zhǎng)江經(jīng)濟(jì)帶工業(yè)綠色發(fā)展水平差異及其分解——基于2004—2013年108個(gè)城市的比較研究[J].軟科學(xué)，2016，30（11）：48-53.

[2] 張永凱，崔佳新.山東省城市工業(yè)綠色發(fā)展水平評(píng)價(jià)[J].蘭州財(cái)經(jīng)大學(xué)學(xué)報(bào)，2019，35（1）：33-41.

[3] 黃郭平，劉子杰.我國(guó)工業(yè)綠色發(fā)展水平之綜合評(píng)價(jià)[J].湖南工業(yè)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2020，25（4）：61-68.

[4] 張彩霞，馬春旺.區(qū)域工業(yè)綠色發(fā)展水平評(píng)價(jià)指標(biāo)體系研究[J].統(tǒng)計(jì)與管理，2021，36（1）：108-113.

[5] 馬越，孫曉雅.省域工業(yè)綠色發(fā)展水平評(píng)價(jià)研究[J].綠色科技，2021，23（6）：277-280.

[6] 吳傳清，黃磊.長(zhǎng)江經(jīng)濟(jì)帶工業(yè)綠色發(fā)展績(jī)效評(píng)估及其協(xié)同效應(yīng)研究[J].中國(guó)地質(zhì)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2018，18（3）：46-55.

[7] 徐成龍，莊貴陽(yáng).供給側(cè)改革驅(qū)動(dòng)中國(guó)工業(yè)綠色發(fā)展的動(dòng)力結(jié)構(gòu)及時(shí)空效應(yīng)[J].地理科學(xué)，2018，38（6）：849-858.

[8] Pittman R W. Multilateral Productivity Comparisons with Undesirable Outputs[J]. The Economic Journal， 1983（93）：883-891.

[9] Schaltegger S， Burritt R. Contemporary Environmental Accounting. Issues， Concepts and Practice[M]. Sheffield：Greenleaf publishing， 2000.

[10] Tone K. A Slacks-based Measure of Efficiency in Data Envelopment Analysis[J]. European Journal of Operational Research， 2001，130（3）：498-509.

[11] 王建民，仇定三，蔣倩穎，等.長(zhǎng)江經(jīng)濟(jì)帶工業(yè)綠色發(fā)展效率測(cè)量與提升路徑研究[J].科技管理研究，2019，39（12）：46-52.

[12] 張仁杰，董會(huì)忠.基于省級(jí)尺度的中國(guó)工業(yè)生態(tài)效率的時(shí)空演變及影響因素[J].經(jīng)濟(jì)地理，2020，40（7）：124-132.

[13] 王喜平，李艷梅.可持續(xù)發(fā)展視角下工業(yè)綠色發(fā)展測(cè)度及其差異演化[J].工業(yè)技術(shù)經(jīng)濟(jì)，2020，39（7）：83-91.

[14] 黃磊，吳傳清.長(zhǎng)江經(jīng)濟(jì)帶城市工業(yè)綠色發(fā)展效率及其空間驅(qū)動(dòng)機(jī)制研究[J].中國(guó)人口·資源與環(huán)境，2019，29（8）：40-49..

[15] 丁顯有，肖雯，田澤.長(zhǎng)三角城市群工業(yè)綠色創(chuàng)新發(fā)展效率及其協(xié)同效應(yīng)研究[J].工業(yè)技術(shù)經(jīng)濟(jì)，2019，38（7）：67-75.

[16] 王韶華，劉曄，張偉，等.基于粗糙集與未確知模型的京津冀工業(yè)綠色發(fā)展的差異性與協(xié)調(diào)性分析[J].運(yùn)籌與管理，2021，30（5）：200-207.

[17] 陳超凡.中國(guó)工業(yè)綠色全要素生產(chǎn)率及其影響因素——基于ML生產(chǎn)率指數(shù)及動(dòng)態(tài)面板模型的實(shí)證研究[J].統(tǒng)計(jì)研究，2016，33（3）：53-62.

[18] 杜永強(qiáng)，遲國(guó)泰.基于指標(biāo)甄別的綠色產(chǎn)業(yè)評(píng)價(jià)指標(biāo)體系構(gòu)建[J].科研管理，2015，36（9）：119-127.

[19] 何小鋼，王自力.能源偏向型技術(shù)進(jìn)步與綠色增長(zhǎng)轉(zhuǎn)型——基于中國(guó)33個(gè)行業(yè)的實(shí)證考察[J].中國(guó)工業(yè)經(jīng)濟(jì)，2015（2）：50-62.

[20] 曾剛，陸琳憶，何金廖.生態(tài)創(chuàng)新對(duì)資源型城市產(chǎn)業(yè)結(jié)構(gòu)與工業(yè)綠色效率的影響[J].資源科學(xué)，2021，41（3）：94-103.

[21] 張惠程.國(guó)內(nèi)外資源型城市綠色轉(zhuǎn)型研究綜述[J].資源與產(chǎn)業(yè)，2018，20（5）：9-15.

[22] 鄧宗兵，宗樹偉，蘇聰文，等.長(zhǎng)江經(jīng)濟(jì)帶生態(tài)文明建設(shè)與新型城鎮(zhèn)化耦合協(xié)調(diào)發(fā)展及動(dòng)力因素研究[J].經(jīng)濟(jì)地理，2019，39（10）：78-86.

[23] 段永蕙，景建邦，張乃明.山西省人口、資源環(huán)境與經(jīng)濟(jì)協(xié)調(diào)發(fā)展分析[J].生態(tài)經(jīng)濟(jì)，2017，33（4）：64-68.

[24] 王毅，丁正山，余茂軍，等.基于耦合模型的現(xiàn)代服務(wù)業(yè)與城市化協(xié)調(diào)關(guān)系量化分析——以江蘇省常熟市為例[J].地理研究，2015，34（1）：94-108.

[25] 張榮天，焦華富.中國(guó)省際城鎮(zhèn)化與生態(tài)環(huán)境的耦合協(xié)調(diào)與優(yōu)化探討[J].干旱區(qū)資源與環(huán)境，2015，29（7）：12-17.

[26] 劉衡，朱鐵輝，辛嶺.農(nóng)業(yè)現(xiàn)代化與城鄉(xiāng)融合的耦合協(xié)調(diào)和驅(qū)動(dòng)因素[J].農(nóng)業(yè)現(xiàn)代化研究，2021，42（6）：1-14.

[27] 陳國(guó)生，蕭烽，黃鑫.湖南農(nóng)村人力資本與農(nóng)業(yè)現(xiàn)代化耦合協(xié)調(diào)發(fā)展[J].經(jīng)濟(jì)地理，2020，40（10）：176-182.

[28] 雷勛平，邱廣華.基于熵權(quán)TOPSIS模型的區(qū)域資源環(huán)境承載力評(píng)價(jià)實(shí)證研究[J].環(huán)境科學(xué)學(xué)報(bào)，2016，36（1）：314-323.

[29] 刁莉娟，姚建，艾怡凝.基于SPA-TOPSIS耦合的工業(yè)綠色發(fā)展綜合評(píng)價(jià)[J].生態(tài)經(jīng)濟(jì)，2020，36（9）：54-57.

[30] 李曉鐘，高嬌嬌.中國(guó)工業(yè)綠色發(fā)展效率時(shí)空差異比較研究[J].杭州電子科技大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2020，16（4）：1-8.

[31] 趙洋.我國(guó)資源型城市產(chǎn)業(yè)綠色轉(zhuǎn)型效率研究——基于地級(jí)資源型城市面板數(shù)據(jù)實(shí)證分析[J].經(jīng)濟(jì)問(wèn)題探索，2019（7）：94-101.

[32] 孫菲，楊曉龍，陳雍淦，等.東北地區(qū)資源型城市綠色發(fā)展效率測(cè)度——基于2008—2017年數(shù)據(jù)的實(shí)證研究[J].中國(guó)石油大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2021，37（3）：18-25.

[33] 李燦，張鳳榮，朱泰峰，等.基于熵權(quán)TOPSIS模型的土地利用績(jī)效評(píng)價(jià)及關(guān)聯(lián)分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2013（5）：217-227.

[34] Mohtadi H. Environment，Growth and Optimal Policy Design[J]. Journal of Public Economics， 1996，63（1）：119-140.

[35] 藏媛.西部地區(qū)綠色發(fā)展評(píng)價(jià)及影響因素研究[D].蘭州：蘭州財(cái)經(jīng)大學(xué)，2019.

[36] 董旭，吳傳清.中國(guó)城市全要素生產(chǎn)率的時(shí)空演變與影響因素研究——來(lái)自35個(gè)主要城市2000—2014年的經(jīng)驗(yàn)證據(jù)[J].學(xué)習(xí)與實(shí)踐，2017（5）：5-16.

[37] 李劍波，鮮學(xué)福.綠色低碳約束下重慶能源DEA效率實(shí)證分析[J].重慶大學(xué)學(xué)報(bào)，2017，40（2）：26-36.

[38] 張勝男.人工智能與能源產(chǎn)業(yè)的融合發(fā)展研究[D].徐州：中國(guó)礦業(yè)大學(xué)， 2020.

[39] 汪東芳，曹建華.互聯(lián)網(wǎng)發(fā)展對(duì)中國(guó)全要素能源效率的影響及網(wǎng)絡(luò)效應(yīng)研究[J].中國(guó)人口·資源與環(huán)境，2019，29（1）：86-95.

責(zé)任編輯：曲 紅

Study on Coupling Coordination Effects of Industrial Green Development

Level and Industrial Green Development Efficiency of Oil and Gas Cities

SUN Fei1， LI Youjun1， GONG Yanfen1， PENG Min2

（1.School of Economics and Management， Northeast Petroleum University， Daqing 163318，? Heilongjiang， China;

2.Institute of Petroleum Economics and Management， Northeast Petroleum University， Daqing 163318， Heilongjiang， China）

Abstract： The paper constructs the respective evaluation index systems in light of both the level and the efficiency of the industrial green development， measures the level and the efficiency by using the entropy-TOPSIS model and SBM model， and studies the effects by applying the coupling coordination model in oil and gas cities. The results show： The coupling coordination effects of industrial green development in various regions are apparent with the gradient increasing pattern over time. The two systems show an overall trend with greater intensity from the year 2004 to 2018 in the oil and gas cities. There is a trend of gap decrease in coupling coordination effects over time in different cities， an obvious type classification and a gap in effects depending on the city concerned. It can be concluded that quality talents， advanced technology， regional differences， environmental regulations， etc. are the primary factors to the coupling coordination effect of two systems. In the final analysis， the paper proposes multiple strategies to boost industrial green development of oil and gas city： optimizing the talent cultivation mechanism， increasing the investment in green technology innovation， building coordinated development mechanism and so on.

Keywords： oil and gas cities;industrial green development level; industrial green development efficiency; coupling coordination degree