王立成 張 軍 陳前前

(山東工商學院經(jīng)濟學院，山東 煙臺 264005；北京物資學院商學院，北京 101149)

商貿(mào)流通業(yè)上市公司發(fā)展指數(shù)綜合評價研究

王立成 張 軍 陳前前

(山東工商學院經(jīng)濟學院，山東 煙臺 264005；北京物資學院商學院，北京 101149)

本文利用主成分分析法，選擇經(jīng)營規(guī)模等7大類內(nèi)含22個單項指標構(gòu)建了我國商貿(mào)流通上市公司發(fā)展狀況的綜合指標體系——商貿(mào)流通業(yè)上市公司發(fā)展指數(shù)，并根據(jù)該指數(shù)綜合評價我國商貿(mào)流通業(yè)上市公司的發(fā)展狀況。研究發(fā)現(xiàn)，2008—2015年期間我國商貿(mào)流通業(yè)上市公司發(fā)展水平整體上表現(xiàn)為“增長—穩(wěn)定—增長”特征，呈現(xiàn)穩(wěn)步上升趨勢；公司運營能力、利潤成長是支撐其發(fā)展的正向因素，而盈利能力是制約行業(yè)發(fā)展的短板因素。

商貿(mào)流通業(yè)；發(fā)展指數(shù)；主成分模型；因子分析

一、引言

商貿(mào)流通業(yè)是連接社會生產(chǎn)和消費的中間環(huán)節(jié)，在國民經(jīng)濟發(fā)展中起到平衡生產(chǎn)和消費的作用，一方面實現(xiàn)產(chǎn)品的最終消費，完成資金和商品循環(huán)流動的最后一環(huán)；另一方面向生產(chǎn)環(huán)節(jié)反饋消費信息，引導(dǎo)生產(chǎn)向著滿足需求的方向發(fā)展。商貿(mào)流通業(yè)的發(fā)展水平是市場經(jīng)濟成熟程度的直接體現(xiàn)，對商貿(mào)流通業(yè)發(fā)展水平及其影響因素的研究不僅有助于提升該行業(yè)的發(fā)展水平，對國民經(jīng)濟的發(fā)展也具有重要的現(xiàn)實意義。

對于我國商貿(mào)流通業(yè)的研究，隨著商貿(mào)流通業(yè)態(tài)的日益豐富而不斷增多。劉曉雪①劉曉雪：《商業(yè)企業(yè)規(guī)模結(jié)構(gòu)啞鈴狀分布究因》，《北京工商大學學報(社會科學版)》2007年第9期。總結(jié)了北京地區(qū)批發(fā)零售行業(yè)企業(yè)發(fā)展的特點，袁武聰?shù)娶谠渎?、夏春玉、曹志來：《中國大型零售企業(yè)規(guī)模擴張方式的選擇》，《中國零售研究》2009年第6期。研究了我國大型零售企業(yè)的擴張方式選擇問題；陳金偉等③陳金偉、張昊：《零售企業(yè)規(guī)模不經(jīng)濟問題研究》，《中國流通經(jīng)濟》2013年第3期。分析了我國零售企業(yè)擴張過程中規(guī)模不經(jīng)濟的主要原因；趙萍④趙萍：《2014年中國流通產(chǎn)業(yè)回顧與2015年展望》，《中國流通經(jīng)濟》2015年第1期。結(jié)了2014年中國流通產(chǎn)業(yè)產(chǎn)生的新特點、闡述了我國零售業(yè)2015年的發(fā)展情況，并對2016年的發(fā)展趨勢進行了展望。洪濤等⑤洪濤、鄭強：《城市流通力的內(nèi)涵及其相應(yīng)指標體系的建立》，《商業(yè)經(jīng)濟與管理》2002年第11期。、孫薇⑥孫薇：《基于因子分析法的地區(qū)流通力比較研究》，《財貿(mào)研究》2005年第4期。、宋則等⑦⑧宋則、張弘：《中國流通現(xiàn)代化評價指標體系研究》，《商業(yè)時代》2003年第11期。宋則、郭冬樂、荊林波：《中國流通理論前沿》，社會科學文獻出版社2006版，第60-83頁。、石憶邵等⑨石憶邵、朱衛(wèi)鋒：《商貿(mào)流通業(yè)競爭力初探：以南通市為例》，《財經(jīng)研究》2004年第5期。，周日星等⑩周日星、蘇為華、張悅等：《商貿(mào)流通業(yè)統(tǒng)計監(jiān)測評價體系研究》，中國市場出版社2006版，第20-33頁。分別從流通力、現(xiàn)代化程度、行業(yè)競爭力等角度對商貿(mào)流通業(yè)進行了評價，通過構(gòu)建指標體系等方法對商貿(mào)流通業(yè)發(fā)展水平進行考察。這些評價方法或側(cè)重于商貿(mào)流通業(yè)發(fā)展的單一方面，或側(cè)重于評價本身，未能從整體上考察商貿(mào)流通業(yè)的發(fā)展情況，進而無法為提高行業(yè)發(fā)展水平提供可靠的依據(jù)。

