李思維 張正勇 孫樹壘

[摘要]2008年中國奶制品污染事件已過去十多年，乳制品業(yè)經(jīng)過這十多年的重新洗牌，呈現(xiàn)出另一番景象。根據(jù)相關(guān)新聞報(bào)道，某些曾經(jīng)受到波及的企業(yè)經(jīng)過這些年的精益求精，目前已成為行業(yè)中的領(lǐng)跑者，然而某些未受到影響反而在當(dāng)時(shí)獲利的企業(yè)卻由于近年來管理決策不當(dāng)，地位早已不復(fù)從前。本文為了解自中國奶制品污染事件以來不同企業(yè)在乳制品供應(yīng)鏈各環(huán)節(jié)的效率差異，構(gòu)建由生產(chǎn)加工環(huán)節(jié)與儲(chǔ)運(yùn)銷售環(huán)節(jié)組成的兩階段數(shù)據(jù)包絡(luò)分析（Data Envelopment Analysis， DEA）的效率測度模型，并運(yùn)用軟件將相關(guān)企業(yè)各環(huán)節(jié)投入產(chǎn)出數(shù)據(jù)代入模型進(jìn)行分析討論。結(jié)果表明，不同企業(yè)在供應(yīng)鏈的不同環(huán)節(jié)各有優(yōu)勢，在各環(huán)節(jié)處于劣勢的企業(yè)應(yīng)向優(yōu)勢企業(yè)學(xué)習(xí)借鑒，并通過減小資源投入或提高技術(shù)能力，不斷改進(jìn)完善自我;乳制品供應(yīng)鏈各環(huán)節(jié)對其整體效率的影響也各不相同，在實(shí)際管理過程中，管理者應(yīng)根據(jù)影響程度的大小，合理設(shè)置權(quán)重以便更好地進(jìn)行分析。

[關(guān)鍵詞]數(shù)據(jù)包絡(luò)分析;乳制品企業(yè);供應(yīng)鏈;技術(shù)效率;績效研究

中圖分類號(hào)：C931.2文獻(xiàn)標(biāo)識(shí)碼：ADOI：10.16465/j.gste.cn431252ts.202003

2008年9月，三鹿奶粉中含有三聚氰胺的新聞報(bào)道引發(fā)社會(huì)各界的強(qiáng)烈關(guān)注，2008年中國奶制品污染事件使得眾多乳制品紛紛下架，乳制品企業(yè)成為眾矢之的，整個(gè)乳制品行業(yè)面臨前所未有的信任危機(jī)。

如今，國內(nèi)的乳制品工業(yè)正由工業(yè)化、機(jī)械化向智能化階段邁進(jìn)。相關(guān)數(shù)據(jù)表明，2018年我國乳制品產(chǎn)量達(dá)2 687.1萬t，主營業(yè)務(wù)收入高達(dá)3 398.9億元，占食品制造行業(yè)收入的18.5%。種種數(shù)據(jù)表明乳制品業(yè)目前已成為中國食品制造業(yè)的第一大行業(yè)。2008—2018年中國乳業(yè)的變化翻天覆地，政府建立了更為嚴(yán)格的市場監(jiān)管體系，乳制品企業(yè)也紛紛從原料采集、生產(chǎn)加工、儲(chǔ)運(yùn)銷售等供應(yīng)鏈的各環(huán)節(jié)著手，對乳制品質(zhì)量進(jìn)行監(jiān)管。

為了解長時(shí)間以來乳制品業(yè)供應(yīng)鏈各環(huán)節(jié)的績效變化，以及企業(yè)間相同供應(yīng)鏈各環(huán)節(jié)的效率差異，本文利用數(shù)據(jù)包絡(luò)分析（Data Envelopment Analysis， DEA）分別對三家企業(yè)連續(xù)多年的供應(yīng)鏈各環(huán)節(jié)進(jìn)行細(xì)致的分析，構(gòu)建兩階段乳制品供應(yīng)鏈績效測量指標(biāo)體系，并建立乳制品企業(yè)兩階段供應(yīng)鏈績效測度模型，將三家企業(yè)各環(huán)節(jié)投入產(chǎn)出數(shù)據(jù)代入其中，進(jìn)行測量與結(jié)果分析。

