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

    重型柴油機(jī)顆粒物分布規(guī)律的試驗(yàn)研究

    2018-08-10 07:11:50周小波胡清華蘇萬(wàn)華
    關(guān)鍵詞:顆粒物

    周小波,胡清華,閆 峰,蘇萬(wàn)華

    ?

    重型柴油機(jī)顆粒物分布規(guī)律的試驗(yàn)研究

    周小波1,2,胡清華1※,閆 峰3,蘇萬(wàn)華2

    (1. 天津大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院,天津 300350; 2. 天津大學(xué)內(nèi)燃機(jī)燃燒學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072;3. 中國(guó)汽車(chē)技術(shù)研究中心有限公司,天津 300300)

    為了同時(shí)控制車(chē)用重型柴油機(jī)的顆粒物排放質(zhì)量和數(shù)目,該文對(duì)穩(wěn)態(tài)工況及瞬態(tài)工況下柴油機(jī)排放顆粒物分布規(guī)律進(jìn)行了試驗(yàn)研究。試驗(yàn)結(jié)果表明:在穩(wěn)態(tài)工況下,隨著負(fù)荷的增加或者轉(zhuǎn)速的提高,積聚態(tài)及核態(tài)顆粒物數(shù)目濃度、中位直徑(count median diameter,CMD)、表面積和質(zhì)量均呈現(xiàn)增大趨勢(shì),且峰值向大粒徑方向偏移。在本研究中,穩(wěn)態(tài)工況全工況總顆粒物數(shù)目濃度為1.5×106~4.5×106個(gè)/cm3,積聚態(tài)顆粒物數(shù)目濃度為2×106~1×107個(gè)/cm3,而核態(tài)顆粒物數(shù)目濃度為1×107~3×107個(gè)/cm3,在總顆粒物數(shù)目濃度中占比為65%~96%。全工況當(dāng)量比均小于0.7,在中低轉(zhuǎn)速,當(dāng)量比對(duì)顆粒物分布影響較為明顯,在高轉(zhuǎn)速尤其是大負(fù)荷條件下,當(dāng)量比的影響減弱。在瞬態(tài)工況下,顆粒物數(shù)目濃度出現(xiàn)了與顆粒物質(zhì)量類(lèi)似的排放尖峰現(xiàn)象,濃度峰值達(dá)到2×108~7×108個(gè)/cm3,比其對(duì)應(yīng)的穩(wěn)態(tài)工況出現(xiàn)的濃度峰值高出1個(gè)數(shù)量級(jí)。而且排放尖峰現(xiàn)象中,積聚態(tài)顆粒物數(shù)目濃度占主要部分,其峰值濃度比穩(wěn)態(tài)工況要高出2個(gè)數(shù)量級(jí),此時(shí)排放尖峰現(xiàn)象中核態(tài)顆粒物數(shù)目濃度也出現(xiàn)明顯增長(zhǎng)。排放尖峰現(xiàn)象對(duì)應(yīng)的粒徑主要集中在積聚態(tài)顆粒物的50~200 nm范圍和核態(tài)顆粒物的6~8 nm及20~50 nm范圍。這主要是因?yàn)楫?dāng)量比在瞬態(tài)過(guò)程出現(xiàn)了超過(guò)臨界當(dāng)量比0.8的情況。研究結(jié)果對(duì)重型柴油機(jī)顆粒物排放控制以及生成機(jī)理有重要指導(dǎo)意義,并可為排放后處理器的匹配計(jì)算提供數(shù)據(jù)支持。

    柴油機(jī);排氣;顆粒;穩(wěn)態(tài)工況;瞬態(tài)工況;分布規(guī)律;核態(tài)和積聚態(tài)顆粒

    0 引 言

    近年來(lái),中國(guó)各主要城市均出現(xiàn)了不同程度的霧霾現(xiàn)象。流行病學(xué)和毒理學(xué)研究表明[1-7],大氣中的超細(xì)顆粒物(顆粒物粒徑D<100 nm)附著大量有毒物質(zhì),通過(guò)呼吸道進(jìn)入人體后會(huì)對(duì)人體呼吸道和心血管等器官造成嚴(yán)重危害,有數(shù)據(jù)證明其與人類(lèi)發(fā)病率和死亡率有明顯關(guān)聯(lián)。中國(guó)約有30%的國(guó)土面積以及超過(guò)8億人口受到霧霾的影響[8]。同時(shí),顆粒物的輻射效應(yīng)僅次于二氧化碳(CO2),是造成溫室效應(yīng)的重要原因之一[9]。柴油機(jī)排放的微粒直徑一般為3 nm~1m,其數(shù)目濃度高達(dá)10×107~10×109個(gè)/cm3,是城市大氣中超細(xì)顆粒的重要來(lái)源。通常,排放微粒按粒徑分為核態(tài)顆粒物(5 nm<D<50 nm)和積聚態(tài)顆粒物(50 nm<D<0.15m),其中積聚態(tài)顆粒物主要由不完全燃燒產(chǎn)生的一次碳粒(D為2 nm左右)發(fā)生團(tuán)聚作用并凝結(jié)有部分HC和硫酸鹽等揮發(fā)和半揮發(fā)組分形成;核態(tài)顆粒物的成因比較復(fù)雜,通常認(rèn)為核態(tài)顆粒物是由燃燒室內(nèi)形成的一次碳粒以及H2SO4或HC等氣態(tài)前體物成核形成的二次顆粒物,并且其數(shù)濃度和粒徑分布特征受稀釋系統(tǒng)參數(shù)(如顆粒在通道內(nèi)的停留時(shí)間、稀釋比例、稀釋空氣的溫度和濕度等)影響較大[10-11]。

    從2000年開(kāi)始,隨著DPF的推廣使用,研究者發(fā)現(xiàn)稱(chēng)重式的PM檢測(cè)方法已經(jīng)不足以準(zhǔn)確反映內(nèi)燃機(jī)顆粒排放水平。因此2001年歐盟組織排放顆粒物測(cè)量計(jì)劃PMP(particle measurement program)開(kāi)始研究更靈敏和準(zhǔn)確的內(nèi)燃機(jī)排放顆粒物檢測(cè)手段。Giechaskiel等[12-13]于2007年和2010年分別提交輕型車(chē)和重型車(chē)顆粒數(shù)檢測(cè)方法的最終報(bào)告并被PMP采用,這種去除顆粒物中揮發(fā)成分后進(jìn)行顆粒物數(shù)量統(tǒng)計(jì)的檢測(cè)顆粒物數(shù)量排放的新方法得到了廣泛的認(rèn)可并被列入法規(guī)考察范圍。自此,顆粒物排放數(shù)目和粒徑分布成為研究?jī)?nèi)燃機(jī)顆粒物排放特性的新熱點(diǎn)。

    目前,國(guó)內(nèi)外針對(duì)柴油機(jī)顆粒物數(shù)目和粒徑分布的研究主要集中在穩(wěn)態(tài)工況下[14-19]。1998年Kittelson等[20]提出典型柴油機(jī)穩(wěn)態(tài)工況顆粒物粒徑分布、數(shù)目、表面積和質(zhì)量。顆粒物數(shù)量和質(zhì)量濃度基本呈現(xiàn)對(duì)數(shù)正態(tài)分布規(guī)律,核態(tài)顆粒物占總數(shù)量的90%以上,積聚態(tài)顆粒物占總質(zhì)量的比例最大,粗態(tài)顆粒物數(shù)量較少,其質(zhì)量約占總質(zhì)量的5~20%。但是近十幾年來(lái)隨著燃燒理論的發(fā)展以及后處理器的使用,粗態(tài)顆粒物明顯減少,研究者把研究目標(biāo)主要集中在核態(tài)和積聚態(tài)顆粒物上。Kazakov等[21]提出碳煙顆粒碰撞理論,認(rèn)為燃燒溫度是決定顆粒物尺寸和分布的主要因素。Desantes等[22]的試驗(yàn)結(jié)果顯示,積聚態(tài)顆粒物與核態(tài)顆粒物之間是此消彼長(zhǎng)的關(guān)系。即積聚態(tài)顆粒物如同“海綿體”一樣來(lái)吸收HC和硫酸鹽等揮發(fā)和半揮發(fā)組分。當(dāng)積聚態(tài)顆粒大幅度增加后,用以促進(jìn)核態(tài)顆粒物生長(zhǎng)的組分減少,導(dǎo)致其濃度較低。Shi等[23-24]在穩(wěn)態(tài)工況下系統(tǒng)地研究了顆粒物的物理、化學(xué)和形態(tài)學(xué)特征。

    在瞬態(tài)工況下,Tan等[25]在輕型柴油機(jī)上研究了生物柴油和普通柴油的不同摻混比例在穩(wěn)態(tài)和瞬態(tài)工況下對(duì)顆粒物數(shù)目排放的影響。試驗(yàn)結(jié)果表明在穩(wěn)態(tài)工況下,隨著生物柴油比例的提高,積聚態(tài)顆粒物數(shù)目排放減少,核態(tài)顆粒物數(shù)目增加。在瞬態(tài)工況下,扭矩瞬變率增大,總顆粒物數(shù)目明顯提高;使用生物柴油時(shí),顆粒物數(shù)目排放明顯高于普通柴油。Wang等[26]研究了在美國(guó)瞬態(tài)測(cè)試循環(huán)中轉(zhuǎn)速變化率及扭矩變化率對(duì)柴油機(jī)顆粒物數(shù)目的影響。研究表明在從怠速開(kāi)始的急加速工況中顆粒物數(shù)目會(huì)大幅增加。在瞬態(tài)工況下,柴油機(jī)顆粒物數(shù)目及粒徑分布特性的研究還比較少,因此,對(duì)比分析穩(wěn)態(tài)工況和瞬態(tài)工況下柴油機(jī)排放顆粒物的生成規(guī)律和影響因素具有重大的現(xiàn)實(shí)意義和科研價(jià)值。

    1 試驗(yàn)裝置與方法

    本文試驗(yàn)臺(tái)架是無(wú)后處理器、滿足歐IV法規(guī)的試驗(yàn)樣機(jī),并通過(guò)了中國(guó)汽車(chē)技術(shù)研究中心的認(rèn)證,主體部分包括某12 L重型柴油機(jī)和電力測(cè)功機(jī)[27],如圖1所示。

    1.數(shù)據(jù)采集系統(tǒng) 2. 柴油機(jī) 3.IVCT系統(tǒng) 4.高壓共軌燃油噴射系統(tǒng) 5.電力測(cè)功系統(tǒng) 6.中冷器 7.EGR系統(tǒng) 8.兩級(jí)渦輪增壓系統(tǒng) 9.VGT系統(tǒng) 10.智能控制器 11.瞬時(shí)流量計(jì)

    發(fā)動(dòng)機(jī)進(jìn)氣系統(tǒng)部分包括高壓級(jí)為可變幾何截面渦輪(variable geometry turbocharger ,VGT)的兩級(jí)增壓系統(tǒng)、高壓廢氣再循環(huán) (exhaust gas recirculation,EGR)系統(tǒng)和課題組自主研發(fā)的進(jìn)氣門(mén)晚關(guān)系統(tǒng)(IVCT系統(tǒng))[28]。測(cè)試儀器包括高響應(yīng)的進(jìn)排放壓力傳感器、進(jìn)氣流量計(jì)、油耗儀、排放采集儀等。試驗(yàn)中使用Cambustion公司的DMS500快速顆粒光譜儀研究分析柴油機(jī)瞬態(tài)工況顆粒物粒徑分布規(guī)律,其測(cè)量范圍是5~1 000 nm,響應(yīng)時(shí)間是100 ms。試驗(yàn)方法參照GB17691中ETC試驗(yàn)規(guī)范執(zhí)行。

    2 結(jié)果與分析

    2.1 穩(wěn)態(tài)工況下顆粒物粒徑分布規(guī)律

    為了研究顆粒物粒徑分布規(guī)律,穩(wěn)態(tài)試驗(yàn)選取歐洲穩(wěn)態(tài)測(cè)試循環(huán)(european steady-state cycle,ESC)十三工況作為主要研究工況,并在穩(wěn)態(tài)工況平面內(nèi)另選了54個(gè)工況點(diǎn)以完成全工況顆粒物排放數(shù)目脈譜。54個(gè)工況點(diǎn)按照轉(zhuǎn)速和扭矩均勻選取,兼顧了轉(zhuǎn)速分布和扭矩分布。

    中位直徑(count median diameter,CMD)是累積百分比為50%時(shí)所對(duì)應(yīng)的粒子直徑,直接關(guān)系到粒譜特征的判斷和表征,是表達(dá)顆粒物粒徑的重要參數(shù)之一。研究表明柴油機(jī)排放微粒的粒徑分布近似于對(duì)數(shù)正態(tài)分布,CMD用式(1)計(jì)算[29]:

    式中d是各分級(jí)的切割直徑,N是分級(jí)對(duì)應(yīng)的粒子直徑。

    2.1.1 負(fù)荷對(duì)穩(wěn)態(tài)工況顆粒物粒徑分布的影響

    圖2為穩(wěn)態(tài)工況下該發(fā)動(dòng)機(jī)為轉(zhuǎn)速1 600 r/min時(shí)負(fù)荷變化對(duì)顆粒物粒徑分布、數(shù)目濃度、表面積、質(zhì)量和CMD的影響。圖2a為負(fù)荷對(duì)積聚態(tài)顆粒物數(shù)目濃度和粒徑分布的影響,隨著負(fù)荷的增大,積聚態(tài)顆粒物數(shù)目濃度和粒徑分布的峰值隨之增大,范圍處于4×106~8×106個(gè)/cm3之間,此時(shí)CMD也隨之增大,范圍在40~50 nm之間。由碳煙顆粒碰撞理論[21]可知,隨著負(fù)荷的增大,缸內(nèi)燃燒溫度增加,碳煙顆粒之間碰撞的頻率提高,顆粒粒徑增大,CMD增大。圖2b為負(fù)荷對(duì)核態(tài)顆粒物數(shù)目濃度和粒徑分布的影響。隨著負(fù)荷的增大,核態(tài)顆粒物數(shù)目濃度和粒徑分布的峰值總體呈現(xiàn)增大趨勢(shì),范圍處于1×107~2×107個(gè)/cm3之間,同時(shí)CMD也隨之增大,范圍在5~9 nm之間,如圖2c所示。核態(tài)顆粒物的數(shù)目濃度相比積聚態(tài)顆粒物要高1個(gè)量級(jí),如圖2d所示,核態(tài)顆粒數(shù)目濃度在總顆粒數(shù)濃度中占比59%~80%,排放顆粒物主要為核態(tài)顆粒物。隨著負(fù)荷的增加,積聚態(tài)顆粒物數(shù)目濃度呈現(xiàn)增大趨勢(shì),其對(duì)應(yīng)的CMD也呈現(xiàn)增大趨勢(shì);核態(tài)顆粒物數(shù)目濃度總體呈現(xiàn)增大趨勢(shì),其對(duì)應(yīng)的CMD也呈現(xiàn)增大趨勢(shì)。

    顆粒物毒性與顆粒物表面積密切相關(guān)[30-31]。假設(shè)顆粒物均為球形規(guī)則體,根據(jù)顆粒物數(shù)目濃度和粒徑分布,可以得到顆粒物表面積分布,如圖2e與圖2f所示。隨著負(fù)荷的增加,積聚態(tài)顆粒物表面積增加,對(duì)應(yīng)峰值在5×1010~18×1010nm2/cm3之間,核態(tài)顆粒物表面積隨負(fù)荷增加整體呈現(xiàn)增大趨勢(shì),其對(duì)應(yīng)峰值在1×109~7×109nm2/cm3之間,比積聚態(tài)顆粒物表面積峰值低1~2個(gè)量級(jí)。而且,相比顆粒物數(shù)目濃度,顆粒物表面積和粒徑分布對(duì)應(yīng)的峰值無(wú)論是積聚態(tài)還是核態(tài)均往大粒徑方向偏離。由于顆粒物密度無(wú)法測(cè)量且形狀不規(guī)則,此處假設(shè)顆粒物密度均勻,核態(tài)與積聚態(tài)密度一致,且均為球形規(guī)則體。圖2g與圖2h為負(fù)荷對(duì)積聚態(tài)和核態(tài)顆粒物質(zhì)量和粒徑分布的影響。顆粒物質(zhì)量和粒徑的分布規(guī)律與顆粒物數(shù)目濃度和粒徑的分布規(guī)律有明顯區(qū)別,隨著負(fù)荷的增加,積聚態(tài)顆粒物質(zhì)量呈明顯增加趨勢(shì),對(duì)應(yīng)峰值在1×104~4×104ng/m3之間,核態(tài)顆粒物質(zhì)量與粒徑分布隨負(fù)荷增加整體呈現(xiàn)增大趨勢(shì),其對(duì)應(yīng)峰值在20~80 ng/m3之間,比積聚態(tài)顆粒物質(zhì)量低2~3個(gè)數(shù)量級(jí)。相比顆粒物數(shù)目濃度和表面積分布,顆粒物質(zhì)量和粒徑分布的對(duì)應(yīng)峰值無(wú)論是積聚態(tài)還是核態(tài)均往大粒徑方向偏離。

    圖2 穩(wěn)態(tài)工況下負(fù)荷對(duì)顆粒物分布的影響

    2.1.2 轉(zhuǎn)速對(duì)穩(wěn)態(tài)工況下顆粒物粒徑分布的影響

    圖3為穩(wěn)態(tài)工況下50%負(fù)荷時(shí)轉(zhuǎn)速變化(1 300、1 600和1 900 r/min)對(duì)顆粒物粒徑與數(shù)目濃度分布、表面積、質(zhì)量和CMD的影響。圖3a為轉(zhuǎn)速對(duì)積聚態(tài)顆粒物數(shù)目濃度和粒徑分布的影響,隨著轉(zhuǎn)速的提高,積聚態(tài)顆粒物數(shù)目濃度和粒徑分布的峰值隨之增大,范圍處于5×106~8×106個(gè)/cm3之間,此時(shí)CMD也隨之增大,范圍在37~47 nm之間。圖3b為轉(zhuǎn)速對(duì)核態(tài)顆粒物數(shù)目濃度和粒徑分布的影響,隨著轉(zhuǎn)速的提高,核態(tài)顆粒物數(shù)目濃度和粒徑分布的峰值隨之增大,范圍處于1×107~2×107個(gè)/cm3之間,同時(shí)CMD也隨之增大,范圍在6~8.5 nm之間,如圖3c所示。核態(tài)顆粒物的數(shù)目濃度相比積聚態(tài)顆粒物要高1個(gè)量級(jí),如圖3d所示,在總顆粒物數(shù)目中,核態(tài)顆粒物占比在64%~75%之間。這主要是由于隨著轉(zhuǎn)速的上升,燃燒過(guò)程中的實(shí)際混合時(shí)間減少,在缸內(nèi)局部缺氧條件下,更易導(dǎo)致碳煙顆粒生成量的增加,這些顆粒以積聚態(tài)為主;而核態(tài)顆粒物其成核作用受排放溫度影響較大,隨著轉(zhuǎn)速的提高,排放溫度從670 K增大到780 K,核態(tài)顆粒物數(shù)目濃度隨著轉(zhuǎn)速的提高而增大,且CMD也隨之增大。因此隨著轉(zhuǎn)速的提高,積聚態(tài)顆粒物數(shù)目濃度上升,其對(duì)應(yīng)的CMD也增大;核態(tài)顆粒物數(shù)目濃度呈現(xiàn)增大趨勢(shì),其對(duì)應(yīng)的CMD也呈現(xiàn)增大趨勢(shì),排放顆粒物主要為核態(tài)顆粒物,且隨轉(zhuǎn)速提高而增加。

    顆粒物表面積分布如圖3e與圖3f所示。顆粒物表面積和粒徑分布規(guī)律與顆粒物數(shù)目和粒徑分布規(guī)律一致,隨著轉(zhuǎn)速的提高,積聚態(tài)顆粒物表面積明顯增加,對(duì)應(yīng)峰值在6×1010~16×1010nm2/cm3之間,跨度較大;核態(tài)顆粒物表面積隨轉(zhuǎn)速增加而增大,其對(duì)應(yīng)峰值在2×109~6×109nm2/cm3之間,比積聚態(tài)顆粒物表面積峰值低1~2個(gè)量級(jí)。而且,相比顆粒物數(shù)目濃度,顆粒物表面積和粒徑分布對(duì)應(yīng)峰值無(wú)論是積聚態(tài)還是核態(tài)均往大粒徑方向偏離。

    圖3g與圖3h為轉(zhuǎn)速對(duì)積聚態(tài)和核態(tài)顆粒物質(zhì)量與粒徑分布的影響。隨著轉(zhuǎn)速的增加,積聚態(tài)顆粒物質(zhì)量呈明顯增加趨勢(shì),對(duì)應(yīng)峰值在7×106~30×106ng/cm3之間,核態(tài)顆粒物質(zhì)量隨轉(zhuǎn)速增加而增大,其對(duì)應(yīng)峰值在2×104~9×104ng/cm3之間,比積聚態(tài)顆粒物質(zhì)量峰值低2個(gè)量級(jí)以上。相比顆粒物數(shù)目濃度和表面積分布,顆粒物質(zhì)量分布對(duì)應(yīng)峰值無(wú)論是積聚態(tài)還是核態(tài)均往大粒徑方向偏離。因此,大粒徑顆粒物是顆粒物質(zhì)量的主要來(lái)源,且積聚態(tài)顆粒物遠(yuǎn)大于核態(tài)顆粒物。

    圖3 穩(wěn)態(tài)工況下轉(zhuǎn)速對(duì)顆粒物分布的影響

    2.1.3 穩(wěn)態(tài)工況顆粒物分布

    為完成全工況顆粒物排放數(shù)目脈譜,穩(wěn)態(tài)試驗(yàn)包括十三工況及其他54個(gè)工況。圖4為穩(wěn)態(tài)工況下全工況平面內(nèi)總顆粒物數(shù)目濃度(圖4a)、核態(tài)顆粒物數(shù)目濃度(圖4b)、積聚態(tài)顆粒物數(shù)目濃度(圖4c)及當(dāng)量比分布規(guī)律(圖4d)。其中總顆粒物數(shù)目濃度在1.5×107~4.5×107個(gè)/cm3,積聚態(tài)顆粒物數(shù)目濃度為2×106~1×107個(gè)/cm3,核態(tài)顆粒物數(shù)目濃度為1×107~3×107個(gè)/cm3,在總顆粒物數(shù)目濃度中占比65%~96%。從整個(gè)穩(wěn)態(tài)工況平面可以看出,全工況當(dāng)量比均小于0.7,在中低轉(zhuǎn)速,當(dāng)量比對(duì)顆粒物分布影響較為明顯,而在高轉(zhuǎn)速尤其是大負(fù)荷條件下,當(dāng)量比的影響減弱。

    圖4 穩(wěn)態(tài)工況下顆粒物數(shù)目濃度及對(duì)應(yīng)的當(dāng)量比分布

    2.2 瞬態(tài)測(cè)試循環(huán)的顆粒物分布

    本文按照歐洲瞬態(tài)測(cè)試循環(huán)(ETC)完成,瞬態(tài)控制策略參考文獻(xiàn)[27]。圖5為在ETC瞬態(tài)測(cè)試循環(huán)中轉(zhuǎn)速和扭矩、當(dāng)量比變化趨勢(shì)和顆粒物分布規(guī)律,包括總顆粒物數(shù)目濃度、核態(tài)顆粒物數(shù)目濃度和積聚態(tài)顆粒物數(shù)目濃度。從顆粒物數(shù)目濃度來(lái)看,瞬態(tài)工況下出現(xiàn)了與顆粒物質(zhì)量類(lèi)似的排放尖峰現(xiàn)象,數(shù)目濃度峰值甚至達(dá)到了2×108~7×108個(gè)/cm3,比其對(duì)應(yīng)的穩(wěn)態(tài)工況出現(xiàn)的濃度峰值高出1個(gè)數(shù)量級(jí)。而且排放尖峰現(xiàn)象中,積聚態(tài)顆粒物占主要部分,這與穩(wěn)態(tài)工況有很大不同,其數(shù)目濃度峰值比穩(wěn)態(tài)工況要高出2個(gè)數(shù)量級(jí),從穩(wěn)態(tài)工況的2×106~1×107個(gè)/cm3跨越式增長(zhǎng)到1×108~8×108個(gè)/cm3;此時(shí)排放尖峰現(xiàn)象中核態(tài)顆粒物數(shù)目濃度也出現(xiàn)明顯的增長(zhǎng),其數(shù)目濃度峰值從穩(wěn)態(tài)工況的1×107~3×107個(gè)/cm3增長(zhǎng)到3×107~3×108個(gè)/cm3。排放尖峰現(xiàn)象主要是因?yàn)楫?dāng)量比在瞬態(tài)過(guò)程出現(xiàn)了對(duì)應(yīng)的峰值,穩(wěn)態(tài)工況下當(dāng)量比沒(méi)有超過(guò)0.7,但是排放尖峰對(duì)應(yīng)的工況出現(xiàn)了超過(guò)臨界當(dāng)量比0.8的情況[27]。

    圖5 瞬態(tài)工況下顆粒物分布

    在城市道路工況中,排放尖峰出現(xiàn)12次,總顆粒物數(shù)目濃度峰值為2×108~8×108個(gè)/cm3,積聚態(tài)顆粒物數(shù)目濃度峰值為1×108~5×108個(gè)/cm3,核聚態(tài)顆粒物數(shù)目濃度峰值為3×107~3×108個(gè)/cm3;在鄉(xiāng)村道路工況中,排放尖峰出現(xiàn)8次,總顆粒物數(shù)目濃度峰值明顯低于城市道路工況,為1×108~3×108個(gè)/cm3,積聚態(tài)顆粒物數(shù)目濃度峰值為1×108~2×108個(gè)/cm3,核聚態(tài)顆粒物數(shù)目濃度峰值為3×107~1×108個(gè)/cm3;在高速道路工況中,排放尖峰沒(méi)有出現(xiàn),總顆粒物數(shù)目濃度、積聚態(tài)顆粒物數(shù)目濃度、核聚態(tài)顆粒物數(shù)目濃度與穩(wěn)態(tài)工況顆粒物粒徑分布規(guī)律相當(dāng)。圖6為在ETC瞬態(tài)測(cè)試循環(huán)中顆粒物粒徑分布的頻譜圖,更準(zhǔn)確地反應(yīng)了各個(gè)顆粒物粒徑對(duì)應(yīng)的顆粒物數(shù)目濃度。

    圖6 瞬態(tài)工況顆粒物粒徑分布頻譜圖

    為了方便分析,把城市道路工況、鄉(xiāng)村道路工況和高速道路工況的前50 s對(duì)應(yīng)的顆粒物粒徑分布頻譜圖放大,與整個(gè)測(cè)試循環(huán)并列展示。

    從整個(gè)測(cè)試循環(huán)可以發(fā)現(xiàn),在顆粒物粒徑大于200 nm時(shí),3個(gè)道路工況并沒(méi)有明顯區(qū)別。在顆粒物粒徑50~200 nm內(nèi),數(shù)目濃度峰值達(dá)到1×108個(gè)/cm3以上的情況(即圖中白色區(qū)域)主要集中在城市道路工況(12次)和鄉(xiāng)村道路工況(8次),高速道路工況沒(méi)有,這也就是前面討論的瞬態(tài)排放尖峰中的積聚態(tài)顆粒物。城市道路工況對(duì)應(yīng)的數(shù)目濃度峰值達(dá)到1×108個(gè)/cm3以上的情況(即圖中白色區(qū)域)持續(xù)時(shí)間長(zhǎng),大多超過(guò)10 s,而鄉(xiāng)村道路工況對(duì)應(yīng)的數(shù)目濃度峰值達(dá)到1×108個(gè)/cm3以上的情況(即圖中白色區(qū)域)持續(xù)時(shí)間較短,大多在5 s以內(nèi)。顆粒物粒徑在20~50 nm時(shí),城市道路工況和鄉(xiāng)村道路工況出現(xiàn)排放尖峰的時(shí)刻,該粒徑段也出現(xiàn)了數(shù)目濃度峰值達(dá)到1×108個(gè)/cm3以上的情況(即圖中白色區(qū)域)。顆粒物粒徑小于20 nm時(shí),城市道路工況中數(shù)目濃度峰值達(dá)到1×108個(gè)/cm3以上的情況(即圖中白色區(qū)域)很少,而鄉(xiāng)村道路工況和高速道路工況較多,且數(shù)目濃度峰值均出現(xiàn)在顆粒物粒徑為6~8 nm時(shí)。因此,在瞬態(tài)工況下排放尖峰現(xiàn)象對(duì)應(yīng)的粒徑主要集中在積聚態(tài)顆粒物的50~200 nm范圍和核態(tài)顆粒物的6~8 nm及20~50 nm范圍。

    3 結(jié) 論

    本文對(duì)穩(wěn)態(tài)工況及瞬態(tài)工況下柴油機(jī)排放顆粒物粒徑分布規(guī)律分別進(jìn)行了分析,結(jié)果如下:

    1)在穩(wěn)態(tài)工況下,隨著負(fù)荷的增加,積聚態(tài)顆粒物數(shù)目濃度呈現(xiàn)上升趨勢(shì),其對(duì)應(yīng)的CMD也呈現(xiàn)增大趨勢(shì),積聚態(tài)顆粒物表面積和質(zhì)量均呈現(xiàn)增大趨勢(shì)且峰值向大粒徑方向偏移;核態(tài)顆粒物數(shù)目濃度總體呈現(xiàn)增大趨勢(shì),其對(duì)應(yīng)的CMD也呈現(xiàn)增大趨勢(shì),核態(tài)顆粒物表面積和質(zhì)量均呈現(xiàn)增大趨勢(shì)且峰值向大粒徑方向偏移。

    2)隨著轉(zhuǎn)速的提高,積聚態(tài)顆粒物數(shù)目濃度上升,其對(duì)應(yīng)的CMD呈現(xiàn)增大趨勢(shì),積聚態(tài)顆粒物表面積和質(zhì)量均呈現(xiàn)增大趨勢(shì)且峰值向大粒徑方向偏移;核態(tài)顆粒物數(shù)目濃度呈現(xiàn)增大趨勢(shì),其對(duì)應(yīng)的CMD呈現(xiàn)增大趨勢(shì),核態(tài)顆粒物表面積和質(zhì)量均呈現(xiàn)增大趨勢(shì)且峰值向大粒徑方向偏移。

    3)穩(wěn)態(tài)工況全工況平面內(nèi)總顆粒物數(shù)目濃度為1.5×107~4.5×107個(gè)/cm3,積聚態(tài)顆粒物數(shù)目濃度為2×106~1×107個(gè)/cm3,核態(tài)顆粒物數(shù)目濃度為1×107~3×107個(gè)/cm3,且在總顆粒物數(shù)目濃度中占比65%~96%之間,因此在穩(wěn)態(tài)工況下,核態(tài)顆粒物數(shù)目是總顆粒物數(shù)目的主要來(lái)源。全工況當(dāng)量比均小于0.7,在中低轉(zhuǎn)速當(dāng)量比對(duì)顆粒物分布影響較為明顯,而在高轉(zhuǎn)速尤其是大負(fù)荷條件下,當(dāng)量比的影響減弱。

    4)在瞬態(tài)工況下,顆粒物數(shù)目濃度出現(xiàn)了與顆粒物質(zhì)量類(lèi)似的排放尖峰現(xiàn)象,濃度峰值達(dá)到2×108~7×108個(gè)/cm3,比其對(duì)應(yīng)的穩(wěn)態(tài)工況出現(xiàn)的濃度峰值高出1個(gè)數(shù)量級(jí)。而且排放尖峰現(xiàn)象中,積聚態(tài)顆粒物數(shù)目濃度占主要部分,其峰值濃度比穩(wěn)態(tài)工況要高出2個(gè)數(shù)量級(jí),從穩(wěn)態(tài)工況的2×106~1×107個(gè)/cm3跨越式增長(zhǎng)到1×108~8×108個(gè)/cm3;排放尖峰現(xiàn)象中核態(tài)顆粒物數(shù)目濃度也出現(xiàn)明顯的增長(zhǎng),其峰值濃度從穩(wěn)態(tài)工況的1×107~3×107個(gè)/cm3增長(zhǎng)到3×107~3×108個(gè)/cm3。排放尖峰現(xiàn)象對(duì)應(yīng)的粒徑主要集中在積聚態(tài)顆粒物的50~200 nm范圍和核態(tài)顆粒物的6~8 nm及20~50 nm范圍。

    [1] Lamarche B, Moorjani S, Cantin B, et al. Small, dense low-density lipoprotein particles as a predictor of the risk of ischemic heart disease in men: prospective results from the Qué bec Cardiovascular Study[J]. Circulation, 1997, 95(1): 69-75.

    [2] Wasko M K, Goodman S B. Emperor’s new clothes: Is particle disease really infected particle disease?[J]. Journal of Orthopaedic Research, 2016, 34(9): 1497-1504.

    [3] 馮仲科,毛海穎,李虹. 環(huán)首都圈植被分布與可吸入顆粒物的空間相關(guān)性[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(1):220-227.

    Feng Zhongke, Mao Haiying, Li Hong. Spatial correlation between vegetation distribution and respirable particulate matter around capital region[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(1): 220-227. (in Chinese with English abstract)

    [4] Suzuki S, Nishikawa A, Kuwana M, et al. Inflammatory myopathy with anti-signal recognition particle antibodies: case series of 100 patients[J]. Orphanet journal of rare diseases, 2015, 10(1): 61.

    [5] 張延峰,宋崇林,成存玉,等. 車(chē)用柴油機(jī)排放顆粒物中有機(jī)組分和無(wú)機(jī)組分的分析[J]. 燃燒科學(xué)與技術(shù),2004,10(3):197-201.

    Zhang Yanfeng, Song Zonglin, Cheng Cuiyu, et al. Analysis of the organic and inorganic components in the emission particulates from diesel engine[J]. Journal of Combustion Science & Technology, 2004, 10(3): 197-201. (in Chinese with English abstract)

    [6] Stewart Goldman M D, Turner C D, Jr D P, et al. Katsouyanni K, Touloumi G, Spix C. Short-term effects of ambient sulphur dioxide and particulate matter on mortality in 12 European cities: results from times series data from the APHEA project[J]. British Medical Journal, 1997, 314(7095): 1658-1653.

    [7] 李林,周啟星. 不同粒徑大氣顆粒物與死亡終點(diǎn)關(guān)系的流行病學(xué)研究回顧[J]. 環(huán)境與職業(yè)醫(yī)學(xué),2015,32(2):168-174.

    Li Lin, Zhou Qixing. Size-fractionated atmospheric particulate matters and mortality: An epidemiological review[J]. Journal of Environmental & Occupational Medicine, 2015, 32(2): 168-174. (in Chinese with English abstract)

    [8] 郝斌. 不同燃料對(duì)柴油機(jī)排放顆粒物的影響研究[D]. 天津:天津大學(xué),2014.

    Hao Bin. Effect of Fuel Identity on the Exhaust Particles from Diesel Engine[D]. Tianjin: Tianjin University, 2014. (in Chinese with English abstract)

    [9] 梅德清,趙翔,王書(shū)龍,等. 柴油機(jī)排放顆粒物的熱重特性分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(16):50-56.

    Mei Deqing, Zhao Xiang, Wang Shulong, et al. Thermogravimetric characteristics analysis of particulate matter of emission of divided diesel[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(16): 50-56. (in Chinese with English abstract)

    [10] 王強(qiáng). 乙烯預(yù)混火焰中碳煙微觀形貌、納觀結(jié)構(gòu)及氧化特性研究[D]. 天津:天津大學(xué),2015.

    Wang Qiang. Studies on the Morphology, Nanostructure and Oxidation Reactivity of Soot in a C2H4/O2/Ar Premixed Flame[D]. Tianjin: Tianjin University, 2015. (in Chinese with English abstract)

    [11] 于文斌. 基于混合與化學(xué)反應(yīng)協(xié)同作用的混合燃燒控制策略的研究[D]. 天津:天津大學(xué),2012.

    Yu Wenbin. Research on Hybrid Combustion Control Strategy Based on Interaction of Mixing and Chemical Factors[D]. Tianjin: Tianjin University, 2012. (in Chinese with English abstract)

    [12] Giechaskiel B, Carriero M, Martini G, et al. Heavy Duty Particle Measurement Programme (PMP): Exploratory Work for the Definition of the Test Protocol[J]. Sae International Journal of Engines, 2009, 2(1): 1528-1546.

    [13] Giechaskiel B, Dilara P, Sandbach E, et al. Particle measurement programme (PMP) light-duty inter-laboratory exercise: comparison of different particle number measurement systems[J]. Measurement Science & Technology, 2008, 19(9): 44-46.

    [14] Benajes J, García A, Monsalve-Serrano J, et al. Gaseous emissions and particle size distribution of dual-mode dual-fuel diesel-gasoline concept from low to full load[J]. Applied Thermal Engineering, 2017, 120: 138-149.

    [15] Barrios C C, Domínguez-Sáez A, Martín C, et al. Effects of animal fat based biodiesel on a TDI diesel engine performance, combustion characteristics and particle number and size distribution emissions[J]. Fuel, 2014, 117: 618-623.

    [16] Barrios C C, Martin C, Dominguez-Saez A, et al. Effects of the addition of oxygenated fuels as additives on combustion characteristics and particle number and size distribution emissions of a TDI diesel engine[J]. Fuel, 2014, 132: 93-100.

    [17] 王忠,孫波,趙洋,等. 小型非道路柴油機(jī)排氣管內(nèi)顆粒的粒徑分布與氧化特性[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(10):41-46.

    Wang Zhong, Sun Bo, Zhao Yang, et al. Characteristics of particle coagulation and oxidation in exhaust pipe of diesel engine[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016. (in Chinese with English abstract)

    [18] 李新令,黃震,王嘉松,等. 柴油機(jī)排氣顆粒濃度和粒徑分布特征試驗(yàn)研究[J]. 內(nèi)燃機(jī)學(xué)報(bào),2007,25(2):113-117.

    Li Xinling, Huang Zhen, Wang Jiasong, et al. Investigation on concentrations and size distribution characteristic of particles from diesel engine[J]. Transactions of Csice, 2007, 25(2): 113-117. (in Chinese with English abstract)

    [19] 孟建,劉軍恒,孫平,等. 雙極荷電凝并降低柴油機(jī)顆粒排放數(shù)量的試驗(yàn)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(14):78-84.

    Meng Jian, Liu Junheng, Sun Ping, et al. Experimental research on decreasing particles emission number of diesel engine by bipolar charged coagulation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(14): 78-84. (in Chinese with English abstract)

    [20] Kittelson D B. Engines and nanoparticles: a review[J]. Journal of Aerosol Science, 1998, 29(5/6): 575-588.

    [21] Kazakov A, Foster D E. Modeling of soot formation during di diesel combustion using a multi-step phenomenological model[J]. Sae Technical Papers, 1998.

    [22] Desantes J M, Bermúdez V, García J M, et al. Effects of current engine strategies on the exhaust aerosol particle size distribution from a Heavy-Duty Diesel Engine[J]. Journal of Aerosol Science, 2005, 36(10): 1251-1276.

    [23] Shi J P, Khan A A, Harrison R M. Measurements of ultrafine particle concentration and size distribution in the urban atmosphere[J]. Science of the Total Environment, 1999, 235(1/2/3): 51-64.

    [24] Shi J P, Evans D E, Khan A A, et al. Sources and concentrationof nanoparticles (<10 nm diameter) in the urban atmosphere[J]. Atmospheric Environment, 2001, 35(7): 1193-1202.

    [25] Tan P Q, Ruan S S, Hu Z Y, et al. Particle number emissions from a light-duty diesel engine with biodiesel fuels under transient-state operating conditions[J]. Applied Energy, 2014, 113(113): 22-31.

    [26] Wang Jian, John Storey, Norberto Domingo, et al. Studies of diesel engine particle emissions during transient operations using an engine exhaust particle sizer[J]. Aerosol Science and Technology, 2006, 40(11): 1002-1015.

    [27] Xiaobo Zhou, Erxi Liu, Dezeng Sun, Wanhua Su, Study on transient emission spikes reduction of a heavy-duty diesel engine equipped with a variable intake valve closing timing mechanism and a two-stage turbocharger[J]. International Journal of Engine Research, 2018(13): 146808741774883.

    [28] 戰(zhàn)強(qiáng). 重型柴油機(jī)可變氣門(mén)系統(tǒng)的開(kāi)發(fā)及實(shí)驗(yàn)研究[D]. 天津:天津大學(xué),2013.

    Zhan Qiang. Design of The Variable Valve Timing System for Heavy Duty Diesel Engine and Its Experimental Study [D]. Tianjin: Tianjin University, 2013. (in Chinese with English abstract)

    [29] 劉志華,葛蘊(yùn)珊,丁艷,等. 柴油機(jī)和 LNG 發(fā)動(dòng)機(jī)排放顆粒物粒徑分布特性研究[J]. 內(nèi)燃機(jī)學(xué)報(bào),2009,27(6):518-522.

    Liu Zhihua, Ge Yunshan, Ding Yan, et al. The size distribution characteristic of PM emitted from diesel engine and LNG engine[J]. Transactions of Csice, 2009. (in Chinese with English abstract)

    [30] Ibald-Mulli A, Wichmann H E, Kreyling W, et al. Epidemiological evidence on health effects of ultrafine particles. J Aerosol Med, 2002, 15: 189-201

    [31] Oberd?ster G, Utell M J. Ultrafine particles in the urban air: To the respiratory track and beyond. Environ Health Perspec, 2002, 110 A: 440-441.

    Experimental study on particle distribution of exhaust emission of heavy-duty diesel engine

    Zhou Xiaobo1,2, Hu Qinghua1※, Yan Feng3, Su Wanhua2

    (1.300350,; 2.300072,; 3.,300300,)

    In recent years, smog has emerged in most of major cities in China. Epidemiological and toxicological studies showed that the ultrafine particles in the atmosphere adhere to a large number of toxic substances, which can cause serious harm to human respiratory and cardiovascular and other organs after they enter the human body through the respiratory tract. and the data indicated that it had a significant correlation with human morbidity and mortality. Particulate emissions from diesel engines are an important source of ultrafine particles in urban atmosphere. The researchers found that the weighing method of particulate matter detection cannot accurately reflect the emission levels of internal combustion engines. This new method for measuring the number of particulate matter emissions after removing volatile matter from the method has been widely recognized and included in the scope of regulatory investigation. However, there are few studies on the number and size distribution of particles in the transient process of diesel engine. In order to simultaneously control the quality and quantity of particulate matter emitted from heavy duty diesel engines, the particle distribution of exhaust particulates from diesel engines under steady and transient conditions was studied. The test bench is a 12 L heavy duty diesel engine and an electric dynamometer. The test bench is a prototype with no post-processor and meets the Euro IV regulations. It has been certified by China Automotive Technology Research Center (CATRC). The engine intake system consists of a two-stage turbocharging system with variable geometry section turbine (VGT), a high-pressure EGR system and an intake valve late closing system (IVCT system, which is developed by the research group independently). Cambustion DMS500 fast particle spectrometer was used to analyze the particle size distribution of diesel engine under transient conditions. The experimental results showed that, in the steady state, with the increase of load or speed, the concentration of accumulated particles showed an upward trend, corresponding to the increase of count median diameter (CMD), the surface area and mass of accumulated particles showed an increasing trend, and the peak value shifted to the direction of large particle size. In this study, the total number concentration of particles during the steady state is 1.5×106-4.5×106/cm3, and accumulation mode particle was 2×106-1×107/cm3, nucleation mode particle was 1×107-3×107/cm3which account for 65%-96%. In the steady state, the number of nucleation mode particles is the main source of the total number of particles. The equivalent ratio of all steady state conditions is less than 0.7. The effect of equivalent ratio on particle distribution is obvious at middle and low rotational speed, but weakens at high rotational speed, especially at high load. The European Transient Cycle (ETC) was used for the transient test. In transient condition, the spikes also appear in the number concentration of particles which similar to that of quality of particulate matter, and the number concentration peak even reaches 2×108-7×108/cm3, which is 2 orders of magnitude higher than that of the corresponding steady state operation. The number concentration of nucleation mode particle increases significantly in the spikes, but the proportion of the number concentration of nucleation mode particle in the total particles is reduced. The number concentration of accumulation mode particle is the main part of the spikes which is different from the steady state condition. The particle size peak is mainly concentrated in the 50-200 nm range at accumulation mode, the 6-8 nm and 20-50 nm at nucleation mode, this mainly because that the equivalent ratio in the transient process appears to exceed the critical equivalent ratio of 0.8 working conditions. The equivalent ratio does not exceed 0.7 in steady state operation, but the condition corresponding to the emission peak appears to exceed the critical equivalent ratio 0.8. The results are of great significance for particulate emission control and generation mechanism of heavy-duty diesel engine, and can provide data support for matching calculation of post-exhaust processor.

    diesel engineers; exhaust gas; particulate matter; steady state condition; transient condition; distribution function; nucleation and accumulation mode

    周小波,胡清華,閆 峰,蘇萬(wàn)華. 重型柴油機(jī)顆粒物分布規(guī)律的試驗(yàn)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(13):62-69.doi:10.11975/j.issn.1002-6819.2018.13.008 http://www.tcsae.org

    Zhou Xiaobo, Hu Qinghua, Yan Feng, Su Wanhua. Experimental study on particle distribution of exhaust emission of heavy-duty diesel engine[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(13): 62-69. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.13.008 http://www.tcsae.org

    2018-03-26

    2018-06-06

    周小波,博士后,主要研究柴油機(jī)排放控制及自動(dòng)駕駛。 Email:xiaobo_zhou@tju.edu.cn

    胡清華,教授,博士生導(dǎo)師。主要研究不確定性建模和多模態(tài)學(xué)習(xí)等。Email:huqinghua@tju.edu.cn

    10.11975/j.issn.1002-6819.2018.13.008

    TK421+.5

    A

    1002-6819(2018)-13-0062-08

    猜你喜歡
    顆粒物
    道路空氣顆粒物污染與騎行流的相關(guān)性研究
    H2O(g)對(duì)富氧燃燒超細(xì)顆粒物生成特性影響
    柴油機(jī)顆粒物實(shí)時(shí)診斷傳感器技術(shù)研究進(jìn)展
    某高中校園內(nèi)顆粒物塵源調(diào)查實(shí)測(cè)與分析
    南平市細(xì)顆粒物潛在來(lái)源分析
    某高校室內(nèi)顆粒物成分及污染特征分析
    固定源細(xì)顆粒物監(jiān)測(cè)技術(shù)現(xiàn)狀分析與思考
    MMT 對(duì)直噴汽油機(jī)顆粒物排放影響的實(shí)驗(yàn)研究
    錯(cuò)流旋轉(zhuǎn)填料床脫除細(xì)顆粒物研究
    多層介質(zhì)阻擋放電處理柴油機(jī)尾氣顆粒物
    美女脱内裤让男人舔精品视频| 熟妇人妻不卡中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 国产有黄有色有爽视频| 乱码一卡2卡4卡精品| 少妇丰满av| 午夜福利在线观看免费完整高清在| kizo精华| 男女无遮挡免费网站观看| 中文字幕免费在线视频6| 大片免费播放器 马上看| 观看美女的网站| 国产视频首页在线观看| 久久99蜜桃精品久久| 中文字幕av成人在线电影| 国产老妇女一区| 成人二区视频| 国产精品成人在线| 亚洲自拍偷在线| 2021天堂中文幕一二区在线观| 蜜臀久久99精品久久宅男| 色网站视频免费| 亚洲欧美日韩东京热| 亚洲精品影视一区二区三区av| 国产黄片美女视频| 久热这里只有精品99| 亚洲国产成人一精品久久久| 亚洲高清免费不卡视频| av免费在线看不卡| 国产午夜福利久久久久久| 高清视频免费观看一区二区| 99热这里只有是精品50| 国产爽快片一区二区三区| 在线免费十八禁| 天美传媒精品一区二区| 丰满人妻一区二区三区视频av| av播播在线观看一区| 亚洲精品乱码久久久久久按摩| 丝瓜视频免费看黄片| 一区二区三区免费毛片| 久久99热这里只有精品18| 欧美日韩视频精品一区| 又大又黄又爽视频免费| av国产免费在线观看| 日韩制服骚丝袜av| 一本色道久久久久久精品综合| 久久精品久久久久久噜噜老黄| 全区人妻精品视频| 欧美丝袜亚洲另类| 亚洲精品影视一区二区三区av| 69av精品久久久久久| av卡一久久| 青春草亚洲视频在线观看| 日韩精品有码人妻一区| 国产日韩欧美在线精品| 亚洲aⅴ乱码一区二区在线播放| 舔av片在线| 日本猛色少妇xxxxx猛交久久| 各种免费的搞黄视频| 99热国产这里只有精品6| 深夜a级毛片| 国产爱豆传媒在线观看| 大陆偷拍与自拍| 十八禁网站网址无遮挡 | 久久热精品热| 成人亚洲精品一区在线观看 | 欧美xxxx性猛交bbbb| 一级毛片久久久久久久久女| 亚洲精品色激情综合| 亚洲国产av新网站| 1000部很黄的大片| 成人综合一区亚洲| 久久午夜福利片| 色视频在线一区二区三区| 国精品久久久久久国模美| 菩萨蛮人人尽说江南好唐韦庄| 日日啪夜夜撸| 夫妻性生交免费视频一级片| 欧美激情久久久久久爽电影| 91狼人影院| 新久久久久国产一级毛片| 精品一区二区三卡| 尤物成人国产欧美一区二区三区| 国产中年淑女户外野战色| 亚洲高清免费不卡视频| 日日摸夜夜添夜夜爱| 91在线精品国自产拍蜜月| av在线天堂中文字幕| 一个人观看的视频www高清免费观看| 久久精品久久精品一区二区三区| 观看美女的网站| 久久久久九九精品影院| av国产久精品久网站免费入址| 国模一区二区三区四区视频| 国精品久久久久久国模美| 色播亚洲综合网| 国产免费视频播放在线视频| 国产成人a∨麻豆精品| 日韩成人av中文字幕在线观看| 欧美性猛交╳xxx乱大交人| 在线观看一区二区三区| 久久99蜜桃精品久久| 国产黄色免费在线视频| 九色成人免费人妻av| 久久久成人免费电影| 亚洲精品成人久久久久久| 三级国产精品欧美在线观看| 黄色一级大片看看| 五月开心婷婷网| 内射极品少妇av片p| 波多野结衣巨乳人妻| 亚洲精华国产精华液的使用体验| 中文字幕久久专区| 99久久精品一区二区三区| 国产在线男女| 久久亚洲国产成人精品v| 伦精品一区二区三区| 小蜜桃在线观看免费完整版高清| 亚洲精品一二三| 天堂中文最新版在线下载 | 精品国产露脸久久av麻豆| 国产精品久久久久久久电影| 免费大片黄手机在线观看| 亚洲国产日韩一区二区| 在线精品无人区一区二区三 | 九九爱精品视频在线观看| 久久久久国产精品人妻一区二区| 亚洲国产精品成人久久小说| 国产精品嫩草影院av在线观看| 天天躁日日操中文字幕| 久久久久久久精品精品| 欧美性感艳星| 国产亚洲91精品色在线| 九九久久精品国产亚洲av麻豆| 亚洲激情五月婷婷啪啪| 人妻夜夜爽99麻豆av| 国产有黄有色有爽视频| 欧美区成人在线视频| 免费看光身美女| 色5月婷婷丁香| 一级毛片aaaaaa免费看小| 日本三级黄在线观看| 国产探花在线观看一区二区| 大又大粗又爽又黄少妇毛片口| 国产 一区 欧美 日韩| 亚洲国产精品国产精品| 国产亚洲av片在线观看秒播厂| av黄色大香蕉| 国产免费又黄又爽又色| 久久久精品免费免费高清| 久久久欧美国产精品| 日韩欧美一区视频在线观看 | 久久精品久久久久久噜噜老黄| 九九久久精品国产亚洲av麻豆| 亚洲自拍偷在线| www.色视频.com| 亚洲成人中文字幕在线播放| 亚洲欧美成人综合另类久久久| 色吧在线观看| 啦啦啦在线观看免费高清www| 久久影院123| 丰满乱子伦码专区| 成人国产av品久久久| 午夜日本视频在线| 欧美日韩视频精品一区| 久久久久久久久大av| 日韩不卡一区二区三区视频在线| 最近最新中文字幕免费大全7| 欧美国产精品一级二级三级 | 欧美成人精品欧美一级黄| 熟女人妻精品中文字幕| 免费在线观看成人毛片| 国产欧美日韩一区二区三区在线 | 久热久热在线精品观看| 亚洲精品乱码久久久v下载方式| 亚洲精品第二区| 成人亚洲精品av一区二区| 国产欧美亚洲国产| 国产 一区精品| 久久久精品欧美日韩精品| 婷婷色综合www| 在线免费观看不下载黄p国产| 国产乱来视频区| 在线亚洲精品国产二区图片欧美 | 国产v大片淫在线免费观看| 国产乱人偷精品视频| 欧美日韩国产mv在线观看视频 | 3wmmmm亚洲av在线观看| 777米奇影视久久| 国语对白做爰xxxⅹ性视频网站| 国产美女午夜福利| 如何舔出高潮| av线在线观看网站| 亚洲av成人精品一区久久| 男女无遮挡免费网站观看| 赤兔流量卡办理| 精品一区二区三卡| 国产又色又爽无遮挡免| 日韩国内少妇激情av| 国产黄片美女视频| 永久网站在线| 成年女人看的毛片在线观看| 永久免费av网站大全| 男女啪啪激烈高潮av片| 久久热精品热| 国产黄色免费在线视频| av免费观看日本| 七月丁香在线播放| 亚洲欧洲日产国产| 校园人妻丝袜中文字幕| 最近中文字幕高清免费大全6| 舔av片在线| 亚洲国产精品成人综合色| 亚洲最大成人手机在线| 一本一本综合久久| 麻豆精品久久久久久蜜桃| 性色av一级| 免费黄频网站在线观看国产| 日韩电影二区| 国产伦理片在线播放av一区| www.av在线官网国产| 午夜福利在线观看免费完整高清在| 国产精品一区二区性色av| 亚洲欧美精品自产自拍| av在线蜜桃| 国产精品偷伦视频观看了| 亚洲av成人精品一二三区| 成人国产av品久久久| 在线天堂最新版资源| 大香蕉久久网| 国产老妇女一区| 大话2 男鬼变身卡| 日产精品乱码卡一卡2卡三| 五月伊人婷婷丁香| 99久久中文字幕三级久久日本| 久久久国产一区二区| 校园人妻丝袜中文字幕| 在线观看免费高清a一片| 最近手机中文字幕大全| 天美传媒精品一区二区| 超碰av人人做人人爽久久| 免费观看无遮挡的男女| 国产亚洲av嫩草精品影院| 成年免费大片在线观看| 久久精品国产a三级三级三级| 久久这里有精品视频免费| 五月开心婷婷网| 三级国产精品欧美在线观看| 边亲边吃奶的免费视频| 看非洲黑人一级黄片| 午夜激情久久久久久久| 国产精品欧美亚洲77777| 黑人欧美特级aaaaaa片| 国产精品久久久久久精品电影小说| 三上悠亚av全集在线观看| 日韩熟女老妇一区二区性免费视频| 天天添夜夜摸| 麻豆av在线久日| 日韩一卡2卡3卡4卡2021年| 激情五月婷婷亚洲| 久久综合国产亚洲精品| 亚洲精品国产av成人精品| 精品国产超薄肉色丝袜足j| 国产日韩欧美亚洲二区| 成人午夜精彩视频在线观看| 精品福利永久在线观看| 在线观看免费午夜福利视频| 综合色丁香网| 最近2019中文字幕mv第一页| 精品亚洲成国产av| 日韩欧美一区视频在线观看| 最近中文字幕2019免费版| 视频在线观看一区二区三区| 日韩中文字幕视频在线看片| 欧美少妇被猛烈插入视频| 午夜福利视频精品| 国产精品久久久久久精品古装| 人妻 亚洲 视频| 久久国产精品大桥未久av| 只有这里有精品99| 色综合欧美亚洲国产小说| 久久97久久精品| 侵犯人妻中文字幕一二三四区| 精品亚洲成国产av| 在线天堂中文资源库| 又大又爽又粗| 操美女的视频在线观看| 啦啦啦在线免费观看视频4| 少妇人妻精品综合一区二区| 日本91视频免费播放| 精品国产乱码久久久久久小说| 精品一区二区三区四区五区乱码 | 欧美xxⅹ黑人| 日本av免费视频播放| 日韩大片免费观看网站| 观看美女的网站| 狂野欧美激情性bbbbbb| 老司机靠b影院| 亚洲视频免费观看视频| 51午夜福利影视在线观看| av片东京热男人的天堂| 国产av国产精品国产| 中文乱码字字幕精品一区二区三区| 飞空精品影院首页| 蜜桃在线观看..| 肉色欧美久久久久久久蜜桃| 婷婷色麻豆天堂久久| 亚洲精品久久久久久婷婷小说| 成年人午夜在线观看视频| 午夜福利视频在线观看免费| 国产女主播在线喷水免费视频网站| 久久久久精品国产欧美久久久 | 亚洲图色成人| 亚洲成av片中文字幕在线观看| 日韩av不卡免费在线播放| 各种免费的搞黄视频| 亚洲色图综合在线观看| 侵犯人妻中文字幕一二三四区| 国产淫语在线视频| 在线看a的网站| 一本—道久久a久久精品蜜桃钙片| 亚洲,一卡二卡三卡| 久久人人爽人人片av| 国产成人午夜福利电影在线观看| 亚洲精品视频女| 国产男女内射视频| 好男人视频免费观看在线| 99久久人妻综合| 99精国产麻豆久久婷婷| 免费观看a级毛片全部| 精品一区在线观看国产| 王馨瑶露胸无遮挡在线观看| xxxhd国产人妻xxx| 99久国产av精品国产电影| 下体分泌物呈黄色| 日韩中文字幕欧美一区二区 | 亚洲伊人色综图| 亚洲成av片中文字幕在线观看| 国产精品国产三级专区第一集| 久久久精品94久久精品| 男人爽女人下面视频在线观看| 亚洲图色成人| 黄色视频不卡| 别揉我奶头~嗯~啊~动态视频 | 在线观看人妻少妇| 国产探花极品一区二区| 丝袜人妻中文字幕| svipshipincom国产片| 亚洲欧美清纯卡通| 国产日韩欧美亚洲二区| 国产欧美亚洲国产| 久久性视频一级片| 欧美久久黑人一区二区| 国产1区2区3区精品| 侵犯人妻中文字幕一二三四区| 久热这里只有精品99| 天堂中文最新版在线下载| 精品一区二区三区四区五区乱码 | 国产一区二区在线观看av| 99久久综合免费| 人人妻,人人澡人人爽秒播 | 新久久久久国产一级毛片| 国产精品二区激情视频| 天天操日日干夜夜撸| 国产一区亚洲一区在线观看| 国产精品久久久久久精品电影小说| 久久国产精品男人的天堂亚洲| 桃花免费在线播放| 国产精品嫩草影院av在线观看| 亚洲精品国产av蜜桃| 国产一级毛片在线| 国产伦人伦偷精品视频| av在线播放精品| 午夜免费观看性视频| 最新的欧美精品一区二区| 又黄又粗又硬又大视频| 18禁国产床啪视频网站| 成人国语在线视频| 亚洲三区欧美一区| 亚洲av电影在线进入| 精品午夜福利在线看| av线在线观看网站| 国产黄色免费在线视频| 校园人妻丝袜中文字幕| videos熟女内射| 亚洲人成电影观看| 日韩免费高清中文字幕av| 久久精品熟女亚洲av麻豆精品| 婷婷色综合大香蕉| 久久久久国产精品人妻一区二区| 国产欧美亚洲国产| 777米奇影视久久| 如何舔出高潮| 亚洲激情五月婷婷啪啪| 亚洲综合精品二区| 亚洲精品美女久久av网站| 久久久国产一区二区| 亚洲精品中文字幕在线视频| 亚洲欧美一区二区三区黑人| 久久精品久久精品一区二区三区| 精品国产一区二区三区四区第35| 女人精品久久久久毛片| 免费黄色在线免费观看| 免费高清在线观看日韩| 一区二区三区四区激情视频| 考比视频在线观看| 人人妻人人澡人人看| 97人妻天天添夜夜摸| 精品国产一区二区三区久久久樱花| 日本猛色少妇xxxxx猛交久久| 三上悠亚av全集在线观看| 日韩av在线免费看完整版不卡| 超碰成人久久| 波野结衣二区三区在线| 亚洲国产精品一区三区| 久久久国产欧美日韩av| 波多野结衣一区麻豆| 亚洲三区欧美一区| 男的添女的下面高潮视频| 亚洲精品av麻豆狂野| 美女午夜性视频免费| 男人操女人黄网站| 少妇人妻精品综合一区二区| 五月开心婷婷网| 精品午夜福利在线看| 中文字幕制服av| av网站免费在线观看视频| 久久综合国产亚洲精品| 男人添女人高潮全过程视频| av福利片在线| 成年人午夜在线观看视频| 亚洲精品一二三| 亚洲熟女精品中文字幕| 免费观看性生交大片5| 天天影视国产精品| 最新的欧美精品一区二区| 人人澡人人妻人| 国产深夜福利视频在线观看| 亚洲成人av在线免费| 国产亚洲欧美精品永久| a级毛片黄视频| 亚洲国产精品999| av免费观看日本| 欧美乱码精品一区二区三区| 日韩一区二区视频免费看| 一区在线观看完整版| 国产成人a∨麻豆精品| 欧美日韩精品网址| 女人精品久久久久毛片| 卡戴珊不雅视频在线播放| 久久精品熟女亚洲av麻豆精品| 国产黄频视频在线观看| 久久韩国三级中文字幕| 久热这里只有精品99| 日韩一卡2卡3卡4卡2021年| 久久精品久久久久久久性| 中文字幕人妻丝袜制服| 精品一区在线观看国产| 日本91视频免费播放| 最新的欧美精品一区二区| 精品国产国语对白av| 国产又色又爽无遮挡免| 一本色道久久久久久精品综合| 一区二区三区激情视频| 男女之事视频高清在线观看 | 亚洲精品在线美女| 波野结衣二区三区在线| 精品第一国产精品| 涩涩av久久男人的天堂| 母亲3免费完整高清在线观看| 丁香六月欧美| 亚洲第一青青草原| www.精华液| 在线观看人妻少妇| 日韩一本色道免费dvd| 亚洲 欧美一区二区三区| 国产免费又黄又爽又色| 哪个播放器可以免费观看大片| 欧美 亚洲 国产 日韩一| 亚洲视频免费观看视频| 日本vs欧美在线观看视频| 秋霞在线观看毛片| 免费高清在线观看日韩| 国产高清国产精品国产三级| 另类精品久久| 国产高清不卡午夜福利| 国产精品一区二区在线不卡| 欧美成人精品欧美一级黄| 国产高清不卡午夜福利| 精品一区在线观看国产| 亚洲精品中文字幕在线视频| 男女免费视频国产| 欧美久久黑人一区二区| 亚洲精品日本国产第一区| 18禁国产床啪视频网站| 国产视频首页在线观看| 国产一卡二卡三卡精品 | 成人18禁高潮啪啪吃奶动态图| 国产野战对白在线观看| 91aial.com中文字幕在线观看| 五月天丁香电影| 高清视频免费观看一区二区| 操美女的视频在线观看| 精品国产超薄肉色丝袜足j| 99国产综合亚洲精品| 成人18禁高潮啪啪吃奶动态图| 国产色婷婷99| 日韩大片免费观看网站| 女人被躁到高潮嗷嗷叫费观| 汤姆久久久久久久影院中文字幕| 国产1区2区3区精品| 欧美老熟妇乱子伦牲交| 国产精品亚洲av一区麻豆 | 成人毛片60女人毛片免费| av在线播放精品| 国产午夜精品一二区理论片| 成人国产av品久久久| 久久 成人 亚洲| 国产色婷婷99| 精品一区二区免费观看| 悠悠久久av| 男女边摸边吃奶| 亚洲精品av麻豆狂野| 欧美激情 高清一区二区三区| 啦啦啦在线观看免费高清www| 悠悠久久av| 免费人妻精品一区二区三区视频| 国产亚洲av高清不卡| 又大又爽又粗| 精品少妇一区二区三区视频日本电影 | 亚洲欧美激情在线| 夫妻性生交免费视频一级片| 亚洲色图综合在线观看| 99国产精品免费福利视频| 国产男女内射视频| 国产在线一区二区三区精| 亚洲精品久久久久久婷婷小说| 一区二区三区乱码不卡18| 黄色视频不卡| 久久热在线av| 国产在线一区二区三区精| 国产黄频视频在线观看| 国产精品嫩草影院av在线观看| 亚洲人成77777在线视频| 成人漫画全彩无遮挡| avwww免费| 成人国产av品久久久| 亚洲成人手机| 国产又色又爽无遮挡免| 午夜激情久久久久久久| 韩国精品一区二区三区| 中文字幕人妻熟女乱码| 丰满饥渴人妻一区二区三| 黄频高清免费视频| 精品国产国语对白av| 亚洲精品美女久久久久99蜜臀 | 精品久久蜜臀av无| 精品亚洲成a人片在线观看| 青草久久国产| 国产亚洲av高清不卡| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 天天添夜夜摸| 国产免费一区二区三区四区乱码| 高清视频免费观看一区二区| 亚洲国产中文字幕在线视频| 国产xxxxx性猛交| 18禁裸乳无遮挡动漫免费视频| 中文字幕制服av| av免费观看日本| 亚洲av综合色区一区| 日韩人妻精品一区2区三区| 热99国产精品久久久久久7| av有码第一页| 久久99热这里只频精品6学生| 亚洲成人国产一区在线观看 | av天堂久久9| 无遮挡黄片免费观看| 精品午夜福利在线看| h视频一区二区三区| av网站在线播放免费| 欧美黑人精品巨大| 男女国产视频网站| 人妻一区二区av| 久久久久精品国产欧美久久久 | 欧美日韩亚洲综合一区二区三区_| 亚洲欧美一区二区三区久久| 黑丝袜美女国产一区| xxx大片免费视频| 亚洲专区中文字幕在线 | 精品第一国产精品| 宅男免费午夜| 国产亚洲最大av| 成人午夜精彩视频在线观看| 色婷婷av一区二区三区视频| 亚洲精品日本国产第一区| 在线精品无人区一区二区三| 亚洲熟女精品中文字幕| 亚洲婷婷狠狠爱综合网| 黄色视频不卡| 高清不卡的av网站| 一本一本久久a久久精品综合妖精| 国产视频首页在线观看| 午夜免费男女啪啪视频观看| 国产精品欧美亚洲77777| 欧美97在线视频| 亚洲国产欧美在线一区| 无遮挡黄片免费观看| 亚洲成国产人片在线观看| 午夜影院在线不卡| 色婷婷av一区二区三区视频| 久久久久久人妻| 最黄视频免费看| 亚洲精品在线美女| 亚洲天堂av无毛|