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

    持續(xù)蒸發(fā)與補(bǔ)給蒸發(fā)過程中水體咸化及同位素分餾的實(shí)驗(yàn)研究

    2015-07-02 03:19:58靜,聰,杏,斌,
    地球化學(xué) 2015年6期
    關(guān)鍵詞:氫氧水化學(xué)衡水

    李 靜, 王 聰, 梁 杏, 馬 斌, 葛 勤

    (1. 中國地質(zhì)大學(xué)(武漢) 環(huán)境學(xué)院, 湖北 武漢 430074; 2. 國土資源部 地裂縫地質(zhì)災(zāi)害重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210018; 3. 湖北省 環(huán)境科學(xué)研究院, 湖北 武漢 430072; 4. 生物地質(zhì)與環(huán)境地質(zhì)國家重點(diǎn)實(shí)驗(yàn)室, 中國地質(zhì)大學(xué)(武漢), 湖北 武漢 430074)

    0 引 言

    華北平原淺層廣泛分布的咸水體[1–2], 是在 12 ka前后干旱的氣候環(huán)境下, 受蒸發(fā)濃縮作用形成, 即大陸鹽化咸水層[3]。查明蒸發(fā)過程中水體的咸化及氫氧同位素分餾機(jī)制, 對于揭示地表水或地下水的成因與演化規(guī)律具有重要意義[4–5]。很多學(xué)者已開展蒸發(fā)實(shí)驗(yàn)研究了水體在咸化過程中水化學(xué)及同位素的演變機(jī)制, 得出了蒸發(fā)過程中的蒸發(fā)結(jié)晶路線、析鹽規(guī)律[6–8]、氫氧同位素的分餾效應(yīng)及濕度和溫度對其的影響[9–11]。自然界中水體處于不斷的水文循環(huán)當(dāng)中, 通常大陸鹽化作用涉及水的(遷移)流動和蒸發(fā)(排泄)兩個(gè)方面。而前人[12–15]對蒸發(fā)過程的實(shí)驗(yàn)研究大多采用持續(xù)蒸發(fā)的實(shí)驗(yàn)方式, 未考慮蒸發(fā)過程中水流補(bǔ)給的影響。在實(shí)際水文循環(huán)系統(tǒng)中,蒸發(fā)排泄和水的側(cè)向徑流補(bǔ)給(或大氣降水補(bǔ)給)同時(shí)存在, 這種條件下的水化學(xué)和同位素演變鮮見報(bào)道,有待進(jìn)一步研究。本文擬利用唐山和衡水地下水樣品,通過設(shè)計(jì)持續(xù)蒸發(fā)和補(bǔ)水蒸發(fā)兩種模式開展蒸發(fā)實(shí)驗(yàn)研究, 旨在查明在有水源補(bǔ)給條件下蒸發(fā)過程中水化學(xué)和同位素富集規(guī)律及其與純蒸發(fā)的差異。

    1 實(shí)驗(yàn)方案及測試

    1.1 持續(xù)蒸發(fā)實(shí)驗(yàn)

    分別采用唐山地下淡水(TDS = 0.7 g/L)和衡水地下微咸水(TDS = 1.9 g/L)進(jìn)行實(shí)驗(yàn), 蒸發(fā)溫度均為60 ℃。實(shí)驗(yàn)設(shè)計(jì): 唐山水樣取450 mL盛放在燒杯中, 利用電熱板加熱。實(shí)驗(yàn)過程分4個(gè)蒸發(fā)階段,利用量筒測量蒸發(fā)前后水體積變化, 即蒸發(fā)量。每階段蒸發(fā)量控制為50 mL, 采取表層水體50 mL進(jìn)行測試, 總蒸發(fā)時(shí)長約216 h。衡水水樣采用了6個(gè)盛放250 mL水樣的燒杯和2個(gè)盛放500 mL水樣的燒杯進(jìn)行加熱, 不同燒杯加熱不同時(shí)長, 即分 8個(gè)蒸發(fā)階段, 測量每個(gè)燒杯的蒸發(fā)量, 采集表層水體進(jìn)行測試, 最大蒸發(fā)時(shí)長約134 h。所取唐山水樣測試水化學(xué)組分, 衡水水樣測試水化學(xué)組分和氫氧同位素。

    1.2 補(bǔ)給蒸發(fā)實(shí)驗(yàn)

    分別把盛放唐山水樣 500 mL的燒杯和衡水水樣300 mL的燒杯放在電熱板加熱。分7個(gè)蒸發(fā)階段,測量每一階段的蒸發(fā)量, 設(shè)計(jì)補(bǔ)水量等于蒸發(fā)量(水體總量保持不變, 除了取樣損失)。每次取樣前, 先使用原始水樣補(bǔ)充蒸發(fā)的水量, 使其混合均勻。唐山水樣蒸發(fā)率控制為20%, 蒸發(fā)總時(shí)長469 h, 每次取樣 50 mL, 測試水化學(xué)組成。衡水水樣蒸發(fā)率約50%, 蒸發(fā)總時(shí)長 191 h, 每次取樣 25 mL, 測試水化學(xué)組成和氫氧同位素。

    1.3 測試項(xiàng)目與方法

    所采水樣水化學(xué)測試包括pH、電導(dǎo)率、陰陽離子。其中pH和電導(dǎo)率采用哈希便攜式水質(zhì)分析參數(shù)儀(儀器型號: HQ30d)測試; 堿度(HCO3–計(jì))采用滴定法測試; 陰陽離子分別采用離子色譜儀(儀器型號: DIONEX-120)和 ICP-AES(儀器型號:IRIS Intrepid Ⅱ X SP 型)測試, 精度分別為 0.1和0.01。

    氫氧同位素的分析儀器為 Mat253氣體同位素質(zhì)譜儀, 氫同位素測試精度優(yōu)于2‰, 氧同位素測試精度優(yōu)于0.2‰。

    2 結(jié)果與討論

    實(shí)驗(yàn)水樣的原始組成及蒸發(fā)水樣的測試結(jié)果如表1和表2所示。其中唐山水樣主要用來分析兩種蒸發(fā)實(shí)驗(yàn)中水化學(xué)演變差異, 衡水水樣主要用來分析同位素富集的差異。

    2.1 化學(xué)組分演化差異

    在持續(xù)蒸發(fā)與補(bǔ)給蒸發(fā)實(shí)驗(yàn)中, 水樣的電導(dǎo)率(EC)和pH隨著蒸發(fā)時(shí)間總體上呈增大的趨勢(圖1)。EC的升高表明鹽分的濃縮, pH的升高是由于水中CO2的釋出造成的。但兩種蒸發(fā)過程中, EC和pH變化曲線表現(xiàn)出明顯分離: 持續(xù)蒸發(fā)條件下, 蒸發(fā)約200 h, 水樣電導(dǎo)率由原始的1107 μS/cm增加到1881 μS/cm, pH由7.8增加到9.0, 濃縮速率逐漸增大。補(bǔ)給蒸發(fā)條件下, 鹽分的濃縮速率明顯較低, 蒸發(fā)約 500 h, 電導(dǎo)率由原始的 1107 μS/cm 增加到1865 μS/cm, pH由7.8增加到8.9。這是因?yàn)樵谟械a(bǔ)給過程中, 蒸發(fā)和稀釋作用同時(shí)存在, 抵消了部分蒸發(fā)效果。因此在干旱氣候環(huán)境中, 水體濃度的高低不僅取決于蒸發(fā)強(qiáng)度, 還與水的徑流(循環(huán))密切相關(guān)。

    表1 唐山原始水樣(TC0)及蒸發(fā)樣品水化學(xué)測試結(jié)果(mg/L)Table 1 Chemistry data (mg/L) for original water (TC0) and evaporation water samples of Tangshan

    表2 衡水原始水樣(HB0)和蒸發(fā)樣品同位素測試結(jié)果(‰)Table 2 Isotopic data (‰) for original water (HB0) and evaporation water samples of Hengshui

    圖1 唐山淡水持續(xù)蒸發(fā)和補(bǔ)給蒸發(fā)模式下電導(dǎo)率和pH變化規(guī)律Fig.1 Diagrams showing electrical conductivity and pH variation trends with time under continuous evaporation and recharge evaporation for Tangshan freshwater samples

    仔細(xì)觀察圖 1可以發(fā)現(xiàn), 持續(xù)蒸發(fā)過程中水樣的電導(dǎo)率在蒸發(fā)初期(0~100 h)存在下降過程, 蒸發(fā)實(shí)驗(yàn)中濃度一般不會出現(xiàn)降低現(xiàn)象。本情況的發(fā)生可能與唐山水樣有關(guān), 實(shí)驗(yàn)水樣的水化學(xué)類型為HCO3·Cl·SO4-Ca·Mg·Na 型水。以 H CO3-、Ca2+為主的水樣, 在蒸發(fā)過程中水溫增加, 極易發(fā)生如下化學(xué)反應(yīng)(1)而生成 CaCO3, 溶解度大大減小, 發(fā)生沉淀。在蒸發(fā)初期當(dāng)該反應(yīng)產(chǎn)生的沉淀量大于其他離子濃縮程度時(shí), 造成了電導(dǎo)率的下降。

    各離子的變化過程可進(jìn)一步說明這一點(diǎn)。如圖2所示, 持續(xù)蒸發(fā)過程中Ca2+和 H CO3-含量在蒸發(fā)初期濃度下降, 一旦 CaCO3沉淀存在于水樣中, 此水環(huán)境中則 Ca2+達(dá)到飽和, 如圖 2所示后期 Ca2+保持穩(wěn)定。H CO3-離子在蒸發(fā)過程中可繼續(xù)分解為 C O32-,可與Mg2+結(jié)合生成MgCO3微溶沉淀, 但從圖2中可以看出Mg2+處于持續(xù)緩慢增加的過程中, 而 H CO3-含量并未持續(xù)下降, 而是相對穩(wěn)定, 表明在蒸發(fā)過程中基本不存在MgCO3的析出或含量甚少。在補(bǔ)水蒸發(fā)過程中, 雖然存在 Ca2+和 H CO3-下降過程, 但降低值較小, 主要與受到淡水補(bǔ)給有關(guān), 電導(dǎo)率并未顯示出下降值。

    其他各離子(S O24-、Mg2+、Na+、Cl–和 K+)與電導(dǎo)率同步變化, 兩種蒸發(fā)過程離子與電導(dǎo)率關(guān)系曲線基本吻合(K+含量較低, 在誤差范圍內(nèi)波動)。說明在實(shí)驗(yàn)設(shè)計(jì)蒸發(fā)程度下, 水中礦物在蒸發(fā)早期以方解石(CaCO3)沉淀為主, 其他鹽分以濃縮作用為主。李亞文等[8]所做的南海蒸發(fā)實(shí)驗(yàn)得出在海水蒸發(fā)了80%后, 才開始析出石膏。在整個(gè)蒸發(fā)過程中, 方解石的動態(tài)溶解-沉淀平衡和鹽分的濃縮作用共存。持續(xù)蒸發(fā)過程中, 早期以方解石析出作用為主, 引起電導(dǎo)率下降, 后期以鹽分的濃縮為主; 在補(bǔ)水蒸發(fā)過程中, 方解石早期的析出作用較小, 如 Ca2+和HCO3-含量高于持續(xù)蒸發(fā)實(shí)驗(yàn), 以鹽分的濃縮為主。

    綜上可知, 持續(xù)蒸發(fā)與補(bǔ)水蒸發(fā)過程中水化學(xué)變化的差異主要為濃縮速率的不同。同時(shí) Ca2+和HCO3-濃度在兩種蒸發(fā)過程中差異明顯, 是由于這兩種離子參與了化學(xué)反應(yīng); 當(dāng)溶液中離子的礦物未飽和時(shí), 如 Na+、K+、 S O24-、Cl–以蒸發(fā)濃縮作用為主。其離子變化主要受水樣鹽度控制, 與是否有水源補(bǔ)給影響不大。兩種蒸發(fā)方式下, 水化學(xué)類型也發(fā)生規(guī)律性的變化。持續(xù)蒸發(fā): 水型變化,HCO3·Cl·SO4-Ca·Mg·Na → HCO3·Cl·SO4-Mg ·Ca·Na →HCO3·Cl·SO4-Mg·Na → Cl·SO4·HCO3-Mg·Na → Cl·SO4-Mg·Na; 補(bǔ) 給 蒸 發(fā) : HCO3·Cl·SO4-Ca·Mg·Na →Cl·HCO3·SO4-Mg·Na·Ca → Cl·HCO3·SO4-Mg·Na →Cl·SO4·HCO3-Mg·Na → Cl·SO4-Mg·Na。水型由復(fù)雜向簡單轉(zhuǎn)化, H CO3-和Ca2+由主要離子逐漸退出水型組成。

    2.2 氫氧同位素的富集過程

    2.2.1 氫氧同位素與剩余水關(guān)系

    持續(xù)蒸發(fā)實(shí)驗(yàn)與前人實(shí)驗(yàn)方法相似[9], 假定溫度和濕度都不變的情況下, 每次揮發(fā)出的水汽即刻分離出系統(tǒng), 剩余水的同位素逐漸富集, 此過程符合瑞利平衡蒸發(fā)模型。瑞利蒸發(fā)模型是英國科學(xué)家瑞利(Rayleigh)討論液體混合物蒸發(fā)過程時(shí)提出來的一種模型[12]。在瑞利過程中, 體系中的兩相物質(zhì)處于一種瞬時(shí)平衡狀態(tài), 其過程可用下式來描述:

    式中: Rl為任一瞬間反應(yīng)物(剩余水體)的同位素比值;R0為蒸發(fā)水體初始同位素比值; f為剩余水體體積比;α 為液-氣間分餾系數(shù)。利用公式 δ = (R樣品/R標(biāo)準(zhǔn)– 1)轉(zhuǎn)化可得δ形式下的瑞利分餾模型:

    可以看出在瑞利分餾模式下, 剩余水體同位素δ與剩余水體積比f呈指數(shù)關(guān)系[13]。

    圖2 唐山淡水不同蒸發(fā)過程中水樣各離子濃度與電導(dǎo)率關(guān)系Fig.2 Diagrams showing relationships between ion content and electrical conductivity under the two evaporation experiments for Tangshan freshwater samples TC–持續(xù)蒸發(fā); TB–補(bǔ)給蒸發(fā)。TC–continuous evaporation; TB–recharge evaporation.

    持續(xù)蒸發(fā)過程中, 剩余水氫氧同位素含量隨著剩余水體積比的減少而呈指數(shù)富集, 這與原有的瑞利平衡分餾模型結(jié)果一致(圖3實(shí)線)。實(shí)驗(yàn)設(shè)計(jì)蒸發(fā)溫度為60 ℃ , 平衡蒸發(fā)模型的分餾系數(shù)α只與溫度相關(guān), 據(jù)文獻(xiàn)[16]可知: α18O= 1.0061, αD= 1.046,由公式(2)可得平衡分餾模型的剩余水同位素富集曲線(圖3中虛線)。由圖3可知, 實(shí)驗(yàn)測試剩余水同位素值高于平衡分餾模型趨勢線, 說明實(shí)驗(yàn)在蒸發(fā)過程中同位素分餾不僅與溫度有關(guān), 還受到濕度的影響, 如存在水-汽間同位素交換等非平衡條件下的動力分餾。其次, 在蒸發(fā)過程中, 水體的氧同位素值還會受到方解石析出的影響。一般巖石及礦物的 δ18O值要比水的δ18O值大得多[16], 如發(fā)生礦物與水的同位素交換反應(yīng), 會明顯增大水溶液的δ18O值。周根陶等[17]指出碳酸鈣緩慢沉淀過程中, 碳酸鈣與溶液中之間有足夠的時(shí)間進(jìn)行同位素交換, 而碳酸鈣快速沉淀, 則繼承溶液中碳酸根離子的同位素組成。本實(shí)驗(yàn)中碳酸鈣持續(xù)沉淀, 可能與水體發(fā)生同位素交換反應(yīng)。但受溫度影響(< 60°)水與巖石礦物發(fā)生氫氧同位素交換的反應(yīng)速度十分緩慢, 一般很難達(dá)到交換平衡[16]。

    圖3 衡水微咸水持續(xù)蒸發(fā)實(shí)驗(yàn)中剩余水氫氧同位素與剩余水體積比關(guān)系Fig.3 Diagrams showing relationships between stable isotopic value in residual water and residual water ratio under continuous evaporation experiment for Hengshui brackish water samples實(shí)線–實(shí)測數(shù)據(jù)趨勢線; 虛線–瑞利平衡模型趨勢線。Solid line – trend line of measured data; Dotted line – trend line of Rayleigh balance mode.

    從數(shù)值來看, 隨著蒸發(fā)率的增大, δ18O 由–9.19‰增加到 20.81‰, 變化值為 30‰; δD 由–62.6‰增加到56.9‰, 變化值為119.5‰, 同等條件下, D的富集速率約是18O的4倍, 說明氫同位素在富集過程中更為敏感, 這是由于H218O和DHO蒸汽壓的差別使他們在蒸發(fā)過程中在水相的富集不成比例[18]。在 δD-δ18O 關(guān)系圖中, 蒸發(fā)線斜率的大小主要取決于相對濕度, 其大小通常介于 4~6之間[19],如圖 5所示, 蒸發(fā)水樣點(diǎn)均落在全球雨水線的右下方, 并呈線性關(guān)系 δD = 4.11δ18O – 25.23 (R2= 0.99),斜率為 4.11, 即為持續(xù)蒸發(fā)實(shí)驗(yàn)條件下的蒸發(fā)線方程(圖 5)。

    2.2.2 氫氧同位素富集與鹽分的關(guān)系

    在持續(xù)蒸發(fā)模式下, 原始水樣由初始 TDS為1.9 g/L增加到11.7 g/L。在含鹽的溶液中, 氫氧同位素的分餾系數(shù)不僅是溫度的函數(shù), 而且也是鹽度的函數(shù), 在離子水合作用影響下, 含鹽溶液中的重同位素水分子的活度系數(shù)與輕同位素水分子的活度系數(shù)之比不等于 1, 從而影響了鹽水蒸發(fā)過程中氫氧同位素的分餾系數(shù), 使鹽水的蒸發(fā)不同于淡水的蒸發(fā)[19]。

    如圖4所示, 在蒸發(fā)過程中, δ18O值隨著陰陽離子含量的升高呈現(xiàn)兩段式富集: 在低鹽水中, 氫氧同位素隨著離子含量升高, 富集速率較快, 當(dāng)鹽度增加到一定程度時(shí), 氫氧同位素富集速率變慢。以Na+為例, 在Na+含量小于500 mg/L時(shí), δ18O與Na+關(guān)系方程為 δ18O = 0.09 × Na – 31.9 (R2= 0.99), 在Na+含量大于 500 mg/L時(shí), δ18O與 Na+關(guān)系方程為δ18O = 0.01 × Na – 4.08 (R2= 0.96)。Mg2+、Cl–、S O24-與 Na+相似, 在蒸發(fā)過程中均與 δ18O成兩段式線性關(guān)系。

    隨著鹽度增加, 水體相對純水18O富集的現(xiàn)象,即為鹽效應(yīng)。Horita[19]研究得出鹽水(礦化度大于海水, 50~150 g/L)中陰離子(Cl–, S O24-)對 D和 δ18O 的分餾具有重要影響, 隨著溫度增加, 鹽效應(yīng)降低。如蒸發(fā)溫度從10 ℃增加到150 ℃, D同位素的鹽效應(yīng)由2.4‰/mol降低到1‰/mol。其次, 堿土金屬(Ca2+、Mg2+)對D和18O產(chǎn)生的鹽效應(yīng)具有方向上的差別[20]。而本實(shí)驗(yàn)中未設(shè)置相應(yīng)條件下淡水蒸發(fā)實(shí)驗(yàn), 無法比對鹽度對氫氧同位素分餾的影響, 但從實(shí)驗(yàn)結(jié)果看陰陽離子與δ18O的關(guān)系一致, 可能是由于該水樣濃度較低(< 2 g/L), 對氫氧同位素的分餾未產(chǎn)生顯著影響。

    圖4 衡水微咸水陰陽離子與δ18O關(guān)系Fig.4 Diagrams showing relationships of ions content vs. δ18O for Hengshui brackish water samples

    圖5 衡水微咸水不同蒸發(fā)實(shí)驗(yàn)?zāi)J较滤碾妼?dǎo)率與氫氧同位素關(guān)系(a)及水的蒸發(fā)線方程(b)Fig.5 Diagrams showing relationships between stable isotopic value for residual water and electrical conductivity (a) and evaporation lines (b)under the two evaporation experiments of Hengshui brackish water samples

    2.2.3 持續(xù)蒸發(fā)與補(bǔ)給蒸發(fā)過程中同位素富集差異

    如圖 5a所示, 在持續(xù)蒸發(fā)與補(bǔ)給蒸發(fā)過程中,隨著蒸發(fā)的進(jìn)行, 水的鹽度(電導(dǎo)率)逐漸增大, δD的變化范圍明顯大于δ18O的變化范圍, 進(jìn)一步說明D同位素在蒸發(fā)過程中富集速率大于氧同位素。對比兩種蒸發(fā)實(shí)驗(yàn)可以發(fā)現(xiàn), 達(dá)到相同的鹽度(電導(dǎo)率), 持續(xù)蒸發(fā)模式下氫氧同位素值相對補(bǔ)給蒸發(fā)模式明顯偏正。這是由于持續(xù)蒸發(fā)情況下, 輕同位素逐漸揮發(fā)出系統(tǒng), 重同位素持續(xù)富集; 而補(bǔ)給蒸發(fā)模式下, 由于存在原水的補(bǔ)給, 同位素不斷受到混合, 減弱了重同位素的富集程度。

    其次, 持續(xù)蒸發(fā)模式下, 同位素隨著鹽度(電導(dǎo)率)的增大持續(xù)增加, δ18O最大值達(dá)20.81‰, 而補(bǔ)給蒸發(fā)模式下, 水樣的電導(dǎo)率在后期緩慢增加, 氧同位素趨于穩(wěn)定, δ18O最大值達(dá)2.9‰。Gonfiantini[20]指出受限盆地在過度蒸發(fā)的條件下δ18O將會增大趨近于一個(gè)穩(wěn)定值, 該值由補(bǔ)給水的流入和相對濕度控制。補(bǔ)給蒸發(fā)模式控制蒸發(fā)率為 50%, 每蒸發(fā)近一半水體積, 則進(jìn)行補(bǔ)給, 在這種蒸發(fā)條件下, 氧同位素的富集由補(bǔ)給水源和氣象因素(溫度、濕度)共同控制。

    由氫氧同位素關(guān)系曲線可知(圖 5b), 不同蒸發(fā)過程, 蒸發(fā)線方程也不同。持續(xù)蒸發(fā)模式蒸發(fā)線斜率為4.11大于補(bǔ)給蒸發(fā)模式蒸發(fā)線斜率3.41。王永森等[13]研究得出自由水體蒸發(fā)線斜率主要受溫度和濕度控制。在溫度不變的情況下, 蒸發(fā)線線斜率隨濕度的增加而增加; 在濕度不變的情況下, 斜率隨溫度的升高而減低。補(bǔ)給蒸發(fā)模式蒸發(fā)線斜率較低主要受同位素混合作用和相對更強(qiáng)的非平衡動力分餾影響。在有補(bǔ)給源存在的情況下, 打破原始的水環(huán)境狀態(tài), 如水-汽的平衡及水溫的變化, 水分子擴(kuò)散運(yùn)動相對強(qiáng)烈等。

    綜上可知, 補(bǔ)給蒸發(fā)方式下水同位素的富集過程更接近實(shí)際情況, 在有其他水補(bǔ)給情況下, 氫氧同位素在蒸發(fā)過程中不會持續(xù)富集, 且蒸發(fā)線斜率小于持續(xù)蒸發(fā)實(shí)驗(yàn)值。

    3 結(jié) 論

    (1) 在持續(xù)蒸發(fā)和補(bǔ)給蒸發(fā)兩種模式下, 水化學(xué)演化的相同點(diǎn)表現(xiàn)為: 逐漸咸化和弱堿性化。水中Ca2+、 H CO3-濃度在初期存在瞬時(shí)降低并保持穩(wěn)定, Mg+、Na+、Cl–、 S O24-濃度逐漸升高, 該水樣水化學(xué)類型由復(fù)雜向簡單型轉(zhuǎn)變; 不同點(diǎn)表現(xiàn)為持續(xù)蒸發(fā)過程中水體濃縮速率較快, 蒸發(fā)初期水體以方解石沉淀作用為主, 引起電導(dǎo)率下降, 后期方解石處于動態(tài)平衡中, 其他組分以濃縮為主; 補(bǔ)給蒸發(fā)過程中水體濃縮速率明顯減慢, 水中Ca2+、H CO3-含量較高, 鹽分濃縮掩蓋了方解石早期析出, 電導(dǎo)率持續(xù)升高。

    (2) 蒸發(fā)過程中, 氫氧同位素逐漸偏正, 氫同位素相對氧同位素更為敏感。氫氧同位素富集與剩余水比例呈指數(shù)關(guān)系, 符合瑞利分餾模型, 同時(shí)也受到界面處水-汽非平衡分餾與水巖同位素交換的影響; 與水中陰陽離子呈兩段式線性關(guān)系, 離子在低濃度段增加時(shí), 氫氧同位素富集較快。

    (3) 得出了兩種蒸發(fā)過程中衡水地區(qū)的蒸發(fā)線方程: 持續(xù)蒸發(fā)條件下, 氫氧同位素的富集隨著電導(dǎo)率的增加持續(xù)偏正, 蒸發(fā)方程為δD = 4.11δ18O –25.23; 補(bǔ)水蒸發(fā)條件下, 氫氧同位素隨著電導(dǎo)率的增加, 富集速率較慢, 在后期相對穩(wěn)定, 同位素分餾程度小于持續(xù)蒸發(fā), 蒸發(fā)線方程為 δD = 3.41δ18O – 28.2。

    :

    [1] 張光輝, 劉春華, 嚴(yán)明疆, 王金哲, 周在明. 環(huán)渤海平原土壤鹽分不同聚形的的水動力學(xué)特征[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2012, 42(6): 1873–1879.Zhang Guang-hui, Liu Chun-hua, Yan Ming-jiang, Wang Jin-zhe, Zhou Zai-ming. Moisture dynamics characteristics of the soils profile of different salinity poly types in the Bohai Sea Plain [J]. J Jilin Univ (Earth Sci Ed), 2012, 42(6):1873–1879 (in Chinese with English abstract).

    [2] 張宗祜, 施德鴻, 任福弘, 殷正宙, 孫繼朝, 張翠云. 論華北平原第四系地下水系統(tǒng)之演化[J]. 中國科學(xué)(D輯), 1997,27(2): 168–173.Zhang Zonghu, Shi Dehong, Ren Fuhong, Yin Zhengzhou,Sun Jichao, Zhang Cuiyun. Evolution of Quaternary groundwater system in North China Plain [J]. Sci China (D), 1997,27(2): 168–173 (in Chinese with English abstract).

    [3] 張人權(quán), 梁杏, 靳孟貴. 末次盛冰期以來河北平原第四系地下水流系統(tǒng)的演變[J]. 地學(xué)前緣, 2013, 20(3): 217–226.Zhang Ren-quan, Liang Xing, Jin Meng-gui. The evolution of groundwater flow systems in the Quaternary of Hebei Plain [J].Earth Sci Front, 2013, 20(3): 217–226 (in Chinese with English abstract).

    [4] Horita J, Wesolowski D J, Cole D R. The activity-composition relationship of oxygen and hydrogen isotopes in aqueous salt solution:Ⅰ. Vapor-liquid water equilibration of single salt solution from 50 to 100 ℃[J]. Geochim Cosmochim Acta,1993, 57 (12): 2797–2817.

    [5] Horita J, Wesolowski D J. Liquid-vapor fractionation of oxygen and hydrogen isotopes of water from the freezing to the critical temperature [J]. Geochim Cosmochim Acta, 1994,58(16): 3425–3437.

    [6] 趙文玉, 王啟山, 吳國平, 吳立波, 孫增峰, 鄭力燕. 引黃水在北大港水庫蓄存期水質(zhì)咸化的研究[J]. 水科學(xué)進(jìn)展,2005, 16(6): 763–766.Zhao Wen-yu, Wang Qi-shan, Wu Guo-ping, Wu Li-bo, Sun Zeng-feng, Zheng Li-yan. Study on salinization of water channeled from Yellow River in Beidagang reservoir during the stored period [J]. Adv Water Sci, 2005, 16(6): 763–766 (in Chinese with English abstract).

    [7] Gremillion P, Wanielista M. Effects of evaporative enrichment on the stable isotope hydrology of a central Florida(USA) river [J]. Hydrol Process, 2000, 14: 1465–1484.

    [8] 李亞文, 韓蔚田. 南海海水 25 ℃等溫蒸發(fā)實(shí)驗(yàn)研究[J]. 地質(zhì)科學(xué), 1995, 30(3): 233–239.Li Ya-wen, Han Wei-tian. Experimental study of 25 ℃ isothermal evaporation of sea water from the South China Sea [J].Sci Geol Sinica, 1995, 30(3): 233–239 (in Chinese with English abstract).

    [9] Hu H Y, Bao W M, Wang T, Qu S M. Experimental study on stable isotopic fractionation of evaporating water under varying temperature [J]. Water Sci Eng, 2009, 2(2): 11–18.

    [10] Zhang X P, Tian L, Liu J M. Fractionation mechanism of stable isotope in evaporating water body [J]. J Geogr Sci, 2005,15(3): 375–384.

    [11] Cappa C D, Hendricks M B, DePaolo D J, Cohen R C. Isotopic fractionation of water during evaporation [J]. J Geophys Res, 2003, 108 (D16): 1–10, doi:10.1029/2003JD003597

    [12] 包為民, 胡海英, 王濤, 瞿思敏. 蒸發(fā)皿中水面蒸發(fā)氫氧同位素分餾的實(shí)驗(yàn)研究[J]. 水科學(xué)進(jìn)展, 2008, 19(6): 780–785.Bao Wei-min, Hu Hai-ying, Wang Tao, Qu Si-min. Experimental study on the fractionation mechanism of hydrogen and oxygen stable isotopes in evaporation from water surface of evaporation pans [J]. Adv Water Sci, 2008, 19( 6): 780–785(in Chinese with English abstract).

    [13] 王永森, 馬振民, 徐征和. 基于瑞利分餾模式的水體蒸發(fā)線斜率模型[J]. 水科學(xué)進(jìn)展, 2011, 22(6): 795–800.Wang Yong-sen, Ma Zhen-min, Xu Zheng-he. Slope of evaporation lines in a model based on Rayleigh fractionation formula [J]. Adv Water Sci, 2011, 22(6): 795–800 (in Chinese with English abstract).

    [14] 王君波, 朱立平, 鞠建廷, 汪勇. 西藏納木錯東部湖水及入湖河流水化學(xué)特征初步研究[J]. 地理科學(xué), 2009, 29(2):288–293.Wang Jun-bo, Zhu Li-ping, Ju Jian-ting, Wang Yong. Water chemistry of Eastern Nam Lake Area and inflowing rivers in Tibet [J]. Sci Geogr Sinica, 2009, 29(2): 288–293 (in Chinese with English abstract).

    [15] 王海雷, 鄭綿平. 青藏高原湖泊水化學(xué)與鹽度的相關(guān)性初步研究[J]. 地質(zhì)學(xué)報(bào), 2010, 84(10): 1517–1522.Wang Hai-lei, Zheng Mian-ping. Preliminary study of the correlation between hydrochemistry and salinity of lakes in the Qinghai-Tibetan Plateau [J]. Acta Geol Sinica, 2010,84(10): 1517–1522 (in Chinese with English abstract).

    [16] 王恒純. 同位素水文地質(zhì)概論[M]. 北京: 地質(zhì)出版社,1991: 1–191.Wang Heng-chun. Introduction to Isotope Hydrogeology [M].Beijing: Geological Publishing House, 1991: 1–191 (in Chinese)

    [17] 周根陶, 鄭永飛. 碳酸鈣水體系氧同位素分餾系數(shù)的低溫實(shí)驗(yàn)研究[J]. 地學(xué)前緣, 2000, 7(2): 321–338.Zhou Gen-tao, Zheng Yong-fei. Experimental studies of oxygen isotope fraction ions between CaCO3and H2O at low temperatures [J]. Earth Sci Front, 2000, 7(2): 321–338.

    [18] 張繼習(xí), 劉耘. 穩(wěn)定同位素分餾蒸汽壓效應(yīng)的計(jì)算方法[J].地球化學(xué), 2014, 43(1): 1–10.Zhang Ji-xi, Liu Yun. Theoretical estimation of vapor pressure isotope effects [J]. Geochimica, 2014, 43(1): 1–10 (in Chinese with English abstract).

    [19] Horita J. Saline waters [M]//Aggarwal P K, Gat J R, Froehlich F O. Isotopes in the Water Cycle. Vienna: IAEA, 2005:271–287.

    [20] 顧慰祖, 龐忠和, 王全九, 宋獻(xiàn)方. 同位素水文學(xué)[M]. 北京: 科學(xué)出版社, 2011: 53–80.Gu Wei-zu, Pang Zhong-he, Wang Quan-jiu, Song Xian-fang.Isotope Hydrology [M]. Beijing: Science Publishing House,2011: 53–80 (in Chinese).

    猜你喜歡
    氫氧水化學(xué)衡水
    幸福湖泉域地下水水化學(xué)特征及其演化機(jī)制
    人民黃河(2025年1期)2025-01-16 00:00:00
    衡水鴻昊企業(yè)有限責(zé)任公司
    肉類研究(2022年7期)2022-08-05 04:47:26
    氫氧燃料電池演示實(shí)驗(yàn)的改進(jìn)
    自制液壓儲氣式氫氧燃料電池
    秦皇島北部地表水和地下水水化學(xué)特征及水質(zhì)評價(jià)
    衡水專場(二)
    大眾文藝(2020年20期)2020-11-05 14:33:52
    連鑄坯氫氧切割應(yīng)用與碳排放研究
    衡水市新近系館陶組地?zé)崴瘜W(xué)特征與分布規(guī)律研究
    筠連煤層氣井產(chǎn)出水化學(xué)特征及意義
    中國煤層氣(2015年5期)2015-08-22 03:26:32
    大推力氫氧火箭發(fā)動機(jī)試驗(yàn)噪聲治理技術(shù)概述
    久久久久久久久免费视频了| 国产午夜精品论理片| 国产激情偷乱视频一区二区| 久久久久免费精品人妻一区二区| 国产精品1区2区在线观看.| 欧美成人午夜精品| av天堂在线播放| 男人舔女人下体高潮全视频| 日韩三级视频一区二区三区| cao死你这个sao货| 欧美人与性动交α欧美精品济南到| 法律面前人人平等表现在哪些方面| 国产精品一区二区免费欧美| 一本大道久久a久久精品| 免费无遮挡裸体视频| 国语自产精品视频在线第100页| 久久国产精品人妻蜜桃| 精品日产1卡2卡| 欧美中文综合在线视频| 99国产精品99久久久久| 午夜激情av网站| 美女大奶头视频| 久久中文看片网| 久久精品综合一区二区三区| 69av精品久久久久久| 国产亚洲av嫩草精品影院| 老鸭窝网址在线观看| 人成视频在线观看免费观看| 国产v大片淫在线免费观看| 国产aⅴ精品一区二区三区波| 亚洲精品av麻豆狂野| 国产伦一二天堂av在线观看| 亚洲精品中文字幕在线视频| 色在线成人网| 长腿黑丝高跟| 欧美在线一区亚洲| 中国美女看黄片| 日本一区二区免费在线视频| 精品无人区乱码1区二区| 怎么达到女性高潮| 亚洲一区二区三区不卡视频| 日本黄色视频三级网站网址| 久久久国产欧美日韩av| 99久久精品热视频| 国产成人影院久久av| 国产视频内射| 啦啦啦观看免费观看视频高清| netflix在线观看网站| 午夜视频精品福利| 欧美黑人精品巨大| 性色av乱码一区二区三区2| 青草久久国产| 色综合欧美亚洲国产小说| 国产成人一区二区三区免费视频网站| 91老司机精品| 变态另类成人亚洲欧美熟女| 在线免费观看的www视频| 国产精品精品国产色婷婷| 免费看十八禁软件| 久久久久性生活片| 在线播放国产精品三级| 国产伦一二天堂av在线观看| 亚洲全国av大片| 午夜a级毛片| 国产av不卡久久| 搡老妇女老女人老熟妇| 午夜福利视频1000在线观看| 亚洲avbb在线观看| 国产伦一二天堂av在线观看| 嫩草影视91久久| 国产又色又爽无遮挡免费看| 搡老岳熟女国产| 欧美另类亚洲清纯唯美| 久久亚洲精品不卡| 亚洲欧美日韩无卡精品| 91字幕亚洲| 成人特级黄色片久久久久久久| 国产亚洲精品一区二区www| 中文字幕高清在线视频| 黄色毛片三级朝国网站| 亚洲国产欧美网| 看黄色毛片网站| 亚洲av中文字字幕乱码综合| 精品国内亚洲2022精品成人| 日韩欧美在线二视频| 欧美成狂野欧美在线观看| 国产一区二区三区在线臀色熟女| 国产男靠女视频免费网站| 老熟妇仑乱视频hdxx| 一级黄色大片毛片| 日日夜夜操网爽| 老司机午夜福利在线观看视频| 日韩大码丰满熟妇| 午夜福利在线观看吧| 亚洲一码二码三码区别大吗| 两个人视频免费观看高清| 最近视频中文字幕2019在线8| 听说在线观看完整版免费高清| 桃红色精品国产亚洲av| 蜜桃久久精品国产亚洲av| 狂野欧美激情性xxxx| 精品国内亚洲2022精品成人| 妹子高潮喷水视频| 国产不卡一卡二| a级毛片a级免费在线| 午夜福利免费观看在线| 国产黄色小视频在线观看| 真人一进一出gif抽搐免费| 国产单亲对白刺激| 日日夜夜操网爽| 在线a可以看的网站| 人人妻,人人澡人人爽秒播| 国产免费av片在线观看野外av| 亚洲熟妇熟女久久| 久久九九热精品免费| 狂野欧美激情性xxxx| 精品国内亚洲2022精品成人| 一二三四在线观看免费中文在| 97碰自拍视频| 岛国视频午夜一区免费看| 中文字幕人妻丝袜一区二区| 19禁男女啪啪无遮挡网站| 亚洲精品久久国产高清桃花| 成年免费大片在线观看| 国内精品久久久久精免费| 国产亚洲精品久久久久久毛片| av中文乱码字幕在线| 五月玫瑰六月丁香| 99久久无色码亚洲精品果冻| 亚洲一区高清亚洲精品| 丁香六月欧美| 亚洲精品在线美女| 香蕉av资源在线| 国产黄色小视频在线观看| 99热6这里只有精品| 两个人看的免费小视频| 又紧又爽又黄一区二区| 午夜a级毛片| 国内毛片毛片毛片毛片毛片| 日韩高清综合在线| www.熟女人妻精品国产| 曰老女人黄片| 亚洲乱码一区二区免费版| 99精品欧美一区二区三区四区| 性色av乱码一区二区三区2| 久久久久性生活片| 亚洲av成人一区二区三| 亚洲全国av大片| 99精品欧美一区二区三区四区| 国产成人精品无人区| 97超级碰碰碰精品色视频在线观看| 99热这里只有是精品50| 久热爱精品视频在线9| 亚洲国产精品sss在线观看| 色综合站精品国产| 日韩有码中文字幕| 啦啦啦观看免费观看视频高清| 精品少妇一区二区三区视频日本电影| 亚洲片人在线观看| 亚洲一码二码三码区别大吗| 欧美中文综合在线视频| 又黄又爽又免费观看的视频| 久久精品夜夜夜夜夜久久蜜豆 | 伊人久久大香线蕉亚洲五| 美女黄网站色视频| 国产日本99.免费观看| 欧美久久黑人一区二区| 麻豆成人av在线观看| 一a级毛片在线观看| 露出奶头的视频| 欧美在线一区亚洲| 亚洲欧美精品综合久久99| 亚洲午夜理论影院| 久久欧美精品欧美久久欧美| 久久久久久大精品| 热99re8久久精品国产| 欧美最黄视频在线播放免费| 三级国产精品欧美在线观看 | 三级男女做爰猛烈吃奶摸视频| 亚洲欧美日韩无卡精品| 天天添夜夜摸| 人妻久久中文字幕网| 人人妻人人看人人澡| 国产高清视频在线播放一区| 国产亚洲欧美在线一区二区| tocl精华| 脱女人内裤的视频| 色哟哟哟哟哟哟| 亚洲av美国av| 狠狠狠狠99中文字幕| 老司机午夜福利在线观看视频| 国产熟女xx| 伊人久久大香线蕉亚洲五| 九色成人免费人妻av| 91国产中文字幕| 国产片内射在线| 欧美日本视频| 婷婷精品国产亚洲av在线| 久久久国产精品麻豆| 免费在线观看黄色视频的| 免费在线观看日本一区| 国产亚洲精品一区二区www| 国产男靠女视频免费网站| 亚洲人成77777在线视频| 午夜精品久久久久久毛片777| 很黄的视频免费| 亚洲av中文字字幕乱码综合| 精品电影一区二区在线| 国产精品自产拍在线观看55亚洲| 色av中文字幕| 亚洲av五月六月丁香网| av欧美777| 老鸭窝网址在线观看| 精品久久蜜臀av无| 亚洲专区国产一区二区| 变态另类丝袜制服| 亚洲人成伊人成综合网2020| 国产亚洲精品一区二区www| 亚洲美女黄片视频| 成人午夜高清在线视频| 国产久久久一区二区三区| 少妇熟女aⅴ在线视频| 啦啦啦免费观看视频1| 免费搜索国产男女视频| 久久国产精品人妻蜜桃| 国产成+人综合+亚洲专区| 亚洲18禁久久av| 午夜精品久久久久久毛片777| 天堂av国产一区二区熟女人妻 | 亚洲熟女毛片儿| 90打野战视频偷拍视频| 久99久视频精品免费| 久久久久久亚洲精品国产蜜桃av| 丁香六月欧美| 国产精品影院久久| 国产亚洲精品久久久久5区| 国产成人欧美在线观看| 国产精品98久久久久久宅男小说| 丰满的人妻完整版| 午夜精品在线福利| 欧美极品一区二区三区四区| 亚洲人成网站在线播放欧美日韩| 欧美成人性av电影在线观看| 亚洲男人的天堂狠狠| 国产精品一区二区精品视频观看| 波多野结衣高清作品| 国产熟女午夜一区二区三区| 亚洲午夜理论影院| 色综合亚洲欧美另类图片| 最近最新免费中文字幕在线| 日日夜夜操网爽| 欧美又色又爽又黄视频| 国产亚洲精品久久久久5区| 熟女少妇亚洲综合色aaa.| 99精品在免费线老司机午夜| 精品少妇一区二区三区视频日本电影| av有码第一页| 久久久久性生活片| 亚洲国产精品合色在线| 久久久久久久久久黄片| 午夜激情福利司机影院| 亚洲自偷自拍图片 自拍| 久久国产乱子伦精品免费另类| 麻豆成人午夜福利视频| 欧美乱码精品一区二区三区| 午夜福利欧美成人| 毛片女人毛片| 激情在线观看视频在线高清| 亚洲成人国产一区在线观看| 免费在线观看亚洲国产| 国产精品一区二区精品视频观看| 黄色女人牲交| 久久久久九九精品影院| 国产精品电影一区二区三区| 亚洲美女黄片视频| 桃红色精品国产亚洲av| 日日摸夜夜添夜夜添小说| 精品电影一区二区在线| 色精品久久人妻99蜜桃| 窝窝影院91人妻| 久久人人精品亚洲av| 人成视频在线观看免费观看| 黑人欧美特级aaaaaa片| 老司机靠b影院| 免费电影在线观看免费观看| 国产成人av教育| 婷婷丁香在线五月| 久久婷婷成人综合色麻豆| 亚洲熟妇熟女久久| 欧美乱妇无乱码| 国产99久久九九免费精品| 亚洲av成人不卡在线观看播放网| 国产成人精品久久二区二区91| 最近最新中文字幕大全电影3| 我的老师免费观看完整版| 极品教师在线免费播放| 一本精品99久久精品77| 每晚都被弄得嗷嗷叫到高潮| 90打野战视频偷拍视频| 久久草成人影院| 亚洲第一欧美日韩一区二区三区| 黄色成人免费大全| 最近视频中文字幕2019在线8| 精品久久久久久成人av| 国产乱人伦免费视频| 久久天堂一区二区三区四区| 午夜福利在线观看吧| 午夜福利成人在线免费观看| 日本五十路高清| 久久久国产成人免费| 黑人巨大精品欧美一区二区mp4| 久久久久久国产a免费观看| 精品第一国产精品| 亚洲激情在线av| 成熟少妇高潮喷水视频| 婷婷精品国产亚洲av在线| 免费在线观看黄色视频的| 不卡一级毛片| 久久性视频一级片| 亚洲中文字幕日韩| 免费看十八禁软件| 国产精品99久久99久久久不卡| 亚洲精品中文字幕一二三四区| 亚洲av片天天在线观看| 亚洲熟妇中文字幕五十中出| 99热这里只有是精品50| 1024香蕉在线观看| 最近最新中文字幕大全电影3| 日韩精品免费视频一区二区三区| 亚洲欧美激情综合另类| 午夜影院日韩av| 一级作爱视频免费观看| 麻豆成人av在线观看| 亚洲,欧美精品.| 日本三级黄在线观看| 日本黄色视频三级网站网址| 久久久精品欧美日韩精品| 熟女少妇亚洲综合色aaa.| 亚洲熟女毛片儿| 一级片免费观看大全| 国产激情偷乱视频一区二区| 国产高清激情床上av| 亚洲成人久久爱视频| 欧美精品啪啪一区二区三区| 十八禁网站免费在线| 久久久久久亚洲精品国产蜜桃av| 999精品在线视频| 免费电影在线观看免费观看| 无人区码免费观看不卡| 久久国产精品影院| 极品教师在线免费播放| 国内精品久久久久精免费| 午夜福利18| 俺也久久电影网| 中文字幕高清在线视频| 少妇被粗大的猛进出69影院| 黄色片一级片一级黄色片| 亚洲欧美日韩东京热| 草草在线视频免费看| 无人区码免费观看不卡| 国产亚洲精品久久久久久毛片| 国产精品爽爽va在线观看网站| 色综合欧美亚洲国产小说| 国产亚洲av高清不卡| 国产一区二区三区在线臀色熟女| 久久人妻av系列| 欧美av亚洲av综合av国产av| 青草久久国产| 国产精品国产高清国产av| 亚洲,欧美精品.| 可以免费在线观看a视频的电影网站| 人人妻人人看人人澡| 午夜免费观看网址| 又爽又黄无遮挡网站| 久久婷婷成人综合色麻豆| 69av精品久久久久久| 天堂av国产一区二区熟女人妻 | 好男人在线观看高清免费视频| 特级一级黄色大片| 中文字幕人妻丝袜一区二区| 蜜桃久久精品国产亚洲av| 欧美+亚洲+日韩+国产| 免费看十八禁软件| 亚洲午夜理论影院| 欧美午夜高清在线| 又爽又黄无遮挡网站| 黄色 视频免费看| 亚洲黑人精品在线| 亚洲成av人片在线播放无| 欧美在线黄色| 久久午夜亚洲精品久久| 欧美一区二区精品小视频在线| 久久精品成人免费网站| 少妇熟女aⅴ在线视频| 久久久久久国产a免费观看| 国产精品永久免费网站| 人人妻,人人澡人人爽秒播| 亚洲狠狠婷婷综合久久图片| 亚洲av成人一区二区三| 99精品欧美一区二区三区四区| 成人三级做爰电影| 久久欧美精品欧美久久欧美| 亚洲欧美日韩高清专用| 精品国内亚洲2022精品成人| 黄色 视频免费看| 亚洲精华国产精华精| 黑人操中国人逼视频| 宅男免费午夜| 欧美日韩亚洲国产一区二区在线观看| 国产精品影院久久| 51午夜福利影视在线观看| 成人18禁高潮啪啪吃奶动态图| 五月伊人婷婷丁香| 久久久精品国产亚洲av高清涩受| 老汉色∧v一级毛片| 亚洲aⅴ乱码一区二区在线播放 | 久久久久久大精品| 午夜久久久久精精品| 搡老岳熟女国产| 国产成人一区二区三区免费视频网站| 精品不卡国产一区二区三区| 欧美在线一区亚洲| 久久久久久久久免费视频了| 18禁观看日本| av超薄肉色丝袜交足视频| 亚洲国产日韩欧美精品在线观看 | 少妇的丰满在线观看| 国产伦一二天堂av在线观看| 90打野战视频偷拍视频| 亚洲人成网站高清观看| 18禁黄网站禁片免费观看直播| 91av网站免费观看| 黄色a级毛片大全视频| av欧美777| 国产日本99.免费观看| 亚洲专区国产一区二区| 亚洲欧美一区二区三区黑人| 在线国产一区二区在线| 1024香蕉在线观看| 欧美另类亚洲清纯唯美| av中文乱码字幕在线| tocl精华| 亚洲国产高清在线一区二区三| 亚洲七黄色美女视频| 天堂av国产一区二区熟女人妻 | 亚洲精品美女久久av网站| 国产久久久一区二区三区| 午夜福利18| 在线观看免费午夜福利视频| 精品国产亚洲在线| 亚洲国产精品成人综合色| 亚洲精品av麻豆狂野| 亚洲av美国av| 最近最新中文字幕大全电影3| 99re在线观看精品视频| 99久久国产精品久久久| 久久国产精品影院| 亚洲欧美一区二区三区黑人| 欧美乱码精品一区二区三区| 亚洲精品国产精品久久久不卡| 国产人伦9x9x在线观看| www日本黄色视频网| bbb黄色大片| 一级黄色大片毛片| 人成视频在线观看免费观看| 好看av亚洲va欧美ⅴa在| 日韩欧美三级三区| 成人18禁在线播放| av天堂在线播放| 亚洲欧美一区二区三区黑人| 久久 成人 亚洲| 日韩精品青青久久久久久| 国产亚洲精品久久久久5区| 99国产精品99久久久久| 久久久久久久久久黄片| 欧美日韩亚洲国产一区二区在线观看| www.自偷自拍.com| 久久久国产成人免费| 国产精品亚洲一级av第二区| av视频在线观看入口| 久久久久久久久中文| 国产成年人精品一区二区| 色尼玛亚洲综合影院| 精华霜和精华液先用哪个| 国产日本99.免费观看| 欧美+亚洲+日韩+国产| www国产在线视频色| 国产精品永久免费网站| 99re在线观看精品视频| 一本精品99久久精品77| 亚洲精品美女久久av网站| 欧美乱色亚洲激情| www.www免费av| 午夜福利高清视频| 午夜福利成人在线免费观看| 国产欧美日韩精品亚洲av| 国产精品亚洲美女久久久| 日本 欧美在线| 午夜福利高清视频| 欧美另类亚洲清纯唯美| 国产黄a三级三级三级人| 欧美性长视频在线观看| 桃红色精品国产亚洲av| 男女午夜视频在线观看| 久久午夜综合久久蜜桃| 欧美性猛交黑人性爽| 在线十欧美十亚洲十日本专区| 青草久久国产| 精品久久久久久,| 久久久久精品国产欧美久久久| 久久精品综合一区二区三区| 免费无遮挡裸体视频| 在线观看日韩欧美| 色噜噜av男人的天堂激情| 高潮久久久久久久久久久不卡| 日韩欧美免费精品| 熟女电影av网| 成年免费大片在线观看| 精品乱码久久久久久99久播| 91在线观看av| 一本久久中文字幕| 国产野战对白在线观看| 观看免费一级毛片| 精品国产乱码久久久久久男人| 麻豆成人av在线观看| 久久这里只有精品19| 少妇被粗大的猛进出69影院| 99久久99久久久精品蜜桃| 母亲3免费完整高清在线观看| 悠悠久久av| 成熟少妇高潮喷水视频| 女生性感内裤真人,穿戴方法视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产99精品国产亚洲性色| 丰满人妻一区二区三区视频av | 狂野欧美激情性xxxx| 国产成年人精品一区二区| 熟妇人妻久久中文字幕3abv| 国产精华一区二区三区| 亚洲国产精品999在线| 九色国产91popny在线| 我要搜黄色片| 人妻夜夜爽99麻豆av| 久久草成人影院| 久久精品影院6| 成人国产综合亚洲| 亚洲人与动物交配视频| 国产主播在线观看一区二区| 色综合婷婷激情| 亚洲五月天丁香| 神马国产精品三级电影在线观看 | 国产精品九九99| 午夜激情av网站| 国内精品久久久久久久电影| 免费av毛片视频| 一本综合久久免费| 亚洲成av人片免费观看| 成人三级黄色视频| 可以在线观看毛片的网站| 欧美黑人巨大hd| 午夜免费成人在线视频| 中文字幕最新亚洲高清| 制服丝袜大香蕉在线| 我的老师免费观看完整版| 精品久久久久久成人av| 亚洲av电影在线进入| 亚洲成人免费电影在线观看| 国产单亲对白刺激| 日韩三级视频一区二区三区| 国产精品1区2区在线观看.| 99久久精品国产亚洲精品| 亚洲男人天堂网一区| 亚洲国产精品合色在线| 三级毛片av免费| 宅男免费午夜| 亚洲欧美日韩东京热| 可以在线观看的亚洲视频| 免费av毛片视频| 欧美乱色亚洲激情| 一级作爱视频免费观看| 国产免费av片在线观看野外av| 女人爽到高潮嗷嗷叫在线视频| 在线观看免费日韩欧美大片| 在线视频色国产色| 亚洲精品av麻豆狂野| 久久精品91无色码中文字幕| av天堂在线播放| 夜夜躁狠狠躁天天躁| 岛国在线免费视频观看| 日韩欧美国产一区二区入口| 一本精品99久久精品77| 狠狠狠狠99中文字幕| 嫩草影视91久久| 男女那种视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | av片东京热男人的天堂| 岛国视频午夜一区免费看| 欧美性长视频在线观看| www.www免费av| 最好的美女福利视频网| 久久国产精品人妻蜜桃| 69av精品久久久久久| 国产av在哪里看| 观看免费一级毛片| 成人午夜高清在线视频| 日韩欧美 国产精品| 亚洲欧洲精品一区二区精品久久久| 黄频高清免费视频| 亚洲 国产 在线| 久久久久久久久免费视频了| 免费搜索国产男女视频| 欧美日本亚洲视频在线播放| 国产激情欧美一区二区| 久久久久九九精品影院| 男女下面进入的视频免费午夜|