國內(nèi)關(guān)于乳制品供應(yīng)鏈的研究大多集中在產(chǎn)品的質(zhì)量安全方面。王磊等[1]基于質(zhì)量安全的視角探討了合作伙伴特性對乳制品供應(yīng)鏈合作關(guān)系穩(wěn)定性的影響;楊琦等[2]應(yīng)用主成分因子分析和二項(xiàng)Logistic回歸分析，對我國乳制品供應(yīng)鏈中銷售環(huán)節(jié)中的質(zhì)量安全影響因素進(jìn)行研究;桂冬冬等[3]提出一種包含原料供應(yīng)商、生產(chǎn)商和消費(fèi)市場的供應(yīng)鏈網(wǎng)絡(luò)均衡模型，研究供應(yīng)鏈各主體的決策行為對食品質(zhì)量安全的影響程度;另有學(xué)者基于供應(yīng)鏈的視角，應(yīng)用博弈論等方法對食品質(zhì)量安全風(fēng)險(xiǎn)進(jìn)行相關(guān)的管控研究[4-6]。而國外關(guān)于食品供應(yīng)鏈的研究，主要集中在供應(yīng)鏈各環(huán)節(jié)的安全管理[7-8]以及整體的可持續(xù)性方面[9-10]，與國內(nèi)相比更加側(cè)重于供應(yīng)鏈的管理以及各環(huán)節(jié)效率的提高。

在乳制品供應(yīng)鏈效率測度方面，國內(nèi)的研究大部分為宏觀性研究，例如魏艷驕等[11]利用中國28個(gè)原奶主產(chǎn)省區(qū)面板數(shù)據(jù)，分析乳制品進(jìn)口對于中國不同規(guī)模奶牛養(yǎng)殖效率的影響;陳吉銘等[12]基于京津冀地區(qū)有代表性的企業(yè)數(shù)據(jù)，對其冷鏈物流效率進(jìn)行研究。國內(nèi)關(guān)于效率測度的研究大多是利用DEA對供應(yīng)鏈整體[13-14]或是分階段進(jìn)行研究[15- 16]，然而在乳制品業(yè)的效率測度方面，則缺少對其供應(yīng)鏈“黑箱”內(nèi)部的研究。在國外，近年來關(guān)于DEA的研究，大都針對DEA模型本身進(jìn)行改進(jìn)[17-18]，或?qū)EA與其他研究方法相結(jié)合進(jìn)行實(shí)例分析[19]，但是關(guān)于乳制品供應(yīng)鏈“黑箱”的研究同樣十分有限。因此，本文基于國內(nèi)外在此方面的共同不足，“打開黑箱”，構(gòu)建乳制品供應(yīng)鏈兩階段DEA模型，并從微觀角度對三家乳制品企業(yè)連續(xù)多年的績效進(jìn)行對比分析，以期為各企業(yè)提供相關(guān)參考，促進(jìn)企業(yè)在激烈的市場競爭中更快更好地發(fā)展，并逐漸獲得更高的關(guān)注度與知名度。

1 模型與數(shù)據(jù)

1.1 模型構(gòu)建

乳制品供應(yīng)鏈由于其產(chǎn)品的特殊性，對各環(huán)節(jié)的把控十分嚴(yán)格，在對供應(yīng)鏈進(jìn)行績效評估時(shí)，由于原料采集階段受牧場規(guī)模與乳牛產(chǎn)奶量不可控因素的影響，針對該階段的績效改進(jìn)短時(shí)間內(nèi)難以實(shí)現(xiàn)。因此，本文只考慮乳制品供應(yīng)鏈的生產(chǎn)加工和儲(chǔ)運(yùn)銷售兩個(gè)階段，并利用人機(jī)料法環(huán)的原則建立相應(yīng)的投入產(chǎn)出指標(biāo)體系。