本文利用主成分分析法，使用商貿(mào)流通行業(yè)上市公司的財務(wù)數(shù)據(jù)構(gòu)建發(fā)展指數(shù)來對商貿(mào)流通業(yè)上市公司的發(fā)展水平進行評價和分析，首先，從指數(shù)層面對業(yè)內(nèi)上市公司的發(fā)展水平進行了分析；其次，對各主成分的值也進行了統(tǒng)計，并與行業(yè)的發(fā)展情況相互印證。本文的創(chuàng)新之處在于，在對原始指標進行重歸類的基礎(chǔ)上，分析了各主成分在一段時間內(nèi)的變化趨勢，展示了各主成分對行業(yè)發(fā)展水平影響的時間特征及其相互關(guān)系。二是將商貿(mào)流通業(yè)細分為若干子行業(yè)，利用發(fā)展指數(shù)，對不同子行業(yè)的發(fā)展情況進行分析。

二、研究設(shè)計

(一)樣本選取

經(jīng)濟體中的所有微觀企業(yè)構(gòu)成了其宏觀經(jīng)濟主體中的企業(yè)，但非上市公司不需要公開披露財務(wù)報表，難以獲得其財務(wù)數(shù)據(jù)，為了研究數(shù)據(jù)可獲得，研究對象選擇了境內(nèi)A股上市公司。首先，本文的樣本選取時間區(qū)間為2008—2015年。這主要是由于我國于2007年實施了新的企業(yè)會計準則，難免會存在會計準則的變遷所帶來的若干不一致現(xiàn)象。為此，基于上述考慮，本研究中的樣本從2008年始，到2015年止。其次，本文選取的樣本公司有如下限制條件：一是上市時間必須在2008年(不含)之前，即挑選符合樣本時間區(qū)間的樣本公司；二是剔除掉上市公司中的ST公司，主要是為了規(guī)避數(shù)據(jù)的異常情況；三是樣本公司必須在2008—2015年區(qū)間內(nèi)未發(fā)生主營業(yè)務(wù)轉(zhuǎn)變或退市，以避免數(shù)據(jù)的不確定性；四是樣本公司數(shù)據(jù)完整，不存在數(shù)據(jù)缺失，以保證數(shù)據(jù)的一致性。本研究從國泰安數(shù)據(jù)庫中獲得了餐飲業(yè)、倉儲業(yè)、道路運輸業(yè)、航空運輸業(yè)、零售業(yè)、批發(fā)業(yè)、水上運輸業(yè)、鐵路運輸業(yè)、住宿業(yè)上市公司的主要財務(wù)指標，得到了172家商貿(mào)流通業(yè)上市公司，共1376個樣本。

(二)評價指標選取

本文主要著眼于企業(yè)內(nèi)部因素對商貿(mào)流通業(yè)上市公司發(fā)展水平的影響，以上市公司公開披露的財務(wù)數(shù)據(jù)作為基礎(chǔ)指標的主要來源。在遵循全面性、科學性、系統(tǒng)性和可操作性等原則的基礎(chǔ)上，充分借鑒現(xiàn)有研究成果，選擇了7個大類，22個小項的財務(wù)指標。具體見表1：

表1 商貿(mào)流通上市公司發(fā)展狀況評價指標體系

(三)研究方法

對各項選取指標進行相應(yīng)的賦權(quán)是構(gòu)建發(fā)展指數(shù)的核心問題。為分析172家商貿(mào)流通業(yè)上市公司在2008-2015年中紛繁復(fù)雜的數(shù)據(jù)，本文試圖采用目前國際社會研究中比較流行的主成分分析法，分析出上述22項指標各自的權(quán)重，最終得到商貿(mào)流通業(yè)上市公司的發(fā)展指數(shù)。

