• <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ù)概述
    麻豆一二三区av精品| 天天影视国产精品| 曰老女人黄片| 久久久久久久午夜电影 | 在线观看免费视频网站a站| 亚洲 国产 在线| 每晚都被弄得嗷嗷叫到高潮| 国产99白浆流出| 少妇粗大呻吟视频| 激情视频va一区二区三区| 欧美亚洲日本最大视频资源| 欧美日韩精品网址| 老司机亚洲免费影院| av免费在线观看网站| 777久久人妻少妇嫩草av网站| 夜夜看夜夜爽夜夜摸 | 久久亚洲精品不卡| 亚洲精品一区av在线观看| 久久香蕉精品热| 黄网站色视频无遮挡免费观看| 久热爱精品视频在线9| 亚洲人成网站在线播放欧美日韩| 久久香蕉精品热| 91字幕亚洲| 男女午夜视频在线观看| 亚洲成人精品中文字幕电影 | 精品免费久久久久久久清纯| 一级,二级,三级黄色视频| 香蕉国产在线看| 很黄的视频免费| 老司机福利观看| 女人爽到高潮嗷嗷叫在线视频| 午夜精品国产一区二区电影| 久久久久久大精品| 久久这里只有精品19| 免费在线观看日本一区| 天天躁夜夜躁狠狠躁躁| 五月开心婷婷网| 波多野结衣一区麻豆| 久久影院123| 变态另类成人亚洲欧美熟女 | 国产精品av久久久久免费| 日韩人妻精品一区2区三区| 在线看a的网站| 久久人妻熟女aⅴ| 亚洲国产欧美一区二区综合| 一区二区三区精品91| 国产熟女午夜一区二区三区| 国产麻豆69| 中文字幕色久视频| 成人18禁高潮啪啪吃奶动态图| 可以在线观看毛片的网站| 香蕉久久夜色| 欧美精品亚洲一区二区| 成人av一区二区三区在线看| 久久国产乱子伦精品免费另类| 在线免费观看的www视频| 在线观看午夜福利视频| 高潮久久久久久久久久久不卡| 日本黄色视频三级网站网址| 19禁男女啪啪无遮挡网站| 99香蕉大伊视频| 国产99久久九九免费精品| 久久精品亚洲精品国产色婷小说| 中文字幕av电影在线播放| 国产午夜精品久久久久久| av中文乱码字幕在线| 纯流量卡能插随身wifi吗| 老司机靠b影院| 久久久久久免费高清国产稀缺| 亚洲精品国产区一区二| 亚洲国产精品合色在线| 国产免费男女视频| 99国产精品99久久久久| www日本在线高清视频| 一级a爱视频在线免费观看| 激情视频va一区二区三区| 波多野结衣高清无吗| 国产精品久久久av美女十八| 亚洲久久久国产精品| 人成视频在线观看免费观看| www.熟女人妻精品国产| 色哟哟哟哟哟哟| 老汉色av国产亚洲站长工具| 极品人妻少妇av视频| 男男h啪啪无遮挡| 97碰自拍视频| 国产精品久久久久久人妻精品电影| a在线观看视频网站| 国产日韩一区二区三区精品不卡| 亚洲 欧美 日韩 在线 免费| 90打野战视频偷拍视频| 国产精品久久电影中文字幕| 亚洲欧美日韩另类电影网站| 一进一出好大好爽视频| 一级作爱视频免费观看| 国产免费男女视频| 欧美黑人精品巨大| 真人一进一出gif抽搐免费| 一级毛片女人18水好多| 午夜影院日韩av| 国产伦人伦偷精品视频| 国产高清国产精品国产三级| 夜夜看夜夜爽夜夜摸 | 国产精华一区二区三区| 亚洲男人的天堂狠狠| 国产精品一区二区在线不卡| 91精品三级在线观看| 在线免费观看的www视频| 国产精品自产拍在线观看55亚洲| 天天添夜夜摸| 一级片免费观看大全| 如日韩欧美国产精品一区二区三区| 久久久国产成人精品二区 | 亚洲 国产 在线| 精品人妻在线不人妻| 日韩三级视频一区二区三区| 高清黄色对白视频在线免费看| 一区在线观看完整版| 丰满人妻熟妇乱又伦精品不卡| 日本黄色视频三级网站网址| 叶爱在线成人免费视频播放| 色综合欧美亚洲国产小说| 9热在线视频观看99| 精品国产美女av久久久久小说| 十八禁人妻一区二区| xxxhd国产人妻xxx| 欧美日本亚洲视频在线播放| 久久久久九九精品影院| 国产熟女xx| 天天影视国产精品| 一区在线观看完整版| 精品卡一卡二卡四卡免费| 亚洲成人精品中文字幕电影 | 韩国av一区二区三区四区| 午夜免费成人在线视频| videosex国产| 久久精品成人免费网站| 美女 人体艺术 gogo| 国产成人欧美在线观看| 久久人妻av系列| 又紧又爽又黄一区二区| 悠悠久久av| 国产不卡一卡二| 在线天堂中文资源库| 狠狠狠狠99中文字幕| 欧美人与性动交α欧美精品济南到| 狂野欧美激情性xxxx| 亚洲 欧美一区二区三区| 久久久久国产精品人妻aⅴ院| 久久人人97超碰香蕉20202| 久久中文看片网| 久久精品国产亚洲av高清一级| 精品熟女少妇八av免费久了| 色综合欧美亚洲国产小说| 亚洲专区字幕在线| 天天躁夜夜躁狠狠躁躁| 美女高潮喷水抽搐中文字幕| 日本a在线网址| 香蕉国产在线看| 色婷婷av一区二区三区视频| 可以在线观看毛片的网站| 欧美乱色亚洲激情| 国产精品久久久av美女十八| 亚洲七黄色美女视频| 国产亚洲欧美98| 国产亚洲精品一区二区www| 精品人妻1区二区| av欧美777| 在线看a的网站| 成人亚洲精品av一区二区 | 日韩视频一区二区在线观看| 一进一出抽搐动态| 久久久久久久久中文| 亚洲专区中文字幕在线| 欧美黑人欧美精品刺激| 久久性视频一级片| 国产成人啪精品午夜网站| 亚洲午夜精品一区,二区,三区| 美女午夜性视频免费| 一级片免费观看大全| 亚洲avbb在线观看| 岛国在线观看网站| 亚洲精品在线观看二区| 午夜福利免费观看在线| 精品无人区乱码1区二区| 一级毛片精品| 999精品在线视频| 长腿黑丝高跟| 一区二区三区精品91| 黑人欧美特级aaaaaa片| 亚洲 欧美 日韩 在线 免费| 久久国产乱子伦精品免费另类| 欧美日韩中文字幕国产精品一区二区三区 | 色婷婷av一区二区三区视频| 免费观看精品视频网站| 精品国产亚洲在线| 国产一卡二卡三卡精品| 老鸭窝网址在线观看| 精品日产1卡2卡| 久久精品亚洲精品国产色婷小说| 最好的美女福利视频网| 一个人免费在线观看的高清视频| 女警被强在线播放| 99精品在免费线老司机午夜| 超碰97精品在线观看| 亚洲精品一二三| 狠狠狠狠99中文字幕| 很黄的视频免费| 丰满的人妻完整版| 一级a爱片免费观看的视频| 亚洲熟妇中文字幕五十中出 | 午夜激情av网站| avwww免费| 亚洲七黄色美女视频| 可以免费在线观看a视频的电影网站| 亚洲av成人一区二区三| 久久久久亚洲av毛片大全| 国产成人精品在线电影| 亚洲欧美激情综合另类| 午夜91福利影院| 免费看a级黄色片| 极品教师在线免费播放| 日韩一卡2卡3卡4卡2021年| a级毛片黄视频| 久久狼人影院| 亚洲精品中文字幕一二三四区| 亚洲人成伊人成综合网2020| 亚洲精品国产一区二区精华液| 美国免费a级毛片| 夜夜看夜夜爽夜夜摸 | 国产色视频综合| 9色porny在线观看| 亚洲熟女毛片儿| 久久精品国产清高在天天线| 日韩大码丰满熟妇| 亚洲在线自拍视频| 妹子高潮喷水视频| 亚洲熟妇中文字幕五十中出 | 国产欧美日韩一区二区三区在线| 国产亚洲欧美在线一区二区| 国产aⅴ精品一区二区三区波| 精品人妻1区二区| 母亲3免费完整高清在线观看| 正在播放国产对白刺激| 久久香蕉精品热| 亚洲,欧美精品.| 免费女性裸体啪啪无遮挡网站| 视频区图区小说| 淫妇啪啪啪对白视频| 999精品在线视频| av有码第一页| 国产精品秋霞免费鲁丝片| 成人永久免费在线观看视频| 波多野结衣av一区二区av| 免费观看人在逋| 亚洲精品国产一区二区精华液| 变态另类成人亚洲欧美熟女 | 成人免费观看视频高清| aaaaa片日本免费| 亚洲专区字幕在线| 午夜福利欧美成人| 国产精品 国内视频| netflix在线观看网站| 国产精品二区激情视频| av天堂久久9| 免费少妇av软件| 国产精品美女特级片免费视频播放器 | 十八禁人妻一区二区| 日韩成人在线观看一区二区三区| 99re在线观看精品视频| 欧美激情高清一区二区三区| 少妇的丰满在线观看| 欧美日韩瑟瑟在线播放| 丁香欧美五月| 日韩av在线大香蕉| 精品日产1卡2卡| 欧美乱色亚洲激情| 村上凉子中文字幕在线| 欧美日韩亚洲综合一区二区三区_| 一边摸一边抽搐一进一小说| 亚洲免费av在线视频| 黄片播放在线免费| 看片在线看免费视频| 精品福利观看| 在线国产一区二区在线| 伊人久久大香线蕉亚洲五| 欧美日韩乱码在线| √禁漫天堂资源中文www| 中文亚洲av片在线观看爽| 国产成年人精品一区二区 | 美女福利国产在线| 色在线成人网| 久久久久久久久久久久大奶| 国产日韩一区二区三区精品不卡| 久久人人97超碰香蕉20202| 男女之事视频高清在线观看| 韩国av一区二区三区四区| 十八禁网站免费在线| 法律面前人人平等表现在哪些方面| av中文乱码字幕在线| 亚洲人成电影免费在线| 欧美日韩亚洲高清精品| 久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 亚洲一区二区三区欧美精品| 变态另类成人亚洲欧美熟女 | 一级片'在线观看视频| 亚洲少妇的诱惑av| 久久久精品国产亚洲av高清涩受| 韩国精品一区二区三区| 18禁黄网站禁片午夜丰满| 欧美日韩黄片免| 久久久久国内视频| 国产视频一区二区在线看| 亚洲精品中文字幕一二三四区| 午夜91福利影院| 亚洲中文av在线| 国产成人av激情在线播放| 久热这里只有精品99| 自拍欧美九色日韩亚洲蝌蚪91| 最好的美女福利视频网| 亚洲成a人片在线一区二区| 日本免费一区二区三区高清不卡 | 操出白浆在线播放| 不卡一级毛片| 一级片免费观看大全| 男人操女人黄网站| 国产一区二区三区综合在线观看| 777久久人妻少妇嫩草av网站| 免费少妇av软件| 一区在线观看完整版| 国产精品一区二区精品视频观看| 精品国产一区二区三区四区第35| 欧美日韩视频精品一区| 免费在线观看日本一区| 日本 av在线| 两性夫妻黄色片| 亚洲av片天天在线观看| 欧美在线黄色| 男人操女人黄网站| 中出人妻视频一区二区| 999精品在线视频| 精品免费久久久久久久清纯| 国产日韩一区二区三区精品不卡| 亚洲人成网站在线播放欧美日韩| 欧美黑人欧美精品刺激| 成人国语在线视频| 黄色成人免费大全| 亚洲熟妇熟女久久| 欧美成人免费av一区二区三区| 亚洲精品中文字幕在线视频| 亚洲国产精品一区二区三区在线| 免费在线观看日本一区| 国产精品香港三级国产av潘金莲| 法律面前人人平等表现在哪些方面| 啦啦啦免费观看视频1| 亚洲第一欧美日韩一区二区三区| 国产在线观看jvid| 国产亚洲欧美精品永久| 日韩欧美一区视频在线观看| 免费av中文字幕在线| 一区福利在线观看| 国产精品 欧美亚洲| 国产高清激情床上av| 婷婷丁香在线五月| 国产高清激情床上av| 国产激情久久老熟女| 国产xxxxx性猛交| 亚洲欧美精品综合久久99| 午夜免费成人在线视频| 欧美黄色淫秽网站| 亚洲第一av免费看| 精品久久久久久成人av| 一进一出好大好爽视频| 91老司机精品| 91精品三级在线观看| 久久热在线av| 18禁国产床啪视频网站| 欧美黄色淫秽网站| 久久久久久久久中文| 啦啦啦 在线观看视频| 99精品在免费线老司机午夜| 久久精品人人爽人人爽视色| 在线国产一区二区在线| 日韩欧美免费精品| 美女国产高潮福利片在线看| 亚洲五月天丁香| 变态另类成人亚洲欧美熟女 | 麻豆国产av国片精品| 不卡一级毛片| 18禁黄网站禁片午夜丰满| 大型黄色视频在线免费观看| 色婷婷av一区二区三区视频| 日本欧美视频一区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲视频免费观看视频| 日韩大尺度精品在线看网址 | 黄色片一级片一级黄色片| 97碰自拍视频| www国产在线视频色| 一级a爱片免费观看的视频| 精品日产1卡2卡| 国产高清激情床上av| 久久久久精品国产欧美久久久| 久久天堂一区二区三区四区| 人妻久久中文字幕网| ponron亚洲| 中文欧美无线码| 免费在线观看亚洲国产| 制服人妻中文乱码| 精品人妻在线不人妻| 91麻豆av在线| 天天躁夜夜躁狠狠躁躁| 如日韩欧美国产精品一区二区三区| 老司机深夜福利视频在线观看| 国产午夜精品久久久久久| av天堂久久9| 亚洲欧美日韩无卡精品| 一区在线观看完整版| www日本在线高清视频| 激情在线观看视频在线高清| 精品乱码久久久久久99久播| 亚洲精品国产一区二区精华液| 黑人巨大精品欧美一区二区蜜桃| 视频在线观看一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 日本精品一区二区三区蜜桃| 久久草成人影院| 国产99白浆流出| 精品国产乱码久久久久久男人| 久久人人爽av亚洲精品天堂| 亚洲国产精品999在线| 女警被强在线播放| 国产99白浆流出| 1024视频免费在线观看| 亚洲专区字幕在线| 亚洲av日韩精品久久久久久密| 黄片播放在线免费| 国产蜜桃级精品一区二区三区| 日本vs欧美在线观看视频| 久久婷婷成人综合色麻豆| 亚洲五月天丁香| 夜夜看夜夜爽夜夜摸 | 免费观看精品视频网站| 成在线人永久免费视频| 老司机亚洲免费影院| 中国美女看黄片| 身体一侧抽搐| 国产亚洲欧美98| 国产一区在线观看成人免费| 丰满迷人的少妇在线观看| 19禁男女啪啪无遮挡网站| 在线观看免费高清a一片| 欧美黄色片欧美黄色片| 黑人猛操日本美女一级片| 国产真人三级小视频在线观看| 最好的美女福利视频网| 欧美日韩黄片免| 亚洲一区二区三区欧美精品| 老司机亚洲免费影院| cao死你这个sao货| 亚洲av五月六月丁香网| 99国产精品免费福利视频| 久久青草综合色| 啦啦啦 在线观看视频| 欧美最黄视频在线播放免费 | 18禁美女被吸乳视频| 色播在线永久视频| www.999成人在线观看| 伦理电影免费视频| 男人舔女人下体高潮全视频| 精品免费久久久久久久清纯| 十八禁人妻一区二区| 精品国产一区二区三区四区第35| av片东京热男人的天堂| 亚洲 国产 在线| 国产成人av教育| 免费在线观看影片大全网站| 亚洲色图 男人天堂 中文字幕| 亚洲欧美激情在线| 一区二区三区精品91| 国产视频一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 韩国av一区二区三区四区| 91精品三级在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产成年人精品一区二区 | 韩国精品一区二区三区| 午夜免费激情av| 激情在线观看视频在线高清| 日韩视频一区二区在线观看| 丁香六月欧美| 午夜免费观看网址| 欧洲精品卡2卡3卡4卡5卡区| 在线观看一区二区三区激情| 午夜精品久久久久久毛片777| 午夜激情av网站| 久久99一区二区三区| 美国免费a级毛片| 黑丝袜美女国产一区| 丁香六月欧美| 成年人免费黄色播放视频| 国产在线精品亚洲第一网站| av视频免费观看在线观看| 熟女少妇亚洲综合色aaa.| 在线观看一区二区三区激情| 国产成人欧美在线观看| 搡老岳熟女国产| 91麻豆av在线| 日韩精品免费视频一区二区三区| 亚洲国产中文字幕在线视频| 亚洲视频免费观看视频| 久久人人爽av亚洲精品天堂| 99国产精品免费福利视频| 国产精品电影一区二区三区| 美女高潮到喷水免费观看| 久久精品国产综合久久久| 国产亚洲精品久久久久5区| 色婷婷久久久亚洲欧美| 国产精品久久视频播放| 亚洲欧美精品综合一区二区三区| a在线观看视频网站| 热re99久久精品国产66热6| 国产av又大| 人人妻,人人澡人人爽秒播| 在线av久久热| 亚洲片人在线观看| 最新美女视频免费是黄的| 国产精品二区激情视频| 久久 成人 亚洲| 国产精品爽爽va在线观看网站 | 一区二区日韩欧美中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区二区三卡| 亚洲第一欧美日韩一区二区三区| 男人舔女人的私密视频| 69精品国产乱码久久久| 久久国产精品影院| 久久久久国内视频| 亚洲中文字幕日韩| 两性夫妻黄色片| 999精品在线视频| 成人国语在线视频| 日本精品一区二区三区蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 久久午夜亚洲精品久久| 啦啦啦免费观看视频1| 国产99久久九九免费精品| 天堂中文最新版在线下载| 国产无遮挡羞羞视频在线观看| 成年女人毛片免费观看观看9| 亚洲精品美女久久久久99蜜臀| 无限看片的www在线观看| 欧美精品一区二区免费开放| 亚洲精品一区av在线观看| 国产日韩一区二区三区精品不卡| 国产亚洲精品久久久久5区| 国产99白浆流出| 自拍欧美九色日韩亚洲蝌蚪91| 黄色视频不卡| 亚洲精品一卡2卡三卡4卡5卡| 老司机午夜十八禁免费视频| 免费在线观看完整版高清| 男男h啪啪无遮挡| 极品教师在线免费播放| 国产精品国产高清国产av| 十八禁人妻一区二区| 老熟妇仑乱视频hdxx| 女人精品久久久久毛片| 啦啦啦免费观看视频1| 国产高清国产精品国产三级| 亚洲午夜理论影院| 色综合站精品国产| 久久精品亚洲熟妇少妇任你| 欧美在线一区亚洲| 国产精品久久电影中文字幕| 国产男靠女视频免费网站| 法律面前人人平等表现在哪些方面| 国产主播在线观看一区二区| 50天的宝宝边吃奶边哭怎么回事| 日韩视频一区二区在线观看| 老司机在亚洲福利影院| tocl精华| 乱人伦中国视频| 熟女少妇亚洲综合色aaa.| 在线永久观看黄色视频| 免费在线观看日本一区| 嫩草影院精品99| 久久久久国产精品人妻aⅴ院| 丝袜人妻中文字幕| 免费在线观看影片大全网站| 久久久国产成人精品二区 | 一区在线观看完整版| 新久久久久国产一级毛片| 成人18禁在线播放| 19禁男女啪啪无遮挡网站| 热99国产精品久久久久久7| 中文欧美无线码| 久久久国产欧美日韩av| 久久午夜综合久久蜜桃| 99国产精品免费福利视频| 韩国av一区二区三区四区| 久久天躁狠狠躁夜夜2o2o| 12—13女人毛片做爰片一| 最新在线观看一区二区三区| 久久天堂一区二区三区四区| 免费观看精品视频网站| 日本免费a在线| 国产亚洲av高清不卡| 99热国产这里只有精品6| 国产真人三级小视频在线观看| 成人18禁高潮啪啪吃奶动态图| 91成年电影在线观看| 777久久人妻少妇嫩草av网站|