孫菲 李友俊 鞏艷芬 彭民

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

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

中圖分類號：F124

文獻標識碼：A

文章編號：1673-5595（2022）02-0027-13

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

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

一、文獻綜述

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

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

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

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

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

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

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

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

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

工業(yè)綠色發(fā)展效率是決定工業(yè)綠色發(fā)展的內(nèi)生動力。從宏觀層面看，工業(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ā)展效率的互補發(fā)展程度，反映工業(yè)綠色發(fā)展水平對工業(yè)綠色發(fā)展效率的支撐作用及工業(yè)綠色發(fā)展效率對工業(yè)綠色發(fā)展水平的引致效應(yīng)。一般采用耦合協(xié)調(diào)度模型測度工業(yè)綠色發(fā)展耦合協(xié)調(diào)效應(yīng)。油氣城市工業(yè)綠色發(fā)展進程中，工業(yè)綠色發(fā)展水平和工業(yè)綠色發(fā)展效率兩大子系統(tǒng)互為前提，且相輔相成，二者與工業(yè)綠色發(fā)展的關(guān)系如圖1所示。

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

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

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

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

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

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

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

（一）評價指標體系構(gòu)建

1.油氣城市工業(yè)綠色發(fā)展水平評價指標體系

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

2.油氣城市工業(yè)綠色發(fā)展效率評價指標體系

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

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

（二）綜合評價模型

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

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

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

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

X=X11…X1n

Xm1…Xmnm×n（1）

（1）數(shù)據(jù)標準化。

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

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

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

負向指標（越小越優(yōu)型）為：

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

（2）計算各項指標的比重。

第j項指標下第i個樣本對象占該指標的比重為：

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

（3）計算第j項指標的熵值。

信息熵值為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）計算熵權(quán)。

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

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

利用熵權(quán)法確定的指標權(quán)重向量wj來計算加權(quán)規(guī)范化決策矩陣，即R=（Rij）m×n通過矩陣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）確定正、負理想解。

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

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

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

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

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

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

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

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

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

（8）計算相對貼近度。

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

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

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

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

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

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

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ī)模報酬不變前提下，可定義生產(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-為投入的松弛變量，表示投入過剩;sg、sb為產(chǎn)出的松弛變量，表示產(chǎn)出不足和非期望產(chǎn)出過多。目標函數(shù)ρ′關(guān)于s-、sg、sb嚴格遞減，且0≤ρ′≤1。對于特定的DMU，當(dāng)ρ′=0時，說明該DMU是無效的;當(dāng)ρ′=1，s-=sg=sb=0時，該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)性越強。

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

在某些情況下，耦合度雖能反映工業(yè)綠色發(fā)展水平和效率的作用強度和方向，但無法準確反映二者的整體功效和耦合協(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ā)展水平與效率這兩個子系統(tǒng)之間的耦合協(xié)調(diào)關(guān)系，故認為各子系統(tǒng)對油氣城市的工業(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，說明兩子系統(tǒng)相互作用、相互影響的效果越顯著，兩子系統(tǒng)之間的耦合協(xié)調(diào)關(guān)系越強;當(dāng)D=1時，兩子系統(tǒng)達到良性互動的最優(yōu)協(xié)調(diào)發(fā)展階段;D=0，表明工業(yè)綠色發(fā)展水平和效率兩子系統(tǒng)之間沒有必然聯(lián)系，對彼此的現(xiàn)狀和未來發(fā)展沒有影響。[38]

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

耦合度和耦合協(xié)調(diào)度能夠表征兩子系統(tǒng)的耦合與協(xié)調(diào)關(guān)系，依據(jù)中值分段原理，對耦合度的等級進行劃分，[38-39]標準如表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ā)展階段低水平耦合頡頏階段初步磨合階段中級磨合階段深度磨合階段高水平耦合階段

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

五、結(jié)論與建議

（一）結(jié)論

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

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

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

（二）建議

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

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

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

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

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

參考文獻：

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

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

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

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

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

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

[7] 徐成龍，莊貴陽.供給側(cè)改革驅(qū)動中國工業(yè)綠色發(fā)展的動力結(jié)構(gòu)及時空效應(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] 王建民，仇定三，蔣倩穎，等.長江經(jīng)濟帶工業(yè)綠色發(fā)展效率測量與提升路徑研究[J].科技管理研究，2019，39（12）：46-52.

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

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

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

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

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

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

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

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

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

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

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

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

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

[25] 張榮天，焦華富.中國省際城鎮(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ū)動因素[J].農(nóng)業(yè)現(xiàn)代化研究，2021，42（6）：1-14.

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

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

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

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

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

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

[33] 李燦，張鳳榮，朱泰峰，等.基于熵權(quán)TOPSIS模型的土地利用績效評價及關(guān)聯(lián)分析[J].農(nóng)業(yè)工程學(xué)報，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ā)展評價及影響因素研究[D].蘭州：蘭州財經(jīng)大學(xué)，2019.

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

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

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

[39] 汪東芳，曹建華.互聯(lián)網(wǎng)發(fā)展對中國全要素能源效率的影響及網(wǎng)絡(luò)效應(yīng)研究[J].中國人口·資源與環(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