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

    凝灰?guī)r成因及儲集空間類型

    2025-03-04 00:00:00楊鳳英張云峰曾琪潘清川范存輝李一葦尹宏張本健
    沉積學報 2025年1期

    摘 要 【目的】川西北地區(qū)二疊系吳家坪組海相碳酸鹽巖中廣泛分布厚度不等的凝灰?guī)r,并具有一定的天然氣儲集能力。川西北大坪剖面吳家坪組剖面厚45.02 m,在吳一段和吳二段中發(fā)育單層厚5~55 cm的凝灰?guī)r,分布穩(wěn)定,為研究該區(qū)凝灰?guī)r成因及儲集空間類型提供了良好的素材。【方法】 基于野外剖面實測,并采集凝灰?guī)r樣品進行全巖礦物含量、地球化學特征、鋯石微量元素、U-Pb年齡等分析對比,探討凝灰?guī)r成因及儲集空間類型?!窘Y(jié)果】吳家坪組凝灰?guī)r以火山塵凝灰?guī)r為主,由粒徑小于0.05 mm的玻屑組成;礦物成分以伊利石、伊蒙混層為主,含少量石英、石膏;吳一段凝灰?guī)r具有Al2O3/TiO2 高的比值(23.36~34.56)、SiO2/Al2O3的低比值(1.84~2.75)及Eu、Ba、Sr、Ti元素明顯虧損等特征,吳二段凝灰?guī)r具有Al2O3/TiO2高的比值(21.140~53.049)、SiO2/Al2O3的低比值(1.829~2.682)及Eu弱虧損、虧損Nb、Ti等高場強元素?!窘Y(jié)論】吳家坪組凝灰?guī)r受到蝕變程度較高,其來源均與峨眉山大火成巖省的多期酸性火山作用有關(guān),形成于板內(nèi)拉張環(huán)境,且吳二段凝灰?guī)r受到殼?;旌系挠绊?。凝灰?guī)r中儲集空間類型包括洞、不規(guī)則溶孔、粒內(nèi)溶孔、有機質(zhì)孔、粒間溶孔、裂縫等,凝灰?guī)r中熱液鋯石測得的U-Pb年齡為132.37±0.53 Ma,其形成受早白堊世熱液、構(gòu)造活動等多種因素影響。

    關(guān)鍵詞 凝灰?guī)r成因;儲集空間類型;峨眉山大火成巖?。粎羌移航M;二疊系;川西北地區(qū)

    第一作者簡介 楊鳳英,女,1998年出生,碩士,儲層地質(zhì)學,E-mail: yangfengying1122@163.com

    通信作者 張云峰,男,副教授,沉積學、儲層地質(zhì)學,E-mail: zhyf@swpu.edu.cn

    中圖分類號 P588.14 文獻標志碼 A

    0 引言

    凝灰?guī)r是火山噴發(fā)過程中產(chǎn)生細粒碎屑物質(zhì),通過風或者水的搬運作用及固結(jié)壓實作用形成的一類火山碎屑巖,因其在大氣中滯留的時間較短(1~3年),且橫向上具有連續(xù)性,是極好的地層對比標志層[1],其縱橫向分布亦可約束火山活動時限及規(guī)模[2?3]。四川盆地及華南地區(qū)晚二疊世凝灰?guī)r廣泛分布[4],成分以伊利石、伊蒙混層、高嶺石等黏土礦物為主,含少量石英、黃鐵礦等[5?6];主要分布于四川盆地的西北部、西部、南部[7?9]、重慶南川[10]等。前人基于鋯石U-Pb年齡、礦物組分、微量元素分析等,認為上二疊統(tǒng)凝灰?guī)r成因與峨眉山大火成巖省酸性火山活動[11?12]、華南板塊西南緣古特提斯二疊紀大陸巖漿弧[13?14]、西伯利亞大火成巖省相關(guān)的基性火山活動[15]等相關(guān),并廣泛用于火山噴發(fā)活動時間及期次的確定。

    受沉積及成巖作用的影響,凝灰?guī)r中普遍存在不同尺度原生孔隙及次生孔隙,前者包括收縮孔、晶間孔等,其發(fā)育程度受巖性及火山噴發(fā)旋回的控制[16?17],后者包括溶蝕孔洞、粒內(nèi)溶孔、粒間溶孔、有機質(zhì)孔、脫?;?、裂縫等,發(fā)育程度受熱液作用、脫?;饔谩⒍嗥谌芪g、破裂作用[18?19]等的影響。國內(nèi)目前已在四川盆地[20]、三塘湖盆地[21]、鄂爾多斯盆地[22?23]等多個地區(qū)發(fā)現(xiàn)凝灰?guī)r儲層及油氣藏。近年來,川西北地區(qū)S101、元壩7、YS1等井的上二疊統(tǒng)凝灰?guī)r發(fā)現(xiàn)規(guī)模凝灰?guī)r孔隙型、溶洞型、裂縫型儲層及工業(yè)氣流[24?26],顯示了較大的儲集能力及天然氣工業(yè)產(chǎn)出價值;與此同時,由于鉆井實物資料的缺乏,該區(qū)凝灰?guī)r成因、儲層形成機制等缺乏系統(tǒng)研究。

    選取出露條件好、特征典型的川西北大坪剖面二疊系吳家坪組凝灰?guī)r作為研究對象,基于剖面實測及密集采樣,開展了巖石學、地球化學等方面的測試,獲取了凝灰?guī)r發(fā)育位置、厚度、礦物類型、微量元素、鋯石U-Pb年齡、不同尺度儲集空間類型等數(shù)據(jù),明確了凝灰?guī)r的巖石類型及礦物組成、成因、不同巖石類型的儲集性等,以期為川西北地區(qū)上二疊統(tǒng)凝灰?guī)r成因、分布、儲集性及天然氣勘探提供支撐。

    1 地質(zhì)背景及剖面特征

    川西北大坪剖面構(gòu)造上位于四川盆地的西北部川北坳陷低緩處的龍門山斷層北段的東緣(圖1a),地理位置位于旺蒼縣雙匯鎮(zhèn)大坪村西河邊公路一側(cè)(圖1b;起點坐標30°40′33.79″ N,104°03′32.39″ E,終點坐標30°40′33.58″ N,104°03′32.85″ E),剖面頂?shù)捉缇€清楚,無植被覆蓋。受中—晚二疊世峨眉地裂運動影響,整個四川盆地及云南部分地區(qū)發(fā)育厚度巨大的玄武巖[11,27?29](圖1a),并在四川盆地表現(xiàn)為中下部鐵玄武巖、頂部酸性火山巖的火山序列[30]。晚二疊世吳家坪階在川西北稱為“吳家坪組”,為一套淺海相碳酸鹽巖夾凝灰?guī)r地層[30],自下而上可分為吳一段—吳三段三個巖性段,碳酸鹽巖中珊瑚、海百合、有孔蟲等普見,具備一定的造礁能力;四川盆地及周緣吳家坪期凝灰?guī)r在各段均有不同程度的分布,主要分布在吳一段—吳二段,呈多個薄層狀產(chǎn)出(圖1c)[31],具多個形成年齡,如賓川地區(qū)峨眉山火山序列酸性熔結(jié)凝灰?guī)r年齡為259.1±0.5 Ma[32]和259.2±0.3 Ma[33],廣元朝天剖面吳家坪組凝灰?guī)r年齡為260.1±2.8 Ma[8]等鋯石U-Pb年齡,指示了峨眉山大火成巖省的火山活動的多階段性。

    川西北大坪剖面二疊系共實測22層,其中第1~7層為茅口組,第8~22層為吳家坪組(厚45.02 m)。剖面呈近南北走向、南西西傾向、傾角40°左右,分布穩(wěn)定。巖性以泥晶灰?guī)r、生屑泥晶灰?guī)r、泥質(zhì)泥晶灰?guī)r、泥巖為主,夾薄層凝灰?guī)r,共識別出12層凝灰?guī)r,單層厚度5~54 cm(圖1d)。

    2 樣品采集及分析方法

    基于野外剖面實測,對該剖面進行了分層,描述了各層巖石類型、沉積構(gòu)造、顏色等。對凝灰?guī)r發(fā)育層位、層厚、顏色等進行了詳細描述,采集凝灰?guī)r樣品22件(DP1~22,采樣位置見圖1d),制作薄片22件,全巖及黏土礦物含量、主量元素、微量元素、掃描電鏡各20件,鋯石U-Pb年齡樣品11件。

    薄片制備執(zhí)行石油天然氣行業(yè)標準《巖石薄片制備方法》(SY/T 5913—2021);顯微成分、組構(gòu)及儲集空間觀察采用偏光顯微鏡(Zeiss AXIO SCOPE.A1);超微成分、組構(gòu)鑒定及能譜分析采用場發(fā)射環(huán)境掃描電鏡(Quanta650FEG);主量元素測試采用電感耦合等離子體發(fā)射光譜儀(PE5300V);全巖及黏土礦物測定采用X射線衍射儀(Ultima IV),電壓為40 kV,電流為40 mV,取樣品磨成粉末放在載物臺進行XRD分析,全巖礦物鑒定與含量計算采用礦物特征衍射峰值和K值法[34],黏土礦物具體步驟為取樣品磨成粉末,用蒸餾水進行脫氣,再進行沉淀和離心法分離小于2 μm的顆粒,接著對黏土礦物樣品在空氣干燥、乙二醇溶劑化24 h和在490 °C下加熱2 h后,進行三次XRD測試[35],X射線衍射數(shù)據(jù)處理使用軟件HighScore在乙二醇曲線上進行,不同物質(zhì)的X 射線吸收指數(shù)不一樣,黏土礦物吸收指數(shù)較低,黏土礦物的相對含量主要使用(001)晶面衍射峰的面積比,伊利石使用1 nm(001)晶面,蒙脫石使用1.7 nm(001)晶面,高嶺石(001)和綠泥石(002)采用0.7 nm疊加峰,計算出來的礦物相對含量誤差約為5%[34]。關(guān)鍵礦物的鑒定依據(jù)單晶礦物卡片庫與實測譜圖對比,一一對應(yīng)就能檢索出樣品的全部物相;微量元素、稀土測試均采用電感耦合等離子體質(zhì)譜儀(Aglient Technologie s7700 Series),測試方法與文獻[36]相同;以上測試在西南石油大學油氣藏地質(zhì)及開發(fā)工程國家重點實驗室和天然氣地質(zhì)四川省重點實驗室完成。

    鋯石U-Pb年齡測試過程包括鋯石樣品挑選、清洗、熱液鋯石分選、測試等步驟:首先將凝灰?guī)r樣品清洗晾干、全巖碎樣和鋯石挑選,挑選出晶形完好和色澤度較好的鋯石用環(huán)氧樹脂固定于樣品靶上,樣品靶拋光后在超純水中超聲清洗至表面可能的污染除去;其次對鋯石進行陰極發(fā)光圖像觀察及拍照,根據(jù)熱液鋯石無明顯震蕩環(huán)帶、晶形通常為半自形—他形、其棱線也不顯著等特征[37?38]挑選出熱液鋯石;最后以NIST610作為外標、91Zr作為內(nèi)標,采用激光剝蝕—電感耦合等離子體質(zhì)譜儀(ResolutionSE+Agilent7900)進行年齡測試;數(shù)據(jù)處理用Iolite程序[39],協(xié)和圖繪制采用Isoplot R程序[40],同位素比值誤差為2σ,年齡采用207Pb/235U 和206Pb/238U 比值和年齡;鋯石U-Pb年齡測試在南京宏創(chuàng)地質(zhì)勘查技術(shù)服務(wù)有限公司微區(qū)分析實驗室完成。

    3 凝灰?guī)r巖石學及分布特征

    通過野外剖面實測,將第8~9層劃分為吳一段,第10~20 層劃分為吳二段,第21~22 層劃分為吳三段,根據(jù)取樣點位置(圖1d),將DP1~6,DP11~18,DP20~22稱為吳一段凝灰?guī)r樣品,將DP7~10,DP19稱為吳二段凝灰?guī)r樣品。

    大坪剖面二疊系吳家坪組中凝灰?guī)r除第8層厚度相對較大(54 cm)外,其他層厚度一般小于20 cm,夾于灰?guī)r地層之間,展布穩(wěn)定。顏色呈灰黑色、土黃色、灰綠色等,且在同一層內(nèi)可見多種顏色變化,單層頂?shù)锥嗯c臨近地層顏色相近,可能是受多種蝕變或成巖作用影響所致(圖2a1~c1,a2~c2)。薄片、掃描電鏡及XRD測試結(jié)果表明,吳一段凝灰?guī)r樣品的礦物組成以黏土為主(69.9%~88.0%)(圖3a),同時存在數(shù)量不等的石英(3.2%~16.3%)、鉀長石(0~7.2%)、黃鐵礦(0~12.5%)、半水石膏(0~8.7%)、黃鉀鐵礬(2.0%~19.4%)、金紅石(0.6%~2.3%)、石膏(0~1.9%)(圖3c)等礦物(表1);黏土礦物以伊利石(22.9%~73.7%)、伊/蒙混層(26.3%~77.1%)為主,伊/蒙混層比中蒙皂石層和伊利石層均為10%和90%(表2);吳二段凝灰?guī)r樣品的礦物組成以黏土為主(32.9%~87.5%)(圖3b),同時存在數(shù)量不等的石英(5.3%~10.0%)、鉀長石(0~3.8%)、斜長石(0~1.8%)、方解石(0~61.4%)、半水石膏(0~2.5%)、黃鉀鐵礬(0~6.2%)、金紅石(0.4%~1.1%)等(表1);黏土礦物以伊利石(0~71.8%)、伊/蒙混層(28.2%~100%)為主,伊/蒙混層比中蒙皂石層和伊利石層均為10% 和90%(表2)。凝灰?guī)r可分為晶屑凝灰?guī)r、玻屑凝灰?guī)r及火山塵凝灰?guī)r,晶屑凝灰?guī)r中晶屑類型主要為長石、石英晶屑等,粒徑一般介于0.2~0.5 mm;玻屑灰?guī)r中玻屑多呈氣泡的浮巖狀,部分呈鐮刀狀,粒度一般介于0.1~0.8 mm;火山塵主要由粒徑小于0.05 mm玻屑組成,可含少量晶屑、玻屑。大坪剖面吳家坪組凝灰?guī)r以火山塵凝灰?guī)r為主,火山塵顯微鏡下整體為深灰色、塵點狀,可見脫?;饔玫牟痪鶆虬邎F以及形成微晶長石、石英等(圖2a3,b3,c3)。礦物學以富含伊/蒙混層和伊利石的黏土為主,含量一般大于70%;長石、石英以晶屑、孔洞充填物、脫玻化作用產(chǎn)物等存在;黃鐵礦、石膏、方解石等以溶蝕孔洞、裂縫中充填物存在(圖2a1~a3);黃鉀鐵礬系黃鐵礦氧化分解后的次生礦物,可能是暴露期形成。

    為分析研究剖面凝灰?guī)r與鄰近區(qū)域分布的共性及差異性,結(jié)合鉆井巖屑鑒定及電性曲線特征,建立了鉆井與野外剖面凝灰?guī)r分布圖(圖4)。可以看出,層位上,凝灰?guī)r均分布在吳一段、吳二段,垂向上發(fā)育6~10層,反映了凝灰?guī)r分布的時限性;厚度上,大坪剖面西側(cè)單井的凝灰?guī)r單層厚度(1~8 m)及累積厚度均較大,與凝灰?guī)r物源來自西南側(cè)相吻合,大坪剖面和鄰近的鐵爐壩剖面凝灰?guī)r單層厚度顯著減??;橫向上,凝灰?guī)r各層整體具有比較好的可對比性。

    4 凝灰?guī)r元素地球化學特征及成因

    4.1 全巖主量元素特征

    大坪剖面凝灰?guī)r全巖主量元素變化較大(表3),吳一段凝灰?guī)r燒失量介于9.010%~14.454%;SiO2含量介于46.185%~55.748%,大多數(shù)在50%以上;Al2O3的含量介于20.26%~26.51%;SiO2/Al2O3 介于1.84~2.75;TiO2 的含量0.723%~1.069%;全堿(Na2O+K2O)含量介于4.802%~6.561%;TiO2/Al2O3 介于0.028 9~0.045 6,Na2O/K2O介于0.008 7~0.056 5。吳二段凝灰?guī)r樣品燒失量介于9.521%~22.354%;SiO2含量介于32.677%~55.352%,大多數(shù)在50%以上;Al2O3的含量介于17.52%~25.47%;SiO2/Al2O3 為1.829~2.682;TiO2的含量介于0.367%~1.204%;全堿(Na2O+K2O)含量介于3.375%~6.221%;TiO2/Al2O3 介于0.019~0.047;Na2O/K2O 介于0.015 5~0.040 3。凝灰?guī)r主量元素TiO2-SiO2圖解[41](圖5)顯示,吳一段、吳二段凝灰?guī)r數(shù)據(jù)均在火成巖區(qū)域,指示大坪剖面凝灰?guī)r主要為火成巖特性。

    4.2 微量及稀土元素特征

    根據(jù)微量元素數(shù)據(jù)繪制了原始地幔標準化微量元素蛛網(wǎng)圖(圖6a)[42?43],可知大部分樣品大離子親石元素處于波谷,有明顯的Ba、Sr、Ti元素負異常。元素Ba、Sr的負虧損與斜長石有關(guān),這是由于巖漿分離結(jié)晶作用過程中斜長石的強分離結(jié)晶作用[44]。吳二段樣品DP9的Sr處于波峰,顯示Sr的正異常和Nb、Ti負異常,DP10顯示Nb、Ti負異常,可能源區(qū)受到地殼、地?;旌献饔?。

    球粒隕石標準化稀土元素配分模式圖(圖6b)[42?43]顯示,吳一段凝灰?guī)r樣品具明顯的Eu負異常(Eu/Eu*值介于0.165~0.422)、無Ce異常。這是由于斜長石能對Eu分配,在巖漿分離結(jié)晶過程中,斜長石的大量晶出將導致殘余熔體中形成明顯Eu 負異常[44],因此凝灰?guī)r中Eu 負異??梢哉J為與火山作用有關(guān)。(La/Sm)N 值介于1.17~5.58,平均為3.46;(La/Yb)N 值介于4.64~26.99,平均為12.97,(Gd/Lu)N值介于1.66~5.16,平均為3.24。稀土元素配分圖顯示吳一段所有樣品的配分模式相似,總體表現(xiàn)為LREE富集、HREE虧損。吳二段樣品顯示具Eu負異常(Eu/Eu*值介于0.41~0.77);(La/Sm)N 值介于0.57~2.45,平均為1.21,其中DP9、DP10樣品的(La/Sm)N 值均小于1;(La/Yb)N 值介于0.38~4.49,平均為1.86,DP9、DP10樣品的(La/Yb)N 值小于1,(Gd/Lu)N值介于0.95~1.37,平均為1.23。因此,吳二段DP9、DP10表現(xiàn)為LREE虧損、HREE平緩,DP8表現(xiàn)與吳一段相似的配分模式。

    4.3 凝灰?guī)r來源

    高場強元素(Nb、Ti、Zr、Ta、Hf、Sc、Th)具有抗蝕變性及Al2O3/TiO2 的比值可以用來解釋巖漿的來源[45],基性、中性、酸性巖漿Al2O3/TiO2值分別介于3~8、8~21、21~70[46]。大坪剖面吳一段凝灰?guī)r樣品Al2O3/TiO2值介于23.36~34.56(圖7a)[6,11,42,47?50],吳二段凝灰?guī)r樣品Al2O3/TiO2值介于39.70~53.04,均高于峨眉山玄武巖,范圍在酸性巖漿21~70內(nèi),與賓川酸性熔結(jié)凝灰?guī)r、峨眉山流紋巖、上寺剖面下部凝灰?guī)r相似,表明大坪剖面吳家坪組凝灰?guī)r可能與酸性火山噴發(fā)有關(guān)。凝灰?guī)r在蝕變和風化過程中,表現(xiàn)為Si的釋放和Al的累積,因此SiO2/Al2O3用以指示火山玻璃轉(zhuǎn)化為黏土礦物過程中受風化和蝕變作用影響的程度[51],大坪剖面吳一段樣品SiO2/Al2O3 含量介于1.84~2.75,吳二段樣品SiO2/Al2O3 介于1.829~2.682,明顯低于賓川酸性熔結(jié)凝灰?guī)r樣品的比值(5.47~6.31)[47],表明其受到蝕變、風化程度較高。

    Zr-Ti判別圖(圖7b)[6,11,42,47?50,52]中,大坪剖面吳一段、吳二段凝灰?guī)r樣品落在板內(nèi)巖漿內(nèi),表明其與板內(nèi)巖漿地球化學特征具有相似性;Zr/Sc-Th/Sc圖解(圖7c)[6,11,42,47?50,53]中,大坪剖面吳一段凝灰?guī)r樣品與賓川酸性熔結(jié)凝灰?guī)r、峨眉山流紋巖、上寺剖面下部凝灰?guī)r沿峨眉山大火成巖省成分變化線分布,吳二段凝灰?guī)r樣品也鄰近峨眉山大火成巖省成分變化線,暗示其與這幾者具有相同的來源;在Nb/Y-Zr/TiO2圖解(圖7d)[6,11,42,47?50,54]中,吳一段凝灰?guī)r樣品點大多落入響巖之間,說明其來源為酸性巖;吳二段凝灰?guī)r樣品落在流紋巖中,其原巖可能為中酸性流紋巖;進一步通過Rb、Sr、Nb、Ta、Nd、Th、Zr、Hf等具有源區(qū)指示意義的微量元素對吳二段樣品進行源區(qū)分析,地殼的Nb/Ta平均值為11,地幔為17.8[55?56];地殼的Rb/Sr大于0.5,地幔值小于0.05,殼?;旌蠟?.05~0.5;地殼的Nd/Th約為3[57],地幔則高于15[58];地殼的Zr/Hf為33,而地幔的為37[55?56]。吳二段凝灰?guī)r樣品中Nb/Ta介于8.59~22.04,平均值為13.95,平均值介于地殼值和地幔值;Rb/Sr介于0.04~1.93,平均值為0.92,平均值更接近地殼值;Nd/Th 介于0.17~1.30,平均值為0.57,平均值與殼源接近;Zr/Hf介于31.12~36.58,平均值為34.69,樣品值既有接近地殼值,也有接近地幔值,表明吳二段樣品為殼?;旌咸卣?。

    在Y+Nb-Rb判別圖中(圖8a),吳一段、吳二段凝灰?guī)r樣品均落在板內(nèi)區(qū)域;在Yb+Ta-Rb 判別圖中(圖8b),吳一段、吳二段凝灰?guī)r樣品也均落在板內(nèi)區(qū)域;在Y-Nb判別圖中(圖8c),吳一段凝灰?guī)r樣品落在板內(nèi)區(qū)域,吳二段落在板內(nèi)與洋脊交界處、洋脊區(qū)域;在Yb-Ta判別圖中(圖8d),吳一段凝灰?guī)r樣品落在板內(nèi)區(qū)域,吳二段凝灰?guī)r樣品落在板內(nèi)與洋脊交界處,由此顯示吳一段、吳二段凝灰?guī)r顯示出與形成于板內(nèi)伸展環(huán)境下花崗巖類巖石相似的特征。

    綜上,大坪剖面吳一段、吳二段凝灰?guī)r源巖區(qū)為中酸性火山巖,但吳一段、吳二段凝灰?guī)r的微量、稀土元素配分模式表示兩者之間有一定的差異性,暗示兩者的成因存在一定的差異。吳一段凝灰?guī)r微量元素、稀土元素均與賓川酸性熔結(jié)凝灰?guī)r、峨眉山流紋巖、上寺剖面下部凝灰?guī)r特征相似,指示了這幾者為同時同源產(chǎn)物,因此認為吳一段凝灰?guī)r可能是峨眉山大火成巖省酸性火山活動的產(chǎn)物。吳二段凝灰?guī)r樣品Zr-Ti圖中顯示與上寺剖面下部凝灰?guī)r、峨眉山大火成巖省相關(guān)巖石具有相似特征(圖7b),且樣品沿峨眉山火成巖成分變化線分布(圖7c),指示該凝灰?guī)r也可能是峨眉山大火成巖省酸性火山活動的產(chǎn)物,形成于板內(nèi)伸展環(huán)境;但其微量元素和稀土元素與吳一段凝灰?guī)r有不同的特征,微量元素表示富集Rb、Th、U,虧損Nb、Ti以及Eu弱虧損,可能受到大陸地殼的混染[60],稀土元素球粒隕石標準化配分模式呈較平坦的左傾型,與洋中脊巖石相似,且在Y-Nb、Yb-Ta判別圖(圖8c,d)中樣品投在板內(nèi)和洋中脊交界處,暗示其源區(qū)可能是DMM源的巖漿經(jīng)分離結(jié)晶作用和受地殼混染的分離結(jié)晶混染作用[61];Nb/Ta、Rb/Sr、Nd/Th、Zr/Hf等值也表明樣品主要為殼?;旌献饔玫某梢颍瑢?yīng)的峨眉山大火成巖省酸性火山作用與殼?;旌嫌嘘P(guān)。據(jù)此推測,吳家坪組吳一段、吳二段凝灰?guī)r是峨眉山大火成巖省不同期火山作用的產(chǎn)物,均形成于板內(nèi)拉張環(huán)境,但吳二段凝灰?guī)r形成可能是DMM源的巖漿受到地殼物質(zhì)的混合作用。

    5 凝灰?guī)r儲集空間類型及成巖作用

    5.1 凝灰?guī)r儲集空間類型

    凝灰?guī)r自身可存在晶間孔等原生孔隙,受成巖作用改造可形成各類溶蝕孔隙[16,18?19,62]。根據(jù)野外觀察、薄片及掃描電鏡等鑒定,在研究剖面凝灰?guī)r中識別出洞、不規(guī)則溶孔、粒內(nèi)溶孔、有機質(zhì)孔、粒間溶孔、裂縫等儲集空間類型。洞主要見于第8層,洞徑3 cm×5 cm~30 cm×47 cm,呈不規(guī)則狀,未充填或黃鐵礦、石英、石膏、鋯石等充填(圖2a1~a3、圖9a),多與裂縫相伴生,形成縫洞復合體;從洞的形態(tài)、充填物判斷其成因與埋藏期來自深部的熱液溶蝕作用相關(guān)。裂縫為流體主要運移通道,形成沿縫分布的大小不等的不規(guī)則擴溶洞,部分被黃鐵礦、石英、石膏等熱液礦物充填。不規(guī)則溶孔一般常與裂縫相伴生,溶孔多被硅質(zhì)、方解石等充填(圖2a3、圖9d,f),發(fā)育較為普遍,第8、14、16、20層凝灰?guī)r中均可見。粒內(nèi)溶孔主要為長石晶屑、伊利石、石膏等粒內(nèi)溶孔,多呈圓形—橢圓形、蜂窩形或不規(guī)則狀,孔徑細小,多呈微米級(圖9g)。有機質(zhì)孔系有機質(zhì)生烴演化過程中體積收縮或者被破壞而形成的孔隙(圖10h,i),孔隙細小,常密集分布。粒間溶孔多形成于黏土、石英長石、石膏、黃鐵礦等晶粒間,多為多邊形或者不規(guī)則狀,為納米—微米級孔隙,礦物邊緣部分具港灣狀溶蝕邊(圖9j~l)。裂縫可見構(gòu)造縫、擴溶縫,前者較細小平直,局部呈密集分布,寬一般小于100 μm(圖2a3~c3、圖9b,d,e);擴溶縫系構(gòu)造縫基礎(chǔ)上溶蝕擴大形成,相較于構(gòu)造縫張開度明顯增大,可至厘米級,縫邊緣不規(guī)則,港灣狀溶蝕邊顯著(圖9b,c,e,f);裂縫未充填或硅質(zhì)、方解石充填。

    5.2 成巖作用對儲集性的影響

    1) 溶蝕作用

    如前所述,大坪剖面凝灰?guī)r儲集空間以次生成因為主,溶蝕作用是其主要成因,凝灰?guī)r內(nèi)部礦物被溶蝕,如石膏、石英等,根據(jù)溶蝕強度的大小,可形成包括洞、不規(guī)則溶孔、粒間溶孔、粒間溶孔以及沿縫發(fā)生擴溶形成擴溶縫等不同尺度的儲集空間類型(圖9),這些孔洞、縫具有較強的儲集性。從充填物類型來看,包括黃鐵礦、石英、方解石、隱晶硅質(zhì)等,指示溶蝕作用以埋藏期為主。

    2) 熱液作用

    通過對11個樣品的10個鋯石顆粒(測試點位為9~18)的測試分析,確定凝灰?guī)r受到熱液作用的影響。稀土元素球粒隕石標準化配分模式表現(xiàn)為虧損輕稀土、重稀土異常富集的左傾型,Ce正異常和Eu負異常(圖10),Th、U值含量相對較高,U的含量介于422~3 293 μg/g,Th的含量介于25.4~115.6 μg/g,Th/U介于0.034~0.083,Th、U含量及Th/U值不能很好地指示熱液鋯石,鋯石顯示較強的Ce正異常、重稀土富集是熱液鋯石另一種不常見的熱液鋯石特征,以澳大利亞Mole 花崗巖為代表[38],與澳大利亞Boggy Plain巖體微量元素具有不同的特征[63]。

    陰極發(fā)光圖像顯示,大坪剖面凝灰?guī)r中鋯石形狀多為柱狀,半自形—他形晶,鋯石不發(fā)光,呈灰白色或黑色,鋯石一般不發(fā)育巖漿的震蕩環(huán)帶或少數(shù)模糊環(huán)帶,局部出現(xiàn)震蕩環(huán)帶可能是由于熱液改造不徹底,且棱線不明顯(圖11a),與巖漿型鋯石明顯震蕩環(huán)帶的特征存在差異,應(yīng)為熱液鋯石特征[37?38],表明其具有熱液鋯石的特性。熱液鋯石的U-Pb同位素平均年齡為132.37±0.53 Ma(圖11b),顯示凝灰?guī)r熱液流體活動為早白堊世(96~137 Ma);該期熱液在四川樂山[64]、滇東北[65]等地峨眉山玄武巖次生杏仁體及不規(guī)則溶蝕孔洞中均有響應(yīng),充填物包括黃鐵礦、銅、綠泥石等。熱液改造是大坪剖面凝灰?guī)r形成規(guī)模性儲層的重要因素,形成的儲集空間類型以不規(guī)則洞為主,洞徑3 mm×5 mm~30 mm×47 mm,常與構(gòu)造縫相伴生,充填物包括黃鐵礦、隱晶硅質(zhì)等;該期溶蝕在川西北葛底壩剖面[31]、鐵爐壩剖面[62]凝灰?guī)r中均有分布,且溶蝕特征類似,經(jīng)熱液溶蝕改造的凝灰?guī)r儲集性能顯著提升;該期熱液亦是四川盆地原油大量裂解生氣時間[66?67],對天然氣運聚成藏具有重要影響。

    3) 構(gòu)造破裂作用

    區(qū)域上多次構(gòu)造運動及其形成的如龍門山構(gòu)造帶、大巴山構(gòu)造帶等多條斷裂帶是大坪剖面凝灰?guī)r中構(gòu)造縫發(fā)育的主要作用力,裂縫多呈北東—南西向展布的特征與主構(gòu)造力基本一致,表明構(gòu)造帶應(yīng)力常對裂縫的展布具有控制作用。薄片中可見平直細小的多期構(gòu)造縫(圖9b~f)。裂縫既是重要的儲集空間類型,又能連通先期形成的孔隙,改善了孔隙連通網(wǎng)絡(luò),并為后期流體改造提供了良好的通道,整體上較大提高了凝灰?guī)r的儲滲性能。

    4) 脫?;饔?/p>

    火山玻璃受到外界環(huán)境的影響轉(zhuǎn)化為晶體,形成微晶石英、斜長石,該過程會伴隨著顆粒體積縮小,形成晶間微孔及脫?;譡68]。在埋藏前或者埋藏期,都可能發(fā)生脫?;饔?。研究區(qū)可見脫?;饔玫牟痪鶆虬邎F(圖2b3)以及脫?;饔脤е率㈩w粒間形成石英間溶孔(圖9j)。脫?;饔眯纬奢^多細小孔隙,進一步改善凝灰?guī)r的儲集性能。

    6 結(jié)論

    (1) 大坪剖面吳家坪組凝灰?guī)r主要分布在吳一段、吳二段,呈薄層狀;巖石類型以火山塵凝灰?guī)r為主,礦物以伊利石、伊/蒙混層等黏土礦物為主,存在數(shù)量不等的石英、鉀長石、斜長石、方解石、黃鐵礦、石膏、黃鉀鐵礬、金紅石等;區(qū)域上厚度向西南地區(qū)增大,發(fā)育期次與鉆井剖面具有較好的可對比性。

    (2) 大坪剖面凝灰?guī)r與峨眉山大火成巖省多期酸性火山活動有關(guān),吳一段和吳二段凝灰?guī)r均形成于板內(nèi)拉張性構(gòu)造環(huán)境,吳二段可能受到地殼地?;旌系挠绊憽?/p>

    (3) 大坪剖面凝灰?guī)r儲集空間類型包括洞、不規(guī)則溶孔、粒內(nèi)溶孔、有機質(zhì)孔、粒間溶孔、裂縫等,裂縫與不同尺度的孔、洞形成性能良好的儲滲體,儲集空間未充填或黃鐵礦、石英、方解石、隱晶硅質(zhì)等充填;其次生儲集空間受多期構(gòu)造活動及熱液流體等的綜合改造,尋找二者聯(lián)合作用區(qū)應(yīng)是川西北地區(qū)二疊系吳家坪組凝灰?guī)r儲層勘探的關(guān)鍵。

    致謝 衷心感謝評審專家和編輯提出的建設(shè)性修改意見,對本文的改進帶來了極大的幫助。

    參考文獻(References)

    [1] K?niger S, Lorenz V, Stollhofen H, et al. Origin, age and stratigraphicsignificance of distal fallout ash tuffs from the Carboniferous-Permian continental Saar-Nahe Basin (SW Germany)[J]. InternationalJournal of Earth Sciences, 2002, 91(2): 341-356.

    [2] He B, Xu Y G, Chung S L, et al. Sedimentary evidence for a rapid,kilometer-scale crustal doming prior to the eruption of theEmeishan flood basalts[J]. Earth and Planetary Science Letters,2003, 213(3/4): 391-405.

    [3] 王曼,鐘玉婷,侯瑩玲,等. 華南地區(qū)二疊紀—三疊紀界線酸性火山灰的源區(qū)與規(guī)模[J]. 巖石學報,2018,34(1):36-48.[WangMan, Zhong Yuting, Hou Yingling, et al. Source and extent of thefelsic volcanic ashes at the Permian-Triassic boundary in SouthChina[J]. Acta Petrologica Sinica, 2018, 34(1): 36-48.]

    [4] 殷鴻福,黃思驥,張克信,等. 華南二疊紀—三疊紀之交的火山活動及其對生物絕滅的影響[J]. 地質(zhì)學報,1989,63(2):169-180.[Yin Hongfu, Huang Siji, Zhang Kexin, et al. Volcanism atthe Permian-Triassic boundary in South China and its effects onmass extinction[J]. Acta Geologica Sinica, 1989, 63(2):169-180.]

    [5] 雒洋冰. 川東川南晚二疊世煤及凝灰?guī)r中微量元素地球化學研究[D]. 北京:中國礦業(yè)大學(北京),2014.[Luo Yangbing.The study of trace elements geochemistry in Late Permian coaland tuff samples from east and south Sichuan province, China[D]. Beijing: China University of Mining and Technology(Beijing),2014.]

    [6] 張晗. 四川廣元晚二疊世凝灰?guī)r地球化學特征及其地質(zhì)意義[D]. 成都:成都理工大學,2020.[Zhang Han. Geochemicalcharacteristics of the Late Permian tuffs in Guangyuan and its geologicalsignificance[D]. Chengdu: Chengdu University of Technology,2020.]

    [7] 向芳,肖倩,喻顯濤,等. 四川盆地元壩地區(qū)上二疊統(tǒng)海相凝灰沉積儲層特征[J]. 石油與天然氣地質(zhì),2022,43(4):889-901.[Xiang Fang, Xiao Qian, Yu Xiantao, et al. Reservoir characteristicsof the Upper Permian marine tuffaceous deposits in Yuanba area,Sichuan Basin[J]. Oil amp; Gas Geology, 2022, 43(4): 889-901.]

    [8] 張晗,黃虎,侯明才. 四川廣元地區(qū)朝天剖面上二疊統(tǒng)吳家坪組凝灰?guī)r成因及其地質(zhì)意義[J]. 地球科學與環(huán)境學報,2020,42(1):36-48.[Zhang Han, Huang Hu, Hou Mingcai. Origin oftuffs from Upper Permian Wujiaping Formation in Chaotian sectionof Guangyuan area, Sichuan, China and its geological significance[J]. Journal of Earth Sciences and Environment, 2020, 42(1): 36-48.]

    [9] Mahoney J J, Coffin M F. Large igneous provinces: Continental,oceanic, and planetary flood volcanism[M]. Washington: AmericanGeophysical Union, 1997: 1-438.

    [10] 田和明,代世峰,李大華,等. 重慶南川晚二疊世凝灰?guī)r的元素地球化學特征[J]. 地質(zhì)論評,2014,60(1):169-177.[TianHeming, Dai Shifeng, Li Dahua, et al. Geochemical features ofthe Late Permian tuff in Nanchuan district, Chongqing, southwesternChina[J]. Geological Review, 2014, 60(1): 169-177.]

    [11] Huang H, Cawood P A, Hou M C, et al. Silicic ash beds bracketEmeishan large igneous province tolt;1m. y. at~260Ma[J]. Lithos,2016, 264: 17-27.

    [12] Huang H, Huyskens M H, Yin Q Z, et al. Eruptive tempo ofEmeishan large igneous province, southwestern China and northernVietnam: Relations to biotic crises and paleoclimate changesaround the Guadalupian-Lopingian boundary[J]. Geology, 2022,50(9): 1083-1087.

    [13] Gao Q L, Zhang N, Xia W C, et al. Origin of volcanic ash bedsacross the Permian-Triassic boundary, Daxiakou, South China:Petrology and U-Pb age, trace elements and Hf-isotope compositionof zircon[J]. Chemical Geology, 2013, 360-361: 41-53.

    [14] Yang J H, Cawood P A, Du Y S, et al. Large igneous provinceand magmatic arc sourced Permian-Triassic volcanogenic sedimentsin China[J]. Sedimentary Geology, 2012, 261-262:120-131.

    [15] Xu L, Lin Y T, Shen W J, et al. Platinum-group elements of theMeishan Permian-Triassic boundary section: Evidence for floodbasaltic volcanism[J]. Chemical Geology, 2007, 246(1/2):55-64.

    [16] 白國帥,蔣有錄,侯帥,等. 龍鳳山地區(qū)火石嶺組火山巖優(yōu)質(zhì)儲層主控因素及成因機制[J]. 地球科學,2022,47(5):1748-1761.[Bai Guoshuai, Jiang Youlu, Hou Shuai, et al. Main controllingfactors and genetic mechanism of high-quality volcanicreservoirs in Huoshiling Formation of Longfengshan area[J].Earth Science, 2022, 47(5): 1748-1761.]

    [17] 馬劍. 馬朗凹陷條湖組含沉積有機質(zhì)凝灰?guī)r致密油成儲—成藏機理[D]. 北京:中國石油大學(北京),2016.[Ma Jian.Mechanism of reservoir formation and hydrocarbon accumulationof sedimentary organic-matter-bearing tuffaceous tight oil inthe Malang Sag, northwest China[D]. Beijing: China Universityof Petroleum (Beijing), 2016.]

    [18] 李壯. 馬朗凹陷卡拉崗組致密凝灰?guī)r儲層特征研究[D]. 青島:中國石油大學(華東),2018.[Li Zhuang. Study on the characteristicsof the tight tuff reservoir in the Kalagang Formation in theMalang Sag[D]. Qingdao: China University of Petroleum (EastChina), 2018.]

    [19] 劉小洪,王委委,馮明友,等. 準噶爾盆地克拉美麗氣田滴西14 井區(qū)石炭系蝕變凝灰?guī)r儲層熱液作用過程及時限[J]. 地球科學,2022,47(5):1694-1710.[Liu Xiaohong, Wang Weiwei,F(xiàn)eng Mingyou, et al. Hydrothermal process and duration of Carboniferousaltered tuff reservoir in well Dixi 14 area of Kelameiligas field (Junggar Basin), NW China[J]. Earth Science, 2022,47(5): 1694-1710.]

    [20] 曾琪,胡欣,屈海洲,等. 川西北部地區(qū)上二疊統(tǒng)吳家坪組的地層劃分及有利儲集巖分布[J]. 天然氣勘探與開發(fā),2020,43(4):33-47.[Zeng Qi, Hu Xin, Qu Haizhou, et al. Division ofstrata and distribution of favorable reservoir rocks of the UpperPermian Wujiaping Formation in the northwestern Sichuan Basin[J]. Natural Gas Exploration and Development, 2020, 43(4):33-47.]

    [21] 陳旋,劉俊田,龍飛,等. 三塘湖盆地二疊系凝灰?guī)r致密油勘探開發(fā)實踐及認識[J]. 中國石油勘探,2019,24(6):771-780.[Chen Xuan, Liu Juntian, Long Fei, et al. Practical experienceand understanding of exploration and development of Permiantight tuff reservoir in the Santanghu Basin[J]. China PetroleumExploration, 2019, 24(6): 771-780.]

    [22] 邱欣衛(wèi),劉池洋,毛光周,等. 鄂爾多斯盆地延長組火山灰沉積物巖石地球化學特征[J]. 地球科學:中國地質(zhì)大學學報,2011,36(1):139-150.[Qiu Xinwei, Liu Chiyang, Mao Guangzhou,et al. Petrological-geochemical characteristics of volcanicash sediments in Yanchang Formation in Ordos Basin[J]. EarthScience: Journal of China University of Geosciences, 2011, 36(1): 139-150.]

    [23] 盧浩. 鄂爾多斯盆地南部長7 段凝灰?guī)r儲層特征及形成機理研究[D]. 北京:中國石油大學(北京),2020.[Lu Hao. Studyon tuff reservoir characteristics and formation mechanism of theChang 7 member in the southern Ordos Basin[D]. Beijing:China University of Petroleum (Beijing), 2020.]

    [24] 曾琪,田云英,鄭超,等. 川西北部上二疊統(tǒng)吳家坪組沉凝灰?guī)r分布及地質(zhì)意義[J]. 天然氣技術(shù)與經(jīng)濟,2020,14(3):20-26.[Zeng Qi, Tian Yunying, Zheng Chao, et al. Distribution ofsedimentary tuff in Wujiaping Formation, northwestern SichuanBasin and Its geological significance[J]. Natural Gas Technologyand Economy, 2020, 14(3): 20-26.]

    [25] 徐敏,梁虹,鄧紹強,等. 四川盆地二疊系火山巖地震響應(yīng)特征及分布規(guī)律[J]. 科學技術(shù)與工程,2019,19(29):27-32.[Xu Min, Liang Hong, Deng Shaoqiang, et al. Seismic responsecharacteristics and distribution of Permian volcanic rocks inSichuan Basin[J]. Science Technology and Engineering, 2019,19(29): 27-32.]

    [26] 陸建林,左宗鑫,師政,等. 四川盆地西部二疊系火山作用特征與天然氣勘探潛力[J]. 天然氣工業(yè),2019,39(2):46-53.[Lu Jianlin, Zuo Zongxin, Shi Zheng, et al. Characteristics ofPermian volcanism in the western Sichuan Basin and its naturalgas exploration potential[J]. Natural Gas Industry, 2019, 39(2):46-53.]

    [27] 羅志立. 峨眉地裂運動觀對川東北大氣區(qū)發(fā)現(xiàn)的指引作用[J]. 新疆石油地質(zhì),2012,33(4):401-407.[Luo Zhili.Guidance function of Emei taphrogenesis viewpoint on discoveryof large gas province in northeastern Sichuan[J]. Xinjiang PetroleumGeology, 2012, 33(4): 401-407.]

    [28] Zhu B, Guo Z J, Zhang S N, et al. Low-Ti gabbroic pluton inDali, SW China: New evidence for back-arc lithospheric meltinginducing early-stage magmatism of the Emeishan large igneousprovince[J]. Journal of the Geological Society, 2021, 178(6):jgs2020-224.

    [29] Zhu B, Guo Z J, Zhang S N, et al. What triggered the earlystageeruption of the Emeishan large igneous province[J]. GSABulletin, 2019, 131(11/12): 1837-1856.

    [30] Xiao L, Xu Y G, Mei H J, et al. Distinct mantle sources of low-Ti and high-Ti basalts from the western Emeishan large igneousprovince, SW China: Implications for plume-lithosphere interaction[J]. Earth and Planetary Science Letters, 2004, 228(3/4):525-546.

    [31] 張云峰,黃泓森,曾琪,等. 川西北葛底壩剖面上二疊統(tǒng)吳家坪組石灰?guī)r微相[J]. 地球科學與環(huán)境學報,2022,44(2):287-297.[Zhang Yunfeng, Huang Hongsen, Zeng Qi, et al. Microfaciesof limestones from Upper Permian Wujiaping Formation atGediba section in the northwestern Sichuan, southwest China[J].Journal of Earth Sciences and Environment, 2022, 44(2):287-297.]

    [32] 鐘玉婷. 晚二疊世火山活動與生物滅絕事件的時間關(guān)系[D].北京:中國科學院大學,2013.[Zhong Yuting. Temporal relationshipbetween Late Permian volcanismsand mass extinctions[D]. Beijing: University of Chinese Academy of Sciences,2013.]

    [33] Zhong Y T, He B, Mundil R, et al. CA-TIMS zircon U-Pb datingof felsic ignimbrite from the Binchuan section: Implications forthe termination age of Emeishan large igneous province[J]. Lithos,2014, 204: 14-19.

    [34] 陳積權(quán),高遠,秦健銘,等. 松遼盆地東緣嫩江組一二段黏土礦物和主量元素地球化學特征及其古氣候意義[J]. 中國煤炭地質(zhì),2017,29(8):17-24.[Chen Jiquan, Gao Yuan, Qin Jianming,et al. Clay mineral and major element geochemical featuresand their paleoclimate significance in Nenjiang Formation1st and 2nd members, eastern margin of Songliao Basin[J]. CoalGeology of China, 2017, 29(8): 17-24.]

    [35] Gao Y, Wang C S, Liu Z F, et al. Clay mineralogy of the middleMingshui Formation (Upper Campanian to Lower Maastrichtian)from the SKIn borehole in the Songliao Basin, NE China: Implicationsfor palaeoclimate and provenance[J]. Palaeogeography,Palaeoclimatology, Palaeoecology, 2013, 385: 162-170.

    [36] 文俊,任光明,龐維華,等. 四川沐川地區(qū)中三疊統(tǒng)凝灰?guī)r鋯石U-Pb 年代學、地球化學特征及地質(zhì)意義[J]. 礦床地質(zhì),2022,41(3):585-597.[Wen Jun, Ren Guangming, Pang Weihua,et al. Zircon U-Pb geochronology, geochemical characteristicsand geological signifi-cance of tuff from Middle Triassic inMuchuan area, Sichuan province[J]. Mineral Deposits, 2022, 41(3): 585-597.]

    [37] 張丁寧,陶中一. 熱液鋯石形成機制、特征及地質(zhì)意義研究[J]. 中國錳業(yè),2017,35(3):31-33.[Zhang Dingning, TaoZhongyi. Formation amp; characteristics of hydrothermal zirconand its geological significance[J]. China’s Manganese Industry,2017, 35(3): 31-33.]

    [38] 蔡宏明,劉桂萍,展新忠. 熱液鋯石鑒定特征及在熱液型金礦床年代學研究中的應(yīng)用[J]. 中國金屬通報,2019(9):270-271.[Cai Hongming, Liu Guiping, Zhan Xinzhong. Characteristicsof hydrothermal zircon identification and its application in geochronologyof hydrothermal gold deposits[J]. China Metal Bulletin,2019(9): 270-271.]

    [39] Paton C, Woodhead J D, Hellstrom J C, et al. Improved laser ablationU-Pb zircon geochronology through robust downhole fractionationcorrection[J]. Geochemistry, Geophysics, Geosystems,2010, 11(2): Q0AA06.

    [40] Vermeesch P. IsoplotR: A free and open toolbox for geochronology[J]. Geoscience Frontiers, 2018, 9(5): 1479-1493.

    [41] Tarney J. Geochemistry of Archaean high-grade gneisses, withimplications as to the origin and evolution of the Precambriancrust[M]. London: Wiley, 1976: 405-417.

    [42] Xu Y G, Chung S L, Shao H, et al. Silicic magmas from theEmeishan large igneous province, southwest China: Petrogenesisand their link with the end-Guadalupian biological crisis[J].Lithos, 2010, 119(1/2): 47-60.

    [43] Sun S S, McDonough W F. Chemical and isotopic systematicsof oceanic basalts: Implications for mantle composition and processes[J]. Geological Society, London, Special Publications,1989, 42(1): 313-345.

    [44] 張宏飛,高山. 地球化學[M]. 北京:地質(zhì)出版社,2012:137-140. [Zhang Hongfei,Gao Shan. Geochemistry[M]. Beijing:Geological Publishing House, 2012: 137-140.]

    [45] Dai S F, Nechaev V P, Chekryzhov I Y, et al. A model for Nb-Zr-REE-Ga enrichment in Lopingian altered alkaline volcanic ashes:Key evidence of H-O isotopes[J]. Lithos, 2018, 302-303:359-369.

    [46] Hayashi K I, Fujisawa H, Holland H D, et al. Geochemistry of ~1. 9 Ga sedimentary rocks from northeastern Labrador, Canada[J]. Geochimica et Cosmochimica Acta, 1997, 61(19): 4115-4137.

    [47] Huang H, Cawood P A, Hou M C, et al. Provenance of LatePermian volcanic ash beds in South China: Implications for theage of Emeishan volcanism and its linkage to climate cooling[J].Lithos, 2018, 314-315: 293-306s.

    [48] Anh T V, Pang K N, Chung S L, et al. The Song Da magmaticsuite revisited: A petrologic, geochemical and Sr-Nd isotopicstudy on picrites, flood basalts and silicic volcanic rocks[J].Journal of Asian Earth Sciences, 2011, 42(6): 1341-1355.

    [49] Cheng L L, Wang Y, Herrin J S, et al. Origin of K-feldsparmegacrysts in rhyolites from the Emeishan large igneousprovince, southwest China[J]. Lithos, 2017, 294-295: 397-411.

    [50] Shellnutt J G, Jahn B M. Formation of the Late Permian Panzhihuaplutonic-hypabyssal-volcanic igneous complex: Implicationsfor the genesis of Fe-Ti oxide deposits and A-type granites ofSW China[J]. Earth and Planetary Science Letters, 2010, 289(3/4): 509-519.

    [51] Hong H L, Fang Q, Wang C W, et al. Clay mineralogy of alteredtephra beds and facies correlation between the Permian-Triassicboundary stratigraphic sets, Guizhou, South China[J]. AppliedClay Science, 2017, 143: 10-21.

    [52] Pearce J A. Trace element characteristics of lavas from destructiveplate boundaries[M]//Thorpe R S. Orogenic andesites andrelated rocks. Chichester: John Wiley and Sons, 1982: 528-548.

    [53] McLennan S M, Hemming S R, McDaniel D K, et al. Geochemicalapproaches to sedimentation, provenance, and tectonics[M].Boulder: Geological Society of America, 1993.

    [54] Winchester J A, Floyd P A. Geochemical discrimination of differentmagma series and their differentiation products using immobileelements[J]. Chemical Geology, 1977, 20: 325-343.

    [55] McDonough W F, Sun S S. The composition of the earth[J].Chemical Geology, 1995, 120(3/4): 223-253.

    [56] Gibbs A K. The continental crust: Its composition and evolution[J]. The Journal of Geology, 1986, 94(4): 632-633.

    [57] 姜春發(fā),王宗起,李錦軼. 中央造山帶開合構(gòu)造[M]. 北京:地質(zhì)出版社,2000:1-107.[Jiang Chunfa, Wang Zongqi, Li Jinyi.Open-close tectonics of the central orogenic belt[M]. Beijing:Geological Publishing House, 2000: 1-107.]

    [58] 肖慶輝,鄧晉福,馬大銓,等. 花崗巖研究思維與方法[M]. 北京:地質(zhì)出版社,2002:1-294.[Xiao Qinghui, Deng Jinfu, MaDaquan, et al. The ways of investigation on granitoids[M]. Beijing:Geological Publishing House, 2002: 1-294.]

    [59] Pearce J A, Harris N B W, Tindle A G. Trace element discriminationdiagrams for the tectonic interpretation of granitic rocks[J].Journal of Petrology, 1984, 25(4): 956-983.

    [60] 劉秀. 天山東段晚古生代火山巖南北對比及其大地構(gòu)造意義[D]. 北京:中國地質(zhì)大學(北京),2020.[Liu Xiu. Late Paleozoicvolcanic rocks tectonic evolution and its north-south comparisonin the eastern Tianshan[D]. Beijing: China University ofGeosciences (Beijing), 2020.]

    [61] 趙振華. 關(guān)于巖石微量元素構(gòu)造環(huán)境判別圖解使用的有關(guān)問題[J]. 大地構(gòu)造與成礦學,2007,31(1):92-103.[Zhao Zhenhua.How to use the trace element diagrams to discriminate tectonicsettings[J]. Geotectonica et Metallogenia, 2007, 31(1):92-103.]

    [62] 李一葦. 川西北地區(qū)上二疊統(tǒng)吳家坪組凝灰?guī)r特征及儲集性[D]. 成都:西南石油大學,2022.[Li Yiwei. Characteristics andreservoir property of tuffs of Upper Permian Wujiaping Formationin northwest Sichuan[D]. Chengdu: Southwest PetroleumUniversity, 2022.]

    [63] 雷瑋琰,施光海,劉迎新. 不同成因鋯石的微量元素特征研究進展[J]. 地學前緣,2013,20(4):273-284.[Lei Weiyan, ShiGuanghai, Liu Yingxin. Research progress on trace elementcharacteristics of zircons of different origins[J]. Earth ScienceFrontiers, 2013, 20(4): 273-284.]

    [64] 鄭江,郗愛華,葛玉輝,等. 峨眉山玄武巖蝕變年齡及其對揚子西緣熱事件的響應(yīng)[J]. 科學通報,2022,67(36):4415-4428.[Zheng Jiang, Xi Aihua, Ge Yuhui, et al. Alteration age of Emeishanbasalt and its response to thermal events on the westernmargin of the Yangtze Block[J]. Chinese Science Bulletin, 2022,67(36): 4415-4428.]

    [65] 朱炳泉,戴橦謨,胡耀國,等. 滇東北峨眉山玄武巖中兩階段自然銅礦化的40Ar/39Ar 與U-Th-Pb 年齡證據(jù)[J]. 地球化學,2005,34(3):235-247. [Zhu Bingquan, Dai Tongmo, HuYaoguo, et al. 40Ar/39Ar and U-Th-Pb dating for native coppermineralizations of two stages from the Emeishan flood basalts innortheastern Yunnan province, China[J]. Geochimica, 2005, 34(3): 235-247.]

    [66] 李建忠,白斌,白瑩,等. 川西北地區(qū)二疊系棲霞組超深層氣藏流體演化過程與成藏模式[J]. 石油勘探與開發(fā),2022,49(4):627-636.[Li Jianzhong, Bai Bin, Bai Ying, et al. Fluid evolutionand hydrocarbon accumulation model of sultra-deep gasreservoirs in Permian Qixia Formation of northwest Sichuan Basin,SW China[J]. Petroleum Exploration and Development,2022, 49(4): 627-636.]

    [67] 管樹巍,姜華,魯雪松,等. 四川盆地中部走滑斷裂系統(tǒng)及其控油氣作用[J]. 石油學報,2022,43(11):1542-1557.[GuanShuwei, Jiang Hua, Lu Xuesong, et al. Strike-slip fault systemand its control on oil amp; gas accumulation in central Sichuan Basin[J]. Acta Petrolei Sinica, 2022, 43(11): 1542-1557.]

    [68] 劉書強,姚孟多,馮亞琴,等. 馬朗凹陷條湖組凝灰?guī)r脫玻化孔影響因素與油氣關(guān)系[J]. 特種油氣藏,2018,25(1):16-19.[Liu Shuqiang, Yao Mengduo, Feng Yaqin, et al. Influencingfactor analysis of devitrification pore and its relationship with hydrocarbonin Tiaohu Formation of Malang Depression[J]. SpecialOil and Gas Reservoirs, 2018, 25(1): 16-19.]

    基金項目:國家自然科學基金項目(41702122);中國石油—西南石油大學創(chuàng)新聯(lián)合體科技合作項目(2020CX010301)

    国产久久久一区二区三区| 精品久久久噜噜| 国产色爽女视频免费观看| 国产免费男女视频| 成人特级黄色片久久久久久久| 国产av在哪里看| 日韩欧美三级三区| 成人无遮挡网站| 极品教师在线视频| 欧美+日韩+精品| 国内精品美女久久久久久| 精品一区二区免费观看| 亚洲av一区综合| bbb黄色大片| 成人美女网站在线观看视频| 1000部很黄的大片| 精品久久国产蜜桃| 日韩 亚洲 欧美在线| 精品一区二区三区av网在线观看| 人妻夜夜爽99麻豆av| 久久久久久伊人网av| 丝袜美腿在线中文| 亚洲第一电影网av| 日本黄色视频三级网站网址| 99久久中文字幕三级久久日本| 日韩av在线大香蕉| 国产精品伦人一区二区| 欧美精品啪啪一区二区三区| 一a级毛片在线观看| 一夜夜www| 一进一出抽搐动态| 18+在线观看网站| 日本免费一区二区三区高清不卡| 少妇的逼好多水| 国内精品宾馆在线| 男女那种视频在线观看| av天堂中文字幕网| 亚洲av免费在线观看| 日韩欧美精品免费久久| 欧美一级a爱片免费观看看| 九九在线视频观看精品| 亚洲av美国av| 久99久视频精品免费| 中文字幕久久专区| 国产三级中文精品| 超碰av人人做人人爽久久| 国内精品久久久久久久电影| 国产伦在线观看视频一区| 美女免费视频网站| 日日干狠狠操夜夜爽| 人妻久久中文字幕网| 欧美激情国产日韩精品一区| 日韩中字成人| 女人十人毛片免费观看3o分钟| 成人鲁丝片一二三区免费| 久久久久免费精品人妻一区二区| 久久精品夜夜夜夜夜久久蜜豆| 一区二区三区四区激情视频 | 亚洲av免费高清在线观看| 九色国产91popny在线| 国产高清不卡午夜福利| 精品午夜福利在线看| 给我免费播放毛片高清在线观看| 精品一区二区三区av网在线观看| 欧美又色又爽又黄视频| av视频在线观看入口| 成人鲁丝片一二三区免费| 搞女人的毛片| 欧美日本亚洲视频在线播放| 舔av片在线| 很黄的视频免费| 中文字幕高清在线视频| 免费看a级黄色片| 亚洲性夜色夜夜综合| 午夜激情福利司机影院| 亚洲欧美日韩东京热| www日本黄色视频网| 国内毛片毛片毛片毛片毛片| 久久草成人影院| 搞女人的毛片| 国产精品,欧美在线| 国内精品宾馆在线| 国产精品不卡视频一区二区| 久久精品国产自在天天线| 色播亚洲综合网| 亚洲不卡免费看| 波野结衣二区三区在线| 十八禁网站免费在线| 麻豆成人午夜福利视频| 国产伦在线观看视频一区| 91av网一区二区| 狠狠狠狠99中文字幕| 欧美极品一区二区三区四区| 一区二区三区激情视频| 久久久久久久午夜电影| 免费看av在线观看网站| 少妇高潮的动态图| 天堂动漫精品| 免费在线观看成人毛片| 亚洲国产精品合色在线| 日韩欧美一区二区三区在线观看| 在现免费观看毛片| 国产亚洲精品综合一区在线观看| eeuss影院久久| 天天躁日日操中文字幕| 国内精品久久久久精免费| 12—13女人毛片做爰片一| 九色成人免费人妻av| 无人区码免费观看不卡| 伦理电影大哥的女人| 黄色女人牲交| 国产淫片久久久久久久久| 看片在线看免费视频| 日日啪夜夜撸| 色精品久久人妻99蜜桃| 日日夜夜操网爽| ponron亚洲| 亚洲精品国产成人久久av| 欧美色欧美亚洲另类二区| 日日摸夜夜添夜夜添av毛片 | 老师上课跳d突然被开到最大视频| 亚洲成人免费电影在线观看| 狂野欧美激情性xxxx在线观看| 国产69精品久久久久777片| 搞女人的毛片| 五月玫瑰六月丁香| 韩国av在线不卡| 黄色配什么色好看| 亚洲自偷自拍三级| 亚洲欧美激情综合另类| 小说图片视频综合网站| 丰满乱子伦码专区| 内地一区二区视频在线| 国产精品久久久久久久电影| 噜噜噜噜噜久久久久久91| 一级av片app| 国产黄片美女视频| 国产一区二区亚洲精品在线观看| 国产精品久久久久久久久免| 熟女电影av网| 最近视频中文字幕2019在线8| 老司机午夜福利在线观看视频| 成年版毛片免费区| 亚洲av日韩精品久久久久久密| 亚洲一区高清亚洲精品| 亚洲电影在线观看av| 亚洲精品粉嫩美女一区| 欧美极品一区二区三区四区| 国产免费一级a男人的天堂| 床上黄色一级片| aaaaa片日本免费| 人妻制服诱惑在线中文字幕| 欧美精品国产亚洲| 极品教师在线视频| 国产免费一级a男人的天堂| 听说在线观看完整版免费高清| 91麻豆av在线| 我的老师免费观看完整版| 狂野欧美激情性xxxx在线观看| 99久国产av精品| 精品不卡国产一区二区三区| 舔av片在线| 亚洲图色成人| 亚洲 国产 在线| 一区二区三区四区激情视频 | 日本在线视频免费播放| 午夜福利在线在线| 一区二区三区免费毛片| 午夜视频国产福利| 五月伊人婷婷丁香| 国产成人av教育| 国产精品乱码一区二三区的特点| 免费不卡的大黄色大毛片视频在线观看 | 成人特级av手机在线观看| 久久99热6这里只有精品| .国产精品久久| 免费高清视频大片| 成人亚洲精品av一区二区| 在线天堂最新版资源| 舔av片在线| 黄色丝袜av网址大全| 国产午夜福利久久久久久| 国产精品自产拍在线观看55亚洲| 亚洲美女搞黄在线观看 | 国产av不卡久久| 欧美中文日本在线观看视频| 亚洲成人免费电影在线观看| av在线天堂中文字幕| 亚洲男人的天堂狠狠| 身体一侧抽搐| 亚洲国产精品成人综合色| 精品久久久久久久久久免费视频| 精品久久久久久,| 午夜福利视频1000在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产91精品成人一区二区三区| 一进一出抽搐动态| 免费观看在线日韩| 成人午夜高清在线视频| 小蜜桃在线观看免费完整版高清| 日韩欧美一区二区三区在线观看| 男人狂女人下面高潮的视频| 12—13女人毛片做爰片一| 岛国在线免费视频观看| 亚洲成av人片在线播放无| 亚洲美女视频黄频| 麻豆av噜噜一区二区三区| 三级男女做爰猛烈吃奶摸视频| 国产一区二区亚洲精品在线观看| 最近最新中文字幕大全电影3| 久久欧美精品欧美久久欧美| 国产精品久久视频播放| 日韩国内少妇激情av| 亚洲国产欧洲综合997久久,| 在线播放国产精品三级| 成人三级黄色视频| 人妻少妇偷人精品九色| 亚洲中文字幕日韩| 国产高清有码在线观看视频| 少妇的逼水好多| 色综合亚洲欧美另类图片| 亚洲成人久久爱视频| 可以在线观看毛片的网站| 在线观看av片永久免费下载| 亚洲精品国产成人久久av| 一级黄片播放器| 内射极品少妇av片p| 久久亚洲精品不卡| 高清毛片免费观看视频网站| 淫妇啪啪啪对白视频| 日本黄大片高清| 国产一区二区激情短视频| 在线免费十八禁| 国产亚洲精品av在线| 深夜a级毛片| 麻豆国产av国片精品| 看十八女毛片水多多多| 极品教师在线免费播放| 精品一区二区三区av网在线观看| 色综合亚洲欧美另类图片| 久久99热6这里只有精品| 国产极品精品免费视频能看的| 一夜夜www| 成人欧美大片| 国产精品一区二区三区四区久久| 国产精品久久电影中文字幕| 美女大奶头视频| 国产欧美日韩一区二区精品| 伦理电影大哥的女人| 国产精品爽爽va在线观看网站| 国内久久婷婷六月综合欲色啪| 国产三级在线视频| 给我免费播放毛片高清在线观看| 国产老妇女一区| 校园人妻丝袜中文字幕| 国产精品电影一区二区三区| 国产精品人妻久久久影院| 欧美另类亚洲清纯唯美| 最好的美女福利视频网| 成人性生交大片免费视频hd| 午夜福利在线观看免费完整高清在 | 国产一级毛片七仙女欲春2| 日本爱情动作片www.在线观看 | 日本免费一区二区三区高清不卡| 久久午夜亚洲精品久久| 熟女电影av网| 一边摸一边抽搐一进一小说| 国产白丝娇喘喷水9色精品| 听说在线观看完整版免费高清| 欧美一级a爱片免费观看看| 国产高清激情床上av| 在线看三级毛片| 国产免费一级a男人的天堂| 中亚洲国语对白在线视频| 国产免费男女视频| 18禁在线播放成人免费| 91在线观看av| or卡值多少钱| 国产激情偷乱视频一区二区| 波多野结衣高清作品| 日韩欧美在线二视频| 国产精品无大码| 三级毛片av免费| 亚洲中文字幕日韩| 99riav亚洲国产免费| 成人高潮视频无遮挡免费网站| 国产一区二区在线av高清观看| 国产av一区在线观看免费| 无遮挡黄片免费观看| 别揉我奶头~嗯~啊~动态视频| 婷婷精品国产亚洲av| 偷拍熟女少妇极品色| 中文字幕久久专区| 免费高清视频大片| 中文字幕久久专区| 亚洲精品色激情综合| 欧美国产日韩亚洲一区| 国内揄拍国产精品人妻在线| 精品久久久久久成人av| 午夜a级毛片| 国产精品伦人一区二区| 成人三级黄色视频| 蜜桃亚洲精品一区二区三区| 91av网一区二区| 麻豆精品久久久久久蜜桃| 色哟哟哟哟哟哟| 国产成人影院久久av| 亚洲欧美日韩东京热| 亚洲性夜色夜夜综合| 亚洲在线自拍视频| 国产午夜精品久久久久久一区二区三区 | 日韩精品中文字幕看吧| 久久精品人妻少妇| 麻豆久久精品国产亚洲av| 听说在线观看完整版免费高清| 在线播放国产精品三级| 亚洲色图av天堂| 哪里可以看免费的av片| 少妇人妻精品综合一区二区 | 免费看a级黄色片| 一区二区三区免费毛片| 亚洲三级黄色毛片| 人人妻人人澡欧美一区二区| 美女xxoo啪啪120秒动态图| 国产又黄又爽又无遮挡在线| av福利片在线观看| 狂野欧美激情性xxxx在线观看| 欧美极品一区二区三区四区| 全区人妻精品视频| 欧美性猛交黑人性爽| 欧美另类亚洲清纯唯美| 成人国产麻豆网| 久久精品国产99精品国产亚洲性色| 欧美又色又爽又黄视频| xxxwww97欧美| 最好的美女福利视频网| 成人二区视频| 午夜福利欧美成人| 蜜桃久久精品国产亚洲av| 少妇被粗大猛烈的视频| 日韩欧美在线二视频| 在线观看一区二区三区| 亚洲人成网站在线播| 18禁黄网站禁片免费观看直播| bbb黄色大片| 国内少妇人妻偷人精品xxx网站| 精品免费久久久久久久清纯| netflix在线观看网站| 国产精品三级大全| 美女被艹到高潮喷水动态| 欧美精品国产亚洲| 99久久精品一区二区三区| 波多野结衣高清作品| 亚洲在线自拍视频| 亚洲精品粉嫩美女一区| 人人妻人人看人人澡| 国产精品三级大全| 免费av观看视频| 国国产精品蜜臀av免费| 国产亚洲精品av在线| 乱人视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 嫩草影院新地址| 天天一区二区日本电影三级| av视频在线观看入口| 黄色配什么色好看| videossex国产| 精品不卡国产一区二区三区| 亚洲最大成人手机在线| 国产高清视频在线播放一区| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼| 欧美色视频一区免费| 免费在线观看影片大全网站| 99精品久久久久人妻精品| 一区二区三区高清视频在线| 国产中年淑女户外野战色| 亚洲精品在线观看二区| 久久久国产成人精品二区| 亚洲自拍偷在线| 国产精品久久久久久久久免| 国产主播在线观看一区二区| 亚洲第一电影网av| 精品午夜福利在线看| 联通29元200g的流量卡| 亚洲无线在线观看| 成人永久免费在线观看视频| 69人妻影院| 伦理电影大哥的女人| 亚洲五月天丁香| 国语自产精品视频在线第100页| 亚洲av第一区精品v没综合| 日韩,欧美,国产一区二区三区 | av.在线天堂| 午夜免费激情av| 国产亚洲精品久久久久久毛片| 欧美成人一区二区免费高清观看| 全区人妻精品视频| 久99久视频精品免费| 久久欧美精品欧美久久欧美| 91精品国产九色| 男人舔女人下体高潮全视频| 色尼玛亚洲综合影院| 国产精品一及| 国产精品电影一区二区三区| 亚洲精品国产成人久久av| 尾随美女入室| 99国产精品一区二区蜜桃av| 变态另类成人亚洲欧美熟女| 欧美色视频一区免费| 大型黄色视频在线免费观看| 国产亚洲欧美98| 欧美日韩精品成人综合77777| 九九爱精品视频在线观看| 热99在线观看视频| 老女人水多毛片| АⅤ资源中文在线天堂| 成人精品一区二区免费| 深爱激情五月婷婷| 日韩精品中文字幕看吧| 男人的好看免费观看在线视频| 国产精品一区二区性色av| 亚洲人与动物交配视频| 亚洲精品日韩av片在线观看| 亚洲美女搞黄在线观看 | 久久午夜亚洲精品久久| 此物有八面人人有两片| 国产精品伦人一区二区| 他把我摸到了高潮在线观看| 99热这里只有是精品50| 狂野欧美激情性xxxx在线观看| 内射极品少妇av片p| 夜夜爽天天搞| 国内精品美女久久久久久| 18禁黄网站禁片午夜丰满| 亚洲美女黄片视频| 内地一区二区视频在线| 欧美精品国产亚洲| 久久久久免费精品人妻一区二区| 99久久无色码亚洲精品果冻| 一区二区三区激情视频| 非洲黑人性xxxx精品又粗又长| 自拍偷自拍亚洲精品老妇| 老司机福利观看| 高清日韩中文字幕在线| 啦啦啦韩国在线观看视频| 久久久精品欧美日韩精品| 特级一级黄色大片| 最后的刺客免费高清国语| 国产免费一级a男人的天堂| 久久久久性生活片| 欧美日韩精品成人综合77777| 99久久成人亚洲精品观看| 88av欧美| 人人妻人人看人人澡| 91在线精品国自产拍蜜月| 看十八女毛片水多多多| 好男人在线观看高清免费视频| 神马国产精品三级电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| www.色视频.com| 久久久精品大字幕| 国产精品女同一区二区软件 | 好男人在线观看高清免费视频| 亚洲人成网站在线播| 免费电影在线观看免费观看| 亚洲中文日韩欧美视频| 日韩av在线大香蕉| 国产免费一级a男人的天堂| 亚洲成a人片在线一区二区| 欧美日韩精品成人综合77777| 一本精品99久久精品77| 国国产精品蜜臀av免费| 日韩精品中文字幕看吧| 好男人在线观看高清免费视频| 成人毛片a级毛片在线播放| 久久精品国产亚洲网站| 少妇丰满av| 国产精品一区二区免费欧美| 黄色视频,在线免费观看| 一个人看的www免费观看视频| 美女cb高潮喷水在线观看| 日韩欧美国产一区二区入口| 内地一区二区视频在线| 亚洲精品粉嫩美女一区| 欧美日韩瑟瑟在线播放| 嫁个100分男人电影在线观看| 国产爱豆传媒在线观看| 亚洲精品456在线播放app | 色综合婷婷激情| 69人妻影院| 亚洲第一区二区三区不卡| 麻豆精品久久久久久蜜桃| 国产伦精品一区二区三区视频9| 99精品在免费线老司机午夜| 免费大片18禁| 九九在线视频观看精品| 亚洲欧美激情综合另类| 99热网站在线观看| 午夜视频国产福利| 性插视频无遮挡在线免费观看| 色哟哟哟哟哟哟| 性欧美人与动物交配| 久久久久久伊人网av| 少妇的逼水好多| 99久久无色码亚洲精品果冻| 国内少妇人妻偷人精品xxx网站| 亚洲欧美精品综合久久99| 欧美黑人欧美精品刺激| 亚洲国产日韩欧美精品在线观看| 午夜福利在线观看吧| 国产高潮美女av| 亚洲欧美精品综合久久99| 天堂动漫精品| 亚洲欧美日韩东京热| 久久午夜亚洲精品久久| 国产 一区精品| 亚洲av成人精品一区久久| 免费大片18禁| 国产极品精品免费视频能看的| 国产亚洲欧美98| 我的女老师完整版在线观看| 午夜日韩欧美国产| 最后的刺客免费高清国语| 男人舔奶头视频| 高清在线国产一区| 欧美一区二区国产精品久久精品| 国产乱人伦免费视频| 18禁黄网站禁片午夜丰满| 精品国产三级普通话版| 国产色婷婷99| 久久久久久久精品吃奶| 日韩人妻高清精品专区| 久久久久久久久久久丰满 | 免费在线观看日本一区| 日本a在线网址| 人妻夜夜爽99麻豆av| 欧美成人免费av一区二区三区| 长腿黑丝高跟| 亚洲av五月六月丁香网| 日韩 亚洲 欧美在线| 婷婷六月久久综合丁香| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久av不卡| 日韩欧美国产一区二区入口| 亚洲天堂国产精品一区在线| 欧美一区二区亚洲| 一级黄色大片毛片| 免费在线观看成人毛片| 人妻少妇偷人精品九色| 又粗又爽又猛毛片免费看| 久久精品国产鲁丝片午夜精品 | 精品福利观看| 亚洲人成伊人成综合网2020| 色吧在线观看| 久久久久免费精品人妻一区二区| 国产亚洲av嫩草精品影院| 国产午夜精品久久久久久一区二区三区 | 亚洲自偷自拍三级| 两性午夜刺激爽爽歪歪视频在线观看| 好男人在线观看高清免费视频| 3wmmmm亚洲av在线观看| 国产伦精品一区二区三区四那| 91在线精品国自产拍蜜月| 12—13女人毛片做爰片一| 九九爱精品视频在线观看| 国产一区二区三区在线臀色熟女| 麻豆成人av在线观看| 在线看三级毛片| 美女高潮喷水抽搐中文字幕| 看免费成人av毛片| 狂野欧美白嫩少妇大欣赏| 91麻豆精品激情在线观看国产| 中文字幕av在线有码专区| 国产成人影院久久av| 成人三级黄色视频| 日日摸夜夜添夜夜添小说| 中国美女看黄片| 无遮挡黄片免费观看| 免费人成在线观看视频色| 亚洲五月天丁香| 亚洲美女黄片视频| 精品国内亚洲2022精品成人| 亚洲aⅴ乱码一区二区在线播放| 22中文网久久字幕| 午夜久久久久精精品| 国产精品精品国产色婷婷| 真实男女啪啪啪动态图| 亚洲人成网站在线播| 亚洲狠狠婷婷综合久久图片| 成人特级黄色片久久久久久久| 无人区码免费观看不卡| 国产精品久久久久久精品电影| 国产精品99久久久久久久久| 久久久色成人| 久久精品久久久久久噜噜老黄 | 国产男人的电影天堂91| 黄色一级大片看看| 国产精品美女特级片免费视频播放器| 久久精品国产清高在天天线| 成人永久免费在线观看视频| 成人性生交大片免费视频hd| 日日撸夜夜添| 亚洲国产精品合色在线| 国内精品久久久久久久电影| 国产精品国产三级国产av玫瑰| 亚洲国产精品合色在线| 久久精品国产亚洲网站| 十八禁网站免费在线| av天堂中文字幕网| 蜜桃亚洲精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 全区人妻精品视频| 亚洲男人的天堂狠狠| 免费看美女性在线毛片视频|