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

    氫抑制條件下的有機質(zhì)熱演化動力學(xué)模型
    ——以三塘湖盆地蘆草溝組為例

    2017-11-09 09:50:29孟元林崔存蕭李新寧陶士振吳晨亮胡安文
    石油與天然氣地質(zhì) 2017年5期
    關(guān)鍵詞:富氫鏡質(zhì)塘湖

    張 磊,孟元林,崔存蕭,李新寧,陶士振,吳晨亮,胡安文,許 丞

    (1.東北石油大學(xué) 地球科學(xué)學(xué)院,黑龍江 大慶 163318; 2.中國石油 吐哈油田分公司 勘探開發(fā)研究院,新疆 哈密 839009;3.中國石油 勘探開發(fā)研究院 ,北京 100083; 4.萊斯大學(xué) 地球、環(huán)境與空間系,休斯頓,77005)

    氫抑制條件下的有機質(zhì)熱演化動力學(xué)模型
    ——以三塘湖盆地蘆草溝組為例

    張 磊1,孟元林1,崔存蕭1,李新寧2,陶士振3,吳晨亮4,胡安文1,許 丞1

    (1.東北石油大學(xué) 地球科學(xué)學(xué)院,黑龍江 大慶 163318; 2.中國石油 吐哈油田分公司 勘探開發(fā)研究院,新疆 哈密 839009;3.中國石油 勘探開發(fā)研究院 ,北京 100083; 4.萊斯大學(xué) 地球、環(huán)境與空間系,休斯頓,77005)

    三塘湖盆地馬朗凹陷蘆草溝組烴源巖含有較多的富氫組分,有機質(zhì)的熱演化受到氫抑制。為了模擬氫抑制背景下的有機質(zhì)熱演化史,在經(jīng)典的Sweeny和Burham Easy%Ro時-溫雙控化學(xué)動力學(xué)模型基礎(chǔ)上,本文考慮氫抑制對有機質(zhì)熱演化的影響,根據(jù)氫指數(shù)調(diào)節(jié)反應(yīng)活化能,建立了描述有機質(zhì)熱演化的時-溫-氫抑制多控化學(xué)動力學(xué)模型。研究結(jié)果表明,蘆草溝組氫抑制對Ro的貢獻值為0.2%左右;由于氫抑制的影響,蘆草溝組烴源巖進入生烴門限和開始大量生烴的時間均被延遲。文中新建立的化學(xué)動力學(xué)模型可以應(yīng)用于各種類型干酪根的有機質(zhì)熱演化史模擬,其精度更高、應(yīng)用范圍更廣。

    氫抑制;化學(xué)動力學(xué);頁巖油;中二疊統(tǒng);三塘湖盆地

    鏡質(zhì)組反射率Ro是烴源巖評價、成巖階段劃分、蓋層封閉性研究、成巖作用數(shù)值模擬、成藏動力學(xué)模擬和盆地模擬最常用的指標[1-7]。傳統(tǒng)的觀點認為,鏡質(zhì)組反射率主要受溫度與時間的影響與控制[8],也有人認為壓力或超壓對Ro有明顯的影響[9-13],Carr(2000),Zou (2001)和肖麗華(2003)考慮超壓對Ro的影響,分別建立了描述有機質(zhì)熱演化的時-溫-壓多控化學(xué)動力學(xué)模型[14-16]。但越來越多的研究表明,有機質(zhì)的熱演化受到氫的抑制作用[17-24],即在富含氫元素的干酪根中有機質(zhì)的熱演化受到氫抑制(Hydrogen Suppression),Ro隨埋深和地溫的增加演化緩慢,Ro的實測值偏低。而形成頁巖油或致密油聚集的烴源巖常常是氫組分含量很高的優(yōu)質(zhì)烴源巖[25],有機質(zhì)的熱演化受到氫抑制,比較典型的是Williston盆地Barkken 頁巖的氫抑制現(xiàn)象[21]。因此,在研究富含藻質(zhì)體和殼質(zhì)組等顯微組分的泥巖有機質(zhì)熱演化時,必須考慮氫抑制對Ro的影響。更重要的是,以前被廣泛應(yīng)用的Sweeny 和Burham的Easy%Ro時-溫雙控模型不再適用,需要加以改進。本文試圖以三塘湖盆地為例,研究有機質(zhì)熱演化被氫抑制的特征,建立氫抑制背景下有機質(zhì)熱演化的化學(xué)動力學(xué)模型。

    1 地質(zhì)背景

    三塘湖盆地位于新疆維吾爾自治區(qū)境內(nèi),總面積為2.3×104km2,可進一步劃分為北部隆起帶、中央坳陷帶和南緣沖斷帶3個二級構(gòu)造單元(圖1a)[26]。條湖凹陷和馬朗凹陷相鄰,位于中央坳陷帶的中部,是三塘湖盆地最有利的勘探區(qū)。三塘湖盆地從下到上發(fā)育上古生界石炭系(C)、二疊系(P),中生界三疊系(T)、侏羅系(J)、白堊系(K)和新生界(Q+R)。其中,上古生界僅發(fā)育中二疊統(tǒng),由蘆草溝組(P2l)和條湖組(P2t)構(gòu)成。蘆草溝組從下到上分為蘆一段(P2l1)、蘆二段(P2l2)和蘆三段(P2l3)(圖1b)。蘆草溝組屬于咸化湖碳酸鹽和泥質(zhì)巖混合沉積,在沉積的同時還伴有火山噴發(fā)和火山灰堆積,主要巖性為泥質(zhì)沉凝灰?guī)r、碳酸鹽質(zhì)沉凝灰?guī)r、沉凝灰?guī)r和泥質(zhì)凝灰質(zhì)碳酸鹽巖[26-27]。

    蘆草溝組富含有機質(zhì),發(fā)育優(yōu)質(zhì)烴源巖,源儲一體,形成了致密油或頁巖油聚集。鉆遇蘆草溝組的48口井中,有28口井測試見油,7口井產(chǎn)出工業(yè)油流。

    2 烴源巖有機質(zhì)熱演化的氫抑制與機理解釋

    2.1 富氫有機質(zhì)的熱演化特征與氫抑制

    蘆草溝組是三塘湖盆地一套極為優(yōu)質(zhì)的烴源巖,該盆地80%的儲量來自這套烴源巖[26-27]。由于凝灰質(zhì)成分對有機質(zhì)的富集作用[28],烴源巖有機質(zhì)極為豐富,有機質(zhì)呈層狀富集(圖2)。馬朗凹陷蘆草溝組烴源巖的TOC主要分布在2.0%~8.0%,條湖凹陷蘆草溝組TOC主要分布1.0%~4.0%。有機質(zhì)類型為主要Ⅰ型和Ⅱ1型,富含氫元素,H/C原子比為1.0~1.5,O/C原子比為0.05~0.15。顯微組分中富氫組分含量很高,殼質(zhì)組和腐泥組的含量主要分布在40%~98%(圖3),富氫組分含量高的烴源巖主要發(fā)育在馬朗凹陷。

    正是由于蘆草溝組(P2l)烴源巖干酪根中含有大量富氫組分,蘆草溝組烴源巖的有機質(zhì)熱演化很慢,受到氫抑制。在深度相同的情況下,蘆草溝組(P2l)烴源巖的Ro低于三疊系(T)、二疊系條湖組(P2t)和石炭系(C)烴源巖的Ro(圖4)。

    值得說明的是,超壓也可以抑制有機質(zhì)的熱演化[9-16]。但是,由于三塘湖盆地蘆草溝組屬于一套由碳酸鹽、凝灰質(zhì)和陸源碎屑混雜堆積形成的細粒沉積,非常致密,欠壓實的程度很弱;而且有機質(zhì)熱演化程度較低,目前處于液態(tài)窗內(nèi),由生烴產(chǎn)生的壓力有限,所以實測地層壓力不高。吐哈油田的實測地層壓力表明[29-30],蘆草溝組的壓力系數(shù)分布在0.5~1.29,主要集中在1.0~1.2的范圍內(nèi)(圖5),基本屬于正常壓力,僅僅少數(shù)井存在弱超壓(壓力系數(shù)≥1.2),超壓對有機質(zhì)熱演化的影響不大,可以忽略。

    圖1 研究區(qū)位置(a)及地層柱狀圖(b)[26]Fig.1 Location and stratigraphic column of the study area[26]

    圖2 馬朗凹陷蘆草溝組烴源巖有機質(zhì)富集分布Fig.2 Distribution of organic matters abundance in the source rocks of the Lucaogou Formation in the Malang Saga.蘆1井,蘆草溝組,埋深3 172.68 m,有機質(zhì)(OM)層狀富集(綠色箭頭),Si為石英,F(xiàn)p為斜長石,Ca為方解石,F(xiàn)k為鉀長石,紅色箭頭為白云石; b.蘆1井,蘆草溝組,埋深3 123.68 m,有機質(zhì)層狀富集(綠色箭頭),F(xiàn)k為鉀長石,Ab為斜長石,Q為石英,Py為黃鐵礦(紅色箭頭)

    圖3 三塘湖盆地蘆草溝組干酪根顯微組分組成Fig.3 Organic macerals of kerogen in the Lucaogou Formation in the Santanghu Basin

    2.2 氫抑制的化學(xué)機理

    許多學(xué)者研究過鏡質(zhì)組反射率被氫抑制的原因。Taylor(1966,1987)最早在無定形有機質(zhì)占優(yōu)勢的烴源巖樣品中,發(fā)現(xiàn)了Ro氫抑制的現(xiàn)象[17-18]。而Hutton和Cook(1980)認為,Ro的氫抑制是富氫鏡質(zhì)體本身造成的[19]。Kalkreuth(1982)的研究結(jié)果則表明[20],氫抑制由干酪根中高含量的殼質(zhì)體引起。Leigh和Charles(1985)認為[21],Wilingston盆地Bakken頁巖Ro的抑制現(xiàn)象與高含量的藻質(zhì)體有關(guān),藻質(zhì)體的含量與鏡質(zhì)組反射率呈負相關(guān)??傊彩桥c富氫顯微組分共生的鏡質(zhì)組,它們的Ro測定值往往比預(yù)料的要低。而且Ro抑制現(xiàn)象多發(fā)生于H/C原子比較高的Ⅰ,Ⅱ型干酪根烴源巖中[22]。三塘湖盆地蘆草溝組烴源巖殼質(zhì)組與Ro亦呈負相關(guān)關(guān)系(圖6),顯示了鏡質(zhì)體反射率被氫抑制的特征。

    圖4 三塘湖盆地鏡質(zhì)組反射率Ro與深度關(guān)系Fig.4 Relation beween Ro and depth in the Santanghu Basin

    圖5 三塘湖盆地實測壓力系數(shù)與深度關(guān)系[29]Fig.5 Measured pressure coefficient vs.depth in the Santanghu Basin[29]

    圖6 三塘湖盆地蘆草溝組干酪根殼質(zhì)組含量與Ro關(guān)系Fig.6 Relation between Ro and exinite of kerogen in the Lucaogou Formation in the Santanghu Basin

    由上可見,抑制有機質(zhì)熱演化的主要原因是干酪根中含有較多的富氫組分。有機質(zhì)熱演化氫抑制的化學(xué)動力學(xué)原理是,干酪根中的富氫顯微組分(藻質(zhì)體、殼質(zhì)體等)含有許多的脂族鏈,在有機質(zhì)熱演化的過程中,這些脂族鏈的C—C鍵斷裂,生成烴類,而C—C鍵斷裂需要的反應(yīng)活化能比貧氫顯微組分生烴時雜原子鍵斷裂所必需的活化能更高。也就是說,富氫鏡質(zhì)組熱演化生烴所需的活化能比貧氫鏡質(zhì)組更高。因此熱演化生烴的速度更慢,在經(jīng)歷同樣時-溫的條件下,富氫鏡質(zhì)組的Ro更低。

    3 氫抑制條件下的有機質(zhì)熱演化化學(xué)動力學(xué)模型

    在Sweeny和Burham(1990)有機質(zhì)熱演化Easy%Ro時-溫雙控化學(xué)動力學(xué)模型的基礎(chǔ)上,本文考慮氫抑制的影響,通過增加熱演化化學(xué)反應(yīng)的活化能,建立了時-溫-氫抑制多控的化學(xué)動力學(xué)模型,描述富氫有機質(zhì)熱演化的過程:

    (1)

    (2)

    Ro%=exp(-1.6+3.7F)

    (3)

    (4)

    式中:wi為干酪根鏡質(zhì)組中組分i的反應(yīng)物濃度,%,i=1,2,…,20;wi0為干酪根鏡質(zhì)組中組分i的初始濃度,%;A為頻率因子,1.13×103/s;t為時間,s;T(t)為t時刻的絕對溫標,K;E為活化能, kcal/mol;Pi為權(quán)系數(shù);λ為氫抑制因子,λ>1,可用下式計算:

    λ=1+β[ (IH-IH0) /IH]

    (5)

    式中:IH為氫指數(shù),mg/g;IH0為不發(fā)生氫抑制的最大氫指數(shù),IH0=150mg/g[23];β為比例因子。

    由公式(5)可見,當IH=IH0時,λ=1,與Sweeny 和Burham(1990)提出的Easy%Ro經(jīng)典模型完全相同;當IHIH0時,可根據(jù)Ro的實測值和模擬值,用“試錯法”求得β和λ,由計算機按照從小到大的順序反復(fù)選取一系列的β,計算出λ,直到鏡質(zhì)組反射率的模擬值達到給定的精度為止:

    (6)

    4 參數(shù)選取與模擬計算

    4.1 氫抑制因子λ的確定

    為了定量研究富氫組分對Ro抑制的程度、確定氫抑制因子λ、模擬馬朗凹陷蘆草溝組有機質(zhì)熱演化的化學(xué)動力學(xué)過程,本文在條湖凹陷和馬朗凹陷分別選取了一口蘆草溝組有機質(zhì)熱演化沒有受到氫抑制的條5井和受到氫抑制的塘參3井,進行了試算。

    首先,以條湖凹陷條5井的實測Ro資料為基礎(chǔ),應(yīng)用Sweeny和Burham(1990)的Easy%Ro經(jīng)典模型,反演出條5井的古大地?zé)崃?圖7),直到鏡質(zhì)組反射率的計算值與實測值相符為止(圖8a);然后,根據(jù)由于這兩口井基底埋深差異引起的大地?zé)崃髡凵淠P蚚31-34],由基底埋藏較淺的條5井的古大地?zé)崃?,計算出基底埋藏較深的塘參3井的古大地?zé)崃?圖7),并以此恢復(fù)塘參3井不同地質(zhì)時期的古地溫,用Sweeny和Burham(1990)的Easy%Ro模型,計算出塘參3井的Ro值(圖8b)。由圖可見,塘參3井Ro實測值小于理論計算值,這是氫抑制有機質(zhì)熱演化的結(jié)果。其差值ΔRo就是氫抑制的貢獻值,ΔRo=0.20%左右。最后,應(yīng)用本文建立的考慮氫抑制的化學(xué)動力學(xué)模型[公式(1)~(6)],由計算機給定一系列的βi值(i=1,2,…,N),得出相應(yīng)的氫抑制因子λi值,反復(fù)試算,直到計算的精度達到給定的精度[公式(6)]。計算結(jié)果表明,當β=0.073 5,氫抑制因子λ=1.044 4時,鏡質(zhì)組反射率的計算值和實測數(shù)據(jù)符合較好(圖8c)。

    4.2 有機質(zhì)熱演化史對比

    應(yīng)用得到的氫抑制因子λ值,模擬了受到氫抑制的塘參3井的有機質(zhì)熱演化史(圖9a)。由圖可見,蘆草溝組(P2l)烴源巖在侏羅紀末期才進入生油門限(Ro≥0.5%);在白堊紀末期進入大量生烴的成熟階段(Ro≥0.7%)。如果塘參3井不存在氫抑制,應(yīng)用Sweeny和Burham(1990)的Easy%Ro經(jīng)典模型恢復(fù)該井的有機質(zhì)熱演化史(圖9b),蘆草溝組(P2l)烴源巖在中二疊統(tǒng)條湖組(P2t)沉積末期就已進入生油門限;在白堊紀中期進入大量生烴的成熟階段。由圖9可知,正是由于鏡質(zhì)組反射率被氫抑制的結(jié)果,使蘆草溝組(P2l)烴源巖進入生烴門限的時間和開始大量生烴的時間均被大大延遲。

    5 結(jié) 論

    1) 三塘湖盆地馬朗凹陷蘆草溝組烴源巖有機質(zhì)富含氫組分,有機質(zhì)熱演化被氫抑制,氫抑制的效應(yīng)為ΔRo=0.2%。

    圖7 三塘湖盆地條5井與塘參3井不同地質(zhì)時期大地?zé)崃鲗Ρ菷ig.7 Heat flows comparison of Well Tiao-5 and Well Tangcan-3 in different geologic periods of the Santanghu Basin

    圖8 三塘湖盆地Ro實測值與模擬值Fig.8 Measured Ro and modeled Ro in the Santanghu Basina.條5井Ro實測值與模擬值;b.塘參3井Ro實測值與Easy%Ro模型計算的理論值;c.塘參3井Ro實測值與氫抑制模型計算的理論值

    圖9 三塘湖盆地地?zé)崾返幕謴?fù)Fig.9 Reconstructed thermal history in the Santanghu Basina.用氫抑制模型恢復(fù)的地?zé)崾?b.用Sweeny和Burham Easy% Ro模型恢復(fù)的地?zé)崾?/p>

    2) 通過調(diào)節(jié)氫抑制因子λ、增加有機質(zhì)熱演化反應(yīng)的活化能,所建立的有機質(zhì)熱演化時-溫-氫抑制多控模型可以更好地模擬富氫有機質(zhì)的熱演化史,而且精度更高、應(yīng)用范圍更廣,可以模擬各種類型有機質(zhì)的熱演化過程。

    3) 由于氫抑制的影響,三塘湖盆地蘆草溝組富含殼質(zhì)組和腐泥組的烴源巖進入生烴門限和開始大量生烴的時間均被推遲。

    [1] 黃飛,辛茂安.SY/T5735-1995中華人民共和國石油天然氣行業(yè)標準并陸相烴源巖地球化學(xué)評價方法[S].北京:石油工業(yè)出版社,1996:1-19.

    Huang Fei and Xin Maoan.Oil and gas industry standard of the People’s Republic of China:The geochemical evaluation method of continental-facies source rocks,SY/T5735-1995 [S].Beijing:Petroleum Industry Press,1996:1-19.

    [2] 應(yīng)鳳祥,何東博,龍玉梅,等.SY/T5477-2003 中華人民共和國石油天然氣行業(yè)標準并碎屑巖成巖階段劃分[S].北京:石油工業(yè)出版社,2003:1-5.

    Ying Fengxiang,He Dongfu,Long Yumei,et al.Oil and gas industry standard of the People’s Republic of China:The division of diagentic stages in clastic rocks,SY/T5477-2003 [S].Beijing:Petroleum Industry Press,2003:1-5.

    [3] 龐雄奇,付廣,萬龍貴,等.蓋層封油氣性綜合定量評價—盆地模擬在蓋層評價中的應(yīng)用[M].北京:地質(zhì)出版社,1993:14-17.

    Pang Xiongqi,F(xiàn)u Guang,Wan Longgui,et al.Caprock sealing comprehensive quantitative assessment of oil and gas basin simulation applied in the evaluation of cap rock[M].Beijing:Geological Publishing House,1993:14-17.

    [4] Hantschel T,Kauerauf I A.Fundamentals of basin and petroleum systems modeling[M].Berlin Heidelberg:Springer-Verlag,2009:169-176.

    [5] 孟元林,劉德來,賀如,等.歧北凹陷超壓背景下的成巖場分析與儲層孔隙度預(yù)測[J].沉積學(xué)報,2005,23(3):389-396.

    Meng Yuanlin,Liu Delai,He Ru,et al.Diagentic field analysis and porosity prediction of the Shaer member(Es2) in overpressure setting in the Qibei depression[J].Acta Sedimentologica Sinica,2005,23(3):389-396.

    [6] 孟元林,王粵川,羅憲嬰,等.渤海灣盆地孔西潛山構(gòu)造帶成藏史數(shù)值模擬[J],地質(zhì)力學(xué)學(xué)報,2005,11(1):11-16.

    Meng Yuanlin,Wang Yuechuan,Luo Xianying,et al.Reservoir forming history modeling of Kongxi buried hill,Bohaiwan Basin[J].Journal of Geomechanics,2005,11(1):11-16.

    [7] 黃文彪,鄧守偉,盧雙舫,等.松遼盆地南部扶余油層致密儲層成巖序列及成藏期次[J].石油與天然氣地質(zhì),2017,38(03):508-516.

    Huang Wenbiao,Deng Shouwei,Lu Shuangfang,et al.Diagenetic sequence and hydrocarbon accumulation phases of the Fuyu layer tight reservoir in the southern Songliao Basin[J].Oil & Gas Geology,2017,38(03):508-516.

    [8] Sweeny J J,Burham A K.Evaluation of a simple model of vitrinite reflectance based on chemical kinetics[J].AAPGBulletin,1990,74:1559-1570.

    [9] McTavish R A.Pressure retardation of vitrinite diagenesis,offshore north—West Europe[J].Nature,1978,271(16):648-650.

    [10] Hao F,Li S T,Dong W H,et al.Abnormal organic matter maturation in the Yinggehai Basin,offshore South China Sea:implications for hydrocarbon expulsion and fluid migration from overpressured systems[J].Journal of Petroleum Geology,1998,21:427-444.

    [11] Huang W L.Experimental study of vitrinite maturation:effects of temperature,time,pressure,water,and hydrogen index[J].Organic Geochemistry,1996,24(2):233-241.

    [12] Dalla T M,Mahlman F R,Ernt W G.Experimental study on the pressure dependence of vitrinite maturation[J].Organic Geochemistry,1997,61:2921-2928.

    [13] 王永詩,邱貽博.濟陽坳陷超壓結(jié)構(gòu)差異性及其控制因素[J].石油與天然氣地質(zhì),2017,38(03):430-437.

    Wang Yongshi,Qiu Yibo.Overpressure structure dissimilarity and its controlling factors in the Jiyang Depression[J].Oil & Gas Geology,2017,38(03):430-437.

    [14] Carr A D.A vitrinite reflectance kinetic model incorporating overpressure retardation[J].Marine and Petroleum Geology,1999,16:355-377.

    [15] Zou Y,Peng P.Overpressure retardation of organic-matter maturation: a kinetic model and its application[J].Marine and Petroleum Geology,2001,18: 707-713.

    [16] 肖麗華,盂元林,張連雪,等.超壓地層中鏡質(zhì)組反射率的計算[J].石油勘探與開發(fā),2005,32(1):14-17.

    Xiao Lihua,Meng Yuanlin,Zhang Lianxue,et al.Vitrinite reflectance modeling in the pressured formations[J].Petroleum Exploration and Development,2005,32(1):14-17.

    [17] Taylor G.H.The electron microscopy of vitrinites.In Coal Sciences (Edited by Gould R.F.) [D].Am.Chemical Society Publications,Washington,D.C.1966,274-283.

    [18] Taylor G H,Liu S Y.Biodegradation in coals and other organic-rich rocks[J].Fuel,1987,66(9):1269-1273.

    [19] Hutton A C,Cook A C.Influence of alginate on the reflectance of vitrinite from Joadja.NSW,and some other coals and oil shales containing alginite[J].Fuel,1980,59(10):711-714.

    [20] Kalkreuth W D.Rank and petrolographic composition of selected Jurassic-Lower Cretaceous coals of British Columbia,Canada[J].Bulletin of Canadian Petroleum Geology,1982,30(2):112-139.

    [21] Leigh C P,Charles E B.Suppression of vitrinite reflectance inamorphous rich kerogen-amajor unrecognized problem[J].Journal of Petroleum Geology,1985,8(1):59-84.

    [22] Hao F,Chen J.The cause and mechanism of vitrinite refletance anomolies[J].Journal of Petroleum Geology,1992,15(4):419-434.

    [23] Lo H B.Correction criteria for the suppression of vitrinite reflectance in hydrogen-rich kerogens: preliminary guidelines[J].Organic Geochemistry,1993,20(6):653-657.

    [24] Dow W G.Kerogen studies and geological interpretations[J].Journal of Geochemical Exploration,1977,7:79-99.

    [25] 鄒才能.非常規(guī)油氣地質(zhì)(第二版)[M].北京:石油工業(yè)出版社,2013:93-169.

    Zou Caineng.Unconventional petroleum geology(second edition)[M].Beijing:Petroleum Industry Press,2011.

    [26] 李新寧,馬強,梁輝,等.三塘湖盆地蘆草溝組二段混積巖致密油地質(zhì)特征及勘探潛力[J].石油勘探與開發(fā),2015,42(6):1-9.

    Li Xinning,Ma Qiang,Liang Hui,et al.Geological characteristics and exploration potential of diamictite tight oil in the second Member of the Permian Lucaogou Formation,Santanghu Basin[J].Petroleum Exploration and Development,2015,42(6):1-9.

    [27] 國建英,鐘寧寧,梁浩,等.三塘湖盆地中二疊統(tǒng)原油的來源及其分布特征[J].地球化學(xué),2012,41(3):266-277.

    Guo Jianying,Zhong Ningning,Liang Hao,et al.Study on the source and distribution of Middle Permian oils in the Santanghu Basin[J].Geochimica,2012,41(3):266-277.

    [28] 孟元林,張磊,李新寧,等.雜積巖地化特征與生烴機理—以三塘湖盆地蘆草溝組為例,第十五屆全國有機地球化學(xué)學(xué)術(shù)會議,青島,2015.

    Meng Yuanlin,Zhang Lei,Li Xinning,et al.Mixed siliciclastic-volcaniclastic-carbonate sediments—taking Lu Caogou Formation of Santanghu Basin as an example,The 15thNational Meeting on Organic Geochemistry in China,Qingdao,2015.

    [29] 李新寧.三塘湖盆地蘆草溝組頁巖油藏地質(zhì)特征與勘探開發(fā)潛力[R].2011,中國石油吐哈油田分公司勘探開發(fā)研究院.

    Li Xinning .Geologic characteristics and exploration and development potential of the Lucaogou shale oil play in Santanghu Basin[R].2011,Exploration and Development Institute of Tuha Oil-field Branch,Petrochina.

    [30] 黃志龍,馬劍,吳紅燭,等.馬朗凹陷蘆草溝組頁巖油流體壓力與初次運移特征[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2012,36(5):7-11

    Huang Zhilong,Ma Jian,Wu Honghu,et al .Fluid pressure and primary migration characteristics of shale oil of Lucaogou formation in Malang sag[J].Journal of China University of Petroleum,2012,36(5):7-11.

    [31] 張菊明,熊亮萍.有限單元法在地?zé)嵫芯恐械膽?yīng)用[M].北京:科學(xué)出版社,1986:15-22.

    Zhang Juming,Xiong Liangping.The applications of finite element in geothermal researches[M].Beijing: Science Press,1986:15-22.

    [32] 肖麗華,張靖,孟元林,等.地?zé)釁?shù)及邊界條件的探討[J].大慶石油學(xué)院學(xué)報,1996,20(2):28-31.

    Xiao Lihua,Zhang Jing,Meng Yuanlin,et al.Study of geothermal parameters and boundary conditions[J].Journal of Daqing Petroleum Institute,1996,20(2): 28-31.

    [33] Meng Y L,Xiao L H,Zhang J.Basin modeling by gravity,magnetics and electrical information and its application[C].//Liu B.T.,Li S.T.eds.Basin analysis,global sedimentary geology and sedimentology.Amsterdam:VSP,1997:97-207.

    [34] 崔軍平,任戰(zhàn)利,李金翔,等.海拉爾盆地呼倫湖凹陷熱演化史恢復(fù)[J].石油與天然氣地質(zhì),2015,36(01):35-42.

    Cui Junping,Ren Zhanli,Li Jinxiang,et al.Reconstruction of geothermal history in Hulunhu Depression,Hailaer Basin[J].Oil & Gas Geology,2015,36(01):35-42.

    A new kinetic model of organic thermal evolution under the condition of hydrogen suppression:A case study from the Lucaogou Formation in the Santanghu Basin

    Zhang Lei1,Meng Yuanlin1,Cui Cunxiao1,Li Xinning2,Tao Shizhen3, Wu Chenliang4,Hu Anwen1, Xu Cheng1

    (1.CollegeofEarthScience,NortheastPetroleumUniversity,Daqing,Heilongjiang163318,China;2.ResearchInstituteofPetroleumExploration&Development,TuhaOilfieldBranch,Petrochina,Hami,Xinjiang839009,China;3.ResearchInstituteofPetroleumExploration&Development,Petrochina,Beijing100083,China;4.DepartmentofEarthScience,RiceUniversity,Houston77005,USA)

    Organic matters of the source rocks in the Lucaogou Formation of the Malang Sag,Santanghu Basin are rich in hydrogen components,which suppresses the thermal evolution of organic matters.To simulate thermal evolution of organic matters,a new kinetic model of organic thermal evolution under the condition of hydrogen suppression was established integrating time,temperature,and hydrogen retardation based on Sweeny and Burham Easy %Rotime and temperature controlling kinetic model.In the new kinetic model,activity energies are adjusted by hydrogen index.The results show that hydrogen suppression contributes to 0.2% ofRovalue in the Lucaogou Formation.Hydrogen suppression also resulted in delays of both initial generation of hydrocarbons and peak generation of hydrocarbons.This new established model is widely applicable to thermal histories for different kerogen types with higher accuracy.

    hydrogen suppression,chemical kinetics,shale oil,Middle Permian,Santanghu Basin

    2016-01-05;

    2017-08-15。

    張磊(1990—),男,博士,石油地質(zhì)。E-mail:zhlangei245@hotmail.com。

    孟元林(1960—),男,教授、博士生導(dǎo)師,石油地質(zhì)。E-mail:qhdmyl@163.com。

    國家自然科學(xué)基金項目(41572135);國家科技重大專項(2016ZX05046-001-006);東北石油大學(xué)研究生創(chuàng)新科研項目(YJSCX2017-004NEPU);黑龍江省科學(xué)基金項目(QC2016049);黑龍江省普通本科高等學(xué)校青年創(chuàng)新人才培養(yǎng)計劃項目(UNPYSCT-2016122)。

    0253-9985(2017)05-0862-07

    10.11743/ogg20170504

    TE122.1

    A

    (編輯 張玉銀)

    猜你喜歡
    富氫鏡質(zhì)塘湖
    燒結(jié)工序降低固體燃耗節(jié)能減碳的措施
    水蒸氣氣化制備富氫燃氣
    利用鏡質(zhì)組反射率鑒定蘭炭與煤粉互混樣的方法解析
    冶金能源(2022年5期)2022-10-14 00:35:16
    三塘湖油田馬中致密油水平井大修技術(shù)研究與實踐
    科學(xué)家(2022年3期)2022-04-11 00:48:50
    大落差管道改性輸送工藝研究
    化工機械(2021年5期)2021-10-28 08:47:10
    超聲處理對顯微組分浮選分離及熱解影響研究
    煤炭工程(2019年11期)2019-12-13 07:38:10
    火紅
    不同變質(zhì)程度煉焦煤的鏡質(zhì)組特性研究
    木屑水蒸氣氣化制取富氫燃氣研究
    三塘湖盆地條湖組烴源巖地化特征及致密油油源對比
    91老司机精品| 久久人人爽av亚洲精品天堂| 午夜免费激情av| 91精品三级在线观看| 亚洲成人国产一区在线观看| 欧美黑人精品巨大| 在线观看www视频免费| 国内精品久久久久精免费| 国产亚洲欧美精品永久| 90打野战视频偷拍视频| 18禁国产床啪视频网站| 人人妻,人人澡人人爽秒播| 欧美激情久久久久久爽电影 | 亚洲专区字幕在线| 在线观看舔阴道视频| 久久草成人影院| 美女高潮喷水抽搐中文字幕| 欧美一区二区精品小视频在线| 变态另类成人亚洲欧美熟女 | 国产亚洲欧美在线一区二区| 天天添夜夜摸| 国产三级黄色录像| 欧美午夜高清在线| 日韩欧美一区视频在线观看| 一级a爱片免费观看的视频| 亚洲精品中文字幕在线视频| 国产亚洲精品一区二区www| 亚洲国产欧美网| 亚洲专区字幕在线| 侵犯人妻中文字幕一二三四区| 亚洲av片天天在线观看| 亚洲人成电影观看| 性色av乱码一区二区三区2| 久久精品国产亚洲av高清一级| 欧美黑人欧美精品刺激| 国产私拍福利视频在线观看| 国产精品98久久久久久宅男小说| 夜夜躁狠狠躁天天躁| 成人亚洲精品av一区二区| 女警被强在线播放| 丝袜人妻中文字幕| 国产高清激情床上av| 99国产精品99久久久久| 精品午夜福利视频在线观看一区| 久久香蕉激情| 天天添夜夜摸| 欧美午夜高清在线| 亚洲精品在线观看二区| 一级a爱视频在线免费观看| 久久 成人 亚洲| 免费久久久久久久精品成人欧美视频| 久9热在线精品视频| 一夜夜www| 长腿黑丝高跟| 国产野战对白在线观看| 亚洲片人在线观看| 亚洲国产毛片av蜜桃av| 午夜两性在线视频| 好看av亚洲va欧美ⅴa在| 级片在线观看| www.精华液| 国产色视频综合| 中国美女看黄片| 成人手机av| 国产精品美女特级片免费视频播放器 | 真人一进一出gif抽搐免费| a级毛片在线看网站| 久久久久久久久久久久大奶| 亚洲人成电影观看| 国产一区二区三区在线臀色熟女| 国产精品久久久久久精品电影 | 亚洲国产看品久久| 国产精品久久久人人做人人爽| 国产精品久久久久久亚洲av鲁大| 黑人操中国人逼视频| 窝窝影院91人妻| 九色国产91popny在线| 久久久久久久久中文| 美国免费a级毛片| 欧美成人午夜精品| 国语自产精品视频在线第100页| 亚洲国产精品久久男人天堂| 亚洲avbb在线观看| 国产97色在线日韩免费| 亚洲成国产人片在线观看| 亚洲第一电影网av| 日日夜夜操网爽| 精品国产国语对白av| 欧美 亚洲 国产 日韩一| 精品久久久久久成人av| a在线观看视频网站| 亚洲黑人精品在线| 欧美在线一区亚洲| cao死你这个sao货| 久久久久久国产a免费观看| 亚洲 欧美一区二区三区| 91麻豆av在线| 三级毛片av免费| 亚洲一卡2卡3卡4卡5卡精品中文| 十八禁人妻一区二区| 欧美成人一区二区免费高清观看 | 中文字幕人成人乱码亚洲影| 成人特级黄色片久久久久久久| 日韩国内少妇激情av| 精品国内亚洲2022精品成人| 日本三级黄在线观看| 少妇的丰满在线观看| 精品无人区乱码1区二区| 真人一进一出gif抽搐免费| 久久精品91无色码中文字幕| 免费高清在线观看日韩| 亚洲第一青青草原| 久久婷婷人人爽人人干人人爱 | 国产熟女午夜一区二区三区| 最近最新免费中文字幕在线| 一夜夜www| 99国产综合亚洲精品| av有码第一页| 亚洲av成人av| 99国产精品免费福利视频| 免费少妇av软件| 国产男靠女视频免费网站| 男女午夜视频在线观看| 97人妻精品一区二区三区麻豆 | 国产精品秋霞免费鲁丝片| 国产成人影院久久av| 如日韩欧美国产精品一区二区三区| 国产麻豆69| 中国美女看黄片| 大香蕉久久成人网| 亚洲免费av在线视频| 久久热在线av| 99国产精品99久久久久| 亚洲国产精品合色在线| 91九色精品人成在线观看| 高清在线国产一区| 法律面前人人平等表现在哪些方面| 亚洲成a人片在线一区二区| 色综合欧美亚洲国产小说| 午夜免费鲁丝| 久久精品亚洲精品国产色婷小说| 9191精品国产免费久久| 久久国产亚洲av麻豆专区| 亚洲一区二区三区不卡视频| 亚洲va日本ⅴa欧美va伊人久久| 欧美乱码精品一区二区三区| 人人澡人人妻人| 日本三级黄在线观看| 午夜成年电影在线免费观看| 在线播放国产精品三级| 免费在线观看日本一区| 亚洲在线自拍视频| 国产1区2区3区精品| 精品高清国产在线一区| 成人精品一区二区免费| 亚洲中文日韩欧美视频| 丝袜美足系列| 99久久99久久久精品蜜桃| 亚洲国产日韩欧美精品在线观看 | 1024视频免费在线观看| 美女 人体艺术 gogo| 欧美久久黑人一区二区| 美女国产高潮福利片在线看| 最好的美女福利视频网| 国产精品久久久久久精品电影 | 亚洲无线在线观看| 欧美色视频一区免费| 国产精品亚洲av一区麻豆| 久久热在线av| av免费在线观看网站| 99久久久亚洲精品蜜臀av| 久久精品aⅴ一区二区三区四区| 一二三四社区在线视频社区8| 久久久久久国产a免费观看| 天天添夜夜摸| 亚洲欧美一区二区三区黑人| 欧美人与性动交α欧美精品济南到| 法律面前人人平等表现在哪些方面| 欧美最黄视频在线播放免费| 又黄又粗又硬又大视频| 麻豆久久精品国产亚洲av| 精品人妻在线不人妻| 亚洲九九香蕉| 亚洲五月婷婷丁香| 精品国产一区二区久久| 一卡2卡三卡四卡精品乱码亚洲| 91老司机精品| www日本在线高清视频| 999久久久国产精品视频| 亚洲激情在线av| 国产成人免费无遮挡视频| 国产aⅴ精品一区二区三区波| 脱女人内裤的视频| 国产精品98久久久久久宅男小说| 一区福利在线观看| 成人av一区二区三区在线看| 亚洲精品在线观看二区| 亚洲熟妇中文字幕五十中出| 日本精品一区二区三区蜜桃| 亚洲国产看品久久| 久久久国产欧美日韩av| 丁香欧美五月| 亚洲精品一区av在线观看| 免费看a级黄色片| 久久亚洲精品不卡| 一区二区日韩欧美中文字幕| 精品乱码久久久久久99久播| 亚洲成av片中文字幕在线观看| 午夜福利免费观看在线| 亚洲av第一区精品v没综合| 免费人成视频x8x8入口观看| 在线观看午夜福利视频| www.www免费av| 99精品久久久久人妻精品| 亚洲成a人片在线一区二区| 亚洲成a人片在线一区二区| 日韩欧美在线二视频| 国产高清videossex| 国产99久久九九免费精品| 丝袜美足系列| 久久久久国产精品人妻aⅴ院| 老汉色∧v一级毛片| 国产精品久久久久久亚洲av鲁大| 成人三级黄色视频| 精品国产美女av久久久久小说| 国产精品免费一区二区三区在线| АⅤ资源中文在线天堂| 国产一区二区三区综合在线观看| 嫩草影院精品99| 亚洲av第一区精品v没综合| 一级,二级,三级黄色视频| 亚洲人成电影观看| 波多野结衣一区麻豆| 日本撒尿小便嘘嘘汇集6| 精品日产1卡2卡| 亚洲av电影不卡..在线观看| 国产在线精品亚洲第一网站| 亚洲av片天天在线观看| 岛国在线观看网站| 亚洲精品在线美女| 亚洲精品在线观看二区| 黄色女人牲交| videosex国产| 最近最新免费中文字幕在线| 一进一出抽搐动态| 色综合欧美亚洲国产小说| 亚洲情色 制服丝袜| 亚洲一区二区三区不卡视频| 亚洲成人免费电影在线观看| 久久欧美精品欧美久久欧美| 欧美黑人精品巨大| АⅤ资源中文在线天堂| 大型av网站在线播放| 成人手机av| 一进一出抽搐动态| 亚洲第一电影网av| 国产成人精品久久二区二区免费| 欧美 亚洲 国产 日韩一| 激情视频va一区二区三区| 婷婷精品国产亚洲av在线| 美女午夜性视频免费| 亚洲精品美女久久av网站| 性色av乱码一区二区三区2| 中亚洲国语对白在线视频| 午夜精品久久久久久毛片777| 99久久精品国产亚洲精品| 日本在线视频免费播放| 国产亚洲精品综合一区在线观看 | 亚洲成人精品中文字幕电影| 欧美黑人精品巨大| 免费高清在线观看日韩| 一级,二级,三级黄色视频| 日本一区二区免费在线视频| 老司机午夜福利在线观看视频| 91九色精品人成在线观看| 女人爽到高潮嗷嗷叫在线视频| 变态另类成人亚洲欧美熟女 | 久久久久国内视频| 久久精品国产99精品国产亚洲性色 | 一区福利在线观看| 国产精品日韩av在线免费观看 | 亚洲一区二区三区不卡视频| 亚洲精品国产区一区二| 国产精品久久久人人做人人爽| 国产三级黄色录像| 日本a在线网址| 久9热在线精品视频| 狂野欧美激情性xxxx| √禁漫天堂资源中文www| 真人做人爱边吃奶动态| 国产精品亚洲美女久久久| 亚洲专区国产一区二区| 欧美日本亚洲视频在线播放| www.999成人在线观看| 俄罗斯特黄特色一大片| av有码第一页| 黄色视频,在线免费观看| 黄片大片在线免费观看| 久久精品国产99精品国产亚洲性色 | 在线观看免费午夜福利视频| 99久久精品国产亚洲精品| 成熟少妇高潮喷水视频| 欧美一级a爱片免费观看看 | 99国产精品免费福利视频| 亚洲,欧美精品.| 老汉色av国产亚洲站长工具| 叶爱在线成人免费视频播放| 制服诱惑二区| 校园春色视频在线观看| 亚洲国产精品合色在线| 日韩成人在线观看一区二区三区| 18禁国产床啪视频网站| 欧美日本中文国产一区发布| 免费高清视频大片| 久久精品人人爽人人爽视色| 午夜福利在线观看吧| 美女 人体艺术 gogo| 极品人妻少妇av视频| bbb黄色大片| 亚洲人成伊人成综合网2020| 麻豆国产av国片精品| 一卡2卡三卡四卡精品乱码亚洲| www.精华液| 高潮久久久久久久久久久不卡| 亚洲精品av麻豆狂野| 亚洲一区二区三区不卡视频| 亚洲三区欧美一区| 国产精品国产高清国产av| 国产成+人综合+亚洲专区| 亚洲情色 制服丝袜| www.熟女人妻精品国产| 久久这里只有精品19| 日日干狠狠操夜夜爽| 亚洲精华国产精华精| 一区在线观看完整版| 在线观看日韩欧美| 高清在线国产一区| 手机成人av网站| 亚洲国产精品sss在线观看| 久久久久久大精品| 欧美乱妇无乱码| 欧美国产日韩亚洲一区| 欧美精品亚洲一区二区| 精品人妻1区二区| 国产三级黄色录像| 亚洲五月婷婷丁香| 欧美国产日韩亚洲一区| 欧美绝顶高潮抽搐喷水| 久久人人97超碰香蕉20202| 90打野战视频偷拍视频| 男人操女人黄网站| 黄片播放在线免费| 亚洲成人免费电影在线观看| 国产91精品成人一区二区三区| 一本久久中文字幕| 最新美女视频免费是黄的| 又紧又爽又黄一区二区| 在线观看一区二区三区| 亚洲精品中文字幕一二三四区| 久久久久国产一级毛片高清牌| netflix在线观看网站| 精品国产超薄肉色丝袜足j| 亚洲熟妇中文字幕五十中出| 香蕉国产在线看| 我的亚洲天堂| 亚洲国产欧美日韩在线播放| 日本精品一区二区三区蜜桃| 一级片免费观看大全| 女人被狂操c到高潮| 超碰成人久久| 美女免费视频网站| 国内毛片毛片毛片毛片毛片| 国产区一区二久久| 美国免费a级毛片| 无遮挡黄片免费观看| 91成年电影在线观看| 日韩一卡2卡3卡4卡2021年| 久久久久久久久中文| 午夜视频精品福利| 18禁美女被吸乳视频| 欧美成狂野欧美在线观看| 国内久久婷婷六月综合欲色啪| 9热在线视频观看99| 日韩欧美一区视频在线观看| 一个人观看的视频www高清免费观看 | 国产成人av教育| 国产精品电影一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 日韩欧美一区视频在线观看| 精品国产一区二区三区四区第35| 国产欧美日韩精品亚洲av| 桃红色精品国产亚洲av| 9色porny在线观看| 成人特级黄色片久久久久久久| 国产一区在线观看成人免费| 夜夜爽天天搞| 人妻丰满熟妇av一区二区三区| 国产免费男女视频| 91字幕亚洲| 免费高清视频大片| 给我免费播放毛片高清在线观看| 女人精品久久久久毛片| 一级毛片精品| 国产精品99久久99久久久不卡| 亚洲片人在线观看| or卡值多少钱| 久久亚洲真实| 大码成人一级视频| 久久精品亚洲熟妇少妇任你| 午夜视频精品福利| 欧美老熟妇乱子伦牲交| 亚洲欧美一区二区三区黑人| 两人在一起打扑克的视频| 亚洲av成人一区二区三| 成人av一区二区三区在线看| 午夜免费观看网址| 国产精品 国内视频| 女人精品久久久久毛片| 麻豆一二三区av精品| 最近最新中文字幕大全电影3 | www.精华液| 亚洲伊人色综图| 亚洲人成电影观看| 女性被躁到高潮视频| 国产精品免费一区二区三区在线| 麻豆av在线久日| 精品久久久久久成人av| 一本大道久久a久久精品| 国产精品免费一区二区三区在线| 免费不卡黄色视频| 黄网站色视频无遮挡免费观看| 极品教师在线免费播放| 久久热在线av| 亚洲国产看品久久| 无人区码免费观看不卡| 在线观看66精品国产| 久久久久久大精品| 最新美女视频免费是黄的| 亚洲国产精品sss在线观看| 别揉我奶头~嗯~啊~动态视频| 久久久久国产精品人妻aⅴ院| 久久久久久大精品| 91大片在线观看| 午夜福利欧美成人| 在线观看免费午夜福利视频| 精品不卡国产一区二区三区| 最近最新免费中文字幕在线| 桃色一区二区三区在线观看| 女人爽到高潮嗷嗷叫在线视频| 在线视频色国产色| 亚洲黑人精品在线| 人妻丰满熟妇av一区二区三区| av天堂在线播放| 久久久久国内视频| 精品久久久久久久人妻蜜臀av | 国产精品久久久久久精品电影 | 国产精品 国内视频| 一级毛片女人18水好多| 乱人伦中国视频| 极品教师在线免费播放| 色尼玛亚洲综合影院| 亚洲成人久久性| 国产精品av久久久久免费| 欧美大码av| 国产av在哪里看| 亚洲国产欧美一区二区综合| 久久久久久久久免费视频了| 亚洲伊人色综图| 99国产精品一区二区蜜桃av| 国产精品永久免费网站| 亚洲 国产 在线| 日本a在线网址| 我的亚洲天堂| 日本欧美视频一区| 久久九九热精品免费| 成人永久免费在线观看视频| 国产私拍福利视频在线观看| 免费不卡黄色视频| 一边摸一边做爽爽视频免费| 国产一区二区三区在线臀色熟女| 国产av一区在线观看免费| 人人妻,人人澡人人爽秒播| 中国美女看黄片| 9191精品国产免费久久| 男女做爰动态图高潮gif福利片 | www.熟女人妻精品国产| 国产精品综合久久久久久久免费 | 禁无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 很黄的视频免费| 一边摸一边做爽爽视频免费| 亚洲中文av在线| 亚洲avbb在线观看| 亚洲 国产 在线| 久久久久国内视频| 香蕉久久夜色| 免费高清视频大片| 天堂影院成人在线观看| 久热这里只有精品99| 国产单亲对白刺激| 一区二区三区高清视频在线| 国产精品 国内视频| 丝袜美足系列| 国产精品一区二区三区四区久久 | 亚洲人成77777在线视频| 人人澡人人妻人| 国产蜜桃级精品一区二区三区| 成人国语在线视频| 久久人人97超碰香蕉20202| 久久九九热精品免费| 9热在线视频观看99| 丝袜美腿诱惑在线| 两个人看的免费小视频| 日韩精品免费视频一区二区三区| 免费不卡黄色视频| 一区二区三区高清视频在线| 国产熟女午夜一区二区三区| 久久久久国产精品人妻aⅴ院| 久久久国产精品麻豆| 最近最新中文字幕大全电影3 | 一级毛片精品| 午夜福利,免费看| 国产亚洲欧美精品永久| 久久欧美精品欧美久久欧美| 嫩草影院精品99| 国产精品久久电影中文字幕| 日韩成人在线观看一区二区三区| 黄频高清免费视频| 国产精品 欧美亚洲| 天天躁狠狠躁夜夜躁狠狠躁| bbb黄色大片| 亚洲精品在线美女| av网站免费在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 十八禁网站免费在线| 午夜福利成人在线免费观看| 国产精品久久久久久亚洲av鲁大| 欧美色欧美亚洲另类二区 | 国产成人一区二区三区免费视频网站| 日韩免费av在线播放| 精品无人区乱码1区二区| 自线自在国产av| 久久这里只有精品19| 精品久久久久久,| 国产激情久久老熟女| 他把我摸到了高潮在线观看| 91大片在线观看| 啦啦啦 在线观看视频| 在线国产一区二区在线| 夜夜爽天天搞| 女同久久另类99精品国产91| 真人一进一出gif抽搐免费| 大型黄色视频在线免费观看| 免费在线观看黄色视频的| 精品国产亚洲在线| 好男人在线观看高清免费视频 | 久热爱精品视频在线9| 1024视频免费在线观看| av免费在线观看网站| 啪啪无遮挡十八禁网站| 久久久久国产精品人妻aⅴ院| 免费女性裸体啪啪无遮挡网站| 可以在线观看毛片的网站| 欧美午夜高清在线| 多毛熟女@视频| 丰满的人妻完整版| 成年女人毛片免费观看观看9| 国产高清视频在线播放一区| 真人一进一出gif抽搐免费| 在线观看www视频免费| 日韩高清综合在线| 成人特级黄色片久久久久久久| 午夜福利视频1000在线观看 | 狠狠狠狠99中文字幕| 国产伦一二天堂av在线观看| 黑丝袜美女国产一区| 国产精品亚洲美女久久久| 日韩欧美免费精品| 一本综合久久免费| 少妇裸体淫交视频免费看高清 | av天堂在线播放| 日本撒尿小便嘘嘘汇集6| 国产野战对白在线观看| 国产单亲对白刺激| 精品国产一区二区三区四区第35| 午夜福利成人在线免费观看| 亚洲精品久久成人aⅴ小说| 在线播放国产精品三级| 美女高潮喷水抽搐中文字幕| 男男h啪啪无遮挡| 久久中文字幕人妻熟女| 韩国av一区二区三区四区| 免费高清在线观看日韩| 咕卡用的链子| 国产精品久久久久久人妻精品电影| 亚洲欧美精品综合一区二区三区| 久久精品91无色码中文字幕| 亚洲国产高清在线一区二区三 | 88av欧美| 国产精品野战在线观看| 天天添夜夜摸| 国产亚洲精品一区二区www| 久久精品国产亚洲av高清一级| 国产精品99久久99久久久不卡| 国产亚洲av嫩草精品影院| 成人av一区二区三区在线看| 国产成人一区二区三区免费视频网站| 国产野战对白在线观看| 午夜福利成人在线免费观看| 精品国产美女av久久久久小说| 人人澡人人妻人| 久久久久久国产a免费观看| 波多野结衣巨乳人妻| 性色av乱码一区二区三区2| 757午夜福利合集在线观看|