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

    多芳烴參數(shù)定量評價烴源巖成熟度的方法
    ——以銀額盆地中生界烴源巖為例

    2022-03-01 03:18:04陳治軍張亞雄王永昌葛宏選高志亮李子梁
    石油實驗地質 2022年1期
    關鍵詞:噻吩烷基烴源

    陳治軍,張亞雄,王永昌,王 希,葛宏選,高志亮,李子梁

    1.長江大學 資源與環(huán)境學院,武漢 430100;2.陜西延長石油(集團)有限責任公司 研究院,西安 710075;3.中國石化 石油勘探開發(fā)研究院,北京 100083;4.中國石油集團 渤海鉆探工程有限公司 定向井分公司,天津 300280

    在烴源巖成烴有效性評價、原油性質分析等方面,成熟度研究是一個重要的環(huán)節(jié)。傳統(tǒng)的成熟度評價方法有鏡質體反射率法(Ro)和巖石熱解最高峰溫法(Tmax)等[1-3]。然而傳統(tǒng)的成熟度評價方法在實際應用中存在著局限性,比如以烴源巖為基礎建立的上述方法無法應用于原油成熟度研究,再比如由于海相烴源巖普遍缺乏鏡質組而不能采用Ro法評價成熟度。

    在有機質熱演化過程中,生物標志化合物中熱穩(wěn)定性差的化合物會向熱穩(wěn)定性好的化合物轉化,使得生物標志化合物成為成熟度評價的有效指標,在烴源巖或原油成熟度研究方面有廣泛的應用[1]。相對于常用于成熟度研究的飽和烴化合物,芳烴化合物具有含量較高、熱穩(wěn)定性較好等特點,在成熟度研究方面更具優(yōu)越性[3-8]。芳烴成熟度參數(shù)建立的理論依據(jù)為:在有機質熱演化過程中,芳烴化合物的甲基化、脫甲基化、甲基重排等作用均受控于熱穩(wěn)定性,伴隨熱演化程度的增加,熱穩(wěn)定性好的化合物的相對含量逐漸增加,而熱穩(wěn)定性差的化合物的相對含量會逐漸減少,相關化合物的含量比值對成熟度有很好的指示作用[1,3]。基于上述原理,前人建立了多種芳烴成熟度參數(shù),這些參數(shù)被廣泛應用于不同地區(qū)、不同類型烴源巖成熟度的研究中,但應用效果不一[9-13]?;诜紵N參數(shù)和成熟度的相關性分析,前人建立了一些利用芳烴參數(shù)計算視鏡質體反射率(Rc)的關系式[7-10,13-14]。但這些定量評價公式都是基于單一的芳烴參數(shù),且在其他區(qū)域運用效果較差,甚至出現(xiàn)計算出的視鏡質體反射率與實際地質情況相矛盾的現(xiàn)象。如利用甲基菲指數(shù)MPI1計算Ro的定量評價關系式最早由RADKE等[7]建立,計算公式得到了廣泛應用,但王輝[6]在利用該公式計算某一地區(qū)湖相烴源巖成熟度時,某一樣品的實測鏡質體反射率(Ro)為1.47%,而計算出的Rc為2.30%或0.40%,數(shù)值相差甚遠。而實際上,一些芳烴成熟度參數(shù)可能并非僅受成熟度控制,生源、沉積環(huán)境等因素對其也有很大的影響,烴源巖的區(qū)域差異性和芳烴參數(shù)影響因素的復雜性是芳烴成熟度參數(shù)及其定量評價關系式在很多區(qū)域應用效果不佳的主要原因。越來越多的學者發(fā)現(xiàn),在運用芳烴參數(shù)定量評價烴源巖成熟度時,需要建立適合于研究區(qū)的定量評價關系式[1,5]。

    在相關地區(qū)利用芳烴參數(shù)對烴源巖成熟度開展定量評價時,發(fā)現(xiàn)利用單一芳烴成熟度參數(shù)建立的計算Rc的公式開展成熟度定量評價時往往具有較大的誤差。筆者探索出了一種全新的利用芳烴參數(shù)定量評價烴源巖成熟度的方法,即利用多元回歸法建立適合研究區(qū)的多芳烴參數(shù)定量評價成熟度的模型,該方法在一定程度上提高了芳烴參數(shù)定量評價烴源巖成熟度的準確度。

    1 芳烴成熟度參數(shù)的理論基礎

    1.1 萘系列參數(shù)

    基于一些烷基取代的萘化合物與天然產(chǎn)物結構的類似性,前人認為萘系列化合物衍生于高等植物及細菌的倍半萜至三萜類化合物[15]。萘系列化合物對有機質成熟度指示作用的理論依據(jù)為:隨著有機質熱演化程度的增加,多甲基萘系列中熱穩(wěn)定性強的異構體含量升高,熱穩(wěn)定性差的異構體含量降低,可用多甲基萘系列化合物的分布特征來研究有機質的熱演化程度[15-16]。萘由2個苯環(huán)組成,2個苯環(huán)共用2個碳原子,由于α位電子云密度高于β位,α位活性大于β位,在甲基取代反應中一般容易得到α位取代物[3,17-18](圖1)。故α位熱穩(wěn)定性比β位差,β位取代萘一旦形成就比較穩(wěn)定,β位取代萘與α位取代萘異構體的比值可以表征烴源巖和原油的成熟度[16]?;谏鲜鲈?,前人建立的烷基萘類成熟度參數(shù)有甲基萘參數(shù)(MNR)、二甲基萘參數(shù)(DNR-1—DNR-5)、三甲基萘參數(shù)(TMNr、TMNrl—TMNr3)等、四甲基萘參數(shù)(TeMNr)、五甲基萘(PMNr)等[3],并確定了上述參數(shù)的表達式(表1)。如RADKE等[7]認為隨著熱演化程度的增加,αα取代型的1,5-二甲基萘(1,5-DMN)經(jīng)甲基遷移和甲基重排重組為2,6-DMN、2,7-DMN等ββ取代型的異構體,參數(shù)DNR-1的表達式為[(2,6-DMN+2,7-DMN)/1,5-DMN],其與Ro有很好的相關性[7]。

    圖1 萘、菲和二苯并噻吩碳骨架型式據(jù)文獻[15,18]修改。

    表1 芳烴成熟度參數(shù)

    然而這些參數(shù)對于不同地區(qū)的應用效果不一,體現(xiàn)出烴源巖的區(qū)域差異性和參數(shù)影響因素的復雜性。如前人認為三甲基萘指數(shù)與成熟度具有持續(xù)的正相關關系,但劉亞洲等[4]運用這些指數(shù)對鄂爾多斯盆地三疊系延長組長7烴源巖成熟度進行評價時,發(fā)現(xiàn)了與實際情況有很多相矛盾之處,指出三甲基萘指數(shù)可能受生源、沉積環(huán)境等其他因素影響較大。因此,有學者指出烷系列化合物在甲基遷移、甲基重排等作用的過程中并不是只受熱動力學作用控制,黏土等礦物的催化作用、母質來源、生物降解等也對烷基萘成熟度參數(shù)有很大的影響[3,18-19]。因此,在運用萘系列參數(shù)評價烴源巖成熟度時,要特別注意其他因素對參數(shù)的影響。

    1.2 烷基菲系列參數(shù)

    菲系列化合物在各類沉積有機質及化石燃料中廣泛分布,且豐度較高,并在中—高成熟階段顯示出最高的豐度,因此在烴源巖或原油成熟度評價方面具有一定的優(yōu)勢,也是應用于成熟度研究較早的芳烴化合物[14,18]。菲系列化合物的碳骨架型式見圖1,由3個苯環(huán)相連而組成,β位取代基較α位取代基穩(wěn)定,隨著有機質熱演化程度的增加,α位取代的異構體的相對含量逐漸降低,β位取代的異構體的相對含量逐漸增加,β位異構體與α位異構體的比值可以表征成熟度[3,20-22]?;谶@一理論基礎,前人建立的烷基菲系列成熟度參數(shù)有甲基菲指數(shù)MPI1—MPI3、甲基菲比值MPR、甲基菲分布分數(shù)F1和F2、二甲基菲指數(shù)DPR和DPR2等[1-2,5,9](表1)。

    這些參數(shù)被廣泛應用于不同地區(qū)烴源巖或原油的成熟度研究中,但在應用中也存在部分參數(shù)并非主要受成熟度控制、對某些參數(shù)隨成熟度的變化規(guī)律存較大的爭議、一些參數(shù)的定量評價應用效果不佳等問題,這可能由烴源巖存在區(qū)域差異性、參數(shù)影響因素存在復雜性、化驗分析存在誤差性等原因所致。如前人[4]利用KVALHEIM等[9]建立的關系式計算出的某一地區(qū)的烴源巖樣品的甲基菲分布分數(shù)F2為0.22~0.28,按照已有的評價標準[23](低熟階段F2<0.27,成熟階段F2=0.27~0.35,高熟階段F2>0.35),這些烴源巖以低成熟烴源巖為主,而實際上這些烴源巖為成熟—高成熟烴源巖,這種矛盾表明對于研究區(qū)甲基菲比值可能受到成熟度以外的其他因素的影響較大?;诖耍叭苏J為應用烷基菲參數(shù)定量評價成熟度,需建立適合于研究區(qū)的關系式[1,6]。

    1.3 烷基二苯并噻吩系列參數(shù)

    二苯并噻吩類化合物廣泛存在于沉積有機質中,由于分子結構對稱(圖1),該類化合物具有較高的穩(wěn)定性[4,11]。二苯并噻吩系列化合物在有機質熱演化過程中也存在著甲基化和甲基重排作用,同樣由于熱穩(wěn)定性的差異,一些熱穩(wěn)定性強的化合物的相對含量越來越大,而一些熱穩(wěn)定性差的化合物的相對含量越來越小,相關化合物的比值可以指示成熟度。這些參數(shù)在不同地區(qū)的應用中有很好的效果,且這些參數(shù)在較廣的成熟度范圍內對成熟度均有較好的響應,在高成熟度下更具持續(xù)、穩(wěn)定的優(yōu)勢[11]?,F(xiàn)有研究已建立了很多烷基二苯并噻吩成熟度參數(shù),如甲基二苯并噻吩參數(shù)MDR[12-13]、甲基二苯并噻吩分布指數(shù)MDBI[14]、二甲基二苯并噻吩參數(shù)DMDBT1、DMDBT2、DMDBT3、DMDBT4、DMDBT5[10-12]等(表1)。如HUGHES[13]認為4-甲基二苯并噻吩(4-MDBT)的熱穩(wěn)定性好于1-甲基二苯并噻吩(1-MDBT),建立了甲基二苯并噻吩參數(shù)MDR,認為隨熱演化程度增加該值應該逐漸增大。魏志彬等[14]提出了甲基二苯并噻吩分布指數(shù)MDBI,并驗證了Ro在0.67%~1.22%范圍內,MDBI與Ro有很好的正相關關系。CHAKHMAKHCHEV等[10]研究發(fā)現(xiàn)2,4-二甲基二苯并噻吩(2,4-DMDBT)和4,6-二甲基二苯并噻吩(4,6-DMDBT)的相對含量隨著熱演化程度的增加而增大,而1,4-二甲基二苯并噻吩(1,4-DMDBT)相對含量隨著熱演化程度的增加而減小,提出了一系列二甲基二苯并噻吩成熟度參數(shù)。

    盡管烷基二苯并噻吩系列成熟度參數(shù)被廣泛應用,但他們的影響因素可能更為復雜。周葉駿等[11]的研究表明,二苯并噻吩類的烷基化程度在咸水湖相、海相、淡水湖相中呈現(xiàn)較為明顯的分異,認為鹽度差異相對于海陸原始有機物組成差異是主導二苯并噻吩類化合物烷基化程度的更為重要的因素。前人建立了一些利用烷基二苯并噻吩系列成熟度參數(shù)計算Rc的關系式,但基本都是針對于中—低成熟度烴源巖,對于高—過成熟度烴源巖還處于探索階段[1]。這主要是因為一些烷基二苯并噻吩系列成熟度參數(shù)與烴源巖實測Ro的關系更為復雜,存在變化不明顯、甚至反轉等現(xiàn)象,使得建立定量評價關系式變得較為困難[11,14]。

    2 方法介紹

    由于芳烴化合物成熟度參數(shù)的影響因素比理論中更為復雜,且烴源巖存在著區(qū)域差異性,使得利用單一參數(shù)定量評價烴源巖成熟度的做法會出現(xiàn)較大誤差。這就需要我們在實際研究中做到以下兩點:一是針對不同地區(qū)、不同類型的烴源巖,套用他區(qū)已建立的關系式可能會達不到理想的應用效果,需要建立適合研究區(qū)的芳烴參數(shù)定量評價關系式;二是在利用芳烴成熟度參數(shù)定量評價成熟度時,要盡量避免利用單一參數(shù),采用多參數(shù)綜合分析可能會減少誤差。

    基于上述考慮,筆者在相關研究中探索出了一種提高芳烴參數(shù)定量評價烴源巖成熟度準確度的方法。該方法基于研究區(qū)一定數(shù)量的烴源巖芳烴色譜—質譜(GC-MS)分析數(shù)據(jù),通過開展各芳烴成熟度參數(shù)與Ro的相關性分析,分別在3大類芳烴成熟度參數(shù)(烷基萘系列、烷基菲系列和烷基二苯并噻吩系列)中找出1個對成熟度最敏感(相關系數(shù)最大)的參數(shù),利用數(shù)值分析軟件對這幾個參數(shù)和Ro進行多元線性回歸,擬合出一個利用多參數(shù)定量計算Rc的關系式。

    該方法的具體實施步驟是:(1)選取研究目標,對具有代表性的烴源巖進行取樣。研究目標區(qū)應為同一盆地或同一凹陷。目標層系選擇參照沉積體系,如可將海相沉積體系、湖相沉積體系和海陸過渡相沉積體系中的1個沉積體系所包含的烴源巖作為研究目標層系。烴源巖樣品具有代表性體現(xiàn)在以下3個方面:一是樣品的數(shù)量多,為開展相關性分析,樣品數(shù)量要求大于3個,當然樣品數(shù)量越多越好;二是成熟度范圍大,最好涵蓋未成熟、低成熟、成熟、高成熟、過成熟等所有熱演化階段;三是樣品在各演化階段均有較為均勻的分布。在獲取烴源巖樣品后,對所有烴源巖樣品開展干酪根Ro測定,對瀝青抽提物的芳烴餾分開展GC-MS分析,也可收集研究目標區(qū)前期的相關資料。(2)計算每個樣品的芳烴成熟度參數(shù),參數(shù)的名稱和計算公式參照表1。開展各芳烴成熟度參數(shù)與Ro的相關性分析,分別在烷基萘類、烷基菲類和烷基二苯并噻吩類3大類芳烴成熟度參數(shù)中找出1個對成熟度最敏感的參數(shù),對成熟度最敏感體現(xiàn)在該參數(shù)在該系列所有參數(shù)中與Ro的相關系數(shù)最大。運用數(shù)值分析軟件(如SPSS等),對所確定的3個成熟度參數(shù)和Ro進行多元回歸,擬合出一個利用多參數(shù)定量計算Rc的多元關系式。(3)對建立的多元關系式開展有效性評價,即利用所建立的模型計算所有樣品的Rc,對Rc和Ro開展相關性分析。如果相關系數(shù)高(如相關系數(shù)R2>0.9)、且相比于單一芳烴成熟度參數(shù)與Ro的相關系數(shù)有大幅度提高,說明定量評價的準確度有明顯的提高,建立的模型適合在研究區(qū)推廣或作為相關研究的參考。如果相關系數(shù)不高(如相關系數(shù)R2<0.9),但相比于單一芳烴成熟度參數(shù)與Ro的相關系數(shù)有大幅度提高,說明選取的目標研究區(qū)不合理(即選擇的區(qū)域太大、或目標層系所包含的沉積體系過多),需要調整縮小研究目標,再次嘗試該方法,直到定量評價的準確度有明顯的提高。如果相關系數(shù)較低(如R2<<0.9)、且相比于單一芳烴成熟度參數(shù)與Ro的相關系數(shù)也沒有大幅度提高,說明芳烴成熟度參數(shù)在研究區(qū)的影響因素較為復雜,不適合開展芳烴成熟度參數(shù)定量評價烴源巖成熟度。

    該方法提供了一種全新的利用芳烴參數(shù)定量評價烴源巖成熟度的研究思路,在很大程度上提高了芳烴參數(shù)定量評價烴源巖成熟度的準確度,研究目標區(qū)所建立的評價模型對該區(qū)相關研究有很好的參考意義。同時,該方法具備很好的可操作性,對于各種勘探研究區(qū)塊均值得推廣。

    3 應用實例

    3.1 烴源巖地球化學特征

    以銀額盆地哈日凹陷中生界湖相烴源巖為例(圖2),驗證該方法的應用效果。哈日凹陷中生代為斷陷湖盆的全面發(fā)展期,主要發(fā)育湖泊沉積體系[24-27]。中生代的沉積地層主要有下白堊統(tǒng)和上白堊統(tǒng),下白堊統(tǒng)自下而上有巴音戈壁組、蘇紅圖組和銀根組,上白堊統(tǒng)僅有烏蘭蘇海組[27-30]。該區(qū)中生界烴源巖較為發(fā)育,烴源巖巖性主要有深灰色—黑色含灰泥巖、灰色—深灰色泥巖、灰色灰質泥巖、灰色白云質泥巖、灰色泥質白云巖等[27-28,31]。

    圖2 銀額盆地哈日凹陷構造位置及相關取樣井

    本次研究的烴源巖樣品來自于H1井—H6井(圖2),數(shù)量為69個,樣品巖性為灰質泥巖、含灰泥巖、含云泥巖等。對所有烴源巖樣品開展有機地球化學測試分析,由長江大學地球化學實驗室完成;對所有烴源巖樣品的瀝青抽提物的芳烴餾分開展色譜—質譜(GC-MS)分析,相關測試分析也是由長江大學地球化學實驗室完成,分析儀器及流程見文獻[27]。

    結果表明,樣品的總有機碳(TOC)含量為0.05%~6.91%,平均為1.47%;熱解生烴潛量(S1+S2)為0.03~77.87 mg/g,平均為6.40 mg/g。從TOC的分布頻率來看(圖3a),有機質豐度較低的非烴源巖(TOC小于0.4%,評價標準據(jù)文獻[32-33],下同)占總樣品的33.3%,差烴源巖(TOC為0.4%~0.6%)占8.7%,中等有機質豐度烴源巖(TOC為0.6%~1.0%)占17.4%,好的烴源巖(TOC為1.0%~2.0%)占15.9%,極好的烴源巖(TOC大于2.0%)占24.6%。從S1+S2的分布頻率同樣也可以看出本次研究樣品有機質豐度的多樣性,部分烴源巖的生烴潛力較大,但較小生烴潛力的烴源巖也有很大的占比(圖3b)。

    圖3 銀額盆地哈日凹陷中生界湖相烴源巖樣品的TOC(a)、S1+S2(b)分布頻率直方圖

    巖石熱解參數(shù)氫指數(shù)(IH)、氧指數(shù)(IO)、最高熱解峰溫(Tmax)對烴源巖的有機質類型有很好的指示作用[27,32-33],分類結果顯示烴源巖有機質類型主要為混合型(Ⅱ1—Ⅱ2型),也有少量有機質類型較好的Ⅰ型和較差的Ⅲ型存在,樣品具有典型的湖相烴源巖特征(圖4)。樣品的Ro為0.52%~2.35%,平均為0.95%。其中低成熟烴源巖(Ro為0.50%~0.70%,評價標準據(jù)文獻[32-33],下同)樣品占總樣品的36%,成熟烴源巖(Ro為0.70%~1.30%)樣品占總樣品的47%,高成熟烴源巖(Ro為1.30%~2.00%)樣品占總樣品的9%,過成熟烴源巖(Ro大于2.00%)樣品占總樣品的8%。

    圖4 銀額盆地哈日凹陷中生界湖相烴源巖樣品有機質類型劃分

    總的來說,本次研究的烴源巖樣品數(shù)量多,有機質豐度具有多樣性,有機質類型豐富,有機質成熟度跨度大,且烴源巖在低成熟、成熟、高成熟、過成熟等熱演化階段有較為均勻的分布,具有很好的代表性,適合運用該方法進行研究。

    3.2 芳烴宏觀組成

    前人研究表明,烴源巖生源特征、沉積環(huán)境、成熟度等均會對芳烴的組成產(chǎn)生影響。如1,2,5-三甲基萘和1,2,5,6-四甲基萘可能來源于高等植物生源的五環(huán)三萜香樹素或樹脂生源的二環(huán)二萜刺柏酸,較高的1,2,5-三甲基萘和1,2,5,6-四甲基萘含量預示著較高的陸源高等植物輸入[5,33-34]。“三芴”系列化合物的相對組成可以指示烴源巖和原油的沉積環(huán)境,陸相淡水或微咸水湖相下的烴源巖/原油中芴的含量較高,鹽湖相和海相碳酸鹽巖烴源巖/原油中硫芴(二苯并噻吩)含量較高,在沼澤相煤或煤成油中氧芴(二苯并呋喃)含量較高[35-36]。芳烴化合物的宏觀組成對成熟度也有一定的指示作用,未成熟—低成熟度原油,芳烴化合物中四環(huán)、五環(huán)化合物含量較高;中—高成熟度的原油,芳烴化合物中二環(huán)或三環(huán)化合物占優(yōu)勢[5,21]。

    圖5 銀額盆地哈日凹陷中生界湖相烴源巖樣品芳烴系列相對含量

    3.3 芳烴參數(shù)與Ro相關性分析

    前人建立的烷基萘成熟度參數(shù)多達12個,由于本次GC-MS分析未能檢出1,8-二甲基萘(1,8-DMN),相關參數(shù)DNR-2、DNR-3、DNR-4和DNR-5未能求取,繪制出其他8個烷基萘類成熟度參數(shù)與Ro的交會圖(圖6a-h),并分別開展線性回歸和相關性分析??梢钥闯?,烷基萘成熟度參數(shù)整體與Ro的相關性差別很大。甲基萘參數(shù)MNR、二甲基萘比值DNR-1、三甲基萘指數(shù)TMNr3和五甲基萘指數(shù)PMNr在本次研究的成熟度范圍內(Ro為0.52%~2.35%)與Ro相關性較差或者幾乎無相關性,相關系數(shù)R2僅為0.000 4~0.110 1(圖6a,b,f,h),表明對于湖相烴源巖,這些參數(shù)可能受沉積環(huán)境、生源等其他因素影響較大。三甲基萘指數(shù)TMNr、TMNr1、TMNr2和四甲基萘指數(shù)TeMNr等參數(shù)雖然在中低成熟度階段(Ro<1.00%)值較為分散,但整體與Ro相關性較好,R2為0.281 1~0.550 3(圖6c-e,g)。烷基萘成熟度參數(shù)中參數(shù)TMNr的R2最大,說明研究區(qū)烷基萘類成熟度參數(shù)中TMNr為對Ro最敏感的參數(shù)。

    同樣,繪制出8個烷基菲成熟度參數(shù)與Ro的交會圖(圖6i-p),并開展線性回歸分析??梢钥闯觯瑓?shù)DPR、DPR2與Ro的相關性較差,因為它們與Ro的相關系數(shù)R2分別僅為0.068 6和0.013 5(圖6o,p)。除了參數(shù)DPR和DPR2,其他參數(shù)與Ro的相關性顯著,R2為0.674 4~0.923 7,表明烷基菲成熟度參數(shù)對于湖相烴源巖有很好的適用性。其中,參數(shù)MPR的R2最大,為0.923 7,預示著研究區(qū)烷基菲類成熟度參數(shù)中MPR對Ro最敏感。

    圖6 銀額盆地哈日凹陷中生界湖相烴源巖芳烴參數(shù)與Ro相關關系

    烷基二苯并噻吩類7個成熟度參數(shù)與Ro的線性回歸分析結果表明,除了參數(shù)MDR、MDBI與Ro的相關性不顯著外(R2分別為0.089 9和0.213 2),其他烷基二苯并噻吩成熟度參數(shù)與Ro的相關性均較好(圖5q-u),R2為0.808 5~0.904 5,其中參數(shù)DMDBT1與Ro的相關系數(shù)R2最大,為0.9045,說明研究區(qū)烷基二苯并噻吩類成熟度參數(shù)中DMDBT1對Ro最敏感。

    3.4 模型的建立與應用效果分析

    基于以上芳烴參數(shù)與Ro的相關性分析,確定的研究區(qū)烷基萘類、烷基菲類和烷基二苯并噻吩類芳烴成熟度參數(shù)中對成熟度最為敏感的參數(shù)分別為TMNr、MPR和DMDBT1。對TMNr、MPR和DMDBT1和Ro進行多元線性回歸,用SPSS軟件擬合出的利用這3個參數(shù)計算視鏡質體反射率(Rc)的關系式:

    Rc=0.316 8×TMNr-0.291 3×MPR+

    0.167 3×DMDBT1+0.395 6

    (1)

    利用公式(1)計算所有樣品的Rc,并對Rc和Ro開展相關性分析,新模型計算出的Rc和Ro之間的相關系數(shù)R2高達0.96(圖7)。

    圖7 銀額盆地哈日凹陷中生界湖相 烴源巖樣品新模型計算的Rc與Ro相關關系

    本應用實例所有樣品的23個芳烴成熟度參數(shù)與Ro的相關系數(shù)R2為0.000 4~0.923 7。雖然有極少數(shù)單一參數(shù)與Ro的相關系數(shù)較高,如MPR與Ro的R2為0.923 7(為最高),DMDBT1與Ro的R2為0.904 5,這幾個參數(shù)在一定程度上能夠滿足烴源巖成熟度評價的精度要求,但地質定量評價研究中的精度越高往往會得到更為客觀的結論。且絕大多數(shù)單一參數(shù)與Ro的相關系數(shù)均較小,如所有單一參數(shù)中與Ro的相關系數(shù)R2小于0.8的占81%。相比于單一芳烴參數(shù)與Ro的相關系數(shù),新模型計算出的Rc與Ro的相關系數(shù)有了大幅度提高。如分別選取樣品中成熟度較低和較高的2個樣品,實測Ro分別為0.56%和2.31%,利用單一芳烴參數(shù)中與Ro的相關系性最好的MPR參數(shù)建立的模型計算出的Rc分別為0.66%和2.14%,而新模型計算出的Rc分別為0.54%和2.23%,可見新模型計算出的Rc精度明顯較高。

    以上表明利用多芳烴參數(shù)建立定量計算Rc的方法具有明顯的優(yōu)越性,體現(xiàn)在相對于單因素定量評價的準確度能有大幅度地提高。本研究應用實例的研究對象為湖相烴源巖,代表了最重要的、最廣泛的一類烴源巖類型,用本方法建立的公式(1)對于銀額盆地中生界烴源巖具有很好的適用性,對于其他區(qū)域的湖相烴源巖的相關研究也能提供參考。

    4 結論

    (1)基于熱穩(wěn)定性的差異性,前人建立了30余個芳烴成熟度參數(shù),但實際上一些參數(shù)可能并非僅受成熟度控制,生源特征、沉積環(huán)境等因素對其也會有很大影響,導致一些參數(shù)在某些研究區(qū)的應用效果并不理想,且前人基于單一參數(shù)建立的計算Rc的模型在其他地區(qū)也并不適用。

    (2)本研究探索出了一種提高芳烴參數(shù)定量評價烴源巖成熟度準確度的方法。該方法基于研究區(qū)一定數(shù)量的烴源巖芳烴色譜—質譜(GC-MS)分析數(shù)據(jù),通過開展各芳烴成熟度參數(shù)與Ro的相關性分析,分別在3大類芳烴成熟度參數(shù)(烷基萘系列、烷基菲系列和烷基二苯并噻吩系列)中找出1個對成熟度最敏感(相關系數(shù)最大)的參數(shù),利用數(shù)值分析軟件對這幾個參數(shù)和Ro進行多元線性回歸,擬合出一個利用多參數(shù)定量計算Rc的關系式。

    (3)本方法在銀額盆地哈日凹陷中生界湖相烴源巖研究中取得了很好的應用效果,用本方法擬合出的關系式計算所有樣品的Rc,結果Rc和Ro相關系數(shù)R2高達0.96,相對于單一芳烴參數(shù),相關系數(shù)有了大幅度的提高,表明該方法具有較好的適用性和較強的可操作性,在相關研究中值得推廣。

    猜你喜歡
    噻吩烷基烴源
    烷基胺插層蒙脫土的陽離子交換容量研究
    東濮凹陷西南部晚古生代—早新生代烴源研究
    探討醫(yī)藥中間體合成中噻吩的應用
    4,7-二噻吩-[2,1,3]苯并硒二唑的合成及其光電性能
    合成化學(2015年10期)2016-01-17 08:56:47
    中國近海湖相優(yōu)質烴源巖形成的主要控制因素
    直接合成法制備載銀稻殼活性炭及其對苯并噻吩的吸附
    應用化工(2014年9期)2014-08-10 14:05:08
    柴達木盆地柴北緣侏羅系烴源巖特征分析
    塔里木盆地下古生界主力烴源巖分布
    五種小麥麩皮烷基酚類化合物體外抗腫瘤作用及初步的機制研究
    噻吩類硫化物在FAU分子篩上吸附的分子模擬研究
    国产成人啪精品午夜网站| 亚洲综合精品二区| 午夜精品国产一区二区电影| 人妻一区二区av| 国产精品 国内视频| 最新在线观看一区二区三区 | 99久久99久久久精品蜜桃| 免费看av在线观看网站| 国产成人午夜福利电影在线观看| 免费高清在线观看日韩| 免费人妻精品一区二区三区视频| 免费久久久久久久精品成人欧美视频| 国产精品久久久久久精品电影小说| 精品亚洲乱码少妇综合久久| 久久国产精品大桥未久av| 欧美人与性动交α欧美精品济南到| 久久久久人妻精品一区果冻| 丰满饥渴人妻一区二区三| 久热这里只有精品99| 国产男女超爽视频在线观看| 久久97久久精品| 一级黄片播放器| 国产免费又黄又爽又色| 在线天堂最新版资源| 少妇人妻 视频| 最近中文字幕高清免费大全6| 国产黄色视频一区二区在线观看| 免费少妇av软件| 最近中文字幕2019免费版| 日韩精品有码人妻一区| 久久国产亚洲av麻豆专区| 亚洲成人一二三区av| a级片在线免费高清观看视频| 久久久久久久久久久久大奶| 美女扒开内裤让男人捅视频| 免费黄色在线免费观看| 宅男免费午夜| 成人黄色视频免费在线看| 欧美久久黑人一区二区| 日韩中文字幕欧美一区二区 | 精品福利永久在线观看| 欧美国产精品一级二级三级| 欧美日韩av久久| 欧美精品高潮呻吟av久久| 最近中文字幕2019免费版| 99久久精品国产亚洲精品| 国产免费视频播放在线视频| 亚洲精品国产av蜜桃| 国产一区二区三区综合在线观看| av福利片在线| 丝袜美足系列| 丝瓜视频免费看黄片| 黄色一级大片看看| 久久精品久久久久久噜噜老黄| 我的亚洲天堂| 人人妻,人人澡人人爽秒播 | 老鸭窝网址在线观看| 午夜激情久久久久久久| 美女国产高潮福利片在线看| 一区二区日韩欧美中文字幕| 久久免费观看电影| 欧美精品av麻豆av| 黄色 视频免费看| 国产精品二区激情视频| 观看av在线不卡| 观看美女的网站| 热re99久久精品国产66热6| 少妇被粗大的猛进出69影院| 香蕉国产在线看| 国产成人av激情在线播放| 一边摸一边抽搐一进一出视频| 色婷婷久久久亚洲欧美| 成人午夜精彩视频在线观看| 欧美亚洲 丝袜 人妻 在线| 国产亚洲最大av| 97精品久久久久久久久久精品| 国产又色又爽无遮挡免| 在线观看免费日韩欧美大片| 亚洲情色 制服丝袜| 韩国高清视频一区二区三区| 欧美激情极品国产一区二区三区| 操美女的视频在线观看| 99久国产av精品国产电影| 亚洲欧美色中文字幕在线| 亚洲第一区二区三区不卡| 国产成人欧美在线观看 | 天天影视国产精品| 婷婷色麻豆天堂久久| av一本久久久久| 国产精品久久久久成人av| 女人高潮潮喷娇喘18禁视频| 免费观看av网站的网址| 成人影院久久| 丝袜在线中文字幕| 欧美日韩综合久久久久久| 七月丁香在线播放| 国产精品一区二区精品视频观看| 宅男免费午夜| 国产熟女午夜一区二区三区| 丝袜美足系列| 久久久精品免费免费高清| 视频在线观看一区二区三区| 久久青草综合色| 亚洲国产欧美日韩在线播放| 最近中文字幕2019免费版| 日韩中文字幕视频在线看片| 国产亚洲av高清不卡| 久久av网站| 国产 一区精品| 日日爽夜夜爽网站| 一级爰片在线观看| 日韩大码丰满熟妇| 少妇猛男粗大的猛烈进出视频| 女的被弄到高潮叫床怎么办| 人体艺术视频欧美日本| www.精华液| av女优亚洲男人天堂| 欧美变态另类bdsm刘玥| 久久久久精品久久久久真实原创| 99久国产av精品国产电影| 久热这里只有精品99| 黄片无遮挡物在线观看| 黄色毛片三级朝国网站| 天天影视国产精品| 制服诱惑二区| av卡一久久| 国产野战对白在线观看| 美女大奶头黄色视频| www.精华液| 国产精品一二三区在线看| 久久鲁丝午夜福利片| 精品国产一区二区久久| 久久久久网色| 久久久久精品人妻al黑| 在线观看三级黄色| 妹子高潮喷水视频| 九草在线视频观看| 亚洲综合色网址| 一二三四在线观看免费中文在| 欧美激情 高清一区二区三区| 亚洲天堂av无毛| 亚洲婷婷狠狠爱综合网| 热99久久久久精品小说推荐| 国产欧美日韩综合在线一区二区| 免费在线观看视频国产中文字幕亚洲 | 日本wwww免费看| 无限看片的www在线观看| 男女无遮挡免费网站观看| 婷婷色麻豆天堂久久| 久久久久久久久免费视频了| av在线老鸭窝| 一级片免费观看大全| 男女高潮啪啪啪动态图| 韩国精品一区二区三区| 日韩一本色道免费dvd| 欧美精品av麻豆av| 国产亚洲av片在线观看秒播厂| 欧美精品一区二区大全| 国产成人a∨麻豆精品| 在线观看人妻少妇| 免费高清在线观看视频在线观看| 午夜福利影视在线免费观看| 在线精品无人区一区二区三| 亚洲一级一片aⅴ在线观看| 丝袜脚勾引网站| 人妻一区二区av| 日韩一卡2卡3卡4卡2021年| 女性被躁到高潮视频| 欧美黄色片欧美黄色片| 日韩中文字幕视频在线看片| 国产精品蜜桃在线观看| 亚洲欧美一区二区三区黑人| 中文字幕av电影在线播放| 久久久国产欧美日韩av| 十八禁人妻一区二区| 啦啦啦在线免费观看视频4| 十八禁高潮呻吟视频| 中文天堂在线官网| netflix在线观看网站| 美女中出高潮动态图| 免费人妻精品一区二区三区视频| 咕卡用的链子| 欧美日本中文国产一区发布| 黑人猛操日本美女一级片| 美女视频免费永久观看网站| 免费观看av网站的网址| 中文字幕av电影在线播放| 欧美97在线视频| av在线观看视频网站免费| 一二三四中文在线观看免费高清| 国产av一区二区精品久久| 欧美日韩福利视频一区二区| 又黄又粗又硬又大视频| 国产一区二区三区综合在线观看| 欧美少妇被猛烈插入视频| 国产精品蜜桃在线观看| 国产成人精品无人区| 91老司机精品| 少妇精品久久久久久久| 免费人妻精品一区二区三区视频| 美国免费a级毛片| 赤兔流量卡办理| 精品久久久久久电影网| 老汉色∧v一级毛片| 夫妻性生交免费视频一级片| 18禁裸乳无遮挡动漫免费视频| 国产 精品1| 亚洲,欧美,日韩| 一本色道久久久久久精品综合| 啦啦啦在线免费观看视频4| 久久精品国产a三级三级三级| 久久国产精品大桥未久av| 婷婷色综合大香蕉| 国产精品99久久99久久久不卡 | 国产精品嫩草影院av在线观看| 亚洲,欧美精品.| 免费不卡黄色视频| 国产一级毛片在线| 国产黄色免费在线视频| 不卡av一区二区三区| 亚洲在久久综合| xxx大片免费视频| 下体分泌物呈黄色| 国产高清不卡午夜福利| 成人18禁高潮啪啪吃奶动态图| 伊人久久国产一区二区| 国产熟女午夜一区二区三区| 国产亚洲av高清不卡| 欧美日韩亚洲综合一区二区三区_| 黑人猛操日本美女一级片| 99久久人妻综合| 十八禁人妻一区二区| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美成人精品一区二区| 日韩av免费高清视频| 久久久久久人人人人人| 亚洲精品视频女| 亚洲婷婷狠狠爱综合网| 尾随美女入室| 国产成人免费观看mmmm| 999久久久国产精品视频| 国产精品久久久久久人妻精品电影 | 97在线人人人人妻| 欧美日韩成人在线一区二区| 国产老妇伦熟女老妇高清| 黄色一级大片看看| 国产片内射在线| 99久久人妻综合| 日韩大码丰满熟妇| 欧美精品亚洲一区二区| 国产熟女欧美一区二区| 少妇猛男粗大的猛烈进出视频| 久久久精品免费免费高清| 99热网站在线观看| 在线亚洲精品国产二区图片欧美| 久久精品亚洲av国产电影网| 亚洲一区中文字幕在线| 人体艺术视频欧美日本| 国产1区2区3区精品| 一区二区日韩欧美中文字幕| 爱豆传媒免费全集在线观看| 国产一级毛片在线| 久久性视频一级片| 一二三四在线观看免费中文在| 色网站视频免费| 久久久久人妻精品一区果冻| 精品一品国产午夜福利视频| 视频区图区小说| 国产精品免费视频内射| 欧美激情高清一区二区三区 | 亚洲成色77777| 亚洲国产中文字幕在线视频| 午夜福利视频在线观看免费| 日本wwww免费看| 丰满饥渴人妻一区二区三| 国产熟女欧美一区二区| 亚洲国产精品一区三区| 成年av动漫网址| 亚洲精品视频女| 久久精品亚洲熟妇少妇任你| 中文字幕另类日韩欧美亚洲嫩草| 王馨瑶露胸无遮挡在线观看| 美女高潮到喷水免费观看| 成人18禁高潮啪啪吃奶动态图| 日本91视频免费播放| 男人操女人黄网站| 一级a爱视频在线免费观看| 99精国产麻豆久久婷婷| 国产成人a∨麻豆精品| 深夜精品福利| 伦理电影免费视频| 考比视频在线观看| 国产激情久久老熟女| 精品国产超薄肉色丝袜足j| 亚洲欧美一区二区三区黑人| 精品卡一卡二卡四卡免费| 亚洲男人天堂网一区| 男女边吃奶边做爰视频| 精品一区二区三卡| 成年女人毛片免费观看观看9 | 搡老岳熟女国产| 99热国产这里只有精品6| 久久久久国产精品人妻一区二区| 亚洲一码二码三码区别大吗| 在线 av 中文字幕| 国产免费福利视频在线观看| 老司机深夜福利视频在线观看 | av在线app专区| 69精品国产乱码久久久| 国产精品蜜桃在线观看| 久久性视频一级片| 国产无遮挡羞羞视频在线观看| 交换朋友夫妻互换小说| 欧美日韩精品网址| 亚洲国产欧美一区二区综合| 久久久久精品性色| 国产成人精品无人区| 夜夜骑夜夜射夜夜干| 成人免费观看视频高清| 日本av手机在线免费观看| 国产精品久久久人人做人人爽| 可以免费在线观看a视频的电影网站 | 日韩人妻精品一区2区三区| 热99国产精品久久久久久7| 国产精品一区二区在线观看99| 天美传媒精品一区二区| 男女床上黄色一级片免费看| 久久综合国产亚洲精品| 中文字幕人妻丝袜一区二区 | 看十八女毛片水多多多| 黄色毛片三级朝国网站| 成人免费观看视频高清| 色播在线永久视频| 欧美最新免费一区二区三区| 欧美乱码精品一区二区三区| 欧美激情高清一区二区三区 | 亚洲精品自拍成人| 国产精品久久久久久精品电影小说| 成年美女黄网站色视频大全免费| 女性被躁到高潮视频| 精品国产乱码久久久久久男人| 王馨瑶露胸无遮挡在线观看| 精品国产一区二区三区久久久樱花| 久久 成人 亚洲| 在线精品无人区一区二区三| 免费观看性生交大片5| 中文字幕人妻丝袜一区二区 | 人体艺术视频欧美日本| 高清欧美精品videossex| 香蕉丝袜av| 天天躁日日躁夜夜躁夜夜| 性高湖久久久久久久久免费观看| 亚洲婷婷狠狠爱综合网| 国产熟女欧美一区二区| 在现免费观看毛片| 亚洲一码二码三码区别大吗| 国产精品99久久99久久久不卡 | 纯流量卡能插随身wifi吗| 哪个播放器可以免费观看大片| 999精品在线视频| 91国产中文字幕| 欧美成人午夜精品| 大香蕉久久成人网| 如何舔出高潮| 国产 一区精品| 蜜桃国产av成人99| 国产精品无大码| 久久久精品区二区三区| 观看av在线不卡| 亚洲国产成人一精品久久久| 国产麻豆69| 国产乱人偷精品视频| 亚洲精品国产区一区二| 只有这里有精品99| 韩国高清视频一区二区三区| 婷婷色综合大香蕉| 亚洲av在线观看美女高潮| 日韩精品有码人妻一区| 精品人妻一区二区三区麻豆| 99久久精品国产亚洲精品| 91成人精品电影| 桃花免费在线播放| 国产成人系列免费观看| 日本猛色少妇xxxxx猛交久久| 在线观看免费视频网站a站| 亚洲欧洲日产国产| 午夜福利一区二区在线看| 中文精品一卡2卡3卡4更新| 国产免费又黄又爽又色| 丰满乱子伦码专区| 亚洲精品一二三| av天堂久久9| 亚洲欧美成人精品一区二区| 免费观看性生交大片5| 国产在线免费精品| 欧美乱码精品一区二区三区| 久久99一区二区三区| a 毛片基地| 视频在线观看一区二区三区| 蜜桃国产av成人99| 国产精品秋霞免费鲁丝片| 亚洲美女搞黄在线观看| 亚洲免费av在线视频| 国产亚洲精品第一综合不卡| 2018国产大陆天天弄谢| av视频免费观看在线观看| 伦理电影免费视频| 不卡视频在线观看欧美| 嫩草影院入口| 狠狠婷婷综合久久久久久88av| 在现免费观看毛片| 亚洲精品日韩在线中文字幕| 久久久久久免费高清国产稀缺| 国产一区二区三区av在线| 中文字幕亚洲精品专区| 99热网站在线观看| 美女中出高潮动态图| 欧美少妇被猛烈插入视频| 18禁国产床啪视频网站| 亚洲情色 制服丝袜| 久久鲁丝午夜福利片| 午夜日韩欧美国产| 满18在线观看网站| 青春草国产在线视频| 亚洲色图 男人天堂 中文字幕| 999久久久国产精品视频| 色婷婷av一区二区三区视频| 国产淫语在线视频| 黑人欧美特级aaaaaa片| 欧美激情极品国产一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品一区二区三区在线| 欧美精品亚洲一区二区| 国产一区二区三区av在线| 久久久久久人妻| 最新在线观看一区二区三区 | 国产在线视频一区二区| 一区二区三区乱码不卡18| 国产亚洲最大av| 国产av一区二区精品久久| 9色porny在线观看| 国产成人精品在线电影| 亚洲精品日本国产第一区| 视频在线观看一区二区三区| 女的被弄到高潮叫床怎么办| 日韩大片免费观看网站| av视频免费观看在线观看| 亚洲av电影在线观看一区二区三区| 亚洲精品美女久久av网站| 美女福利国产在线| 国产福利在线免费观看视频| 一边摸一边做爽爽视频免费| 国产精品一国产av| 国产乱人偷精品视频| 午夜激情久久久久久久| 亚洲国产精品999| 日本av手机在线免费观看| 久久女婷五月综合色啪小说| 亚洲国产精品成人久久小说| 亚洲精品视频女| 悠悠久久av| 日韩制服骚丝袜av| h视频一区二区三区| 中国国产av一级| 女人久久www免费人成看片| 亚洲自偷自拍图片 自拍| 美女大奶头黄色视频| 啦啦啦啦在线视频资源| 国产精品秋霞免费鲁丝片| 国产免费又黄又爽又色| 久久人人爽av亚洲精品天堂| 90打野战视频偷拍视频| 成年人午夜在线观看视频| 最新的欧美精品一区二区| 在线亚洲精品国产二区图片欧美| 大码成人一级视频| 人人妻人人爽人人添夜夜欢视频| 精品国产国语对白av| 精品福利永久在线观看| 亚洲国产欧美日韩在线播放| 狂野欧美激情性xxxx| 1024香蕉在线观看| av卡一久久| 国产女主播在线喷水免费视频网站| 性色av一级| 精品久久久久久电影网| 亚洲精品久久成人aⅴ小说| 午夜免费鲁丝| 叶爱在线成人免费视频播放| 亚洲欧美日韩另类电影网站| 国产极品粉嫩免费观看在线| 国产免费又黄又爽又色| 午夜免费鲁丝| 精品免费久久久久久久清纯 | 捣出白浆h1v1| 又黄又粗又硬又大视频| 在线观看免费日韩欧美大片| 涩涩av久久男人的天堂| 中文字幕最新亚洲高清| 午夜福利在线免费观看网站| 曰老女人黄片| 一区在线观看完整版| 少妇精品久久久久久久| 亚洲情色 制服丝袜| 色网站视频免费| 久久久久国产精品人妻一区二区| 亚洲欧洲国产日韩| 日韩欧美一区视频在线观看| av在线播放精品| 男女免费视频国产| 热re99久久精品国产66热6| 桃花免费在线播放| 熟女av电影| 欧美日韩国产mv在线观看视频| 色精品久久人妻99蜜桃| 亚洲专区中文字幕在线 | 亚洲成av片中文字幕在线观看| 一区二区三区乱码不卡18| 18在线观看网站| 高清在线视频一区二区三区| 国产精品一区二区在线不卡| 99精国产麻豆久久婷婷| 日韩中文字幕欧美一区二区 | 欧美人与性动交α欧美软件| 成人手机av| 丝袜喷水一区| 久久人人爽av亚洲精品天堂| 亚洲第一区二区三区不卡| 亚洲精品在线美女| 久久久久视频综合| 亚洲在久久综合| 国产精品熟女久久久久浪| 亚洲熟女精品中文字幕| 日韩视频在线欧美| 日日爽夜夜爽网站| 亚洲国产欧美一区二区综合| 亚洲图色成人| 久久久久人妻精品一区果冻| 国产有黄有色有爽视频| 免费看av在线观看网站| 在线观看三级黄色| 成年人免费黄色播放视频| 成年人午夜在线观看视频| 99香蕉大伊视频| 欧美在线黄色| 久久青草综合色| 纵有疾风起免费观看全集完整版| 久久久久久久国产电影| 丰满乱子伦码专区| 久久国产精品大桥未久av| 国产日韩欧美视频二区| 国产成人免费观看mmmm| 男女床上黄色一级片免费看| 在线 av 中文字幕| 国产极品粉嫩免费观看在线| 赤兔流量卡办理| 久久精品国产综合久久久| 1024视频免费在线观看| 精品亚洲成a人片在线观看| 久久午夜综合久久蜜桃| 国产伦人伦偷精品视频| 飞空精品影院首页| 国产成人啪精品午夜网站| 成人亚洲精品一区在线观看| 一级黄片播放器| 最近最新中文字幕免费大全7| 国产精品二区激情视频| 美女国产高潮福利片在线看| 亚洲欧美成人综合另类久久久| 久热这里只有精品99| 99久久99久久久精品蜜桃| 一区二区三区精品91| 午夜影院在线不卡| 一二三四中文在线观看免费高清| 日韩av在线免费看完整版不卡| 涩涩av久久男人的天堂| 精品少妇久久久久久888优播| 精品福利永久在线观看| 亚洲国产精品国产精品| 丰满迷人的少妇在线观看| 亚洲欧美一区二区三区黑人| 亚洲专区中文字幕在线 | 午夜福利免费观看在线| 亚洲三区欧美一区| 亚洲国产欧美在线一区| 久久久久久久久久久免费av| 建设人人有责人人尽责人人享有的| 国产精品人妻久久久影院| 国产成人精品福利久久| 久久精品亚洲熟妇少妇任你| 国产精品国产三级专区第一集| 老汉色∧v一级毛片| 国产97色在线日韩免费| 欧美激情 高清一区二区三区| 免费日韩欧美在线观看| 亚洲美女视频黄频| 午夜福利,免费看| 五月天丁香电影| 久久这里只有精品19| 丝袜人妻中文字幕| 欧美av亚洲av综合av国产av | 国产女主播在线喷水免费视频网站| 老汉色∧v一级毛片| 女人精品久久久久毛片| 考比视频在线观看| 王馨瑶露胸无遮挡在线观看| 如日韩欧美国产精品一区二区三区| 国产亚洲av片在线观看秒播厂| 激情五月婷婷亚洲| 国产精品.久久久| 下体分泌物呈黄色| 亚洲免费av在线视频| 国产在线免费精品| 亚洲熟女毛片儿| 国产男女内射视频|