在生產(chǎn)加工階段，員工的精細(xì)操作和設(shè)備的有效運(yùn)行是保證產(chǎn)品質(zhì)量安全的充分必要條件，員工投入數(shù)量的多少與固定資產(chǎn)投入的大小是影響乳制品最終產(chǎn)量的因素之一。因此，將生產(chǎn)技術(shù)人員與固定設(shè)備資產(chǎn)作為該階段的投入指標(biāo)，將乳制品產(chǎn)量作為產(chǎn)出指標(biāo)。

在儲(chǔ)運(yùn)銷售階段，除了生產(chǎn)加工階段的乳制品產(chǎn)量投入，還存在大量的為產(chǎn)品銷售輔助產(chǎn)生的費(fèi)用，比如廣告營銷費(fèi)、裝卸運(yùn)輸費(fèi)、職工薪酬等?？紤]到指標(biāo)的精簡問題，將這些費(fèi)用統(tǒng)稱為銷售費(fèi)用，儲(chǔ)運(yùn)銷售階段的投入指標(biāo)為乳制品產(chǎn)量（生產(chǎn)加工階段產(chǎn)出指標(biāo)）、銷售費(fèi)用與銷售人員，將利潤總額作為最終的產(chǎn)出指標(biāo)。

綜上所述，本文所建立的乳制品供應(yīng)鏈兩階段指標(biāo)體系如圖1所示。

由圖1可知，乳制品產(chǎn)量既是上一階段的產(chǎn)出指標(biāo)，又是下一階段的投入指標(biāo)，本文將這類指標(biāo)統(tǒng)稱為中間產(chǎn)品。值得注意的是，由于相關(guān)企業(yè)年報(bào)中乳制品產(chǎn)量這一數(shù)據(jù)公布不完全，因此采用乳制品營業(yè)成本這一指標(biāo)來反映乳制品產(chǎn)量的多少，乳制品產(chǎn)量越大營業(yè)成本越高，從而這一指標(biāo)在構(gòu)建的DEA模型中可以很好地代替乳制品產(chǎn)量指標(biāo)。

在建立指標(biāo)體系的基礎(chǔ)上，進(jìn)一步構(gòu)建兩階段乳制品供應(yīng)鏈結(jié)構(gòu)模型（如圖2所示），圖2中每個(gè)子階段都對應(yīng)著各自的投入產(chǎn)出指標(biāo)。

顯而易見，上述模型結(jié)構(gòu)既包含并聯(lián)又包含串聯(lián)，屬于復(fù)合型結(jié)構(gòu)。以往有相關(guān)學(xué)者針對不同類型結(jié)構(gòu)的生產(chǎn)系統(tǒng)提出了相應(yīng)的DEA模型，本文在此基礎(chǔ)上根據(jù)圖2分別建立乳制品供應(yīng)鏈兩階段的DEA模型，采用規(guī)模報(bào)酬可變的BC2模式，并作相應(yīng)假設(shè)如下。

生產(chǎn)加工階段和儲(chǔ)運(yùn)銷售階段的第i個(gè)決策單元投入指標(biāo)權(quán)重值分別為αi，βi，乳制品供應(yīng)鏈兩階段的效率分別為θ1，θ2。

1.2 相關(guān)數(shù)據(jù)及來源

本文在上述DEA模型結(jié)構(gòu)下，對不同乳制品企業(yè)的供應(yīng)鏈績效做對比研究，文中所有數(shù)據(jù)來源于各企業(yè)不同時(shí)間段年報(bào)，整理后分別得到三家乳制品企業(yè)投入產(chǎn)出數(shù)據(jù)如表1～表3所示。本文將利用表1～表3的數(shù)據(jù)，分別對三家處于不同梯級(jí)的乳制品企業(yè)進(jìn)行績效對比研究。

2 模型運(yùn)算及相關(guān)分析

在乳制品供應(yīng)鏈中，任何環(huán)節(jié)對于乳制品企業(yè)來說都同等重要，因此本文假設(shè)生產(chǎn)加工環(huán)節(jié)與儲(chǔ)運(yùn)銷售環(huán)節(jié)的權(quán)重相同，并通過DEAP2.1軟件對兩個(gè)環(huán)節(jié)分別進(jìn)行效率測度分析。

