• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    配額-交易機(jī)制下動(dòng)態(tài)批量生產(chǎn)和減排投資策略研究

    2016-05-16 02:28:46帝,劍,
    中國(guó)管理科學(xué) 2016年4期
    關(guān)鍵詞:總成本配額生產(chǎn)商

    黃 帝, 陳 劍, 周 泓

    (1. 北京交通大學(xué)經(jīng)濟(jì)管理學(xué)院,北京 100044;2. 清華大學(xué)經(jīng)濟(jì)管理學(xué)院,北京 100084;3. 北京航空航天大學(xué)經(jīng)濟(jì)管理學(xué)院,北京 100191)

    ?

    配額-交易機(jī)制下動(dòng)態(tài)批量生產(chǎn)和減排投資策略研究

    黃 帝1, 2, 陳 劍2, 周 泓3

    (1. 北京交通大學(xué)經(jīng)濟(jì)管理學(xué)院,北京 100044;2. 清華大學(xué)經(jīng)濟(jì)管理學(xué)院,北京 100084;3. 北京航空航天大學(xué)經(jīng)濟(jì)管理學(xué)院,北京 100191)

    隨著我國(guó)碳排放交易市場(chǎng)的建立和發(fā)展,在碳排放約束下逐步降低單位產(chǎn)出的碳排放水平成為企業(yè)生產(chǎn)經(jīng)營(yíng)管理中的中長(zhǎng)期約束性目標(biāo)。本文在一個(gè)多周期決策模型中研究了配額-交易機(jī)制下企業(yè)的最優(yōu)動(dòng)態(tài)批量生產(chǎn)、碳排放權(quán)交易和減排投資聯(lián)合決策問(wèn)題。生產(chǎn)商在整個(gè)決策周期期初決定是否進(jìn)行減排投資以及投資規(guī)模,根據(jù)每個(gè)周期的生產(chǎn)計(jì)劃決定減排設(shè)備的運(yùn)行計(jì)劃。根據(jù)節(jié)能減排技術(shù)的特點(diǎn),本文假設(shè)生產(chǎn)商運(yùn)行減排設(shè)備時(shí)不僅降低了產(chǎn)品的單位生產(chǎn)碳排放量,而且降低了產(chǎn)品的單位生產(chǎn)成本。本文基于廣義Benders分解法對(duì)模型進(jìn)行了最優(yōu)性分析,得到了最優(yōu)生產(chǎn)決策和最優(yōu)減排投資決策的一些基本性質(zhì),并通過(guò)數(shù)值實(shí)驗(yàn)分析了碳排放配額和碳排放權(quán)價(jià)格對(duì)生產(chǎn)商總成本、總排放以及減排投資決策的影響。本文的數(shù)值實(shí)驗(yàn)分析結(jié)果發(fā)現(xiàn):(1)當(dāng)碳交易市場(chǎng)上的碳排放權(quán)充足時(shí),減少碳排放配額或改變碳排放配額的分配方式并不能影響生產(chǎn)商的碳排放水平;(2)碳排放權(quán)價(jià)格是影響生產(chǎn)商的碳排放水平和減排投資規(guī)模的關(guān)鍵因素;(3)隨著碳排放權(quán)價(jià)格的上升,即使擁有足夠的碳排放配額,生產(chǎn)商仍會(huì)不斷提高減排投資規(guī)模以獲得減排收益。研究結(jié)果對(duì)碳排放交易體系下生產(chǎn)企業(yè)進(jìn)行減排技術(shù)投資具有較強(qiáng)的管理啟示。

    動(dòng)態(tài)批量;減排投資;配額-交易;廣義Benders分解

    1 引言

    近年來(lái),環(huán)境保護(hù)和經(jīng)濟(jì)發(fā)展之間的矛盾愈加突出。一方面,全球氣候變暖使人們認(rèn)識(shí)到經(jīng)濟(jì)、社會(huì)和環(huán)境協(xié)調(diào)發(fā)展的重要性,以低能耗、低污染為理念的低碳經(jīng)濟(jì)成為人類(lèi)社會(huì)可持續(xù)發(fā)展的必然選擇。另一方面,隨著經(jīng)濟(jì)的發(fā)展,全球能源需求仍在繼續(xù)增加。國(guó)際能源署(IEA)預(yù)測(cè),全球能源需求量到2035年將至少增長(zhǎng)三分之一。工業(yè)是國(guó)民經(jīng)濟(jì)的重要組成部分,是推動(dòng)經(jīng)濟(jì)增長(zhǎng)的主要?jiǎng)恿?,也是能源消費(fèi)及二氧化碳等溫室氣體排放的主要領(lǐng)域。為實(shí)現(xiàn)經(jīng)濟(jì)持續(xù)健康發(fā)展,降低工業(yè)生產(chǎn)中的二氧化碳排放量成為國(guó)民經(jīng)濟(jì)發(fā)展中的長(zhǎng)期約束性指標(biāo),也成為當(dāng)前學(xué)術(shù)研究的熱點(diǎn)。

    碳排放規(guī)則是指國(guó)際社會(huì)為應(yīng)對(duì)全球氣候變暖給人類(lèi)生存和發(fā)展帶來(lái)的嚴(yán)峻挑戰(zhàn)而采取的一系列控制二氧化碳等溫室氣體排放的政策措施。配額—交易機(jī)制作為一種二氧化碳排放總量控制手段,被中國(guó)、美國(guó)、澳大利亞和歐盟等國(guó)家和地區(qū)廣泛使用。它通過(guò)總量控制和市場(chǎng)調(diào)節(jié)相結(jié)合的手段,引導(dǎo)和鼓勵(lì)企業(yè)在節(jié)能減排技術(shù)和可再生能源領(lǐng)域進(jìn)行投資,實(shí)現(xiàn)經(jīng)濟(jì)社會(huì)的總體減排。在碳排放權(quán)交易市場(chǎng)中,碳排放權(quán)被賦予特定的產(chǎn)權(quán),可以在市場(chǎng)上進(jìn)行轉(zhuǎn)讓、交換等交易活動(dòng),通過(guò)碳排放權(quán)價(jià)格的調(diào)節(jié)來(lái)實(shí)現(xiàn)環(huán)境資源的合理配置。1997年12月制定的國(guó)際公約《京都議定書(shū)》提出了旨在減排的三個(gè)靈活合作機(jī)制:聯(lián)合履約、清潔發(fā)展和碳排放權(quán)交易機(jī)制。此后,美國(guó)、澳大利亞、歐盟和中國(guó)等國(guó)家和地區(qū)都相繼建立起了碳排放權(quán)交易市場(chǎng)。以配額—交易機(jī)制為核心的歐洲碳排放權(quán)交易體系(EU-ETS)自2005年至今已運(yùn)行三個(gè)階段,在減少二氧化碳等溫室氣體排放方面取得了初步成效。我國(guó)自2012年以來(lái)在北京、天津、上海、重慶、深圳、廣東和湖北等七個(gè)省市逐步開(kāi)展碳排放交易試點(diǎn)工作。截止2014年6月19日,七個(gè)碳排放權(quán)交易試點(diǎn)省市已全部開(kāi)始實(shí)際交易。目前,廣東碳交易市場(chǎng)已發(fā)展成為僅次于歐盟、韓國(guó)的全球第三大碳交易市場(chǎng)。

    碳排放規(guī)則的實(shí)施給企業(yè)帶來(lái)了諸多管理挑戰(zhàn),逐步降低碳排放水平成為企業(yè)生產(chǎn)經(jīng)營(yíng)管理中的中長(zhǎng)期約束性指標(biāo)。本文在一個(gè)多周期決策模型中研究配額—交易機(jī)制下企業(yè)的最優(yōu)生產(chǎn)計(jì)劃和減排計(jì)劃,分析減排投資成本、碳排放價(jià)格、產(chǎn)品需求等參數(shù)對(duì)企業(yè)生產(chǎn)和減排決策的影響,著力解決以下幾個(gè)企業(yè)運(yùn)營(yíng)管理中的關(guān)鍵問(wèn)題:(1) 企業(yè)何時(shí)投資減排設(shè)備?(2) 企業(yè)何時(shí)運(yùn)行減排設(shè)備?(3) 企業(yè)如何安排生產(chǎn)計(jì)劃?這些都是在企業(yè)管理實(shí)踐中迫切需要解決,但目前又尚未進(jìn)行深入研究的管理問(wèn)題。對(duì)于這些問(wèn)題的研究無(wú)論是從理論層面還是從企業(yè)運(yùn)作層面都具有非常重要的意義

    目前對(duì)于低碳背景下企業(yè)生產(chǎn)和減排投資問(wèn)題的研究主要集中在兩個(gè)方向。第一個(gè)研究方向是碳排放約束下的批量生產(chǎn)問(wèn)題研究。Hua Guowei等[1]研究了配額—交易機(jī)制下碳排放權(quán)交易價(jià)格和碳排放權(quán)配額對(duì)企業(yè)的經(jīng)濟(jì)訂貨批量、總排放和總成本的影響。Benjaafar等[2]在動(dòng)態(tài)批量生產(chǎn)模型中研究了碳稅、配額—交易、配額—抵消等不同的碳政策下單個(gè)企業(yè)及供應(yīng)鏈的多周期生產(chǎn)庫(kù)存決策問(wèn)題。他們還通過(guò)數(shù)值實(shí)驗(yàn)分析了碳排放權(quán)配額、碳政策、節(jié)能減排技術(shù)和供應(yīng)鏈聯(lián)合減排等因素對(duì)企業(yè)總排放和總成本的影響。Absi等[3]研究了動(dòng)態(tài)批量模型在周期性單位生產(chǎn)碳排放約束、累積周期單位生產(chǎn)碳排放約束、滾動(dòng)周期單位生產(chǎn)碳排放約束和全周期單位生產(chǎn)碳排放約束等約束條件下的算法時(shí)間復(fù)雜度。他們證明只有周期性單位生產(chǎn)碳排放約束下的動(dòng)態(tài)批量模型存在多項(xiàng)式時(shí)間算法。胡海菊和李勇建[4]利用系統(tǒng)動(dòng)力學(xué)仿真工具對(duì)考慮再制造和產(chǎn)品需求可替代的動(dòng)態(tài)批量生產(chǎn)問(wèn)題進(jìn)行了建模分析。戴道明[5]研究了產(chǎn)品動(dòng)態(tài)定價(jià)與批量生產(chǎn)的聯(lián)合決策問(wèn)題,分析了不同碳政策情形下碳排放配額、碳排放交易價(jià)格和碳稅對(duì)企業(yè)利潤(rùn)和碳排放水平的影響。

    第二個(gè)研究方向是企業(yè)碳減排投資決策問(wèn)題研究。王明喜等[6]指出,企業(yè)在生產(chǎn)環(huán)節(jié)可以通過(guò)減排設(shè)備投資、減排研發(fā)投資和能源投入結(jié)構(gòu)調(diào)整三個(gè)方面來(lái)降低其單位產(chǎn)值的二氧化碳排放量。趙道致和張學(xué)強(qiáng)[7]研究了考慮總成本和總排放量?jī)蓚€(gè)優(yōu)化目標(biāo)的供應(yīng)鏈網(wǎng)絡(luò)設(shè)計(jì)和優(yōu)化問(wèn)題。他們假設(shè)生產(chǎn)商有多種可供選擇的技術(shù)手段,分別對(duì)應(yīng)不同的單位產(chǎn)品碳減排量和減排投資額。謝鑫鵬和趙道致[8-9]研究了單周期兩級(jí)供應(yīng)鏈的定價(jià)和減排決策,但他們?cè)谘芯恐屑僭O(shè)減排投資對(duì)產(chǎn)品的單位生產(chǎn)成本沒(méi)有影響。杜少甫等[10]認(rèn)為企業(yè)的碳減排來(lái)自于對(duì)排放物的凈化處理,他們分別研究了確定的凈化水平和可變的凈化水平下企業(yè)的生產(chǎn)決策和碳排放權(quán)交易決策。此后,杜少甫等[10]的工作被何大義和馬洪云[11]和侯玉梅等[12]分別擴(kuò)展到不確定性需求情形下和雙寡頭企業(yè)情形下。陳伯成和李英杰[13]對(duì)杜少甫等[10]的工作進(jìn)行了補(bǔ)充研究,豐富了對(duì)企業(yè)最優(yōu)決策結(jié)果的分析和參數(shù)敏感性分析。Islegen等[14]對(duì)碳捕捉與存儲(chǔ)技術(shù)的投資決策進(jìn)行了經(jīng)濟(jì)性分析。Jiang Yan等[15]對(duì)比分析了不同碳政策下企業(yè)的總利潤(rùn)、總排放量、減排投資額及投資時(shí)機(jī)等決策。駱瑞玲等[16]在一個(gè)兩級(jí)供應(yīng)鏈中研究了配額—交易機(jī)制下供應(yīng)鏈減排技術(shù)投資決策問(wèn)題,分析了不同決策情形下消費(fèi)者碳足跡敏感系數(shù)、碳排放配額和碳減排成本系數(shù)對(duì)供應(yīng)鏈碳排放總量、供應(yīng)鏈利潤(rùn)、碳減排投資水平等的影響。王明喜等[17]給出了碳排放約束下我國(guó)企業(yè)在生產(chǎn)技術(shù)、清潔能源和節(jié)能管理三條減排路徑上的最優(yōu)減排投資決策,分析了“強(qiáng)制減排”和“自愿減排”兩種措施下不同的碳排放配額分配方式和企業(yè)間碳交易模式對(duì)企業(yè)實(shí)施最優(yōu)減排投資的影響。陳伯成等[18]假設(shè)政府對(duì)企業(yè)排放凈化量進(jìn)行補(bǔ)貼、對(duì)企業(yè)排放量進(jìn)行懲罰,給出了不同碳政策情形下的最優(yōu)凈化水平和最優(yōu)產(chǎn)量。周艷菊等[19]在一個(gè)兩級(jí)供應(yīng)鏈中研究了減排成本分擔(dān)契約對(duì)產(chǎn)品訂貨量和銷(xiāo)售價(jià)格的影響,以及消費(fèi)者低碳認(rèn)知和低碳情感變化對(duì)最優(yōu)決策的影響。

    本文與前述文獻(xiàn)的主要區(qū)別在于:1)研究碳排放約束下企業(yè)的多周期生產(chǎn)、碳排放交易和減排聯(lián)合決策問(wèn)題。而已有文獻(xiàn)或者只考慮碳排放約束下企業(yè)的多周期生產(chǎn)決策問(wèn)題[2-3,5],或者只考慮單周期的生產(chǎn)和減排聯(lián)合決策問(wèn)題[10-13,16,18-19],少有兩者兼顧。2)結(jié)合國(guó)內(nèi)實(shí)踐,主要研究企業(yè)的減排投資通過(guò)降低單位產(chǎn)出的耗能量來(lái)實(shí)現(xiàn)降低單位產(chǎn)出的碳排放量,即技術(shù)減排路徑[17]。而已有文獻(xiàn)認(rèn)為企業(yè)的減排投資僅僅降低單位碳排放量,不影響產(chǎn)品的單位生產(chǎn)成本(包含能耗成本)[8-13,16,18-19]。3)考慮企業(yè)在某一周期的減排投資會(huì)在后續(xù)周期發(fā)揮減排效果[17,20]。而大多數(shù)文獻(xiàn)在研究中假設(shè)企業(yè)的減排投資僅在當(dāng)前決策周期產(chǎn)生減排效果[10-13,16,18-19]。

    鑒于此,本文在確定需求環(huán)境下研究多周期情形下的動(dòng)態(tài)批量生產(chǎn)、碳排放權(quán)交易和減排投資聯(lián)合決策問(wèn)題。生產(chǎn)商在整個(gè)決策周期期初決定是否進(jìn)行減排投資和投資規(guī)模,根據(jù)每個(gè)周期的生產(chǎn)計(jì)劃決定減排設(shè)備的運(yùn)行計(jì)劃。根據(jù)節(jié)能減排技術(shù)的特點(diǎn),本文假設(shè)生產(chǎn)商運(yùn)行減排設(shè)備時(shí)不僅降低了產(chǎn)品的單位生產(chǎn)碳排放量,而且降低了產(chǎn)品的單位生產(chǎn)成本。本文與Benjaafar等[2]的研究不同之處在于,本文在模型中額外考慮了生產(chǎn)商的減排投資決策,而且基于廣義Benders分解得到了生產(chǎn)和減排投資決策的最優(yōu)性質(zhì)。分析發(fā)現(xiàn),配額—交易機(jī)制下同時(shí)考慮碳排放權(quán)交易和減排投資時(shí)的最優(yōu)生產(chǎn)決策與經(jīng)典批量生產(chǎn)問(wèn)題的最優(yōu)決策相似,即若生產(chǎn)商在某個(gè)周期生產(chǎn)產(chǎn)品,則其生產(chǎn)量必定滿(mǎn)足前后若干個(gè)周期的總需求。此外,本文的數(shù)值實(shí)驗(yàn)還得到一些與Benjaafar等[2]的數(shù)值實(shí)驗(yàn)不同的觀察結(jié)果:(1)當(dāng)碳排放配額固定時(shí),隨著碳價(jià)格的升高生產(chǎn)商的碳排放水平不斷下降,而總成本先增加后減少。(2)當(dāng)碳排放權(quán)價(jià)格固定時(shí),生產(chǎn)商的減排投資規(guī)模不受碳排放配額的影響;當(dāng)碳排放配額固定時(shí),生產(chǎn)商的減排投資規(guī)模隨碳排放權(quán)價(jià)格的增加不斷上升。(3)當(dāng)碳排放權(quán)價(jià)格與碳排放配額呈負(fù)相關(guān)關(guān)系時(shí),隨著碳排放權(quán)價(jià)格的上升,生產(chǎn)商的碳排放水平不斷下降,而總成本先減少后增加。(4)碳排放配額按周期平均分配與一次性分配的結(jié)果相同,即碳排放權(quán)的初始分配方式并不影響生產(chǎn)商的決策。

    2 問(wèn)題描述和建模

    2.1 問(wèn)題描述和符號(hào)定義

    在一個(gè)含有T個(gè)周期的決策問(wèn)題中,生產(chǎn)商在第t(t=1,…,T)個(gè)周期的市場(chǎng)需求為dt,固定生產(chǎn)成本為St。第t個(gè)周期的生產(chǎn)量為Qt,產(chǎn)品的單位生產(chǎn)成本為ct,yt為一個(gè)二進(jìn)制變量,第t個(gè)周期生產(chǎn)時(shí)為1,不生產(chǎn)時(shí)為0。第t個(gè)周期期末的庫(kù)存水平為It,單位庫(kù)存成本為ht。第t個(gè)周期期末的延期交貨量為Bt,延期交貨的單位成本為bt。

    2.2 模型的數(shù)學(xué)表示

    生產(chǎn)商的批量生產(chǎn)、碳排放權(quán)交易和減排投資聯(lián)合決策模型如下:

    (1)

    s.t.It-1+Qt+Bt-Bt-1-It=dt,t=1,…,T

    (2)

    (3)

    (4)

    yt∈{0,1},t=1,…,T

    (5)

    0≤xt≤min(Qt,a),t=1,…,T

    (6)

    (7)

    Qt,It,Bt≥0,t=1,…,T

    (8)

    目標(biāo)函數(shù)(1)中第一項(xiàng)是生產(chǎn)啟動(dòng)成本;第二項(xiàng)是生產(chǎn)成本;第三項(xiàng)是延期交貨成本;第四項(xiàng)是庫(kù)存成本;第五項(xiàng)是碳排放權(quán)交易成本(負(fù)值表示收益);第六項(xiàng)是減排設(shè)備運(yùn)行成本;第七項(xiàng)是運(yùn)行減排設(shè)備節(jié)約的生產(chǎn)成本;最后一項(xiàng)是減排設(shè)備投資成本。約束條件(2)是庫(kù)存平衡方程;約束條件(3)是整個(gè)決策周期內(nèi)的碳排放平衡約束;約束條件(4)表示產(chǎn)量Qt和yt之間的邏輯關(guān)系;約束條件(5) 表示yt是0-1變量;約束條件(6)表示每周期的減排處理量不得超過(guò)其投資規(guī)模和該周期的生產(chǎn)量;約束條件(7)和(8)是變量的非負(fù)性約束。

    3 模型的廣義Benders分解

    碳排放約束下的動(dòng)態(tài)批量模型及其擴(kuò)展模型一般都是NP難問(wèn)題[3],一般的處理方法是將原問(wèn)題進(jìn)行分解,以降低求解規(guī)模和求解難度。魯奎等[21-22]提出了一種約束松弛的模型分解方法:對(duì)耦合約束條件進(jìn)行拉格朗日松弛,然后將原問(wèn)題直接拆分為兩個(gè)獨(dú)立的子問(wèn)題。但其缺點(diǎn)是被松弛的約束條件對(duì)應(yīng)的拉格朗日乘子被拆分到不同的子問(wèn)題中,影響了解的質(zhì)量。有鑒于此,本文使用廣義Benders分解法對(duì)原問(wèn)題進(jìn)行模型分解,以分析最優(yōu)解的性質(zhì)。

    由于碳排放權(quán)的市場(chǎng)價(jià)格在整個(gè)決策周期內(nèi)是固定不變的,不存在投機(jī)套利的機(jī)會(huì),因此生產(chǎn)商的最優(yōu)碳排放權(quán)交易數(shù)量應(yīng)等于其短缺或剩余的碳排放權(quán)數(shù)量。故約束條件(3)等價(jià)于:

    (9)

    將約束條件(9)代入到(1)中,可得

    (10)

    s.t. (2), (4), (5), (6), (8)

    廣義Benders分解法[23]的基本思想是,對(duì)于給定的復(fù)雜變量的每一組值,容易找出剩下的優(yōu)化問(wèn)題(子問(wèn)題)的最優(yōu)解,在這些最優(yōu)解中使主問(wèn)題目標(biāo)函數(shù)值最大的那一組最優(yōu)解及其對(duì)應(yīng)的復(fù)雜變量值即為原問(wèn)題的最優(yōu)解。其中,復(fù)雜變量是指當(dāng)其給定時(shí)會(huì)使原問(wèn)題的優(yōu)化變得非常容易的那些變量。根據(jù)廣義Benders分解法,將問(wèn)題(10)分解為一個(gè)主問(wèn)題和一個(gè)子問(wèn)題,其中復(fù)雜變量定義為生產(chǎn)商的減排投資決策變量a。

    (11)

    s.t. (6)

    其中,ZQ是下列子問(wèn)題的最優(yōu)值:

    (12)

    s.t. (2), (4), (5), (8)

    子問(wèn)題(12)是一個(gè)考慮碳排放成本的允許缺貨的動(dòng)態(tài)批量生產(chǎn)模型,它是Wager和Whitin[24]提出的動(dòng)態(tài)經(jīng)濟(jì)批量訂貨模型的擴(kuò)展形式[25],可以在O(T2)時(shí)間內(nèi)找出問(wèn)題的最優(yōu)解[26]。子問(wèn)題(12)以其最優(yōu)目標(biāo)函數(shù)值ZQ反饋到主問(wèn)題(11)的目標(biāo)函數(shù)中,以其最優(yōu)解Qt反饋到主問(wèn)題(11)的約束條件(6)中。

    值得注意的是,子問(wèn)題(12)的最優(yōu)解不受復(fù)雜變量a取值的影響。因此,子問(wèn)題(12)的最優(yōu)解即為原問(wèn)題(1)的最優(yōu)解。子問(wèn)題(12)的目標(biāo)函數(shù)是線(xiàn)性的,約束條件(2), (4), (5), (8)組成的可行域是封閉的凸集,所以其最優(yōu)解必在該凸集的頂點(diǎn)處達(dá)到。

    定理1 子問(wèn)題(12)的最優(yōu)解滿(mǎn)足如下性質(zhì):

    (1)It-1Qt=It-1Bt=ItBt=QtBt=0,t=1,…,T;

    (2)同樣使用反證法:若最優(yōu)解不滿(mǎn)足(2),易知其必定不滿(mǎn)足(1),則其不是最優(yōu)解。

    證畢。

    定理1說(shuō)明,生產(chǎn)商的最優(yōu)生產(chǎn)策略應(yīng)該遵循以下原則:

    (1) 庫(kù)存水平為0時(shí)才進(jìn)行生產(chǎn);

    (2)有庫(kù)存時(shí)不允許缺貨;

    (3) 若在某個(gè)周期進(jìn)行生產(chǎn),則該周期不應(yīng)存在缺貨,即產(chǎn)量至少應(yīng)該滿(mǎn)足當(dāng)期需求和該周期之前的所有缺貨;

    (4) 一次生產(chǎn)的產(chǎn)量應(yīng)該等于該周期前后若干個(gè)周期的總需求。

    定義1 如果Qt>0,則稱(chēng)周期t為生產(chǎn)周期;如果It>0,則稱(chēng)周期t為庫(kù)存周期;如果Bt>0,則稱(chēng)周期t為缺貨周期。特別地,如果It=0且Bt=0,則稱(chēng)周期t為再生周期。

    在定理1給出的最優(yōu)生產(chǎn)策略中,任意兩個(gè)連續(xù)的再生周期之間存在且僅存在一個(gè)生產(chǎn)周期以及若干個(gè)庫(kù)存周期和缺貨周期。由此,可以給出子問(wèn)題(12)的動(dòng)態(tài)規(guī)劃遞推公式[27]:

    (13)

    其中,E(j)表示從周期1到周期j的最小總成本,周期i和j表示兩個(gè)連續(xù)的再生周期(i

    一般來(lái)說(shuō),生產(chǎn)商獲得碳排放權(quán)的來(lái)源有三種方式:政府分配配額、碳排放權(quán)市場(chǎng)交易和生產(chǎn)商投資減排技術(shù)[8,10,11,12]。其中,碳排放權(quán)配額是由政府管理部門(mén)分配的外生變量,其它兩項(xiàng)是由生產(chǎn)商決定的決策變量。當(dāng)問(wèn)題(1)的最優(yōu)生產(chǎn)決策變量由子問(wèn)題(12)確定后,主問(wèn)題(11)中生產(chǎn)商的最優(yōu)減排投資和運(yùn)行決策滿(mǎn)足如下性質(zhì):

    生產(chǎn)商在某個(gè)決策周期內(nèi)進(jìn)行減排設(shè)備投資決策時(shí),需要在減排投資和運(yùn)行總成本與減排總收益之間進(jìn)行權(quán)衡。其中,減排投資和運(yùn)行總成本包括減排投資成本和每個(gè)周期的減排設(shè)備運(yùn)行成本,減排總收益包括每個(gè)周期運(yùn)行減排設(shè)備節(jié)約的碳排放成本和能耗成本。定理3的管理意義是,當(dāng)減排投資的邊際收益小于邊際成本時(shí),生產(chǎn)商不進(jìn)行減排投資。

    4 數(shù)值實(shí)驗(yàn)與分析

    為了分析配額—交易機(jī)制下碳排放權(quán)價(jià)格pc和碳排放權(quán)配額C對(duì)企業(yè)生產(chǎn)和減排投資決策的影響,在數(shù)值實(shí)驗(yàn)中分別對(duì)其取不同的值,其中pc={10, 15, 20, 25, 30, 35},C={9000, 9500, 10000, 10500, 11000, 11500}。下列數(shù)值實(shí)驗(yàn)結(jié)果均通過(guò)對(duì)隨機(jī)生成的25組需求數(shù)據(jù)進(jìn)行運(yùn)算并求其平均值得到。

    圖1 碳排放配額對(duì)總成本和碳排放量的影響

    (1)固定碳排放權(quán)的市場(chǎng)價(jià)格pc=10,逐漸變化碳排放權(quán)配額C發(fā)現(xiàn),生產(chǎn)商的碳排放水平不受碳排放配額的影響,而總成本隨著碳排放權(quán)配額的減少不斷上升(圖1)。這是因?yàn)閺?10)中可以看出, 碳排放權(quán)配額C不影響生產(chǎn)商的決策, 與生產(chǎn)商的總成本之間是遞減關(guān)系。這個(gè)結(jié)果說(shuō)明當(dāng)碳交易市場(chǎng)上的碳排放權(quán)充足時(shí),減少碳排放配額并不能降低生產(chǎn)商的碳排放水平。

    圖2 碳價(jià)格對(duì)總成本和碳排放量的影響

    (2)固定碳排放配額C=10000,逐漸變化碳排放權(quán)的市場(chǎng)價(jià)格pc發(fā)現(xiàn),生產(chǎn)商的碳排放水平隨著碳排放權(quán)價(jià)格的增加不斷下降,而總成本隨著碳排放價(jià)格的增加先增加后減少(圖2)。這是因?yàn)殡S著碳排放權(quán)價(jià)格的上升,生產(chǎn)商的碳排放水平不斷下降,減排投資規(guī)模不斷上升(表2),生產(chǎn)商由碳排放權(quán)不足(支付碳排放權(quán)買(mǎi)入成本)逐漸轉(zhuǎn)變?yōu)樘寂欧艡?quán)剩余(賣(mài)出碳排放權(quán)獲得收益)。從圖2中可以看出,當(dāng)碳排放權(quán)價(jià)格大于15時(shí),生產(chǎn)商的總成本和碳排放水平隨著碳價(jià)格的升高同時(shí)降低。

    (3) 當(dāng)碳排放權(quán)價(jià)格固定時(shí),生產(chǎn)商的減排投資規(guī)模不受碳排放配額的影響(表1);當(dāng)碳排放配額固定時(shí),生產(chǎn)商的減排投資規(guī)模隨碳排放權(quán)價(jià)格的增加不斷上升(表2)。表1-2給出了兩種情況下生產(chǎn)商的實(shí)際生產(chǎn)量,減排處理量和最大減排處理能力的統(tǒng)計(jì)平均值。從表1可以看出,約96.4%的生產(chǎn)量經(jīng)過(guò)了減排處理,但同時(shí)約6.4%的減排能力被閑置。從表2可以看出,雖然在pc=20時(shí)生產(chǎn)商的減排投資規(guī)模已基本實(shí)現(xiàn)對(duì)生產(chǎn)量的100%減排處理,但是隨著碳價(jià)格的增加,生產(chǎn)商的減排投資規(guī)模仍在繼續(xù)增加,使得減排能力閑置率由5.9%增加至8.4%。這說(shuō)明過(guò)低的碳排放權(quán)價(jià)格會(huì)抑制減排投資,而過(guò)高的碳排放權(quán)價(jià)格則會(huì)引起減排投資過(guò)度,減排能力閑置。

    表1 不同碳排放配額下的實(shí)際生產(chǎn)量,減排處理量和最大減排處理能力(統(tǒng)計(jì)平均值)

    注:“實(shí)際生產(chǎn)量”指整個(gè)決策周期內(nèi)生產(chǎn)商累計(jì)生產(chǎn)的數(shù)量,它等于總需求量;“減排處理量”指整個(gè)決策周期內(nèi)減排設(shè)備累計(jì)運(yùn)行處理的數(shù)量;“最大減排處理能力”指減排設(shè)備滿(mǎn)負(fù)荷運(yùn)行時(shí)處理的數(shù)量,它等于減排投資規(guī)模乘以生產(chǎn)周期個(gè)數(shù)。

    表2 不同碳排放權(quán)價(jià)格下的實(shí)際生產(chǎn)量,減排處理量和最大減排處理能力(統(tǒng)計(jì)平均值)

    注:“實(shí)際生產(chǎn)量”, “減排處理量”和“最大減排處理能力”的定義同表1。

    (4) 當(dāng)碳排放權(quán)價(jià)格與碳排放配額呈負(fù)相關(guān)關(guān)系(如pc=125-0.009C)時(shí),生產(chǎn)商的碳排放水平隨著碳排放權(quán)價(jià)格的增加不斷下降(圖3),生產(chǎn)商的減排投資規(guī)模隨著碳排放權(quán)價(jià)格的增加不斷上升(參見(jiàn)表2)。生產(chǎn)商的總成本隨著碳排放權(quán)價(jià)格的增加先減少后增加(圖3)。從圖3中可以看出,當(dāng)碳排放權(quán)價(jià)格小于20時(shí),生產(chǎn)商的總成本和碳排放水平隨著碳價(jià)格的升高同時(shí)降低。將pc=125-0.009C代入(10)中可以看出,這是因?yàn)樯a(chǎn)商的總成本函數(shù)(14)中最后一項(xiàng)是關(guān)于pc的凸二次函數(shù):

    (14)

    s.t. (2), (4), (5), (6), (8)

    圖3 碳價(jià)格與碳配額負(fù)相關(guān)時(shí)碳價(jià)格對(duì)總成本和碳排放量的影響

    (5)當(dāng)碳排放配額按周期平均分配給生產(chǎn)商時(shí),與一次性分配碳排放配額的結(jié)果相同,生產(chǎn)商的碳排放水平和減排投資規(guī)模不受碳排放配額的影響(參見(jiàn)圖1),而總成本隨著碳排放權(quán)配額的減少而不斷上升(參見(jiàn)表1)。

    假設(shè)生產(chǎn)商每個(gè)周期得到CT的碳排放配額,則問(wèn)題(1)中的約束條件(3)調(diào)整為

    (15)

    同樣,生產(chǎn)商每周期的最優(yōu)碳排放權(quán)交易數(shù)量應(yīng)等于其短缺或剩余的碳排放權(quán)數(shù)量。故約束條件(15)等價(jià)于:

    (16)

    代入到(1)中,生產(chǎn)商的總成本函數(shù)為:

    (17)

    s.t. (2), (4), (5), (6), (8)

    從(17)中可以看出,碳排放權(quán)的周期配額CT同樣不影響生產(chǎn)商的決策。也就是說(shuō),當(dāng)存在穩(wěn)定的碳交易市場(chǎng)時(shí),碳排放權(quán)配額的分配方式并不影響生產(chǎn)商的決策。

    5 結(jié)語(yǔ)

    本文研究了一個(gè)多周期動(dòng)態(tài)需求下企業(yè)生產(chǎn)和減排投資聯(lián)合優(yōu)化決策問(wèn)題,提出了一個(gè)考慮碳排放總量約束、碳排放交易和減排投資的動(dòng)態(tài)批量?jī)?yōu)化模型?;趶V義Benders分解法對(duì)模型進(jìn)行了分解分析,得到了最優(yōu)生產(chǎn)決策和最優(yōu)減排投資決策的一些性質(zhì),并通過(guò)數(shù)值實(shí)驗(yàn)分析了碳排放配額和碳排放權(quán)價(jià)格對(duì)生產(chǎn)商總成本和總排放以及減排投資決策的影響。本文的研究結(jié)果揭示了配額—交易機(jī)制的幾個(gè)特點(diǎn):(1)當(dāng)碳交易市場(chǎng)上的碳排放權(quán)充足時(shí),減少碳排放配額或改變碳排放配額的分配方式并不能影響生產(chǎn)商的碳排放水平;(2)碳排放權(quán)價(jià)格是影響生產(chǎn)商的碳排放水平和減排投資規(guī)模的關(guān)鍵因素;(3)隨著碳排放權(quán)價(jià)格的上升,即使擁有足夠的碳排放配額,生產(chǎn)商仍會(huì)不斷提高減排投資規(guī)模以獲得減排收益。因此,對(duì)于給定的碳排放配額,總存在一個(gè)碳排放權(quán)價(jià)格“拐點(diǎn)”,在此拐點(diǎn)價(jià)格水平之上提高碳排放權(quán)價(jià)格可以同時(shí)降低生產(chǎn)商總成本和碳排放水平。(4)若減少企業(yè)的碳排放配額會(huì)引起碳排放權(quán)價(jià)格的上漲,則企業(yè)從增加減排投資規(guī)模中獲得的減排收益最終會(huì)被迅速增加的碳排放成本所抵消。因此也存在一個(gè)碳排放權(quán)價(jià)格“拐點(diǎn)”, 在此拐點(diǎn)價(jià)格水平之下提高碳排放權(quán)價(jià)格可以同時(shí)降低生產(chǎn)商總成本和碳排放水平。這些都是在已有文獻(xiàn)中尚未涉及,但又對(duì)配額—交易機(jī)制的應(yīng)用實(shí)踐具有重要意義的結(jié)論。

    [1] Hua Guowei, Cheng T C E, Wang Shouyang. Managing carbon footprints in inventory management[J]. International Journal of Production Economics, 2011, 132(2): 178-185.

    [2] Benjaafar S, Li Yanzhi, Daskin M S. Carbon footprint and the management of supply chains: Insights from simple models[J]. IEEE transactions on automation science and engineering, 2013, 10(1): 99-116.

    [3] Absi N, Dauzère-Pérès S, Kedad-Sidhoum S, et al. Lot sizing with carbon emission constraints[J]. European Journal of Operational Research, 2013, 227(1): 55-61.

    [4] 胡海菊, 李勇建. 考慮再制造和產(chǎn)品需求可替代的短生命周期產(chǎn)品動(dòng)態(tài)批量生產(chǎn)計(jì)劃問(wèn)題[J]. 系統(tǒng)工程理論與實(shí)踐, 2007, (12): 76-84.

    [5] 戴道明.碳排放約束下動(dòng)態(tài)定價(jià)與批量的聯(lián)合決策[J]. 合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 36(8): 996-1001.

    [6] 王明喜, 王明榮, 汪壽陽(yáng), 等. 最優(yōu)減排策略及其實(shí)施的理論分析[J]. 管理評(píng)論, 2010, 22(6): 42-47.

    [7] 趙道致, 張學(xué)強(qiáng). 面向碳減排投資優(yōu)化的低碳供應(yīng)鏈網(wǎng)絡(luò)設(shè)計(jì)及優(yōu)化研究[J]. 物流技術(shù), 2013, 32(3): 215-218.

    [8] 謝鑫鵬, 趙道致. 基于CDM的兩級(jí)低碳供應(yīng)鏈企業(yè)產(chǎn)品定價(jià)與減排決策機(jī)制研究[J]. 軟科學(xué), 2013, 27(5): 80-85.

    [9] 謝鑫鵬, 趙道致. 低碳供應(yīng)鏈企業(yè)減排合作策略研究[J]. 管理科學(xué), 2013, 26(3): 108-119.

    [10] 杜少甫, 董駿峰, 梁樑, 等. 考慮排放許可與交易的生產(chǎn)優(yōu)化[J]. 中國(guó)管理科學(xué), 2009, 17(3): 81-86.

    [11] 何大義, 馬洪云. 碳排放約束下企業(yè)生產(chǎn)與存儲(chǔ)策略研究[J]. 資源與產(chǎn)業(yè), 2011, 13(2): 63-68.

    [12] 侯玉梅, 潘登, 梁聰智. 碳排放交易下雙寡頭企業(yè)生產(chǎn)與減排研究[J]. 商業(yè)研究,2013, 55(1): 176-182.

    [13] 陳伯成,李英杰.“考慮排放許可交易的生產(chǎn)優(yōu)化”的補(bǔ)充研究[J]. 中國(guó)管理科學(xué), 2014, 22(9):141-148.

    [14] Islegen O, Reichelstein S J. Carbon capture by fossil fuel power plants: An economic analysis[J]. Management Science, 2011, 57(1): 21-39.

    [15] Jiang Yan, Klabjan D. Optimal emissions reduction investment under green house gas emissions regulations[R]. Working paper, Northwestern University, 2012.

    [16] 駱瑞玲,范體軍,夏海洋. 碳排放交易政策下供應(yīng)鏈碳減排技術(shù)投資的博弈分析[J]. 中國(guó)管理科學(xué), 2014, 22(11): 44-53.

    [17] 王明喜,鮑勤,湯鈴,等. 碳排放約束下的企業(yè)最優(yōu)減排投資行為[J]. 管理科學(xué)學(xué)報(bào), 2015, 18(6): 41-57.

    [18] 陳伯成,李英杰,閆學(xué)為. 補(bǔ)貼及懲罰模式下的碳排放許可交易生產(chǎn)優(yōu)化[J]. 中國(guó)管理科學(xué), 2014, 22(SI): 774-781.

    [19] 周艷菊,黃雨晴,陳曉紅,等. 促進(jìn)低碳產(chǎn)品需求的供應(yīng)鏈減排成本分擔(dān)模型[J]. 中國(guó)管理科學(xué), 2015, 23(7): 85-93.

    [20] Vogt-Schilb A, Meunier G, Hallegatte S. Should marginal abatement costs differ across sectors: The effect of low-carbon capital accumulation[A]. Working Paper, World Bank, 2013.

    [21] 魯奎, 楊昌輝, 戴道明. 能力受限批量問(wèn)題的啟發(fā)式算法與CPLEX仿真優(yōu)化[J]. 系統(tǒng)仿真學(xué)報(bào), 2008, 20(3): 6365-6371.

    [22] 魯奎, 楊昌輝, 戴道明. 運(yùn)輸能力受限與費(fèi)用時(shí)變批量問(wèn)題的拉格朗日松弛啟發(fā)式算法[J]. 系統(tǒng)工程理論與實(shí)踐, 2008,(10): 47-52.

    [23] Geoffrion A M. Generalized benders decomposition[J]. Journal of Optimization Theory and Application. 1972, 10(4):237-260.

    [24] Wagner H M, Whitin T M. Dynamic version of the economic lot size model[J]. Management Science, 1958, 5(1): 89-96.

    [25] Brahimi N, Dauzère-Pérès S, Najid N M, et al. Single item lot sizing problems[J]. European Journal of Operational Research, 2006, 168(1): 1-16.

    [26] Aggarwal A, Park J K. Improved algorithms for economic lot size problem[J]. Operations Research, 1993, 41(3):549-571.

    [27] Gade D, Kücükyavuz S. Formulations for dynamic lot sizing with service levels[J]. Naval Research Logistics, 2013, 60(2):87-101.

    Optimal Production and Emissions Reduction Investment Policies in a Dynamic Lot Sizing Model Under Cap-and-trade

    HUANG Di1,2, CHEN Jian2, ZHOU Hong3

    (1.School of Economics and Management, Beijing Jiaotong University, Beijing 100044, China; 2.School of Economics and Management, Tsinghua University, Beijing 100084, China; 3.School of Economics and Management, Beihang University, Beijing 100191, China)

    With the establishment and development of China’s emissions trading market, the reduction of carbon emissions intensity has become a medium and long term mandatory target in companies’ production operations management. A dynamic lot sizing model is proposed to investigate the optimal production, emissions trading and abatement investment decisions of a manufacturer under cap-and-trade. At the beginning of the planning horizon, the manufacturer determines the optimal abatement investment decisions, including the timing of investment in an abatement technology and the optimal abatement capacity. In each period, the manufacturer decides the optimal emission abatement size according to the production planning. It is assumed that the abatement technology can reduce both carbon emissions and energy consumption in the production process. Based on the generalized Benders decomposition approach, the dynamic model is analyzed and some properties of the optimal solution are derived. Numerical experiments are conducted to examine the effects of emission cap and emission allowance price on the manufacturer’s total cost, total emissions level, and the optimal abatement investment decisions. It is found that: (1) when the supply of emission allowances in the market is abundant, reducing the carbon cap or changing the initial allocation rule of free emission allowances will not reduce the manufacturer’s carbon emissions level; (2) carbon emission allowance price is the key factor to control the manufacturer’s carbon emissions level and abatement capacity; (3) with the increase of carbon price, the manufacturer will expand the abatement capacity even if it has already held enough allowances to cover all its emissions. Our research provides useful managerial insights for manufacturing firms to make investment in carbon emissions abatement under the cap-and-trade system.

    dynamic lot sizing; abatement investment; cap-and-trade; generalized benders decomposition

    1003-207(2016)04-0129-09

    10.16381/j.cnki.issn1003-207x.2016.04.015

    2014-06-23;

    2015-12-16

    國(guó)家自然科學(xué)基金資助項(xiàng)目(71232007, 71301082, 71471007, 71390334);中央高?;究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金資助項(xiàng)目(2014RC032)

    黃帝(1984-),男(漢族),河南南陽(yáng)人,北京交通大學(xué)經(jīng)濟(jì)管理學(xué)院講師,研究方向:低碳供應(yīng)鏈運(yùn)作管理,E-mail:huangd@sem.tsinghua.edu.cn.

    F253.4

    A

    猜你喜歡
    總成本配額生產(chǎn)商
    2020年中國(guó)棉花種植成本調(diào)查
    碳減排量及碳配額的區(qū)別
    新疆鋼鐵(2021年1期)2021-10-14 08:45:42
    魚(yú)粉:秘魯A季配額低于預(yù)期,內(nèi)外盤(pán)短期大幅上漲
    魚(yú)粉:秘魯A季配額公布,國(guó)內(nèi)外魚(yú)粉價(jià)格反彈
    數(shù)據(jù)驅(qū)動(dòng)下的庫(kù)存優(yōu)化模型研究
    線(xiàn)性盈虧平衡分析在TBM隧洞工程中的應(yīng)用
    生產(chǎn)商名錄
    生產(chǎn)商名錄
    生產(chǎn)商名錄
    關(guān)于煤化工生產(chǎn)企業(yè)成本管控的思考
    90打野战视频偷拍视频| 大香蕉久久成人网| 精品一区二区三区四区五区乱码| 久久九九热精品免费| 桃花免费在线播放| 中文字幕高清在线视频| 国产一级毛片在线| 另类亚洲欧美激情| 又大又爽又粗| 国产精品 欧美亚洲| 亚洲专区字幕在线| 亚洲国产毛片av蜜桃av| 丝袜美腿诱惑在线| 国产精品亚洲av一区麻豆| 免费在线观看日本一区| 国产欧美日韩一区二区三区在线| 国产精品二区激情视频| 亚洲天堂av无毛| 熟女少妇亚洲综合色aaa.| 男男h啪啪无遮挡| 另类精品久久| 日韩三级视频一区二区三区| 中文字幕人妻熟女乱码| 亚洲精品av麻豆狂野| 国产亚洲一区二区精品| www.999成人在线观看| 中国美女看黄片| 午夜成年电影在线免费观看| 欧美日韩精品网址| 久久久精品免费免费高清| 狠狠狠狠99中文字幕| 国产精品亚洲av一区麻豆| 女人爽到高潮嗷嗷叫在线视频| 丝瓜视频免费看黄片| www.av在线官网国产| 亚洲熟女精品中文字幕| av在线播放精品| 美女视频免费永久观看网站| 19禁男女啪啪无遮挡网站| 青草久久国产| 天堂俺去俺来也www色官网| av视频免费观看在线观看| 纵有疾风起免费观看全集完整版| 欧美少妇被猛烈插入视频| 欧美日韩av久久| 国产精品香港三级国产av潘金莲| 国产伦人伦偷精品视频| 亚洲国产中文字幕在线视频| 精品福利观看| 80岁老熟妇乱子伦牲交| 久久精品亚洲av国产电影网| 久久热在线av| 纯流量卡能插随身wifi吗| 黄色视频不卡| 免费在线观看黄色视频的| 日日摸夜夜添夜夜添小说| 丝袜在线中文字幕| 十八禁网站网址无遮挡| 三级毛片av免费| 少妇猛男粗大的猛烈进出视频| 999久久久国产精品视频| 女人久久www免费人成看片| 99香蕉大伊视频| 国产av国产精品国产| 动漫黄色视频在线观看| www.av在线官网国产| 久久精品久久久久久噜噜老黄| 久久国产亚洲av麻豆专区| 中文字幕制服av| 国产在线免费精品| 日韩大片免费观看网站| 国产精品久久久久成人av| 亚洲第一av免费看| 18在线观看网站| 欧美大码av| 秋霞在线观看毛片| 大香蕉久久网| 国产区一区二久久| 老司机亚洲免费影院| 精品一区二区三卡| 黄片大片在线免费观看| 一级毛片女人18水好多| 国产精品影院久久| 人成视频在线观看免费观看| 黑人欧美特级aaaaaa片| 国产深夜福利视频在线观看| 高清黄色对白视频在线免费看| 免费久久久久久久精品成人欧美视频| 国产亚洲一区二区精品| 成在线人永久免费视频| 一级片'在线观看视频| a级毛片黄视频| 两个人免费观看高清视频| 国产精品自产拍在线观看55亚洲 | 亚洲av片天天在线观看| 国产成人系列免费观看| 又大又爽又粗| 看免费av毛片| 亚洲欧美日韩另类电影网站| 午夜视频精品福利| 午夜视频精品福利| 欧美在线黄色| 一个人免费在线观看的高清视频 | 国产亚洲精品第一综合不卡| 国产免费av片在线观看野外av| 亚洲欧美精品自产自拍| 国产成人免费观看mmmm| 国产精品一区二区免费欧美 | 丝袜美腿诱惑在线| 国产成人一区二区三区免费视频网站| 欧美国产精品一级二级三级| 亚洲欧美清纯卡通| 老司机靠b影院| 黄片小视频在线播放| 亚洲av欧美aⅴ国产| 午夜老司机福利片| 午夜激情久久久久久久| 欧美另类一区| 亚洲成人手机| 一级a爱视频在线免费观看| 精品亚洲乱码少妇综合久久| 国产伦理片在线播放av一区| 久久99一区二区三区| 亚洲欧美成人综合另类久久久| 欧美少妇被猛烈插入视频| 国产在视频线精品| 亚洲熟女毛片儿| 丝袜喷水一区| 国产三级黄色录像| 美女主播在线视频| 国产精品欧美亚洲77777| 国产日韩欧美视频二区| 叶爱在线成人免费视频播放| 国产人伦9x9x在线观看| 午夜精品国产一区二区电影| 日本a在线网址| 欧美国产精品一级二级三级| 午夜精品国产一区二区电影| 日本a在线网址| 国产高清国产精品国产三级| 乱人伦中国视频| 男女下面插进去视频免费观看| 男女高潮啪啪啪动态图| 日韩欧美国产一区二区入口| 1024视频免费在线观看| 国产免费一区二区三区四区乱码| 精品人妻在线不人妻| 后天国语完整版免费观看| 亚洲中文日韩欧美视频| 国产主播在线观看一区二区| 精品国产一区二区久久| 狂野欧美激情性bbbbbb| 在线精品无人区一区二区三| 亚洲免费av在线视频| 在线 av 中文字幕| 19禁男女啪啪无遮挡网站| 日韩大片免费观看网站| 婷婷成人精品国产| 国产伦人伦偷精品视频| 啦啦啦中文免费视频观看日本| 国产亚洲av片在线观看秒播厂| 纵有疾风起免费观看全集完整版| 国产成人a∨麻豆精品| 日韩视频在线欧美| 亚洲五月婷婷丁香| 我要看黄色一级片免费的| 国产又色又爽无遮挡免| 99久久精品国产亚洲精品| 亚洲一区中文字幕在线| 亚洲av日韩精品久久久久久密| 老司机午夜福利在线观看视频 | 国产一级毛片在线| 男女下面插进去视频免费观看| 黄色a级毛片大全视频| 精品亚洲成a人片在线观看| 人成视频在线观看免费观看| 午夜日韩欧美国产| 精品人妻1区二区| 男女床上黄色一级片免费看| 亚洲国产欧美在线一区| 在线 av 中文字幕| 少妇被粗大的猛进出69影院| 久9热在线精品视频| 免费在线观看视频国产中文字幕亚洲 | 老汉色av国产亚洲站长工具| 中文字幕色久视频| 999久久久国产精品视频| 99国产精品一区二区蜜桃av | 欧美少妇被猛烈插入视频| 午夜精品国产一区二区电影| 免费高清在线观看视频在线观看| 欧美国产精品va在线观看不卡| 亚洲第一青青草原| 久久这里只有精品19| 九色亚洲精品在线播放| 成人亚洲精品一区在线观看| 两个人看的免费小视频| 国产真人三级小视频在线观看| tocl精华| 老司机影院成人| 十八禁网站网址无遮挡| 国产高清videossex| 久久久久久久国产电影| 桃红色精品国产亚洲av| 99国产精品一区二区三区| 亚洲男人天堂网一区| 免费不卡黄色视频| 97精品久久久久久久久久精品| 91成年电影在线观看| 法律面前人人平等表现在哪些方面 | 母亲3免费完整高清在线观看| 91大片在线观看| 日韩免费高清中文字幕av| 日韩欧美国产一区二区入口| 久久精品国产亚洲av高清一级| 久久久精品94久久精品| 精品少妇一区二区三区视频日本电影| kizo精华| 亚洲欧美成人综合另类久久久| 超碰成人久久| av欧美777| 最近最新免费中文字幕在线| 男女高潮啪啪啪动态图| 一本久久精品| 国产99久久九九免费精品| 亚洲三区欧美一区| 日韩欧美一区二区三区在线观看 | 伦理电影免费视频| 水蜜桃什么品种好| 80岁老熟妇乱子伦牲交| 欧美精品一区二区大全| 老司机午夜十八禁免费视频| 亚洲成人手机| av天堂在线播放| 久久人妻福利社区极品人妻图片| 女人精品久久久久毛片| 欧美国产精品va在线观看不卡| 国产av一区二区精品久久| 国内毛片毛片毛片毛片毛片| 麻豆乱淫一区二区| 男女无遮挡免费网站观看| 国产三级黄色录像| 国产极品粉嫩免费观看在线| 成人18禁高潮啪啪吃奶动态图| 久热爱精品视频在线9| 欧美激情 高清一区二区三区| 国产成人影院久久av| av在线播放精品| 亚洲国产成人一精品久久久| 在线精品无人区一区二区三| 香蕉国产在线看| 亚洲色图综合在线观看| 国产高清视频在线播放一区 | 午夜激情久久久久久久| 精品福利观看| 人人妻人人爽人人添夜夜欢视频| tocl精华| 黄色毛片三级朝国网站| 亚洲成人免费av在线播放| 丝袜人妻中文字幕| 9热在线视频观看99| 国产男女超爽视频在线观看| 三上悠亚av全集在线观看| 电影成人av| 人妻人人澡人人爽人人| 性高湖久久久久久久久免费观看| 国产高清videossex| 午夜影院在线不卡| 操出白浆在线播放| 亚洲三区欧美一区| 亚洲一区二区三区欧美精品| 国产老妇伦熟女老妇高清| 国产亚洲欧美在线一区二区| 久久热在线av| 1024香蕉在线观看| 久久久精品94久久精品| 国产在线观看jvid| 一进一出抽搐动态| 999久久久国产精品视频| 国产黄频视频在线观看| 亚洲一区中文字幕在线| 成年美女黄网站色视频大全免费| 丰满少妇做爰视频| 老司机影院毛片| 我的亚洲天堂| 午夜福利,免费看| 国产成人啪精品午夜网站| 男女午夜视频在线观看| 中文字幕人妻熟女乱码| 欧美 日韩 精品 国产| 一级片免费观看大全| 搡老熟女国产l中国老女人| 国产在线视频一区二区| 50天的宝宝边吃奶边哭怎么回事| 久久久欧美国产精品| 一级片'在线观看视频| 男女午夜视频在线观看| 国产成人免费观看mmmm| 99国产极品粉嫩在线观看| 免费观看a级毛片全部| 精品少妇一区二区三区视频日本电影| 一级片免费观看大全| 亚洲自偷自拍图片 自拍| 国产亚洲欧美在线一区二区| 国产成人免费无遮挡视频| 欧美老熟妇乱子伦牲交| 欧美另类一区| 国产精品香港三级国产av潘金莲| 老司机福利观看| 97在线人人人人妻| 日韩一卡2卡3卡4卡2021年| 亚洲国产成人一精品久久久| 欧美日韩一级在线毛片| 欧美 日韩 精品 国产| 国产精品偷伦视频观看了| 精品福利观看| 亚洲七黄色美女视频| 99九九在线精品视频| 国产亚洲欧美在线一区二区| 91精品伊人久久大香线蕉| 91精品三级在线观看| 法律面前人人平等表现在哪些方面 | 国产精品成人在线| 精品福利观看| 老鸭窝网址在线观看| 亚洲自偷自拍图片 自拍| 久久久国产精品麻豆| 亚洲欧美日韩高清在线视频 | 亚洲一区二区三区欧美精品| 亚洲精品成人av观看孕妇| 亚洲av日韩在线播放| 欧美乱码精品一区二区三区| videosex国产| 大片免费播放器 马上看| 美女高潮到喷水免费观看| 黄色怎么调成土黄色| 欧美激情 高清一区二区三区| 久久久久精品国产欧美久久久 | 日韩一卡2卡3卡4卡2021年| 在线亚洲精品国产二区图片欧美| 亚洲色图 男人天堂 中文字幕| 国产精品1区2区在线观看. | 国产精品久久久av美女十八| 国产精品秋霞免费鲁丝片| 18禁观看日本| 国产麻豆69| 丝瓜视频免费看黄片| 欧美成狂野欧美在线观看| 久久人人爽人人片av| 亚洲黑人精品在线| 国产精品熟女久久久久浪| 亚洲九九香蕉| 又黄又粗又硬又大视频| 高清黄色对白视频在线免费看| 午夜福利影视在线免费观看| 亚洲九九香蕉| 国产精品久久久久成人av| 久久女婷五月综合色啪小说| 99久久99久久久精品蜜桃| 欧美日韩中文字幕国产精品一区二区三区 | 大码成人一级视频| 国产真人三级小视频在线观看| 久久久久视频综合| 精品视频人人做人人爽| 欧美另类亚洲清纯唯美| 日韩欧美一区视频在线观看| 久久久精品94久久精品| 午夜福利,免费看| 久久人人爽av亚洲精品天堂| 亚洲av美国av| av国产精品久久久久影院| 欧美黑人欧美精品刺激| 亚洲精品av麻豆狂野| 美女高潮到喷水免费观看| 国产成人免费无遮挡视频| 性高湖久久久久久久久免费观看| 欧美乱码精品一区二区三区| 狠狠狠狠99中文字幕| 欧美国产精品va在线观看不卡| 国产亚洲精品久久久久5区| 两性午夜刺激爽爽歪歪视频在线观看 | 一级毛片精品| videosex国产| 亚洲欧美清纯卡通| 精品福利永久在线观看| 国产精品国产三级国产专区5o| 这个男人来自地球电影免费观看| 啦啦啦啦在线视频资源| 欧美日韩一级在线毛片| 十八禁高潮呻吟视频| 首页视频小说图片口味搜索| 成年动漫av网址| 免费av中文字幕在线| 在线永久观看黄色视频| 91老司机精品| 亚洲精品国产精品久久久不卡| 啦啦啦视频在线资源免费观看| 五月开心婷婷网| 亚洲av成人一区二区三| 成年美女黄网站色视频大全免费| 九色亚洲精品在线播放| 欧美成狂野欧美在线观看| 女人被躁到高潮嗷嗷叫费观| 免费不卡黄色视频| 国产淫语在线视频| 精品第一国产精品| www.av在线官网国产| 亚洲中文日韩欧美视频| 51午夜福利影视在线观看| 久久久久网色| 超色免费av| 欧美亚洲日本最大视频资源| 国产成人一区二区三区免费视频网站| av在线老鸭窝| 国产精品欧美亚洲77777| 日韩欧美一区二区三区在线观看 | 成年动漫av网址| 美女扒开内裤让男人捅视频| 他把我摸到了高潮在线观看 | 在线观看www视频免费| 久久香蕉激情| 国产亚洲av片在线观看秒播厂| 一本—道久久a久久精品蜜桃钙片| 巨乳人妻的诱惑在线观看| 美女扒开内裤让男人捅视频| kizo精华| 久久国产亚洲av麻豆专区| 精品国产一区二区三区久久久樱花| 精品久久久精品久久久| 欧美精品人与动牲交sv欧美| 动漫黄色视频在线观看| 欧美日韩一级在线毛片| 极品少妇高潮喷水抽搐| 日韩中文字幕视频在线看片| 久久天堂一区二区三区四区| 淫妇啪啪啪对白视频 | 欧美大码av| 在线精品无人区一区二区三| 王馨瑶露胸无遮挡在线观看| 99国产精品一区二区蜜桃av | a 毛片基地| 久久精品成人免费网站| 免费在线观看影片大全网站| 日本撒尿小便嘘嘘汇集6| 啦啦啦免费观看视频1| 伦理电影免费视频| 国产三级黄色录像| 国产精品熟女久久久久浪| 两性夫妻黄色片| 欧美人与性动交α欧美软件| 母亲3免费完整高清在线观看| 欧美成人午夜精品| 成人黄色视频免费在线看| 女人被躁到高潮嗷嗷叫费观| 一本色道久久久久久精品综合| 精品人妻在线不人妻| 精品第一国产精品| 午夜福利在线免费观看网站| 欧美精品啪啪一区二区三区 | 无限看片的www在线观看| 母亲3免费完整高清在线观看| 欧美成人午夜精品| 久久久精品国产亚洲av高清涩受| 男女无遮挡免费网站观看| 两个人看的免费小视频| 老司机影院毛片| 亚洲午夜精品一区,二区,三区| 国产黄频视频在线观看| av天堂在线播放| 亚洲国产中文字幕在线视频| 亚洲精品美女久久av网站| 天天躁日日躁夜夜躁夜夜| 母亲3免费完整高清在线观看| 欧美日韩黄片免| 免费不卡黄色视频| 在线永久观看黄色视频| 99国产精品99久久久久| 人成视频在线观看免费观看| 亚洲国产欧美日韩在线播放| 可以免费在线观看a视频的电影网站| 在线永久观看黄色视频| 午夜福利视频精品| 大型av网站在线播放| 久久久久久免费高清国产稀缺| 久久久精品国产亚洲av高清涩受| 色老头精品视频在线观看| 久久性视频一级片| 91成人精品电影| 伦理电影免费视频| 人人妻人人澡人人爽人人夜夜| 亚洲国产精品成人久久小说| 69精品国产乱码久久久| 国产精品久久久av美女十八| 久久久久久久久免费视频了| 欧美激情高清一区二区三区| 亚洲五月色婷婷综合| 国产成人av教育| 巨乳人妻的诱惑在线观看| 两个人看的免费小视频| 国内毛片毛片毛片毛片毛片| 亚洲中文字幕日韩| 亚洲国产av影院在线观看| 久久精品aⅴ一区二区三区四区| 精品国内亚洲2022精品成人 | 午夜免费观看性视频| 一区二区三区四区激情视频| 在线看a的网站| 妹子高潮喷水视频| 日韩大片免费观看网站| 免费在线观看完整版高清| 在线观看舔阴道视频| 最近最新免费中文字幕在线| www.自偷自拍.com| 多毛熟女@视频| 亚洲精品粉嫩美女一区| 欧美性长视频在线观看| 国产一区二区三区av在线| 一个人免费看片子| 黄色a级毛片大全视频| 国内毛片毛片毛片毛片毛片| 韩国高清视频一区二区三区| 女人精品久久久久毛片| 国产伦理片在线播放av一区| 人妻人人澡人人爽人人| 国产黄频视频在线观看| 少妇精品久久久久久久| 人人妻人人添人人爽欧美一区卜| 亚洲,欧美精品.| 91麻豆av在线| 欧美日本中文国产一区发布| 色老头精品视频在线观看| 嫁个100分男人电影在线观看| 人人妻,人人澡人人爽秒播| 久久久久国产一级毛片高清牌| av视频免费观看在线观看| 国产成人a∨麻豆精品| a级毛片在线看网站| 成人国产一区最新在线观看| 男女国产视频网站| 99国产综合亚洲精品| 又紧又爽又黄一区二区| 日韩一卡2卡3卡4卡2021年| 国产又色又爽无遮挡免| 女人爽到高潮嗷嗷叫在线视频| 桃红色精品国产亚洲av| 欧美av亚洲av综合av国产av| 人人澡人人妻人| 亚洲七黄色美女视频| 欧美日韩亚洲高清精品| 国产亚洲欧美精品永久| 久久天堂一区二区三区四区| 国产在视频线精品| 国产精品 国内视频| 久久精品国产a三级三级三级| 成人影院久久| 俄罗斯特黄特色一大片| 动漫黄色视频在线观看| 中文精品一卡2卡3卡4更新| 亚洲精品粉嫩美女一区| 777久久人妻少妇嫩草av网站| 久久久国产成人免费| 欧美在线黄色| 91精品伊人久久大香线蕉| 这个男人来自地球电影免费观看| 午夜免费成人在线视频| 久久午夜综合久久蜜桃| 国内毛片毛片毛片毛片毛片| 丝袜在线中文字幕| 国产精品欧美亚洲77777| 黑丝袜美女国产一区| 在线观看www视频免费| 中文精品一卡2卡3卡4更新| 黄色视频不卡| 又紧又爽又黄一区二区| 国产精品亚洲av一区麻豆| 一本一本久久a久久精品综合妖精| 亚洲自偷自拍图片 自拍| 99精品久久久久人妻精品| 日韩欧美一区视频在线观看| 丰满少妇做爰视频| 两性夫妻黄色片| 欧美在线黄色| 中文字幕av电影在线播放| 亚洲国产欧美一区二区综合| 乱人伦中国视频| svipshipincom国产片| 天天添夜夜摸| 日韩中文字幕欧美一区二区| 狂野欧美激情性xxxx| 日韩三级视频一区二区三区| 97精品久久久久久久久久精品| 成人国语在线视频| 亚洲精品粉嫩美女一区| 秋霞在线观看毛片| av又黄又爽大尺度在线免费看| 亚洲中文av在线| 亚洲三区欧美一区| 老司机福利观看| 亚洲人成77777在线视频| 999久久久精品免费观看国产| 最新在线观看一区二区三区| 91麻豆精品激情在线观看国产 | 天堂中文最新版在线下载| 亚洲精品国产av蜜桃| 久久午夜综合久久蜜桃| 男女免费视频国产| 亚洲 国产 在线| 久久综合国产亚洲精品| 精品福利永久在线观看| 男女午夜视频在线观看| 日日夜夜操网爽| 免费高清在线观看视频在线观看| 欧美97在线视频| 亚洲成人手机| 成人黄色视频免费在线看|