三、發(fā)展指數(shù)編制

(一)適用性檢驗

在使用主成分分析法確定主要影響因子之前，需要首先使用Stata軟件對原始數(shù)據(jù)指標進行KMO(Kaiser-Meyer-Olkin)檢驗和巴特利球形(Bartlett's Testof Sphericity)檢驗，以此來辨別所選取樣本數(shù)據(jù)是否可以進行有效的因子分析過程。針對商貿(mào)流通業(yè)上市公司數(shù)據(jù)的KMO檢驗和巴特利球形檢驗結(jié)果顯示：KMO 統(tǒng)計量的具體數(shù)值為0.693，大于0.5，則KMO檢驗有效性成立；巴特利特球形檢驗的數(shù)值為2191.918，顯著性P值為0.000，小于0.05，則巴特利特球形檢驗有效性成立。因此本文所選取的樣本數(shù)據(jù)可以滿足主成分分析的要求。

表2 商貿(mào)流通業(yè)上市公司數(shù)據(jù)的主成分特征值和方差貢獻率值

(二)數(shù)據(jù)的縮尾處理和標準化處理

為了使得數(shù)據(jù)分析更加符合實際情況，需要對數(shù)據(jù)進行微小的調(diào)整和處理?？s尾處理是將數(shù)據(jù)兩側(cè)各進行了2.5%的縮尾，從而避免了極端數(shù)值的不準確性對主成分的影響；標準化處理是將所有數(shù)據(jù)進行均值為0和方差為1的標準簡化，從而有效防止由于各個指標的量綱不同對分析結(jié)果帶來不良影響。

(三)主成分因素的選擇

表2列示了各成分的特征值和對方差的解釋程度，前7個成分的特征值大于1，累計貢獻率也達到了72.8%，因此，提取這7個成分來代替原始的財務(wù)指標，作為主成分來衡量商貿(mào)流通業(yè)上市公司的發(fā)展水平。

(四)各主成分的經(jīng)濟意義

表3報告了旋轉(zhuǎn)后的因子載荷矩陣，反映各原始指標對主成分的載荷大小，即對各主成分的解釋程度。

表3 旋轉(zhuǎn)后的因子載荷矩陣

根據(jù)7個主成分中相關(guān)指標顯示的主要經(jīng)濟內(nèi)容，將其分別進行了因子的重新命名。速動資產(chǎn)比率、流動資產(chǎn)比率、現(xiàn)金流動負債比和資產(chǎn)負債比率等各項指標，主要反映了商貿(mào)流通業(yè)上市公司的財務(wù)穩(wěn)健性，故將主成分f1定義為穩(wěn)健因子；凈資產(chǎn)收益率、每股收益和資產(chǎn)凈利率等各項指標，主要反映了公司的盈利能力，故將主成分f2定義為盈利因子；總資產(chǎn)周轉(zhuǎn)率、應(yīng)付賬款周轉(zhuǎn)率、期間費用率和銷售毛利率等各項指標，主要反映了公司的銷售和運營管理能力，故將主成分f3定義為運營因子；資產(chǎn)總額、營業(yè)收入和凈利潤等各項指標，主要反映了公司的規(guī)模因素，故將主成分f4定義為規(guī)模因子；資產(chǎn)總額增長率和營業(yè)收入增長率等各項指標，主要反映了公司的規(guī)模增長情況，故將主成分f5定義為規(guī)模成長因子；凈利潤增長率和營業(yè)利潤增長率等各項指標，主要反映了公司的利潤增長情況，故將主成分f6定義為利潤成長因子；應(yīng)收賬款周轉(zhuǎn)率、營業(yè)收入現(xiàn)金含量和非經(jīng)常損益比率等各項指標，主要反映了公司的現(xiàn)金管理水平，故將主成分f7定義為財務(wù)質(zhì)量因子。

(五)主成分得分矩陣

表4中的數(shù)字為主成分得分的系數(shù)矩陣，即為計算各主成分的系數(shù)。通過標準化后的原始數(shù)據(jù)矩陣乘以主成分得分系數(shù)矩陣，從而可以確定7個主成分的得分函數(shù)。

表4 主成分得分的系數(shù)矩陣

(六)建立主成分分析模型

依據(jù)上述7個主成分的得分，可以具體計算出各主成分對商貿(mào)流通業(yè)上市公司發(fā)展指數(shù)的貢獻率。并將此7個貢獻率作為主成分因子的系數(shù)，生成商貿(mào)流通業(yè)上市公司的主成分模型如下：

原始數(shù)據(jù)代入模型的公式中，從而得到172家上市公司在各年份的主成分f值結(jié)果。

(七)發(fā)展指數(shù)的修正

經(jīng)過上述主成分f值的計算過程，需要采用對數(shù)變換的方法使負數(shù)的f值變?yōu)檎龜?shù)，從而便于對各個年份和各個上市公司之間進行比較。具體對發(fā)展指數(shù)的修正公式如下：

