• <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ù)的測量
    美女午夜性视频免费| 蜜桃久久精品国产亚洲av| 国产精品电影一区二区三区| 日本黄色视频三级网站网址| 国产91精品成人一区二区三区| 精品免费久久久久久久清纯| 中出人妻视频一区二区| 天天一区二区日本电影三级| 男人舔女人的私密视频| 老汉色av国产亚洲站长工具| 国产精品美女特级片免费视频播放器 | 两性午夜刺激爽爽歪歪视频在线观看| 日韩欧美 国产精品| 国产91精品成人一区二区三区| 一个人观看的视频www高清免费观看 | 日本熟妇午夜| 国产不卡一卡二| 成在线人永久免费视频| 免费在线观看亚洲国产| 亚洲一区二区三区不卡视频| 少妇人妻一区二区三区视频| 亚洲av日韩精品久久久久久密| 久久午夜亚洲精品久久| 窝窝影院91人妻| 999精品在线视频| 色吧在线观看| 欧美成狂野欧美在线观看| 午夜久久久久精精品| 成人高潮视频无遮挡免费网站| av天堂在线播放| 亚洲中文av在线| 精品一区二区三区视频在线 | 全区人妻精品视频| 一二三四社区在线视频社区8| 一边摸一边抽搐一进一小说| 真实男女啪啪啪动态图| 日本黄色视频三级网站网址| 久久香蕉国产精品| 亚洲欧美激情综合另类| 99国产精品一区二区三区| 日本五十路高清| 久久久久久大精品| 亚洲人与动物交配视频| 中文资源天堂在线| 巨乳人妻的诱惑在线观看| 免费大片18禁| 欧美色欧美亚洲另类二区| 人人妻人人看人人澡| 伊人久久大香线蕉亚洲五| 亚洲精品粉嫩美女一区| 午夜激情福利司机影院| 国产精品野战在线观看| 午夜福利高清视频| 国产精品精品国产色婷婷| 一进一出好大好爽视频| 99riav亚洲国产免费| 丰满人妻熟妇乱又伦精品不卡| 老司机深夜福利视频在线观看| 国产野战对白在线观看| 最近最新中文字幕大全免费视频| 欧美国产日韩亚洲一区| 国产极品精品免费视频能看的| 18禁国产床啪视频网站| 99精品久久久久人妻精品| 啦啦啦观看免费观看视频高清| 99视频精品全部免费 在线 | 欧美乱码精品一区二区三区| 国产三级在线视频| 在线观看美女被高潮喷水网站 | 性色av乱码一区二区三区2| 99久久精品国产亚洲精品| 成人精品一区二区免费| 欧美又色又爽又黄视频| 搡老妇女老女人老熟妇| 在线播放国产精品三级| 人妻夜夜爽99麻豆av| АⅤ资源中文在线天堂| 中文字幕精品亚洲无线码一区| 成人一区二区视频在线观看| 国产av在哪里看| 亚洲性夜色夜夜综合| xxxwww97欧美| av天堂中文字幕网| 国产高清视频在线观看网站| 成人av在线播放网站| 一级毛片精品| 亚洲欧美日韩无卡精品| 老鸭窝网址在线观看| 亚洲精品乱码久久久v下载方式 | 长腿黑丝高跟| 日本撒尿小便嘘嘘汇集6| 香蕉久久夜色| 国产主播在线观看一区二区| 精品欧美国产一区二区三| 丁香六月欧美| av国产免费在线观看| 岛国在线免费视频观看| 国产aⅴ精品一区二区三区波| 国产精品久久久久久久电影 | 成人av在线播放网站| 午夜日韩欧美国产| 国产成人系列免费观看| 在线观看午夜福利视频| 国产欧美日韩精品一区二区| 精品人妻1区二区| 香蕉丝袜av| 色av中文字幕| 国产午夜福利久久久久久| 长腿黑丝高跟| aaaaa片日本免费| 一本综合久久免费| 国产主播在线观看一区二区| 欧美色视频一区免费| 天堂av国产一区二区熟女人妻| 国产av在哪里看| 亚洲欧美一区二区三区黑人| 制服人妻中文乱码| 国产一级毛片七仙女欲春2| 免费高清视频大片| 成年免费大片在线观看| 国内精品一区二区在线观看| 国产精品一区二区三区四区久久| 很黄的视频免费| 亚洲专区字幕在线| 草草在线视频免费看| 99久久无色码亚洲精品果冻| 国产aⅴ精品一区二区三区波| 久久香蕉精品热| 久久天堂一区二区三区四区| 精品一区二区三区视频在线 | 日韩精品中文字幕看吧| 欧美午夜高清在线| 午夜精品在线福利| 在线观看舔阴道视频| 国产av在哪里看| 欧美日韩乱码在线| 国产高清videossex| 在线a可以看的网站| 久久99热这里只有精品18| 又大又爽又粗| 精品久久久久久久人妻蜜臀av| 欧美+亚洲+日韩+国产| 午夜福利免费观看在线| 午夜福利成人在线免费观看| 亚洲中文字幕日韩| 久久久久久久午夜电影| 狂野欧美激情性xxxx| 亚洲午夜精品一区,二区,三区| 巨乳人妻的诱惑在线观看| 一个人观看的视频www高清免费观看 | 精品久久久久久久久久久久久| 又紧又爽又黄一区二区| 久久久久久国产a免费观看| 国产欧美日韩一区二区精品| 99国产精品99久久久久| 精品不卡国产一区二区三区| 国产成人系列免费观看| 嫁个100分男人电影在线观看| 精品一区二区三区av网在线观看| 日韩欧美 国产精品| 精品久久久久久久毛片微露脸| 少妇的丰满在线观看| 亚洲最大成人中文| 亚洲,欧美精品.| 别揉我奶头~嗯~啊~动态视频| 亚洲中文日韩欧美视频| 国产精品久久久人人做人人爽| 亚洲欧美日韩高清专用| 国产毛片a区久久久久| 88av欧美| 成年女人毛片免费观看观看9| 变态另类成人亚洲欧美熟女| 久久久精品欧美日韩精品| 特大巨黑吊av在线直播| 黄色 视频免费看| 我要搜黄色片| 人妻丰满熟妇av一区二区三区| av天堂在线播放| 母亲3免费完整高清在线观看| 久久久久久久久免费视频了| 美女午夜性视频免费| 欧美乱妇无乱码| 久久久久免费精品人妻一区二区| 久久婷婷人人爽人人干人人爱| 一级作爱视频免费观看| 国产精品99久久99久久久不卡| 九色成人免费人妻av| 少妇的丰满在线观看| 熟女电影av网| 一二三四社区在线视频社区8| 神马国产精品三级电影在线观看| 三级国产精品欧美在线观看 | 后天国语完整版免费观看| 久久草成人影院| 国产伦在线观看视频一区| 最好的美女福利视频网| 欧美日韩中文字幕国产精品一区二区三区| 操出白浆在线播放| 小蜜桃在线观看免费完整版高清| 久久久久亚洲av毛片大全| 丁香欧美五月| 日韩av在线大香蕉| 波多野结衣巨乳人妻| 久久亚洲精品不卡| 欧美成狂野欧美在线观看| 国产免费男女视频| 亚洲精品美女久久久久99蜜臀| 国产成人aa在线观看| 男女那种视频在线观看| 伊人久久大香线蕉亚洲五| 色播亚洲综合网| 啦啦啦韩国在线观看视频| 美女大奶头视频| 精品欧美国产一区二区三| 日本熟妇午夜| 久久久精品欧美日韩精品| 国产人伦9x9x在线观看| 五月伊人婷婷丁香| 国产成人一区二区三区免费视频网站| 亚洲中文字幕日韩| 久久亚洲真实| 久久精品亚洲精品国产色婷小说| 亚洲18禁久久av| 女人被狂操c到高潮| 午夜福利在线观看免费完整高清在 | 97碰自拍视频| 最近在线观看免费完整版| 不卡av一区二区三区| 变态另类丝袜制服| 国产精品久久电影中文字幕| 久久精品人妻少妇| 国产精华一区二区三区| 成人国产综合亚洲| 99国产综合亚洲精品| 国产成人av激情在线播放| 国产精品久久久av美女十八| 国产成人av教育| 国产精品99久久99久久久不卡| 欧美性猛交黑人性爽| 欧美国产日韩亚洲一区| 国产真人三级小视频在线观看| 亚洲成av人片免费观看| 婷婷精品国产亚洲av| 亚洲性夜色夜夜综合| 91老司机精品| 国产私拍福利视频在线观看| 精品日产1卡2卡| 欧美+亚洲+日韩+国产| 99riav亚洲国产免费| 成人永久免费在线观看视频| 露出奶头的视频| 制服人妻中文乱码| 国产精品av视频在线免费观看| 久久久久性生活片| 亚洲 国产 在线| 午夜免费激情av| 亚洲欧美日韩东京热| 19禁男女啪啪无遮挡网站| 国产久久久一区二区三区| 婷婷精品国产亚洲av在线| 性色avwww在线观看| 亚洲一区二区三区不卡视频| 亚洲av电影在线进入| 男女那种视频在线观看| 国产欧美日韩一区二区三| 狂野欧美激情性xxxx| 日韩人妻高清精品专区| av女优亚洲男人天堂 | 国产美女午夜福利| 1000部很黄的大片| 香蕉国产在线看| 国产精品99久久99久久久不卡| 久久中文字幕人妻熟女| 午夜日韩欧美国产| 亚洲av日韩精品久久久久久密| 日本一本二区三区精品| 国产精品久久久人人做人人爽| 国产成人欧美在线观看| 久久人人精品亚洲av| 精品欧美国产一区二区三| 国产精品1区2区在线观看.| 噜噜噜噜噜久久久久久91| 一本精品99久久精品77| 亚洲美女视频黄频| 精品久久久久久成人av| 日韩成人在线观看一区二区三区| 亚洲国产色片| 人妻夜夜爽99麻豆av| 欧美xxxx黑人xx丫x性爽| 亚洲专区字幕在线| 欧美性猛交╳xxx乱大交人| 男女那种视频在线观看| 国产欧美日韩一区二区三| 精品午夜福利视频在线观看一区| 美女大奶头视频| 久久这里只有精品19| 国产极品精品免费视频能看的| 亚洲五月天丁香| 在线十欧美十亚洲十日本专区| 亚洲av日韩精品久久久久久密| 欧美一级a爱片免费观看看| av欧美777| 欧美中文日本在线观看视频| 国产男靠女视频免费网站| 好男人在线观看高清免费视频| 亚洲欧美精品综合一区二区三区| 久久中文字幕一级| 黑人操中国人逼视频| 免费高清视频大片| 国产高清三级在线| 国产91精品成人一区二区三区| 亚洲男人的天堂狠狠| 村上凉子中文字幕在线| xxxwww97欧美| 亚洲乱码一区二区免费版| 操出白浆在线播放| 久久精品人妻少妇| 午夜激情欧美在线| 国产乱人视频| 国产成人一区二区三区免费视频网站| 成人国产一区最新在线观看| 午夜激情福利司机影院| 99国产极品粉嫩在线观看| 精品熟女少妇八av免费久了| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美日韩高清专用| 国产高潮美女av| 久久午夜亚洲精品久久| 黄色 视频免费看| 国产亚洲精品久久久com| 国产精品一区二区精品视频观看| 精品久久久久久久末码| 国产精品一区二区精品视频观看| 国产精品久久视频播放| 一个人看的www免费观看视频| 性色av乱码一区二区三区2| 搡老熟女国产l中国老女人| 女生性感内裤真人,穿戴方法视频| 天堂网av新在线| 国产精品一及| 国产精品一区二区免费欧美| 波多野结衣高清作品| 国产人伦9x9x在线观看| 精品乱码久久久久久99久播| 国产成人av激情在线播放| 欧美黄色片欧美黄色片| 国内精品久久久久精免费| 岛国视频午夜一区免费看| 国产又黄又爽又无遮挡在线| 99精品欧美一区二区三区四区| 一本精品99久久精品77| 美女高潮喷水抽搐中文字幕| 免费观看精品视频网站| 变态另类丝袜制服| 亚洲中文日韩欧美视频| 俄罗斯特黄特色一大片| 亚洲自偷自拍图片 自拍| 久久久精品欧美日韩精品| 中文字幕av在线有码专区| 色综合站精品国产| 一级a爱片免费观看的视频| 男插女下体视频免费在线播放| 国产乱人伦免费视频| 亚洲七黄色美女视频| 国产成+人综合+亚洲专区| 黄色视频,在线免费观看| 久久精品国产综合久久久| 嫩草影院精品99| 成人精品一区二区免费| 听说在线观看完整版免费高清| 9191精品国产免费久久| 亚洲无线在线观看| 亚洲一区二区三区不卡视频| 日本a在线网址| 午夜日韩欧美国产| 亚洲av美国av| 日本黄色片子视频| 精品国产亚洲在线| 国产亚洲精品久久久com| 动漫黄色视频在线观看| 一本一本综合久久| 成人国产一区最新在线观看| 夜夜夜夜夜久久久久| 亚洲自拍偷在线| 精品久久蜜臀av无| 久久精品aⅴ一区二区三区四区| 欧洲精品卡2卡3卡4卡5卡区| 美女黄网站色视频| 日日夜夜操网爽| 19禁男女啪啪无遮挡网站| 免费一级毛片在线播放高清视频| 欧美日韩中文字幕国产精品一区二区三区| 在线观看舔阴道视频| 偷拍熟女少妇极品色| 国产亚洲精品一区二区www| 极品教师在线免费播放| 免费看光身美女| 91久久精品国产一区二区成人 | 每晚都被弄得嗷嗷叫到高潮| 美女 人体艺术 gogo| 女人高潮潮喷娇喘18禁视频| 成人三级做爰电影| 久久精品国产清高在天天线| 精品一区二区三区av网在线观看| 国产成人福利小说| 啦啦啦韩国在线观看视频| 久久久久久国产a免费观看| 免费在线观看日本一区| 中文字幕人妻丝袜一区二区| 小蜜桃在线观看免费完整版高清| 免费观看人在逋| 免费看光身美女| 免费人成视频x8x8入口观看| 日本成人三级电影网站| 国产淫片久久久久久久久 | 麻豆一二三区av精品| 舔av片在线| 亚洲av成人一区二区三| 欧美日韩乱码在线| 成人欧美大片| 久久中文字幕人妻熟女| 亚洲av五月六月丁香网| 十八禁人妻一区二区| 久久精品91蜜桃| 性色avwww在线观看| 色哟哟哟哟哟哟| 狂野欧美白嫩少妇大欣赏| 国产男靠女视频免费网站| 99久国产av精品| 久久婷婷人人爽人人干人人爱| 久久人妻av系列| 99在线人妻在线中文字幕| 色综合婷婷激情| 哪里可以看免费的av片| 国产亚洲精品综合一区在线观看| 国内毛片毛片毛片毛片毛片| 欧美不卡视频在线免费观看| 毛片女人毛片| av天堂在线播放| 精品午夜福利视频在线观看一区| 宅男免费午夜| 亚洲 欧美一区二区三区| 成在线人永久免费视频| 国产乱人伦免费视频| 禁无遮挡网站| 午夜福利在线在线| 亚洲欧美精品综合久久99| 他把我摸到了高潮在线观看| 怎么达到女性高潮| www.999成人在线观看| 精品国内亚洲2022精品成人| 国产一区二区三区在线臀色熟女| 国产精品一区二区三区四区免费观看 | 岛国在线免费视频观看| 久久久久久久久中文| 午夜免费观看网址| 国产激情偷乱视频一区二区| 舔av片在线| 国产精品国产高清国产av| 中国美女看黄片| 国内精品美女久久久久久| 免费在线观看成人毛片| 国产主播在线观看一区二区| 国产成人av激情在线播放| 久久久久久久午夜电影| 一个人看视频在线观看www免费 | 观看免费一级毛片| 亚洲中文av在线| 日韩欧美在线乱码| 免费搜索国产男女视频| 麻豆久久精品国产亚洲av| 别揉我奶头~嗯~啊~动态视频| 夜夜看夜夜爽夜夜摸| 99精品久久久久人妻精品| 国产乱人视频| 老司机福利观看| 草草在线视频免费看| 人妻丰满熟妇av一区二区三区| 综合色av麻豆| h日本视频在线播放| 亚洲成av人片免费观看| 2021天堂中文幕一二区在线观| 亚洲成a人片在线一区二区| 国内揄拍国产精品人妻在线| 久久久精品大字幕| 国产欧美日韩精品亚洲av| 亚洲无线观看免费| 亚洲国产精品成人综合色| 成人特级黄色片久久久久久久| 午夜精品久久久久久毛片777| 国产av在哪里看| 亚洲欧美精品综合久久99| 国产高清激情床上av| 成熟少妇高潮喷水视频| 香蕉av资源在线| 波多野结衣巨乳人妻| 免费看美女性在线毛片视频| 精品福利观看| 欧美黑人巨大hd| 国内精品美女久久久久久| 国产日本99.免费观看| 男人舔女人下体高潮全视频| 亚洲欧美日韩高清在线视频| 国产一区二区三区在线臀色熟女| 亚洲欧美精品综合一区二区三区| 一进一出抽搐动态| 中文字幕人成人乱码亚洲影| 精品一区二区三区av网在线观看| 久久精品综合一区二区三区| 超碰成人久久| 亚洲色图 男人天堂 中文字幕| 日韩精品青青久久久久久| 婷婷丁香在线五月| 亚洲18禁久久av| 又黄又粗又硬又大视频| 久久九九热精品免费| 婷婷精品国产亚洲av在线| 色尼玛亚洲综合影院| 成年女人永久免费观看视频| 免费在线观看视频国产中文字幕亚洲| 久久亚洲真实| 变态另类成人亚洲欧美熟女| 十八禁人妻一区二区| 日韩中文字幕欧美一区二区| 中文字幕精品亚洲无线码一区| 热99在线观看视频| 国产欧美日韩精品一区二区| 69av精品久久久久久| 一区二区三区高清视频在线| 欧美中文日本在线观看视频| 五月玫瑰六月丁香| 搞女人的毛片| 毛片女人毛片| 亚洲成人中文字幕在线播放| 免费观看精品视频网站| 亚洲七黄色美女视频| 1000部很黄的大片| 欧美极品一区二区三区四区| 日本精品一区二区三区蜜桃| 天天一区二区日本电影三级| 国产淫片久久久久久久久 | 韩国av一区二区三区四区| 精品久久久久久久末码| 嫩草影视91久久| 高清在线国产一区| 国产精品精品国产色婷婷| 一级黄色大片毛片| av片东京热男人的天堂| 99热只有精品国产| 久久精品亚洲精品国产色婷小说| 久久午夜亚洲精品久久| 啦啦啦免费观看视频1| 亚洲欧美日韩高清在线视频| 国产成人精品久久二区二区91| 黄色 视频免费看| 91麻豆av在线| 天堂√8在线中文| 成人高潮视频无遮挡免费网站| 久久午夜亚洲精品久久| 亚洲国产看品久久| 成人精品一区二区免费| 久久天躁狠狠躁夜夜2o2o| 99国产极品粉嫩在线观看| 最近最新中文字幕大全电影3| 丰满人妻一区二区三区视频av | 亚洲中文字幕日韩| 日本黄色片子视频| 亚洲18禁久久av| 久久久久国产精品人妻aⅴ院| 99在线视频只有这里精品首页| 少妇丰满av| 亚洲中文日韩欧美视频| 成人无遮挡网站| 久久久久亚洲av毛片大全| 最近最新中文字幕大全电影3| 国内毛片毛片毛片毛片毛片| 色哟哟哟哟哟哟| 亚洲男人的天堂狠狠| 久久久久久久午夜电影| 在线观看免费视频日本深夜| 久久久国产成人精品二区| 久久精品夜夜夜夜夜久久蜜豆| 色综合婷婷激情| 国产精华一区二区三区| 丁香欧美五月| 亚洲av中文字字幕乱码综合| 中亚洲国语对白在线视频| 99国产极品粉嫩在线观看| 亚洲人成网站高清观看| 国产成人精品无人区| 精品电影一区二区在线| 午夜久久久久精精品| 一个人免费在线观看电影 | 精品日产1卡2卡| 国产精品九九99| 黄色日韩在线| 成人鲁丝片一二三区免费| 九色成人免费人妻av| 麻豆成人av在线观看| 美女高潮的动态| 日韩三级视频一区二区三区| 亚洲最大成人中文| 日韩有码中文字幕| 国产免费av片在线观看野外av| 成人永久免费在线观看视频| 国产欧美日韩一区二区精品| 99久久精品热视频| 国产免费男女视频| 午夜精品久久久久久毛片777| 日韩欧美一区二区三区在线观看| 在线观看舔阴道视频| 香蕉国产在线看| 久久精品国产亚洲av香蕉五月| 国产亚洲精品av在线| 久久久色成人|