• <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ù)概述
    757午夜福利合集在线观看| 天堂av国产一区二区熟女人妻| 如何舔出高潮| 色av中文字幕| 日韩欧美国产在线观看| 国产 一区 欧美 日韩| 国内久久婷婷六月综合欲色啪| 久久久久精品国产欧美久久久| 国产久久久一区二区三区| 成人亚洲精品av一区二区| 亚洲美女搞黄在线观看 | 国产精品av视频在线免费观看| 午夜a级毛片| 深夜精品福利| 国产一区二区在线av高清观看| 国产黄a三级三级三级人| 欧美性猛交╳xxx乱大交人| 国产精品日韩av在线免费观看| 久久久久久久久大av| 国产91精品成人一区二区三区| 精品午夜福利在线看| 色综合站精品国产| 欧美性感艳星| 一区二区三区免费毛片| 91在线精品国自产拍蜜月| 露出奶头的视频| 日日干狠狠操夜夜爽| 日韩欧美精品免费久久 | 国产乱人伦免费视频| 在线天堂最新版资源| 国产男靠女视频免费网站| 色在线成人网| avwww免费| 精品乱码久久久久久99久播| 90打野战视频偷拍视频| 日本 欧美在线| 精品熟女少妇八av免费久了| 在线免费观看不下载黄p国产 | 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区三区av网在线观看| 小蜜桃在线观看免费完整版高清| 看黄色毛片网站| 激情在线观看视频在线高清| 久久香蕉精品热| 久久精品人妻少妇| 国产不卡一卡二| 丰满人妻一区二区三区视频av| 我要看日韩黄色一级片| 日韩欧美精品免费久久 | 自拍偷自拍亚洲精品老妇| 在线观看一区二区三区| 免费av毛片视频| 国产精品自产拍在线观看55亚洲| 99久久无色码亚洲精品果冻| 国产亚洲av嫩草精品影院| 一进一出抽搐动态| bbb黄色大片| 色综合欧美亚洲国产小说| 精品久久久久久久久久免费视频| 午夜日韩欧美国产| 国产精品国产高清国产av| 1000部很黄的大片| 少妇人妻一区二区三区视频| 亚洲av.av天堂| 九色成人免费人妻av| 岛国在线免费视频观看| 麻豆国产av国片精品| 久久久国产成人精品二区| 99久久成人亚洲精品观看| 免费在线观看日本一区| 国产伦在线观看视频一区| 在线观看66精品国产| 久久久久久九九精品二区国产| 成人国产一区最新在线观看| av天堂在线播放| 亚洲综合色惰| 欧美性猛交黑人性爽| 啪啪无遮挡十八禁网站| 日韩亚洲欧美综合| 首页视频小说图片口味搜索| 国产aⅴ精品一区二区三区波| 免费观看精品视频网站| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久久久久电影| 国产野战对白在线观看| 黄色女人牲交| 国产成人福利小说| 日本免费a在线| 久久久久九九精品影院| 一进一出抽搐动态| 日本三级黄在线观看| 看免费av毛片| 少妇高潮的动态图| 永久网站在线| 丁香欧美五月| 日韩欧美三级三区| 偷拍熟女少妇极品色| 午夜老司机福利剧场| 一进一出抽搐动态| 亚洲av不卡在线观看| 亚洲人成网站高清观看| 国产亚洲精品av在线| 看十八女毛片水多多多| 亚洲美女黄片视频| 日本精品一区二区三区蜜桃| 一级毛片久久久久久久久女| 美女cb高潮喷水在线观看| 国产成人欧美在线观看| 成人美女网站在线观看视频| 波多野结衣巨乳人妻| 一个人看视频在线观看www免费| 九九在线视频观看精品| av中文乱码字幕在线| 在线免费观看的www视频| 色播亚洲综合网| 三级毛片av免费| 直男gayav资源| 婷婷亚洲欧美| 日韩免费av在线播放| 在线观看一区二区三区| av在线观看视频网站免费| 久久亚洲精品不卡| 久久热精品热| 久久久成人免费电影| 久久久久久大精品| 午夜老司机福利剧场| 亚洲av电影不卡..在线观看| 精品久久久久久久人妻蜜臀av| av天堂中文字幕网| 日韩欧美精品v在线| 欧美一区二区国产精品久久精品| 给我免费播放毛片高清在线观看| 久久精品国产99精品国产亚洲性色| 亚洲第一电影网av| 床上黄色一级片| 男人的好看免费观看在线视频| 亚洲精品456在线播放app | 欧美日韩国产亚洲二区| 欧美+亚洲+日韩+国产| 精品午夜福利在线看| 一区二区三区四区激情视频 | 美女被艹到高潮喷水动态| h日本视频在线播放| 久久久久亚洲av毛片大全| 国产精品一及| 91在线精品国自产拍蜜月| 免费在线观看亚洲国产| 看免费av毛片| 一个人看的www免费观看视频| 午夜福利18| 色综合站精品国产| 性色avwww在线观看| 女人被狂操c到高潮| 国产精品人妻久久久久久| 精品一区二区三区人妻视频| 国产蜜桃级精品一区二区三区| 观看美女的网站| 国产成人av教育| 国产高清三级在线| 欧美绝顶高潮抽搐喷水| 日韩人妻高清精品专区| 他把我摸到了高潮在线观看| 日韩欧美免费精品| 国产真实乱freesex| 亚洲第一电影网av| 亚洲欧美日韩东京热| 日韩人妻高清精品专区| 黄片小视频在线播放| 嫩草影院入口| 精品久久久久久成人av| 看片在线看免费视频| 国产在线男女| 一个人看视频在线观看www免费| 日韩免费av在线播放| 久久99热这里只有精品18| 国产精品影院久久| 亚洲精品色激情综合| 国产91精品成人一区二区三区| 国产三级黄色录像| 俺也久久电影网| 成人亚洲精品av一区二区| 日日摸夜夜添夜夜添av毛片 | 午夜福利在线在线| 国产成人福利小说| 国产真实乱freesex| 国产色婷婷99| 日本成人三级电影网站| 国内揄拍国产精品人妻在线| 俄罗斯特黄特色一大片| 午夜精品久久久久久毛片777| 真实男女啪啪啪动态图| 欧美最新免费一区二区三区 | 天天躁日日操中文字幕| 免费黄网站久久成人精品 | 岛国在线免费视频观看| 99热6这里只有精品| 看免费av毛片| 女同久久另类99精品国产91| 日韩欧美精品v在线| 午夜亚洲福利在线播放| 午夜激情福利司机影院| 欧美最黄视频在线播放免费| 变态另类成人亚洲欧美熟女| 少妇被粗大猛烈的视频| 性色av乱码一区二区三区2| 国产成人aa在线观看| 高潮久久久久久久久久久不卡| 变态另类成人亚洲欧美熟女| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 午夜精品在线福利| 丰满人妻熟妇乱又伦精品不卡| 黄色一级大片看看| 婷婷色综合大香蕉| 国产黄色小视频在线观看| 国产单亲对白刺激| 国产精品一区二区三区四区免费观看 | 久久久久久久久久黄片| 丝袜美腿在线中文| 性欧美人与动物交配| 看黄色毛片网站| 亚洲精品亚洲一区二区| 久久精品综合一区二区三区| 99在线视频只有这里精品首页| 九色国产91popny在线| 在线观看免费视频日本深夜| 国产主播在线观看一区二区| 亚洲黑人精品在线| 久久精品国产99精品国产亚洲性色| 9191精品国产免费久久| 又粗又爽又猛毛片免费看| 波多野结衣高清作品| 性插视频无遮挡在线免费观看| 国产成人啪精品午夜网站| 亚洲最大成人中文| 精品久久久久久久久亚洲 | 精品久久久久久久末码| 亚洲第一电影网av| 亚洲经典国产精华液单 | 国产亚洲欧美98| 小蜜桃在线观看免费完整版高清| 成人av一区二区三区在线看| 国产91精品成人一区二区三区| 欧美黑人巨大hd| 99久久成人亚洲精品观看| 午夜两性在线视频| 久久久久久久久大av| 欧美一级a爱片免费观看看| 少妇高潮的动态图| 一a级毛片在线观看| 两个人视频免费观看高清| 精品日产1卡2卡| 成人国产综合亚洲| 日本 av在线| 男女视频在线观看网站免费| 欧美国产日韩亚洲一区| 中文亚洲av片在线观看爽| 久久久久久久久大av| 18禁裸乳无遮挡免费网站照片| 九色国产91popny在线| 国产一级毛片七仙女欲春2| 久久精品影院6| 欧美xxxx性猛交bbbb| 99精品在免费线老司机午夜| 国产精品一区二区性色av| 波多野结衣高清作品| 亚洲国产色片| 日韩欧美在线乱码| 精品人妻1区二区| 极品教师在线免费播放| 亚洲一区高清亚洲精品| 欧美乱色亚洲激情| 午夜影院日韩av| 嫩草影院入口| 一本精品99久久精品77| 精品国产亚洲在线| 亚洲av免费在线观看| 久久久色成人| 亚洲欧美日韩高清在线视频| 亚洲 欧美 日韩 在线 免费| 一二三四社区在线视频社区8| 亚洲人成电影免费在线| 日韩大尺度精品在线看网址| 每晚都被弄得嗷嗷叫到高潮| 国产主播在线观看一区二区| АⅤ资源中文在线天堂| 成人av在线播放网站| 赤兔流量卡办理| 国产男靠女视频免费网站| 美女xxoo啪啪120秒动态图 | 村上凉子中文字幕在线| 男女之事视频高清在线观看| 欧美另类亚洲清纯唯美| 又紧又爽又黄一区二区| 亚洲熟妇中文字幕五十中出| 精品欧美国产一区二区三| 啦啦啦观看免费观看视频高清| 男人舔奶头视频| 天美传媒精品一区二区| 久久伊人香网站| 欧美区成人在线视频| 久久久久久久久久黄片| 精品人妻熟女av久视频| 99久久九九国产精品国产免费| av欧美777| 真人一进一出gif抽搐免费| 中文资源天堂在线| 无遮挡黄片免费观看| 亚洲国产精品成人综合色| 天堂影院成人在线观看| 午夜福利高清视频| a级毛片免费高清观看在线播放| 99热这里只有是精品在线观看 | 成人午夜高清在线视频| 国产在视频线在精品| 国内精品美女久久久久久| 国产精品自产拍在线观看55亚洲| 国产v大片淫在线免费观看| 中文亚洲av片在线观看爽| 国产一区二区亚洲精品在线观看| 97超级碰碰碰精品色视频在线观看| 欧美+日韩+精品| 久久精品国产亚洲av香蕉五月| 欧美成狂野欧美在线观看| 日日摸夜夜添夜夜添av毛片 | 天堂动漫精品| 99热6这里只有精品| 国产伦在线观看视频一区| 在线免费观看的www视频| 非洲黑人性xxxx精品又粗又长| 国产黄片美女视频| www.999成人在线观看| 老司机午夜十八禁免费视频| 亚洲,欧美,日韩| 男女做爰动态图高潮gif福利片| 亚洲人成电影免费在线| 亚洲欧美日韩东京热| 亚洲,欧美精品.| 很黄的视频免费| 亚洲色图av天堂| 精品久久国产蜜桃| 欧美日韩亚洲国产一区二区在线观看| 夜夜看夜夜爽夜夜摸| 国产一区二区在线av高清观看| 一级av片app| 国产视频一区二区在线看| 久久久久久大精品| 日本三级黄在线观看| 国产精品一区二区性色av| 乱人视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 我的老师免费观看完整版| 国内毛片毛片毛片毛片毛片| 久久草成人影院| 一区二区三区激情视频| 欧美性猛交黑人性爽| 亚洲av不卡在线观看| 亚洲第一区二区三区不卡| 国产黄色小视频在线观看| 高潮久久久久久久久久久不卡| 亚洲国产精品999在线| 午夜免费激情av| 亚洲午夜理论影院| 亚洲精品成人久久久久久| 亚洲欧美日韩高清专用| 好看av亚洲va欧美ⅴa在| 久久婷婷人人爽人人干人人爱| 精品欧美国产一区二区三| 国产精品乱码一区二三区的特点| 天堂动漫精品| 精品午夜福利视频在线观看一区| 国产一级毛片七仙女欲春2| 国产精品女同一区二区软件 | 免费观看精品视频网站| 日韩成人在线观看一区二区三区| 亚洲最大成人手机在线| 国产伦在线观看视频一区| 国产91精品成人一区二区三区| 99久久精品国产亚洲精品| 精品不卡国产一区二区三区| 综合色av麻豆| 99久久成人亚洲精品观看| 婷婷精品国产亚洲av在线| www.www免费av| 91狼人影院| 男女床上黄色一级片免费看| 欧美丝袜亚洲另类 | 国内精品美女久久久久久| 舔av片在线| 国产精品一区二区三区四区免费观看 | 午夜精品久久久久久毛片777| 看十八女毛片水多多多| 国产成人福利小说| 精品99又大又爽又粗少妇毛片 | 自拍偷自拍亚洲精品老妇| 美女cb高潮喷水在线观看| 一本综合久久免费| 亚洲av美国av| 99热这里只有精品一区| 国产三级中文精品| 亚洲精华国产精华精| 91九色精品人成在线观看| 99热只有精品国产| 中出人妻视频一区二区| 国产爱豆传媒在线观看| 欧美黑人巨大hd| 国产午夜精品论理片| www.999成人在线观看| 亚洲久久久久久中文字幕| 国产三级中文精品| 性色avwww在线观看| 中文资源天堂在线| 亚洲av美国av| 成熟少妇高潮喷水视频| 好男人电影高清在线观看| 久久久久国内视频| 亚洲,欧美,日韩| АⅤ资源中文在线天堂| 国产成人啪精品午夜网站| 亚洲乱码一区二区免费版| 久99久视频精品免费| 国内精品一区二区在线观看| 国产亚洲欧美在线一区二区| 99热这里只有精品一区| 99国产极品粉嫩在线观看| 亚洲成av人片免费观看| 在线十欧美十亚洲十日本专区| 九九热线精品视视频播放| 国产综合懂色| 少妇丰满av| 97超级碰碰碰精品色视频在线观看| 色哟哟哟哟哟哟| 中亚洲国语对白在线视频| 内射极品少妇av片p| 精品久久国产蜜桃| 国产在视频线在精品| 日韩中字成人| 国内精品久久久久久久电影| 成年女人看的毛片在线观看| 51午夜福利影视在线观看| 国产精品人妻久久久久久| 男女床上黄色一级片免费看| 99热这里只有精品一区| 国产午夜精品久久久久久一区二区三区 | 最新在线观看一区二区三区| 赤兔流量卡办理| 国产精品av视频在线免费观看| x7x7x7水蜜桃| 国产极品精品免费视频能看的| 一区二区三区免费毛片| 欧美激情久久久久久爽电影| 国产精品爽爽va在线观看网站| 欧美最黄视频在线播放免费| 宅男免费午夜| 内地一区二区视频在线| 亚洲精品日韩av片在线观看| 51国产日韩欧美| 精品人妻偷拍中文字幕| 午夜福利视频1000在线观看| 欧美日韩综合久久久久久 | 国产国拍精品亚洲av在线观看| 国产又黄又爽又无遮挡在线| 在线观看一区二区三区| 夜夜夜夜夜久久久久| 成人毛片a级毛片在线播放| 欧美乱妇无乱码| 中文亚洲av片在线观看爽| 国产精品久久久久久人妻精品电影| 久久精品国产亚洲av涩爱 | 一级作爱视频免费观看| 免费av不卡在线播放| 久久久久久久精品吃奶| 国产精品不卡视频一区二区 | 啦啦啦观看免费观看视频高清| 日韩欧美一区二区三区在线观看| 精品日产1卡2卡| 丁香欧美五月| 免费看日本二区| 亚洲乱码一区二区免费版| 亚洲人与动物交配视频| 熟女电影av网| 日韩欧美精品免费久久 | 国模一区二区三区四区视频| 看片在线看免费视频| 亚洲熟妇熟女久久| 国产又黄又爽又无遮挡在线| 国产免费av片在线观看野外av| 欧美极品一区二区三区四区| 2021天堂中文幕一二区在线观| 免费在线观看亚洲国产| 国产精品综合久久久久久久免费| 毛片女人毛片| 亚洲成人免费电影在线观看| 中亚洲国语对白在线视频| 国产精品1区2区在线观看.| 色播亚洲综合网| 黄色一级大片看看| 两个人视频免费观看高清| 国产探花极品一区二区| 免费观看人在逋| 看黄色毛片网站| 久久久久久国产a免费观看| 免费人成在线观看视频色| 国产av不卡久久| 国产av麻豆久久久久久久| 天堂动漫精品| 夜夜躁狠狠躁天天躁| 丰满人妻一区二区三区视频av| 欧美国产日韩亚洲一区| 精品国产三级普通话版| 精品99又大又爽又粗少妇毛片 | 久久久久久久久久黄片| 国内久久婷婷六月综合欲色啪| avwww免费| 日本在线视频免费播放| 搡女人真爽免费视频火全软件 | 国产爱豆传媒在线观看| 亚洲国产色片| 日韩免费av在线播放| 精品熟女少妇八av免费久了| 日本免费一区二区三区高清不卡| 国产人妻一区二区三区在| 欧美一区二区精品小视频在线| 欧美色视频一区免费| 两个人的视频大全免费| 国产精品嫩草影院av在线观看 | 老司机深夜福利视频在线观看| 偷拍熟女少妇极品色| 18+在线观看网站| 乱人视频在线观看| 欧美日韩国产亚洲二区| 欧美日韩乱码在线| 亚洲最大成人中文| а√天堂www在线а√下载| 熟妇人妻久久中文字幕3abv| 国产伦人伦偷精品视频| 特级一级黄色大片| 国产伦人伦偷精品视频| h日本视频在线播放| 夜夜爽天天搞| 国产视频一区二区在线看| 国产高清激情床上av| 99热这里只有精品一区| www日本黄色视频网| 精品久久久久久成人av| 一进一出好大好爽视频| 乱码一卡2卡4卡精品| 国产精品一及| av在线老鸭窝| 一二三四社区在线视频社区8| 精品国内亚洲2022精品成人| 精品久久久久久久久av| 午夜免费男女啪啪视频观看 | 欧美最新免费一区二区三区 | 99热这里只有是精品50| 亚洲熟妇中文字幕五十中出| 可以在线观看毛片的网站| 精品久久国产蜜桃| 18+在线观看网站| 99国产综合亚洲精品| 757午夜福利合集在线观看| 99久久无色码亚洲精品果冻| 国产免费男女视频| 久久久久久大精品| 国产成人a区在线观看| 国产精品影院久久| 成人午夜高清在线视频| 精品一区二区三区av网在线观看| 久久精品人妻少妇| 亚洲七黄色美女视频| 99热这里只有精品一区| 91久久精品电影网| 一级毛片久久久久久久久女| 免费看日本二区| 中文亚洲av片在线观看爽| 一本一本综合久久| 神马国产精品三级电影在线观看| 一级黄片播放器| 2021天堂中文幕一二区在线观| 精品午夜福利视频在线观看一区| 99久久精品热视频| 久久久久久久亚洲中文字幕 | 国产av麻豆久久久久久久| 亚洲最大成人av| 亚洲国产精品合色在线| 日韩高清综合在线| 日韩欧美 国产精品| 麻豆av噜噜一区二区三区| av欧美777| 精品一区二区三区av网在线观看| 成人特级av手机在线观看| 亚洲,欧美,日韩| 久久亚洲真实| 国产极品精品免费视频能看的| 99精品久久久久人妻精品| 精品人妻1区二区| 亚洲av不卡在线观看| 黄片小视频在线播放| 国产在线男女| 亚洲av第一区精品v没综合| 99视频精品全部免费 在线| 亚洲欧美日韩高清专用| 精品人妻1区二区| 蜜桃久久精品国产亚洲av| 国产91精品成人一区二区三区| 国产欧美日韩精品一区二区| 黄色视频,在线免费观看| 精品福利观看| 1024手机看黄色片| 日韩欧美精品v在线| 天堂影院成人在线观看| 国产三级黄色录像| 国产精品乱码一区二三区的特点| 18禁裸乳无遮挡免费网站照片| 91字幕亚洲|