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

    煤及油母中常見(jiàn)C—X(X=N,S)的解離能

    2017-10-16 08:13:28李璐樊紅軍胡浩權(quán)
    化工學(xué)報(bào) 2017年10期

    李璐,樊紅軍,胡浩權(quán)

    (1中國(guó)科學(xué)院大連化學(xué)物理研究所分子反應(yīng)動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116023;2大連理工大學(xué)化工學(xué)院,精細(xì)化工國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116024)

    煤及油母中常見(jiàn)C—X(X=N,S)的解離能

    李璐1,樊紅軍1,胡浩權(quán)2

    (1中國(guó)科學(xué)院大連化學(xué)物理研究所分子反應(yīng)動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116023;2大連理工大學(xué)化工學(xué)院,精細(xì)化工國(guó)家重點(diǎn)實(shí)驗(yàn)室,遼寧 大連 116024)

    煤及油頁(yè)巖中除了碳、氫和氧原子是主要的組成元素外,氮和硫等雜原子大都以C—X(X=N,S)的鍵合形式存在,在其結(jié)構(gòu)及轉(zhuǎn)化利用中同樣發(fā)揮重要的作用。認(rèn)識(shí)C—X的解離能,有助于建立氮、硫熱解過(guò)程中的遷移模型,豐富對(duì)煤及油母中常見(jiàn)化學(xué)鍵性質(zhì)的認(rèn)識(shí),對(duì)發(fā)展高效清潔的能源利用技術(shù)至關(guān)重要。利用雙雜化密度泛函方法,系統(tǒng)研究了煤及油母中典型 C—X鍵的解離能(BDE)范圍。研究結(jié)果表明,煤及油母中常見(jiàn) C—N和 C—S 鍵的 BDE 值范圍分別是 154.1~55.7 kcal·mol?1和 83.0~56.6 kcal·mol?1。在熱解過(guò)程中,苯硫類自由基以及苯胺類自由基會(huì)在初期產(chǎn)生,其次才是巰基、胺類自由基等側(cè)鏈取代自由基脫落。C—S鍵的整體BDE值范圍比其他類型化學(xué)鍵更低,各類化學(xué)鍵最低BDE值的高低順序符合O—H > C—H > C—C > C—N > C—S > C—O的規(guī)律,其中,只有當(dāng)有PhO·生成時(shí),C—S > C—O,否則C—O > C—S。

    煤;油母;解離能;焓;熱力學(xué);計(jì)算化學(xué)

    Abstract:The carbon,hydrogen,oxygen are the majority elements in the structure of coal and oil shale,while the nitrogen and sulfur atoms,mostly in the form of C—X (X=N,S),still play a nonnegligible role in the structure and application of coal and oil shale.The homolytic bond dissociation enthalpies (BDE) of C—X for model compounds that are representative of the functionalities present in coal and kerogen were computed by using a double-hybrid method mPW2PLYP.The studies will be helpful for people to construct the transfer model of N and S in pyrolysis,and enrich the knowledge of chemical reaction of typical bonds in coal and kerogen.The BDE for C—N and C—S cover a range from 154.1—55.7 kcal·mol?1and 83.0—56.6 kcal·mol?1,respectively.The results suggest that the thiophenyl radical or phenylamino radical are the most favorable intermediates in the early pyrolytic stage.With the increase of temperature,the loss of sulfydryl and amidogen groups become feasible,and the ring cleavage reaction of five-membered ring (pyrrole and thiophene) through homolysis of C—N and C—S bond can also occur.The phenyl radical is the most difficult to form,and the directly bond dissociation of C—N bond in pyridine with six-membered ring is extremely unfavorable in the pyrolysis of coal and kerogen.Comparedto other types of bonds,C—S bond has the lowest general BDE range.The order of the lowest BDE for various bonds is O—H > C—H > C—C > C—N > C—S > C—O.The bond dissociation of C—O to generate PhO·radical is easier than to generate PhS· radical through homolysis of C—S,otherwise the BDE of C—O is higher than C—S.

    Key words:coal;kerogen;bond dissociation enthalpies;enthalpy;thermodynamics;computational chemistry

    引 言

    能源是人類活動(dòng)的物質(zhì)基礎(chǔ),人類社會(huì)的前進(jìn)與發(fā)展離不開(kāi)優(yōu)質(zhì)能源的開(kāi)采使用和先進(jìn)能源技術(shù)的開(kāi)發(fā)利用。我國(guó)缺油、少氣和煤炭相對(duì)豐富的能源結(jié)構(gòu)決定了煤在我國(guó)能源消費(fèi)結(jié)構(gòu)中的主導(dǎo)地位。油頁(yè)巖是重要的石油補(bǔ)充能源,近年來(lái)得到廣泛關(guān)注[1-2]。低溫?zé)峤馐且环N煤及油頁(yè)巖常用的加工利用方式,但是由于煤及油頁(yè)巖本身的成分特點(diǎn)和現(xiàn)有工藝技術(shù)的缺陷,目前的熱解技術(shù)普遍存在油氣收率低、重質(zhì)組分含量高等技術(shù)難題,并且加工過(guò)程中由組分中的氮和硫轉(zhuǎn)化生成的NOx、H2S和SO2等被普遍認(rèn)為是導(dǎo)致酸雨的重要原因[3],N2O則能夠?qū)Τ粞鯇赢a(chǎn)生破壞[4],這些都會(huì)對(duì)環(huán)境造成極大的破壞。對(duì)煤及油頁(yè)巖的油母結(jié)構(gòu)中常見(jiàn)C—X(X=N,S)化學(xué)鍵在熱場(chǎng)中的解離行為的深入認(rèn)識(shí)有助于建立它們?cè)跓峤膺^(guò)程中的硫、氮遷移模型[5],對(duì)于發(fā)展高效清潔的能源利用技術(shù)具有重要的指導(dǎo)意義。

    由于煤及油母初始結(jié)構(gòu)的復(fù)雜性和可變性,選擇能代表其特征官能團(tuán)性質(zhì)的模型化合物進(jìn)行研究,有助于從分子水平理解它們的化學(xué)性質(zhì)。Ling等研究了含氮原子的喹啉和異喹啉[6]以及含硫原子的噻吩[7]和苯硫酚[8]的熱解反應(yīng)機(jī)理,并指出分子內(nèi)氫轉(zhuǎn)移反應(yīng)在模型化合物的熱解過(guò)程中對(duì)產(chǎn)物的形成和分配起到了非常重要的作用,苯硫酚熱解比噻吩熱解更容易生成H2S。Dubnikova等[9]詳細(xì)研究了鄰位喹啉自由基和鄰位異喹啉自由基之間的異構(gòu)化過(guò)程以及它們各自后續(xù)的熱解反應(yīng)路徑,結(jié)果發(fā)現(xiàn)兩種自由基在熱解反應(yīng)初期會(huì)經(jīng)歷分子內(nèi)異構(gòu)化反應(yīng)迅速達(dá)到平衡,使得兩種自由基的最終熱解產(chǎn)物及分布完全一致,并且計(jì)算得出的穩(wěn)定產(chǎn)物分布與實(shí)驗(yàn)結(jié)果有良好的一致性。Liang等[10]利用ReaxFF分子動(dòng)力學(xué)模擬研究了木質(zhì)素?zé)峤夂图託錈峤庵械牧蜻w移機(jī)理,指出單純熱解中硫原子會(huì)以雜環(huán)結(jié)構(gòu)形式存在于熱解產(chǎn)物中,而在有H2條件下的加氫熱解中,H2會(huì)促使硫原子從含硫雜環(huán)形態(tài)轉(zhuǎn)變成苯硫酚/醇形態(tài),從而更易于轉(zhuǎn)變成H2S,以達(dá)到脫除硫的目的。并且,在脫硫過(guò)程中,合理控制溫度可達(dá)到最大限度脫硫的目的。

    除了上述對(duì)含氮或硫模型化合物的具體產(chǎn)物生成機(jī)理研究之外,對(duì)于C—X(X=N,S)鍵解離能(bond dissociation enthalpies,BDE)的研究也引起了廣泛關(guān)注。Khrapkovskii等[11]考察了多種方法水平下硝基苯和取代硝基苯類化合物的C—NO2鍵的BDE值,提出B3LYP/6-31G(d,p)方法的計(jì)算結(jié)果能很好地吻合實(shí)驗(yàn)值,苯環(huán)上的吸電子取代基會(huì)使C—NO2鍵的 BDE值降低。Yu等[12]評(píng)價(jià)了多種密度泛函方法(density functional theory,DFT)和基組在計(jì)算有機(jī)反應(yīng)中常見(jiàn)的C—SO2R鍵解離BDE值時(shí)的表現(xiàn),指出 M06-2X/6-31G(d)的計(jì)算結(jié)果最為可靠。Yao等[13]則通過(guò)對(duì)多種化學(xué)鍵型(C—H、O—H、C—C、C—N、C—S等)BDE值的計(jì)算,得出MPW1P86方法在多種DFT方法中表現(xiàn)最為可靠。由此可見(jiàn),雖然已有文獻(xiàn)關(guān)注過(guò)C—S和C—N鍵的反應(yīng)性質(zhì)[14-22],但研究所涉及到的煤及油母結(jié)構(gòu)中的C—N和C—S鍵類型并不完善,并且未對(duì)其解離行為規(guī)律進(jìn)行系統(tǒng)研究。此外,各文獻(xiàn)選擇的最優(yōu)理論方法并不一致,有關(guān)煤及油母中化學(xué)鍵的 BDE也并未在同一計(jì)算水平下得到,以致不能直接采用已有數(shù)值進(jìn)行相互之間的比較。鑒于熱場(chǎng)中的鍵解離是反應(yīng)的第一步,因此本文選取合理的模型化合物,對(duì)煤及油母中 C—X(X=N,S)的BDE值進(jìn)行系統(tǒng)研究,其結(jié)果勢(shì)必會(huì)對(duì)理解N和S在后續(xù)過(guò)程中的遷移途徑有積極的推進(jìn)作用,并能補(bǔ)充對(duì)整個(gè)煤及油母中各個(gè)類型化學(xué)鍵熱反應(yīng)性質(zhì)的認(rèn)識(shí)。

    1 計(jì)算方法

    分子A—B的鍵解離能(BDE)被定義為在標(biāo)準(zhǔn)狀況下,氣相反應(yīng)式(1)中分子A—B鍵均裂生成A·和B·兩個(gè)自由基所需的能量[23]

    BDE的值只與反應(yīng)物和生成物的相對(duì)焓值有關(guān),可通過(guò)式(2)的計(jì)算獲得

    其中,H(A—B)為分子A—B的標(biāo)準(zhǔn)摩爾生成焓,H(A·)和H(B·)分別為自由基 A·和 B·的標(biāo)準(zhǔn)摩爾生成焓。吡啶等雜環(huán)開(kāi)環(huán)的鍵解離能則通過(guò)鍵斷裂后形成雙自由基與母體分子的焓值差進(jìn)行計(jì)算,以吡啶為例,其開(kāi)環(huán)的C—N鍵的BDE值由反應(yīng)式(3)的能量決定

    本文中所有的計(jì)算都是在Gaussian 09量化軟件包下完成,所有分子的優(yōu)化都是在沒(méi)有對(duì)稱性限制下進(jìn)行,每一個(gè)優(yōu)化結(jié)構(gòu)經(jīng)過(guò)頻率分析確認(rèn)沒(méi)有虛頻。根據(jù)作者前期對(duì)解離能計(jì)算方法評(píng)估的研究結(jié)果[24],本文中選取雙雜化密度泛函mPW2PLYP[25]方法在 cc-pVDZ[26]基組下進(jìn)行構(gòu)型優(yōu)化和頻率計(jì)算,并對(duì)優(yōu)化結(jié)構(gòu)在 cc-pVTZ[27]基組下進(jìn)行了單點(diǎn)能量計(jì)算。所有的熱力學(xué)校正數(shù)據(jù)均由頻率計(jì)算的結(jié)果獲得。為了使計(jì)算結(jié)果更準(zhǔn)確,計(jì)算中包括了基組重疊誤差(basis set superposition error,BSSE)校正[28-29]。本文中所有能量單位均為 kcal·mol?1(1 kcal·mol?1=4.184 kJ·mol?1)。

    文中選擇了煤及油母中常見(jiàn)C—N和C—S鍵的模型化合物進(jìn)行研究,所選擇的化合物結(jié)構(gòu)見(jiàn)圖1。主要包括含氮或含硫側(cè)鏈以及橋鍵的芳烴化合物、含氮或含硫雜環(huán)化合物,如吡啶、吡咯和噻吩等。

    圖1 所選模型化合物Fig.1 Model compounds under study in this work

    2 實(shí)驗(yàn)結(jié)果與討論

    2.1 C—N鍵的解離能范圍

    表1為計(jì)算得到的不同化合物中 C—N鍵的BDE值。C—N鍵的 BDE值范圍為 154.1~55.7 kcal·mol?1。其中,吡啶環(huán)的C—N鍵斷裂的解離能最高,達(dá)到 154.1 kcal·mol?1,而 PhCH2—NHPh 鍵斷裂的BDE值最低。與先前研究報(bào)道的C—C鍵型和C—O鍵型的BDE值規(guī)律相似[30],生成Ph·自由基的 C—N鍵解離能量較高,其 BDE值范圍為82.1~99.8 kcal·mol?1。有 PhCH2·自由基和類芐基的o,m,p-C5H4NCH2·自由基生成的 C—N鍵的斷裂所需能量最低,其 BDE值范圍為 55.7~71.2 kcal·mol?1。化合物中甲基取代的胺類自由基解離出CH3·自由基的 C—N 鍵 BDE值范圍為 63.2~73.3 kcal·mol?1,其中處于橋鍵位置的 N—CH3鍵解離比苯環(huán)側(cè)鏈取代的N—CH3鍵解離所需能量更低,因此更易發(fā)生。由吡啶環(huán)替代芐胺中苯環(huán)得到的o,m,p-C5H4NCH2—NH2中的C—N鍵的BDE值分別比PhCH2—NH2高 1.8、0.6 和 1.1 kcal·mol?1。值得注意的是,不同于吡啶六元環(huán)開(kāi)環(huán)需要非常高的能量值,吡咯的 C—N鍵解離所需能量為 84.5 kcal·mol?1。各種 C—N 鍵的 BDE 值的順序見(jiàn)圖2。

    表1 不同化合物中C—N鍵的解離能Table 1 BDE of C—N bonds/kcal·mol?1

    2.2 C—S鍵的解離能范圍

    表2中列出了計(jì)算得到的不同化合物中 C—S鍵 BDE值。C—S鍵的BDE值范圍為 83.0~56.6 kcal·mol?1。其中,Ph—SH 鍵斷裂產(chǎn)生的 BDE 值最高,而PhCH2—SH鍵斷裂產(chǎn)生的BDE值最低。與C—N鍵型的BDE值規(guī)律相似,生成Ph·自由基的C—S鍵解離能量較高,其BDE值范圍為80.0~83.0 kcal·mol?1。有 PhCH2·自由基生成的 PhCH2—S 鍵的斷裂所需能量最低,其 BDE值范圍為 56.6~67.4 kcal·mol?1。顯而易見(jiàn),這些 C—S 鍵的 BDE 值明顯低于相應(yīng)的C—N鍵BDE值。噻吩的C—S鍵斷裂所需能量為 80.8 kcal·mol?1,該值與吡咯的 C—N 鍵BDE值較為接近。鄰位取代的吡啶環(huán)上發(fā)生甲硫基解離所需能量明顯低于間位和對(duì)位取代,該現(xiàn)象與Hayes等[31]計(jì)算的甲基在吡啶環(huán)上不同取代位的解離能規(guī)律一致。各種C—S鍵的BDE值的順序見(jiàn)圖3。

    圖 2 C—N 鍵的解離能區(qū)間(50~155 kcal·mol?1)Fig.2 C—N BDE range from 50 to 155 kcal·mol?1

    表2 不同化合物中C—S鍵的解離能Table 2 BDE of C—S bonds/kcal·mol?1

    2.3 煤及油母中常見(jiàn)化學(xué)鍵的解離能范圍

    圖 3 C—S 鍵的解離能區(qū)間(50~85 kcal·mol?1)Fig.3 C—S BDE range from 50 to 85 kcal·mol?1

    圖4 煤及油母中常見(jiàn)鍵的解離能區(qū)間(括號(hào)中的數(shù)值為所計(jì)算鍵型的數(shù)量)Fig.4 BDE range of different bond types in coal and kerogen (data in parentheses is number of bond type calculated)

    結(jié)合先前的研究結(jié)果[30],可以得到煤及油母中幾種常見(jiàn)化學(xué)鍵的BDE值范圍(圖4)。C—H、C—C、C—O、O—H、C—N和C—S鍵的BDE值范圍分別為 111.4~81.2 kcal·mol?1、102.9~62.8 kcal·mol?1、107.6~52.6 kcal·mol?1、111.2~86.6 kcal·mol?1、154.1~55.7 kcal·mol?1和 83.0~56.6 kcal·mol?1。可以看出,每一種鍵型的 BDE值都是一個(gè)很大的范圍,并且不同鍵型之間的 BDE值存在很大程度的重疊。有H·自由基生成的O—H和C—H鍵的BDE值普遍較高,C—C、C—O、C—N和 C—S鍵的最低 BDE值均來(lái)自于 PhCH2—YPh(Y=CH2、O、NH 和 S)。C—C、C—N 和 C—S鍵的高能量區(qū)域均由生成苯基和吡啶雜環(huán)的芳環(huán)自由基的鍵解離反應(yīng)組成,而在C—O鍵的BDE值區(qū)間內(nèi),芳香羧酸類及芳香脂類化合物中羧基基團(tuán)中C—O鍵斷裂所需能量最高,其次才是有芳環(huán)自由基生成相關(guān)的C—O解離。從圖4可以看出,每種鍵型的最低BDE值的高低順序?yàn)镺—H > C—H >C—C > C—N > C—S > C—O,結(jié)合先前研究成果,將同種類型化學(xué)鍵斷裂的C—O和C—S解離能進(jìn)行比較發(fā)現(xiàn),有PhS·自由基生成的C—S鍵的BDE值總是比有PhO·生成的C—O鍵的解離能高,即C—SPh > C—OPh,這得益于PhO·自由基中的O可以與苯環(huán)形成良好的共軛,使得PhO·自由基具有更高的穩(wěn)定性,從而有更低的 BDE值。而在其他同種類型的BDE值比較中,C—S鍵的BDE值均低于對(duì)應(yīng)的C—O鍵的BDE值。Shi等[32]提出脂肪醚和脂肪醇類化合物中的C—O鍵鍵長(zhǎng)遠(yuǎn)小于脂肪硫醚和硫醇類化合物中的C—S鍵鍵長(zhǎng)是導(dǎo)致這一現(xiàn)象的直接原因。

    2.4 熱解中化學(xué)鍵的斷裂順序

    鍵解離是熱解過(guò)程的第一步,解離能的大小直接關(guān)系煤及油母中化學(xué)鍵斷裂的先后順序,從而影響后續(xù)的整個(gè)熱解過(guò)程,因此對(duì)煤及油母中常見(jiàn)C—X(X=N,S)化學(xué)鍵解離能的系統(tǒng)研究不僅有利于了解N,S原子的遷移途徑,而且對(duì)煤及油母結(jié)構(gòu)中化學(xué)鍵以及熱解反應(yīng)過(guò)程的認(rèn)識(shí)具有重要的作用。

    綜合先前對(duì)煤中常見(jiàn)4種類型(C—H、C—C、C—O和O—H)化學(xué)鍵BDE值的研究結(jié)果[30]可知,煤及油母熱解過(guò)程中應(yīng)首先斷裂由苯氧類自由基、芐基類自由基、苯硫類自由基以及苯胺類自由基相互連接而形成的弱橋鍵。隨著溫度的逐漸升高,芳環(huán)的側(cè)鏈取代基開(kāi)始解離,巰基相比其他小自由基更易解離,隨后是胺類自由基、烷氧類自由基和烷基類自由基的生成,而氫自由基的直接生成較為困難。同時(shí),在此溫度區(qū)間內(nèi),五元環(huán)的吡咯和噻吩也有可能會(huì)發(fā)生開(kāi)環(huán)反應(yīng)。此外,生成苯自由基或類似的芳環(huán)自由基同樣需要較高的能量,因此這類反應(yīng)只能在高溫下發(fā)生。由于六元環(huán)的吡啶直接開(kāi)環(huán)所需的鍵解離能量過(guò)高,可以推測(cè)該六元芳香環(huán)并不是通過(guò)直接鍵斷裂的開(kāi)環(huán)反應(yīng)參與熱解反應(yīng)過(guò)程。

    3 結(jié) 論

    (1)煤及油母中常見(jiàn)C—N和C—S鍵的BDE值范圍分別是 154.1~55.7 kcal·mol?1和 83.0~56.6 kcal·mol?1,與 C—H、C—C、C—O 和 O—H 鍵的BDE值范圍間同樣存在很大程度的重疊。

    (2)煤及油母中的N和S傾向于首先形成苯硫類自由基以及苯胺類自由基,巰基比胺類自由基更易形成。

    (3)煤及油母中常見(jiàn)鍵型的最低 BDE值的高低順序符合O—H > C—H > C—C > C—N > C—S >C—O的規(guī)律。當(dāng)有PhO·生成時(shí),C—S > C—O,否則 C—O > C—S。

    [1]柏靜儒,林衛(wèi)生,潘朔,等.油頁(yè)巖低溫?zé)峤膺^(guò)程中輕質(zhì)氣體的析出特性[J].化工學(xué)報(bào),2015,66(3): 1104-1110.BAI J R,LIN W S,PAN S,et al.Characteristics of light gases evolution during oil shale pyrolysis[J].CIESC Journal,2015,66(3):1104-1110.

    [2]DYNI J R.Geology and resources of some world oil shale oil deposits[J].Shale,2003,20(3): 193-252.

    [3]GAFFNEY J S,STREIT G E,SPALL W D,et al.Beyond acid rain.Do soluble oxidants and organic toxins interact with SO2and NOxto increase ecosystem effects [J].Environmental Science Technology,1987,21(6): 519-524.

    [4]余岳溪,高正陽(yáng),季鵬,等.煤焦異相還原N2O的反應(yīng)機(jī)理[J].化工學(xué)報(bào),2017,68(1): 369-374.YU Y X,GAO Z Y,JI P,et al.Heterogeneous reduction reaction of N2O by char[J].CIESC Journal,2017,68(1): 369-374.

    [5]謝建軍,楊學(xué)民,呂雪松,等.煤熱解過(guò)程中硫氮及遷移規(guī)律研究進(jìn)展[J].化工進(jìn)展,2004,23(11): 1214-1218.XIE J J,YANG X M,Lü X S,et al.Progress on transformation behavior of sulfur and nitrogen during coal pyrolysis[J].Chemical Industry and Engineering Progress,2004,23(11): 1214-1218.

    [6]LING L,ZHANG R G,WANG B J,et al.Pyrolysis mechanisms of quinoline and isoquinoline with density functional theory[J].Chinese Journal of Chemical Engineering,2009,17(5): 805-813.

    [7]LING L,ZHANG R G,WANG B J,et al.Density functional theory study on the pyrolysis mechanism of thiophene in coal[J].Journal of Molecular Structure (THEOCHEM),2009,905(1/2/3): 8-12.

    [8]LING L,ZHANG R G,WANG B J,et al.DFT study on the sulfur migration during benzenethiol pyrolysis in coal[J].Journal of Molecular Structure (THEOCHEM),2009,952(1/2/3): 31-35.

    [9]DUBNIKOVA F,TAMBURU C,LIFSHITZ A.A deep insight into the details of the interisomerization and decomposition mechanism ofo-quinolyl ando-isoquinolyl radicals.Quantum chemical calculations and computer modeling[J].The Journal of Physical Chemistry A,2016,120(38): 7538-7547.

    [10]LI G Y,WANG F,WANG J P,et al.ReaxFF and DFT study on the sulfur transformation mechanism during the oxidation process of lignite[J].Fuel,2016,181: 238-247.

    [11]KHRAPKOVSKII G M,SHARIPOV D D,SHAMOV A G,et al.Theoretical study of substituents effect on C—NO2bond strength in mono substituted nitrobenzenes[J].Computational and Theoretical Chemistry,2013,1017: 7-13.

    [12]YU H Z,FU F,ZHANG L,et al.Accurate predictions of C—SO2R bond dissociation enthalpies using density functional theorymethods[J].Physical Chemistry Chemical Physics,2014,16:20964-20970.

    [13]YAO X Q,HOU X J,JIAO H J,et al.Accurate calculations of bond dissociation enthalpies with density functional methods[J].The Journal of Physical Chemistry A,2003,107(46): 9991-9996.

    [14]DEL GIACCO T,LANZALUNGA O,LAPI A,et al.Photosensitized oxidation of aryl benzyl sulfoxides.Evidence for nucleophilic assistance to the C—S bond cleavage of aryl benzyl sulfoxide radical cations[J].Journal of Organic Chemistry,2015,80(4): 2310-2318.

    [15]KHRAPKOVSKII G M,SHARIPOV D D,SHAMOV A G,et al.Enthalpies of formation of mono substituted nitrobenzenes: a quantum chemistry study[J].Computational and Theoretical Chemistry,2013,1011: 37-43.

    [16]FAYET G,ROUTUREAU P,MINISINI B.Decomposition mechanisms of trinitroalkyl compounds: a theoretical study from aliphatic to aromatic nitro compounds[J].Physical Chemistry Chemical Physics,2014,16(14): 6614-6622.

    [17]MAROCHKIN I I,DOROFEEVA O V.Amide bond dissociation enthalpies: effect of substitution on N—C bond strength[J].Computational and Theoretical Chemistry,2012,991: 182-191.

    [18]JOHNSON E R,CLARKIN O J,DILABIO G A.Density functional theory based model calculations for accurate bond dissociation enthalpies(3): A single approach for X—H,X—X,and X—Y (X,Y=C,N,O,S,halogen) bonds[J].The Journal of Physical Chemistry A,2003,107(46): 9953-9963.

    [19]SONG X L,PARISH C A.Pyrolysis mechanisms of thiophene and methylthiophene in asphaltenes[J].The Journal of Physical Chemistry A,2011,115(13): 2882-2891.

    [20]ZHENG W R,FU Y,GUO Q X.G3//BMK and its application to calculation of bond dissociation enthalpies[J].Journal of Chemical Theory and Computation,2008,4(8): 1324-1331.

    [21]TAO X X,TANG L F,XIE M H,et al.Dielectric properties analysis of sulfur-containing models in coal and energy evaluation of their sulfur-containing bond dissociation in microwave field[J].Fuel,2016,181: 1027-1033.

    [22]TANG L F,WANG S W,GUO J F,et al.Exploration on the removal mechanism of sulfur ether model compounds for coal by microwave irradiation with peroxyacetic acid[J].Fuel Processing Technology,2017,159: 442-447.

    [23]BARCKHOLTZ C,BARCKHOLTZ T A,HADAD C M.C–H and N—H bond dissociation energies of small aromatic hydrocarbons[J].Journal of the American Chemical Society,1999,121(3): 491-500.

    [24]LI L,FAN H J,HU H Q.Assessment of contemporary theoretical methods for bond dissociation enthalpies[J].Chinese Journal of Chemical Physics,2016,29(4): 453-461.

    [25]SCHWABE T,GRIMME S.Towards chemical accuracy for the thermodynamics of large molecules: new hybrid density functionals including non-local correlation effects[J].Physical Chemistry Chemical Physics,2006,8(38): 4398-4401.

    [26]DUNNING J T H.Gaussian basis sets for use in correlated molecular calculations(Ⅰ): The atoms boron through neon and hydrogen[J].The Journal of Chemical Physics,1989,90(2): 1007-1023.

    [27]KENDALL R A,DUNNING J T H,HARRISON R J.Electron affinities of the first-row atoms revisited.Systematic basis sets and wave functions[J].The Journal of Chemical Physics,1992,96(9):6796-6806.

    [28]BOYS S F,BERNARDI F.Calculation of small molecular interactions by differences of separate total energies—some procedures with reduced errors[J].Molecular Physics,1970,19(4):553.

    [29]SIMON S,DURAN M,DANNENBERG J J.How does basis set superposition error change the potential surfaces for hydrogen bonded dimers[J].Journal of Chemical Physics,1996,105(24): 11024-11031.

    [30]LI L,FAN H J,HU H Q.A theoretical study on bond dissociation enthalpies of coal based model compounds[J].Fuel,2015,153: 70-77.

    [31]HAYES C J,HADAD C M.Combustion pathways of the alkylated heteroaromatics: bond dissociation enthalpies and alkyl group fragmentations[J].The Journal of Physical Chemistry A,2009,1133(45): 12370-12379.

    [32]SHI J,HU X R,LIANG S.A computational study of C—S bond dissociation enthalpies in petroleum chemistry[J].Heteroatom Chemistry,2011,22(2): 97-105.

    Bond dissociation enthalpies of C—X (X=N,S) in coal and kerogen

    LI Lu1,FAN Hongjun1,HU Haoquan2
    (1State Key Laboratory of Molecular Reaction Dynamics,Dalian Institute of Chemical Physics,Chinese Academy of Sciences,Dalian116023,Liaoning,China;2State Key Laboratory of Fine Chemicals,Institute of Coal Chemical Engineering,School of Chemical Engineering,Dalian University of Technology,Dalian116024,Liaoning,China)

    TQ 021.2

    A

    0438—1157(2017)10—3900—06

    10.11949/j.issn.0438-1157.20170377

    2017-04-10收到初稿,2017-06-23收到修改稿。

    聯(lián)系人:樊紅軍,胡浩權(quán)。

    李璐(1990—),女,博士后。

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFB0600301)。

    Received date:2017-04-10.

    Corresponding author:Prof.FAN Hongjun,fanhj@dicp.ac.cn;Prof.HU Haoquan,hhu@dlut.edu.cn

    Foundation item:supported by the National Key R&D Program of China (2016YFB0600301).

    国产精品国产三级国产专区5o| 国产男女超爽视频在线观看| 国产黄频视频在线观看| 男女国产视频网站| 热99国产精品久久久久久7| 国产无遮挡羞羞视频在线观看| 日韩视频在线欧美| 男女之事视频高清在线观看 | 亚洲精品久久午夜乱码| 亚洲精品中文字幕在线视频| 日韩,欧美,国产一区二区三区| 久久久亚洲精品成人影院| 狂野欧美激情性xxxx| 国产一区二区 视频在线| 下体分泌物呈黄色| 国产精品久久久人人做人人爽| 精品免费久久久久久久清纯 | 久久精品aⅴ一区二区三区四区| 99精国产麻豆久久婷婷| 少妇人妻久久综合中文| 精品国产国语对白av| 亚洲av综合色区一区| 人妻 亚洲 视频| 肉色欧美久久久久久久蜜桃| 视频区图区小说| 国产日韩一区二区三区精品不卡| 少妇 在线观看| 午夜福利视频在线观看免费| 久久精品人人爽人人爽视色| 国产熟女午夜一区二区三区| 久久 成人 亚洲| 婷婷色av中文字幕| 精品人妻熟女毛片av久久网站| 汤姆久久久久久久影院中文字幕| 国产高清videossex| 老司机亚洲免费影院| 十八禁高潮呻吟视频| 女人久久www免费人成看片| 青草久久国产| 看免费av毛片| 91字幕亚洲| 在线精品无人区一区二区三| 久久精品亚洲av国产电影网| 国产成人av激情在线播放| 又紧又爽又黄一区二区| 另类亚洲欧美激情| 国产黄频视频在线观看| 亚洲精品久久午夜乱码| 久久久久久久国产电影| 国产男女超爽视频在线观看| 爱豆传媒免费全集在线观看| 欧美精品亚洲一区二区| 亚洲精品一二三| 亚洲av国产av综合av卡| 黄片播放在线免费| 热99国产精品久久久久久7| 亚洲精品在线美女| 一区二区日韩欧美中文字幕| 日韩,欧美,国产一区二区三区| 国产精品一国产av| 亚洲中文日韩欧美视频| 亚洲国产毛片av蜜桃av| 久久99一区二区三区| 首页视频小说图片口味搜索 | 亚洲黑人精品在线| 精品一区在线观看国产| 国产精品一二三区在线看| 国产精品久久久人人做人人爽| 亚洲黑人精品在线| 黄色a级毛片大全视频| 麻豆乱淫一区二区| 男女床上黄色一级片免费看| 99香蕉大伊视频| 国产99久久九九免费精品| 韩国高清视频一区二区三区| 国产在线视频一区二区| 国产精品秋霞免费鲁丝片| 啦啦啦视频在线资源免费观看| 老汉色∧v一级毛片| 亚洲av美国av| 99热国产这里只有精品6| 亚洲,一卡二卡三卡| 999精品在线视频| 精品人妻熟女毛片av久久网站| 后天国语完整版免费观看| 国产精品免费大片| 亚洲国产精品一区三区| 精品国产乱码久久久久久男人| 久久久久久久大尺度免费视频| 美女中出高潮动态图| 亚洲,欧美精品.| h视频一区二区三区| 久久天堂一区二区三区四区| 成人免费观看视频高清| 免费女性裸体啪啪无遮挡网站| 1024视频免费在线观看| 老司机午夜十八禁免费视频| 欧美大码av| 欧美性长视频在线观看| 国产亚洲av高清不卡| 少妇粗大呻吟视频| 一区二区日韩欧美中文字幕| 亚洲成色77777| a级毛片在线看网站| 午夜老司机福利片| 欧美精品av麻豆av| 婷婷色av中文字幕| 人人妻人人澡人人爽人人夜夜| 又大又黄又爽视频免费| 老司机影院成人| 精品一区二区三卡| 又大又爽又粗| av在线播放精品| 久久久久久久大尺度免费视频| 五月开心婷婷网| 99久久精品国产亚洲精品| 99九九在线精品视频| 极品少妇高潮喷水抽搐| 1024视频免费在线观看| 国产成人a∨麻豆精品| 女人高潮潮喷娇喘18禁视频| 在线看a的网站| 欧美少妇被猛烈插入视频| 欧美成人午夜精品| 亚洲人成77777在线视频| 国产欧美日韩综合在线一区二区| 久久久久精品国产欧美久久久 | 一本大道久久a久久精品| 亚洲国产毛片av蜜桃av| 在线观看免费日韩欧美大片| 男男h啪啪无遮挡| 精品人妻1区二区| 成在线人永久免费视频| 多毛熟女@视频| 日韩一区二区三区影片| 成在线人永久免费视频| 国产精品免费视频内射| 十八禁网站网址无遮挡| 色婷婷av一区二区三区视频| 国产黄色视频一区二区在线观看| 久久久久精品人妻al黑| 老司机午夜十八禁免费视频| 一级a爱视频在线免费观看| av视频免费观看在线观看| 久久性视频一级片| 99久久99久久久精品蜜桃| 在线观看一区二区三区激情| 日韩av不卡免费在线播放| 久久综合国产亚洲精品| 久久久久久人人人人人| 男女之事视频高清在线观看 | 十八禁高潮呻吟视频| 中文字幕人妻丝袜制服| 男女高潮啪啪啪动态图| 亚洲成人免费电影在线观看 | 欧美 亚洲 国产 日韩一| 搡老岳熟女国产| 满18在线观看网站| 看免费av毛片| 美女午夜性视频免费| 伊人亚洲综合成人网| 男女高潮啪啪啪动态图| 黑人巨大精品欧美一区二区蜜桃| 黄频高清免费视频| av天堂久久9| 成年av动漫网址| 男女之事视频高清在线观看 | 亚洲第一av免费看| 如日韩欧美国产精品一区二区三区| 免费一级毛片在线播放高清视频 | 免费看av在线观看网站| videosex国产| 国产精品免费大片| 美女视频免费永久观看网站| 亚洲九九香蕉| 中国国产av一级| 丝袜在线中文字幕| 看十八女毛片水多多多| 午夜久久久在线观看| 秋霞在线观看毛片| 99精国产麻豆久久婷婷| 丰满人妻熟妇乱又伦精品不卡| 日韩中文字幕欧美一区二区 | 曰老女人黄片| 大香蕉久久网| 午夜久久久在线观看| 婷婷成人精品国产| 大话2 男鬼变身卡| 国产成人免费无遮挡视频| 免费少妇av软件| 麻豆av在线久日| 日日摸夜夜添夜夜爱| 丁香六月天网| 亚洲美女黄色视频免费看| 天天躁日日躁夜夜躁夜夜| 午夜福利乱码中文字幕| 亚洲国产av影院在线观看| 亚洲午夜精品一区,二区,三区| 久久亚洲国产成人精品v| 一区二区av电影网| 美女国产高潮福利片在线看| 国产三级黄色录像| 亚洲国产中文字幕在线视频| 亚洲图色成人| 国产精品国产三级专区第一集| 丝袜美足系列| 欧美老熟妇乱子伦牲交| 婷婷丁香在线五月| 婷婷成人精品国产| 日韩av不卡免费在线播放| 免费一级毛片在线播放高清视频 | 欧美在线黄色| 一区二区av电影网| 欧美亚洲日本最大视频资源| 一级片免费观看大全| 久久 成人 亚洲| 欧美变态另类bdsm刘玥| 久久精品久久久久久噜噜老黄| 日本欧美国产在线视频| www.av在线官网国产| 中文字幕色久视频| 成年av动漫网址| 亚洲精品自拍成人| 国产在视频线精品| 国产在线视频一区二区| 夫妻午夜视频| 亚洲成人免费av在线播放| 亚洲激情五月婷婷啪啪| 精品久久久精品久久久| 十八禁网站网址无遮挡| 久久久久精品国产欧美久久久 | 考比视频在线观看| 亚洲av成人精品一二三区| 下体分泌物呈黄色| 中文字幕另类日韩欧美亚洲嫩草| 青青草视频在线视频观看| videos熟女内射| 国产精品偷伦视频观看了| 欧美日韩国产mv在线观看视频| 国产精品久久久人人做人人爽| 欧美日韩一级在线毛片| 国产日韩欧美视频二区| 美女大奶头黄色视频| 中文字幕高清在线视频| 久久亚洲国产成人精品v| 日本vs欧美在线观看视频| 久久国产精品影院| 国产福利在线免费观看视频| 色播在线永久视频| 久久精品国产亚洲av涩爱| 少妇人妻久久综合中文| 午夜视频精品福利| 最新的欧美精品一区二区| 我要看黄色一级片免费的| 亚洲第一青青草原| 夫妻午夜视频| av在线app专区| 后天国语完整版免费观看| 中文乱码字字幕精品一区二区三区| 成人免费观看视频高清| 黄色怎么调成土黄色| 一区福利在线观看| 亚洲免费av在线视频| 久久精品久久久久久久性| 亚洲国产精品一区二区三区在线| 国产免费一区二区三区四区乱码| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产一区二区精华液| 亚洲国产成人一精品久久久| 午夜av观看不卡| 在线观看免费高清a一片| 熟女av电影| 一边亲一边摸免费视频| 在线观看人妻少妇| 99香蕉大伊视频| 日韩一区二区三区影片| 亚洲国产欧美一区二区综合| 婷婷色综合www| 国产欧美亚洲国产| 搡老乐熟女国产| av在线老鸭窝| 国产成人一区二区三区免费视频网站 | 婷婷色av中文字幕| 十八禁人妻一区二区| 一本综合久久免费| 亚洲av电影在线进入| 日本av手机在线免费观看| 日韩欧美一区视频在线观看| 99热网站在线观看| 久久久久久久大尺度免费视频| 一级,二级,三级黄色视频| 午夜91福利影院| 久久久久精品人妻al黑| 欧美国产精品一级二级三级| 美女中出高潮动态图| 日韩一卡2卡3卡4卡2021年| 久久久久久免费高清国产稀缺| 大话2 男鬼变身卡| 欧美日韩国产mv在线观看视频| 人人妻人人爽人人添夜夜欢视频| 欧美少妇被猛烈插入视频| 涩涩av久久男人的天堂| 中文字幕人妻熟女乱码| 久久久久久免费高清国产稀缺| 亚洲一码二码三码区别大吗| 亚洲精品国产av成人精品| 国产成人精品在线电影| 脱女人内裤的视频| 99国产精品免费福利视频| √禁漫天堂资源中文www| 久久久国产欧美日韩av| 操美女的视频在线观看| 亚洲图色成人| 午夜老司机福利片| 成人国产一区最新在线观看 | 少妇精品久久久久久久| 十分钟在线观看高清视频www| 午夜精品国产一区二区电影| 中文字幕人妻熟女乱码| 性色av一级| 亚洲av欧美aⅴ国产| 一本大道久久a久久精品| 日韩熟女老妇一区二区性免费视频| 青青草视频在线视频观看| 国产精品久久久人人做人人爽| 性少妇av在线| 三上悠亚av全集在线观看| 一级,二级,三级黄色视频| 老司机亚洲免费影院| 免费高清在线观看日韩| 男男h啪啪无遮挡| 国产欧美日韩精品亚洲av| 男女国产视频网站| 国产精品久久久人人做人人爽| 亚洲精品乱久久久久久| 欧美av亚洲av综合av国产av| 亚洲成国产人片在线观看| 国产深夜福利视频在线观看| 精品国产一区二区久久| 亚洲专区中文字幕在线| 欧美人与性动交α欧美精品济南到| 一二三四在线观看免费中文在| 欧美激情高清一区二区三区| 91九色精品人成在线观看| 一二三四社区在线视频社区8| 免费久久久久久久精品成人欧美视频| 巨乳人妻的诱惑在线观看| 国产成人91sexporn| 午夜视频精品福利| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇粗大呻吟视频| 一区二区av电影网| 精品亚洲乱码少妇综合久久| 精品国产一区二区久久| 久久人人97超碰香蕉20202| 中文字幕人妻熟女乱码| 三上悠亚av全集在线观看| 欧美日韩福利视频一区二区| 蜜桃国产av成人99| 99精品久久久久人妻精品| 国产一级毛片在线| 一区二区三区精品91| 男人爽女人下面视频在线观看| 黑丝袜美女国产一区| 国产1区2区3区精品| 欧美激情 高清一区二区三区| 国产亚洲欧美精品永久| 午夜免费观看性视频| 2018国产大陆天天弄谢| 丝袜喷水一区| 麻豆国产av国片精品| 一区二区三区激情视频| 国产成人精品久久二区二区91| 午夜91福利影院| 啦啦啦 在线观看视频| 国产人伦9x9x在线观看| 考比视频在线观看| 久久影院123| 久久精品aⅴ一区二区三区四区| 国产精品人妻久久久影院| 999久久久国产精品视频| 日韩大码丰满熟妇| 老司机深夜福利视频在线观看 | 老鸭窝网址在线观看| 又黄又粗又硬又大视频| 夫妻性生交免费视频一级片| 国产99久久九九免费精品| 18禁国产床啪视频网站| 国产精品 欧美亚洲| 97在线人人人人妻| 菩萨蛮人人尽说江南好唐韦庄| 另类精品久久| 男人操女人黄网站| 久久久久国产精品人妻一区二区| 男的添女的下面高潮视频| 欧美人与性动交α欧美软件| 久久人妻福利社区极品人妻图片 | 国产亚洲欧美在线一区二区| av国产精品久久久久影院| 欧美变态另类bdsm刘玥| √禁漫天堂资源中文www| 精品人妻在线不人妻| 色综合欧美亚洲国产小说| 午夜免费成人在线视频| 久久这里只有精品19| 久久久久久久大尺度免费视频| 国产熟女午夜一区二区三区| 欧美日韩亚洲高清精品| 成人手机av| 亚洲精品成人av观看孕妇| 男的添女的下面高潮视频| 免费日韩欧美在线观看| 人成视频在线观看免费观看| 亚洲午夜精品一区,二区,三区| 丝袜美腿诱惑在线| 午夜91福利影院| 19禁男女啪啪无遮挡网站| 人人妻人人添人人爽欧美一区卜| 欧美人与善性xxx| 可以免费在线观看a视频的电影网站| 日韩制服丝袜自拍偷拍| 成年人黄色毛片网站| 侵犯人妻中文字幕一二三四区| 高清黄色对白视频在线免费看| 亚洲国产欧美网| 国产欧美亚洲国产| 啦啦啦中文免费视频观看日本| 国产在线视频一区二区| 久久精品久久精品一区二区三区| 亚洲,欧美,日韩| 国产视频首页在线观看| 免费在线观看完整版高清| 中国国产av一级| 赤兔流量卡办理| 精品第一国产精品| 国产91精品成人一区二区三区 | 亚洲精品美女久久av网站| 久久国产精品人妻蜜桃| 亚洲精品美女久久久久99蜜臀 | 一二三四社区在线视频社区8| 9色porny在线观看| 热99久久久久精品小说推荐| 夫妻午夜视频| 久久久精品免费免费高清| 久久久精品94久久精品| 久久亚洲国产成人精品v| 日本av手机在线免费观看| 国产成人精品在线电影| 欧美精品啪啪一区二区三区 | 免费观看a级毛片全部| 国产高清视频在线播放一区 | 亚洲国产最新在线播放| 永久免费av网站大全| 九色亚洲精品在线播放| 国产一区二区三区av在线| 午夜精品国产一区二区电影| 首页视频小说图片口味搜索 | 国产高清国产精品国产三级| 在线天堂中文资源库| 99久久99久久久精品蜜桃| 波野结衣二区三区在线| www日本在线高清视频| 1024视频免费在线观看| 晚上一个人看的免费电影| 菩萨蛮人人尽说江南好唐韦庄| 丁香六月天网| 欧美另类一区| 亚洲成av片中文字幕在线观看| 大片免费播放器 马上看| 美女脱内裤让男人舔精品视频| 亚洲一区中文字幕在线| 日韩av不卡免费在线播放| 国产精品一二三区在线看| 欧美人与性动交α欧美精品济南到| 999精品在线视频| 一区在线观看完整版| 亚洲伊人色综图| 中文精品一卡2卡3卡4更新| 国产伦理片在线播放av一区| 人人妻人人添人人爽欧美一区卜| 午夜福利免费观看在线| 婷婷色综合www| 色精品久久人妻99蜜桃| 亚洲国产欧美日韩在线播放| 国产成人91sexporn| 亚洲第一av免费看| 在线av久久热| 国产av一区二区精品久久| 最新在线观看一区二区三区 | cao死你这个sao货| 免费看av在线观看网站| 精品久久久久久电影网| 精品国产一区二区三区四区第35| 国产一区二区 视频在线| 国产三级黄色录像| 欧美人与性动交α欧美软件| 男女边摸边吃奶| 妹子高潮喷水视频| 精品熟女少妇八av免费久了| 亚洲精品一区蜜桃| 又粗又硬又长又爽又黄的视频| 熟女av电影| 亚洲成人手机| 麻豆av在线久日| 在线亚洲精品国产二区图片欧美| 日韩大片免费观看网站| 国产av一区二区精品久久| 人人妻人人添人人爽欧美一区卜| 久久国产精品大桥未久av| 成年人午夜在线观看视频| 欧美亚洲 丝袜 人妻 在线| 久久人人97超碰香蕉20202| 2018国产大陆天天弄谢| 一本一本久久a久久精品综合妖精| 久久国产亚洲av麻豆专区| 肉色欧美久久久久久久蜜桃| 又紧又爽又黄一区二区| 在线观看免费午夜福利视频| 男女午夜视频在线观看| 免费av中文字幕在线| 精品少妇黑人巨大在线播放| 老司机深夜福利视频在线观看 | 9色porny在线观看| 99国产综合亚洲精品| 97人妻天天添夜夜摸| 在线观看免费日韩欧美大片| 新久久久久国产一级毛片| 国产精品免费大片| 欧美人与善性xxx| 亚洲精品国产色婷婷电影| 欧美黑人精品巨大| 人人妻人人添人人爽欧美一区卜| 精品久久久久久久毛片微露脸 | 午夜免费鲁丝| 欧美xxⅹ黑人| 亚洲精品中文字幕在线视频| 777米奇影视久久| 在线精品无人区一区二区三| 一区二区三区乱码不卡18| 又大又爽又粗| 人人妻,人人澡人人爽秒播 | 国产xxxxx性猛交| 国产成人av激情在线播放| 色网站视频免费| 秋霞在线观看毛片| 午夜老司机福利片| 久久亚洲精品不卡| 黄色怎么调成土黄色| 亚洲九九香蕉| 久久久精品区二区三区| 日韩大片免费观看网站| 中文字幕人妻丝袜一区二区| netflix在线观看网站| 国产精品 国内视频| 国产高清视频在线播放一区 | 日日夜夜操网爽| 国产伦人伦偷精品视频| 又紧又爽又黄一区二区| 国产国语露脸激情在线看| 日韩欧美一区视频在线观看| 国产免费一区二区三区四区乱码| 久久免费观看电影| 少妇粗大呻吟视频| 美女扒开内裤让男人捅视频| 国产精品二区激情视频| netflix在线观看网站| 亚洲激情五月婷婷啪啪| 国产野战对白在线观看| 久久人人爽av亚洲精品天堂| 性色av乱码一区二区三区2| 一级黄片播放器| bbb黄色大片| 国产真人三级小视频在线观看| 国产成人av激情在线播放| 香蕉丝袜av| 国产亚洲一区二区精品| 两人在一起打扑克的视频| 国产伦理片在线播放av一区| 欧美国产精品va在线观看不卡| 一级,二级,三级黄色视频| 成人国产一区最新在线观看 | 另类亚洲欧美激情| 欧美精品av麻豆av| 久久久欧美国产精品| 国产男女内射视频| 久久久国产欧美日韩av| 亚洲精品一二三| 1024视频免费在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产成人一区二区三区免费视频网站 | 制服人妻中文乱码| 国产午夜精品一二区理论片| 亚洲精品久久午夜乱码| 精品国产一区二区三区久久久樱花| 中文字幕精品免费在线观看视频| 日韩人妻精品一区2区三区| 黄色视频在线播放观看不卡| 国产免费福利视频在线观看| 精品国产超薄肉色丝袜足j| 成人黄色视频免费在线看| 国产欧美日韩一区二区三区在线| 一级毛片 在线播放| 在线观看人妻少妇| 国产精品人妻久久久影院| 无遮挡黄片免费观看| 免费在线观看黄色视频的| 天堂俺去俺来也www色官网| 国产福利在线免费观看视频| 精品久久久精品久久久| 51午夜福利影视在线观看| 亚洲国产精品成人久久小说| 麻豆国产av国片精品| 亚洲国产看品久久| 欧美日韩视频高清一区二区三区二| 人妻一区二区av|