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

    馬里亞納弧前蛇紋巖泥火山無機(jī)成因甲烷–氫氣水合物形成條件與穩(wěn)定帶發(fā)育特征

    2022-04-28 00:50:26湯加麗曹運(yùn)誠陳多福
    地球化學(xué) 2022年2期
    關(guān)鍵詞:泥火山馬里亞納水合物

    湯加麗, 曹運(yùn)誠, 陳多福

    馬里亞納弧前蛇紋巖泥火山無機(jī)成因甲烷–氫氣水合物形成條件與穩(wěn)定帶發(fā)育特征

    湯加麗, 曹運(yùn)誠*, 陳多福

    (上海海洋大學(xué) 海洋科學(xué)學(xué)院 上海深淵科學(xué)工程技術(shù)研究中心, 上海 201306)

    海洋水合物的成藏不僅需要合適的溫度、壓力條件, 而且需要充足的氣源。一般認(rèn)為大洋區(qū)海底沉積物缺乏豐富的有機(jī)質(zhì), 不具備天然氣水合物發(fā)育的氣源條件。然而最近研究表明洋殼廣泛的蛇紋巖化作用可以產(chǎn)生大量的無機(jī)成因甲烷, 并找到了與蛇紋巖化有關(guān)的水合物發(fā)育證據(jù)。蛇紋巖化過程中不僅有甲烷生成, 還生成大量的氫氣, 很可能形成甲烷–氫氣水合物。本文根據(jù)IODP 366航次鉆探資料, 應(yīng)用甲烷–氫氣水合物的熱力學(xué)模型, 計(jì)算了馬里亞納弧前三個(gè)蛇紋巖泥火山鉆探實(shí)測(cè)獲得的不同氫氣含量條件下的甲烷–氫氣水合物穩(wěn)定帶分布特征。結(jié)果顯示氫氣比例越高, 計(jì)算獲得的甲烷–氫氣水合物穩(wěn)定帶底界深度越淺。研究的6個(gè)馬里亞納弧前蛇紋巖泥火山站位中, IODP 1491、1492、1496、1498四個(gè)站位可能具備甲烷–氫氣水合物發(fā)育溫壓和氣體成分組成條件; IODP 1493、1497站位幾乎不具備甲烷–氫氣水合物的溫壓和氣體組成條件。

    甲烷–氫氣水合物; 蛇紋巖泥火山; 無機(jī)成因甲烷; 水合物穩(wěn)定帶; 馬里亞納弧前

    0 引 言

    天然氣水合物是以烴類為主的氣體分子與水分子在合適的溫度、壓力下形成的冰晶狀固體化合物, 常見的氣體分子以甲烷為主。天然氣水合物是潛力巨大的能源(Sloan and Koh, 2008; Johnson, 2011), 并且具有重要的氣候和環(huán)境等效應(yīng)(Kennett et al., 2003; Maslin et al., 2010)。

    天然氣水合物的成藏需要合適的溫度、壓力等熱力學(xué)條件以及充足的氣源。大陸邊緣海洋環(huán)境生物生產(chǎn)力高, 海底沉積物中有大量有機(jī)質(zhì), 這些有機(jī)質(zhì)通過微生物或熱解作用轉(zhuǎn)變?yōu)樘烊粴? 為水合物提供充足的甲烷氣源(Paull and Dillon, 2001; Zheng et al., 2020), 是海洋天然氣水合物的主要分布區(qū)(Davie and Buffett, 2003; Chen et al., 2006)。但大洋區(qū)海底沉積物有機(jī)質(zhì)匱乏, 缺乏必要的氣源供給, 一直被認(rèn)為不適合天然氣水合物發(fā)育。近年深海調(diào)查研究發(fā)現(xiàn), 上地幔與洋殼的超基性巖、基性巖與水之間的相互作用可產(chǎn)生廣泛的蛇紋巖化, 并形成氫氣及甲烷等烴類氣體(Mével, 2003; Mccollom and Bach, 2009)。在超基性巖中的橄欖石和輝石等與水相互作用轉(zhuǎn)變?yōu)樯呒y石的過程中, 同時(shí)形成水鎂石、氫氣和磁鐵礦等(Mg1.8Fe0.2SiO4+1.37H2O→ 0.5Mg3Si2O5(OH)4+0.3Mg(OH)2+0.067Fe3O4+0.067H2; Mével, 2003; Evans et al., 2013)。在含碳體系中將會(huì)產(chǎn)生費(fèi)托反應(yīng)(CO2+(2+1)H2→CH2n+2+2H2O)或薩巴蒂爾反應(yīng)(CO2+4H2→CH4+2H2O), 由氫與碳在磁鐵礦催化下形成以甲烷為主的短鏈烷烴和其他有機(jī)化合物(McCollom, 1999; Proskurowski et al., 2008; 汪小妹等, 2010; 黃瑞芳等, 2015; Brovarone et al., 2017; Merdith et al., 2020)。因此, 蛇紋巖化形成的無機(jī)成因甲烷等烴類氣體可以為水合物形成提供潛在的氣源。

    此外, 對(duì)北大西洋和北冰洋之間Fram海峽超慢速擴(kuò)張脊和馬里亞納弧前蛇紋巖泥火山區(qū)的溫壓條件的研究表明, 超慢速擴(kuò)張脊和弧前蛇紋巖泥火山區(qū)具有甲烷水合物穩(wěn)定存在的熱力學(xué)條件(湯加麗等, 2020)。Fram海峽超慢速擴(kuò)張脊地震剖面上, 也發(fā)現(xiàn)了蛇紋巖化無機(jī)成因甲烷水合物發(fā)育的典型地球物理證據(jù)——似海底反射層(BSR, bottom simulating reflector), 深部蛇紋巖化成因的甲烷氣體可以通過拆離斷層向海底運(yùn)移, 直接在BSR之下聚集, 為水合物生成提供充足的甲烷氣源(Johnson et al., 2015)。

    在蛇紋巖化過程中, 不僅生成甲烷等烴類氣體, 還生成了大量的氫氣。馬里亞納弧前蛇紋巖泥火山區(qū)(圖1)大洋鉆探采集的頂空氣顯示了甲烷和氫氣共同發(fā)育, CH4/H2比值介于0.002~17.429(圖2)。一方面, 氫氣是潛在高效清潔能源; 另一方面, 氫氣的存在會(huì)改變水合物形成的整體熱力學(xué)狀態(tài), 包括整體壓力、甲烷分壓及甲烷溶解度。理論和實(shí)驗(yàn)研究表明氫氣和甲烷的混合氣在適宜的溫度、壓力等熱力學(xué)條件下可以形成甲烷–氫氣水合物(Zhang et al., 2000; Pang et al., 2012; Matsumoto et al., 2014)。因此, 研究蛇紋巖化有關(guān)的天然氣水合物生成, 有必要考慮氫氣影響。

    研究自然界甲烷–氫氣水合物, 首先要確定甲烷?氫氣水合物的形成條件和穩(wěn)定性, 因而有必要開展大洋區(qū)蛇紋巖化發(fā)育區(qū)甲烷–氫氣水合物形成條件和穩(wěn)定帶發(fā)育特征的研究。本文以馬里亞納俯沖帶弧前蛇紋巖泥火山的大洋鉆探資料為基礎(chǔ), 計(jì)算了馬里亞納弧前海底蛇紋巖泥火山區(qū)的甲烷–氫氣水合物形成條件和穩(wěn)定帶發(fā)育深度, 探討該海域蛇紋巖泥火山區(qū)甲烷–氫氣水合物發(fā)育的熱力學(xué)穩(wěn)定性, 可以為甲烷–氫氣水合物資源潛力評(píng)估提供參考。

    1 研究區(qū)地質(zhì)背景

    在西太平洋馬里亞納海溝俯沖帶弧前海底, 發(fā)育有規(guī)模巨大的蛇紋巖泥火山群, 在13°N~20°N之間距海溝軸線約30~100 km的海域分布有14個(gè)大型海底蛇紋巖泥火山, 泥火山直徑10~50 km, 高度0.5~2 km, 有大量小型蛇紋巖泥火山發(fā)育(Salisbury et al., 2002; Wheat et al., 2008; Fryer, 2012)。大洋鉆探(ODP)125和195航次及最近國際大洋發(fā)現(xiàn)計(jì)劃(IODP)366航次對(duì)馬里亞納弧前泥火山進(jìn)行了鉆探, 運(yùn)用Alvin號(hào)、Jason2和Shinkai6500等深潛器進(jìn)行了多個(gè)航次的海底調(diào)查, 發(fā)現(xiàn)蛇紋巖泥火山頂部正發(fā)育有滲漏流體活動(dòng), 分析結(jié)果顯示流體富氫氣和甲烷, 低溫(~2 ℃)、強(qiáng)堿性(pH值最高達(dá)12.5)(Fryer et al., 1990, 2018; Salisbury et al., 2002; Mottl et al., 2004; Wheat et al., 2008)。IODP 366航次對(duì)Yinazao, Fantangis?a, Asùt Tesoro三個(gè)海底蛇紋巖泥火山進(jìn)行了鉆探(圖1), 各站位水深在1243~4492 m之間、海底溫度約1.55~3.91 ℃、地溫梯度為11.7~26.5 ℃/km,采集的蛇紋巖泥火山沉積物巖芯頂空氣中的氣體以甲烷和氫氣為主, 并且不同站位和深度的氣體組成均存在差異(圖2; Fryer et al., 2018)。

    2 馬里亞納弧前蛇紋巖泥火山甲烷–氫氣水合物

    2.1 甲烷–氫氣水合物形成熱力學(xué)

    在海底沉積層中滿足一定溫度和壓力等熱力學(xué)條件使天然氣水合物能穩(wěn)定存在的范圍稱之為天然氣水合物穩(wěn)定帶(簡(jiǎn)稱“穩(wěn)定帶”), 其最深位置為天然氣水合物穩(wěn)定帶底界(BHSZ, base of hydrate stability zone)。確定穩(wěn)定帶底界是水合物調(diào)查研究的重要內(nèi)容和評(píng)估水合物資源前景的前提。穩(wěn)定帶底界深度由實(shí)際的地溫與天然氣水合物穩(wěn)定存在的三相平衡溫度交點(diǎn)確定。

    本文應(yīng)用甲烷–氫氣水合物相平衡模型(Pang et al., 2012), 確定氫氣與甲烷水合物形成溫度和壓力關(guān)系。天然氣水合物的三相平衡溫度3由壓力和氣體成分共同控制, 可以根據(jù)Pang et al. (2012)計(jì)算:

    (1)

    其中是壓力(Pa),H2是氫氣比例(%),3為壓力和氫氣比例H2對(duì)應(yīng)的三相平衡溫度(K),(,H2)為給定壓力和氫氣比例求解三相平衡溫度的函數(shù)。式(1)中的在海底沉積物中依靜水壓力式(2)計(jì)算:

    圖1 馬里亞納海溝區(qū)域地圖及Yinazao、Fantangis?a、Asùt Tesoro蛇紋巖泥火山分布位置(Fryer et al., 2018)

    CH4、H2含量來源于IODP 366航次報(bào)告(Fryer et al., 2018)。

    (2)

    其中為壓力(Pa),為海底之下的深度(m),為水深(m),為海水的密度(kg/m3),為重力加速度(9.8 m/s2)。圖3為計(jì)算的氫氣成分為1%~90%時(shí)甲烷–氫氣水合物、水和游離氣三相平衡共存的溫壓關(guān)系。

    海底之下沉積層中的溫度由海底溫度和地溫梯度確定:

    (3)

    其中是海底之下給定深度沉積層的溫度(K),0是為海底溫度(K),是地溫梯度(℃/km), 由鉆孔實(shí)測(cè)確定。

    對(duì)比實(shí)際的地溫和甲烷–氫氣水合物三相平衡溫度, 如果地溫超過三相平衡溫度(>3), 甲烷–氫氣水合物不能穩(wěn)定存在; 所處的地溫低于三相平衡溫度(<3), 甲烷–氫氣水合物可以穩(wěn)定存在; 當(dāng)?shù)販嘏c三相平衡溫度相等時(shí)(=3), 相應(yīng)的深度為甲烷–氫氣水合物穩(wěn)定帶底界深度(BHSZ)。

    圖3 甲烷–氫氣水合物、水和游離氣三相平衡溫壓關(guān)系(Pang et al., 2012)

    2.2 馬里亞納弧前蛇紋巖泥火山水合物穩(wěn)定帶底界分布

    根據(jù)IODP 366航次位于馬里亞納弧前三個(gè)蛇紋巖泥火山的6個(gè)典型站位IODP 1491、1492、1493、1496、1497、1498測(cè)定的海底之下不同氫氣比例(1%~90%), 計(jì)算了甲烷–氫氣水合物三相平衡溫度(圖4)。由于氫氣比例為99%時(shí), 計(jì)算的各站位的三相平衡溫度為負(fù)值, 各站位均不滿足甲烷–氫氣水合物存在的熱力學(xué)條件, 因此圖4只顯示了氫氣比例為1%~90%的結(jié)果。結(jié)果表明6個(gè)站位均在給定壓力時(shí), 氫氣比例越高, 相應(yīng)的三相平衡溫度越低。其中當(dāng)氫氣比例小于60%時(shí), 氫氣比例若增大10%, 相應(yīng)的甲烷–氫氣水合物的三相平衡溫度將會(huì)降低約1~2 ℃; 當(dāng)氫氣比例大于60%時(shí), 氫氣比例增加對(duì)甲烷–氫氣水合物的三相平衡溫度的降低更為顯著。

    根據(jù)地溫和甲烷–氫氣水合物的三相平衡溫度確定了其穩(wěn)定帶底界深度(表1), 氫氣比例增加, 將使甲烷–氫氣水合物的穩(wěn)定帶底界深度顯著變淺。如IODP 1492站位甲烷–氫氣水合物穩(wěn)定帶底界深度受到氫氣比例影響最為顯著, 氫氣比例增加10%, 穩(wěn)定帶底界深度變淺約5%。IODP 1491、1492、1493、1496、1497、1498站位的氫氣比例分別大于90%、90%、90%、65%、90%、90%時(shí), 甲烷–氫氣水合物穩(wěn)定帶底界深度甚至為0。因此, 確定甲烷?氫氣水合物穩(wěn)定的條件還需要考慮氣體成分的因素。根據(jù)不同成分計(jì)算的甲烷–氫氣水合物穩(wěn)定帶底界深度, 擬合了各站位氫氣比例與甲烷–氫氣水合物穩(wěn)定帶底界深度的關(guān)系曲線(圖5)。通過不同氫氣成分的甲烷–氫氣水合物的三相平衡深度曲線, 確定甲烷–氫氣水合物在各深度能穩(wěn)定存在時(shí)的最大氫氣比例。顯然只有氫氣成分低于該擬合線的值, 才具備甲烷–氫氣水合物穩(wěn)定存在的熱力學(xué)條件。因此能同時(shí)滿足甲烷–氫氣水合物穩(wěn)定存在的溫壓條件和氣體成分條件為甲烷–氫氣水合物穩(wěn)定域(圖5中灰色區(qū)域)。

    圖4 IODP 1491、1492、1493、1496、1497和1498站位溫度與甲烷–氫氣水合物三相平衡溫度確定天然氣水合物穩(wěn)定帶底界

    表1 不同氫氣比例計(jì)算的甲烷–氫氣水合物穩(wěn)定帶底界深度(m)

    注: “—”表示水合物在海底不穩(wěn)定。

    黑色點(diǎn)是不同氫氣比例在各個(gè)站位海底水深和溫度及地溫條件下所計(jì)算的甲烷–氫氣水合物三相平衡點(diǎn)的深度, 即穩(wěn)定帶底界深度; 紅線是各個(gè)站位不同氫氣比例計(jì)算的甲烷–氫氣水合物穩(wěn)定帶底界深度擬合氫氣成分–水合物穩(wěn)定帶締結(jié)深度曲線; 藍(lán)色點(diǎn)是各站位不同采樣深度頂空氣樣品對(duì)應(yīng)的氫氣比例; 灰色區(qū)域代表甲烷–氫氣水合物熱力學(xué)穩(wěn)定域。

    IODP 366航次在馬里亞納弧前蛇紋巖泥火山不同深度沉積物進(jìn)行頂空氣采樣, 不同深度和不同站位的氫氣比例存在一定的差異。計(jì)算的甲烷–氫氣水合物穩(wěn)定域由于不同的氣體成分在空間上存在顯著變化。可以歸納分為3類: ①IODP 1491和1492站位: 頂空氣測(cè)定的氫氣成分值主要在甲烷–氫氣水合物穩(wěn)定域之內(nèi), IODP 1492僅一個(gè)值位于穩(wěn)定域之外; ②IODP 1496和1498站位: 部分成分值位于甲烷–氫氣水合物穩(wěn)定域, 部分位于甲烷–氫氣水合物穩(wěn)定域之外, 這兩個(gè)站位不完全具備甲烷–氫氣水合物發(fā)育的熱力學(xué)溫壓和氣體成分條件; ③IODP 1497和1493站位: 頂空氣成分主要以氫氣為主, 處于甲烷–氫氣水合物穩(wěn)定域之外, 其中IODP 1493頂空氣僅有一個(gè)樣品的氫氣成分位于甲烷–氫氣水合物穩(wěn)定帶內(nèi), 剩余均位于甲烷–氫氣水合物非穩(wěn)定域。因此, IODP 1491和1492站位具備甲烷–氫氣水合物發(fā)育的熱力學(xué)溫壓和氣體成分條件, IODP 1496和1498站位可能具備甲烷–氫氣水合物發(fā)育溫壓和氣體成分條件, 而IODP 1497和1493站位基本不具備甲烷–氫氣水合物發(fā)育溫壓和氣體成分條件。

    3 討 論

    IODP 366航次的頂空氣采樣證實(shí)馬里亞納弧前區(qū)域蛇紋巖化泥火山存在氫氣和甲烷。氫氣和甲烷在適宜的溫度、壓力條件下可形成甲烷–氫氣水合物(Zhang et al., 2000; Pang et al., 2012; Matsumoto et al., 2014)。氫氣的存在將會(huì)改變水合物的熱力學(xué)條件, 從而影響蛇紋巖泥火山沉積物中水合物的發(fā)育特征。根據(jù)甲烷–氫氣水合物模型(Pang et al., 2012)計(jì)算不同氫氣比例的甲烷–氫氣水合物溫度壓力相平衡條件, 顯示隨著氫氣比例的增大, 等溫條件下生成水合物所需的壓力將更大, 說明氫氣的加入, 使生成水合物的溫度壓力條件要求更高(圖3)。馬里亞納弧前蛇紋巖泥火山鉆探獲得的氣體組成、海底深度、溫度和地溫梯度顯示該區(qū)域存在甲烷–氫氣水合物生成的熱力學(xué)條件。鑒于氫氣含量可能較高, 氫氣相對(duì)甲烷不利于水合物穩(wěn)定, 甲烷含量可能是控制水合物能否發(fā)育的關(guān)鍵因素。

    氫氣不僅改變水合物生成的溫壓條件, 還可以改變游離氣生成水合物的甲烷分壓或溶解氣生成水合物–水二相甲烷溶解度。為評(píng)估甲烷和氫氣生成水合物所需甲烷的變化, 選取6個(gè)溫度來確定每個(gè)溫度條件下生成甲烷–氫氣水合物的氫氣比例與相應(yīng)的甲烷分壓關(guān)系(圖6), 其中甲烷的分壓利用生成甲烷–氫氣水合物的壓力(溫度控制)與甲烷氣體的比例確定。結(jié)果表明, 隨著溫度的升高, 甲烷–氫氣水合物生成所需要的甲烷分壓增大; 隨氫氣比例增大, 生成水合物所需的甲烷分壓變化很小, 與純甲烷水合物所需要的壓力非常接近(偏差小于10%)。

    圖6 等溫條件下甲烷–氫氣水合物甲烷分壓與氫氣比例關(guān)系

    IODP航次1492A站位實(shí)測(cè)頂空氣甲烷濃度隨深度逐漸增大, 通過線性擬合頂空氣甲烷濃度與深度關(guān)系, 推測(cè)穩(wěn)定帶內(nèi)頂空氣甲烷濃度具有類似大陸邊緣水合物發(fā)育層實(shí)測(cè)的頂空氣甲烷含量特征(湯加麗等, 2020)。由于甲烷–氫氣水合物中的甲烷含量(甲烷分壓)與純甲烷水合物的類似, 因此IODP 1492A站位具有形成甲烷水合物的甲烷含量條件, 在該站位蛇紋巖泥火山沉積物中可能發(fā)育甲烷–氫氣水合物。

    4 結(jié) 論

    本文利用馬里亞納弧前蛇紋巖泥火山區(qū)多個(gè)IODP站位的溫壓條件和頂空氣成分, 分析了甲烷和氫氣生成天然氣水合物的熱力學(xué)條件。結(jié)果顯示氫氣的含量控制了甲烷–氫氣水合物的熱力學(xué)穩(wěn)定性及穩(wěn)定帶底界深度。氫氣含量越高, 甲烷–氫氣水合物的熱力學(xué)穩(wěn)定性越低, 計(jì)算的甲烷–氫氣水合物穩(wěn)定帶底界越淺。其中馬里亞納弧前蛇紋巖區(qū)IODP 1491、1492、1496和1498可能具備甲烷–氫氣水合物發(fā)育的溫壓和氣體成分條件, 而IODP 1493、1497以氫氣為主, 可能不具備甲烷–氫氣水合物生成條件。

    致謝:感謝中國石油大學(xué)(北京)陳光進(jìn)課題組對(duì)本論文的幫助, 感謝兩位匿名審稿人提出的寶貴修改建議。

    黃瑞芳, 孫衛(wèi)東, 丁興, 劉吉強(qiáng), 詹文歡. 2015. 橄欖巖蛇紋石化過程中氫氣和烷烴的形成. 巖石學(xué)報(bào), 31(7): 1901–1907.

    湯加麗, 曹運(yùn)誠, 陳多福. 2020. 大洋蛇紋巖化無機(jī)成因甲烷水合物穩(wěn)定帶底界模擬. 海洋地質(zhì)與第四紀(jì)地質(zhì), 40(4): 107–115.

    汪小妹, 曾志剛, 歐陽荷根, 殷學(xué)博, 王曉媛, 陳帥, 張國良, 武力. 2010. 大洋橄欖巖的蛇紋巖石化研究進(jìn)展評(píng)述. 地球科學(xué)進(jìn)展, 25(6): 605–616.

    Brovarone A V, Martinez I, Elmaleh A, Compagnoni R, Chaduteau C, Ferraris C, Esteve I. 2017. Massive production of abiotic methane during subduction evidenced in metamorphosed ophicarbonates from the Italian Alps., 8, 14134.

    Chen D F, Su Z, Cathles L. 2006. Types of gas hydrates in marine environments and their thermodynamic characteristics., 17: 723–737.

    Davie M K, Buffett B A. 2003. Sources of methane for marine gas hydrate: Inferences from a comparison of observations and numerical models., 206(1): 51–63.

    Evans B W, Hattori K, Baronnet A. 2013. Serpentinite: What, why, where?, 9(2): 99–106.

    Fryer P. 2012. Serpentinite mud volcanism: Observations, processes, and implications., 4: 345–373.

    Fryer P, Pearce J A, Stokking L B. 1990. Proceedings of the Ocean Drilling Program. Initial Reports, vol 125. College Station, TX: 367–380.

    Fryer P, Wheat C G, Williams T, Expedition 366 Scientists. 2018. Expedition 366 preliminary report: Mariana Convergent Margin and South Chamorro Seamount. International Ocean Discovery Program.

    Johnson A H. 2011. Global resource potential of gas hydrate — A new calculation.11(2): 1–4.

    Johnson J E, Mienert J, Plaza-Faverola A, Vadakkepuliyambatta S, Knies J, Bunz S, Andreassen K, Ferre B. 2015. Abioticmethane from ultraslow-spreading ridges can charge Arctic gas hydrates., 43(5): 371–374.

    Kennett J P, Cannariato K G, Hendy I L, Behl R J. 2003. The clathrate gun hypothesis. American Geophysical Union, Washington: DC: 105–107.

    Maslin M, Owen M, Betts R, Day S, Dunkley Jones T, Ridgwell A. 2010. Gas hydrates: Past and future geohazard?, 368: 2369–2393.

    Matsumoto Y, Grim R G, Khan N M, Sugahara T, Ohgaki K, Sloan E D, Koh C A, Sum A K. 2014. Investigating the thermodynamic stabilities of hydrogen and methane binary gas hydrates., 118(7): 3783–3788.

    McCollom T M. 1999. Methanogenesis as a potential source of chemical energy for primary biomass production by autotrophic organisms in hydrothermal systems on Europa., 104: 30729–30742.

    McCollom T M, Bach W. 2009. Thermodynamic constraints on hydrogen generation during serpentinization of ultramafic rocks., 73(3): 856–875.

    Merdith A S, Garc?-a D P, Isabelle D, Muriel A, Nicky W M, Nicolas C. 2020. Pulsated global hydrogen and methane fluxat mid-ocean ridges driven by Pangea breakup,,, 21(4), e2019GC008869.

    Mével C. 2003. Serpentinization of abyssal peridotites at mid-ocean ridges., 335(10–11): 825–852.

    Mottl M J, Wheat C G, Fryer P, Gharib J, Martin J B. 2004. Chemistry of springs across the Mariana forearc shows progressive devolatilization of the subducting plate., 68(23): 4915–4933.

    Pang J Y, Ng H-J, Zuo J L, Zhang D, Ma Q L, Chen G J. 2012. Hydrogen gas hydrate — Measurements and predictions., 316: 6–10.

    Paull C K, Dillon W P. 2001. Natural gas hydrates: Occurrence, distribution, and detection. American Geophysical Union, Washington: DC: 1–315.

    Proskurowski G, Lilley M, Seewald J S, Früh-Green G, Olson E J, Lupton J, Sylva S, Kelley D. 2008. Abiogenic hydrocarbon production at lost city hydrothermal field., 319: 604–607.

    Salisbury M H, Shinohara M, Richter C, Shipboard Scientific Party. 2002. Proceedings of the Ocean Drilling Program. Initial Reports, vol 195. College Station, TX: 1–63.

    Sloan E D, Koh C A. 2008. Clathrate hydrates of natural gases (3rd edition). New Nork: CRC press: 1–701.

    Wheat C G, Fryer P, Fisher A T, Hulme S, Jannasch H, Mottl M J, Becker K. 2008. Borehole observations of fluid flow from South Chamorro Seamount, an active serpentinitemud volcano in the Mariana forearc., 267(3–4): 401–409.

    Zhang S X, Chen G J, Ma C F, Yang L Y, Guo T M. 2000. Hydrate formation of hydrogen + hydrocarbon gas mixtures., 45(5): 908–911.

    Zheng Z H, Cao Y C, Xu W Y, Chen D F. 2020. A numerical model to estimate the effects of variable sedimentation rates on methane hydrate formation — Application to the ODP Site 997 on Blake ridge, southeastern north American continental margin.:, 125, e2019JB018851.

    The abiotic methane-hydrogen hydrate stability zone in Mariana forearc serpentinite mud volcanoes

    TANG Jiali, CAO Yuncheng*, CHEN Duofu

    (Shanghai Engineering Research Center of Hadal Science and Technology, College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China)

    Although thermodynamic conditions for hydrate stability are satisfied throughout most of the world’s oceans, most marine hydrates have been discovered along continental margins where methane is supplied by the conversion of organic carbon within sediments. Abiotic methane has recently been recognized as an additional sub-seafloor gas source generated during the serpentinization of ultramafic rocks. A large amount of hydrogen is also produced during serpentinization, which forms hydrogen gas hydrate, implying that thermodynamic conditionsand gas resources are quite different from that of previously recovered methane hydrate. In this study, we estimated the methane-hydrogen hydrate formation conditions and the base depth of methane-hydrogen hydrate at six sites drilled at the Mariana forearc serpentinite mud volcano for gas mixtures with concentrations of hydrogen from 1% to 90%. Our calculations show that the base of the methane-hydrogen hydrate stability zone would be shallow and the partial pressure of methane needed for methane-hydrogen hydrate formation is constant at a given temperature and pressure when the ratio of hydrogen is high. Additionally, the methane-hydrogen hydrate formation conditions can be satisfied at four sites: IODP 1491, 1492, 1496, and 1498. Hydrogen was the dominant gas at IODP sites 1493 and 1497, where methane-hydrogen hydrates could not be formed.

    methane-hydrogen hydrate; serpentinite mud volcano; abiotic methane; hydrate stability zone; Mariana forearc

    P67

    A

    0379-1726(2022)02-0194-08

    10.19700/j.0379-1726.2022.02.003

    2020-03-14;

    2020-07-27

    國家自然科學(xué)基金項(xiàng)目(41776050, 42176069, 41776080)資助。

    湯加麗(1994–), 女, 碩士研究生, 海洋地質(zhì)專業(yè)。E-mail: 2363589516@qq.com

    曹運(yùn)誠(1983–), 男, 副研究員, 主要從事天然氣水合物數(shù)值模擬研究。E-mail: yccao@shou.edu.cn

    猜你喜歡
    泥火山馬里亞納水合物
    深海魅影:馬里亞納海溝里的奇怪生物
    馬里亞納海溝的奇怪生物
    軍事文摘(2021年20期)2021-11-10 01:58:54
    氣井用水合物自生熱解堵劑解堵效果數(shù)值模擬
    軍器圖說 馬里亞納海戰(zhàn)
    熱水吞吐開采水合物藏?cái)?shù)值模擬研究
    天然氣水合物保壓轉(zhuǎn)移的壓力特性
    我國海域天然氣水合物試采成功
    親近泥火山
    水潤(rùn)的泥火山
    準(zhǔn)噶爾盆地南緣泥火山與油氣苗
    在线观看一区二区三区| 搞女人的毛片| 嫩草影视91久久| 搞女人的毛片| 老司机深夜福利视频在线观看| 国产伦在线观看视频一区| 熟妇人妻久久中文字幕3abv| 伊人久久精品亚洲午夜| 熟妇人妻久久中文字幕3abv| 99国产精品一区二区三区| 成人高潮视频无遮挡免费网站| 看黄色毛片网站| 成人永久免费在线观看视频| 亚洲午夜理论影院| 亚洲av二区三区四区| 亚洲精品在线美女| 好男人电影高清在线观看| 宅男免费午夜| 天天躁日日操中文字幕| 国产 一区 欧美 日韩| 内射极品少妇av片p| 午夜福利高清视频| 精品久久国产蜜桃| 午夜激情欧美在线| 波多野结衣巨乳人妻| 两个人的视频大全免费| www.999成人在线观看| 最后的刺客免费高清国语| 首页视频小说图片口味搜索| 日韩精品中文字幕看吧| 最好的美女福利视频网| 欧美日韩乱码在线| 亚洲av电影不卡..在线观看| 欧美乱色亚洲激情| www日本黄色视频网| 久久久精品欧美日韩精品| 午夜精品久久久久久毛片777| 又黄又爽又刺激的免费视频.| 精品不卡国产一区二区三区| 黄片小视频在线播放| 中亚洲国语对白在线视频| 美女高潮的动态| 天堂av国产一区二区熟女人妻| 国产一级毛片七仙女欲春2| 国产精品一及| 亚洲成av人片免费观看| 久久久精品大字幕| 91午夜精品亚洲一区二区三区 | 搡女人真爽免费视频火全软件 | 美女免费视频网站| 欧美激情在线99| 久久久国产成人精品二区| 午夜免费男女啪啪视频观看 | 夜夜看夜夜爽夜夜摸| 婷婷六月久久综合丁香| 青草久久国产| 少妇的逼水好多| 精品人妻偷拍中文字幕| 免费看美女性在线毛片视频| 欧美另类亚洲清纯唯美| 人妻制服诱惑在线中文字幕| 成人三级黄色视频| 免费人成视频x8x8入口观看| 亚洲欧美精品综合久久99| 国产伦一二天堂av在线观看| 色在线成人网| 免费人成视频x8x8入口观看| 国产大屁股一区二区在线视频| 午夜两性在线视频| 97超视频在线观看视频| 久久久久国内视频| 99久久精品一区二区三区| 我的老师免费观看完整版| 亚州av有码| 国产亚洲av嫩草精品影院| 国产综合懂色| 日本免费一区二区三区高清不卡| 久久久久久国产a免费观看| 亚洲人成电影免费在线| 欧美最黄视频在线播放免费| 午夜福利在线在线| 亚洲av一区综合| 一级av片app| 国产精品亚洲美女久久久| 人人妻人人澡欧美一区二区| 免费看美女性在线毛片视频| 久久精品国产亚洲av天美| 一个人看的www免费观看视频| 国产一级毛片七仙女欲春2| 日本 av在线| 精品福利观看| 丰满的人妻完整版| 午夜精品久久久久久毛片777| 99热精品在线国产| 1000部很黄的大片| 露出奶头的视频| 亚洲成a人片在线一区二区| 啦啦啦韩国在线观看视频| 好男人电影高清在线观看| 国产精品爽爽va在线观看网站| 一进一出抽搐动态| 99国产精品一区二区三区| eeuss影院久久| 欧美激情在线99| 在线播放国产精品三级| 色哟哟哟哟哟哟| 变态另类丝袜制服| 首页视频小说图片口味搜索| 在线观看免费视频日本深夜| 国产一区二区激情短视频| 亚洲人与动物交配视频| 内射极品少妇av片p| 久久精品夜夜夜夜夜久久蜜豆| 免费搜索国产男女视频| 久久久久久久亚洲中文字幕 | 国产v大片淫在线免费观看| 婷婷亚洲欧美| 香蕉av资源在线| 男插女下体视频免费在线播放| 国产精品电影一区二区三区| 久久精品人妻少妇| 最近中文字幕高清免费大全6 | 男女之事视频高清在线观看| 制服丝袜大香蕉在线| 亚洲av一区综合| 欧美乱妇无乱码| 91麻豆av在线| 日韩大尺度精品在线看网址| 欧美激情久久久久久爽电影| 久久热精品热| 少妇丰满av| 亚州av有码| 成人午夜高清在线视频| 熟女人妻精品中文字幕| 国产 一区 欧美 日韩| 国产精品一区二区三区四区免费观看 | 国产一区二区亚洲精品在线观看| x7x7x7水蜜桃| 亚洲专区国产一区二区| 国产一区二区激情短视频| xxxwww97欧美| 白带黄色成豆腐渣| 亚洲av成人av| 久久精品国产99精品国产亚洲性色| 一个人看的www免费观看视频| 在线免费观看的www视频| 国产麻豆成人av免费视频| 国产亚洲精品av在线| 97人妻精品一区二区三区麻豆| 久久久久久久精品吃奶| 窝窝影院91人妻| 国产野战对白在线观看| 桃红色精品国产亚洲av| 国产熟女xx| 日韩中字成人| 亚洲欧美清纯卡通| 一本综合久久免费| 国产91精品成人一区二区三区| 精品欧美国产一区二区三| 91在线观看av| 久久亚洲精品不卡| 亚洲人成电影免费在线| 久久国产精品影院| 国产国拍精品亚洲av在线观看| 美女 人体艺术 gogo| 黄色一级大片看看| 精品乱码久久久久久99久播| 女人被狂操c到高潮| 欧美日韩综合久久久久久 | 嫩草影院精品99| 国产v大片淫在线免费观看| 日韩中文字幕欧美一区二区| 九九热线精品视视频播放| 国产成人欧美在线观看| 午夜两性在线视频| 欧美日韩乱码在线| 日韩大尺度精品在线看网址| 男人舔女人下体高潮全视频| 日韩 亚洲 欧美在线| 亚洲一区二区三区不卡视频| 免费电影在线观看免费观看| 两人在一起打扑克的视频| 网址你懂的国产日韩在线| 99热只有精品国产| 又爽又黄无遮挡网站| 国产老妇女一区| 国产一区二区在线av高清观看| 又爽又黄无遮挡网站| 精品不卡国产一区二区三区| 免费大片18禁| 美女被艹到高潮喷水动态| 亚洲乱码一区二区免费版| 国产精品久久视频播放| 91狼人影院| 国产精品av视频在线免费观看| 一本一本综合久久| 国产精品98久久久久久宅男小说| 淫秽高清视频在线观看| 乱码一卡2卡4卡精品| 波多野结衣高清无吗| 女人十人毛片免费观看3o分钟| 精品久久久久久久久久免费视频| 亚洲精品成人久久久久久| 美女cb高潮喷水在线观看| 天堂av国产一区二区熟女人妻| 精品久久久久久久久亚洲 | 真实男女啪啪啪动态图| 嫩草影视91久久| 色播亚洲综合网| 脱女人内裤的视频| 精品国产三级普通话版| 很黄的视频免费| 免费av不卡在线播放| 有码 亚洲区| 老司机午夜福利在线观看视频| 一区二区三区激情视频| 亚洲七黄色美女视频| 久久国产精品人妻蜜桃| 午夜免费成人在线视频| 久久性视频一级片| 久久久久国产精品人妻aⅴ院| 国产欧美日韩精品亚洲av| 欧美色欧美亚洲另类二区| 在线a可以看的网站| 久久中文看片网| 嫩草影院入口| 日韩欧美在线乱码| 婷婷色综合大香蕉| 欧美激情在线99| 国产又黄又爽又无遮挡在线| 丁香欧美五月| 免费在线观看亚洲国产| 午夜福利在线观看免费完整高清在 | 国产大屁股一区二区在线视频| 人人妻,人人澡人人爽秒播| 欧美xxxx黑人xx丫x性爽| 波多野结衣巨乳人妻| 他把我摸到了高潮在线观看| 国产v大片淫在线免费观看| 亚洲 欧美 日韩 在线 免费| 免费人成视频x8x8入口观看| 亚洲一区二区三区不卡视频| 国产黄a三级三级三级人| 亚洲欧美精品综合久久99| 99久久久亚洲精品蜜臀av| 一区二区三区四区激情视频 | 狠狠狠狠99中文字幕| 精品一区二区三区av网在线观看| 男女下面进入的视频免费午夜| 日韩有码中文字幕| 欧美三级亚洲精品| 毛片一级片免费看久久久久 | 国语自产精品视频在线第100页| 国产黄a三级三级三级人| 色尼玛亚洲综合影院| 日韩欧美国产一区二区入口| 在现免费观看毛片| 中文亚洲av片在线观看爽| 国内毛片毛片毛片毛片毛片| 精品久久久久久久人妻蜜臀av| 国产精品一及| 又爽又黄无遮挡网站| 国产精品国产高清国产av| 91午夜精品亚洲一区二区三区 | 一区福利在线观看| 国产精品精品国产色婷婷| 色视频www国产| 中文资源天堂在线| 久久人人精品亚洲av| 午夜亚洲福利在线播放| 麻豆av噜噜一区二区三区| 欧美激情国产日韩精品一区| 村上凉子中文字幕在线| 嫩草影院精品99| 又爽又黄无遮挡网站| 久久精品国产亚洲av涩爱 | 国内毛片毛片毛片毛片毛片| 美女高潮的动态| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区四区激情视频 | 国产精品久久视频播放| 黄色配什么色好看| 天美传媒精品一区二区| 国产亚洲精品综合一区在线观看| 九色成人免费人妻av| 亚洲精品乱码久久久v下载方式| 久久欧美精品欧美久久欧美| 国产精品美女特级片免费视频播放器| 亚洲 欧美 日韩 在线 免费| 97热精品久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 听说在线观看完整版免费高清| 人妻丰满熟妇av一区二区三区| 草草在线视频免费看| 国产亚洲欧美98| 日本免费一区二区三区高清不卡| 成年人黄色毛片网站| 亚洲国产精品合色在线| 最近中文字幕高清免费大全6 | 国产黄片美女视频| 在线观看美女被高潮喷水网站 | 五月玫瑰六月丁香| av黄色大香蕉| www日本黄色视频网| x7x7x7水蜜桃| 亚洲精品在线观看二区| 欧美日韩福利视频一区二区| 麻豆久久精品国产亚洲av| 一区二区三区免费毛片| 成人国产综合亚洲| 一区二区三区四区激情视频 | 天天躁日日操中文字幕| 亚洲欧美日韩东京热| 国产精品爽爽va在线观看网站| 脱女人内裤的视频| 欧美精品国产亚洲| 久久久久久久亚洲中文字幕 | 成人精品一区二区免费| АⅤ资源中文在线天堂| 日本一本二区三区精品| 亚洲18禁久久av| 黄色丝袜av网址大全| 亚洲精品成人久久久久久| 麻豆av噜噜一区二区三区| 亚洲专区中文字幕在线| 国产成人av教育| av在线蜜桃| 国模一区二区三区四区视频| 婷婷亚洲欧美| 欧美在线一区亚洲| 在线十欧美十亚洲十日本专区| 国产精品女同一区二区软件 | 亚洲av成人不卡在线观看播放网| 精品熟女少妇八av免费久了| 每晚都被弄得嗷嗷叫到高潮| 免费观看的影片在线观看| av福利片在线观看| 欧美黑人巨大hd| 久久久久国产精品人妻aⅴ院| 国产精品女同一区二区软件 | 亚洲精品一区av在线观看| 亚洲最大成人手机在线| 精品国产亚洲在线| 一本久久中文字幕| 夜夜看夜夜爽夜夜摸| 国模一区二区三区四区视频| 小说图片视频综合网站| 日韩欧美精品免费久久 | 欧美日韩福利视频一区二区| 老熟妇仑乱视频hdxx| 午夜精品在线福利| 亚洲精品在线美女| 久久久色成人| 日韩欧美 国产精品| 毛片女人毛片| av在线蜜桃| 国产精品嫩草影院av在线观看 | av中文乱码字幕在线| 国产麻豆成人av免费视频| 在线观看av片永久免费下载| 热99re8久久精品国产| 国产精品综合久久久久久久免费| 亚洲成人久久爱视频| 日韩欧美 国产精品| 一本一本综合久久| 久久热精品热| 午夜激情福利司机影院| 色在线成人网| 如何舔出高潮| 久久精品夜夜夜夜夜久久蜜豆| 精品99又大又爽又粗少妇毛片 | 夜夜爽天天搞| 国产精品国产高清国产av| 亚洲国产欧洲综合997久久,| 国产成人啪精品午夜网站| 精品久久久久久久久亚洲 | 亚洲国产色片| 成年女人看的毛片在线观看| 性插视频无遮挡在线免费观看| 婷婷亚洲欧美| 一区二区三区免费毛片| 亚洲第一区二区三区不卡| a在线观看视频网站| 在线观看免费视频日本深夜| 久久久色成人| 三级国产精品欧美在线观看| 能在线免费观看的黄片| 波多野结衣高清无吗| 少妇人妻精品综合一区二区 | 桃红色精品国产亚洲av| 国内精品美女久久久久久| 最近在线观看免费完整版| 国产午夜福利久久久久久| 美女cb高潮喷水在线观看| 不卡一级毛片| 99在线视频只有这里精品首页| 国产免费一级a男人的天堂| 麻豆国产av国片精品| 两个人的视频大全免费| 无人区码免费观看不卡| 最好的美女福利视频网| 极品教师在线视频| 热99在线观看视频| 国产又黄又爽又无遮挡在线| 一二三四社区在线视频社区8| 脱女人内裤的视频| 九九热线精品视视频播放| 中文资源天堂在线| 最近中文字幕高清免费大全6 | 日本黄大片高清| 欧美色视频一区免费| 国产精品综合久久久久久久免费| x7x7x7水蜜桃| 3wmmmm亚洲av在线观看| 高清在线国产一区| 久久久久久久久大av| 乱码一卡2卡4卡精品| 亚洲av免费在线观看| 成人欧美大片| 我要看日韩黄色一级片| 国产精品一区二区三区四区免费观看 | 亚洲成人中文字幕在线播放| 欧美极品一区二区三区四区| 久久婷婷人人爽人人干人人爱| 别揉我奶头 嗯啊视频| 国产精品,欧美在线| 国产精品一区二区免费欧美| 麻豆国产97在线/欧美| 亚洲成av人片在线播放无| 免费看美女性在线毛片视频| 精品福利观看| 久久性视频一级片| 精品人妻偷拍中文字幕| 久久草成人影院| 男女下面进入的视频免费午夜| 精品人妻一区二区三区麻豆 | 亚洲成人久久性| bbb黄色大片| 自拍偷自拍亚洲精品老妇| 蜜桃亚洲精品一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 99久久成人亚洲精品观看| 黄片小视频在线播放| 国产探花在线观看一区二区| 欧美性猛交黑人性爽| 老司机深夜福利视频在线观看| 69av精品久久久久久| 在线天堂最新版资源| 麻豆成人午夜福利视频| 午夜福利在线在线| 能在线免费观看的黄片| 全区人妻精品视频| 1024手机看黄色片| 人人妻人人澡欧美一区二区| 中文字幕免费在线视频6| 久久中文看片网| 啦啦啦韩国在线观看视频| 午夜福利免费观看在线| 日韩精品中文字幕看吧| 午夜免费激情av| 日韩人妻高清精品专区| 国产精品伦人一区二区| 久久这里只有精品中国| 免费人成在线观看视频色| 欧美日韩中文字幕国产精品一区二区三区| 欧美乱色亚洲激情| 最近最新免费中文字幕在线| 亚洲一区二区三区色噜噜| 国内精品久久久久久久电影| 亚洲成人中文字幕在线播放| 亚洲天堂国产精品一区在线| 亚洲国产欧美人成| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人影院久久av| 91久久精品电影网| xxxwww97欧美| 色综合亚洲欧美另类图片| 欧美一级a爱片免费观看看| 亚洲狠狠婷婷综合久久图片| 99久久九九国产精品国产免费| 午夜免费激情av| 成人国产一区最新在线观看| 中出人妻视频一区二区| 88av欧美| 91在线观看av| 日本 欧美在线| 91久久精品电影网| 亚洲精品456在线播放app | 日韩高清综合在线| 国产高潮美女av| 亚洲成a人片在线一区二区| 亚洲美女视频黄频| 精品国内亚洲2022精品成人| 亚洲性夜色夜夜综合| 国产白丝娇喘喷水9色精品| 嫁个100分男人电影在线观看| 日本免费一区二区三区高清不卡| 成人亚洲精品av一区二区| 亚洲av.av天堂| 国产激情偷乱视频一区二区| 51午夜福利影视在线观看| 三级男女做爰猛烈吃奶摸视频| 精品熟女少妇八av免费久了| 丁香欧美五月| av在线老鸭窝| 我要搜黄色片| 欧美区成人在线视频| 美女 人体艺术 gogo| 午夜福利欧美成人| 亚洲欧美日韩高清在线视频| 18禁在线播放成人免费| av天堂中文字幕网| 精品无人区乱码1区二区| av欧美777| 亚洲三级黄色毛片| 岛国在线免费视频观看| 亚洲 欧美 日韩 在线 免费| 欧美黄色淫秽网站| 国产私拍福利视频在线观看| 国产精品人妻久久久久久| 最近最新免费中文字幕在线| 精品人妻偷拍中文字幕| 亚洲国产精品久久男人天堂| 少妇的逼好多水| 一个人看视频在线观看www免费| 美女cb高潮喷水在线观看| 婷婷精品国产亚洲av在线| www.熟女人妻精品国产| 成人特级av手机在线观看| 亚洲欧美日韩东京热| 国产欧美日韩一区二区精品| 久久国产乱子免费精品| 人妻丰满熟妇av一区二区三区| 亚洲av一区综合| 九九在线视频观看精品| 国产视频内射| 国产91精品成人一区二区三区| 韩国av一区二区三区四区| 最近最新免费中文字幕在线| 成人国产一区最新在线观看| 黄色丝袜av网址大全| 亚洲无线在线观看| 女同久久另类99精品国产91| 亚洲欧美激情综合另类| 九九久久精品国产亚洲av麻豆| 99久久精品一区二区三区| 丝袜美腿在线中文| 又粗又爽又猛毛片免费看| 亚洲久久久久久中文字幕| 精品久久久久久久久av| 亚洲精品一卡2卡三卡4卡5卡| 熟女人妻精品中文字幕| 国产免费av片在线观看野外av| 精品久久久久久久人妻蜜臀av| 男女视频在线观看网站免费| 国产精品久久久久久亚洲av鲁大| 99久久成人亚洲精品观看| 欧美性猛交╳xxx乱大交人| .国产精品久久| 99热6这里只有精品| 亚洲综合色惰| 亚洲午夜理论影院| 国产亚洲精品久久久com| 日韩成人在线观看一区二区三区| 中文字幕久久专区| 国产高清有码在线观看视频| 国产一区二区三区视频了| 欧美日韩亚洲国产一区二区在线观看| 变态另类丝袜制服| 黄色日韩在线| 三级毛片av免费| www.色视频.com| 国产熟女xx| 久久性视频一级片| 久久精品国产自在天天线| 国产欧美日韩精品一区二区| 一个人观看的视频www高清免费观看| 日韩免费av在线播放| 免费黄网站久久成人精品 | 在线观看av片永久免费下载| 亚洲精品乱码久久久v下载方式| 久久久精品欧美日韩精品| 99视频精品全部免费 在线| 男人狂女人下面高潮的视频| 久久人人爽人人爽人人片va | 亚洲av电影在线进入| 三级男女做爰猛烈吃奶摸视频| 欧美色欧美亚洲另类二区| 亚洲自偷自拍三级| 国内精品一区二区在线观看| 日韩 亚洲 欧美在线| 精品一区二区免费观看| 日日摸夜夜添夜夜添av毛片 | 成人欧美大片| 最近最新免费中文字幕在线| 亚洲真实伦在线观看| 日韩大尺度精品在线看网址| 色播亚洲综合网| 最好的美女福利视频网| 亚洲在线观看片| 亚洲aⅴ乱码一区二区在线播放| 97碰自拍视频| 精品人妻一区二区三区麻豆 | 国产不卡一卡二| 中文字幕高清在线视频| www.www免费av| 亚洲欧美日韩无卡精品| 中文在线观看免费www的网站| 欧美日韩乱码在线| 亚洲 国产 在线| 中文字幕人妻熟人妻熟丝袜美| 亚洲av电影在线进入| 精品一区二区三区视频在线观看免费| 亚洲男人的天堂狠狠| 中文字幕熟女人妻在线|