• <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ù)概述
    欧美亚洲日本最大视频资源| 看十八女毛片水多多多| 汤姆久久久久久久影院中文字幕| 十分钟在线观看高清视频www| 激情五月婷婷亚洲| 一边摸一边做爽爽视频免费| 国产成人精品一,二区| www.av在线官网国产| av在线观看视频网站免费| 日韩中文字幕视频在线看片| 99久国产av精品国产电影| 欧美成人精品欧美一级黄| 精品人妻在线不人妻| 成人亚洲精品一区在线观看| 国产一区二区在线观看av| 久久青草综合色| 一区二区三区激情视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品日韩在线中文字幕| 大码成人一级视频| 国产在视频线精品| 巨乳人妻的诱惑在线观看| 欧美最新免费一区二区三区| 女人精品久久久久毛片| 欧美人与性动交α欧美软件| 国产一区有黄有色的免费视频| 人人妻人人澡人人爽人人夜夜| 国产熟女午夜一区二区三区| 久久国产亚洲av麻豆专区| 久久久久久免费高清国产稀缺| 亚洲人成网站在线观看播放| 免费高清在线观看视频在线观看| 国产综合精华液| 曰老女人黄片| 巨乳人妻的诱惑在线观看| 日本免费在线观看一区| 国产精品人妻久久久影院| 久久精品国产综合久久久| 女的被弄到高潮叫床怎么办| 成人免费观看视频高清| 亚洲综合色惰| 99久久精品国产国产毛片| 亚洲国产精品国产精品| 亚洲欧美精品自产自拍| 日本猛色少妇xxxxx猛交久久| 又黄又粗又硬又大视频| 亚洲精品一区蜜桃| 国产成人精品久久久久久| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产欧美网| 伊人久久大香线蕉亚洲五| 九草在线视频观看| 99久久中文字幕三级久久日本| 亚洲国产欧美网| 老熟女久久久| 成年美女黄网站色视频大全免费| 亚洲内射少妇av| 五月天丁香电影| 午夜影院在线不卡| 满18在线观看网站| 一级片'在线观看视频| 国产精品嫩草影院av在线观看| 一个人免费看片子| a级毛片在线看网站| 高清av免费在线| 国产老妇伦熟女老妇高清| 午夜影院在线不卡| 亚洲av国产av综合av卡| 精品第一国产精品| 九九爱精品视频在线观看| 成人免费观看视频高清| 久久久久国产一级毛片高清牌| 午夜福利影视在线免费观看| 亚洲第一av免费看| 亚洲av日韩在线播放| 亚洲av福利一区| 女性被躁到高潮视频| 日本猛色少妇xxxxx猛交久久| a 毛片基地| 久久这里有精品视频免费| 男女下面插进去视频免费观看| 久久久久视频综合| 一二三四在线观看免费中文在| 少妇的逼水好多| 麻豆av在线久日| 中文字幕最新亚洲高清| 国产精品蜜桃在线观看| 99精国产麻豆久久婷婷| 日本猛色少妇xxxxx猛交久久| 国产精品99久久99久久久不卡 | 精品福利永久在线观看| 国产激情久久老熟女| 欧美日韩精品网址| 精品少妇黑人巨大在线播放| 国产精品秋霞免费鲁丝片| av卡一久久| 亚洲图色成人| 精品国产乱码久久久久久男人| 在线观看免费日韩欧美大片| 欧美日韩精品成人综合77777| 制服诱惑二区| 日本vs欧美在线观看视频| 桃花免费在线播放| 国产人伦9x9x在线观看 | 黄色一级大片看看| 伊人亚洲综合成人网| 精品国产一区二区久久| 精品一区在线观看国产| 久久精品国产综合久久久| 日韩中文字幕欧美一区二区 | 人人澡人人妻人| 99久国产av精品国产电影| 精品99又大又爽又粗少妇毛片| 国产一区二区三区综合在线观看| 成人国语在线视频| 久久人人爽人人片av| 国产毛片在线视频| 一级毛片 在线播放| 春色校园在线视频观看| 亚洲精华国产精华液的使用体验| 日韩不卡一区二区三区视频在线| 在线 av 中文字幕| 中文字幕人妻丝袜一区二区 | 国产1区2区3区精品| 亚洲图色成人| 日韩av在线免费看完整版不卡| 国产一区二区 视频在线| 亚洲国产欧美日韩在线播放| 久久国产精品大桥未久av| 18禁裸乳无遮挡动漫免费视频| 男人爽女人下面视频在线观看| 另类亚洲欧美激情| 欧美bdsm另类| 曰老女人黄片| 欧美+日韩+精品| 26uuu在线亚洲综合色| 午夜老司机福利剧场| 欧美在线黄色| 亚洲天堂av无毛| 91精品伊人久久大香线蕉| 天天躁狠狠躁夜夜躁狠狠躁| 最新的欧美精品一区二区| 久久久久久久精品精品| 亚洲第一av免费看| 欧美日韩精品成人综合77777| 亚洲国产最新在线播放| 宅男免费午夜| 日日撸夜夜添| 天堂俺去俺来也www色官网| 国产在线免费精品| 十分钟在线观看高清视频www| 久久久久精品性色| 精品酒店卫生间| 亚洲综合色惰| 99热全是精品| 国产av码专区亚洲av| 有码 亚洲区| 超碰97精品在线观看| 成人亚洲精品一区在线观看| 国产视频首页在线观看| 中文字幕制服av| 你懂的网址亚洲精品在线观看| 久久久精品94久久精品| 日日啪夜夜爽| 91国产中文字幕| 好男人视频免费观看在线| 日本欧美视频一区| av福利片在线| 亚洲欧美一区二区三区国产| 午夜福利在线免费观看网站| 多毛熟女@视频| 成人二区视频| 伊人亚洲综合成人网| xxxhd国产人妻xxx| h视频一区二区三区| 中文欧美无线码| 亚洲婷婷狠狠爱综合网| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜激情av网站| 在线精品无人区一区二区三| 成人国语在线视频| 女的被弄到高潮叫床怎么办| 亚洲激情五月婷婷啪啪| 春色校园在线视频观看| 看免费成人av毛片| 多毛熟女@视频| 国产欧美日韩一区二区三区在线| 日韩人妻精品一区2区三区| 综合色丁香网| 亚洲精品久久成人aⅴ小说| 日本av免费视频播放| 一级毛片我不卡| 国产午夜精品一二区理论片| 精品国产一区二区久久| 在线观看www视频免费| 亚洲人成77777在线视频| 高清av免费在线| 黄色 视频免费看| 亚洲第一区二区三区不卡| 中文字幕最新亚洲高清| 亚洲第一av免费看| av免费观看日本| www日本在线高清视频| 国产毛片在线视频| 有码 亚洲区| 久久久久久久国产电影| 国产熟女欧美一区二区| 青春草亚洲视频在线观看| av线在线观看网站| 国产黄频视频在线观看| 丝袜脚勾引网站| 国产欧美亚洲国产| 精品一区在线观看国产| 在线观看一区二区三区激情| 国产极品粉嫩免费观看在线| 日韩中文字幕欧美一区二区 | 老熟女久久久| 欧美日韩成人在线一区二区| 欧美国产精品va在线观看不卡| 欧美黄色片欧美黄色片| www.精华液| 熟妇人妻不卡中文字幕| 黄片无遮挡物在线观看| 亚洲一区中文字幕在线| 国产av精品麻豆| 久久影院123| 国产日韩欧美亚洲二区| 日韩伦理黄色片| 亚洲欧美一区二区三区国产| 午夜福利在线观看免费完整高清在| 国产综合精华液| 在线观看一区二区三区激情| 极品人妻少妇av视频| 国产日韩欧美视频二区| 男女免费视频国产| 午夜福利视频精品| 亚洲久久久国产精品| 日日撸夜夜添| 久久人人97超碰香蕉20202| 男女啪啪激烈高潮av片| 精品久久久精品久久久| 伦精品一区二区三区| 国产一区二区在线观看av| 99热全是精品| 亚洲欧美精品自产自拍| 日韩制服骚丝袜av| 高清av免费在线| 国产精品秋霞免费鲁丝片| 亚洲成国产人片在线观看| 午夜激情久久久久久久| 在线观看国产h片| 亚洲情色 制服丝袜| 91在线精品国自产拍蜜月| 久久精品国产综合久久久| 久久午夜福利片| 九色亚洲精品在线播放| 啦啦啦啦在线视频资源| 黑人巨大精品欧美一区二区蜜桃| 天天躁夜夜躁狠狠躁躁| 日本欧美视频一区| 欧美精品人与动牲交sv欧美| av一本久久久久| 在线免费观看不下载黄p国产| 精品人妻在线不人妻| 日韩人妻精品一区2区三区| 国产激情久久老熟女| 黄色视频在线播放观看不卡| 欧美黄色片欧美黄色片| 亚洲美女视频黄频| 久久国内精品自在自线图片| 精品人妻在线不人妻| 久久久亚洲精品成人影院| 亚洲欧洲日产国产| 精品国产超薄肉色丝袜足j| 亚洲,欧美,日韩| 成人手机av| av在线播放精品| 考比视频在线观看| 日韩av免费高清视频| 亚洲伊人色综图| 街头女战士在线观看网站| www日本在线高清视频| 亚洲中文av在线| 国产高清国产精品国产三级| h视频一区二区三区| 亚洲av.av天堂| 晚上一个人看的免费电影| 老鸭窝网址在线观看| 久久午夜福利片| 99国产精品免费福利视频| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利在线免费观看网站| 久久 成人 亚洲| 久久亚洲国产成人精品v| av又黄又爽大尺度在线免费看| 中文字幕亚洲精品专区| 亚洲熟女精品中文字幕| 国产精品国产三级国产专区5o| 久久人妻熟女aⅴ| 国产成人免费无遮挡视频| 欧美激情 高清一区二区三区| 大片电影免费在线观看免费| 巨乳人妻的诱惑在线观看| 国产麻豆69| 亚洲精品在线美女| 久久久久久久久久久久大奶| 免费日韩欧美在线观看| 午夜福利在线观看免费完整高清在| 国产成人a∨麻豆精品| 人妻一区二区av| 午夜久久久在线观看| 在线观看人妻少妇| 女人精品久久久久毛片| 精品人妻一区二区三区麻豆| 一区二区av电影网| 三上悠亚av全集在线观看| 国产黄频视频在线观看| 黄网站色视频无遮挡免费观看| 国产福利在线免费观看视频| 久久国内精品自在自线图片| 精品国产乱码久久久久久男人| 国产1区2区3区精品| 欧美+日韩+精品| 国产视频首页在线观看| 亚洲视频免费观看视频| 搡老乐熟女国产| h视频一区二区三区| 中文字幕色久视频| 亚洲欧美色中文字幕在线| 中文字幕人妻熟女乱码| 黄色 视频免费看| 日产精品乱码卡一卡2卡三| 欧美日韩av久久| 天堂中文最新版在线下载| 午夜91福利影院| 久久久久国产网址| 久久精品久久精品一区二区三区| 欧美日本中文国产一区发布| 亚洲 欧美一区二区三区| 亚洲第一青青草原| 精品第一国产精品| a级片在线免费高清观看视频| 高清av免费在线| 在线观看一区二区三区激情| 精品久久蜜臀av无| 久久久久久久亚洲中文字幕| 国产精品亚洲av一区麻豆 | 老司机影院毛片| 啦啦啦中文免费视频观看日本| 久久av网站| 男女国产视频网站| 另类亚洲欧美激情| 青春草亚洲视频在线观看| 国产午夜精品一二区理论片| av天堂久久9| 少妇被粗大的猛进出69影院| 中文字幕最新亚洲高清| 国产伦理片在线播放av一区| 欧美日韩视频高清一区二区三区二| 国产不卡av网站在线观看| 国产精品欧美亚洲77777| 日韩视频在线欧美| 亚洲内射少妇av| 亚洲国产看品久久| 国产一区亚洲一区在线观看| 国产成人精品一,二区| 亚洲成国产人片在线观看| 国产日韩欧美亚洲二区| 少妇的丰满在线观看| 精品亚洲成国产av| 亚洲国产欧美在线一区| 男人添女人高潮全过程视频| 国产男女内射视频| 久热这里只有精品99| 欧美中文综合在线视频| 少妇被粗大的猛进出69影院| 久久影院123| a级片在线免费高清观看视频| 成人影院久久| 婷婷成人精品国产| 五月伊人婷婷丁香| 亚洲精品国产av蜜桃| 久久精品国产自在天天线| 欧美日韩亚洲高清精品| 国产精品无大码| 欧美人与性动交α欧美软件| 黄色一级大片看看| 中国国产av一级| 精品人妻一区二区三区麻豆| 波多野结衣一区麻豆| 亚洲精品日本国产第一区| 欧美日本中文国产一区发布| 激情五月婷婷亚洲| 国产成人精品一,二区| 男人操女人黄网站| 97精品久久久久久久久久精品| 国产免费现黄频在线看| 国产视频首页在线观看| 国产xxxxx性猛交| 欧美日韩精品成人综合77777| 9热在线视频观看99| 女性生殖器流出的白浆| 亚洲情色 制服丝袜| 日韩视频在线欧美| 1024视频免费在线观看| 国产精品久久久久成人av| 捣出白浆h1v1| 成年女人毛片免费观看观看9 | 午夜福利乱码中文字幕| 国产成人精品无人区| 国精品久久久久久国模美| 久久婷婷青草| 五月伊人婷婷丁香| 国产成人aa在线观看| 国产男女超爽视频在线观看| 色哟哟·www| 精品亚洲成a人片在线观看| 亚洲美女黄色视频免费看| 亚洲四区av| videosex国产| 精品第一国产精品| 在线观看免费视频网站a站| 日韩精品有码人妻一区| www日本在线高清视频| 久久久久久久久久久免费av| 边亲边吃奶的免费视频| 麻豆乱淫一区二区| 日本午夜av视频| 成人二区视频| 少妇人妻久久综合中文| 青春草视频在线免费观看| 欧美少妇被猛烈插入视频| 亚洲精华国产精华液的使用体验| av在线播放精品| 亚洲av电影在线进入| 又黄又粗又硬又大视频| 中文天堂在线官网| 国产精品一二三区在线看| 美女主播在线视频| 午夜福利视频在线观看免费| 97在线视频观看| 深夜精品福利| 好男人视频免费观看在线| 成年美女黄网站色视频大全免费| 久久久久国产精品人妻一区二区| 亚洲熟女精品中文字幕| 国产日韩欧美亚洲二区| √禁漫天堂资源中文www| 日韩在线高清观看一区二区三区| 国精品久久久久久国模美| 黑丝袜美女国产一区| 老司机影院成人| 18在线观看网站| 18禁动态无遮挡网站| 亚洲国产av新网站| 国产xxxxx性猛交| 久久ye,这里只有精品| 欧美+日韩+精品| 亚洲成人手机| 成人毛片60女人毛片免费| 免费黄色在线免费观看| 少妇被粗大猛烈的视频| 精品亚洲成a人片在线观看| 亚洲内射少妇av| 国产伦理片在线播放av一区| 亚洲婷婷狠狠爱综合网| 精品福利永久在线观看| 久久久久久免费高清国产稀缺| 亚洲国产最新在线播放| 国产福利在线免费观看视频| 啦啦啦啦在线视频资源| av片东京热男人的天堂| 秋霞伦理黄片| av国产精品久久久久影院| 叶爱在线成人免费视频播放| 水蜜桃什么品种好| av福利片在线| 欧美在线黄色| 黄片播放在线免费| 99国产综合亚洲精品| 亚洲国产精品成人久久小说| 久久99蜜桃精品久久| 成年av动漫网址| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品一二三| 啦啦啦视频在线资源免费观看| 丝袜美腿诱惑在线| 久久精品夜色国产| 少妇人妻 视频| 精品国产一区二区三区四区第35| 亚洲情色 制服丝袜| 亚洲人成77777在线视频| 曰老女人黄片| av免费观看日本| 波多野结衣av一区二区av| 国产亚洲最大av| 亚洲中文av在线| 女性被躁到高潮视频| 黄片无遮挡物在线观看| 日本91视频免费播放| 国产精品免费大片| 看免费成人av毛片| 男人爽女人下面视频在线观看| 国产精品一国产av| 亚洲少妇的诱惑av| 国产高清国产精品国产三级| 夫妻性生交免费视频一级片| 97精品久久久久久久久久精品| 99热国产这里只有精品6| 国产xxxxx性猛交| 综合色丁香网| 日本-黄色视频高清免费观看| 国产日韩欧美视频二区| 性少妇av在线| 欧美变态另类bdsm刘玥| 嫩草影院入口| 电影成人av| 精品99又大又爽又粗少妇毛片| 国产亚洲欧美精品永久| 久久久久国产网址| 美女国产高潮福利片在线看| 成人黄色视频免费在线看| 亚洲欧美色中文字幕在线| 亚洲av中文av极速乱| 三级国产精品片| 亚洲熟女精品中文字幕| 成人二区视频| 如日韩欧美国产精品一区二区三区| 美女国产高潮福利片在线看| 免费观看在线日韩| 夜夜骑夜夜射夜夜干| 日韩熟女老妇一区二区性免费视频| 日韩成人av中文字幕在线观看| 高清视频免费观看一区二区| 亚洲精品久久午夜乱码| 欧美日本中文国产一区发布| 久久97久久精品| 国产免费现黄频在线看| 9191精品国产免费久久| 人体艺术视频欧美日本| 高清欧美精品videossex| 国产免费福利视频在线观看| www日本在线高清视频| 欧美精品一区二区大全| 欧美日韩精品网址| 国产黄色免费在线视频| 黄网站色视频无遮挡免费观看| 国产黄频视频在线观看| 久热久热在线精品观看| 99热国产这里只有精品6| 老熟女久久久| 中国三级夫妇交换| 美女中出高潮动态图| 夫妻午夜视频| 午夜免费男女啪啪视频观看| 国产一区二区 视频在线| 欧美人与善性xxx| 亚洲情色 制服丝袜| 精品人妻熟女毛片av久久网站| 成人漫画全彩无遮挡| 亚洲欧美成人精品一区二区| 久久午夜综合久久蜜桃| 欧美黄色片欧美黄色片| 91在线精品国自产拍蜜月| 国产免费现黄频在线看| 一区二区三区乱码不卡18| 永久网站在线| 久久久久国产精品人妻一区二区| 免费观看a级毛片全部| 免费黄网站久久成人精品| 美女国产视频在线观看| tube8黄色片| 国产精品欧美亚洲77777| 尾随美女入室| 久久精品久久精品一区二区三区| 久久女婷五月综合色啪小说| 婷婷色综合www| 男女边摸边吃奶| 妹子高潮喷水视频| 日韩制服骚丝袜av| 国产xxxxx性猛交| 成人漫画全彩无遮挡| 欧美激情极品国产一区二区三区| 久久 成人 亚洲| 肉色欧美久久久久久久蜜桃| 亚洲欧美中文字幕日韩二区| 亚洲精品日韩在线中文字幕| 国产亚洲精品第一综合不卡| 久久综合国产亚洲精品| 边亲边吃奶的免费视频| 久久精品国产鲁丝片午夜精品| 国产精品女同一区二区软件| 肉色欧美久久久久久久蜜桃| 久久久久久人人人人人| 男女边吃奶边做爰视频| 国产精品av久久久久免费| 亚洲精品在线美女| 亚洲内射少妇av| 欧美日韩综合久久久久久| 国产一级毛片在线| 欧美少妇被猛烈插入视频| 国产熟女欧美一区二区| 18禁动态无遮挡网站| 国产视频首页在线观看| 永久网站在线| 免费黄网站久久成人精品| 成人漫画全彩无遮挡| 午夜福利网站1000一区二区三区| 久久久久久人妻| 9色porny在线观看| 亚洲国产精品一区二区三区在线| 黑人巨大精品欧美一区二区蜜桃| 久久久久久久久免费视频了| 亚洲av欧美aⅴ国产| 国产亚洲精品第一综合不卡|