DEV=ef×100

其中，DEV為上市公司的發(fā)展指數(shù)，e為自然對數(shù)，f為主成分的f值。

上述公式中乘以100的作用是，可以將原有的較小取值的f值擴張，方便對發(fā)展指數(shù)的比較。

就企業(yè)發(fā)展來看，DEV值越大越好，越大的發(fā)展指數(shù)表明在整個行業(yè)中該企業(yè)的發(fā)展水平越高。由于樣本數(shù)據(jù)的選取區(qū)間為2008-2015年，因此本文在計算出各個年度以及各個上市公司的DEV值之后，采用當年度樣本公司DEV值的平均值作為行業(yè)發(fā)展水平的衡量指標。

四、發(fā)展指數(shù)的綜合分析

(一)發(fā)展指數(shù)的總體變化趨勢分析

通過計算2008—2015年商貿(mào)服務(wù)業(yè)上市公司的DEV值，可以得到9個行業(yè)的行業(yè)發(fā)展指數(shù)，并與商貿(mào)流通業(yè)總指數(shù)同列于表5。

表5 2008—2015年商貿(mào)流通業(yè)上市公司發(fā)展指數(shù)細分情況表

從表5中可以看到，我國商貿(mào)流通業(yè)上市公司的發(fā)展指數(shù)總體上呈現(xiàn)出增長趨勢，可分為三個階段：第一階段是2008—2010年，總指數(shù)有大幅度的快速增長。各細分行業(yè)中，除住宿業(yè)、餐飲業(yè)和水上運輸業(yè)下降外，其余細分行業(yè)的發(fā)展指數(shù)都有增長，尤其是航空運輸業(yè)、鐵路運輸業(yè)和倉儲業(yè)增長幅度最大，遠遠高于其他細分行業(yè)。第二階段是2010—2014年，總指數(shù)的基本保持平穩(wěn)，變動幅度微小。各細分行業(yè)中，倉儲業(yè)發(fā)展指數(shù)表現(xiàn)出明顯的下降趨勢；水上運輸業(yè)穩(wěn)步增長，其余細分行業(yè)發(fā)展指數(shù)則基本保持穩(wěn)定。第三階段是2014—2015年，總指數(shù)又出現(xiàn)小幅度的增長，達到了最大值143。多數(shù)細分行業(yè)的發(fā)展指數(shù)都變現(xiàn)為上升，只有倉儲業(yè)、鐵路運輸業(yè)是下降的，其中鐵路運輸業(yè)的下降趨勢較為明顯，下降幅度也很大。多數(shù)細分行業(yè)發(fā)展指數(shù)的變化趨勢與總指數(shù)的趨勢一致，表現(xiàn)為“增長—穩(wěn)定—增長”的發(fā)展趨勢。

其中最重要的原因是，2008年以后，我國為應(yīng)對美國次貸危機帶來的挑戰(zhàn)，由政府主導(dǎo)進行了大規(guī)模的投資活動，使我國經(jīng)濟增速在短暫急劇下降后，表現(xiàn)出了快速回升的趨勢，商貿(mào)流通指數(shù)的變化與投資的效用恰好相符。從指數(shù)的變化上看，2009—2010年，總指數(shù)大幅度的躍升主要歸因于航空運輸業(yè)、鐵路運輸業(yè)和倉儲業(yè)指數(shù)的快速增長，這三個行業(yè)正是4萬億資金的主要投向。四萬億資金及其帶動的地方政府和民營資金所進行的大規(guī)?；A(chǔ)設(shè)施建設(shè)明顯促進了上述行業(yè)的發(fā)展，使其發(fā)展水平有了大幅度的提升。