2.1 生產(chǎn)加工環(huán)節(jié)績效分析

分別將三家企業(yè)生產(chǎn)加工環(huán)節(jié)的表中數(shù)據(jù)代入式（1），并選擇軟件中的VRS模型對所有決策單元的生產(chǎn)加工環(huán)節(jié)進(jìn)行評價(jià)分析，不同企業(yè)該環(huán)節(jié)的分析數(shù)據(jù)分別如表4～表6所示。

在乳制品供應(yīng)鏈生產(chǎn)加工階段，由表中技術(shù)效率值可知，A企業(yè)中決策單元（Decision Making Units，DMU）實(shí)現(xiàn)DEA有效的包括DMU1和DMU4，B企業(yè)中DMU4、DMU5、DMU7滿足DEA有效，C企業(yè)DEA有效的決策單元為DMU3、DMU4、DMU6和DMU11，此外其余的決策單元技術(shù)效率值小于1，說明存在投入冗余或產(chǎn)出不足的情形。在該環(huán)節(jié)造成企業(yè)間效率差異的原因主要是企業(yè)間純技術(shù)效率與規(guī)模效率的差異，A企業(yè)是由于規(guī)模效率逐年降低所導(dǎo)致，B企業(yè)則是由于純技術(shù)效率的后勁不足致使企業(yè)近年的技術(shù)效率停滯不前，C企業(yè)在該階段相比其他企業(yè)效率值更加穩(wěn)定，雖然偶有波動(dòng)，但總體呈現(xiàn)上升趨勢。因此，A、B企業(yè)應(yīng)根據(jù)自身不足合理進(jìn)行調(diào)整;減少投入，合理配置資源;加強(qiáng)管理，提高技術(shù)能力。C企業(yè)在保持現(xiàn)有狀態(tài)的前提下，應(yīng)追求進(jìn)一步穩(wěn)步上升，扭轉(zhuǎn)效率低下的情形。

為進(jìn)一步了解三家企業(yè)在該環(huán)節(jié)的績效演變趨勢以及對比情況，繪制企業(yè)隨時(shí)間變化的技術(shù)效率趨勢圖，如圖3所示。

由圖3可知，A、B企業(yè)在生產(chǎn)加工環(huán)節(jié)的技術(shù)效率隨時(shí)間變化波動(dòng)較大，而技術(shù)效率是純技術(shù)效率與規(guī)模效率的乘積，因此A、B企業(yè)在該環(huán)節(jié)長時(shí)間的低效狀態(tài)是由這兩方面的原因造成的。對于A企業(yè)而言，技術(shù)效率無效是由規(guī)模效率無效引起，B企業(yè)則是純技術(shù)效率無效所導(dǎo)致。圖3顯示B企業(yè)整體效率高于A企業(yè)，這說明規(guī)模無效所引起的技術(shù)無效大于純技術(shù)效率無效引起的技術(shù)無效。值得注意的是，C企業(yè)技術(shù)效率隨時(shí)間變化的趨勢較為平緩，且達(dá)到有效狀態(tài)的決策單元要大于另外兩家企業(yè)。因此，在生產(chǎn)加工環(huán)節(jié)，A、B企業(yè)應(yīng)該向C企業(yè)看齊，合理地改進(jìn)自身不足，并提升技術(shù)效率，C企業(yè)也應(yīng)常存憂患意識(shí)，以A、B企業(yè)為警示，追求技術(shù)效率進(jìn)一步穩(wěn)定、有效發(fā)展。

2.2 儲(chǔ)運(yùn)銷售環(huán)節(jié)績效分析

將不同企業(yè)儲(chǔ)運(yùn)銷售環(huán)節(jié)的投入產(chǎn)出數(shù)據(jù)代入式（2），利用軟件求解，最終得到該環(huán)節(jié)結(jié)果如表7～表9所示。

