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

    山東近海陸源耗氧有機物生物可利用性及其降解動力學(xué)研究?

    2017-06-05 15:08:52馬云鵬戴愛泉李克強王修林
    關(guān)鍵詞:耗氧陸源常數(shù)

    馬云鵬, 戴愛泉,2, 李克強??, 王修林

    (1.中國海洋大學(xué)海洋化學(xué)理論與工程技術(shù)教育部重點實驗室,化學(xué)化工學(xué)院,山東 青島 266100; 2.青島市環(huán)境保護科學(xué)研究院,山東 青島 266000)

    山東近海陸源耗氧有機物生物可利用性及其降解動力學(xué)研究?

    馬云鵬1, 戴愛泉1,2, 李克強1??, 王修林1

    (1.中國海洋大學(xué)海洋化學(xué)理論與工程技術(shù)教育部重點實驗室,化學(xué)化工學(xué)院,山東 青島 266100; 2.青島市環(huán)境保護科學(xué)研究院,山東 青島 266000)

    針對因不同陸源耗氧有機物組成和生物可利用性的差異,單一表觀降解模式難以滿足近海水環(huán)境管理中對耗氧有機物濃度準(zhǔn)確模擬計算的問題,選擇山東省近海典型陸源污染來源,通過實驗室受控培養(yǎng)實驗,研究不同來源耗氧有機物(以COD表征)的生物可利用性及其降解動力學(xué)規(guī)律。結(jié)果顯示,耗氧有機物的降解基本符合準(zhǔn)一級動力學(xué)規(guī)律,即先快后慢最后趨于穩(wěn)定的降解規(guī)律,且不同陸源耗氧有機物降解速率常數(shù)及生物可利用性存在顯著差異(P<0.01)。其中,紡織印染廠來源耗氧有機物的生物可利用性最大,約為92%,其次為污水處理廠和小清河來源,而虞河和啤酒廠來源最??;啤酒廠來源耗氧有機物降解速率常數(shù)最高,約為0.40 d-1,其次為污水處理廠來源,而小清河、虞河以及印染廠來源最低。根據(jù)降解速率常數(shù)的不同,可以將不同來源的耗氧有機物分為四類:食品工業(yè)源(啤酒廠)生物可利用組分,城鎮(zhèn)生活源(污水處理廠)生物可利用組分,河流綜合源(小清河、虞河和印染廠)生物可利用組分以及難降解組分。污水處理廠來源耗氧有機物培養(yǎng)實驗結(jié)果表明,耗氧有機物的降解速率常數(shù)隨溫度的升高及攪拌速率的增強而增大,生物可利用性則隨水動力條件的增強而增大。本文研究結(jié)果可為提高水環(huán)境管理中水質(zhì)模擬計算準(zhǔn)確性提供實驗支持。

    降解動力學(xué);生物可利用性;耗氧有機物(COD);降解速率常數(shù);山東近海

    近年來,山東省近海海域,如萊州灣、乳山灣、膠州灣等半封閉海灣耗氧有機物污染嚴(yán)重,導(dǎo)致季節(jié)性缺氧或低氧區(qū)域的形成[1-2],對海區(qū)初級生產(chǎn)力和浮游植物的分布產(chǎn)生直接影響[3]。此外,山東近海海域赤潮、綠潮頻發(fā)[4-7],可能與耗氧有機物能夠為赤潮生物的生長提供碳源有關(guān)[8]。然而,隨著山東省社會經(jīng)濟的高速發(fā)展,山東近海海域耗氧有機物污染尚未得到有效緩解,特別是河口及膠州灣等半封閉海灣灣頂區(qū)域以化學(xué)需氧量(COD)表征的耗氧有機物超標(biāo)海域面積占比呈明顯升高趨勢[9]。實施陸源COD等污染物排海總量控制管理是近海水環(huán)境質(zhì)量改善的有效措施,也是國家“水污染防治行動計劃”的重要內(nèi)容,在技術(shù)支撐上需要COD污染物水質(zhì)模擬計算準(zhǔn)確性的提高,進而可提高COD污染物總量控制工作的效率性。

    陸源污染排放是自然水體中耗氧有機物的主要來源之一[8,10]。對于陸源耗氧有機物生物可利用性的研究始于1927年關(guān)于生活污水中有機污染物降解規(guī)律的研究[11]。耗氧有機物的生物可利用性主要通過培養(yǎng)實驗手段研究。Skopintsev和Ovtshinnikova在有氧條件下將生活污水與天然水混合,模擬研究了自然水體中陸源耗氧有機物的微生物降解情況[12]。Krause等的研究結(jié)果表明,耗氧有機物組分復(fù)雜,其中蛋白質(zhì)和碳水化合物易被微生物利用,而腐殖質(zhì)組分較難分解[13-14]。Richards等在對Nitinat湖和黑海的有機物降解動力學(xué)的研究中提出耗氧有機物組分中不穩(wěn)定的有機化合物最易被生物分解利用,因此降解實驗開始時降解速率較大,而后期較穩(wěn)定的有機成分的降解帶來的影響增加而導(dǎo)致降解速率減小[15-16]。隨著人類活動影響的增大,對于耗氧有機物生物可利用性的研究多結(jié)合在對于工業(yè)污水[17-19]、飲用水[20-21]及生活污水[22-24]處理技術(shù)的研究中,并開始應(yīng)用于廢水生物處理的數(shù)學(xué)模式[25]及各類海洋水質(zhì)模型[26]中。

    對于水體中耗氧有機物的降解動力學(xué)過程,多數(shù)學(xué)者認(rèn)為其符合一級反應(yīng)動力學(xué)方程[27-28]。一級動力學(xué)方程中的表觀降解速率常數(shù),可以反映耗氧有機污染物在水體中自凈降解速度的快慢[29]。降解系數(shù)亦可通過培養(yǎng)實驗[30-31]確定。但亦有學(xué)者通過人工添加與自然水體中成分近似的有機試劑對耗氧有機物的降解動力學(xué)過程進行研究,并提出了可以運用米氏(Michaelis-Menten)方程[32]及Monod方程[33-35]描述底物耗氧有機物變化過程。耗氧有機物作為混合物,其生物可利用性及降解動力學(xué)主要與組成有關(guān),不同成分間可能會相互抑制或促進,導(dǎo)致耗氧有機物的表觀降解速率發(fā)生變化。例如,Reardon等將由不同有機物混合而成的耗氧有機物生物降解動力學(xué)進行比較,結(jié)果顯示其生物可利用性和降解速率有顯著不同[36]。同時,耗氧有機物的生物可利用性及降解動力學(xué)還受溫度、鹽度、污染物(底物)起始濃度、微生物種類和數(shù)量、pH及水動力條件等環(huán)境因素的影響[22-23,37-38]。宋學(xué)宏等研究了陽澄湖河蟹養(yǎng)殖區(qū)水體中耗氧有機物生化降解規(guī)律,結(jié)果顯示水溫與降解系數(shù)密切相關(guān)[39]。郭棟鵬等對黃海南部海水耗氧有機物降解動力學(xué)的研究,指出初始濃度和溫度是影響其降解速率的重要因素,而氯度、pH和懸浮泥沙含量影響不顯著[40]。

    隨著社會經(jīng)濟的高速發(fā)展,陸源污染物結(jié)構(gòu)發(fā)生了顯著變化,生活污水和工業(yè)廢水占比增大[41],導(dǎo)致陸源耗氧有機物組成和性質(zhì)發(fā)生變化。不同耗氧有機物的組成不同,其生物可利用性和降解動力學(xué)規(guī)律可能也存在差異[34]。目前,相關(guān)的研究難以滿足近海水環(huán)境管理對耗氧有機物濃度準(zhǔn)確模擬計算的要求。因此,本文通過選擇5個山東省近海典型陸源污染源,通過實驗室受控培養(yǎng)實驗,以COD為表征指標(biāo),運用準(zhǔn)一級動力學(xué)方程研究耗氧有機物生物可利用性及其降解動力學(xué)規(guī)律,為提高水環(huán)境管理中水質(zhì)模擬計算準(zhǔn)確性提供實驗支持。

    1 材料與方法

    1.1 取樣

    山東省陸源耗氧有機物來源主要可以分為農(nóng)業(yè)源、工業(yè)源和城鎮(zhèn)生活源,不同來源的陸源耗氧有機物組成差別較大。為使污染物來源能夠基本涵蓋主要的入海污染物來源,選取小清河下游、虞河下游、青島麥島公共污水處理廠排污口、麗源印染廠排污口及青島第二啤酒廠排污口作為水樣采集點(見圖1)。其中,小清河代表綜合源,虞河代表農(nóng)業(yè)源,印染廠和啤酒廠分別代表各自行業(yè)工業(yè)源,直排海公共污水處理廠代表城鎮(zhèn)生活源(見表1)。應(yīng)當(dāng)指出,山東省沿海地區(qū)工業(yè)企業(yè)種類繁多,耗氧有機物來源以紡織、印染、食品、機械制造、造紙等行業(yè)為主,本研究選取食品行業(yè)的啤酒廠和紡織印染行業(yè)的紡織印染廠作為典型工業(yè)源。海水取自萊州灣(初始COD濃度1.7 mg/L)和膠州灣(初始COD濃度1.1 mg/L)近岸海域。河流樣品均在斷面中間位置采樣且均取自表層。樣品冷藏保存,并迅速轉(zhuǎn)移至實驗室。

    (WWTF、YR、WQH、YH和PJ分別代表青島麥島污水處理廠、紡織印染廠、小清河、虞河和青島第二啤酒廠來源.WWTF, YR, WQH, YH and PJ separately stand for Qingdao Maidao wastewater treatment plant, printing and dyeing mill, Xiaoqing river, Yu river and Qingdao No.2 brewery)

    圖1 陸源`

    Fig.1 The sampling stations of typical terrigenous

    oxygen-consumed organic matter and the basin of rivers

    1.2 培養(yǎng)實驗

    耗氧有機物的濃度以COD表征,采用堿性高錳酸鉀法測定,避免在酸性介質(zhì)中氯離子的干擾。

    針對不同來源的耗氧有機物在海域中的生物可利用性和降解動力學(xué)實驗,參考排污口混合區(qū)COD的濃度設(shè)置不同濃度梯度(見表1),將不同來源的水樣與海水混合,其中紡織印染廠、虞河及小清河三組培養(yǎng)實驗采用萊州灣海水,啤酒廠及污水處理廠培養(yǎng)實驗采用膠州灣海水,以不添加陸源耗氧有機物的海水樣品為空白對照組。

    陸源耗氧有機物主要來自于生活源,青島市麥島污水處理廠作為典型城鎮(zhèn)生活源,環(huán)境因素對其耗氧有機物的影響具有代表性。為此,以青島市麥島污水處理廠為例,設(shè)置了不同溫度、水動力條件下的培養(yǎng)實驗。其中,根據(jù)山東近海水溫變化范圍,分別設(shè)置12、16、20和26℃ 4個溫度梯度;采用攪拌轉(zhuǎn)速表征水動力條件,在相同初始濃度及溫度等條件下,設(shè)置0、200、391、579和810 r/min 5個轉(zhuǎn)速梯度。

    表1 各采樣點編號、采樣位置和初始濃度設(shè)置

    分別于2013年8、10月、2014年3、5月開展不同季節(jié)(溫度)培養(yǎng)實驗,于2014年5月開展不同污染物來源、不同初始污染物濃度及不同水動力條件培養(yǎng)實驗。培養(yǎng)實驗裝置為15 L敞口對苯二甲酸乙二酯(PET)桶,培養(yǎng)體積12 L,每個實驗組設(shè)置兩個平行樣。每天定時搖動,起到大氣復(fù)氧作用,并保持微小顆粒呈懸浮狀態(tài)。使用日光燈模擬光照,光照強度為50 W·m-2,光暗比為12∶12。使用水浴控制實驗體系溫度。實驗期間,實驗組鹽度平均為(17.48±6.03)。在實驗開始階段,每天上午9∶00時從培養(yǎng)體系中取樣測定,實驗后期逐漸延長采樣周期。

    1.3 數(shù)據(jù)處理方法

    根據(jù)文獻研究結(jié)果,有機物降解動力學(xué)過程在較短時間內(nèi)(約6~8 d)符合一級反應(yīng)動力學(xué)方程[37,39]。但在較完整的COD降解周期內(nèi)存在難降解的有機成分[30],因此,應(yīng)采用能夠體現(xiàn)難降解COD組分的準(zhǔn)一級動力學(xué)方程描述COD降解動力學(xué)過程[42-43]:

    COD(t)=BCOD·exp(-k·t)+RCOD。

    (1)

    其中:COD、BCOD和RCOD分別表示耗氧有機物、易降解和難降解有機物濃度,mg/L;k表示降解速率常數(shù),d-1;t表示時間。使用origin8.0軟件對各組已扣除海水COD空白值的結(jié)果進行非線性擬合。

    根據(jù)方程(2)計算耗氧有機物生物可利用性(BCOD%):

    (2)

    2 實驗結(jié)果

    2.1 不同來源耗氧有機物的培養(yǎng)實驗結(jié)果

    在20 ℃實驗條件下,不同來源不同初始濃度組的COD濃度變化及擬合曲線如圖2所示。結(jié)果表明,所有實驗組COD濃度變化均呈現(xiàn)先快后慢的下降趨勢,其降解動力學(xué)過程均符合方程1所述動力學(xué)過程。

    分別將5個不同來源的耗氧有機物實驗組中對應(yīng)不同初始濃度COD(0)的參數(shù)k和BCOD%進行拉格朗日插值到對應(yīng)初始濃度為3.0、4.0、5.0、6.0、7.0 mg/L,然后進行雙因素顯著性分析。結(jié)果表明,對于相同來源、不同COD初始濃度的實驗組的降解速率常數(shù)k和BCOD%值,組間均無顯著性差異(P >0.05)。說明不同COD初始濃度下,k值和BCOD%值基本恒定,可取不同初始濃度的平均值和偏差表示各COD污染來源的k值和BCOD%值。

    在20 ℃、無攪拌條件下各不同耗氧有機物來源的實驗組的k值和BCOD%值如圖3所示。

    結(jié)果表明,不同來源的耗氧有機物的降解速率常數(shù)和生物可利用性不同。其中,啤酒廠來源耗氧有機物降解速率常數(shù)k最大,約為(0.40±0.026) d-1,其次為污水處理廠,約為(0.28±0.018) d-1,紡織印染廠、小清河和虞河k值最小,約為0.042~0.072 d-1。紡織印染廠來源的耗氧有機物的生物可利用性最大,BCOD%值為(92±2.9)%,其次分別為小清河和污水處理廠,BCOD%值分別為(70±6.4)%和(58±1.2)%,虞河和啤酒廠最小,BCOD%分別為(24±5.3)%和(25±8.3)%。

    2.2 不同溫度條件下的培養(yǎng)實驗結(jié)果

    不同溫度條件下麥島污水處理廠COD濃度隨時間變化情況及擬合曲線如圖4所示。不同溫度條件下的COD降解動力學(xué)過程均符合方程(1)所示先快后慢的降解過程。各組COD降解速率常數(shù)k和BCOD%值如圖5所示。結(jié)果表明,COD降解速率常數(shù)k隨溫度增加而增大,存在顯著性差異(P<0.05)。 其中,溫度為12℃時,降解速率常數(shù)k值最小,為(0.14±0.032)d-1,溫度為26℃時,k值最大,約為(0.48±0.056)d-1。k隨水溫的升高而升高,表明當(dāng)環(huán)境溫度升高時,COD污染物降解速率會加快,與張世坤等的研究結(jié)果一致[44]。而BCOD%在不同溫度下組間差異不顯著(P>0.1)。

    (其中,WWTF、YR、WQH、YH和PJ分別代表污水處理廠、紡織印染廠、小清河、虞河和啤酒廠來源;編號1~5表示COD初始濃度由低到高的不同組別。(WWTF, YR, WQH, YH and PJ separately stand for Qingdao Maidao wastewater treatment plant, printing and dyeing mill, Xiaoqing river, Yu river and Qingdao No.2 brewery; ID from No.1 to No.5 are different groups as initial concentration of COD.)

    圖2 耗氧有機物(COD)降解實驗結(jié)果(點)及擬合結(jié)果(線)

    Fig.2 The degradation experiment’s result (points) and fitting result (lines) of oxygen-consumed organic matter

    圖3 不同來源耗氧有機物降解速率常數(shù)(k)和生物可利用性(BCOD%)Fig.3 The degradation rate constants (k) and bioavailability (BCOD%) of oxygen-consumed organic matter from different sources

    2.3 不同水動力條件下的培養(yǎng)實驗結(jié)果

    針對污水處理廠排污口來源耗氧有機物,不同水動力條件下COD濃度隨時間變化情況及擬合曲線如圖6所示。結(jié)果表明,隨著攪拌速率增加,降解速率加快。其中,COD降解速率常數(shù)k隨攪拌速率增加呈上升趨勢(見圖7),與Devi R.等對COD降解動力學(xué)的研究結(jié)果一致[23]。攪拌速率對BCOD%值的影響與k相同,污水處理廠來源的耗氧有機物的BCOD%值呈上升趨勢,最大生物可利用性變化幅度約5.2%。

    圖4 不同溫度條件下耗氧有機物(COD)降解實驗結(jié)果Fig.4 The degradation experiment results of oxygen- consumed organic matter in different temperature

    圖5 溫度對降解速率常數(shù)(k)及生物 可利用性(BCOD%)的影響Fig.5 The effect of temperature on degradation rate constants (k) and bioavailability (BCOD%)

    圖6 不同攪拌速率條件下污水處理廠 來源耗氧有機物(COD)降解曲線擬合Fig.6 The degradation curves of oxygen-consumed organic matter from WWTF in different stirring conditions

    圖7 攪拌速率對耗氧有機物降解速率 常數(shù)(k)和生物可利用性(BCOD%)的影響Fig.7 The effect of stirring rate conditions on oxygen-consumed organic matter′s degradation rate constants (k) and bioavailability (BCOD%)

    3 討論

    3.1 不同來源對耗氧有機物生物可利用性及降解速率常數(shù)的影響

    對于降解速率常數(shù)k,以綜合源為主的小清河和以農(nóng)業(yè)面源為主的虞河的耗氧有機物k值較小;來自城鎮(zhèn)生活源為主的青島市麥島污水處理廠耗氧有機物的k值相對于農(nóng)業(yè)和綜合源明顯升高。這與來自入海河流和公共污水處理廠的有機污染物組成密切相關(guān)。河流中的有機污染物多來自于農(nóng)田、林地等沖刷帶來的腐殖質(zhì)及其它各類有機質(zhì)的分解產(chǎn)物,而腐殖質(zhì)等通常為難降解有機質(zhì),因此k值較小。而城鎮(zhèn)生活污水中的有機質(zhì)多為淀粉、糖類及蛋白質(zhì)等易降解物質(zhì),因此k值相比河流的k值明顯增大。對于不同直排海企業(yè)的耗氧有機物降解速率常數(shù)k差異較大,平均為(0.24±0.16) d-1,這與工業(yè)企業(yè)的生產(chǎn)工藝、產(chǎn)品性質(zhì)等因素有密切關(guān)系。其中,青島第二啤酒廠的耗氧有機物的降解速率常數(shù)是所有實驗組別中最高的;而來自麗源印染廠的耗氧有機物k值較小,與河流農(nóng)業(yè)源和綜合源的k值近似。顯著性分析表明,不同污染來源的k值組間差異性不同,紡織印染廠(YR)和小清河(XQH)的k值無顯著性差異(P>0.05),而其它各組間均存在顯著性差異(P<0.01,表2)。

    生物可利用性(BCOD%)與降解速率常數(shù)不同,農(nóng)業(yè)面源為主的虞河和綜合源為主的小清河的BCOD%值相差較大,這是由于小清河流域同時容納了生活源和工業(yè)源污染來源的耗氧有機物。顯著性分析表明,不同耗氧有機物來源組間差異性不同(見表2),啤酒廠(PJ)和虞河(YH)來源的耗氧有機物的BCOD%值無顯著性差異(P>0.05),而其它來源組間差異性均顯著(P<0.01)。

    表2 不同來源的耗氧有機物組間降解速率常數(shù)(k)和生物可利用性(BCOD%)顯著性分析

    Note: ++ 數(shù)據(jù)在P<0.01 條件下有顯著性差異 There are significant difference whenp<0.01+ 數(shù)據(jù)在P<0.05條件下有顯著性差異 There are significant difference whenp<0.05—數(shù)據(jù)在P<0.05條件下無顯著性差異There are no significant difference whenp<0.05

    使用SPSS20.0的聚類分析功能對5個不同來源的耗氧有機物降解速率常數(shù)k值進行聚類分析,采用k-means cluster聚類方法,結(jié)果如圖8所示。

    圖8 按降解速率常數(shù)(k)值聚類分析示意圖Fig.8 The cluster analysis diagram based on the degradation rate constants (k)

    結(jié)果表明,在溫度為20℃,初始COD濃度為7.0 mg/L時,根據(jù)降解速率常數(shù)k,可將耗氧有機物生物可利用組分分為3類(見表3):啤酒廠食品工業(yè)源耗氧有機物(BCOD1),污水處理廠城鎮(zhèn)生活源耗氧有機物(BCOD2),小清河、虞河、紡織印染廠綜合源耗氧有機物(BCOD3)。而不同來源耗氧有機物難降解組分為單獨一類(RCOD)。虞河和小清河屬于農(nóng)業(yè)面源污染,有機物經(jīng)長距離輸送,生物可降解性和生物可利用性降低,且組成上以難降解的腐殖質(zhì)居多,故k值較低;而同為工業(yè)源,由于生產(chǎn)對象及生產(chǎn)工藝的不同,啤酒廠污水中主要是氨基酸、糖類等易降解有機物,而印染廠則主要是芳香、雜環(huán)類等難降解有機物,結(jié)果啤酒廠COD的k值遠(yuǎn)高于印染廠;作為城鎮(zhèn)生活源的污水處理廠,其中的污染物多為蛋白質(zhì)、胺類等易降解的有機物,其COD的k值僅次于啤酒廠。

    3.2 溫度對降解速率常數(shù)和生物可利用性的影響

    來源于青島市麥島污水處理廠的耗氧有機物,在實驗設(shè)定的水溫范圍內(nèi),降解速率常數(shù)k與溫度T(℃)的關(guān)系可以用如下關(guān)系式表示:

    kT=k20(℃)θT-20。

    (4)

    根據(jù)實驗結(jié)果,得=1.079,與Schroopfer(1.056)和Moore(1.065)的研究結(jié)果基本一致。

    根據(jù)Arrheniux方程[45]:

    (5)

    圖9 lnk與關(guān)系示意圖Fig.9 The relationship between lnk and 1/T

    (6)

    計算得到耗氧有機物降解過程的表觀活化能E為58.10 kJ/mol。比陳松等[46]的研究結(jié)果27.00 kJ/mol明顯增大,推測是由于陳松等的實驗時間約為6d,而本研究中的實驗時間相對較長(20d),有更多的難降解耗氧有機物組分參與反應(yīng),導(dǎo)致其表觀活化能結(jié)果升高。

    表3 陸源耗氧有機物聚類分析結(jié)果及參數(shù)范圍

    3.3 攪拌速率對降解速率常數(shù)和生物可利用性的影響

    COD的降解實際上是一個微生物的耗氧降解過程。針對麥島污水處理廠耗氧有機物,在溫度20℃條件下,隨著攪拌速率的增加,一方面,實驗體系中的溶解氧含量增加,另一方面,在攪拌的條件下,水體中的有機物與微生物接觸得更加充分,導(dǎo)致COD降解速率增大,生物可利用性提高[47]。若將不同攪拌速率近似看作不同水動力條件,可以預(yù)測,耗氧有機物降解速率常數(shù)k和生物可利用性BCOD%值隨湍流速率增加均有緩慢上升的趨勢。

    4 結(jié)語

    不同來源的耗氧有機物降解速率常數(shù)k和生物可利用性BCOD%值存在顯著差異,按照降解速率常數(shù)聚類分析結(jié)果,可將不同來源的耗氧有機物分為4類(BCOD1、BCOD2、BCOD3和RCOD)。陸源耗氧有機物降解速率常數(shù)k和生物可利用性BCOD%受溫度、攪拌速率和鹽度等環(huán)境因素影響,其中,k隨溫度升高而顯著增加,k和BCOD%均隨攪拌速率增強而顯著增加。研究結(jié)果可應(yīng)用于水質(zhì)模型的修正中,為提高水環(huán)境管理中水質(zhì)模擬計算準(zhǔn)確性提供實驗支持。

    [1] 夏斌, 張曉理, 崔毅, 等. 夏季萊州灣及附近水域理化環(huán)境及營養(yǎng)現(xiàn)狀評價[J]. 漁業(yè)科學(xué)進展, 2009, 30(3): 103-111. XIA Bin, ZHANG Xiaoli, CUI Yi. The evaluation of physicochemical environment and nutritional status for Laizhou Bay in summer [J]. Marine Fishery Development, 2009, 30(3): 103-111.

    [2] 劉軍, 臧家業(yè), 冉祥濱, 等. 乳山灣外低氧海域沉積物中有機碳、氮、磷及其形態(tài)特征分析[J]. 海洋科學(xué), 2012, 36(7): 70-78. LIU Jun, ZHANG Jiaye, RAN Xiangbin, et al. Characteristics of different forms of nitrogen and phosphorus and organic carbon in the sediments of low-oxygen zones in adjacent Rushan Bay [J]. Marine Sciences, 2012, 36(7): 70-78.

    [3] 楊波, 王保棟, 韋欽勝, 等. 低氧環(huán)境對沉積物中生源要素生物地球化學(xué)循環(huán)的影響[J]. 海洋科學(xué), 2012, 36(5): 124-129. YANG Bo, WANG Baodong, WEI Qinsheng, et al. Impact of hypoxia on biogeochemical cycling of biogenic elements in sediment [J]. Marine Sciences, 2012, 36(5): 124-129.

    [4] 劉霜, 張繼民, 張洪亮, 等. 青島近海赤潮災(zāi)害分級與時空分布及赤潮生物的變化[J]. 水生態(tài)學(xué)雜志, 2014(4): 43-47. LIU Shuang, ZHANG Jimin, ZHANG Hongliang, et al. Spatial-temporal Distribution and Hazard Grading of Red Tide and Variation of Red-tide Species in Qingdao Adjacent Sea Area [J]. Reservoir Fisheries, 2014(4): 43-47.

    [5] 王妍, 董志軍, 劉東艷, 等. 煙臺近海浮游植物的時空變化特征[J]. 海洋通報, 2013, 32(4):408-420. WANG Yan, DONG Zhijun LIU Dongyan, er al.Variation of spatial and temporal distributions of phytoplankton community in coastal waters of Yantai [J]. Marine Science Bulletin, 2013, 32(4): 408-420.

    [6] 喬方利,馬德毅,朱明遠(yuǎn), 等. 2008年黃海滸苔爆發(fā)的基本狀況與科學(xué)應(yīng)對措施[J]. 海洋科學(xué)進展,2008, 26(3):409-410. QIAO Fangli, MA Deyi, ZHU Mingyuan, et al.The scientific countermeasures and state for enteromorpha occurrence of Yellow Sea in 2008 [J]. Advances in Marine Science, 2008,26(3): 409-410.

    [7] Sun S, Wang F, Li C et al. Emerging challenges: Massive green algae blooms in the Yellow Sea[J]. Nature Precedings, 2008, 2266(1): 1-5.

    [8] 王顥, 石曉勇, 張傳松, 等. 2004年春季東海赤潮高發(fā)區(qū)COD分布及其與赤潮關(guān)系的初步研究[J]. 海洋科學(xué), 2008, 32(12): 82-86. WANG Hao, SHI Xiaoyong, ZHANG Chuansong, et al. Primary studies on the COD distribution and relationship with HAB for high frequency HAB occurrence areas of the East China Sea in spring 2004[J]. Marine Sciences, 2008, 32(12): 82-86.

    [9] 中華人民共和國環(huán)境保護部. 近岸海域環(huán)境質(zhì)量公報(2005-2013)[C]. 2006-2014. Chinese Ministry of Environmental Protectio. Environmental quality bulletin of offshore area (2005-2013)[C]. 2006-2014.

    [10] 王修林, 崔正國, 李克強, 等. 渤海COD入海通量估算及其分配容量優(yōu)化研究[J]. 海洋環(huán)境科學(xué), 2009, 28(5):497-500. WANG Xiulin, CUI Zhengguo, LI Keqiang, et al. Study on estimation of flux of COD and allocation capacity in Bohai Sea [J]. Marine Environmental Science, 2009, 28(5): 497-500.

    [11] Theriault E J, Mohlman F W. The oxygen demand of polluted waters [J]. Journal, 1928, 19(6):731-732.

    [12] Skopintsev B A, Ovtshinnikova J S. Study of oxidation processes in the polluted waters under aetobic conditions[J]. Midrobioloclya, 1934, 3: 138-147.

    [13] Krause H R. Biochemische Untersuchungen uber den post-mortalen Abbau vontotem Plankton unter aeroben und anaeroben Bedingungen[J]. Arch Hydrobiol, 1959, Suppl(24): 297-337.

    [14] Skopintsev B A, Krylova L P. The results of some aspects of organic matter dynamics study in the natural waters[J]. Tr Vses Gidrobiol Ova, 1955 (6): 35-45.

    [15] Richards F A, Cline I D, Broenkow W W, er al. Some consequences of the decomposition of organic matter of lake Nitinat, an anoxic fjord. Limnol[J]. Oceanogr, 1965, 10(Suppl.): 185-201.

    [16] Skopintsev B A. Formation of the Present Chemical Composition in the Black Sea Waters[M]. Leningrad: Gidrometeoizdat, 1975.

    [17] Poots V, McKay G, Healy J. Removal of basic dye from effluent using wood as an adsorbent[J]. Journal (Water Pollution Control Federation), 1978, 50(5): 926-935.

    [18] Raji C, Shubha K, Anirudhan T. Use of chemically modified sawdust in the removal of Pb (II) ions from aqueous media[J]. Indian Journal of Environmental Health, 1997, 39(3): 230-238.

    [19] Zeeman G, Sanders W, Wang K Y, et al. Anaerobic treatment of complex wastewater and waste activated sludge—application of an upflow anaerobic solid removal (UASR) reactor for the removal and pre-hydrolysis of suspended COD[J]. Water Science and Technology, 1997, 35(10): 121-128.

    [20] Suffet I H, Brenner L, Coyle J T, et al. Evaluation of the capability of granular activated carbon and XAD-2 resin to remove trace organics from treated drinking water[J]. Environmental Science & Technology, 1978, 12(12): 1315-1322.

    [21] Chen P H. Adsorption of organic compounds in water using a synthetic adsorbent[J]. Environment International, 1997, 23(1): 63-73.

    [22] 韓相奎, 周春生. SBR 法過程動力學(xué)研究[J]. 環(huán)境科學(xué)學(xué)報, 1996, 16(4): 394-399. HAN Xiangkui, ZHOU Chunsheng. Kinetics of SBR process [J]. Journal of Environmental Sciences, 1996, 16(4): 394-399.

    [23] Devi R, Dahiya R. COD and BOD removal from domestic wastewater generated in decentralised sectors[J]. Bioresource technology, 2008, 99(2): 344-349.

    [24] Devi R, Singh V, Kumar A. COD and BOD reduction from coffee processing wastewater using Avacado peel carbon[J]. Bioresource technology, 2008, 99(6): 1853-1860.

    [25] 顧夏聲. 廢水生物處理數(shù)學(xué)模式[M]. 北京: 清華大學(xué)出版社, 1993. GU Xiasheng.Mathematical Model for Biological Treatment of Wastewater [M]. Beijing: Tsinghua University Press, 1993.

    [26] Zhao X X, Wang X L, Shi X Y, et al. Environmental capacity of chemical oxygen demand in the Bohai Sea: Modeling and calculation[J]. Chinese Journal of Oceanology and Limnology, 2011, 1: 46-52.

    [27] 張玉清. 水污染動力學(xué)和水污染控制[M]. 北京: 化學(xué)工業(yè)出版社, 2007. ZHANG Yuqing. Kinetics and Control of Water Pollution[M]. Beijing: Chemical Industry Press, 2007.

    [28] 慕金波, 酒濟明. 河流中有機物降解系數(shù)的室內(nèi)模擬實驗研究[J]. 山東科學(xué), 1997, 10(2): 50-55. MU Jinbo, JIU Jiming. A laboratory simulation study for the degradation coefficient of organic matter in rivers [J]. Shandong Science, 1997, 10(2): 50-55.

    [29] 楊亮秀, 楊常亮, 劉衛(wèi)紅, 等. 滇池流域盤龍江化學(xué)需氧量與總氮降解系數(shù)的測定研究[C]∥. [s.1.]: 2011中國環(huán)境科學(xué)學(xué)會學(xué)術(shù)年會論文集(第一卷). 2011. YANG Liang xiu, YANG Changliang, LIU Weihong, et al. A research for the degradation coefficient of COD and total nitrogen in the Dianchi Lake basin [C]∥. [s.1.]: The Annual Meeting’s Proceedings of Chinese Society of Environmental Science in 2011 (Vol. 1). 2011.

    [30] 慕金波, 韓言柱. 南四湖港口物擴散區(qū) COD 降解系數(shù)的研究[J]. 中國環(huán)境監(jiān)測, 1997, 13(4): 47-50. MU Jinbo, HAN Yanzhu. Study of COD degradation coefficients in the port pollution regions of Nansi Lake [J]. Environmental Monitoring in China, 1997, 13(4): 47-50.

    [31] Leivisk? T, Nurmesniemi H, P?yki?R, et al. Effect of biological wastewater treatment on the molecular weight distribution of solublpounds and on the reduction of BOD, COD and P in pulp and paper mill effluent[J]. Water Research, 2008, 42(14): 3952-3960.

    [32] 陳誦英, 陳平, 李永旺,等. 催化反應(yīng)動力學(xué)[M]. 北京: 化學(xué)工業(yè)出版社, 2007. CHEN Songying, CHEN Ping, LI Yongwang,et al. Catalytic Reaction Kinetics[M]. Beijing: Chemical Industry Press, 2007.

    [33] 祖波, 張代鈞, 周富春. 強化廢水生物脫氮新技術(shù)—厭氧氨氧化、甲烷化、反硝化耦合的機理與動力學(xué)研究[M]. 北京:中國環(huán)境科學(xué)出版社, 2009. ZU Bo, ZHANG Daijun, ZHOU Fuchun. Strengthen the New Technology of Nitrogen’s Biological Removal from Wastewater [M]. Beijing: China Environmental Science Press, 2009.

    [34] Paris D F, Steen W C, Baughman G L, et al. Second-order model to predict microbial degradation of organic compounds in natural waters[J]. Applied and Environmental Microbiology, 1981, 41(3): 603-609.

    [35] Alvarez P J, Anid P J, Vogel T M. Kinetics of aerobic biodegradation of benzene and toluene in sandy aquifer material[J]. Biodegradation, 1991, 2(1): 43-51.

    [36] Reardon K F, Mosteller D C, Rogers J D B. Biodegradation kinetics of benzene, toluene, and phenol as single and mixed substrates for pseudomonas putida F1[J]. Biotechnology and Bioengineering, 2000, 69(4): 385-400.

    [37] 季民, 孫志偉, 王澤良, 等. 污水排海有機物的生化降解動力學(xué)系數(shù)測定及水質(zhì)模擬[J]. 中國給水排水, 1999, 15(11): 62-65. JI Min, SUN Zhiwei, WANG Zeliang, er al. The biochemical degradation kinetic coefficient and water quality simulating of organic matter in wastewater discharged into sea [J].China Water and Wastewater, 1999, 15(11): 62-65.

    [38] 崔有為, 王淑瑩, 甘湘慶, 等. 生物處理含鹽污水的鹽抑制動力學(xué)[J]. 環(huán)境污染治理技術(shù)與設(shè)備, 2005, 6(5): 38-41. CUI Youwei, WANG Shuying, GAN Xiangqing, et al. Salt inhibition kinetics in biological treatment of saline wastewater [J]. Techniques and Equipment for Environmental Pollution Control, 2005, 6 (5): 38-41.

    [39] 宋學(xué)宏, 邴旭文, 孫麗萍, 等. 陽澄湖養(yǎng)殖水體 COD 降解動力學(xué)研究[J]. 安徽農(nóng)業(yè)大學(xué)學(xué)報, 2010(2): 328-332. SONG Xuehong, BING Xuwen, SUN Liping, et al. A research for the degradation kinetics of COD in aquatic water of Yangcheng Lake [J]. Anhui Agricultural University, 2010(2): 328-332.

    [40] 郭棟鵬. 黃海南部海域排海尾水中污染物降解規(guī)律研究[D].太原: 太原理工大學(xué), 2006. GUO Dongpeng. A research for the degradation rule of the pollution in sail water discharged into the south Yellow Sea [D].Taiyuan: Taiyuan University of Technology, 2006.

    [41] 楊東方, 耿曉, 曲延峰等. 膠州灣水體重金屬汞的分布和重力特性[J]. 海洋開發(fā)與管理, 2014(12): 71-77. YANG Dongfang, GENG Xiao, QU Yanfeng, et al. The distribution and gravity property of Hg in Jiaozhou Bay [J]. Development and Management of Ocean, 2014(12): 71-77.

    [42] Cubas S, Foresti E, Rodrigues J, et al. Influence of liquid-phase mass transfer on the performance of a stirred anaerobic sequencing batch reactor containing immobilized biomass[J]. Biochemical Engineering Journal, 2004, 17(2): 99-105.

    [43] Pinho S C, Ratusznei S M, Rodrigues J A D, et al. Influence of the agitation rate on the treatment of partially soluble wastewater in anaerobic sequencing batch biofilm reactor[J]. Water Research, 2004, 38(19): 4117-4124.

    [44] 張世坤, 張建軍, 田依林, 等. 黃河花園口典型污染物自凈降解規(guī)律研究[J]. 人民黃河, 2006, 28(4): 46-47. ZHANG Shikun, ZHANG Jianjun, TIAN Yilin, et al. A research for the self-purification degradation rule of the typical pollution in Yellow River Huayuankou [J]. Yellow River of People, 2006, 28(4): 46-47.

    [45] 鐘巍. 基于 Arrhenius定理的化學(xué)動力學(xué)數(shù)值計算法[J]. 化學(xué)研究, 2011, 22(3): 56-60. ZHONG Wei. Arrhenius theory-based numerical algorithm for chemical kinetics [J]. Chemical Research, 2011, 22(3): 56-60.

    [46] 陳松, 廖文卓, 許愛玉, 等. 廈門市排海污水COD的降解動力學(xué)[J]. 海洋學(xué)報 (中文版), 1991(3): 401-406. CHEN Song, LIAO Wenzhuo, XU Aiyu, et al. A research for degradaion kinetics of COD in the wastewater discharged into sea from Xiamen City [J]. Acta of Oceanologica Sinica (Chinese Version), 1991(3): 401-406.

    [47] Yesilada O, Sik S, Sam M. Biodegradation of olive oil mill wastewaterby Coriolus versicolor and Funalia trogii:effects of agitation, initial COD concentration, inoculum size and immobilization[J]. World Journal of Microbiology & Biotechnology, 1997, 14(14):37-42.

    責(zé)任編輯 徐 環(huán)

    A Research for the Bioavailability and Degradation Kinetics of the Oxygen-ConsumedOrganic Matter from Coastal Sea of Shandong Province

    MA Yun-Peng1, DAI Ai-Quan1,2, LI Ke-Qiang1, WANG Xiu-Lin1

    (1.Key Laboratory of Marine Chemistry Theory and Technology, Ministry of Education, College of Chemistry and Chemical Engineering, Ocean University of China, Qingdao 266100, China; 2.Qingdao Research Academy of Environmental Sciences, Qingdao 266000, China)

    It is difficult to simulate well on the concentration distribution of oxygen-consumed organic matter(characterized by COD)based on a single overall degradation modelwith respect to the bulk organic matterincoastal water,because of the differences in composition and bioavailability of terrigenous COD.Thetypical COD sources in Shandong provincewere chosen to study bioavailability and degradation kinetics of COD by laboratory culture. The results show thatthe degradation of COD basically followed the first-order kinetic law which is fast followed slow and tending to be sable finally; andthe degradation rate constants and bioavailability of different terrigenous COD exist significant differences(P<0.01).The bioavailability of COD fromprinting and dyeing mill source is maximum, about 92%, following by wastewater treatment plant and Xiaoqing river sources, whilstthe brewery and Yu river sources are minimum; thedegradation rate constantsof COD fromthe brewery source is maximum,about 0.40 d-1following by wastewater treatment plant source, whilstXiaoqing river, Yu river and the printing and dyeing mill sources are minimum. Based on the differences of degradation rate constants, COD from different sources can be divided into four types: food industry source’s (brewery) COD, domestic pollution source’s (wastewater treatment plant) COD, river integrated pollution source’s (Xiaoqing river, Yu river and the printing and dyeing mill) COD and refractory COD. The results of the culture experiment of COD from wastewater treatment plant show that the degradation rate constants of COD increase with rising of temperature and stirringrate, and the bioavailability increases with rising of stirring rate. The results of the study can provide support for increasing accuracy of water quality simulations in water environment management.

    degradation kinetics; bioavailability; oxygen-consumed organic matter (COD); degradation rate constant; coastal sea of Shandong

    國家自然科學(xué)基金委員會-山東省人民政府聯(lián)合資助海洋科學(xué)研究中心項目(U1406403)資助 Supported by NSFC-Shandong Joint Fund for Marine Science Research Centers (U1406403)

    2016-03-11;

    2016-05-15

    馬云鵬(1991-),男,碩士生。E-mail:353067464@qq.com

    ?? 通訊作者:E-mail: likeqiang@ouc.edu.cn

    P734.5

    A

    1672-5174(2017)07-079-10

    10.16441/j.cnki.hdxb.20160066

    馬云鵬,戴愛泉,李克強,等.山東近海陸源耗氧有機物生物可利用性及其降解動力學(xué)研究[J].中國海洋大學(xué)學(xué)報(自然科學(xué)版),2017, 47(7): 79-88

    MA Yun-Peng, DAI Ai-Quan, LI Ke-Qiang, et al.A research for the bioavailability and degradation kinetics of the oxygen-consumed organic matter from coastal sea of Shandong Province[J].Periodical of Ocean University of China, 2017, 47(7): 79-88.

    猜你喜歡
    耗氧陸源常數(shù)
    放學(xué)
    關(guān)于Landau常數(shù)和Euler-Mascheroni常數(shù)的漸近展開式以及Stirling級數(shù)的系數(shù)
    心梗時大力咳嗽 能救命?
    蝦夷扇貝對溫度、饑餓脅迫的應(yīng)答機制淺析
    人體代謝耗氧模擬裝置的改進設(shè)計
    陸源有機碳對萊州灣浮游動物能量貢獻的初步研究
    控制陸源污染,保護海洋環(huán)境
    求友聲的變奏
    西湖(2017年5期)2017-05-12 16:40:09
    幾個常數(shù)項級數(shù)的和
    萬有引力常數(shù)的測量
    亚洲国产精品久久男人天堂| 精品国内亚洲2022精品成人| 中亚洲国语对白在线视频| 欧美在线黄色| 国产午夜福利久久久久久| 日韩国内少妇激情av| 国内久久婷婷六月综合欲色啪| 欧美黄色淫秽网站| 这个男人来自地球电影免费观看| 国产一区二区三区在线臀色熟女| 黄色女人牲交| 国产野战对白在线观看| 国产一区在线观看成人免费| 男女下面进入的视频免费午夜| 12—13女人毛片做爰片一| 午夜日韩欧美国产| 久久热在线av| 久久草成人影院| 成人性生交大片免费视频hd| 中亚洲国语对白在线视频| 最近最新中文字幕大全免费视频| 免费看美女性在线毛片视频| 国产一区二区三区在线臀色熟女| 久久久久久久午夜电影| 十八禁网站免费在线| 欧美黄色片欧美黄色片| 最新中文字幕久久久久 | 老司机深夜福利视频在线观看| 窝窝影院91人妻| 色综合欧美亚洲国产小说| 两性夫妻黄色片| 久久亚洲真实| 最近最新免费中文字幕在线| 亚洲av成人av| av欧美777| 此物有八面人人有两片| 老汉色∧v一级毛片| 亚洲午夜精品一区,二区,三区| 国产成人av激情在线播放| 免费看日本二区| 国产v大片淫在线免费观看| 国产精品亚洲av一区麻豆| 国产视频内射| 可以在线观看的亚洲视频| 亚洲五月天丁香| 免费看美女性在线毛片视频| 亚洲av电影在线进入| 久久精品国产亚洲av香蕉五月| 久久精品国产清高在天天线| 亚洲国产看品久久| 精品国产三级普通话版| 人人妻,人人澡人人爽秒播| 国产精品综合久久久久久久免费| 18禁观看日本| 免费av毛片视频| 国产成人系列免费观看| 黄色成人免费大全| 国产精品亚洲一级av第二区| 激情在线观看视频在线高清| 网址你懂的国产日韩在线| 国产亚洲精品综合一区在线观看| 男女那种视频在线观看| 天堂影院成人在线观看| www.精华液| 婷婷丁香在线五月| 在线看三级毛片| 欧美3d第一页| 91麻豆av在线| 免费在线观看成人毛片| 亚洲五月天丁香| 久久精品91无色码中文字幕| 最近最新中文字幕大全免费视频| 亚洲国产欧美人成| 床上黄色一级片| 99精品欧美一区二区三区四区| 国产私拍福利视频在线观看| 少妇丰满av| 亚洲精品一区av在线观看| 伊人久久大香线蕉亚洲五| 身体一侧抽搐| 亚洲成人久久性| 在线十欧美十亚洲十日本专区| 99久久99久久久精品蜜桃| 国产久久久一区二区三区| 久久这里只有精品19| 欧美日韩福利视频一区二区| 中亚洲国语对白在线视频| 国模一区二区三区四区视频 | 熟女人妻精品中文字幕| 欧美3d第一页| svipshipincom国产片| 啦啦啦韩国在线观看视频| 国产熟女xx| 美女被艹到高潮喷水动态| 精品久久久久久久久久久久久| 久久人妻av系列| 国产aⅴ精品一区二区三区波| 香蕉久久夜色| 日韩中文字幕欧美一区二区| 日本与韩国留学比较| 国产不卡一卡二| 一二三四社区在线视频社区8| 男女床上黄色一级片免费看| 免费电影在线观看免费观看| 精品国产美女av久久久久小说| 草草在线视频免费看| 免费电影在线观看免费观看| 久久中文看片网| 男女床上黄色一级片免费看| 中文字幕熟女人妻在线| 亚洲精品美女久久av网站| 国产精品爽爽va在线观看网站| 久久久精品欧美日韩精品| 国产伦人伦偷精品视频| 啪啪无遮挡十八禁网站| 中亚洲国语对白在线视频| 美女cb高潮喷水在线观看 | 真实男女啪啪啪动态图| 国产 一区 欧美 日韩| 国产在线精品亚洲第一网站| 亚洲国产中文字幕在线视频| 丰满人妻熟妇乱又伦精品不卡| 最新美女视频免费是黄的| 两个人看的免费小视频| 午夜久久久久精精品| xxx96com| 看片在线看免费视频| 久久久久久久久中文| 欧美日本视频| 脱女人内裤的视频| 夜夜躁狠狠躁天天躁| 韩国av一区二区三区四区| 久久精品亚洲精品国产色婷小说| 天天一区二区日本电影三级| www日本在线高清视频| 亚洲va日本ⅴa欧美va伊人久久| 成年版毛片免费区| 国语自产精品视频在线第100页| 亚洲av五月六月丁香网| 亚洲欧美日韩东京热| 在线a可以看的网站| 欧美一区二区精品小视频在线| 波多野结衣高清无吗| 日韩成人在线观看一区二区三区| 国产精品久久久人人做人人爽| 亚洲激情在线av| 午夜福利免费观看在线| 制服丝袜大香蕉在线| 91av网站免费观看| 日本免费a在线| 国产精品久久久久久人妻精品电影| 中国美女看黄片| 中文字幕人成人乱码亚洲影| 精品福利观看| 国产精品久久久久久久电影 | 亚洲国产日韩欧美精品在线观看 | 久久久久久久久中文| 国产免费男女视频| 成人国产综合亚洲| 精品不卡国产一区二区三区| 又粗又爽又猛毛片免费看| 久久欧美精品欧美久久欧美| 村上凉子中文字幕在线| 99re在线观看精品视频| 国产av不卡久久| 久久香蕉国产精品| 哪里可以看免费的av片| 琪琪午夜伦伦电影理论片6080| 男女之事视频高清在线观看| 俺也久久电影网| 桃红色精品国产亚洲av| 麻豆成人午夜福利视频| 在线免费观看不下载黄p国产 | 狂野欧美白嫩少妇大欣赏| 精品一区二区三区视频在线 | 狂野欧美白嫩少妇大欣赏| 老司机午夜十八禁免费视频| 日韩欧美在线二视频| 国产人伦9x9x在线观看| 日本三级黄在线观看| 日韩欧美国产在线观看| 国产 一区 欧美 日韩| 国产av麻豆久久久久久久| 日韩国内少妇激情av| 一二三四在线观看免费中文在| 一区二区三区高清视频在线| 国产野战对白在线观看| 久久亚洲精品不卡| 久久久水蜜桃国产精品网| or卡值多少钱| 岛国在线免费视频观看| 欧美成人一区二区免费高清观看 | 免费看光身美女| 一二三四在线观看免费中文在| 成年女人看的毛片在线观看| 大型黄色视频在线免费观看| 一级a爱片免费观看的视频| 99精品欧美一区二区三区四区| 国产精品日韩av在线免费观看| 欧美成人一区二区免费高清观看 | 中文亚洲av片在线观看爽| 日韩免费av在线播放| 国产精品久久久av美女十八| 波多野结衣高清作品| 制服丝袜大香蕉在线| 亚洲精品久久国产高清桃花| 欧美日韩一级在线毛片| 亚洲欧美精品综合久久99| xxx96com| 91在线观看av| 岛国在线观看网站| 人人妻,人人澡人人爽秒播| 91麻豆精品激情在线观看国产| 久久精品综合一区二区三区| 国产精品久久久久久人妻精品电影| 亚洲色图 男人天堂 中文字幕| 国产成人aa在线观看| 老司机福利观看| 中文字幕久久专区| 日本熟妇午夜| 此物有八面人人有两片| 国产免费av片在线观看野外av| 免费高清视频大片| 国产精品久久视频播放| 99精品久久久久人妻精品| 国产激情久久老熟女| 午夜日韩欧美国产| 亚洲中文字幕日韩| 国产免费av片在线观看野外av| 欧美黄色片欧美黄色片| 精品国产乱码久久久久久男人| 欧美成人一区二区免费高清观看 | 国产精品久久久久久久电影 | 久久欧美精品欧美久久欧美| 男人舔女人下体高潮全视频| 九色成人免费人妻av| 日本a在线网址| 美女cb高潮喷水在线观看 | 91字幕亚洲| 午夜a级毛片| 90打野战视频偷拍视频| 中国美女看黄片| 国产成人欧美在线观看| 日韩欧美国产在线观看| 搞女人的毛片| 不卡一级毛片| 成人三级黄色视频| 美女高潮喷水抽搐中文字幕| 国产精品爽爽va在线观看网站| 麻豆av在线久日| 99热6这里只有精品| 欧美日韩亚洲国产一区二区在线观看| 日本在线视频免费播放| 免费在线观看视频国产中文字幕亚洲| 九九在线视频观看精品| 午夜a级毛片| 亚洲九九香蕉| 两性夫妻黄色片| 人妻久久中文字幕网| 中文字幕久久专区| 国产精品久久久久久久电影 | 日本五十路高清| 在线免费观看的www视频| cao死你这个sao货| 免费在线观看成人毛片| 国产高清三级在线| 国语自产精品视频在线第100页| 99国产极品粉嫩在线观看| 久久久成人免费电影| 精品免费久久久久久久清纯| 在线观看免费午夜福利视频| 国内久久婷婷六月综合欲色啪| 老司机深夜福利视频在线观看| 国产三级在线视频| 国产一区二区在线av高清观看| 欧美成人一区二区免费高清观看 | 91字幕亚洲| 久久久久久国产a免费观看| 国产精华一区二区三区| 午夜福利成人在线免费观看| 日日摸夜夜添夜夜添小说| 国产一区二区三区视频了| 国产高清激情床上av| 精品免费久久久久久久清纯| x7x7x7水蜜桃| 99热精品在线国产| 夜夜看夜夜爽夜夜摸| 日本 av在线| 亚洲九九香蕉| 日韩有码中文字幕| 高潮久久久久久久久久久不卡| 久久久久国产精品人妻aⅴ院| av天堂在线播放| 久久精品影院6| 18禁裸乳无遮挡免费网站照片| 一卡2卡三卡四卡精品乱码亚洲| 国产午夜精品论理片| 日本五十路高清| 中文字幕人妻丝袜一区二区| 99久久无色码亚洲精品果冻| 欧美色欧美亚洲另类二区| 真人做人爱边吃奶动态| 999久久久精品免费观看国产| 亚洲激情在线av| 亚洲av熟女| 国产97色在线日韩免费| 亚洲中文日韩欧美视频| 男人舔女人的私密视频| 亚洲18禁久久av| 国产亚洲av嫩草精品影院| 久久午夜综合久久蜜桃| 变态另类成人亚洲欧美熟女| 亚洲片人在线观看| 啪啪无遮挡十八禁网站| 成年免费大片在线观看| 国产av麻豆久久久久久久| 免费在线观看影片大全网站| 岛国视频午夜一区免费看| 色精品久久人妻99蜜桃| 真实男女啪啪啪动态图| 亚洲精品在线美女| 黄片小视频在线播放| 亚洲五月婷婷丁香| 亚洲欧洲精品一区二区精品久久久| 手机成人av网站| 男人舔女人下体高潮全视频| 欧美在线一区亚洲| 美女高潮的动态| 欧美日韩国产亚洲二区| 欧美av亚洲av综合av国产av| 一边摸一边抽搐一进一小说| 精品久久久久久久人妻蜜臀av| 男人舔女人下体高潮全视频| 久久精品91无色码中文字幕| 亚洲在线观看片| 国产爱豆传媒在线观看| 国产高清视频在线观看网站| 毛片女人毛片| 国产欧美日韩精品一区二区| 伊人久久大香线蕉亚洲五| 观看美女的网站| av片东京热男人的天堂| 在线a可以看的网站| 99在线人妻在线中文字幕| xxxwww97欧美| 国产黄片美女视频| 欧美黄色片欧美黄色片| 午夜激情福利司机影院| 国产精品一区二区精品视频观看| 国产av不卡久久| 欧美丝袜亚洲另类 | 国产亚洲精品综合一区在线观看| 亚洲国产精品sss在线观看| 两性夫妻黄色片| 不卡av一区二区三区| 色吧在线观看| 禁无遮挡网站| 怎么达到女性高潮| 欧美av亚洲av综合av国产av| 久久久久九九精品影院| 亚洲午夜理论影院| 热99在线观看视频| tocl精华| 不卡av一区二区三区| 床上黄色一级片| 欧美丝袜亚洲另类 | 在线观看免费视频日本深夜| 男女下面进入的视频免费午夜| 精品久久久久久久末码| 亚洲avbb在线观看| 我的老师免费观看完整版| 久久久久久大精品| e午夜精品久久久久久久| 亚洲国产精品合色在线| 国产欧美日韩一区二区三| 国产欧美日韩一区二区精品| 一级毛片高清免费大全| 无遮挡黄片免费观看| 日韩av在线大香蕉| 一进一出好大好爽视频| 亚洲精品久久国产高清桃花| 久久热在线av| 日韩欧美国产一区二区入口| 18禁美女被吸乳视频| 国产成人av激情在线播放| 99国产极品粉嫩在线观看| 欧美zozozo另类| 高清在线国产一区| 麻豆久久精品国产亚洲av| 国产成人福利小说| 欧美极品一区二区三区四区| 国产精品一区二区精品视频观看| 国产精品九九99| 国产精品亚洲av一区麻豆| 一级黄色大片毛片| 国产一区二区在线av高清观看| 亚洲人成伊人成综合网2020| 特级一级黄色大片| 亚洲美女视频黄频| 这个男人来自地球电影免费观看| 精品国产乱子伦一区二区三区| 国产精品久久久久久亚洲av鲁大| а√天堂www在线а√下载| 99久久成人亚洲精品观看| 国产精品综合久久久久久久免费| 18禁美女被吸乳视频| xxx96com| 国产成人系列免费观看| 成人性生交大片免费视频hd| 亚洲一区二区三区不卡视频| 在线永久观看黄色视频| x7x7x7水蜜桃| 欧美日韩精品网址| 久久久久久国产a免费观看| 九色国产91popny在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成av人片免费观看| 网址你懂的国产日韩在线| 综合色av麻豆| 亚洲人成网站高清观看| 国产不卡一卡二| 午夜激情福利司机影院| 身体一侧抽搐| 日本一二三区视频观看| 国产视频一区二区在线看| 亚洲国产精品成人综合色| 欧美乱码精品一区二区三区| 男女床上黄色一级片免费看| 亚洲avbb在线观看| 亚洲国产精品999在线| 可以在线观看的亚洲视频| 网址你懂的国产日韩在线| 噜噜噜噜噜久久久久久91| 国产在线精品亚洲第一网站| 亚洲av美国av| 中文字幕av在线有码专区| 叶爱在线成人免费视频播放| 欧美日韩精品网址| 一本一本综合久久| 亚洲真实伦在线观看| 中出人妻视频一区二区| 91麻豆精品激情在线观看国产| 人妻丰满熟妇av一区二区三区| 看片在线看免费视频| 精品国产三级普通话版| 99热这里只有是精品50| 色综合站精品国产| 免费看光身美女| 成人无遮挡网站| 两个人的视频大全免费| 国产69精品久久久久777片 | 久久久国产成人精品二区| 中文字幕人成人乱码亚洲影| 九色国产91popny在线| 国产一区二区激情短视频| 一本综合久久免费| 麻豆av在线久日| 两个人的视频大全免费| 国产精品亚洲av一区麻豆| 久久这里只有精品中国| 免费观看精品视频网站| 99国产综合亚洲精品| 亚洲av五月六月丁香网| 国产欧美日韩一区二区精品| 无遮挡黄片免费观看| 国产精品永久免费网站| 国产精品九九99| 亚洲自偷自拍图片 自拍| 亚洲五月天丁香| 1000部很黄的大片| 久久精品国产99精品国产亚洲性色| 国产免费av片在线观看野外av| 日韩大尺度精品在线看网址| 91av网站免费观看| 中文亚洲av片在线观看爽| 国产成人aa在线观看| 搡老岳熟女国产| 欧美大码av| 欧美性猛交黑人性爽| 国内精品一区二区在线观看| 午夜福利视频1000在线观看| 国产av在哪里看| 欧美中文综合在线视频| 床上黄色一级片| 欧美绝顶高潮抽搐喷水| 免费观看的影片在线观看| 久久欧美精品欧美久久欧美| 亚洲精品乱码久久久v下载方式 | 久9热在线精品视频| 亚洲乱码一区二区免费版| svipshipincom国产片| 亚洲熟妇熟女久久| 两性午夜刺激爽爽歪歪视频在线观看| 男人的好看免费观看在线视频| 国产三级中文精品| 欧美中文综合在线视频| 麻豆国产97在线/欧美| 天天一区二区日本电影三级| 午夜精品在线福利| 亚洲国产精品sss在线观看| 中文在线观看免费www的网站| 亚洲av日韩精品久久久久久密| 搡老熟女国产l中国老女人| 听说在线观看完整版免费高清| 国产成人精品久久二区二区免费| 亚洲精品国产精品久久久不卡| 啦啦啦观看免费观看视频高清| 国产91精品成人一区二区三区| 精品午夜福利视频在线观看一区| 久久这里只有精品中国| 国产伦精品一区二区三区四那| 欧美成人性av电影在线观看| 精品国产超薄肉色丝袜足j| 国产视频内射| 免费电影在线观看免费观看| 在线观看舔阴道视频| 欧美xxxx黑人xx丫x性爽| 老熟妇仑乱视频hdxx| 欧美一级毛片孕妇| 免费高清视频大片| 国内精品一区二区在线观看| 麻豆国产av国片精品| 一级a爱片免费观看的视频| 久久久久久久久免费视频了| 在线观看美女被高潮喷水网站 | 中文在线观看免费www的网站| 噜噜噜噜噜久久久久久91| av福利片在线观看| 高清毛片免费观看视频网站| 亚洲国产高清在线一区二区三| 美女cb高潮喷水在线观看 | 亚洲精华国产精华精| 亚洲欧美精品综合一区二区三区| 成人无遮挡网站| 动漫黄色视频在线观看| 欧美成狂野欧美在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲人与动物交配视频| 久久久久久人人人人人| 桃红色精品国产亚洲av| 成人一区二区视频在线观看| 亚洲,欧美精品.| 最近最新中文字幕大全电影3| 精品国产三级普通话版| 天堂影院成人在线观看| 久久久成人免费电影| 亚洲男人的天堂狠狠| 九九久久精品国产亚洲av麻豆 | 精品日产1卡2卡| 国产成人系列免费观看| 亚洲av电影不卡..在线观看| 久久亚洲精品不卡| 国产精品永久免费网站| 久久亚洲精品不卡| 欧美国产日韩亚洲一区| 亚洲av美国av| 99久久精品一区二区三区| 亚洲国产欧美人成| 免费一级毛片在线播放高清视频| 中文字幕人妻丝袜一区二区| 日韩欧美国产在线观看| 最新中文字幕久久久久 | 中文字幕av在线有码专区| 哪里可以看免费的av片| av黄色大香蕉| 国产蜜桃级精品一区二区三区| 色综合欧美亚洲国产小说| 久久国产乱子伦精品免费另类| 亚洲欧美激情综合另类| 国产乱人视频| 精品一区二区三区视频在线 | 亚洲欧美日韩东京热| 成人午夜高清在线视频| 在线十欧美十亚洲十日本专区| 观看美女的网站| 不卡av一区二区三区| 欧美成人性av电影在线观看| 搞女人的毛片| 91在线精品国自产拍蜜月 | 国产精品野战在线观看| АⅤ资源中文在线天堂| 欧美最黄视频在线播放免费| 欧美日韩中文字幕国产精品一区二区三区| 国产成人av激情在线播放| 99精品欧美一区二区三区四区| 国内精品久久久久精免费| 天堂动漫精品| 亚洲在线自拍视频| 国产欧美日韩精品亚洲av| 亚洲 欧美一区二区三区| 色综合站精品国产| 国产成人精品久久二区二区免费| 国产精品自产拍在线观看55亚洲| 又黄又粗又硬又大视频| 亚洲国产精品成人综合色| 一级黄色大片毛片| 国产欧美日韩一区二区精品| 久久久久亚洲av毛片大全| 欧美成人免费av一区二区三区| 国产成人系列免费观看| www.熟女人妻精品国产| 床上黄色一级片| 又黄又粗又硬又大视频| www.www免费av| 亚洲专区中文字幕在线| 99久国产av精品| 国产精品,欧美在线| 看片在线看免费视频| а√天堂www在线а√下载| 国产精品99久久99久久久不卡| 亚洲色图 男人天堂 中文字幕| 男人舔女人的私密视频| 少妇的丰满在线观看| 国产伦精品一区二区三区视频9 | 在线观看午夜福利视频|