2010年以后，投資增速下降，各細分行業(yè)的發(fā)展也趨于穩(wěn)定，其中對投資依賴程度較高的航空、鐵路和倉儲等行業(yè)指數(shù)還出現(xiàn)了下滑的趨勢，但幅度較小。2014年以后，受航空運輸業(yè)和批發(fā)業(yè)的影響，商貿(mào)流通業(yè)上市公司的發(fā)展指數(shù)又有小幅提升。

發(fā)展指數(shù)較高的細分行業(yè)是航空運輸業(yè)和鐵路運輸業(yè)，受四萬億投資的影響，其增長速度也遠遠高于其余細分行業(yè)；較低的是餐飲業(yè)和住宿業(yè)，發(fā)展指數(shù)基本保持在60左右；其余細分行業(yè)則處于中間水平。這一現(xiàn)象也與我國流通領(lǐng)域的發(fā)展現(xiàn)狀基本吻合，以運輸業(yè)為例，航空運輸和鐵路運輸?shù)陌l(fā)展要明顯好于水上運輸，餐飲和住宿業(yè)也沒有形成行業(yè)發(fā)展的龍頭企業(yè)，發(fā)展水平普遍較低。以快餐為例，雖然存在眾多的地方性連鎖企業(yè)，但其經(jīng)營規(guī)模和發(fā)展速度都遠遠不及行業(yè)內(nèi)的國際巨頭。2012年，中央八項規(guī)定出臺以后，餐飲和住宿業(yè)的發(fā)展形勢持續(xù)低迷，同時面臨著市場萎縮、業(yè)務(wù)增長乏力、營業(yè)成本上升等難題，使其發(fā)展指數(shù)持續(xù)下降。

表6 2008—2015年各主成分平均得分情況

商貿(mào)流通業(yè)上市公司2008—2015年各主成分的得分見表6。數(shù)據(jù)顯示：2009—2010年是發(fā)展指數(shù)增長最快的時期，主要原因是受到了f1、f2、f5和f7四個因子的影響，2010年以后，商貿(mào)流通業(yè)上市公司發(fā)展水平保持穩(wěn)定，發(fā)展指數(shù)保持在132左右，一方面受益于規(guī)模增長和運營能力提升的因素，另一方面也受到市場不利因素的影響。2012年以后，各主成分的變化幅度明顯面小，各主成分的變化趨勢也有趨同的特點。綜合以上因素，我國商貿(mào)流通業(yè)上市公司的發(fā)展指數(shù)在2008—2015年期間，呈現(xiàn)出“增長—穩(wěn)定—增長”的趨勢，受宏觀經(jīng)濟影響的趨勢較為明顯。

(二)各主成分因子分析

各主成分對成長指數(shù)的影響分為兩大類，第一類包括f1、f2、f5和f7，第二類包括f3、f4和f6。第一類先上升后下降，第二類先下降后上升，其最終的作用結(jié)果是總指數(shù)“增長—穩(wěn)定—增長”的變化趨勢。

圖1 各主成分變化趨勢圖

f1、f4和f5受投資影響較為明顯，2010年前后，隨著大規(guī)模資金的投入，上市公司的財務(wù)穩(wěn)健性普遍增強，資產(chǎn)和收入等指標也大幅度提升，規(guī)模增長較快，在七個主成分中，f4(規(guī)模因子)的增長幅度也是最大的。但隨著資產(chǎn)的不斷形成，流動資金占比降低，上市公司的財務(wù)穩(wěn)健性也迅速下降，直至2012年，才出現(xiàn)了小幅改善。f5(規(guī)模成長因子)在2009—2013年之間也處于較高的水平上，2013年以后規(guī)模增長的速度明顯下降。f2(盈利因子)的波動幅度最大，在2010年達到最大值后，持續(xù)下降，到2015年已經(jīng)成為七個主成分中最低的一個，f2的下降反映了我國商貿(mào)流通業(yè)上市公司盈利能力的下滑，這也是商貿(mào)流通業(yè)上市公司發(fā)展中面臨的最大問題。f3和f6持續(xù)上升，是行業(yè)指數(shù)保持穩(wěn)定和實現(xiàn)增長的主要正向因素，對指數(shù)增長的貢獻僅次于規(guī)模因子。f3在2010—2012年有大幅的增長，反應(yīng)了上市公司運營能力的快速提升。雖然在2012年以后又有小幅下降，但一直保持在較高水平，且2015年有所提高。f7呈現(xiàn)出波浪式的變化趨勢，一方面上市公司通過盈余管理操控利潤的行為和動機不斷減少，企業(yè)財務(wù)質(zhì)量提升；另一方面，市場中的不利因素增多，使上市公司應(yīng)收賬款周轉(zhuǎn)率和營業(yè)收入現(xiàn)金含量下降，降低了上市公司抵御財務(wù)風險和市場環(huán)境周期性變化的能力。