為了解不同企業(yè)在該環(huán)節(jié)技術(shù)效率的變化情況，繪制企業(yè)隨時(shí)間變化的技術(shù)效率趨勢圖，如圖4所示。在儲(chǔ)運(yùn)銷售環(huán)節(jié)，A企業(yè)達(dá)到有效狀態(tài)的決策單元包括DMU7、DMU9和DMU10，B企業(yè)包括DMU9、DMU10，并且A、B企業(yè)的技術(shù)效率隨著時(shí)間逐漸增強(qiáng)，雖稍有波動(dòng)，但整體呈現(xiàn)上升趨勢。這說明A、B企業(yè)在該環(huán)節(jié)的資源配置合理，管理技術(shù)能力達(dá)標(biāo)。而C企業(yè)僅DMU1處于有效狀態(tài)，由圖4可知，C企業(yè)在2013年前后，技術(shù)效率呈現(xiàn)先降低后升高的變化，而導(dǎo)致技術(shù)效率降低的原因是純技術(shù)效率降低，因此C企業(yè)應(yīng)該進(jìn)一步提高自身的技術(shù)管理能力，例如冷鏈運(yùn)輸技術(shù)、生產(chǎn)加工技術(shù)等，以追求技術(shù)效率的穩(wěn)步增長。

由圖4可知，在儲(chǔ)運(yùn)銷售環(huán)節(jié)，A、B企業(yè)的技術(shù)效率隨時(shí)間遞增，且A企業(yè)的整體效率高于B企業(yè)，有效決策單元的個(gè)數(shù)也多于B企業(yè)，而C企業(yè)在2013年之前，其技術(shù)效率一直在減小，2013年之后技術(shù)效率則逐漸遞增，這是因其純技術(shù)效率變化造成的，表明C企業(yè)管理者在2013年做了相應(yīng)的技術(shù)調(diào)整，促使技術(shù)效率不斷升高。因此，在實(shí)際管理過程中，A、B企業(yè)在該環(huán)節(jié)應(yīng)盡量保持現(xiàn)狀，并追求技術(shù)效率的連續(xù)穩(wěn)定提升，C企業(yè)也應(yīng)繼續(xù)改進(jìn)自身的技術(shù)能力，使其技術(shù)效率達(dá)到有效狀態(tài)。

此外，為了解本文研究乳制品供應(yīng)鏈中哪一階段對其整體效率影響較大時(shí)，將圖1中所有投入作為整體投入，所有產(chǎn)出作為整體產(chǎn)出，并同樣采用DEAP2.1軟件對三家企業(yè)進(jìn)行整體績效分析，得到不同企業(yè)整體技術(shù)效率值如表10所示。

根據(jù)乳制品企業(yè)供應(yīng)鏈整體技術(shù)效率分析可知，A企業(yè)相對于另外兩家企業(yè)的供應(yīng)鏈整體效率更加有效突出，大部分決策單元均實(shí)現(xiàn)了DEA有效。B企業(yè)隨著時(shí)間推移，整體效率一直在提高，但在某些時(shí)刻還是略低于A企業(yè)。C企業(yè)在早年間也能連續(xù)保持DEA有效的狀態(tài)，然而后勁不足，近些年未能維持自身的有效狀態(tài)。

綜上所述，對比各企業(yè)供應(yīng)鏈整體效率與各個(gè)環(huán)節(jié)的效率表現(xiàn)，發(fā)現(xiàn)兩者并不相同，進(jìn)一步分析各企業(yè)生產(chǎn)加工環(huán)節(jié)和儲(chǔ)運(yùn)銷售環(huán)節(jié)的技術(shù)效率與企業(yè)整體技術(shù)效率的關(guān)系可知，對于A、B兩家企業(yè)而言，儲(chǔ)運(yùn)銷售環(huán)節(jié)的效率情況對整體效率的影響最大，而對于C企業(yè)來說，生產(chǎn)加工環(huán)節(jié)的效率情況對整體效率影響最大。

