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

    ·OH自由基介導2,4-二氯酚光降解機理及產(chǎn)物毒性變化特征的理論研究

    2016-12-12 01:40:20高艷蓬姬越蒙李桂英安太成
    生態(tài)毒理學報 2016年2期
    關鍵詞:光降解中間體毒性

    高艷蓬,姬越蒙,李桂英,安太成,,

    1. 中國科學院廣州地球化學研究所 有機地球化學國家重點實驗室和廣東省環(huán)境資源利用與保護重點實驗室,廣州 510640 2. 廣東工業(yè)大學 環(huán)境健康與污染控制研究院和環(huán)境科學與工程學院,廣州 510006

    ?

    ·OH自由基介導2,4-二氯酚光降解機理及產(chǎn)物毒性變化特征的理論研究

    高艷蓬1,姬越蒙2,李桂英2,安太成1,2,

    1. 中國科學院廣州地球化學研究所 有機地球化學國家重點實驗室和廣東省環(huán)境資源利用與保護重點實驗室,廣州 510640 2. 廣東工業(yè)大學 環(huán)境健康與污染控制研究院和環(huán)境科學與工程學院,廣州 510006

    由于2,4-二氯酚(2,4-dichlorophenol, 簡稱2,4-DCP)具有“三致”作用及潛在的生態(tài)風險,其光化學降解機理及產(chǎn)物的生態(tài)毒理效應目前成為人們關注的焦點。本研究采用理論化學計算手段,系統(tǒng)闡明了2,4-DCP在水環(huán)境中·OH介導的間接光化學轉化機理、動力學和轉化產(chǎn)物的生態(tài)毒性變化特征。結果表明2,4-DCP很容易被·OH氧化降解,其降解主要通過·OH-加成和H-遷移路徑進行。在低溫條件下·OH-加成路徑將占主導,主要形成4,6-二氯苯-1,3-二酚;而當溫度超過313 K時,·OH提取酚羥基上的氫原子為主要降解途徑,主要形成2-氯苯-對苯醌。尤其在活性物種濃度較低的高溫環(huán)境中,H-遷移路徑有可能生成二噁英。計算毒理學結果表明:H-遷移產(chǎn)物的水生毒性超過·OH-加成產(chǎn)物,甚至超出母體2,4-DCP一個毒性等級。因此,我們建議在以后的環(huán)境監(jiān)測以及風險評估過程中,關注2,4-DCP及其轉化產(chǎn)物特別是H-遷移路徑的轉化產(chǎn)物的生態(tài)毒理學問題。

    2,4-二氯酚;光降解機理;轉化產(chǎn)物;理論化學計算;生態(tài)毒性;毒性變化特征

    氯酚類化合物(chlorophenols, 簡稱CPs)是一類典型的難降解毒害有機污染物。CPs作為一類重要的化工原料,已被廣泛地用于農(nóng)藥、殺蟲劑和防腐劑生產(chǎn)以及染料等工業(yè)[1]。隨著大量使用和生產(chǎn)泄露,這類化合物已成為目前煉油、造紙、塑料化工等工業(yè)廢水以及飲用水氯化消毒過程中的主要污染物,在河流和土壤中頻繁檢出。調(diào)查研究表明2,4-二氯酚(2,4-dichlorophenol, 簡稱2,4-DCP)是水環(huán)境中分布最為廣泛、含量最為豐富的CPs[2]。在中國眾多流域和土壤中均檢測到2,4-DCP[3],部分水體中2,4-DCP的濃度甚至高達0.1~1 mg·L-1[4]。同時相關研究發(fā)現(xiàn)2,4-DCP具有潛在的“致癌、致畸、致突變”三致作用[5],能夠通過食物鏈在動物體內(nèi)不斷積累,干擾人體的內(nèi)分泌系統(tǒng)[6]。因此,2,4-DCP已被列為“中國環(huán)境優(yōu)先控制污染物黑名單”以及美國環(huán)保局的“環(huán)境優(yōu)控污染物”。

    在自然水體中,光化學過程是有機污染物轉化的重要過程[7,8]。自然水體中的有機質(zhì)、硝酸鹽或亞硝酸鹽、Fe(III)等經(jīng)光敏化或光化學作用等能產(chǎn)生大量的活性物種如·OH、1O2和H2O2等[9-11]。而在這些活性物種中,·OH尤其具有極強的氧化活性,能夠介導氧化降解自然水體中絕大部分的有機污染物,進而對其遷移轉化和環(huán)境歸趨具有非常重要的作用[11-12]。因此·OH介導的光化學轉化過程在一定程度上決定了環(huán)境水體中有機污染物的環(huán)境命運。但是,近期研究發(fā)現(xiàn)該反應過程中有可能形成仍具有一定毒性的轉化中間產(chǎn)物[13-15]。例如本課題組前期的研究發(fā)現(xiàn)防腐劑三氯生在·OH介導的間接光化學轉化過程中生成了致癌物質(zhì)2,4-DCP[16]。然而目前的研究大多集中在探索提高水環(huán)境中2,4-DCP的光降解效率[17],而對于其光降解機理研究不多見。盡管有研究檢測到2,4-DCP在光降解過程中生成了不同的轉化產(chǎn)物[18],但是對這些轉化產(chǎn)物的形成機理及生態(tài)毒性的關注目前還較少。

    理論化學計算作為一種新的研究手段越來越多地應用于環(huán)境科學領域[15-16,19-21]。它不僅可以從分子尺度進行考慮,提供一些實驗研究中無法觀測到的重要反應機理的信息如過渡態(tài)和反應能壘等,進而得到有機污染物氧化降解的微觀機理與環(huán)境歸趨等;還發(fā)揮著彌補基礎數(shù)據(jù)的缺失、降低測試費用、減少動物實驗等重要作用[22],從理論上輔助解決當前實驗研究有機污染物光化學轉化機理與毒性風險評估的難題。此外,隨著大量有機污染物在環(huán)境中的不斷涌現(xiàn),如何快速準確地了解這些有機污染物在天然水環(huán)境中的氧化降解機制、動力學和降解產(chǎn)物等環(huán)境問題就顯得尤為重要。因此,本研究主要采用理論計算手段,研究了2,4-DCP的·OH介導間接光化學降解過程與機理、動力學及其產(chǎn)物毒性變化特征,利用計算化學手段從分子層面闡明該污染物的光降解產(chǎn)物與環(huán)境歸趨,這將有助于了解該類污染物的環(huán)境轉化過程,正確評估其污染狀況及其健康風險。

    1 理論計算方法(Computational methods)

    本文所涉及的量子化學計算采用Gaussian 09程序完成[23]。采用B3LYP方法,以6-31+G(d,p)為基組,對反應物、產(chǎn)物和過渡態(tài)的構型進行了幾何參數(shù)全優(yōu)化。并在相同水平下計算了反應過程中所有駐點的振動頻率以確定過渡態(tài),確定過渡態(tài)有且只有一個虛頻,并采用內(nèi)稟坐標(IRC)方法確認了過渡態(tài)的真實性。溶劑化的計算是在B3LYP/6-31+G(d,p)水平上優(yōu)化得到結構,采用Gaussian 09推薦的SMD虛擬溶劑模型[24]在B3LYP/6-311++G(3df,2p)水平下計算得到單點能及溶劑化能。反應速率常數(shù)采用過渡態(tài)理論模型,并考慮籠效應、擴散、標準態(tài)校正等情況下計算得到。

    本研究選取魚為典型的水生物種[25],用以評估2,4-DCP及其轉化產(chǎn)物的水生毒性。2,4-DCP及其轉化產(chǎn)物的急性毒性和慢性毒性將采用美國國家環(huán)保局的美國國家環(huán)保局開發(fā)的“生態(tài)-結構活性關系 (ECOSAR)”[26]計算得到。急性毒性是通過魚暴露96 h的半致死濃度LC50來表示,而慢性毒性值(ChV)為無觀察效應濃度(NOECs)和最低可觀測效應濃度(LOEC)的幾何平均值,并且對于所屬多種類別的目標化合物僅列出了最保守的毒性效應值。

    圖1 2,4-DCP和·OH的初始反應路徑圖及其反應焓(ΔH) (單位: kcal·mol-1)Fig. 1 All the possible pathways in the initial reaction of 2,4-DCP with ·OH, as well as the reaction enthalpies (ΔH) (unit: kcal·mol-1)

    2 結果與討論(Results and discussion)

    2.1 ·OH自由基介導的2,4-DCP光化學降解機理

    水體環(huán)境中·OH介導2,4-DCP的光化學轉化可能有3類反應路徑(如圖1所示):(i) ·OH-加成路徑:·OH進攻2,4-DCP的苯環(huán)碳原子(Radd1-6)形成不同的羥基加合物(·DCP-OH),或者脫去氯原子形成4-Cl-鄰苯二酚(CAT)和2-Cl氫醌(HQ);(ii) H-遷移路徑,·OH提取2,4-DCP的羥基氫(RabsO),形成脫氫中間體(RabsO(-H)·);(iii) 單電子轉移路徑:·OH將2,4-DCP上的電子轉移至·OH(Rset)。這些路徑中所涉及到的過渡態(tài)(TS)結構均列于圖2中。計算得到2,4-DCP與·OH反應各路徑的反應焓值(ΔH)和勢能面圖分別列于圖1和圖3中,其中利用反應物(2,4-DCP+·OH)的能量之和為參考零點。

    圖3 在B3LYP/6-311++G(3df,2p)//B3LYP/6-31+G(d,p)水平下,計算得到的2,4-DCP和·OH反應的勢能面圖 (單位: kcal·mol-1)Fig. 3 Schematic free energy diagram for the reaction between 2,4-DCP and ·OH at the B3LYP/6-311++G(3df,2p)//B3LYP/6-31+G(d,p) level (Unit: kcal·mol-1)

    圖2 在B3LYP/6-31+G(d,p)水平下,優(yōu)化得到各路徑的過渡態(tài)(TS)結構,鍵長單位為?Fig. 2 Optimized geometries of transition states of the initial reaction of 2,4-DCP with ·OH at the B3LYP/6-31+G(d,p) level.Bond lengths are in angstroms.

    從圖1列出的ΔH值可以看出,除了路徑Rset為吸熱反應之外(H = 6.91 kcal·mol-1),其他所有·OH-加成和H-遷移路徑均為放熱反應(G = -13.94~-33.01 kcal·mol-1)。這說明在2,4-DCP的·OH介導光降解過程中,單電子轉移路徑與放熱的·OH-加成和H-遷移路徑相比是難以發(fā)生。這與我們以及其他課題組相關的已有污染物電子轉移反應機理的文獻相一致[20-21,27],即·OH與分子形態(tài)化合物發(fā)生電子轉移反應的活性較低,對水體有機污染物的光化學降解過程貢獻幾乎可以忽略。

    從圖3列出的反應勢能面圖可以看出,·OH-加成路徑(Radd5和Radd3)具有較低的能壘,分別是8.84和8.91 kcal·mol-1,比其他·OH-加成路徑(Radd1,2,4,6)低了約1~2 kcal·mol-1。從這一結果初步推測,前者的兩條路徑比其他·OH-加成路徑略微地占一定優(yōu)勢。對于H-遷移路徑而言,路徑RabsO的能壘為8.93 kcal·mol-1,僅比能壘最小的·OH-加成路徑Radd5高出0.09 kcal·mol-1。由此可以看出,各反應路徑之間的能壘差別較小,僅僅從能量的角度無法判斷出主要的反應路徑及其反應中間體,有必要進一步開展反應動力學的計算。

    2.2 ·OH自由基介導的2,4-DCP光化學降解動力學

    為了進一步模擬真實水體環(huán)境中·OH-加成和H-遷移反應路徑對2,4-DCP光化學降解反應的貢獻,更好地理解2,4-DCP在水體環(huán)境中的壽命與環(huán)境歸趨,本文對273~313 K溫度區(qū)間內(nèi)的反應動力學進行了研究。表1分別列出了在該溫度區(qū)內(nèi)所有反應路徑和總反應的速率常數(shù)值(ktotal)。從表中可以看出:總反應速率ktotal值約為109~1010(mol·L-1)-1·s-1,這說明·OH介導的2,4-DCP轉化反應在整個溫度區(qū)間內(nèi)均屬于擴散控制反應過程(< 4109(mol·L-1)-1·s-1)[28]。同時從表1可以看出:所有路徑的反應速率常數(shù)以及總反應速率均隨著溫度的升高而增加。例如,隨著溫度從273 K升高到313 K,總反應速率ktotal值從2.62×109(mol·L-1)-1·s-1逐漸升高到9.82×1010(mol·L-1)-1·s-1。因此,溫度升高將有利于促進·OH-介導2,4-DCP的光化學降解過程。另外,我們計算得到的298 K下·OH與2,4-DCP的反應速率常數(shù)為6.76×109(mol·L-1)-1·s-1。該速率常數(shù)數(shù)值與已有文獻報導的實驗值處于一個數(shù)量級范圍內(nèi)。例如競爭動力學實驗方法測得該反應的速率常數(shù)為5.14×109(mol·L-1)-1·s-1[29]和7.0×109(mol·L-1)-1·s-1[30];輻射實驗技術測得的數(shù)值為(7.4±0.5)×109(mol·L-1)-1·s-1[31]。這些實驗結果均很好地印證了我們計算結果的可靠性。

    另外,本文也給出了273~313 K溫度區(qū)間各反應路徑以及總反應的阿倫尼烏斯方程,具體數(shù)據(jù)列于表2中。從表可以獲得該溫度區(qū)間內(nèi)的指前因子和活化能,以及任意溫度下各反應路徑的速率常數(shù)等。例如在273~313 K溫度區(qū)間內(nèi),計算得到2,4-DCP與·OH反應的活化能為5.63 kcal·mol-1。較小的活化能進一步驗證了在水體環(huán)境中2,4-DCP很容易被·OH氧化降解。因此非常有必要繼續(xù)對2,4-DCP光化學轉化所形成的中間產(chǎn)物進行跟蹤研究,進一步了解其生態(tài)毒性效應變化特征。

    表1 在273~313 K的溫度區(qū)間內(nèi)2,4-DCP與·OH反應的速率常數(shù)值((mol·L-1)-1·s-1)

    2.3 轉化中間體與轉化產(chǎn)物

    為了具體了解各反應路徑對·OH介導2,4-DCP光化學降解的貢獻大小以及定量預測所生成的主要降解中間產(chǎn)物,我們對所有反應路徑的分支比(Γ)與溫度的依賴關系做了進一步計算,結果列于圖4中。從圖中可以看出:·OH-加成路徑Radd3和Radd5在273 K下的分支比分別為31.4%和34.1%,比其他的加成路徑高出24%~32%,同時比H-遷移反應路徑RabsO高出了11%以上。這表明在該溫度條件下,2,4-DCP發(fā)生·OH介導的光化學降解主要是通過·OH-加成路徑Radd3和Radd5兩個途徑進行。然而,隨著溫度的升高,路徑Radd3和Radd5對總反應的貢獻將逐漸降低,當溫度為313 K時H-遷移反應路徑RabsO的分支比增加為22.1%,超過·OH-加成路徑Radd3 (21.2%) 和Radd5 (21.7%)對總反應的貢獻。因此,得到如下結論:2,4-DCP的光化學降解在低溫條件下將主要以·OH-加成路徑Radd3和Radd5為主,主要生成羥基加合物·DCP-OH5和·DCP-OH3。

    從以上研究可以看出,在2,4-DCP的·OH-介導的間接光化學降解的初始反應步驟中,將主要形成·OH-加合物(·DCP-OH3和·DCP-OH5)和脫氫中間體(·DCP(-H))。這些轉化中間體均為自由基,具有較高的反應活性,可以進一步發(fā)生一系列的后續(xù)反應,進而形成轉化降解產(chǎn)物。因此,本研究將分別以主要的·OH-加合物(·DCP-OH5)和脫氫中間體(·DCP(-H))為代表,進一步系統(tǒng)了解這些轉化中間體在水環(huán)境中的后續(xù)反應。

    圖5列出了主要中間體·DCP-OH5的后續(xù)反應路徑圖。從圖中可以看出·DCP-OH5很容易與·OH進一步反應形成二羥基化合物(DCP-OH5-OH4),而該過程為無能壘的放熱反應,放熱量為52.15 kcal·mol-1。這些能量將有助于克服DCP-OH5-OH4的進一步脫水反應的能壘(56.42 kcal·mol-1),形成轉化產(chǎn)物4,6-二氯苯-1,3-二酚(PC1)。在實際水體環(huán)境中,O2也是一種非常重要的氧化劑,而且有可能參與主要中間體的后續(xù)反應,因此我們也考慮了O2參與轉化中間體的后續(xù)反應。圖5列出的中間體·DCP-OH5與O2后續(xù)反應的路徑圖可以看出,中間體·DCP-OH5容易與O2發(fā)生反應,只需經(jīng)過能量為6.90 kcal·mol-1的過渡態(tài)從而形成轉化產(chǎn)物PC1。在以前的實驗研究中也有檢測到該產(chǎn)物PC1[18],再次證實了我們理論計算機理的可靠性。而且該步驟的能壘比初始步驟(Radd5)的能壘還低了1.94 kcal·mol-1,說明在2,4-DCP的·OH-介導光化學降解反應過程中初始步為整個反應的決速步驟。綜上所述,我們認為不論是在·OH充裕的水環(huán)境還是有氧條件下,PC1都將是·OH-加成反應路徑上的主要轉化產(chǎn)物。

    表2 在273~313 K的溫度區(qū)間內(nèi)2,4-DCP各反應路徑以及總反應的阿倫尼烏斯方程

    圖4 ·OH介導2,4-DCP光降解過程中各路徑對總反應的貢獻(Γ%)隨溫度的變化情況Fig. 4 Calculated branching ratios (Γ) of the main transformation pathways for the reaction of 2,4-DCP with ·OH within the temperature range of 273 to 313 K

    圖5 主要中間體·DCP-OH5的后續(xù)反應勢路徑圖 (單位: kcal·mol-1)Fig. 5 Schematic diagram of the subsequent pathways of ·DCP-OH5 intermediates (Unit: kcal·mol-1)

    圖6 脫氫中間體(·DCP(-H))的后續(xù)反應勢路徑圖 (單位: kcal·mol-1)Fig. 6 Schematic diagram of the subsequent pathways of ·DCP(-H) (Unit: kcal·mol-1)

    圖6列出了脫氫中間體(·DCP(-H))的后續(xù)反應路徑圖。從圖中可以看出·DCP(-H)很容易與·OH結合形成羥基化合物(·DCP(-H)-OH4),而且該過程是無能壘的放熱反應,放熱量為43.06 kcal·mol-1。該能量將完全能夠克服DCP-OH5-OH4的后續(xù)反應能壘(15.10 kcal·mol-1),形成轉化產(chǎn)物2-氯苯-對苯醌(PC2)。然而在有氧環(huán)境中,O2一般也會參與脫氫中間體(·DCP(-H))的后續(xù)反應。從圖6中可以看出,中間體·DCP-OH5與O2發(fā)生反應,需經(jīng)過能量為11.83 kcal·mol-1的過渡態(tài)形成轉化過氧化物(DCP(-H)-OH4),且該過程為吸熱過程(8.83 kcal·mol-1)。此外過氧化物脫去一分子·OCl形成轉化產(chǎn)物PC2的過程中,需要較高的能壘20.85 kcal·mol-1。因此O2介導的脫氫中間體(·DCP(-H))的后續(xù)轉化路徑將不如前者·OH參與的反應容易進行。即在·OH濃度較高的的水環(huán)境中,PC2將可能是H-遷移反應路徑上的主要轉化產(chǎn)物。然而,在活性物種濃度比較低的水體環(huán)境中,脫氫中間體(·DCP(-H))將有可能進行自身的雙分子二聚反應。前人的文獻研究發(fā)現(xiàn),該過程有可能會進一步轉化生成二噁英,例如多氯代二噁英(PCDD)和多氯代二苯并呋喃(PCDF)等[32-33]。

    表3 2,4-DCP及其光降解產(chǎn)物對魚類的水生毒性(mg·L-1)

    2.4 2,4-DCP及其主要光降解產(chǎn)物的水生毒性研究

    美國國家環(huán)保局開發(fā)的ECOSAR程序在有機污染物的毒性預測方面表現(xiàn)出了較好的優(yōu)勢,可以成功預測了多種化學品及其轉化產(chǎn)物的生態(tài)毒性[34]。因此,本文采用該程序預測了2,4-DCP及其光降解產(chǎn)物對魚類的急性和慢性毒性,結果列在了表3中。表中的毒性效應濃度值越小,則說明相應的水生毒性越大。從表中可知:2,4-DCP的急性和慢性毒性的效應濃度分別為6.96和0.79 mg·L-1。

    根據(jù)歐盟標準(described in Annex VI of Directive 67/548/EEC) 和中國新化學物質(zhì)風險評估方案(HJ/T 154-2004)可知,2,4-DCP對魚的急性和慢性毒性均屬于有毒等級(1< LC50<10和0.1< ChV <1 mg·L-1)。由此可見2,4-DCP是一種毒性污染物,會對魚類產(chǎn)生一定的危害,因此有必要進一步關注其轉化降解產(chǎn)物對水生生物的危害作用。

    表3分別列出了2,4-DCP分別經(jīng)·OH-加成路徑Radd5和H-遷移反應路徑RabsO所形成轉化產(chǎn)物的水生毒性。從中可以看出:·OH-加成產(chǎn)物PC1和H-遷移產(chǎn)物PC2的LC50值分別為6.17和0.10 mg·L-1,均低于母體2,4-DCP的LC50值(6.96 mg·L-1),尤其是H-遷移產(chǎn)物的數(shù)值比母體2,4-DCP低了將近70倍。這一結果表明:2,4-DCP在·OH介導的光降解過程中,不論經(jīng)由·OH-加成路徑還是H-遷移路徑,均使得其水生毒性有所增加,特別是其H-遷移過程的毒性增加更為顯著。而對于慢性毒性而言,·OH-加成產(chǎn)物的ChV值為3.23 mg·L-1,比母體的ChV值高出4倍,這表明經(jīng)·OH-加成路徑轉化過程時水生毒性有所減小。依據(jù)毒性分類,·OH-加成產(chǎn)物PC1歸屬為有害等級,比母體污染物2,4-DCP低了一個毒性等級。然而H-遷移產(chǎn)物卻呈現(xiàn)出相反的趨勢,即H-遷移產(chǎn)物PC2的ChV值僅為0.01 mg·L-1,比母體2,4-DCP高出了一個毒性等級,歸屬為非常毒的等級。

    綜上所述,2,4-DCP在水環(huán)境中·OH-介導的間接光降解過程對其環(huán)境歸趨起著重要的作用。計算結果表明該過程的活化能僅為5.63 kcal·mol-1,說明該過程在實際水體中比較容易發(fā)生。在298 K下反應速率常數(shù)為6.76×109(mol·L-1)-1·s-1,這與實驗結果(5.14×109~7.0×109) (mol·L-1)-1·s-1符合的較好,說明我們本文計算結果具有一定的可靠性。值得指出的是在這類污染物的環(huán)境監(jiān)測與風險評估過程中,除了關注這類污染物本身之外,也應當重視這些光降解產(chǎn)物,特別是其H-遷移轉化產(chǎn)物的水生生態(tài)風險值得重點關注。

    [1] 羅義, 紀靚靚, 蘇燕, 等. 2,4-二氯苯酚誘導鯽魚活性氧(ROS)的產(chǎn)生及其分子致毒機制[J]. 環(huán)境科學學報, 2007, 27(1): 129-134

    Luo Y, Ji L L, Su Y, et al. Reactive oxygen species generation and the molecular toxic mechanism induced by 2,4- dichlorophenol in Carassius auratus [J]. Acta Scientiae Circumstantiae, 2007, 27(1): 129-134 (in Chinese)

    [2] House W A, Leach D, Long J L A, et al. Micro-organic compounds in the Humber rivers [J]. Science of the TotalEnvironment, 1997, 194: 357-371

    [3] Gao J J, Liu L H, Liu X R, et al. Levels and spatial distribution of chlorophenols 2,4-dichlorophenol, 2,4,6-trichlorophenol, and pentachlorophenol in surface water of China [J]. Chemosphere, 2008, 71(6): 1181-1187

    [4] 李碧, 陳軼, 孔殿超. 2,4-二氯苯酚降解研究進展[J]. 廣州化工, 2013, 41(5): 40-42

    Li B, Chen T, Kong D C. Advance research on the degradation of 2,4-dichlorophenol [J]. Guangzhou Chemical Industry, 2013, 41(5): 40-42 (in Chinese)

    [5] 李學德, 岳永德, 花日茂, 等. 水中2,4-二氯苯酚的光催化降解研究[J]. 農(nóng)業(yè)環(huán)境保護, 2002, 21(2): 156-158, 168

    Li X D, Yue Y D, Hua R M, et al. Photocatalytic degradation of 2,4 - dicholorophenal in aqueous solution [J]. Agro-Enviromental Protection, 2002, 21(2): 156-158, 168 (in Chinese)

    [6] Czaplicka M. Sources and transformations of chlorophenols in the natural environment [J]. Science of the Total Environment, 2004, 322(1-3): 21-39

    [7] Buerge I J, Buser H R, Muller M D, et al. Behavior of the polycyclic musks HHCB and AHTN in lakes, two potential anthropogenic markers for domestic wastewater in surface waters [J]. Environmental Science & Technology, 2003, 37(24): 5636-5644

    [8] Sanchez-Prado L, Lourido M, Lores M, et al. Study of the photoinduced degradation of polycyclic musk compounds by solid-phase microextraction and gas chromatography/mass spectrometry [J]. Rapid Communications in Mass Spectrometry, 2004, 18(11): 1186-1192

    [9] Fang G D, Zhu C Y, Dionysiou D D, et al. Mechanism of hydroxyl radical generation from biochar suspensions: Implications to diethyl phthalate degradation [J]. Bioresource Technology, 2015, 176: 210-217

    [10] Vione D, Falletti G, Maurino V, et al. Sources and sinks of hydroxyl radicals upon irradiation of natural water samples [J]. Environmental Science & Technology, 2006, 40(12): 3775-3781

    [11] Dong M M, Rosario-Ortiz F L. Photochemical formation of hydroxyl radical from effluent organic matter [J]. Environmental Science & Technology, 2012, 46(7): 3788-3794

    [12] Wenk J, von Gunten U,Canonica S. Effect of dissolved organic matter on the transformation of contaminants induced by excited triplet states and the hydroxyl radical [J]. Environmental Science & Technology, 2011, 45(4): 1334-1340

    [13] Gao Y, An T, Ji Y, et al. Eco-toxicity and human estrogenic exposure risks from OH-initiated photochemical transformation of four phthalates in water: A computational study [J]. Environmental Pollution, 2015, 206: 510-517

    [14] Gao Y, Ji Y, Li G, et al. Mechanism, kinetics and toxicity assessment of OH-initiated transformation of triclosan in aquatic environments [J]. Water Research, 2014, 49: 360-370

    [15] Gao Y P, An T C, Fang H S, et al. Computational consideration on advanced oxidation degradation of phenolic preservative, methylparaben, in water: Mechanisms, kinetics, and toxicity assessments [J]. Journal of Hazardous Materials, 2014, 278: 417-425

    [16] Gao Y P, Ji Y M, Li G Y, et al. Mechanism, kinetics and toxicity assessment of OH-initiated transformation of triclosan in aquatic environments [J]. Water Research, 2014, 49: 360-370

    [17] Karci A, Arslan-Alaton I, Olmez-Hanci T, et al. Degradation and detoxification of industrially important phenol derivatives in water by direct UV-C photolysis and H2O2/UV-C process: A comparative study [J]. Chemical Engineering Journal, 2013, 224: 4-9

    [18] Hirvonen A, Trapido M, Hentunen J, et al. Formation of hydroxylated and dimeric intermediates during oxidation of chlorinated phenols in aqueous solution [J]. Chemosphere, 2000, 41(8): 1211-1218

    [19] Zhou J, Chen J W, Liang C H, et al. Quantum chemical investigation on the mechanism and kinetics of PBDE photooxidation by center dot OH: A case study for BDE-15 [J]. Environmental Science & Technology, 2011, 45(11): 4839-4845

    [20] Fang H S, Gao Y P, Li G Y, et al. Advanced oxidation kinetics and mechanism of preservative propylparaben degradation in aqueous suspension of TiO2and risk assessment of its degradation products [J]. Environmental Science & Technology, 2013, 47(6): 2704-2712

    [21] An T C, Gao Y P, Li G Y, et al. Kinetics and mechanism of (OH)-O-center dot mediated degradation of dimethyl phthalate in aqueous solution: Experimental and theoretical studies [J]. Environmental Science & Technology, 2014, 48(1): 641-648

    [22] 陳景文, 李雪花, 于海瀛, 等. 面向毒害有機物生態(tài)風險評價的(Q)SAR技術: 進展與展望[J]. 中國科學 (B輯:化學), 2008, 38(6): 461-474

    [23] Frisch M J, Trucks G W, Schlegel H B, et al. Gaussian 09, Revision D.01 [CP]. Wallingford: Gaussian, Inc., 2009

    [24] Nicolaescu A R, Wiest O, Kamat P V. Mechanistic pathways of the hydroxyl radical reactions of quinoline. 2. Computational analysis of hydroxyl radical attack at C atoms [J]. Journal of Physical Chemistry A, 2005, 109(12): 2829-2835

    [25] Zhao C Y, Boriani E, Chana A, et al. A new hybrid system of QSAR models for predicting bioconcentration factors (BCF) [J]. Chemosphere, 2008, 73(11): 1701-1707

    [26] ECOSAR [CP]. http://www.epa.gov/oppt/newchems/tools/21ecosar.htm. 2014

    [27] Fang X W, Schuchmann H P, von Sonntag C. The reaction of the OH radical with pentafluoro-, pentachloro-, pentabromo- and 2,4,6-triiodophenol in water: Electron transfer vs. addition to the ring [J]. Journal of the Chemical Society-Perkin Transactions, 2000(7): 1391-1398

    [28] Cramer C J, Truhlar D G. Implicit solvation models: Equilibria, structure, spectra, and dynamics [J]. Chemical Reviews, 1999, 99(8): 2161-2200

    [29] Benitez F J, Beltran-Heredia J, Acero J L, et al. Oxidation of several chlorophenolic derivatives by UV irradiation and hydroxyl radicals [J]. Journal of Chemical Technology and Biotechnology, 2001, 76(3): 312-320

    [30] Zimbron J A, Reardon K F. Hydroxyl free radical reactivity toward aqueous chlorinated phenols [J]. Water Research, 2005, 39 (5): 865-869

    [31] Peller J, Kamat P V. Radiolytic transformations of chlorinated phenols and chlorinated phenoxyacetic acids [J]. Journal of Physical Chemistry A, 2005, 109(42): 9528-9535

    [32] Xu F, Yu W N, Zhou Q, et al. Mechanism and direct kinetic study of the polychlorinated dibenzo-p-dioxin and dibenzofuran formations from the radical/radical cross-condensation of 2,4-dichlorophenoxy with 2-chlorophenoxy and 2,4,6-trichlorophenoxy [J]. Environmental Science & Technology, 2011, 45(2): 643-650

    [33] Xu F, Wang H, Zhang Q Z, et al. Kinetic properties for the complete series reactions of chlorophenols with OH radicals-relevance for dioxin formation [J]. Environmental Science & Technology, 2010, 44(4): 1399-1404

    [34] Buth J M, Arnold W A, McNeill K. Unexpected products and reaction mechanisms of the aqueous chlorination of cimetidine [J]. Environmental Science & Technology, 2007, 41(17): 6228-6233

    [35] Veith G, Broderius S. Structure-Toxicity Relationships for Industrial Chemicals Causing Type (II) Narcosis Syndrome [M]// Kaiser K E. ed. QSAR in Environmental Toxicology - II. Springer, 1987

    [36] Holcombe G W, Phipps G L, Fiandt J T. Effects of phenol, 2,4-dimethylphenol, 2,4-dichlorophenol, and pentachlorophenol on embryo, larval, and early-juvenile fathead minnows (Pimephales promelas) [J]. Archives of Environmental Contamination and Toxicology, 1982, 11(1): 73-78

    Theoretical Investigation on the Degradation Mechanism and Toxicity Evolution during Photochemical Transformation of 2,4-Dichlorophenol Mediated by·OH

    Gao Yanpeng1, JiYuemeng2, Li Guiying2, An Taicheng1,2,*

    1. State Key Laboratory of Organic Geochemistry and Guangdong Key Laboratory of Environmental Resources Utilization and Protection, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China 2. Institute of Environmental Health and Pollution Control and School of Environmental Science and Engineering, Guangdong University of Technology, Guangzhou 510006, China

    Received 20 December 2015 accepted 01 February 2016

    Recently, the photochemical degradation mechanisms of 2,4-dichlorophenol (2,4-DCP) and the consequent toxicities of its intermedates cause a worldwide concern due to its potential carcinogenesis, teratogenesis, and mutagenesis as well as ecotoxicological effects. Herein, the ·OH-initiated indirect photochemical transformation mechanism, environmental fate and aquatic toxicity of 2,4-DCP as well as its degradation intermediates in water environment were theoretically calculated in detail. The results show that 2,4-DCP can be degraded readily via ·OH-addition and H-abstraction pathways. The ·OH-addition was the predominant pathway at low temperature (<313 K), mainly producing 4,6-dichlorobenzene-1,3-diol; whereas the phenolic hydroxyl H atom was predominantly abstracted by ·OH at high temperature, mainly forming 2-chlorobenzene-p-benzenequinone. What’s more, dioxin would be also formed via the H-abstraction pathway in the high temperature water environment with limited concentration of reactive species. Furthermore, the computational toxicology results indicate that the transformation products from H-abstraction pathways possess an increased aquatic toxicity to fish compared with ·OH-addition products, and even exceeds one toxic level than 2,4-DCP. Therefore, more attentions should be paid to photochemical products as well as the original 2,4-DCP, particularly those products from H-abstraction pathways.

    2,4-dichlorophenol; photochemical transformation mechanisms; transformation products; theoretical chemical calculation; eco-toxicity; toxicity evolution

    10.7524/AJE.1673-5897.20151220001

    國家杰出青年科學基金項目(41425015);第57批博士后面上資助項目(2015M572375)

    高艷蓬(1987-),女,博士,在站博士后,研究方向為新興有機污染物的遷移轉化與風險評估的理論研究,E-mail: gaoyanpeng@gig.ac.cn

    *通訊作者(Corresponding author), E-mail: antc99@gdut.edu.cn

    2015-12-20 錄用日期:2016-02-01

    1673-5897(2016)2-179-09

    X171.5

    A

    簡介:安太成(1972—),男,環(huán)境科學博士,研究員,博導。主要從事新興毒害有機污染物(EOCs)的環(huán)境地球化學過程及其環(huán)境污染控制原理方面研究,發(fā)表論文250余篇。

    高艷蓬, 姬越蒙, 李桂英, 等. ·OH自由基介導2,4-二氯酚光降解機理及產(chǎn)物毒性變化特征的理論研究[J]. 生態(tài)毒理學報,2016, 11(2): 179-187

    Gao Y P, Ji Y M, Li G Y, et al. Theoretical investigation on the degradation mechanism and toxicity evolution during photochemical transformation of 2, 4-dichlorophenol mediated by ·OH [J]. Asian Journal of Ecotoxicology, 2016, 11(2): 179-187 (in Chinese)

    猜你喜歡
    光降解中間體毒性
    建立A注射液中間體中肉桂酸含量測定方法
    激發(fā)態(tài)和瞬態(tài)中間體的光譜探測與調(diào)控
    動物之最——毒性誰最強
    水體中布洛芬的間接光降解作用機理研究
    RGD肽段連接的近紅外量子點對小鼠的毒性作用
    PM2.5中煤煙聚集物最具毒性
    吸入麻醉藥的作用和毒性分析
    基于β-二亞胺配體的鋁氧硼六元環(huán)化合物和其中間體的合成、表征及其反應性研究
    水中磺胺類抗生素的光降解及富里酸對其光降解的影響
    LC-MS/MS測定馬來酸桂哌齊特原料合成中間體殘留
    欧美日本视频| 91精品伊人久久大香线蕉| 日产精品乱码卡一卡2卡三| 高清视频免费观看一区二区 | 久热久热在线精品观看| 性色avwww在线观看| 亚洲精品乱码久久久久久按摩| 国产伦一二天堂av在线观看| 蜜桃久久精品国产亚洲av| 国产精品熟女久久久久浪| 我要搜黄色片| 日日啪夜夜撸| 91久久精品电影网| 欧美人与善性xxx| 久热久热在线精品观看| 老司机福利观看| 极品教师在线视频| 麻豆精品久久久久久蜜桃| 国产成人a∨麻豆精品| 精品久久国产蜜桃| 亚洲欧洲国产日韩| 欧美三级亚洲精品| 国产乱人偷精品视频| 男人的好看免费观看在线视频| 国产精品三级大全| 少妇人妻一区二区三区视频| 看黄色毛片网站| 69av精品久久久久久| 插逼视频在线观看| 亚洲精品乱码久久久v下载方式| 国产精品国产三级国产专区5o | 久久久欧美国产精品| 国产在线一区二区三区精 | 国产欧美日韩精品一区二区| 男人狂女人下面高潮的视频| 亚洲最大成人手机在线| 一夜夜www| 九九久久精品国产亚洲av麻豆| 亚洲经典国产精华液单| 亚洲国产精品合色在线| 人人妻人人看人人澡| a级一级毛片免费在线观看| 老师上课跳d突然被开到最大视频| 能在线免费观看的黄片| 亚洲综合精品二区| 亚洲欧美清纯卡通| 日日干狠狠操夜夜爽| 久久久久久大精品| 国产精品,欧美在线| 精华霜和精华液先用哪个| 1000部很黄的大片| 黄色日韩在线| 91精品一卡2卡3卡4卡| 丝袜喷水一区| 久久久国产成人精品二区| 亚洲精品国产av成人精品| 国产高清国产精品国产三级 | 神马国产精品三级电影在线观看| 欧美日本视频| 国产精品永久免费网站| 美女高潮的动态| 免费看av在线观看网站| 国内精品宾馆在线| 视频中文字幕在线观看| 丰满人妻一区二区三区视频av| 只有这里有精品99| av福利片在线观看| www.av在线官网国产| 国产精品电影一区二区三区| 极品教师在线视频| 人妻少妇偷人精品九色| 久久午夜福利片| 久久综合国产亚洲精品| 我的女老师完整版在线观看| 亚洲真实伦在线观看| 国产午夜福利久久久久久| 九草在线视频观看| 亚洲性久久影院| 不卡视频在线观看欧美| 天堂av国产一区二区熟女人妻| 亚洲精品一区蜜桃| 成人亚洲精品av一区二区| 国产精品久久久久久精品电影小说 | 欧美成人午夜免费资源| 啦啦啦啦在线视频资源| 国产真实乱freesex| 性色avwww在线观看| 丝袜美腿在线中文| 两个人视频免费观看高清| 淫秽高清视频在线观看| 长腿黑丝高跟| 亚洲综合色惰| 亚洲美女视频黄频| 日本色播在线视频| 亚洲精品aⅴ在线观看| 夜夜看夜夜爽夜夜摸| 亚洲经典国产精华液单| 国产伦在线观看视频一区| 黄色欧美视频在线观看| 欧美日本视频| 在线a可以看的网站| 97人妻精品一区二区三区麻豆| 嫩草影院精品99| 韩国高清视频一区二区三区| 亚洲在线自拍视频| 日韩欧美在线乱码| 看免费成人av毛片| 国语自产精品视频在线第100页| 99视频精品全部免费 在线| 你懂的网址亚洲精品在线观看 | 97超视频在线观看视频| 美女国产视频在线观看| 久久精品久久久久久噜噜老黄 | 亚洲美女搞黄在线观看| 免费在线观看成人毛片| 国产视频内射| 美女高潮的动态| 干丝袜人妻中文字幕| 国产一级毛片在线| 少妇人妻精品综合一区二区| 国产精品久久视频播放| 午夜亚洲福利在线播放| 久久精品综合一区二区三区| 亚洲av成人av| 91久久精品国产一区二区三区| 国产黄a三级三级三级人| 日韩强制内射视频| 免费观看人在逋| 我要看日韩黄色一级片| 淫秽高清视频在线观看| 久久久精品94久久精品| 男女国产视频网站| 在线观看美女被高潮喷水网站| 成人特级av手机在线观看| 我要搜黄色片| 久久午夜福利片| 亚洲高清免费不卡视频| 久久国内精品自在自线图片| 禁无遮挡网站| 久久久国产成人免费| 一级av片app| 亚洲美女搞黄在线观看| 波多野结衣巨乳人妻| 国产伦在线观看视频一区| 白带黄色成豆腐渣| 美女脱内裤让男人舔精品视频| 国产av码专区亚洲av| 三级男女做爰猛烈吃奶摸视频| 小说图片视频综合网站| 两个人的视频大全免费| 神马国产精品三级电影在线观看| 欧美又色又爽又黄视频| 男女下面进入的视频免费午夜| 少妇被粗大猛烈的视频| 国产成人免费观看mmmm| 亚洲18禁久久av| 青青草视频在线视频观看| 久久午夜福利片| 久久久国产成人免费| 亚洲欧美一区二区三区国产| 国产单亲对白刺激| 在现免费观看毛片| 在线免费观看不下载黄p国产| 麻豆av噜噜一区二区三区| 午夜激情欧美在线| 午夜激情欧美在线| 亚洲aⅴ乱码一区二区在线播放| 一本久久精品| 最后的刺客免费高清国语| 男女边吃奶边做爰视频| 99九九线精品视频在线观看视频| 国产精品爽爽va在线观看网站| 久久国产乱子免费精品| 久久欧美精品欧美久久欧美| 久久久久国产网址| 日本黄色视频三级网站网址| 在线观看66精品国产| 中文亚洲av片在线观看爽| 日本-黄色视频高清免费观看| 午夜福利高清视频| 国产精品久久久久久久电影| 国产精品综合久久久久久久免费| h日本视频在线播放| 国产毛片a区久久久久| 国产白丝娇喘喷水9色精品| 成年版毛片免费区| h日本视频在线播放| 黄色日韩在线| 免费av不卡在线播放| 夜夜看夜夜爽夜夜摸| 五月玫瑰六月丁香| 听说在线观看完整版免费高清| 亚洲aⅴ乱码一区二区在线播放| 久久午夜福利片| 男人的好看免费观看在线视频| 2021少妇久久久久久久久久久| 99热精品在线国产| 午夜福利在线在线| 欧美一区二区亚洲| 国产高清国产精品国产三级 | 亚洲久久久久久中文字幕| 亚洲图色成人| 黄色一级大片看看| 亚洲欧美日韩无卡精品| 淫秽高清视频在线观看| 午夜爱爱视频在线播放| 亚洲图色成人| 在线观看66精品国产| 99热网站在线观看| 好男人在线观看高清免费视频| 中文字幕av在线有码专区| 免费看a级黄色片| 亚洲av熟女| 亚洲精品乱码久久久v下载方式| 久久精品91蜜桃| 综合色av麻豆| 成人二区视频| 男女那种视频在线观看| 国产国拍精品亚洲av在线观看| 亚洲成人久久爱视频| 老师上课跳d突然被开到最大视频| 麻豆久久精品国产亚洲av| 草草在线视频免费看| 免费搜索国产男女视频| 美女高潮的动态| 午夜免费男女啪啪视频观看| 欧美最新免费一区二区三区| 村上凉子中文字幕在线| 亚洲欧洲日产国产| 麻豆av噜噜一区二区三区| 天天一区二区日本电影三级| 日韩欧美精品v在线| 成人特级av手机在线观看| 午夜激情福利司机影院| 97超碰精品成人国产| 日韩人妻高清精品专区| 亚洲av.av天堂| 美女xxoo啪啪120秒动态图| 国产老妇伦熟女老妇高清| 国产精品永久免费网站| 国产色婷婷99| 天堂av国产一区二区熟女人妻| 亚洲成人av在线免费| 汤姆久久久久久久影院中文字幕 | 国产精品永久免费网站| 热99re8久久精品国产| 中文字幕免费在线视频6| 日本一本二区三区精品| 一级毛片aaaaaa免费看小| 伦精品一区二区三区| 精品国产露脸久久av麻豆 | 国内精品一区二区在线观看| 亚洲成人中文字幕在线播放| 日韩av在线大香蕉| 最近最新中文字幕大全电影3| 天堂中文最新版在线下载 | 亚洲精品456在线播放app| 欧美xxxx性猛交bbbb| a级一级毛片免费在线观看| 亚洲中文字幕一区二区三区有码在线看| 免费不卡的大黄色大毛片视频在线观看 | 国产三级在线视频| 亚洲国产精品专区欧美| 嫩草影院新地址| 简卡轻食公司| 国产成人a区在线观看| 国产在视频线在精品| 国产精品99久久久久久久久| 三级国产精品欧美在线观看| 亚洲av.av天堂| 国产伦一二天堂av在线观看| 亚洲精品,欧美精品| 亚洲,欧美,日韩| 1000部很黄的大片| 亚洲乱码一区二区免费版| 91在线精品国自产拍蜜月| 日本三级黄在线观看| av.在线天堂| 在线免费观看的www视频| 国产片特级美女逼逼视频| 看黄色毛片网站| 国产午夜精品久久久久久一区二区三区| 偷拍熟女少妇极品色| 晚上一个人看的免费电影| 啦啦啦观看免费观看视频高清| 久久久久久九九精品二区国产| 国产亚洲最大av| 亚洲人与动物交配视频| 久久久久久久国产电影| 在线观看66精品国产| 久久精品综合一区二区三区| 国产亚洲av嫩草精品影院| av专区在线播放| av在线观看视频网站免费| 久久韩国三级中文字幕| 国产日韩欧美在线精品| 国产av在哪里看| 舔av片在线| 简卡轻食公司| 日日摸夜夜添夜夜添av毛片| 日产精品乱码卡一卡2卡三| 我的女老师完整版在线观看| 国产免费视频播放在线视频 | 国产精品一区二区性色av| av免费观看日本| 免费观看性生交大片5| 麻豆乱淫一区二区| a级一级毛片免费在线观看| 身体一侧抽搐| 国产不卡一卡二| 色播亚洲综合网| 国产色爽女视频免费观看| 免费看a级黄色片| 国产精品精品国产色婷婷| 三级国产精品片| 亚洲精品乱码久久久v下载方式| 国产精品电影一区二区三区| 国产精品一区二区三区四区免费观看| www.av在线官网国产| 边亲边吃奶的免费视频| 欧美日韩精品成人综合77777| 亚洲av成人av| 亚洲aⅴ乱码一区二区在线播放| 好男人在线观看高清免费视频| 一级av片app| 天堂av国产一区二区熟女人妻| 国产黄片美女视频| 91av网一区二区| 亚洲综合精品二区| 看片在线看免费视频| 日韩精品青青久久久久久| 亚洲自偷自拍三级| 久久久久网色| 最近最新中文字幕大全电影3| 国产免费一级a男人的天堂| 国产黄a三级三级三级人| 好男人在线观看高清免费视频| 国产精品麻豆人妻色哟哟久久 | 色5月婷婷丁香| 91aial.com中文字幕在线观看| 舔av片在线| 久久久国产成人精品二区| av在线亚洲专区| 91aial.com中文字幕在线观看| 欧美性猛交╳xxx乱大交人| 秋霞在线观看毛片| 亚洲成av人片在线播放无| 欧美一区二区精品小视频在线| 国产一区二区在线av高清观看| 成人性生交大片免费视频hd| 欧美97在线视频| 高清午夜精品一区二区三区| 亚洲怡红院男人天堂| 日本免费a在线| 国产亚洲一区二区精品| 可以在线观看毛片的网站| 大香蕉97超碰在线| 国内精品宾馆在线| 久久精品国产鲁丝片午夜精品| av福利片在线观看| 免费播放大片免费观看视频在线观看 | 中文字幕av成人在线电影| 精品久久国产蜜桃| 亚洲人与动物交配视频| 国产亚洲精品av在线| 亚洲精品亚洲一区二区| 日韩强制内射视频| 久久99热6这里只有精品| 黄色一级大片看看| 欧美精品国产亚洲| 一区二区三区四区激情视频| 韩国高清视频一区二区三区| 亚洲人成网站高清观看| 欧美激情久久久久久爽电影| 22中文网久久字幕| 高清日韩中文字幕在线| 成人午夜高清在线视频| 免费在线观看成人毛片| 日韩视频在线欧美| 国产午夜福利久久久久久| 搡女人真爽免费视频火全软件| 久久久精品欧美日韩精品| 欧美xxxx性猛交bbbb| 午夜视频国产福利| 久久国产乱子免费精品| 一边摸一边抽搐一进一小说| 尤物成人国产欧美一区二区三区| 小说图片视频综合网站| 长腿黑丝高跟| 少妇熟女欧美另类| 国产白丝娇喘喷水9色精品| 色综合亚洲欧美另类图片| 男人和女人高潮做爰伦理| 波多野结衣巨乳人妻| av.在线天堂| 国产一区二区在线av高清观看| 国产乱来视频区| 国产伦精品一区二区三区四那| 亚洲电影在线观看av| 在线观看一区二区三区| 一级毛片电影观看 | 久久精品夜色国产| 精品久久久久久久人妻蜜臀av| 国产精品综合久久久久久久免费| 男女边吃奶边做爰视频| 一本一本综合久久| 中文字幕久久专区| 精品人妻偷拍中文字幕| 床上黄色一级片| 水蜜桃什么品种好| 内地一区二区视频在线| 免费看光身美女| 亚洲av成人精品一二三区| 日本av手机在线免费观看| 99热这里只有精品一区| 看十八女毛片水多多多| 国产国拍精品亚洲av在线观看| 最近手机中文字幕大全| 观看免费一级毛片| 国产老妇女一区| 日韩欧美三级三区| 狂野欧美白嫩少妇大欣赏| 真实男女啪啪啪动态图| 大香蕉久久网| 久久久久精品久久久久真实原创| 国产成人a∨麻豆精品| 蜜桃久久精品国产亚洲av| 岛国毛片在线播放| 日韩制服骚丝袜av| 久久久国产成人免费| 伦精品一区二区三区| 日本wwww免费看| 最近2019中文字幕mv第一页| 免费看光身美女| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 能在线免费观看的黄片| 亚洲欧美成人精品一区二区| 草草在线视频免费看| 亚洲人与动物交配视频| 黄色配什么色好看| 国产成人福利小说| 亚洲国产精品国产精品| 汤姆久久久久久久影院中文字幕 | 国内精品美女久久久久久| 国产一区二区亚洲精品在线观看| 一级av片app| 我的老师免费观看完整版| 视频中文字幕在线观看| 天堂网av新在线| 久久久久精品久久久久真实原创| 在线免费十八禁| 久久久久久久久久久免费av| 综合色av麻豆| 一级爰片在线观看| 日韩精品有码人妻一区| 麻豆成人午夜福利视频| 天天一区二区日本电影三级| 国产成人91sexporn| 国产精品一区www在线观看| 久久久久久久久久成人| 别揉我奶头 嗯啊视频| 中文乱码字字幕精品一区二区三区 | av在线蜜桃| 中国国产av一级| 欧美成人精品欧美一级黄| 高清午夜精品一区二区三区| 九色成人免费人妻av| 国产精品永久免费网站| 黄片无遮挡物在线观看| 一边摸一边抽搐一进一小说| 伦精品一区二区三区| 精品久久久久久久久av| 精品99又大又爽又粗少妇毛片| 村上凉子中文字幕在线| 一级毛片我不卡| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品人妻久久久影院| 一级av片app| 中文字幕av成人在线电影| 午夜免费激情av| 不卡视频在线观看欧美| 丰满乱子伦码专区| 亚洲av成人精品一区久久| 日韩欧美精品v在线| 久久久国产成人免费| 精品国产一区二区三区久久久樱花 | 中文在线观看免费www的网站| 波多野结衣高清无吗| 国语对白做爰xxxⅹ性视频网站| 欧美变态另类bdsm刘玥| 免费电影在线观看免费观看| 国产真实乱freesex| 免费看光身美女| 97在线视频观看| 老女人水多毛片| 日本午夜av视频| 最近的中文字幕免费完整| 国产老妇女一区| 两个人视频免费观看高清| 看黄色毛片网站| 亚洲精品色激情综合| 内地一区二区视频在线| 午夜福利在线在线| 久久久色成人| 国产毛片a区久久久久| 噜噜噜噜噜久久久久久91| 久久久久久久久久久丰满| 成人国产麻豆网| 99久国产av精品| 欧美又色又爽又黄视频| 亚洲欧美清纯卡通| 国产黄片美女视频| 一个人免费在线观看电影| 秋霞伦理黄片| 亚洲熟妇中文字幕五十中出| 久久久久久久亚洲中文字幕| 久久久久精品久久久久真实原创| 免费黄色在线免费观看| 在线天堂最新版资源| 99国产精品一区二区蜜桃av| 免费av观看视频| 99在线人妻在线中文字幕| 国产在线一区二区三区精 | 国产精品一及| 高清毛片免费看| 日本色播在线视频| 水蜜桃什么品种好| 国产极品精品免费视频能看的| 99国产精品一区二区蜜桃av| 亚洲国产日韩欧美精品在线观看| 美女国产视频在线观看| 建设人人有责人人尽责人人享有的 | av国产久精品久网站免费入址| 国产伦理片在线播放av一区| 国产亚洲91精品色在线| 国产精品永久免费网站| 少妇丰满av| 免费看日本二区| 只有这里有精品99| 免费在线观看成人毛片| 国产精品综合久久久久久久免费| 亚洲精品国产av成人精品| 国产黄片美女视频| 岛国在线免费视频观看| 国产又色又爽无遮挡免| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av男天堂| 欧美成人午夜免费资源| 亚洲成人精品中文字幕电影| 水蜜桃什么品种好| 亚洲av日韩在线播放| 丰满少妇做爰视频| 国产高清不卡午夜福利| 色播亚洲综合网| 看免费成人av毛片| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华液的使用体验| 精品久久久久久久久av| 国产成人精品一,二区| av免费在线看不卡| 亚洲欧美精品综合久久99| 国产av在哪里看| 久久人人爽人人片av| 亚洲精品乱码久久久久久按摩| 日本欧美国产在线视频| 熟女电影av网| 97热精品久久久久久| 伦理电影大哥的女人| 校园人妻丝袜中文字幕| 夜夜爽夜夜爽视频| 国产精品国产三级国产av玫瑰| 久久久久久久久久久丰满| 久久99热这里只有精品18| 蜜桃久久精品国产亚洲av| 国产不卡一卡二| 国产91av在线免费观看| www日本黄色视频网| 如何舔出高潮| 国产一区二区三区av在线| 极品教师在线视频| 精品人妻视频免费看| 久久久国产成人精品二区| 日本色播在线视频| 精品99又大又爽又粗少妇毛片| 国产不卡一卡二| 中文字幕亚洲精品专区| 久久国产乱子免费精品| 一卡2卡三卡四卡精品乱码亚洲| 久久这里只有精品中国| 欧美一级a爱片免费观看看| 长腿黑丝高跟| 亚洲国产日韩欧美精品在线观看| 一夜夜www| 夫妻性生交免费视频一级片| 18禁裸乳无遮挡免费网站照片| 久久精品91蜜桃| 简卡轻食公司| 自拍偷自拍亚洲精品老妇| 久久草成人影院| 亚洲怡红院男人天堂| 午夜精品一区二区三区免费看| 亚洲最大成人av| 看免费成人av毛片| 男插女下体视频免费在线播放| 91久久精品电影网| 99热全是精品| eeuss影院久久| 乱系列少妇在线播放| 国产亚洲5aaaaa淫片| 国产精品永久免费网站| 成年女人永久免费观看视频| 日韩国内少妇激情av| 毛片女人毛片| 91aial.com中文字幕在线观看| 亚洲国产精品国产精品|