• <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ù)的測量
    中文欧美无线码| 女人被狂操c到高潮| 精品一区二区三卡| 1000部很黄的大片| 日韩三级伦理在线观看| 美女cb高潮喷水在线观看| 18禁动态无遮挡网站| 国产精品不卡视频一区二区| 亚洲av中文av极速乱| 国产精品一区二区三区四区久久| 99久久精品国产国产毛片| 久久精品国产亚洲av天美| 亚洲欧美日韩东京热| 亚洲激情五月婷婷啪啪| 久久精品综合一区二区三区| 麻豆精品久久久久久蜜桃| 亚洲三级黄色毛片| 高清午夜精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 国产视频内射| 搞女人的毛片| 精品久久久久久久久亚洲| 青青草视频在线视频观看| 亚洲精品成人av观看孕妇| 免费少妇av软件| 99视频精品全部免费 在线| 男人舔女人下体高潮全视频| 久久午夜福利片| 国产v大片淫在线免费观看| 中文资源天堂在线| 国产 一区精品| 18禁在线播放成人免费| 欧美成人精品欧美一级黄| 国产成人一区二区在线| 五月伊人婷婷丁香| 美女内射精品一级片tv| 天堂影院成人在线观看| 国产av码专区亚洲av| 亚洲怡红院男人天堂| 国产精品国产三级专区第一集| 日本猛色少妇xxxxx猛交久久| 国产精品嫩草影院av在线观看| 最近最新中文字幕大全电影3| 欧美区成人在线视频| 菩萨蛮人人尽说江南好唐韦庄| 少妇猛男粗大的猛烈进出视频 | 久久鲁丝午夜福利片| 国产精品福利在线免费观看| 啦啦啦中文免费视频观看日本| 国产中年淑女户外野战色| av专区在线播放| 麻豆久久精品国产亚洲av| 爱豆传媒免费全集在线观看| 春色校园在线视频观看| 又黄又爽又刺激的免费视频.| 九九在线视频观看精品| 天堂√8在线中文| 中文乱码字字幕精品一区二区三区 | 97人妻精品一区二区三区麻豆| a级毛片免费高清观看在线播放| 国产成人aa在线观看| 日韩欧美一区视频在线观看 | 91精品国产九色| 免费看光身美女| 大话2 男鬼变身卡| 黄片无遮挡物在线观看| 少妇被粗大猛烈的视频| 精品久久久久久电影网| 国产欧美另类精品又又久久亚洲欧美| 久久久精品94久久精品| 欧美zozozo另类| 在线 av 中文字幕| 国产精品蜜桃在线观看| 久久97久久精品| 在线 av 中文字幕| 亚洲人成网站在线播| 赤兔流量卡办理| 亚洲电影在线观看av| 少妇人妻一区二区三区视频| 亚洲国产欧美在线一区| av女优亚洲男人天堂| 精品国内亚洲2022精品成人| 街头女战士在线观看网站| 偷拍熟女少妇极品色| 国产一级毛片在线| 久久久久国产网址| 五月玫瑰六月丁香| 人妻系列 视频| 色5月婷婷丁香| 亚洲欧美日韩无卡精品| videossex国产| 天堂影院成人在线观看| 日韩av在线免费看完整版不卡| 国产精品爽爽va在线观看网站| 亚洲精品成人av观看孕妇| 尾随美女入室| 美女国产视频在线观看| 最近2019中文字幕mv第一页| 午夜免费男女啪啪视频观看| 热99在线观看视频| 成年人午夜在线观看视频 | 亚洲伊人久久精品综合| 国产人妻一区二区三区在| 国产在线男女| 男女下面进入的视频免费午夜| 亚洲婷婷狠狠爱综合网| 精品久久久久久久人妻蜜臀av| 国产极品天堂在线| 精品久久久久久久末码| 久久久精品免费免费高清| 神马国产精品三级电影在线观看| 六月丁香七月| 波多野结衣巨乳人妻| 91aial.com中文字幕在线观看| 日日干狠狠操夜夜爽| 亚洲图色成人| 在线天堂最新版资源| 中文精品一卡2卡3卡4更新| 久99久视频精品免费| 日日摸夜夜添夜夜添av毛片| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲av嫩草精品影院| 亚洲三级黄色毛片| 日本与韩国留学比较| 亚洲aⅴ乱码一区二区在线播放| 日本三级黄在线观看| 欧美日韩综合久久久久久| 中文精品一卡2卡3卡4更新| 一级黄片播放器| 伊人久久精品亚洲午夜| 亚洲国产欧美在线一区| 少妇人妻精品综合一区二区| 一个人看视频在线观看www免费| 国产伦精品一区二区三区视频9| 日本一二三区视频观看| 国产在线男女| 国产成人91sexporn| 成年人午夜在线观看视频 | 精品久久久久久久久av| 中文字幕制服av| 国产老妇女一区| 高清av免费在线| 日韩一区二区三区影片| 美女脱内裤让男人舔精品视频| 亚洲熟妇中文字幕五十中出| 啦啦啦啦在线视频资源| 九九在线视频观看精品| 少妇裸体淫交视频免费看高清| 亚洲av电影不卡..在线观看| 国产午夜福利久久久久久| 97在线视频观看| 婷婷色av中文字幕| 亚洲欧美精品自产自拍| 日本免费a在线| 极品少妇高潮喷水抽搐| 国产黄色免费在线视频| 亚洲国产最新在线播放| 亚洲精华国产精华液的使用体验| 亚洲综合色惰| 亚洲一区高清亚洲精品| 日本-黄色视频高清免费观看| 国产视频首页在线观看| 国产乱人偷精品视频| 成人毛片60女人毛片免费| 舔av片在线| 久久精品国产亚洲av天美| av在线蜜桃| 欧美性猛交╳xxx乱大交人| 国产久久久一区二区三区| 国产免费福利视频在线观看| 舔av片在线| 99久久人妻综合| 人妻夜夜爽99麻豆av| 91午夜精品亚洲一区二区三区| 舔av片在线| 日本三级黄在线观看| 边亲边吃奶的免费视频| 六月丁香七月| av在线亚洲专区| 天堂俺去俺来也www色官网 | 91久久精品国产一区二区成人| 又黄又爽又刺激的免费视频.| 国产极品天堂在线| 日本一二三区视频观看| 国产伦理片在线播放av一区| 国产精品一区二区三区四区久久| 91精品国产九色| 亚洲第一区二区三区不卡| 综合色av麻豆| 国产亚洲av片在线观看秒播厂 | 久久久久网色| 网址你懂的国产日韩在线| 天堂√8在线中文| 免费看日本二区| 国产av码专区亚洲av| 搡老乐熟女国产| 亚洲不卡免费看| 中文天堂在线官网| 欧美性感艳星| 大又大粗又爽又黄少妇毛片口| 免费看a级黄色片| 精品久久久久久久末码| 欧美激情国产日韩精品一区| 亚洲欧洲日产国产| 久久久a久久爽久久v久久| 大话2 男鬼变身卡| 国产国拍精品亚洲av在线观看| 亚洲av电影在线观看一区二区三区 | 一级片'在线观看视频| 最近手机中文字幕大全| 欧美高清性xxxxhd video| 亚洲欧美日韩东京热| 免费看a级黄色片| 听说在线观看完整版免费高清| 免费在线观看成人毛片| 午夜免费男女啪啪视频观看| 国产永久视频网站| 欧美精品国产亚洲| 亚洲四区av| 精品久久久噜噜| 国产免费福利视频在线观看| 免费电影在线观看免费观看| 国产人妻一区二区三区在| 岛国毛片在线播放| 精品一区在线观看国产| 一级爰片在线观看| 搡女人真爽免费视频火全软件| 能在线免费观看的黄片| 99视频精品全部免费 在线| 一级二级三级毛片免费看| 国产 一区 欧美 日韩| 亚洲自拍偷在线| 男人和女人高潮做爰伦理| 精品少妇黑人巨大在线播放| 综合色av麻豆| 亚洲三级黄色毛片| 国产 亚洲一区二区三区 | 最近2019中文字幕mv第一页| 亚洲第一区二区三区不卡| 最近中文字幕高清免费大全6| 精品一区二区三卡| 日韩大片免费观看网站| 久久草成人影院| 国产精品综合久久久久久久免费| 午夜免费激情av| av在线亚洲专区| 国产精品一区二区三区四区久久| 午夜久久久久精精品| 国产三级在线视频| 国产精品一区www在线观看| 亚洲经典国产精华液单| 最近中文字幕高清免费大全6| 成人欧美大片| av免费在线看不卡| 亚洲人成网站在线播| 国产黄色免费在线视频| 日韩亚洲欧美综合| 波野结衣二区三区在线| 最近视频中文字幕2019在线8| 国产中年淑女户外野战色| 老司机影院毛片| 国产一区亚洲一区在线观看| 午夜激情福利司机影院| 插逼视频在线观看| 美女内射精品一级片tv| 国产成人精品福利久久| 综合色丁香网| 久久亚洲国产成人精品v| 国产男人的电影天堂91| 男女那种视频在线观看| 亚洲伊人久久精品综合| 高清毛片免费看| 国产亚洲5aaaaa淫片| 一本一本综合久久| 成年女人看的毛片在线观看| 亚洲va在线va天堂va国产| 特级一级黄色大片| 最近视频中文字幕2019在线8| 日韩视频在线欧美| 最近最新中文字幕免费大全7| 蜜桃亚洲精品一区二区三区| 久久99精品国语久久久| 久久精品夜色国产| 天天躁夜夜躁狠狠久久av| 国语对白做爰xxxⅹ性视频网站| 免费大片黄手机在线观看| 国产高清国产精品国产三级 | 高清在线视频一区二区三区| 国内揄拍国产精品人妻在线| 中文字幕制服av| av在线观看视频网站免费| 精品欧美国产一区二区三| 国产黄片美女视频| 日韩av在线大香蕉| 日本熟妇午夜| 国产黄色小视频在线观看| 看黄色毛片网站| 亚洲自偷自拍三级| 国产日韩欧美在线精品| 国产成人免费观看mmmm| 内地一区二区视频在线| 哪个播放器可以免费观看大片| 久久精品国产亚洲网站| 国产极品天堂在线| 日韩一区二区三区影片| 国产乱人视频| 日韩 亚洲 欧美在线| 欧美三级亚洲精品| 插阴视频在线观看视频| 精品国产一区二区三区久久久樱花 | 高清午夜精品一区二区三区| 18禁裸乳无遮挡免费网站照片| 女人久久www免费人成看片| 99九九线精品视频在线观看视频| 欧美日本视频| 成人漫画全彩无遮挡| 午夜视频国产福利| 久久久久久久午夜电影| 日韩欧美一区视频在线观看 | 日韩国内少妇激情av| 欧美另类一区| 免费看不卡的av| 精品久久国产蜜桃| 80岁老熟妇乱子伦牲交| 色综合站精品国产| a级一级毛片免费在线观看| 国产男人的电影天堂91| 青春草视频在线免费观看| 亚洲一级一片aⅴ在线观看| 亚洲av在线观看美女高潮| 黄片wwwwww| 欧美丝袜亚洲另类| 成人综合一区亚洲| 精品午夜福利在线看| 亚洲av在线观看美女高潮| 日韩中字成人| 亚洲久久久久久中文字幕| 久久久午夜欧美精品| 卡戴珊不雅视频在线播放| 国产成人午夜福利电影在线观看| 久久精品综合一区二区三区| 久久精品国产鲁丝片午夜精品| 男女边摸边吃奶| 亚洲国产色片| 国产一级毛片七仙女欲春2| 嫩草影院入口| 九九在线视频观看精品| 婷婷色综合www| 熟女电影av网| 一本一本综合久久| 成人毛片a级毛片在线播放| 国内精品宾馆在线| 永久网站在线| 简卡轻食公司| 午夜日本视频在线| 国产亚洲一区二区精品| 2022亚洲国产成人精品| 色视频www国产| 久久久午夜欧美精品| 亚洲人与动物交配视频| 国产精品国产三级专区第一集| 免费黄色在线免费观看| 神马国产精品三级电影在线观看| 色视频www国产| 国产三级在线视频| 99久国产av精品| 精品人妻视频免费看| 亚洲精品中文字幕在线视频 | 亚洲最大成人手机在线| 永久免费av网站大全| 一级毛片我不卡| 精品少妇黑人巨大在线播放| 又粗又硬又长又爽又黄的视频| a级毛色黄片| 大陆偷拍与自拍| 精品久久国产蜜桃| 婷婷色av中文字幕| 肉色欧美久久久久久久蜜桃 | 国产黄片视频在线免费观看| 欧美97在线视频| 美女主播在线视频| 亚洲四区av| 亚洲欧美日韩东京热| 高清av免费在线| 亚洲国产成人一精品久久久| 午夜福利网站1000一区二区三区| 午夜亚洲福利在线播放| 天堂影院成人在线观看| 国产午夜精品久久久久久一区二区三区| 亚洲欧美中文字幕日韩二区| 久久久午夜欧美精品| 国产精品国产三级国产专区5o| 人人妻人人澡欧美一区二区| 亚洲成人一二三区av| 国产 一区精品| 中文字幕亚洲精品专区| 九九久久精品国产亚洲av麻豆| 精品久久久久久久久久久久久| 国产精品av视频在线免费观看| 日日啪夜夜爽| 久久久久久久久久黄片| 欧美潮喷喷水| 狂野欧美激情性xxxx在线观看| 尾随美女入室| 日日摸夜夜添夜夜爱| 亚洲成色77777| 国语对白做爰xxxⅹ性视频网站| av.在线天堂| 老司机影院成人| 午夜福利高清视频| 日本黄色片子视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲国产精品国产精品| 亚洲18禁久久av| 国产欧美日韩精品一区二区| 色综合站精品国产| 99re6热这里在线精品视频| 久久6这里有精品| 老司机影院毛片| 毛片女人毛片| 国产黄色免费在线视频| 晚上一个人看的免费电影| 免费播放大片免费观看视频在线观看| 最近最新中文字幕大全电影3| 国产亚洲91精品色在线| 免费av观看视频| 欧美日韩视频高清一区二区三区二| 看免费成人av毛片| 亚洲精品国产av成人精品| 国产免费一级a男人的天堂| 欧美xxxx性猛交bbbb| 国产精品福利在线免费观看| 内地一区二区视频在线| 国产成年人精品一区二区| 国产探花在线观看一区二区| 99九九线精品视频在线观看视频| 成年版毛片免费区| 2021少妇久久久久久久久久久| 午夜福利高清视频| 成人二区视频| 国产成人免费观看mmmm| 成人一区二区视频在线观看| 一级二级三级毛片免费看| 欧美激情久久久久久爽电影| 日本一二三区视频观看| 午夜日本视频在线| 欧美极品一区二区三区四区| 亚洲av免费在线观看| 99久久精品一区二区三区| 深爱激情五月婷婷| 免费av不卡在线播放| 一本一本综合久久| 2021天堂中文幕一二区在线观| 亚洲精品国产av蜜桃| 亚洲一级一片aⅴ在线观看| av专区在线播放| 超碰av人人做人人爽久久| 欧美不卡视频在线免费观看| 国产精品嫩草影院av在线观看| 国产黄色免费在线视频| 亚洲人成网站在线播| 日韩欧美三级三区| 亚洲精华国产精华液的使用体验| 国产一区二区亚洲精品在线观看| av专区在线播放| 免费大片18禁| 三级男女做爰猛烈吃奶摸视频| 国产亚洲精品久久久com| 男女边摸边吃奶| 极品教师在线视频| 爱豆传媒免费全集在线观看| 18+在线观看网站| 女人十人毛片免费观看3o分钟| 你懂的网址亚洲精品在线观看| 精品人妻熟女av久视频| eeuss影院久久| 十八禁国产超污无遮挡网站| 舔av片在线| 久久久久久久久久久免费av| 国产淫片久久久久久久久| 只有这里有精品99| 久久精品国产亚洲av天美| 国产成人精品婷婷| 亚洲在久久综合| 国语对白做爰xxxⅹ性视频网站| 精品人妻偷拍中文字幕| 国产一级毛片在线| 亚洲精品乱久久久久久| 国产亚洲精品av在线| 国产成人a∨麻豆精品| 麻豆久久精品国产亚洲av| 麻豆av噜噜一区二区三区| 日韩不卡一区二区三区视频在线| 精品久久久噜噜| 99九九线精品视频在线观看视频| 三级男女做爰猛烈吃奶摸视频| 久99久视频精品免费| 哪个播放器可以免费观看大片| 国产亚洲av片在线观看秒播厂 | 午夜精品一区二区三区免费看| 舔av片在线| 波野结衣二区三区在线| 欧美xxxx性猛交bbbb| 欧美bdsm另类| 啦啦啦韩国在线观看视频| 在线观看美女被高潮喷水网站| 三级男女做爰猛烈吃奶摸视频| 国产美女午夜福利| 又大又黄又爽视频免费| 精品国内亚洲2022精品成人| 日日干狠狠操夜夜爽| 国产免费福利视频在线观看| 又爽又黄无遮挡网站| 亚洲一级一片aⅴ在线观看| 91精品一卡2卡3卡4卡| 久久午夜福利片| 91午夜精品亚洲一区二区三区| 国产91av在线免费观看| 又粗又硬又长又爽又黄的视频| 精品少妇黑人巨大在线播放| 亚洲四区av| 国产一区二区三区av在线| 老师上课跳d突然被开到最大视频| 淫秽高清视频在线观看| 色网站视频免费| 51国产日韩欧美| 少妇丰满av| 一级毛片aaaaaa免费看小| av在线播放精品| 国产免费福利视频在线观看| 爱豆传媒免费全集在线观看| 欧美+日韩+精品| 国产亚洲av嫩草精品影院| 亚洲欧美精品专区久久| 国产精品麻豆人妻色哟哟久久 | 高清午夜精品一区二区三区| 色视频www国产| 97超视频在线观看视频| 欧美+日韩+精品| 毛片女人毛片| av国产免费在线观看| 国产成人精品一,二区| 搡老妇女老女人老熟妇| 两个人的视频大全免费| 久久久久久久国产电影| 国产精品日韩av在线免费观看| 久久久久久久久久久免费av| 久久精品夜色国产| 亚洲综合色惰| 亚洲性久久影院| 高清欧美精品videossex| 自拍偷自拍亚洲精品老妇| av在线老鸭窝| 久久精品久久精品一区二区三区| 纵有疾风起免费观看全集完整版 | 久久99蜜桃精品久久| 男女那种视频在线观看| 日本爱情动作片www.在线观看| 成人毛片60女人毛片免费| 狠狠精品人妻久久久久久综合| 精品久久久久久久久av| 黄色日韩在线| 99热全是精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 九草在线视频观看| 熟女人妻精品中文字幕| 亚洲伊人久久精品综合| 国产伦在线观看视频一区| 亚洲国产最新在线播放| 亚洲成人精品中文字幕电影| 免费观看精品视频网站| 毛片一级片免费看久久久久| 亚洲av中文av极速乱| 在线观看人妻少妇| 国产黄片视频在线免费观看| av网站免费在线观看视频 | 免费观看a级毛片全部| 成人一区二区视频在线观看| 嫩草影院入口| 国产三级在线视频| 国产色婷婷99| 六月丁香七月| 黄片无遮挡物在线观看| 99久久人妻综合| 天堂√8在线中文| 亚洲性久久影院| 欧美成人a在线观看| 人体艺术视频欧美日本| 国产精品嫩草影院av在线观看| 夜夜看夜夜爽夜夜摸| 亚洲国产色片| 只有这里有精品99| 成人亚洲精品av一区二区| av免费在线看不卡| 可以在线观看毛片的网站| 中国国产av一级| 伦理电影大哥的女人| 毛片女人毛片| 午夜福利在线在线| 成人亚洲精品av一区二区| 日韩欧美三级三区| 中文字幕av在线有码专区| 日本免费a在线| 99视频精品全部免费 在线| 亚洲成人久久爱视频| 91aial.com中文字幕在线观看| 日本三级黄在线观看| 国产亚洲91精品色在线| 久久久久免费精品人妻一区二区| 久久韩国三级中文字幕| 精品一区二区三区人妻视频| 亚洲精品成人av观看孕妇| 男人舔女人下体高潮全视频| 天堂网av新在线| 一级二级三级毛片免费看|