因此，乳制品企業(yè)管理者在制定相關(guān)政策時(shí)，需要從供應(yīng)鏈的各個(gè)環(huán)節(jié)出發(fā)去考慮并解決問題，而不能只關(guān)注供應(yīng)鏈的始終，僅從起始投入與最終產(chǎn)出來研究供應(yīng)鏈效率問題。管理者在進(jìn)行研究時(shí)，應(yīng)根據(jù)供應(yīng)鏈各個(gè)環(huán)節(jié)效率對整體效率的影響大小，設(shè)定不同的權(quán)重進(jìn)行分析。以本文研究企業(yè)為例，A、B企業(yè)應(yīng)設(shè)置生產(chǎn)加工環(huán)節(jié)權(quán)重大于儲(chǔ)運(yùn)銷售環(huán)節(jié)權(quán)重，C企業(yè)反之。

3 結(jié) 論

本文為研究長時(shí)間以來不同乳制品企業(yè)關(guān)于供應(yīng)鏈各階段的效率差異，構(gòu)建由生產(chǎn)加工和儲(chǔ)運(yùn)銷售兩個(gè)環(huán)節(jié)構(gòu)成的兩階段DEA模型，并將整理得出的三家乳制品企業(yè)各環(huán)節(jié)投入產(chǎn)出數(shù)據(jù)代入軟件進(jìn)行研究求解，最后得出的結(jié)論總結(jié)為以下幾個(gè)方面。

（1）在生產(chǎn)加工環(huán)節(jié)，A、B兩個(gè)企業(yè)的技術(shù)效率隨時(shí)間波動(dòng)較大，且都呈現(xiàn)出早期技術(shù)效率大于中后期技術(shù)效率的情形，原因卻并不相同。對于A企業(yè)而言，是因規(guī)模效率隨時(shí)間一直降低而使得技術(shù)效率越來越低，B企業(yè)則是由純技術(shù)效率的降低所導(dǎo)致。因此，A企業(yè)應(yīng)該減少投入，合理配置資源以實(shí)現(xiàn)最大產(chǎn)出;B企業(yè)應(yīng)進(jìn)一步提高企業(yè)的管理、技術(shù)能力以實(shí)現(xiàn)DEA有效。C企業(yè)的技術(shù)效率隨時(shí)間波動(dòng)不大，說明C企業(yè)在純技術(shù)效率與規(guī)模效率的管理上都十分適宜，值得其他企業(yè)借鑒與學(xué)習(xí)。

（2）在儲(chǔ)運(yùn)銷售環(huán)節(jié)，A、B企業(yè)隨著時(shí)間增長，技術(shù)效率不斷升高，說明無論是純技術(shù)效率還是規(guī)模效率都隨時(shí)間不斷地改進(jìn)，從而使得DEA有效的決策單元越來越多。然而C企業(yè)的技術(shù)效率在2013年之前波動(dòng)不斷，整體呈現(xiàn)下降趨勢，2013年之后技術(shù)效率則一直在升高。這是由其純技術(shù)效率的波動(dòng)所導(dǎo)致，因此C企業(yè)在未來的發(fā)展過程中要注重提高企業(yè)在該環(huán)節(jié)的技術(shù)支撐能力。

（3）分析企業(yè)各環(huán)節(jié)的技術(shù)效率對供應(yīng)鏈整體技術(shù)效率的影響可知，在A、B企業(yè)，儲(chǔ)運(yùn)銷售環(huán)節(jié)對于整體技術(shù)效率的影響要遠(yuǎn)大于生產(chǎn)加工環(huán)節(jié)，C企業(yè)與之相反，生產(chǎn)加工環(huán)節(jié)對于整體技術(shù)的效率提高的影響所占比重更大。因此，在實(shí)際管理過程中，要全面考慮供應(yīng)鏈的各個(gè)環(huán)節(jié)，并根據(jù)不同企業(yè)各環(huán)節(jié)對整體效率影響的不同，設(shè)置不同的權(quán)重，以便更好地進(jìn)行研究分析。

參考文獻(xiàn)