五、結(jié)論與討論

(一)研究結(jié)論

第一，商貿(mào)流通業(yè)上市公司發(fā)展趨勢表現(xiàn)為“增長-穩(wěn)定-增長”。一方面行業(yè)的發(fā)展受到宏觀經(jīng)濟發(fā)展水平和經(jīng)濟政策的影響與制約，另一方面他也是宏觀經(jīng)濟健康發(fā)展的催化劑，可以起到促進宏觀經(jīng)濟發(fā)展的作用。商貿(mào)流通業(yè)發(fā)展既是國民經(jīng)濟的重要組成部分，也會隨著產(chǎn)業(yè)結(jié)構(gòu)和消費結(jié)構(gòu)的升級換代而帶來自身的發(fā)展機會。

第二，“四萬億投資”顯著刺激了我國航空運輸和鐵路運輸行業(yè)的發(fā)展水平，2009—2010年間這兩個行業(yè)的發(fā)展指數(shù)出現(xiàn)了躍升，行業(yè)規(guī)模增長較快，財務(wù)穩(wěn)健性普遍增強，資產(chǎn)和收入等指標得到改善。但是刺激政策也在一定程度上加劇了商貿(mào)流通業(yè)內(nèi)部發(fā)展水平的分化。

第三，運營因子和利潤成長因子是促進我國商貿(mào)流通業(yè)上市公司最重要的兩個因素，上市公司應(yīng)維持和進一步改善運營能力，實現(xiàn)運營因子的穩(wěn)定增長，并保持利潤的持續(xù)增長。這需要商貿(mào)流通行業(yè)上市公司積極應(yīng)對消費者消費習慣改變的趨勢，積極籌劃，引領(lǐng)消費趨勢，主動對接“互聯(lián)網(wǎng)+”。

最后，盈利能力因子是制約我國商貿(mào)流通業(yè)上市公司發(fā)展水平的短板因素，用凈資產(chǎn)收益率、每股收益等相對指標來反映的盈利能力因子在2010年后呈現(xiàn)持續(xù)下降趨勢，商貿(mào)流通上市公司未來必須著力補齊該因素。

(二)進一步探討

使用財務(wù)指標構(gòu)建的發(fā)展指數(shù)雖然可以反映公司經(jīng)營的結(jié)果，但并不能完全反映產(chǎn)生這種結(jié)果的動因。財務(wù)指標可以反映企業(yè)的規(guī)模、增長的潛力、財務(wù)風險的高低，但并不全面，還有行業(yè)發(fā)展前景、企業(yè)文化、企業(yè)戰(zhàn)略選擇等多方面非財務(wù)的因素。但限于數(shù)據(jù)的可獲得性，本文并未將上述指標納入評價體系之中，此問題的研究仍有待進一步完善。

(責任編輯：曹守新)

F32

A

1003-4145[2017]10-0146-06

2017-07-10

王立成(1966—)，男，遼寧朝陽人，山東工商學院經(jīng)濟學院教授，研究方向：產(chǎn)業(yè)經(jīng)濟與公司治理。張 軍(1977—)，男，山東臨沂人，北京物資學院商學院副教授，管理學博士，主要研究方向：會計理論與企業(yè)發(fā)展。陳前前(1988—)，女，山東濟南人，北京物資學院商學院講師，會計學博士，主要研究方向：會計理論與企業(yè)發(fā)展。

本文系教育部人文社會科學研究基金項目(項目編號：15YJC790144)、北京社科基金研究基地項目(項目編號：16JDGLA016)、北京物資學院科研攻關(guān)團隊階段性研究成果。