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

    反演純剪切模型中β的算法研究

    2011-10-17 08:39:42宋海斌
    物探化探計(jì)算技術(shù) 2011年1期
    關(guān)鍵詞:陸坡巖石圈測線

    張 炯,宋海斌

    (1.遼寧省冶金地質(zhì)勘查局 地質(zhì)勘查研究院,鞍山 114038;

    2.中國科學(xué)院地質(zhì)與地球物理研究所 油氣資源研究院重點(diǎn)實(shí)驗(yàn)室,北京 100029)

    反演純剪切模型中β的算法研究

    張 炯1,宋海斌2

    (1.遼寧省冶金地質(zhì)勘查局 地質(zhì)勘查研究院,鞍山 114038;

    2.中國科學(xué)院地質(zhì)與地球物理研究所 油氣資源研究院重點(diǎn)實(shí)驗(yàn)室,北京 100029)

    以往反演伸展因子的算法多是基于單井資料進(jìn)行擬合計(jì)算得到,很少把構(gòu)造與熱演化史結(jié)合研究。這里利用M cKenzie在1978年提出的純剪切模型中計(jì)算沉降和熱流的公式,基于一維反演伸展因子的思想,采用M atlab軟件GU I功能,編制多點(diǎn)反演伸展因子界面,大大減少在對偽井?dāng)?shù)據(jù)反演伸展因子時的工作量。以南海中北部SO49-25測線為例,計(jì)算了該剖面陸架和陸坡區(qū)的伸展因子和古熱流值,并驗(yàn)證了算法的可行性。

    純剪切模型;伸展因子;反演

    0 前言

    隨著大量資料積累及電子技術(shù)的進(jìn)步和發(fā)展,通過地質(zhì)過程定量分析或模擬研究,加深了對地質(zhì)過程的認(rèn)識,因此盆地定量模擬研究對解決地質(zhì)問題的定量化提供了依據(jù)。其中地史和熱史的數(shù)值模擬最終目的,是為油氣生、排、運(yùn)聚、成藏模擬創(chuàng)造條件。盆地構(gòu)造~地史和熱演化史,是油氣資源評價和油氣成藏動力學(xué)研究的重要內(nèi)容之一。

    二十世紀(jì)七十年代以來,許多學(xué)者對巖石圈拉張、盆地和大陸邊緣的形成演化,進(jìn)行了廣泛深入的研究,提出了多種地球動力學(xué)定量模型。其中以M cKenzie[4]提出的巖石圈純剪切模型最為經(jīng)典,為以后定量研究盆地與張裂大陸邊緣形成演化機(jī)制奠定了基礎(chǔ)。數(shù)值模擬方法必須與觀測資料相結(jié)合,只有這樣模擬結(jié)果才具有地質(zhì)意義,才能為解釋地質(zhì)現(xiàn)象提供科學(xué)依據(jù),換言之,數(shù)值模型需要以觀測數(shù)據(jù)作為約束條件[1~4]。

    1 純剪切模型

    M cKenzie[4]提出了巖石圈均勻拉張模型,又稱純剪切模型。該模型描述了巖石圈對拉伸作用的基本響應(yīng),建立了拉伸系數(shù)(β)和盆地?zé)崃髦g的定量關(guān)系,奠定了拉伸盆地定量模擬的理論基礎(chǔ)。在盆地模擬和沉積盆地?zé)崾坊謴?fù)領(lǐng)域中,具有劃時代意義。構(gòu)造~熱演化模擬方法的發(fā)展,實(shí)際上即為拉伸盆地定量模擬的地質(zhì)地球物理模型的發(fā)展。M cKenzie提出的定量模型屬于瞬時均勻純剪切模型,如下頁圖1所示,即:拉伸時間短(<20M a),拉伸期間的熱擴(kuò)散忽略不計(jì),地溫梯度瞬時增高,巖石圈變形機(jī)制為純剪切變形,拉伸過程中地殼和巖石圈地幔均勻變薄。

    M cKenzie認(rèn)為,裂陷盆地總體沉降由初始沉降和熱沉降二部份組成:①伸展期間由斷層控制的沉降,取決于地殼的初始厚度和伸展系數(shù)β值的大小;②伸展之后的熱沉降,由巖石圈冷卻所引起,取決于伸展量。并認(rèn)為巖石圈伸展是均勻的,斷層控制的沉降是瞬間完成的,而熱沉降速率是隨時間呈指數(shù)衰減的。M cKenzie提出的純剪切模型的基本公式以及各個參數(shù)的意義為:

    圖1 純剪切模型圖(M ckenzie,1978)Fig.1 Pure shearmodel(M cKenzie,1978)

    (1)初始沉降公式:

    (2)熱沉降公式:

    (3)溫度異常采用傅里葉級數(shù)表達(dá):

    其中 an為實(shí)常數(shù);z為深度。(4)熱流值公式:

    以上公式(1)~式(5)中的參數(shù)意義:β是伸展因子;yL是初始的巖石圈厚度;yc是初始的地殼厚度;ρm*是0℃時地幔的密度,通常取值3 330 kg/m3;ρ*c是0℃時地殼的密度,通常取值2 800 kg/m3;ρs是沉積物或水的密度,通常取值1 030 kg/m3;αv是熱擴(kuò)散系數(shù),通常取值3.28×10-5℃;Tm是軟流圈的溫度,通常取值1 333℃;κ是熱膨脹系數(shù),通常取值0.007 5 cm2/s。

    研究分析M cKenzie均勻伸展模型的理論結(jié)果,認(rèn)為:

    (1)M cKenzie均勻伸展模型是一維瞬時拉張模型,其總沉積量由二部份組成:斷層控制的初始沉降和沉降速率隨指數(shù)衰減的熱沉降。熱沉降隨伸展因子β的增大而增大,β為無窮大,代表洋中脊的情況。

    (2)伸展因子β為巖石圈伸展前厚度與伸展后厚度之比,因此β越大,巖石圈伸展減薄量越大,初始沉降也就越大,因此初始熱流異常也就越大。

    (3)M cKenzie模型模擬的是地表熱流,巖石圈拉張瞬間完成時熱流達(dá)到最大值,隨即呈指數(shù)衰減。伸展因子越大,初始熱流越高,在經(jīng)過200M a冷卻后,熱流趨于穩(wěn)態(tài)。

    通過了解M cKenzie均勻伸展模型的機(jī)制,可以得出在利用M cKenzie均勻伸展模型反演計(jì)算伸展因子的計(jì)算過程中,初始沉降量取決于地殼與巖石圈厚度之比(yc/yL)和伸展因子(β)。如果對于某個盆地yc/yL已知,則理論上由斷層控制的初始沉降可用來估算伸展量。由于其它因素造成了裂谷期厚度的變化,使得這種方法可信度降低。因此在擬合伸展因子的計(jì)算過程中,很少利用考慮與初始沉降的關(guān)系。因?yàn)镸 cKenzie模型預(yù)測的熱沉降只取決于伸展因子,影響裂后熱沉降的因素相對較少,所以通過M cKenzie模型計(jì)算熱沉降的公式,可以估算伸展因子。擬合方法是把裂后沉積厚度經(jīng)過沉積物壓實(shí)校正,在沉積物負(fù)載校正,水深校正,以及重力均衡校正之后,得到構(gòu)造沉降史。使得裂后沉降部份與理論模型預(yù)測的熱沉降曲線進(jìn)行對比,得到最佳擬合的伸展因子。

    因此作者利用M cKenzie理論模型一維擬合伸展因子的思想,編制了多點(diǎn)擬合伸展因子的帶界面的模擬程序。

    2 伸展因子估算方法

    A llen等[1]總結(jié)了多種估算拉伸因子的方法,作者在本文中討論的是利用沉降分析估算β。首先利用地震剖面和鉆井資料,通過回剝法(Steck ler和W atts[5])可以得到水載構(gòu)造沉降,即觀測沉降。從純剪切模型的初始沉降公式中可以看到,初始沉降取決于初始地殼厚度與巖石圈厚度之比和拉伸因子。因此,在已知初始地殼與巖石圈厚度的情況下,理論上可以通過斷層控制的初始沉降估算拉伸因子。但由于裂谷期厚度的劇烈變化,使得估算結(jié)果變得不可靠。常用的方法是給定不同的拉伸因子β,利用純剪切模型的熱沉降公式,計(jì)算理論的水載熱沉降曲線,通過比較模型預(yù)測的沉降與觀測沉降,以確定最佳擬合的β。

    基于均勻伸展模型正演模擬研究中的計(jì)算公式,編輯擬合伸展因子的反演算法,步驟如下:

    (1)已知m個時期的基底的埋藏史數(shù)據(jù),并對數(shù)據(jù)進(jìn)行轉(zhuǎn)換,得到構(gòu)造沉降史數(shù)據(jù),其中包括裂谷期沉降(Ys)和裂后期沉降(Yt(t))。

    (2)已知裂后期開始的時間,確定裂后期熱沉降史數(shù)據(jù)有n個,從構(gòu)造沉降史數(shù)據(jù)中減去裂谷期沉降(Ys)的貢獻(xiàn),即得到裂后期熱沉降量(Yt(t)-Ys)。

    (3)基于M cKenzie均勻伸展模型中的計(jì)算熱沉降量的公式:),以β=1.1為初始值,計(jì)算yt(t)。然后通過迭代運(yùn)算,使得為最小,即得到擬合因子β,其中err是單個構(gòu)造沉降量的擬合計(jì)算的計(jì)算誤差,單位為“m”。

    (4)當(dāng)β>6.1時,伸展因子取6.1。

    3 多點(diǎn)伸展因子擬合界面研究

    采用m atlab語言進(jìn)行編程,編制反演多點(diǎn)伸展因子擬合的程序,并利用m atlab軟件的GU I相關(guān)功能,實(shí)現(xiàn)PC機(jī)可視化編程,編輯獨(dú)立的計(jì)算機(jī)界面,用于模型的數(shù)值模擬工作。在PC機(jī)中有m atlab動態(tài)鏈接庫的前提下,可以實(shí)現(xiàn)程序脫離m atlab軟件環(huán)境運(yùn)行[6~9]。

    圖2為主界面,“顯示路徑”鍵:可以顯示參數(shù)輸入的路徑;“運(yùn)行”鍵:即輸出運(yùn)算結(jié)果。

    圖3為時間參數(shù)文件界面,在“tim e.par”文件中,第一行中“8”代表本次參與運(yùn)算的地層時間有八個;第二行分別顯示各個時間點(diǎn);第三行中“23.3”代表裂后熱沉降開始的時間。

    圖4為基本參數(shù)文件界面,在“param eters.par”文件中,“125000”代表原始巖石圈厚度,單位為“m”;“31200”代表原始地殼厚度,單位為“m”;“3330”代表原始地幔密度,單位為“kg/m3”;“2800”代表原始地殼密度,單位為“kg/m3”;“1030”代表海水密度,單位為“kg/m3”;“3”代表熱導(dǎo)率值,單位“J/C/m/s”;“3.28E-5”代表熱擴(kuò)散系數(shù),單位“℃”;“1330”代表軟流圈頂界溫度,單位為“℃”;“1.0E-6”代表熱膨脹系數(shù),單位“cm2/s”。

    圖2 設(shè)置參數(shù)路徑Fig.2 Set the param eterspath

    圖3 時間參數(shù)文件Fig.3 The file of tim e param eters

    圖4 計(jì)算基本參數(shù)文件Fig.4 The file of calculating the elem entary param eters

    圖5(見下頁)為沉降數(shù)據(jù)參數(shù)文件界面,在“subsidence.par”文件中,第一行數(shù)據(jù)中“8”代表有八個層位數(shù)據(jù),單位“m”,“120”代表有一百二十個點(diǎn)數(shù)據(jù)。

    圖6(見下頁)為密度參數(shù)文件界面,在“rho_Sed.par”中,第一行數(shù)據(jù)“8”對應(yīng)八個沉積物平均密度,單位“kg/m3”,“120”代表有一百二十個點(diǎn)。

    圖7和圖8為輸出文件界面,其中“heat_flux.dat”文件和“beta_inv.dat”文件是輸出結(jié)果文件。“heat_flux.dat”文件中,第一行數(shù)據(jù)“120”是代表本次參與擬合伸展因子的點(diǎn)數(shù)為一百二十個,“6”代表是輸出六個時間點(diǎn)的熱流;“beta_inv.dat”文件中,第一行數(shù)據(jù)“120”是代表本次參與擬合伸展因子的結(jié)果有一百二十個;從第二行數(shù)據(jù)開始,第一列數(shù)據(jù)顯示為伸展因子擬合結(jié)果,第二列數(shù)據(jù)顯示為計(jì)算擬合誤差。

    圖5 沉降數(shù)據(jù)參數(shù)文件Fig.5 The file of subsidence data

    圖6 沉積物平均密度參數(shù)文件Fig.6 The file of deposition average density

    4 SO 49-25測線實(shí)例研究

    作者利用在2D-M ove軟件中分別輸出各個時期的基底沉降,已經(jīng)完成地震剖面的埋藏史恢復(fù)工作,只需提取垂直方向各個時期的沉降數(shù)據(jù),對沉降數(shù)據(jù)做A iry均衡校正,轉(zhuǎn)化為基底的構(gòu)造沉降。根據(jù)有效的多點(diǎn)伸展因子擬合的算法,計(jì)算不同區(qū)域的伸展因子的變化,并計(jì)算當(dāng)中存在的誤差分析。圖9(見下頁)中顯示的計(jì)算誤差,是擬合伸展因子過程中的最小err[5]。

    圖7 輸出的熱流結(jié)果文件Fig.7 The file of exporting the heat flow resu lt

    圖8 輸出的伸展因子及擬合誤差文件Fig.8 The file of exporting the stretching factor and fitting error

    25測線北部陸架部份的伸展因子擬合的結(jié)果顯示(見下頁圖9),在測線0 km~20 km處伸展量相對較小,伸展因子在1.7~2.5之間,而且計(jì)算誤差多在20~90之間,分析認(rèn)為是因該區(qū)域位于珠二坳陷,裂后沉降較大,這說明伸展因子在擬合過程中有一定的合理性。在20 km~80 km處伸展量相對較大,伸展因子變化較小導(dǎo)致,在1.3~2.0之間,而且計(jì)算誤差在20~110之間,分析認(rèn)為是因該區(qū)域位于一統(tǒng)暗沙隆起區(qū),相對裂后沉降較小導(dǎo)致。在計(jì)算西沙海槽伸展因子結(jié)果時,因發(fā)生裂后沉降結(jié)果較大,導(dǎo)致伸展因子結(jié)果在2.0~6.1之間變化較大,且計(jì)算誤差在50~250之間變化,導(dǎo)致結(jié)果可信度不高。

    SO49-25測線北部陸坡區(qū)熱流分析(見后面圖10):測線穿過北部陸架~西北海槽~北部陸坡區(qū),西沙海槽和陸坡區(qū)熱流值較大,北部陸架區(qū)0 km~30 km熱流值總體呈增大趨勢;在30 km~40 km處,熱流有減小趨勢,分析認(rèn)為該區(qū)域位于一統(tǒng)暗沙隆起區(qū),使得熱流值較小;在40 km~80 km處,總體為增大趨勢,符合靠近洋盆熱流增大的趨勢。

    圖9 SO49~25北部陸架區(qū)和西沙海槽~陸坡區(qū)伸展因子分布圖Fig.9 The stretching factor distribution in the northern continental shelf and X isha trough area of the SO49-25 line

    對各個時期的熱流值分布的認(rèn)識:

    (1)23.8M a屬于裂后期開始的時間,熱流值最大,可達(dá)到80mW/m2,均在60mW/m2以上。

    (2)16.5 M a熱流逐漸減小,最大值為70mW/m2,均值在58mW/m2以上。

    (3)10.5 M a熱流繼續(xù)減小,最大值為68mW/m2,均在55mW/m2以上。

    (4)現(xiàn)今熱流仍然在減小,最大值為60mW/m2,均在50mW/m2以上。

    在以上熱流計(jì)算中,暫時沒有考慮放射性元素的影響。熱流Q=qh+q+Δq,q代表地殼淺層、深地殼和地幔熱源的熱流;Δq代表異常熱流;qh是來自放射性衰減的熱流。

    自23.8M a以來,測線所在研究區(qū)處于熱沉降階段,熱流從裂陷谷期末期逐漸下降。

    根據(jù)模擬計(jì)算,17測線陸坡區(qū)現(xiàn)今熱流在50mW/m2~62 mW/m2之間,加上放射性生熱貢獻(xiàn)(簡單取28 mW/m2),在78 mW/m2~90mW/m2之間;18測線陸架陸坡部份現(xiàn)今熱流主要在40mW/m2~60mW/m2之間,加上放射性生熱貢獻(xiàn),在68 mW/m2~88 mW/m2之間;25測線陸坡部份在40mW/m2~58 mW/m2之間,加上放射性生熱貢獻(xiàn),在68mW/m2~86mW/m2之間。自東向西,拉伸因子總體減小,熱流相應(yīng)降低。

    5 結(jié)論

    作者在本文重點(diǎn)研究了M atlab軟件的GU I功能,結(jié)合M cKenzie提出的純剪切模型,主要完成以下幾點(diǎn)工作:①編制多點(diǎn)反演伸展因子的帶界面的程序,減少了在對偽井?dāng)?shù)據(jù)反演伸展因子時的工作量;②以南海中北部SO49-25測線為例,計(jì)算了該剖面陸架和陸坡區(qū)的伸展因子和古熱流值,把構(gòu)造史和熱演化史有效地結(jié)合起來,為下一步討論油氣資源的分布提供參考。

    致謝 作者在本文所用SO航次的地震數(shù)據(jù)是由德國BGR的Franke博士提供,在此向Franke博士表示感謝。

    [1]ALLEN P A,ALLEN J R.Basin analysis,p rincip les and app lications[M].Oxford:B lackwell Scientific Publications,1990.

    [2]JARV ISG T,MCKENZIED P.Sedim entary basin form ationw ith finite extension rates[J].Earth and Planetary Science Letters,1980,48(1):42.

    [3]KEEN C E,DEHLER SA.Stretching and subsidence:R ifting of con jugatem argins in theNorth A tlantic region[J].Tectonics,1993,12:1219.

    [4]MCKENZIE D.Som e rem arks on the developm ent of sedim entary basins[J].Earth and Planetary Science Letters,1978,40(1):25.

    [5]STECKLERM S,WATTSA B.Subsidence of the atlan tic-type continentalm argin off New Yo rk[J].Earth and p lanetary sciences letters,1978,41:1.

    [6]李顯宏.M atlab 7.x界面設(shè)計(jì)與編譯技巧[M].北京:電子工業(yè)出版社,2006.

    [7]王正林,劉明.精通MATLAB7[M].北京:電子工業(yè)出版社,2006.

    [8]楊長青,胥澤銀.基于MATLAB的面積計(jì)算方法[J].物探化探計(jì)算技術(shù),2004,26(2):177.

    [9]李曉昌.在MATLAB平臺上實(shí)現(xiàn)可控源音頻大地電磁反演數(shù)據(jù)三維可視化顯示[J].物探化探計(jì)算技術(shù),2007,(增刊):69.

    TU 432

    A

    1001—1749(2011)01—0101—06

    國家基礎(chǔ)研究發(fā)展規(guī)劃項(xiàng)目資助(2007CB411704)

    2010-08-09 改回日期:2010-11-02

    張炯(1984-),男,碩士,從事應(yīng)用地球物理研究工作。

    猜你喜歡
    陸坡巖石圈測線
    極地海洋多波束測量測線布設(shè)系統(tǒng)設(shè)計(jì)及實(shí)現(xiàn)
    基于動態(tài)規(guī)劃的多波束測線布設(shè)模型
    楚科奇陸坡流的海表時空特征
    第四章 堅(jiān)硬的巖石圈
    應(yīng)用陸坡形態(tài)—遷移軌跡組合法尋找大型富砂深水扇
    巖石圈磁場異常變化與巖石圈結(jié)構(gòu)的關(guān)系
    地震研究(2017年3期)2017-11-06 21:54:14
    2014年魯?shù)?—5級地震相關(guān)斷裂的巖石圈磁異常分析
    地震研究(2017年3期)2017-11-06 01:58:51
    層序地層格架內(nèi)白云凹陷北坡珠江組沉積相分析
    Rio Muni盆地第四紀(jì)陸坡地震地貌學(xué)①②
    蘆山7.0級地震前后巖石圈磁場異常變化研究
    地震研究(2014年1期)2014-02-27 09:29:41
    美女 人体艺术 gogo| 好男人电影高清在线观看| 亚洲专区国产一区二区| 嫩草影视91久久| 亚洲精品日韩av片在线观看| 宅男免费午夜| 亚洲熟妇中文字幕五十中出| 国产伦人伦偷精品视频| 午夜精品在线福利| x7x7x7水蜜桃| 国产单亲对白刺激| 俺也久久电影网| 成年人黄色毛片网站| 午夜精品久久久久久毛片777| 美女高潮的动态| 99久国产av精品| 国产大屁股一区二区在线视频| 亚洲第一电影网av| 在线国产一区二区在线| 又爽又黄无遮挡网站| 日本在线视频免费播放| 国产成人啪精品午夜网站| 国产精华一区二区三区| 久9热在线精品视频| 99热6这里只有精品| 精品久久久久久成人av| 搡老熟女国产l中国老女人| 日本在线视频免费播放| 日韩大尺度精品在线看网址| 精品一区二区三区av网在线观看| 日日摸夜夜添夜夜添av毛片 | 国产野战对白在线观看| 精品午夜福利在线看| 国产男靠女视频免费网站| 久久精品91蜜桃| 激情在线观看视频在线高清| 亚洲第一区二区三区不卡| 国产野战对白在线观看| 女生性感内裤真人,穿戴方法视频| 狠狠狠狠99中文字幕| 成人特级av手机在线观看| 黄色女人牲交| 他把我摸到了高潮在线观看| 午夜福利免费观看在线| 桃红色精品国产亚洲av| 国产精品伦人一区二区| av专区在线播放| 国产色爽女视频免费观看| 精品99又大又爽又粗少妇毛片 | 成人国产综合亚洲| 国产色爽女视频免费观看| 丝袜美腿在线中文| 国产伦一二天堂av在线观看| 伦理电影大哥的女人| 成人特级av手机在线观看| 欧美高清成人免费视频www| 黄色视频,在线免费观看| 国产视频一区二区在线看| 青草久久国产| 51国产日韩欧美| 九色成人免费人妻av| 小说图片视频综合网站| 成人鲁丝片一二三区免费| 国产三级黄色录像| 两个人的视频大全免费| 国产伦一二天堂av在线观看| 看片在线看免费视频| 窝窝影院91人妻| 国产精品久久视频播放| 国产精品日韩av在线免费观看| 日韩欧美一区二区三区在线观看| 亚洲 国产 在线| 亚洲人成网站在线播放欧美日韩| netflix在线观看网站| 国产午夜精品论理片| 国产麻豆成人av免费视频| 亚洲av不卡在线观看| 国产激情偷乱视频一区二区| 少妇高潮的动态图| 男人舔奶头视频| 男女下面进入的视频免费午夜| 国产国拍精品亚洲av在线观看| 日本一二三区视频观看| 蜜桃久久精品国产亚洲av| 中文字幕av在线有码专区| .国产精品久久| 欧美又色又爽又黄视频| 亚洲av第一区精品v没综合| 亚洲国产色片| 最近最新中文字幕大全电影3| 日韩欧美在线乱码| 国产乱人视频| 国产精品三级大全| 亚洲最大成人手机在线| 精品午夜福利视频在线观看一区| 中文字幕av在线有码专区| 99久久九九国产精品国产免费| 好男人电影高清在线观看| 香蕉av资源在线| 久久精品影院6| 老鸭窝网址在线观看| 在线看三级毛片| 亚洲第一欧美日韩一区二区三区| 免费在线观看亚洲国产| 久久久久精品国产欧美久久久| 小蜜桃在线观看免费完整版高清| 国产爱豆传媒在线观看| 中文亚洲av片在线观看爽| 一级a爱片免费观看的视频| 国产一区二区亚洲精品在线观看| 1024手机看黄色片| 中文字幕久久专区| .国产精品久久| 欧美日本亚洲视频在线播放| 国产高清激情床上av| 久久精品综合一区二区三区| 无人区码免费观看不卡| 午夜视频国产福利| 久久精品综合一区二区三区| 性色avwww在线观看| 老司机午夜福利在线观看视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | a级一级毛片免费在线观看| 99久久九九国产精品国产免费| 午夜福利高清视频| 少妇人妻一区二区三区视频| 亚洲最大成人av| 国产亚洲av嫩草精品影院| 国产亚洲欧美98| 乱码一卡2卡4卡精品| 欧美最新免费一区二区三区 | 婷婷色综合大香蕉| 国产精品爽爽va在线观看网站| 69av精品久久久久久| 最新在线观看一区二区三区| 波多野结衣巨乳人妻| 中文字幕久久专区| 中文字幕高清在线视频| 国产高清三级在线| 国产亚洲精品av在线| 人妻丰满熟妇av一区二区三区| 欧美最黄视频在线播放免费| 久久久久久久久久成人| 日本黄色片子视频| 欧美乱妇无乱码| 国产精品日韩av在线免费观看| 亚洲av电影在线进入| 国产精品自产拍在线观看55亚洲| 欧美黑人欧美精品刺激| 精品人妻1区二区| 精品国内亚洲2022精品成人| 在线a可以看的网站| 搞女人的毛片| 韩国av一区二区三区四区| a级一级毛片免费在线观看| 国产激情偷乱视频一区二区| 国产精品野战在线观看| 国产激情偷乱视频一区二区| av在线老鸭窝| 日日夜夜操网爽| 丰满人妻熟妇乱又伦精品不卡| 亚洲美女搞黄在线观看 | 精品一区二区三区人妻视频| 欧美绝顶高潮抽搐喷水| 亚洲最大成人手机在线| 亚洲av免费高清在线观看| 两个人视频免费观看高清| 性色av乱码一区二区三区2| 亚洲狠狠婷婷综合久久图片| 亚洲性夜色夜夜综合| 90打野战视频偷拍视频| 色哟哟·www| 少妇的逼好多水| 五月玫瑰六月丁香| 一区福利在线观看| 国产精品,欧美在线| 女生性感内裤真人,穿戴方法视频| 最后的刺客免费高清国语| 国产亚洲欧美98| 国产在线男女| 神马国产精品三级电影在线观看| 国产乱人视频| a级毛片免费高清观看在线播放| 久久久国产成人免费| 午夜亚洲福利在线播放| 亚洲人成网站在线播| 欧美极品一区二区三区四区| 婷婷亚洲欧美| 又粗又爽又猛毛片免费看| h日本视频在线播放| 国产午夜精品久久久久久一区二区三区 | 丰满的人妻完整版| 亚洲精品粉嫩美女一区| 亚洲成a人片在线一区二区| 18+在线观看网站| 麻豆久久精品国产亚洲av| 久久久成人免费电影| 久99久视频精品免费| 国产伦精品一区二区三区四那| 国内精品美女久久久久久| 99热只有精品国产| 欧美成人免费av一区二区三区| 国内精品一区二区在线观看| 久久精品久久久久久噜噜老黄 | 国产一区二区激情短视频| 欧美色视频一区免费| 久久国产乱子伦精品免费另类| 亚洲成人中文字幕在线播放| 18禁在线播放成人免费| 国产真实乱freesex| 成人高潮视频无遮挡免费网站| 91av网一区二区| 两个人视频免费观看高清| 亚洲男人的天堂狠狠| 三级男女做爰猛烈吃奶摸视频| 男女做爰动态图高潮gif福利片| 亚洲精华国产精华精| 最新在线观看一区二区三区| 国产亚洲精品综合一区在线观看| 亚洲电影在线观看av| 男女床上黄色一级片免费看| 午夜福利视频1000在线观看| 亚洲av成人精品一区久久| 我的老师免费观看完整版| 久久久精品欧美日韩精品| 精品无人区乱码1区二区| 别揉我奶头~嗯~啊~动态视频| 午夜视频国产福利| 18禁黄网站禁片免费观看直播| 国产精品爽爽va在线观看网站| 黄色丝袜av网址大全| 午夜两性在线视频| 国产一区二区亚洲精品在线观看| 别揉我奶头 嗯啊视频| 三级男女做爰猛烈吃奶摸视频| 色综合亚洲欧美另类图片| 99久久精品国产亚洲精品| 91狼人影院| 午夜影院日韩av| 噜噜噜噜噜久久久久久91| 高清在线国产一区| 精品一区二区三区人妻视频| 国产一区二区三区在线臀色熟女| 级片在线观看| 精品人妻偷拍中文字幕| 国产精品精品国产色婷婷| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩黄片免| 国产av不卡久久| 69av精品久久久久久| 亚洲内射少妇av| 国产精品国产高清国产av| 国语自产精品视频在线第100页| 一个人看的www免费观看视频| 精品一区二区三区视频在线观看免费| av天堂在线播放| 91av网一区二区| 国产精品免费一区二区三区在线| 无遮挡黄片免费观看| 最新中文字幕久久久久| 日韩欧美精品v在线| 日韩欧美一区二区三区在线观看| 有码 亚洲区| 欧美性感艳星| 久久九九热精品免费| 亚洲国产色片| 波多野结衣高清无吗| av在线观看视频网站免费| 午夜福利18| 搡老岳熟女国产| 亚洲精品日韩av片在线观看| 久久久成人免费电影| 99久久精品热视频| 亚洲国产精品999在线| 波野结衣二区三区在线| 91字幕亚洲| 久久国产乱子伦精品免费另类| 一二三四社区在线视频社区8| 性插视频无遮挡在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 中文字幕人成人乱码亚洲影| 色播亚洲综合网| 日本在线视频免费播放| 亚洲,欧美,日韩| 亚洲av成人av| 一级av片app| 欧美+日韩+精品| 757午夜福利合集在线观看| 久久久久久久午夜电影| 蜜桃亚洲精品一区二区三区| 国产一级毛片七仙女欲春2| avwww免费| 狂野欧美白嫩少妇大欣赏| 久久久久久久久久成人| 国产欧美日韩一区二区三| 国产精品嫩草影院av在线观看 | 国产精品,欧美在线| 免费在线观看亚洲国产| 九九在线视频观看精品| 亚洲无线在线观看| 成人美女网站在线观看视频| 少妇人妻一区二区三区视频| 男人和女人高潮做爰伦理| 两性午夜刺激爽爽歪歪视频在线观看| 在线免费观看的www视频| 国产成人aa在线观看| 国产精品99久久久久久久久| 精品久久久久久久久av| 国产久久久一区二区三区| 999久久久精品免费观看国产| 变态另类丝袜制服| 国产aⅴ精品一区二区三区波| av在线蜜桃| 午夜福利免费观看在线| 国产视频内射| 成人av一区二区三区在线看| 乱人视频在线观看| 亚洲天堂国产精品一区在线| 男女下面进入的视频免费午夜| 久9热在线精品视频| 特级一级黄色大片| 丰满的人妻完整版| 日本在线视频免费播放| 午夜免费成人在线视频| 国产精品自产拍在线观看55亚洲| 亚洲人成网站高清观看| 97人妻精品一区二区三区麻豆| 两个人视频免费观看高清| 99热精品在线国产| 午夜福利欧美成人| 毛片女人毛片| 在线十欧美十亚洲十日本专区| 精品乱码久久久久久99久播| 久久久国产成人精品二区| 色哟哟哟哟哟哟| 日日摸夜夜添夜夜添av毛片 | 欧美日韩综合久久久久久 | 搡女人真爽免费视频火全软件 | 成人欧美大片| 国产三级在线视频| x7x7x7水蜜桃| 久久国产精品影院| 搡女人真爽免费视频火全软件 | 麻豆成人午夜福利视频| 精品不卡国产一区二区三区| 制服丝袜大香蕉在线| 欧美潮喷喷水| 国产高清激情床上av| 一进一出抽搐动态| 色综合婷婷激情| 日韩中字成人| 日韩有码中文字幕| 又黄又爽又免费观看的视频| 免费大片18禁| 夜夜躁狠狠躁天天躁| av中文乱码字幕在线| 老司机午夜十八禁免费视频| 永久网站在线| 国产精品久久久久久人妻精品电影| 免费看光身美女| 亚洲自偷自拍三级| 日本三级黄在线观看| 亚洲片人在线观看| 国产日本99.免费观看| 网址你懂的国产日韩在线| 啪啪无遮挡十八禁网站| 91狼人影院| 国产精品一及| 级片在线观看| 久久精品国产自在天天线| 禁无遮挡网站| 一区二区三区激情视频| 欧美在线一区亚洲| 精品久久久久久久人妻蜜臀av| 有码 亚洲区| 丰满人妻熟妇乱又伦精品不卡| 国产精品综合久久久久久久免费| 欧美色视频一区免费| 免费无遮挡裸体视频| 亚洲av成人av| 国产伦精品一区二区三区四那| 久久99热这里只有精品18| 青草久久国产| 成人高潮视频无遮挡免费网站| 国产一区二区在线av高清观看| 男女那种视频在线观看| 亚洲美女视频黄频| 性欧美人与动物交配| 精品人妻视频免费看| 麻豆国产av国片精品| 成人欧美大片| 淫秽高清视频在线观看| 99热这里只有是精品50| 两个人的视频大全免费| 很黄的视频免费| 亚洲电影在线观看av| 亚洲欧美日韩高清专用| 女同久久另类99精品国产91| 淫秽高清视频在线观看| 韩国av一区二区三区四区| 一二三四社区在线视频社区8| 亚洲精品乱码久久久v下载方式| 国产一区二区在线观看日韩| 精品人妻熟女av久视频| 国产三级在线视频| av女优亚洲男人天堂| 91久久精品电影网| 国产精品久久视频播放| 国内精品一区二区在线观看| 一个人免费在线观看的高清视频| 国产成+人综合+亚洲专区| 精品久久久久久久久av| 国产精品一区二区免费欧美| 18禁黄网站禁片免费观看直播| 精品久久久久久成人av| 亚洲精品日韩av片在线观看| 国产美女午夜福利| av在线观看视频网站免费| 精品一区二区三区视频在线| 免费av不卡在线播放| a在线观看视频网站| 国产午夜精品论理片| 欧美成人性av电影在线观看| 亚洲av中文字字幕乱码综合| av在线天堂中文字幕| 久久久精品大字幕| 亚洲国产精品999在线| 精品免费久久久久久久清纯| 日本三级黄在线观看| 少妇人妻一区二区三区视频| aaaaa片日本免费| 天堂影院成人在线观看| 精品一区二区免费观看| 久久久精品大字幕| 国产精品爽爽va在线观看网站| 欧美日韩国产亚洲二区| 成人午夜高清在线视频| 午夜两性在线视频| 俺也久久电影网| 亚洲成人免费电影在线观看| 日韩 亚洲 欧美在线| 怎么达到女性高潮| 久久热精品热| 欧美丝袜亚洲另类 | 国产日本99.免费观看| 国产成+人综合+亚洲专区| 午夜福利免费观看在线| 国产成人a区在线观看| 搡女人真爽免费视频火全软件 | 国产真实乱freesex| 老熟妇乱子伦视频在线观看| 亚洲欧美清纯卡通| 黄色视频,在线免费观看| 亚洲中文日韩欧美视频| 天堂动漫精品| 亚洲精品乱码久久久v下载方式| 亚洲av美国av| 级片在线观看| 国产主播在线观看一区二区| 怎么达到女性高潮| 国产高清视频在线观看网站| 欧美高清成人免费视频www| 亚洲内射少妇av| 99热精品在线国产| a级毛片a级免费在线| 99精品在免费线老司机午夜| 亚洲av成人av| 91狼人影院| 国产一区二区三区在线臀色熟女| 欧美潮喷喷水| 三级国产精品欧美在线观看| 国产精品精品国产色婷婷| 麻豆一二三区av精品| 欧美一区二区亚洲| 亚洲一区二区三区不卡视频| 亚洲自偷自拍三级| 久久久久免费精品人妻一区二区| 舔av片在线| 免费av观看视频| 搡女人真爽免费视频火全软件 | 日韩 亚洲 欧美在线| 亚洲中文日韩欧美视频| 久久久久久大精品| av黄色大香蕉| 人人妻,人人澡人人爽秒播| 亚洲国产精品成人综合色| 久久久久久久久中文| 舔av片在线| 国产亚洲精品av在线| 亚洲国产欧美人成| 亚洲中文日韩欧美视频| 淫秽高清视频在线观看| 在线观看舔阴道视频| 超碰av人人做人人爽久久| 韩国av一区二区三区四区| 俄罗斯特黄特色一大片| 91狼人影院| 日日干狠狠操夜夜爽| 欧美不卡视频在线免费观看| 91在线精品国自产拍蜜月| 老司机深夜福利视频在线观看| 舔av片在线| 久久久久久久亚洲中文字幕 | 日本撒尿小便嘘嘘汇集6| 亚洲精品一卡2卡三卡4卡5卡| 91久久精品国产一区二区成人| 国产熟女xx| 欧美性猛交╳xxx乱大交人| 日日夜夜操网爽| 深夜精品福利| 精品午夜福利视频在线观看一区| 国产真实伦视频高清在线观看 | 日本三级黄在线观看| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 免费黄网站久久成人精品 | 亚洲综合色惰| 狂野欧美白嫩少妇大欣赏| 日本熟妇午夜| 美女高潮喷水抽搐中文字幕| 日本黄色片子视频| 在线a可以看的网站| а√天堂www在线а√下载| 日韩有码中文字幕| 午夜激情福利司机影院| 久久久久精品国产欧美久久久| 精华霜和精华液先用哪个| 日韩人妻高清精品专区| 18禁在线播放成人免费| 国产色婷婷99| 亚洲乱码一区二区免费版| 久久久久久久精品吃奶| 日本一本二区三区精品| 人人妻,人人澡人人爽秒播| 国模一区二区三区四区视频| 99热精品在线国产| 男女下面进入的视频免费午夜| 久久久色成人| 怎么达到女性高潮| 亚洲精品日韩av片在线观看| 在现免费观看毛片| 1000部很黄的大片| 毛片一级片免费看久久久久 | 亚洲成av人片免费观看| 91麻豆av在线| 1000部很黄的大片| 夜夜爽天天搞| 日韩欧美在线乱码| 亚洲性夜色夜夜综合| 国内久久婷婷六月综合欲色啪| 国产精品影院久久| 性色avwww在线观看| 搡老熟女国产l中国老女人| 人妻制服诱惑在线中文字幕| av国产免费在线观看| 国产成人aa在线观看| 国产精品亚洲一级av第二区| 欧美成人性av电影在线观看| 夜夜夜夜夜久久久久| 亚洲av中文字字幕乱码综合| 成人三级黄色视频| 国产中年淑女户外野战色| 看十八女毛片水多多多| 国产精品女同一区二区软件 | 哪里可以看免费的av片| 91麻豆av在线| 在线观看av片永久免费下载| 亚洲精品久久国产高清桃花| netflix在线观看网站| 国产精品自产拍在线观看55亚洲| 久久久久久久亚洲中文字幕 | 亚洲性夜色夜夜综合| 精品人妻一区二区三区麻豆 | 88av欧美| 欧美乱妇无乱码| 亚洲色图av天堂| 精品午夜福利在线看| 人妻久久中文字幕网| 国产乱人视频| 久久精品国产99精品国产亚洲性色| 又粗又爽又猛毛片免费看| 久久久久久久精品吃奶| 亚洲午夜理论影院| 精品国产三级普通话版| 老司机午夜十八禁免费视频| 亚洲国产精品999在线| 久久久久性生活片| 精品一区二区三区人妻视频| 女人十人毛片免费观看3o分钟| 国产色爽女视频免费观看| 国产av不卡久久| 中文字幕人妻熟人妻熟丝袜美| 3wmmmm亚洲av在线观看| 桃色一区二区三区在线观看| 亚洲 欧美 日韩 在线 免费| 麻豆国产av国片精品| 九九在线视频观看精品| 日韩av在线大香蕉| 国产成+人综合+亚洲专区| 一个人观看的视频www高清免费观看| netflix在线观看网站| 久久久久久久午夜电影| 人妻丰满熟妇av一区二区三区| 日本 av在线| 男人和女人高潮做爰伦理| 日韩欧美国产一区二区入口| 男插女下体视频免费在线播放| 老熟妇乱子伦视频在线观看| 超碰av人人做人人爽久久| 国产伦精品一区二区三区四那| 欧美性猛交╳xxx乱大交人| av专区在线播放| 美女黄网站色视频| www.熟女人妻精品国产| 亚洲片人在线观看|