[1] 王磊，李翠霞，王澤民.合作伙伴特性對乳制品供應(yīng)鏈合作關(guān)系穩(wěn)定性的影響——基于質(zhì)量安全視角的實(shí)證研究[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2019（7）：104-114.

[2] 楊琦，裴磊，魏旭明.中國乳制品銷售環(huán)節(jié)安全監(jiān)管狀況影響因素研究——基于主成分因子分析和二項(xiàng)Logistic回歸[J].南京工業(yè)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2019，18（3）：90-101.

[3] 桂冬冬，周燦，潘喜悅，等.多層次食品供應(yīng)鏈網(wǎng)絡(luò)均衡研究[J].糧食科技與經(jīng)濟(jì)，2017，42（6）：30-34.

[4] 晚春東，王雨佳.供應(yīng)鏈視角下食品質(zhì)量安全風(fēng)險(xiǎn)調(diào)控投資研究[J].科技管理研究，2019，39（5）：215-221.

[5] 劉永勝，王荷麗，徐廣姝.食品供應(yīng)鏈安全風(fēng)險(xiǎn)博弈分析[J].經(jīng)濟(jì)問題，2018（1）：57-64+90.

[6] 何靜，楊翼.物聯(lián)網(wǎng)環(huán)境下的乳制品供應(yīng)鏈質(zhì)量安全風(fēng)險(xiǎn)管理研究[J].中國乳品工業(yè)，2019，47（2）：43-47.

[7] Kirilova E G，Vaklieva-Bancheva N G.Environmentally friendly management of dairy supply chain for designing a green products portfolio[J].Journal of Cleaner Production，2017（167）：493-504.

[8] Njage P M K，Opiyo B，Wangoh J，et al.Scale of production and implementation of food safety programs influence the performance of current food safety management systems：Case of dairy processors[J].Food Control，2018（85）： 85-97.

[9] Herzog A，Winckler C，Zollitsch W.In pursuit of sustainability in dairy farming：A review of interdependent effects of animal welfare improvement and environmental impact mitigation[J].Agriculture，Ecosystems & Environment，2018（267）：174-187.

[10] Kazancoglu Y，Ozkan-Ozen Y D，Ozbiltekin M. Minimizing losses in milk supply chain with sustainability：An example from an emerging economy[J].Resources， Conservation and Recycling，2018（139）：270-279.

[11] 魏艷驕，朱晶.乳制品進(jìn)口對中國不同規(guī)模奶牛養(yǎng)殖效率的影響[J].資源科學(xué)，2019，41（8）：1475-1487.

[12] 陳吉銘，劉芳，何忠偉，等.京津冀乳制品冷鏈物流效率與影響因素研究[J].中國畜牧雜志，2017，53（11）：123-128+143.

[13] 劉剛橋，師建華，張慶平.供給側(cè)改革背景下廣東省物流業(yè) 全要素生產(chǎn)率實(shí)證研究—基于DEA-Malmquist模型[J].物流技術(shù)，2019，38（9）：90-95.

[14] 方煒，楊步.基于DEA方法的企業(yè)綠色供應(yīng)鏈運(yùn)營效率評價(jià)研究[J].工業(yè)技術(shù)經(jīng)濟(jì)，2017，36（12）：19-26.

[15] 羅蓉，袁碧蓉.利率市場化對我國商業(yè)銀行效率的影響研究—基于籌資、投資兩階段DEA模型效率測度[J].湘潭大學(xué)學(xué)報(bào)（哲學(xué)社會(huì)科學(xué)版），2017，41（2）：71-77.

[16] 孟維站，李春艷，石曉冬.中國高技術(shù)產(chǎn)業(yè)創(chuàng)新效率分階段分析—基于三階段DEA模型[J].宏觀經(jīng)濟(jì)研究，2019 （2）：78-91.

[17] Izadikhah M，Saen R F.Assessing sustainability of supply chains by chance-constrained two-stage DEA model in the presence of undesirable factors[J].Computers & Operations Research，2018（100）：343-367.

[18] Cecchini L，Venanzi S，Pierri A，et al.Environmental efficiency analysis and estimation of CO2 abatement costs in dairy cattle farms in Umbria （Italy）：A SBM-DEA model with undesirable output[J].Journal of Cleaner Production，2018（197）：895-907.