基于系統(tǒng)動力學模型的長江經濟帶綠色發(fā)展?jié)摿ρ芯?/h1>
2023-03-14 10:33:14于涵如
統(tǒng)計理論與實踐 2023年1期
袁 靖 馬 騰 于涵如 彭 燕
（山東工商學院 統(tǒng)計學院，山東 煙臺 264005）
一、引言
面對全球生態(tài)問題和經濟復蘇壓力，綠色發(fā)展已經成為世界發(fā)展的重要趨勢。黨的二十大報告指出未來中國要大力推動綠色發(fā)展。由此可見，對中國各省（區(qū)、市）未來綠色發(fā)展?jié)摿M行綜合評價及最優(yōu)情景設定具有重要研究意義。
長江流域是中華民族重要的起源地，長江經濟帶以長江流域為基底，是中國生態(tài)地位突出且發(fā)展?jié)摿薮蟮慕洕鷧^(qū)域，走生態(tài)優(yōu)先、綠色發(fā)展道路是其發(fā)展的應有之義。長江經濟帶劃分為上游、中游、下游，下游地區(qū)包括上海、江蘇、浙江、安徽4省（市），面積35.03萬平方公里，占長江經濟帶總面積的17.1%；中游地區(qū)包括江西、湖北、湖南3省，面積56.46萬平方公里，占長江經濟帶的27.5%；上游地區(qū)包括重慶、四川、貴州、云南4?。ㄊ校?，面積113.74萬平方公里，占長江經濟帶的55.4%。隨著長江經濟帶的不斷開發(fā)，資源環(huán)境生態(tài)已經成為制約長江經濟帶發(fā)展的重要瓶頸，因而對其未來經濟發(fā)展?jié)摿M行綜合評價尤為重要。
國際學術界多圍繞人類社會發(fā)展和自然生態(tài)保護之間的關系來研究綠色發(fā)展，整體看，研究內容多包括綠色能源優(yōu)化、實現綠色發(fā)展以及綠色經濟轉型發(fā)展策略等方面。從綠色發(fā)展評價測度看，研究多側重基于不同方法構建綠色發(fā)展評價指標體系測度研究對象的綠色發(fā)展水平。隨著研究的不斷深入，測度方法呈現多樣化，包括指標賦權法、熵值法、數據包絡分析、神經網絡模型、能值法、網絡層次模型、生態(tài)足跡法等。近年來，學者對中國綠色發(fā)展評價的研究較多，包括對中國省際綠色發(fā)展水平及影響因素進行研究，以及對黃河流域、京津冀等經濟帶的綜合綠色發(fā)展水平進行評價。
學者對中國綠色發(fā)展的內涵、測度等研究日趨完善，但是對中國綠色發(fā)展?jié)摿Φ难芯窟€不夠，現有的綠色發(fā)展評價研究側重如何識別并快速提升綠色發(fā)展短板，不具備發(fā)掘深層次促進或制約區(qū)域發(fā)展?jié)摿Φ哪芰Γ诰G色發(fā)展路徑優(yōu)化方面存在先天性缺陷，而關于綠色發(fā)展實現路徑的研究則側重于某些指標的優(yōu)化，忽視了區(qū)域發(fā)展的系統(tǒng)性和綜合性?？茖W判斷綠色發(fā)展?jié)摿Σ⒅D化為現實競爭力才是國家及地區(qū)決策部門需要關注的重點，本文將綠色發(fā)展指數與系統(tǒng)動力學模型相結合，提出綠色發(fā)展?jié)摿Φ母拍詈投垦芯糠桨?，以長江經濟帶11?。ㄊ校┑娜娼洕l(fā)展數據為基礎，一方面能夠更加綜合地評判長江經濟帶綠色發(fā)展?jié)摿?，為未來長江經濟帶綠色發(fā)展提供有益借鑒；另一方面通過分析綠色發(fā)展系統(tǒng)內部指標間的制約機制，探索提升綠色發(fā)展水平的有效方案，擴展了綠色發(fā)展的研究范圍。
二、構建長江經濟帶11省（市）綠色發(fā)展評價指標體系
（一）綠色發(fā)展評價指標體系選取原則
1.全面性原則。構建綠色發(fā)展評價指標體系的目的在于綜合反映區(qū)域綠色發(fā)展水平，綠色發(fā)展的內涵包括經濟、社會、資源、生態(tài)環(huán)境等多方面因素，因此需要保證指標選取的全面性，將綠色發(fā)展相關變量盡可能多地納入指標體系。
2.差異化原則。各地區(qū)自然資源稟賦、發(fā)展戰(zhàn)略定位以及發(fā)展現狀不盡相同，導致各地區(qū)綠色發(fā)展的關注重點有所差別，構建區(qū)域綠色發(fā)展評價指標體系需要結合區(qū)域發(fā)展特點選取差異化的評價指標。
3.可操作性原則。可操作性主要表現在指標量化與數據獲取的方便可靠、計算結果的可比等方面。另外，本文結合了綠色發(fā)展指數與系統(tǒng)動力學模型兩種模型工具，因此還考慮了選取指標參與構建綠色發(fā)展系統(tǒng)動力學模型的融合性。
（二）構建指標體系
本文在保證指標的全面性、差異性、可操作性基礎上，同時兼顧代表性、簡潔性，評價指標由資源利用、環(huán)境治理、環(huán)境質量、生態(tài)保護、增長質量、綠色生活6個方面體現（見表1）。
表1 長江經濟帶綠色發(fā)展評價指標體系
（續(xù)表）
（三）指標權重確定
常用的指標賦權方法有德爾菲法、層次分析法、熵值法、主成分分析法、因子分析法等，其中，熵值法是利用信息熵來判斷指標有效性和價值的一種客觀賦權方法。采用熵值法確定指標權重，能夠深刻反映指標信息熵值的效用價值，既可以克服主觀賦權法無法避免的隨機性、臆斷性問題，還可以有效解決多指標變量間的信息重疊問題。
本文采用熵值法確定指標權重，計算步驟如下：
1.指標標準化
由于各指標量綱不一致，不能直接進行比較，正向指標參考式（1），逆向指標參考式（2）：
式中，Xij為評價指標；X*ij為評價指標標準化值；i為指標代碼編號；j為年份。
2.計算指標值在所屬一級指標中的占比
3.計算信息熵與信息熵冗余度
式中，K是與樣本數量相關的常數；m為樣本量；ei為第i項指標的信息熵；di為信息熵冗余度。
4.計算指標權重
熵值法計算指標權重的本質是比較指標的信息熵效用價值，因此采用信息熵冗余度的占比來代替指標權重，參考式（7），最終采用式（8）合成全部二級指標權重。
式中，wi表示第i項指標在綠色發(fā)展整體系統(tǒng)中的權重；w'i表示第i項指標在所屬一級指標中的權重；ej為對應的一級指標權重。
三、長江經濟帶11?。ㄊ校┚G色發(fā)展指數評價
（一）數據來源
經濟社會發(fā)展相關數據來源于相關年份長江經濟帶11?。ㄊ校┙y(tǒng)計年鑒，其中經濟數據以2013年為不變價進行換算；能源利用數據來源于相關年份《中國能源統(tǒng)計年鑒》；廢水、廢氣、固體廢物排放數據來源于相關年份長江經濟帶11?。ㄊ校┉h(huán)境狀況公報，其他生態(tài)環(huán)境相關數據來源于相關年份《中國環(huán)境統(tǒng)計年鑒》。
（二）計算綠色發(fā)展指數
采用綜合加權法計算綠色發(fā)展指數：
式中，Gj表示j年綠色發(fā)展指數。
綠色發(fā)展指數從多方面綜合反映了區(qū)域綠色發(fā)展水平，同時也是進行綠色發(fā)展?jié)摿α炕u價的基礎。按照上述計算過程得到中國長江經濟帶11?。ㄊ校?013—2018年綠色發(fā)展指數，見表2。
表2 長江經濟帶11?。ㄊ校┚G色發(fā)展指數
（續(xù)表）
本文構建的綠色發(fā)展指標體系包含資源利用、環(huán)境治理、環(huán)境質量、生態(tài)保護、增長質量、綠色生活6個方面。其中資源利用重點關注了水資源、土地資源、能源資源的利用效率以及可回收利用資源的再利用情況；環(huán)境治理是綠色發(fā)展的核心領域之一，反映政府和社會對環(huán)境污染治理的資金投入力度；環(huán)境質量是環(huán)境治理的結果體現；生態(tài)保護體現了自然生態(tài)質量與生態(tài)保護投入情況；增長質量主要表現為經濟總量增長以及經濟發(fā)展的可持續(xù)性；綠色生活主要關注居民生活方式的綠色性。
從長江經濟帶11?。ㄊ校┑木G色發(fā)展指數計算結果，可以得到以下結論：（1）2013年11?。ㄊ校┑木G色發(fā)展指數存在一定差異，但經過近年來綠色發(fā)展理念的貫徹執(zhí)行，到2018年11省（市）的綠色發(fā)展之間的差異明顯縮小，11?。ㄊ校┚G色發(fā)展指數都超過0.90，說明近年來長江經濟帶11?。ㄊ校┳⒅貐f(xié)同發(fā)展。（2）總體看，針對資源利用、環(huán)境治理、環(huán)境質量、生態(tài)保護、增長質量、綠色生活6個方面指標，我們可以看到11省（市）的環(huán)境治理和增長質量兩個一級指標計算結果增長顯著，尤其是上游4個省（市），2013年重慶環(huán)境治理和增長質量兩個指標計算結果均為0.58，2018年兩個指標分別增長至0.84和0.80；四川兩個指標計算結果由2013年的0.61和0.57，增長至2018年的0.89和0.80；云南兩個指標計算結果由2013年的0.64和0.56，增長至2018年的0.84和0.80；貴州兩個指標計算結果由2013年的0.50和0.59，增長至2018年的0.93和0.80。（3）從資源利用角度看，下游與中游7個?。ㄊ校┍壬嫌?個?。ㄊ校┲笜烁?，上游的重慶、四川、貴州和云南資源豐富，但是未充分開發(fā)利用。（4）環(huán)境治理方面，貴州、江西和湖北近年來表現突出。近年來長江兩岸防護林體系建設造林超千萬畝，通過國家濕地公園、綠色發(fā)展示范區(qū)以及退化濕地恢復等建設，長江經濟帶重點湖庫及流域水質優(yōu)良斷面比例顯著提高，水源涵養(yǎng)功能明顯增強。禁漁退漁在長江流域重點水域取得階段性成效，18萬漁民全部提前退捕上岸。禁漁以來，長江水生物資源逐漸恢復，如作為長江流域生態(tài)環(huán)境“指標生物”的江豚在湖北宜昌長江三江段頻現。（5）由于環(huán)境治理成效顯著，長江經濟帶11?。ㄊ校?013—2018年環(huán)境質量大幅提升。（6）從生態(tài)保護角度看，2013—2016年上海生態(tài)保護措施不到位，但是2017年之后出臺相關政策，進而有了大幅提升。（7）從綠色生活角度看，其他省市都有小幅上升，上海2018年反而下降，說明上海需要特別注意居民生活方式的綠色改善。
經濟的綠色發(fā)展水平是衡量綠色發(fā)展的重要測度指標，今后長江經濟帶11?。ㄊ校┙洕l(fā)展的重要內容仍然是促進經濟結構優(yōu)化升級，推動經濟發(fā)展綠色轉型升級，培育新的經濟增長點，拓展綠色發(fā)展空間，不斷提高經濟發(fā)展的質量和效益。
四、綠色發(fā)展系統(tǒng)動力學模型構建
將長江經濟帶區(qū)域發(fā)展系統(tǒng)劃分為經濟、社會、資源、生態(tài)環(huán)境4個子系統(tǒng)，通過分析子系統(tǒng)和變量間的因果關系與數學關系，開發(fā)設計出長江經濟帶11?。ㄊ校?013—2060年綠色發(fā)展仿真模型并進行檢驗，檢驗結果表明模型運行穩(wěn)健，與系統(tǒng)真實運行情況基本相符，能夠用于長江經濟帶11?。ㄊ校┚G色發(fā)展仿真模擬。
（一）系統(tǒng)邊界及基本假設
2020年9月，習近平主席在第七十五屆聯合國大會一般性辯論上提出，中國將采取更加有力的政策和措施，二氧化碳排放力爭于2030年前達到峰值，努力爭取2060年前實現碳中和。因此本文設定模型有效時間為2013—2060年，其中2013—2018年為模型設定與檢驗區(qū)間，2018—2060年為模型預測區(qū)間；設定2018—2030年為中期預測區(qū)間，2030—2060年為長期預測區(qū)間。設定時間間隔為1年。
一般而言，系統(tǒng)邊界內包含與研究問題有關的全部對象，但由于區(qū)域綠色發(fā)展系統(tǒng)強調經濟、社會、資源、生態(tài)環(huán)境4大子系統(tǒng)的整體與協(xié)調發(fā)展，涉及的指標數量眾多，較難逐一列舉。本文從各子系統(tǒng)的內涵出發(fā)，分別選擇與綠色發(fā)展關聯度較高的變量納入模型：（1）經濟子系統(tǒng)中選擇GDP、產業(yè)結構、固定資產投資、財政支出、研發(fā)支出、居民收入等變量；（2）社會子系統(tǒng)中選擇人口、城鎮(zhèn)化率、就業(yè)、醫(yī)療衛(wèi)生水平、公共交通情況等變量；（3）資源子系統(tǒng)中選擇能耗、水耗、土地資源消耗、碳排放等變量；（4）生態(tài)環(huán)境子系統(tǒng)中選擇水、大氣、固體廢物等污染物排放水平及生態(tài)環(huán)境質量、污染治理投入等變量。
（二）系統(tǒng)動力學流圖
因果關系圖只能反映指標間的相關性和反饋過程，需要進一步繪制系統(tǒng)動力學流圖來實現對各要素運行情況的仿真模擬。長江經濟帶11省（市）綠色發(fā)展系統(tǒng)由經濟、社會、資源、生態(tài)環(huán)境4個子系統(tǒng)構成，得到系統(tǒng)動力學流圖①系統(tǒng)動力學流圖可咨詢作者或《統(tǒng)計理論與實踐》編輯部。。
（三）模型方程式
模型方程式定量化描述指標間因果關系。綠色發(fā)展系統(tǒng)動力學模型中涉及的變量類型有水平變量、速率變量、輔助變量以及常量，各變量之間的定量關系通過模型方程予以描述，共建立95個模型方程。
五、長江經濟帶11?。ㄊ校┚G色發(fā)展?jié)摿η榫胺治?/h2>（一）綠色發(fā)展情景設定
情景分析是政策制定的基礎，通過改變不同控制指標的參數，設定多個發(fā)展情景并進行模擬仿真，分析不同參數變化情況對綠色發(fā)展系統(tǒng)的影響，從而為提升區(qū)域綠色發(fā)展?jié)摿μ峁┖侠斫ㄗh。
情景分析一般分為單要素模擬與多要素綜合模擬，單要素模擬主要分析單個控制指標變動對系統(tǒng)的影響，多要素模擬分析多個控制指標同時變動對系統(tǒng)的影響?？紤]到單一指標調整對綠色發(fā)展系統(tǒng)影響較小，從產業(yè)結構優(yōu)化、能源結構優(yōu)化、固定資產投資管控、科技創(chuàng)新資金投入、污染治理資金投入等角度選取多個指標，設定不同的組合方案進行多要素模擬。具體調控指標包括第二產業(yè)增加值占比、第三產業(yè)增加值占比、工業(yè)增加值占比、GDP增速、煤炭消費占比、天然氣消費占比、財政支出占比、環(huán)境污染治理投資占比、工業(yè)固體廢物產污系數、工業(yè)源COD排污系數、工業(yè)固體廢物綜合利用率以及研發(fā)支出增長率，這些指標在系統(tǒng)動力學模型中均以表函數表示。
P1（基礎模式，不改變現有參數）：作為基礎對照組，情景方案設置延續(xù)現狀發(fā)展路徑，各項調控指標參數保持現有變化趨勢，不進行調整。
P2（產業(yè)結構優(yōu)化模式）：第二產業(yè)增加值占比減少10%，第三產業(yè)增加值占比增加10%，工業(yè)增加值占比減少10%，GDP增速變?yōu)?.2倍。
P3（能源結構優(yōu)化模式）：煤炭消費占比減少15%，天然氣消費占比增加5%，財政支出占比變?yōu)?.2倍。
P4（技術進步模式）：研發(fā)支出增長率變?yōu)?.2倍。
P5（環(huán)境治理模式）：環(huán)境污染治理投資占比和工業(yè)固體廢物綜合利用率變?yōu)?.2倍，工業(yè)固體廢物產污系數和工業(yè)源COD排污系數變?yōu)?.8倍。
P6（綜合模式）：同時進行P2、P3、P4、P5的4種調整方案。
由于長江經濟帶11?。ㄊ校┠Ｐ蛿祿跏贾荡嬖诓町?，因而僅列出上海市綠色發(fā)展?jié)摿η榫澳M調控指標的參數設置，如表3所示。其余各?。ㄊ校┡c之相似，但依賴于具體模型數據初始值。
表3 上海市系統(tǒng)動力學模型情景模擬調控指標的參數設置結果
（二）長江經濟帶11省（市）綠色發(fā)展?jié)摿η榫澳M對比分析
將長江經濟帶11?。ㄊ校┚G色發(fā)展?jié)摿?種情景模擬結果進行對比，中期設置為2018—2030年，長期設置為2030—2060年，結果見表4、表5。
據表4、表5的模擬結果，可以看到：（1）長江經濟帶11?。ㄊ校╅L期的綠色發(fā)展指數變化幅度均顯著于中期變化幅度，上海長期綠色發(fā)展指數變化幅度均達到138以上，貴州長期綠色發(fā)展指數變化幅度在第六種情景模式中也達到78.256。（2）不論是中期變化幅度還是長期變化幅度，上海的變化幅度均位列第一，接著是上游4的重慶、四川、貴州、云南，其次是下游的浙江、安徽和中游的湖南，變化幅度最小的是中游的江西、湖北和下游的江蘇。
表4 2018—2030年長江經濟帶11?。ㄊ校┚G色發(fā)展指數變化幅度模擬結果
表5 2030—2060年長江經濟帶11?。ㄊ校┚G色發(fā)展指數變化幅度模擬結果
根據11?。ㄊ校┚G色發(fā)展指數差異，可以得出以下結論：（1）上海作為長三角一體化發(fā)展及長江經濟帶發(fā)展的龍頭，地處流域出?？冢仁且I流域整體轉型發(fā)展的橋頭堡，也是這一區(qū)域與海外市場的重要樞紐。因此，上海未來的綠色發(fā)展?jié)摿薮蟆＃?）長江上游4個?。ㄊ校ㄖ貞c、四川、貴州、云南）位于中國西南部，是長江上游地區(qū)的重要生態(tài)屏障和經濟、金融、科創(chuàng)、航運和商貿物流中心，是長江經濟帶上游成渝城市群的核心城市，其生態(tài)環(huán)境位居全國前列，未來的綠色發(fā)展?jié)摿σ埠艽?。?）近年來長江下游和中游的江蘇、浙江、安徽、湖北、湖南和江西工業(yè)規(guī)?？偭看?、增速快，重化工業(yè)占比偏大，第二產業(yè)發(fā)展對區(qū)域資源環(huán)境造成的負荷仍然較高，煤炭消費量占比較大，清潔能源消費占比不高，綜合看，其未來綠色發(fā)展?jié)摿π枰M一步挖掘。
根據表4、表5列出的6種情景模擬結果差異對比，得出以下結論：
（1）作為基礎對照方案，該情景下綠色發(fā)展指數穩(wěn)定提升，綠色發(fā)展水平存在一定上升空間，綠色發(fā)展?jié)摿^高。但是對于長江經濟帶11?。ㄊ校?，基礎模式下的綠色發(fā)展指數提升幅度均是6種模式下最小的，因而不是最優(yōu)發(fā)展方案。
（2）對比產業(yè)結構優(yōu)化模式、能源結構優(yōu)化模式、技術進步模式和環(huán)境治理模式，產業(yè)結構優(yōu)化模式下長江經濟帶11?。ㄊ校┚G色發(fā)展指數變化幅度最大，其次是技術進步模式，再次是環(huán)境治理模式，最后是能源結構優(yōu)化模式。
（3）產業(yè)結構優(yōu)化模式情景設置是以提高第三產業(yè)和高新技術產業(yè)產值占比，降低工業(yè)增加值占比為主要出發(fā)點，并協(xié)同帶動經濟增長。從結果看，該情景下產業(yè)結構明顯優(yōu)化，經濟增長動能增強，主要是因為資源利用效率快速提升，工業(yè)污染物排放量明顯降低，對綠色發(fā)展?jié)摿μ嵘哂忻黠@促進效果。
（4）技術進步模式下，科技創(chuàng)新能力提升對經濟增長質量、清潔生產水平、環(huán)境污染治理等多方面產生正向促進作用，與前述發(fā)展方案相比較，該情景下各類指標指數均有一定幅度增長。
（5）環(huán)境治理模式下情景設置旨在通過加強環(huán)境污染治理投資、工業(yè)污染管控水平來降低污染物排放，綠色發(fā)展水平大幅度提升。
（6）能源結構優(yōu)化模式情景設定目的在于優(yōu)化能源結構、提升能源利用效率。從結果看，碳排放削減效果明顯，符合情景設定目的，對綠色發(fā)展綜合指數提升具有一定促進效果，但指數增長幅度不大，表明僅通過優(yōu)化能源利用方式已較難快速提升區(qū)域綠色發(fā)展?jié)摿Γ摲桨高m合作為其他方案的補充。
（7）綜合模式情景從資金管控和政策引導兩個方面進行控制，并綜合考慮到產業(yè)結構、能源結構、污染治理、科技創(chuàng)新等方面的因素進行調整。該情景方案下，長江經濟帶11省（市）綠色發(fā)展指數增長幅度最大，采用多項措施的組合方案對提升區(qū)域綠色發(fā)展質量效果最突出，是提升綠色發(fā)展?jié)摿Φ淖顑?yōu)選擇。
總結上述研究結果認為，加強產業(yè)結構轉型和生態(tài)環(huán)境治理投入是提升長江經濟帶11?。ㄊ校┚G色發(fā)展?jié)摿Φ闹攸c所在，應當以節(jié)能降耗減污為政策導向，通過優(yōu)化產業(yè)結構、提高環(huán)境污染治理資金投入，并配合實施能源結構調整、加大技術創(chuàng)新等一系列組合措施，最大程度促進長江經濟帶11省（市）高質量綠色發(fā)展。
六、結論
本文通過構建綠色發(fā)展評價指標體系，得到了反映綠色發(fā)展綜合水平的綠色發(fā)展指數，并設計了包括經濟、社會、資源、生態(tài)環(huán)境4個子系統(tǒng)的長江經濟帶11?。ㄊ校┑木G色發(fā)展系統(tǒng)動力學模型，模型檢驗的可靠性和穩(wěn)定性較好，能夠反映綠色發(fā)展各關鍵指標變化情況。在此基礎上，從優(yōu)化產業(yè)結構、能源清潔利用、促進科技創(chuàng)新、強化污染治理等方面設定不同的組合方案，仿真模擬了長江經濟帶11省（市）2018—2060年綠色發(fā)展變化趨勢。研究結果表明：（1）根據綠色發(fā)展指數的計算結果，近年來長江經濟帶11?。ㄊ校┚G色發(fā)展均取得明顯成效，且各一級指標指數均有提升。（2）長江經濟帶11省（市）未來發(fā)展，仍需要持續(xù)推動產業(yè)結構轉型升級。（3）科技創(chuàng)新為經濟社會發(fā)展提供重要支撐，不僅能夠提升勞動生產率，還能提升資源環(huán)境產出率。長江經濟帶11?。ㄊ校碛休^好的科技創(chuàng)新土壤，人才科教資源十分豐富，應大力增加綠色技術研發(fā)投入，促進綠色產業(yè)發(fā)展，進一步凸顯綠色技術在優(yōu)化產業(yè)結構、提升資源利用效率、促進污染減排等方面的作用。（4）未來長江經濟帶11?。ㄊ校^續(xù)加大環(huán)境污染治理資金投入力度，削減污染物排放總量。（5）長江經濟帶11省（市）可以根據自身實際情況提升未來綠色發(fā)展?jié)摿??！?/p>
猜你喜歡
一條江的嬗變長江經濟帶綠色發(fā)展之路法人(2021年2期)2021-03-02 07:12:46 長江之頭青年歌聲(2020年11期)2020-11-24 06:57:28 今日農業(yè)(2019年11期)2019-08-13 00:49:02 陜西呼應長江經濟帶當代陜西(2018年9期)2018-11-18 07:26:13 長江之歌(外二首)天津詩人(2017年2期)2017-11-29 01:24:30 長江圖(外二首)天津詩人(2017年2期)2017-11-29 01:24:05 領導決策信息(2017年17期)2017-06-21 09:51:22 冷藏保溫車發(fā)展?jié)摿Ρ患ぐl(fā)專用汽車(2016年5期)2016-03-01 04:14:39 汽車拆解市場的發(fā)展?jié)摿⒅饾u釋放中國資源綜合利用(2016年7期)2016-02-03 03:00:14 長江經濟帶與漢江生態(tài)經濟帶如何協(xié)調融合學習月刊(2015年9期)2015-07-09 05:33:46
袁 靖 馬 騰 于涵如 彭 燕
（山東工商學院 統(tǒng)計學院，山東 煙臺 264005）
一、引言
面對全球生態(tài)問題和經濟復蘇壓力，綠色發(fā)展已經成為世界發(fā)展的重要趨勢。黨的二十大報告指出未來中國要大力推動綠色發(fā)展。由此可見，對中國各省（區(qū)、市）未來綠色發(fā)展?jié)摿M行綜合評價及最優(yōu)情景設定具有重要研究意義。
長江流域是中華民族重要的起源地，長江經濟帶以長江流域為基底，是中國生態(tài)地位突出且發(fā)展?jié)摿薮蟮慕洕鷧^(qū)域，走生態(tài)優(yōu)先、綠色發(fā)展道路是其發(fā)展的應有之義。長江經濟帶劃分為上游、中游、下游，下游地區(qū)包括上海、江蘇、浙江、安徽4省（市），面積35.03萬平方公里，占長江經濟帶總面積的17.1%；中游地區(qū)包括江西、湖北、湖南3省，面積56.46萬平方公里，占長江經濟帶的27.5%；上游地區(qū)包括重慶、四川、貴州、云南4?。ㄊ校?，面積113.74萬平方公里，占長江經濟帶的55.4%。隨著長江經濟帶的不斷開發(fā)，資源環(huán)境生態(tài)已經成為制約長江經濟帶發(fā)展的重要瓶頸，因而對其未來經濟發(fā)展?jié)摿M行綜合評價尤為重要。
國際學術界多圍繞人類社會發(fā)展和自然生態(tài)保護之間的關系來研究綠色發(fā)展，整體看，研究內容多包括綠色能源優(yōu)化、實現綠色發(fā)展以及綠色經濟轉型發(fā)展策略等方面。從綠色發(fā)展評價測度看，研究多側重基于不同方法構建綠色發(fā)展評價指標體系測度研究對象的綠色發(fā)展水平。隨著研究的不斷深入，測度方法呈現多樣化，包括指標賦權法、熵值法、數據包絡分析、神經網絡模型、能值法、網絡層次模型、生態(tài)足跡法等。近年來，學者對中國綠色發(fā)展評價的研究較多，包括對中國省際綠色發(fā)展水平及影響因素進行研究，以及對黃河流域、京津冀等經濟帶的綜合綠色發(fā)展水平進行評價。
學者對中國綠色發(fā)展的內涵、測度等研究日趨完善，但是對中國綠色發(fā)展?jié)摿Φ难芯窟€不夠，現有的綠色發(fā)展評價研究側重如何識別并快速提升綠色發(fā)展短板，不具備發(fā)掘深層次促進或制約區(qū)域發(fā)展?jié)摿Φ哪芰Γ诰G色發(fā)展路徑優(yōu)化方面存在先天性缺陷，而關于綠色發(fā)展實現路徑的研究則側重于某些指標的優(yōu)化，忽視了區(qū)域發(fā)展的系統(tǒng)性和綜合性?？茖W判斷綠色發(fā)展?jié)摿Σ⒅D化為現實競爭力才是國家及地區(qū)決策部門需要關注的重點，本文將綠色發(fā)展指數與系統(tǒng)動力學模型相結合，提出綠色發(fā)展?jié)摿Φ母拍詈投垦芯糠桨?，以長江經濟帶11?。ㄊ校┑娜娼洕l(fā)展數據為基礎，一方面能夠更加綜合地評判長江經濟帶綠色發(fā)展?jié)摿?，為未來長江經濟帶綠色發(fā)展提供有益借鑒；另一方面通過分析綠色發(fā)展系統(tǒng)內部指標間的制約機制，探索提升綠色發(fā)展水平的有效方案，擴展了綠色發(fā)展的研究范圍。
二、構建長江經濟帶11省（市）綠色發(fā)展評價指標體系
（一）綠色發(fā)展評價指標體系選取原則
1.全面性原則。構建綠色發(fā)展評價指標體系的目的在于綜合反映區(qū)域綠色發(fā)展水平，綠色發(fā)展的內涵包括經濟、社會、資源、生態(tài)環(huán)境等多方面因素，因此需要保證指標選取的全面性，將綠色發(fā)展相關變量盡可能多地納入指標體系。
2.差異化原則。各地區(qū)自然資源稟賦、發(fā)展戰(zhàn)略定位以及發(fā)展現狀不盡相同，導致各地區(qū)綠色發(fā)展的關注重點有所差別，構建區(qū)域綠色發(fā)展評價指標體系需要結合區(qū)域發(fā)展特點選取差異化的評價指標。
3.可操作性原則。可操作性主要表現在指標量化與數據獲取的方便可靠、計算結果的可比等方面。另外，本文結合了綠色發(fā)展指數與系統(tǒng)動力學模型兩種模型工具，因此還考慮了選取指標參與構建綠色發(fā)展系統(tǒng)動力學模型的融合性。
（二）構建指標體系
本文在保證指標的全面性、差異性、可操作性基礎上，同時兼顧代表性、簡潔性，評價指標由資源利用、環(huán)境治理、環(huán)境質量、生態(tài)保護、增長質量、綠色生活6個方面體現（見表1）。
表1 長江經濟帶綠色發(fā)展評價指標體系
（續(xù)表）
（三）指標權重確定
常用的指標賦權方法有德爾菲法、層次分析法、熵值法、主成分分析法、因子分析法等，其中，熵值法是利用信息熵來判斷指標有效性和價值的一種客觀賦權方法。采用熵值法確定指標權重，能夠深刻反映指標信息熵值的效用價值，既可以克服主觀賦權法無法避免的隨機性、臆斷性問題，還可以有效解決多指標變量間的信息重疊問題。
本文采用熵值法確定指標權重，計算步驟如下：
1.指標標準化
由于各指標量綱不一致，不能直接進行比較，正向指標參考式（1），逆向指標參考式（2）：
式中，Xij為評價指標；X*ij為評價指標標準化值；i為指標代碼編號；j為年份。
2.計算指標值在所屬一級指標中的占比
3.計算信息熵與信息熵冗余度
式中，K是與樣本數量相關的常數；m為樣本量；ei為第i項指標的信息熵；di為信息熵冗余度。
4.計算指標權重
熵值法計算指標權重的本質是比較指標的信息熵效用價值，因此采用信息熵冗余度的占比來代替指標權重，參考式（7），最終采用式（8）合成全部二級指標權重。
式中，wi表示第i項指標在綠色發(fā)展整體系統(tǒng)中的權重；w'i表示第i項指標在所屬一級指標中的權重；ej為對應的一級指標權重。
三、長江經濟帶11?。ㄊ校┚G色發(fā)展指數評價
（一）數據來源
經濟社會發(fā)展相關數據來源于相關年份長江經濟帶11?。ㄊ校┙y(tǒng)計年鑒，其中經濟數據以2013年為不變價進行換算；能源利用數據來源于相關年份《中國能源統(tǒng)計年鑒》；廢水、廢氣、固體廢物排放數據來源于相關年份長江經濟帶11?。ㄊ校┉h(huán)境狀況公報，其他生態(tài)環(huán)境相關數據來源于相關年份《中國環(huán)境統(tǒng)計年鑒》。
（二）計算綠色發(fā)展指數
采用綜合加權法計算綠色發(fā)展指數：
式中，Gj表示j年綠色發(fā)展指數。
綠色發(fā)展指數從多方面綜合反映了區(qū)域綠色發(fā)展水平，同時也是進行綠色發(fā)展?jié)摿α炕u價的基礎。按照上述計算過程得到中國長江經濟帶11?。ㄊ校?013—2018年綠色發(fā)展指數，見表2。
表2 長江經濟帶11?。ㄊ校┚G色發(fā)展指數
（續(xù)表）
本文構建的綠色發(fā)展指標體系包含資源利用、環(huán)境治理、環(huán)境質量、生態(tài)保護、增長質量、綠色生活6個方面。其中資源利用重點關注了水資源、土地資源、能源資源的利用效率以及可回收利用資源的再利用情況；環(huán)境治理是綠色發(fā)展的核心領域之一，反映政府和社會對環(huán)境污染治理的資金投入力度；環(huán)境質量是環(huán)境治理的結果體現；生態(tài)保護體現了自然生態(tài)質量與生態(tài)保護投入情況；增長質量主要表現為經濟總量增長以及經濟發(fā)展的可持續(xù)性；綠色生活主要關注居民生活方式的綠色性。
從長江經濟帶11?。ㄊ校┑木G色發(fā)展指數計算結果，可以得到以下結論：（1）2013年11?。ㄊ校┑木G色發(fā)展指數存在一定差異，但經過近年來綠色發(fā)展理念的貫徹執(zhí)行，到2018年11省（市）的綠色發(fā)展之間的差異明顯縮小，11?。ㄊ校┚G色發(fā)展指數都超過0.90，說明近年來長江經濟帶11?。ㄊ校┳⒅貐f(xié)同發(fā)展。（2）總體看，針對資源利用、環(huán)境治理、環(huán)境質量、生態(tài)保護、增長質量、綠色生活6個方面指標，我們可以看到11省（市）的環(huán)境治理和增長質量兩個一級指標計算結果增長顯著，尤其是上游4個省（市），2013年重慶環(huán)境治理和增長質量兩個指標計算結果均為0.58，2018年兩個指標分別增長至0.84和0.80；四川兩個指標計算結果由2013年的0.61和0.57，增長至2018年的0.89和0.80；云南兩個指標計算結果由2013年的0.64和0.56，增長至2018年的0.84和0.80；貴州兩個指標計算結果由2013年的0.50和0.59，增長至2018年的0.93和0.80。（3）從資源利用角度看，下游與中游7個?。ㄊ校┍壬嫌?個?。ㄊ校┲笜烁?，上游的重慶、四川、貴州和云南資源豐富，但是未充分開發(fā)利用。（4）環(huán)境治理方面，貴州、江西和湖北近年來表現突出。近年來長江兩岸防護林體系建設造林超千萬畝，通過國家濕地公園、綠色發(fā)展示范區(qū)以及退化濕地恢復等建設，長江經濟帶重點湖庫及流域水質優(yōu)良斷面比例顯著提高，水源涵養(yǎng)功能明顯增強。禁漁退漁在長江流域重點水域取得階段性成效，18萬漁民全部提前退捕上岸。禁漁以來，長江水生物資源逐漸恢復，如作為長江流域生態(tài)環(huán)境“指標生物”的江豚在湖北宜昌長江三江段頻現。（5）由于環(huán)境治理成效顯著，長江經濟帶11?。ㄊ校?013—2018年環(huán)境質量大幅提升。（6）從生態(tài)保護角度看，2013—2016年上海生態(tài)保護措施不到位，但是2017年之后出臺相關政策，進而有了大幅提升。（7）從綠色生活角度看，其他省市都有小幅上升，上海2018年反而下降，說明上海需要特別注意居民生活方式的綠色改善。
經濟的綠色發(fā)展水平是衡量綠色發(fā)展的重要測度指標，今后長江經濟帶11?。ㄊ校┙洕l(fā)展的重要內容仍然是促進經濟結構優(yōu)化升級，推動經濟發(fā)展綠色轉型升級，培育新的經濟增長點，拓展綠色發(fā)展空間，不斷提高經濟發(fā)展的質量和效益。
四、綠色發(fā)展系統(tǒng)動力學模型構建
將長江經濟帶區(qū)域發(fā)展系統(tǒng)劃分為經濟、社會、資源、生態(tài)環(huán)境4個子系統(tǒng)，通過分析子系統(tǒng)和變量間的因果關系與數學關系，開發(fā)設計出長江經濟帶11?。ㄊ校?013—2060年綠色發(fā)展仿真模型并進行檢驗，檢驗結果表明模型運行穩(wěn)健，與系統(tǒng)真實運行情況基本相符，能夠用于長江經濟帶11?。ㄊ校┚G色發(fā)展仿真模擬。
（一）系統(tǒng)邊界及基本假設
2020年9月，習近平主席在第七十五屆聯合國大會一般性辯論上提出，中國將采取更加有力的政策和措施，二氧化碳排放力爭于2030年前達到峰值，努力爭取2060年前實現碳中和。因此本文設定模型有效時間為2013—2060年，其中2013—2018年為模型設定與檢驗區(qū)間，2018—2060年為模型預測區(qū)間；設定2018—2030年為中期預測區(qū)間，2030—2060年為長期預測區(qū)間。設定時間間隔為1年。
一般而言，系統(tǒng)邊界內包含與研究問題有關的全部對象，但由于區(qū)域綠色發(fā)展系統(tǒng)強調經濟、社會、資源、生態(tài)環(huán)境4大子系統(tǒng)的整體與協(xié)調發(fā)展，涉及的指標數量眾多，較難逐一列舉。本文從各子系統(tǒng)的內涵出發(fā)，分別選擇與綠色發(fā)展關聯度較高的變量納入模型：（1）經濟子系統(tǒng)中選擇GDP、產業(yè)結構、固定資產投資、財政支出、研發(fā)支出、居民收入等變量；（2）社會子系統(tǒng)中選擇人口、城鎮(zhèn)化率、就業(yè)、醫(yī)療衛(wèi)生水平、公共交通情況等變量；（3）資源子系統(tǒng)中選擇能耗、水耗、土地資源消耗、碳排放等變量；（4）生態(tài)環(huán)境子系統(tǒng)中選擇水、大氣、固體廢物等污染物排放水平及生態(tài)環(huán)境質量、污染治理投入等變量。
（二）系統(tǒng)動力學流圖
因果關系圖只能反映指標間的相關性和反饋過程，需要進一步繪制系統(tǒng)動力學流圖來實現對各要素運行情況的仿真模擬。長江經濟帶11省（市）綠色發(fā)展系統(tǒng)由經濟、社會、資源、生態(tài)環(huán)境4個子系統(tǒng)構成，得到系統(tǒng)動力學流圖①系統(tǒng)動力學流圖可咨詢作者或《統(tǒng)計理論與實踐》編輯部。。
（三）模型方程式
模型方程式定量化描述指標間因果關系。綠色發(fā)展系統(tǒng)動力學模型中涉及的變量類型有水平變量、速率變量、輔助變量以及常量，各變量之間的定量關系通過模型方程予以描述，共建立95個模型方程。
五、長江經濟帶11?。ㄊ校┚G色發(fā)展?jié)摿η榫胺治?/h2>（一）綠色發(fā)展情景設定
情景分析是政策制定的基礎，通過改變不同控制指標的參數，設定多個發(fā)展情景并進行模擬仿真，分析不同參數變化情況對綠色發(fā)展系統(tǒng)的影響，從而為提升區(qū)域綠色發(fā)展?jié)摿μ峁┖侠斫ㄗh。
情景分析一般分為單要素模擬與多要素綜合模擬，單要素模擬主要分析單個控制指標變動對系統(tǒng)的影響，多要素模擬分析多個控制指標同時變動對系統(tǒng)的影響?？紤]到單一指標調整對綠色發(fā)展系統(tǒng)影響較小，從產業(yè)結構優(yōu)化、能源結構優(yōu)化、固定資產投資管控、科技創(chuàng)新資金投入、污染治理資金投入等角度選取多個指標，設定不同的組合方案進行多要素模擬。具體調控指標包括第二產業(yè)增加值占比、第三產業(yè)增加值占比、工業(yè)增加值占比、GDP增速、煤炭消費占比、天然氣消費占比、財政支出占比、環(huán)境污染治理投資占比、工業(yè)固體廢物產污系數、工業(yè)源COD排污系數、工業(yè)固體廢物綜合利用率以及研發(fā)支出增長率，這些指標在系統(tǒng)動力學模型中均以表函數表示。
P1（基礎模式，不改變現有參數）：作為基礎對照組，情景方案設置延續(xù)現狀發(fā)展路徑，各項調控指標參數保持現有變化趨勢，不進行調整。
P2（產業(yè)結構優(yōu)化模式）：第二產業(yè)增加值占比減少10%，第三產業(yè)增加值占比增加10%，工業(yè)增加值占比減少10%，GDP增速變?yōu)?.2倍。
P3（能源結構優(yōu)化模式）：煤炭消費占比減少15%，天然氣消費占比增加5%，財政支出占比變?yōu)?.2倍。
P4（技術進步模式）：研發(fā)支出增長率變?yōu)?.2倍。
P5（環(huán)境治理模式）：環(huán)境污染治理投資占比和工業(yè)固體廢物綜合利用率變?yōu)?.2倍，工業(yè)固體廢物產污系數和工業(yè)源COD排污系數變?yōu)?.8倍。
P6（綜合模式）：同時進行P2、P3、P4、P5的4種調整方案。
由于長江經濟帶11?。ㄊ校┠Ｐ蛿祿跏贾荡嬖诓町?，因而僅列出上海市綠色發(fā)展?jié)摿η榫澳M調控指標的參數設置，如表3所示。其余各?。ㄊ校┡c之相似，但依賴于具體模型數據初始值。
表3 上海市系統(tǒng)動力學模型情景模擬調控指標的參數設置結果
（二）長江經濟帶11省（市）綠色發(fā)展?jié)摿η榫澳M對比分析
將長江經濟帶11?。ㄊ校┚G色發(fā)展?jié)摿?種情景模擬結果進行對比，中期設置為2018—2030年，長期設置為2030—2060年，結果見表4、表5。
據表4、表5的模擬結果，可以看到：（1）長江經濟帶11?。ㄊ校╅L期的綠色發(fā)展指數變化幅度均顯著于中期變化幅度，上海長期綠色發(fā)展指數變化幅度均達到138以上，貴州長期綠色發(fā)展指數變化幅度在第六種情景模式中也達到78.256。（2）不論是中期變化幅度還是長期變化幅度，上海的變化幅度均位列第一，接著是上游4的重慶、四川、貴州、云南，其次是下游的浙江、安徽和中游的湖南，變化幅度最小的是中游的江西、湖北和下游的江蘇。
表4 2018—2030年長江經濟帶11?。ㄊ校┚G色發(fā)展指數變化幅度模擬結果
表5 2030—2060年長江經濟帶11?。ㄊ校┚G色發(fā)展指數變化幅度模擬結果
根據11?。ㄊ校┚G色發(fā)展指數差異，可以得出以下結論：（1）上海作為長三角一體化發(fā)展及長江經濟帶發(fā)展的龍頭，地處流域出?？冢仁且I流域整體轉型發(fā)展的橋頭堡，也是這一區(qū)域與海外市場的重要樞紐。因此，上海未來的綠色發(fā)展?jié)摿薮蟆＃?）長江上游4個?。ㄊ校ㄖ貞c、四川、貴州、云南）位于中國西南部，是長江上游地區(qū)的重要生態(tài)屏障和經濟、金融、科創(chuàng)、航運和商貿物流中心，是長江經濟帶上游成渝城市群的核心城市，其生態(tài)環(huán)境位居全國前列，未來的綠色發(fā)展?jié)摿σ埠艽?。?）近年來長江下游和中游的江蘇、浙江、安徽、湖北、湖南和江西工業(yè)規(guī)?？偭看?、增速快，重化工業(yè)占比偏大，第二產業(yè)發(fā)展對區(qū)域資源環(huán)境造成的負荷仍然較高，煤炭消費量占比較大，清潔能源消費占比不高，綜合看，其未來綠色發(fā)展?jié)摿π枰M一步挖掘。
根據表4、表5列出的6種情景模擬結果差異對比，得出以下結論：
（1）作為基礎對照方案，該情景下綠色發(fā)展指數穩(wěn)定提升，綠色發(fā)展水平存在一定上升空間，綠色發(fā)展?jié)摿^高。但是對于長江經濟帶11?。ㄊ校?，基礎模式下的綠色發(fā)展指數提升幅度均是6種模式下最小的，因而不是最優(yōu)發(fā)展方案。
（2）對比產業(yè)結構優(yōu)化模式、能源結構優(yōu)化模式、技術進步模式和環(huán)境治理模式，產業(yè)結構優(yōu)化模式下長江經濟帶11?。ㄊ校┚G色發(fā)展指數變化幅度最大，其次是技術進步模式，再次是環(huán)境治理模式，最后是能源結構優(yōu)化模式。
（3）產業(yè)結構優(yōu)化模式情景設置是以提高第三產業(yè)和高新技術產業(yè)產值占比，降低工業(yè)增加值占比為主要出發(fā)點，并協(xié)同帶動經濟增長。從結果看，該情景下產業(yè)結構明顯優(yōu)化，經濟增長動能增強，主要是因為資源利用效率快速提升，工業(yè)污染物排放量明顯降低，對綠色發(fā)展?jié)摿μ嵘哂忻黠@促進效果。
（4）技術進步模式下，科技創(chuàng)新能力提升對經濟增長質量、清潔生產水平、環(huán)境污染治理等多方面產生正向促進作用，與前述發(fā)展方案相比較，該情景下各類指標指數均有一定幅度增長。
（5）環(huán)境治理模式下情景設置旨在通過加強環(huán)境污染治理投資、工業(yè)污染管控水平來降低污染物排放，綠色發(fā)展水平大幅度提升。
（6）能源結構優(yōu)化模式情景設定目的在于優(yōu)化能源結構、提升能源利用效率。從結果看，碳排放削減效果明顯，符合情景設定目的，對綠色發(fā)展綜合指數提升具有一定促進效果，但指數增長幅度不大，表明僅通過優(yōu)化能源利用方式已較難快速提升區(qū)域綠色發(fā)展?jié)摿Γ摲桨高m合作為其他方案的補充。
（7）綜合模式情景從資金管控和政策引導兩個方面進行控制，并綜合考慮到產業(yè)結構、能源結構、污染治理、科技創(chuàng)新等方面的因素進行調整。該情景方案下，長江經濟帶11省（市）綠色發(fā)展指數增長幅度最大，采用多項措施的組合方案對提升區(qū)域綠色發(fā)展質量效果最突出，是提升綠色發(fā)展?jié)摿Φ淖顑?yōu)選擇。
總結上述研究結果認為，加強產業(yè)結構轉型和生態(tài)環(huán)境治理投入是提升長江經濟帶11?。ㄊ校┚G色發(fā)展?jié)摿Φ闹攸c所在，應當以節(jié)能降耗減污為政策導向，通過優(yōu)化產業(yè)結構、提高環(huán)境污染治理資金投入，并配合實施能源結構調整、加大技術創(chuàng)新等一系列組合措施，最大程度促進長江經濟帶11省（市）高質量綠色發(fā)展。
六、結論
本文通過構建綠色發(fā)展評價指標體系，得到了反映綠色發(fā)展綜合水平的綠色發(fā)展指數，并設計了包括經濟、社會、資源、生態(tài)環(huán)境4個子系統(tǒng)的長江經濟帶11?。ㄊ校┑木G色發(fā)展系統(tǒng)動力學模型，模型檢驗的可靠性和穩(wěn)定性較好，能夠反映綠色發(fā)展各關鍵指標變化情況。在此基礎上，從優(yōu)化產業(yè)結構、能源清潔利用、促進科技創(chuàng)新、強化污染治理等方面設定不同的組合方案，仿真模擬了長江經濟帶11省（市）2018—2060年綠色發(fā)展變化趨勢。研究結果表明：（1）根據綠色發(fā)展指數的計算結果，近年來長江經濟帶11?。ㄊ校┚G色發(fā)展均取得明顯成效，且各一級指標指數均有提升。（2）長江經濟帶11省（市）未來發(fā)展，仍需要持續(xù)推動產業(yè)結構轉型升級。（3）科技創(chuàng)新為經濟社會發(fā)展提供重要支撐，不僅能夠提升勞動生產率，還能提升資源環(huán)境產出率。長江經濟帶11?。ㄊ校碛休^好的科技創(chuàng)新土壤，人才科教資源十分豐富，應大力增加綠色技術研發(fā)投入，促進綠色產業(yè)發(fā)展，進一步凸顯綠色技術在優(yōu)化產業(yè)結構、提升資源利用效率、促進污染減排等方面的作用。（4）未來長江經濟帶11?。ㄊ校^續(xù)加大環(huán)境污染治理資金投入力度，削減污染物排放總量。（5）長江經濟帶11省（市）可以根據自身實際情況提升未來綠色發(fā)展?jié)摿??！?/p>
猜你喜歡
法人(2021年2期)2021-03-02 07:12:46
青年歌聲(2020年11期)2020-11-24 06:57:28
當代陜西(2018年9期)2018-11-18 07:26:13
天津詩人(2017年2期)2017-11-29 01:24:30
天津詩人(2017年2期)2017-11-29 01:24:05
專用汽車(2016年5期)2016-03-01 04:14:39
中國資源綜合利用(2016年7期)2016-02-03 03:00:14
學習月刊(2015年9期)2015-07-09 05:33:46
