• <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
    狂野欧美白嫩少妇大欣赏| 动漫黄色视频在线观看| 国产一区在线观看成人免费| 好看av亚洲va欧美ⅴa在| 国产视频一区二区在线看| 好男人电影高清在线观看| 精品国产美女av久久久久小说| 日本a在线网址| 99久久精品国产亚洲精品| 成人欧美大片| 男人舔奶头视频| 免费电影在线观看免费观看| 欧美xxxx黑人xx丫x性爽| 精品日产1卡2卡| 日日夜夜操网爽| 免费av不卡在线播放| 巨乳人妻的诱惑在线观看| 成在线人永久免费视频| 国产69精品久久久久777片 | 日韩欧美精品v在线| 欧美一区二区国产精品久久精品| 在线免费观看不下载黄p国产 | 女生性感内裤真人,穿戴方法视频| 黄片大片在线免费观看| aaaaa片日本免费| 99久久综合精品五月天人人| 亚洲在线自拍视频| 亚洲色图av天堂| 巨乳人妻的诱惑在线观看| 男女之事视频高清在线观看| 亚洲avbb在线观看| 不卡av一区二区三区| 午夜免费成人在线视频| 免费大片18禁| 成人三级做爰电影| 一个人看的www免费观看视频| 国产又色又爽无遮挡免费看| 51午夜福利影视在线观看| 国产97色在线日韩免费| 天天一区二区日本电影三级| 男女下面进入的视频免费午夜| 麻豆国产97在线/欧美| 国产一区二区激情短视频| 国产午夜精品久久久久久| 丁香欧美五月| 亚洲国产精品999在线| 亚洲欧美日韩无卡精品| 在线观看一区二区三区| 丰满人妻一区二区三区视频av | 欧美+亚洲+日韩+国产| www日本黄色视频网| 国产精品免费一区二区三区在线| 亚洲精品国产精品久久久不卡| 亚洲av电影在线进入| 亚洲狠狠婷婷综合久久图片| 久久久久久久精品吃奶| 国产麻豆成人av免费视频| 免费观看精品视频网站| 亚洲专区中文字幕在线| 制服丝袜大香蕉在线| 99久久精品国产亚洲精品| 国产亚洲精品av在线| 亚洲精品一区av在线观看| 嫩草影院精品99| 国产日本99.免费观看| 亚洲中文日韩欧美视频| 午夜福利在线在线| 最近视频中文字幕2019在线8| 18禁黄网站禁片午夜丰满| 99久久精品一区二区三区| 我的老师免费观看完整版| 亚洲av成人av| 一进一出抽搐动态| 一个人看视频在线观看www免费 | 日本 欧美在线| svipshipincom国产片| h日本视频在线播放| 男女做爰动态图高潮gif福利片| 国产成人精品久久二区二区91| 91麻豆av在线| 全区人妻精品视频| 琪琪午夜伦伦电影理论片6080| 午夜激情福利司机影院| 日本熟妇午夜| 久久精品国产亚洲av香蕉五月| 啪啪无遮挡十八禁网站| 国产在线精品亚洲第一网站| 亚洲国产高清在线一区二区三| 欧美黑人欧美精品刺激| 亚洲人成电影免费在线| 黑人操中国人逼视频| 最好的美女福利视频网| 夜夜爽天天搞| 级片在线观看| 此物有八面人人有两片| 亚洲七黄色美女视频| 欧美极品一区二区三区四区| 全区人妻精品视频| 噜噜噜噜噜久久久久久91| e午夜精品久久久久久久| 久久久久免费精品人妻一区二区| 久久人人精品亚洲av| 免费观看人在逋| 中出人妻视频一区二区| 九九在线视频观看精品| 国产亚洲欧美在线一区二区| 国产v大片淫在线免费观看| 成人国产综合亚洲| 日韩欧美三级三区| 精品欧美国产一区二区三| 亚洲成a人片在线一区二区| 亚洲国产精品久久男人天堂| 手机成人av网站| 国内揄拍国产精品人妻在线| 国产淫片久久久久久久久 | 一级作爱视频免费观看| 国产精品亚洲美女久久久| 亚洲成人中文字幕在线播放| www.www免费av| 亚洲在线自拍视频| 亚洲国产精品sss在线观看| 一个人观看的视频www高清免费观看 | 99精品在免费线老司机午夜| 亚洲国产精品sss在线观看| 久久亚洲真实| 黄色女人牲交| 亚洲欧美日韩高清专用| 中文资源天堂在线| 国产精品,欧美在线| 真人一进一出gif抽搐免费| 亚洲av中文字字幕乱码综合| 欧美av亚洲av综合av国产av| 久久精品91无色码中文字幕| 亚洲av中文字字幕乱码综合| 国产高清激情床上av| 桃红色精品国产亚洲av| 美女被艹到高潮喷水动态| 最好的美女福利视频网| 成人无遮挡网站| 午夜精品一区二区三区免费看| 日韩欧美 国产精品| 国产一区二区三区在线臀色熟女| 亚洲第一电影网av| 三级国产精品欧美在线观看 | 天天一区二区日本电影三级| 久久天堂一区二区三区四区| 国产麻豆成人av免费视频| 欧美绝顶高潮抽搐喷水| 久久久久久久久久黄片| 国产精品亚洲av一区麻豆| h日本视频在线播放| 深夜精品福利| 又爽又黄无遮挡网站| 国产精品99久久99久久久不卡| 五月玫瑰六月丁香| 久久久久国产一级毛片高清牌| 变态另类成人亚洲欧美熟女| 久久久久久大精品| 国产人伦9x9x在线观看| 国产三级在线视频| 午夜精品在线福利| 好看av亚洲va欧美ⅴa在| 国产高清激情床上av| 精品欧美国产一区二区三| 91九色精品人成在线观看| av欧美777| 青草久久国产| 欧美高清成人免费视频www| 成人三级黄色视频| 俺也久久电影网| 亚洲av中文字字幕乱码综合| 午夜福利欧美成人| 热99re8久久精品国产| 在线十欧美十亚洲十日本专区| 丰满的人妻完整版| 首页视频小说图片口味搜索| 国内精品美女久久久久久| 精品国产乱码久久久久久男人| 色尼玛亚洲综合影院| 12—13女人毛片做爰片一| 少妇裸体淫交视频免费看高清| 成人18禁在线播放| 亚洲真实伦在线观看| 欧美在线一区亚洲| 成人一区二区视频在线观看| 亚洲片人在线观看| 很黄的视频免费| 亚洲av五月六月丁香网| 综合色av麻豆| 天堂影院成人在线观看| 亚洲 国产 在线| 91在线精品国自产拍蜜月 | 亚洲熟女毛片儿| 欧美日韩中文字幕国产精品一区二区三区| www.熟女人妻精品国产| 国产三级在线视频| 免费观看的影片在线观看| 黑人操中国人逼视频| 久久99热这里只有精品18| bbb黄色大片| 色老头精品视频在线观看| 亚洲精品美女久久久久99蜜臀| 国产成人aa在线观看| 亚洲av中文字字幕乱码综合| 成人av一区二区三区在线看| 欧美+亚洲+日韩+国产| 日本免费一区二区三区高清不卡| 免费观看人在逋| 观看免费一级毛片| 美女免费视频网站| 在线观看免费午夜福利视频| 老司机在亚洲福利影院| 国内毛片毛片毛片毛片毛片| 最新美女视频免费是黄的| 国产精品影院久久| 久久午夜亚洲精品久久| 老司机深夜福利视频在线观看| 午夜免费观看网址| 人妻久久中文字幕网| 三级国产精品欧美在线观看 | 小说图片视频综合网站| 综合色av麻豆| 精品国产乱子伦一区二区三区| 村上凉子中文字幕在线| 亚洲成人久久爱视频| h日本视频在线播放| 视频区欧美日本亚洲| 国产精品永久免费网站| 精品久久久久久,| 美女cb高潮喷水在线观看 | 免费看十八禁软件| 国产免费男女视频| 哪里可以看免费的av片| 精品日产1卡2卡| 国产精品综合久久久久久久免费| 久久久久久久精品吃奶| 国产伦在线观看视频一区| 美女黄网站色视频| 老司机午夜福利在线观看视频| 国产99白浆流出| 亚洲七黄色美女视频| 国产精品电影一区二区三区| 国产私拍福利视频在线观看| 少妇的丰满在线观看| 国内精品美女久久久久久| 欧美成人性av电影在线观看| 欧美日本视频| 久久久久久久午夜电影| 啦啦啦观看免费观看视频高清| 中文字幕人成人乱码亚洲影| 色视频www国产| 色哟哟哟哟哟哟| 欧美zozozo另类| 波多野结衣巨乳人妻| 香蕉av资源在线| 女人被狂操c到高潮| 色哟哟哟哟哟哟| 亚洲国产精品成人综合色| 丁香六月欧美| 五月玫瑰六月丁香| 欧洲精品卡2卡3卡4卡5卡区| 天天躁狠狠躁夜夜躁狠狠躁| 99久久国产精品久久久| 麻豆久久精品国产亚洲av| 99国产精品一区二区蜜桃av| 成人永久免费在线观看视频| 哪里可以看免费的av片| 男女那种视频在线观看| 99精品久久久久人妻精品| 麻豆一二三区av精品| 国产精品久久久久久精品电影| 97超级碰碰碰精品色视频在线观看| 午夜两性在线视频| 欧美xxxx黑人xx丫x性爽| 久久精品亚洲精品国产色婷小说| 亚洲av电影在线进入| 国产精品一区二区免费欧美| 亚洲激情在线av| 久久精品国产99精品国产亚洲性色| 12—13女人毛片做爰片一| 国产精品影院久久| 午夜精品久久久久久毛片777| 亚洲av五月六月丁香网| 日日干狠狠操夜夜爽| 中文字幕高清在线视频| 欧美午夜高清在线| 欧美中文综合在线视频| 亚洲国产精品999在线| 好男人在线观看高清免费视频| 欧美极品一区二区三区四区| 男插女下体视频免费在线播放| 性欧美人与动物交配| 亚洲av中文字字幕乱码综合| 97碰自拍视频| 日韩av在线大香蕉| 国产成人av激情在线播放| 巨乳人妻的诱惑在线观看| 1024手机看黄色片| 熟妇人妻久久中文字幕3abv| 国产精品av视频在线免费观看| 性色avwww在线观看| 国产欧美日韩精品一区二区| 国产真人三级小视频在线观看| 禁无遮挡网站| 人人妻人人看人人澡| 色在线成人网| 日本免费a在线| 亚洲国产精品久久男人天堂| 国产精品久久久久久久电影 | 男人舔女人下体高潮全视频| 99riav亚洲国产免费| 啦啦啦观看免费观看视频高清| 亚洲专区中文字幕在线| 天天添夜夜摸| 狂野欧美激情性xxxx| 国产乱人伦免费视频| 亚洲欧美日韩无卡精品| 床上黄色一级片| 热99re8久久精品国产| 久久久久精品国产欧美久久久| 欧美日韩综合久久久久久 | 丰满人妻一区二区三区视频av | 国产亚洲精品久久久com| 国产av麻豆久久久久久久| 最好的美女福利视频网| 免费看十八禁软件| 亚洲国产中文字幕在线视频| 国产精品一区二区精品视频观看| 婷婷六月久久综合丁香| 在线观看免费午夜福利视频| 少妇熟女aⅴ在线视频| 日本 av在线| 国产欧美日韩一区二区精品| 91九色精品人成在线观看| 男女做爰动态图高潮gif福利片| 午夜精品在线福利| 亚洲无线观看免费| 久久精品国产亚洲av香蕉五月| 亚洲国产色片| 又黄又粗又硬又大视频| 欧美午夜高清在线| 18禁黄网站禁片午夜丰满| 老熟妇乱子伦视频在线观看| 十八禁网站免费在线| 欧美成人免费av一区二区三区| 亚洲精品在线美女| 麻豆av在线久日| 国语自产精品视频在线第100页| 日日干狠狠操夜夜爽| 成人永久免费在线观看视频| 精品免费久久久久久久清纯| 757午夜福利合集在线观看| 一个人免费在线观看电影 | a级毛片在线看网站| 亚洲最大成人中文| 国产精品久久久久久人妻精品电影| 亚洲最大成人中文| 亚洲黑人精品在线| 三级男女做爰猛烈吃奶摸视频| 欧美在线一区亚洲| 亚洲专区国产一区二区| a在线观看视频网站| 曰老女人黄片| 欧美xxxx黑人xx丫x性爽| 亚洲精品一卡2卡三卡4卡5卡| 日韩av在线大香蕉| 国产精品一区二区三区四区免费观看 | 91av网站免费观看| 熟妇人妻久久中文字幕3abv| 我要搜黄色片| 又黄又爽又免费观看的视频| 久久亚洲真实| 熟女电影av网| 好男人在线观看高清免费视频| 欧美日韩黄片免| 久久午夜亚洲精品久久| 国产爱豆传媒在线观看| 很黄的视频免费| 日本五十路高清| 日本撒尿小便嘘嘘汇集6| 欧美一区二区精品小视频在线| 真人一进一出gif抽搐免费| 国产一区二区三区视频了| 午夜精品在线福利| 国产爱豆传媒在线观看| 非洲黑人性xxxx精品又粗又长| 此物有八面人人有两片| 亚洲一区高清亚洲精品| 美女午夜性视频免费| 老司机福利观看| 最新中文字幕久久久久 | 婷婷精品国产亚洲av在线| 精品一区二区三区四区五区乱码| 午夜福利在线在线| 成年女人毛片免费观看观看9| 97人妻精品一区二区三区麻豆| 亚洲欧美精品综合久久99| www日本黄色视频网| 校园春色视频在线观看| 国产欧美日韩一区二区三| 国产 一区 欧美 日韩| 美女高潮喷水抽搐中文字幕| 美女被艹到高潮喷水动态| 欧美日韩中文字幕国产精品一区二区三区| 欧美日本视频| 欧美zozozo另类| 国产精品99久久99久久久不卡| 午夜福利在线观看吧| 又紧又爽又黄一区二区| 国产成人影院久久av| 九九在线视频观看精品| 午夜精品久久久久久毛片777| ponron亚洲| 很黄的视频免费| 亚洲成人精品中文字幕电影| 成人亚洲精品av一区二区| 男女下面进入的视频免费午夜| 日本三级黄在线观看| 精品一区二区三区四区五区乱码| 日本黄色片子视频| 一区二区三区高清视频在线| 非洲黑人性xxxx精品又粗又长| 2021天堂中文幕一二区在线观| 亚洲成a人片在线一区二区| 亚洲国产精品999在线| 一卡2卡三卡四卡精品乱码亚洲| 男女下面进入的视频免费午夜| 国产精品永久免费网站| 天堂动漫精品| 欧美高清成人免费视频www| 制服丝袜大香蕉在线| 曰老女人黄片| 毛片女人毛片| 999久久久国产精品视频| 日本撒尿小便嘘嘘汇集6| 少妇人妻一区二区三区视频| 中出人妻视频一区二区| 国产欧美日韩一区二区三| 国产成人av激情在线播放| 亚洲中文av在线| av女优亚洲男人天堂 | 国产欧美日韩一区二区精品| 婷婷精品国产亚洲av在线| 亚洲精品456在线播放app | 一本精品99久久精品77| 麻豆国产97在线/欧美| 亚洲国产日韩欧美精品在线观看 | www国产在线视频色| 亚洲片人在线观看| 亚洲在线自拍视频| 国产一区二区三区在线臀色熟女| 久久久久国产精品人妻aⅴ院| 最近最新免费中文字幕在线| 黄色 视频免费看| 在线免费观看不下载黄p国产 | 美女扒开内裤让男人捅视频| 国产精品,欧美在线| 国产精品电影一区二区三区| 免费观看人在逋| 婷婷精品国产亚洲av| 成年人黄色毛片网站| 日本免费一区二区三区高清不卡| www国产在线视频色| 久久婷婷人人爽人人干人人爱| 最近最新中文字幕大全电影3| 麻豆国产97在线/欧美| 亚洲国产色片| 麻豆成人av在线观看| 1024手机看黄色片| 亚洲美女黄片视频| 亚洲第一电影网av| 麻豆国产97在线/欧美| 国产伦在线观看视频一区| 母亲3免费完整高清在线观看| 精品不卡国产一区二区三区| 一进一出抽搐动态| 亚洲国产看品久久| 日本一本二区三区精品| 午夜福利免费观看在线| 国产精品一区二区三区四区免费观看 | 一个人免费在线观看电影 | 精品熟女少妇八av免费久了| 欧美av亚洲av综合av国产av| 制服人妻中文乱码| 精品久久久久久久毛片微露脸| 亚洲国产精品合色在线| 欧美成狂野欧美在线观看| 又粗又爽又猛毛片免费看| 国产单亲对白刺激| 亚洲,欧美精品.| 色综合婷婷激情| 欧美在线一区亚洲| 12—13女人毛片做爰片一| 欧美黄色淫秽网站| 1024香蕉在线观看| 国产伦精品一区二区三区四那| 91av网站免费观看| 中文亚洲av片在线观看爽| 热99re8久久精品国产| 国产不卡一卡二| 91老司机精品| 少妇裸体淫交视频免费看高清| 一级黄色大片毛片| 美女 人体艺术 gogo| 欧美在线黄色| 久久久久久久久免费视频了| 欧美性猛交╳xxx乱大交人| 成年版毛片免费区| 精品熟女少妇八av免费久了| 美女午夜性视频免费| 淫妇啪啪啪对白视频| or卡值多少钱| 日本成人三级电影网站| 亚洲人成伊人成综合网2020| 日韩欧美一区二区三区在线观看| 成人欧美大片| 欧美在线一区亚洲| 国产精品综合久久久久久久免费| 成人亚洲精品av一区二区| 天堂动漫精品| 免费看光身美女| 欧美成狂野欧美在线观看| 无遮挡黄片免费观看| 精品久久久久久久久久久久久| 国产精品日韩av在线免费观看| 757午夜福利合集在线观看| 精品欧美国产一区二区三| 久久久国产欧美日韩av| 日本一本二区三区精品| 制服丝袜大香蕉在线| 成人特级黄色片久久久久久久| 精品一区二区三区四区五区乱码| 久久99热这里只有精品18| 小蜜桃在线观看免费完整版高清| 久久久久久久精品吃奶| 午夜福利高清视频| 欧美一级毛片孕妇| 十八禁网站免费在线| 亚洲va日本ⅴa欧美va伊人久久| 黄片大片在线免费观看| 久久久国产欧美日韩av| 亚洲av成人精品一区久久| 少妇裸体淫交视频免费看高清| 日韩中文字幕欧美一区二区| 欧美3d第一页| 中文字幕人妻丝袜一区二区| 少妇熟女aⅴ在线视频| 麻豆国产av国片精品| 99re在线观看精品视频| 窝窝影院91人妻| 69av精品久久久久久| 国内少妇人妻偷人精品xxx网站 | 91av网一区二区| 国产在线精品亚洲第一网站| 一夜夜www| 麻豆久久精品国产亚洲av| 中出人妻视频一区二区| 亚洲熟妇熟女久久| 欧美日韩瑟瑟在线播放| 99精品久久久久人妻精品| 国内揄拍国产精品人妻在线| 中文字幕高清在线视频| 99riav亚洲国产免费| 亚洲乱码一区二区免费版| 精品久久久久久久久久免费视频| 亚洲欧美一区二区三区黑人| 国产精品综合久久久久久久免费| 法律面前人人平等表现在哪些方面| 欧美性猛交黑人性爽| 一二三四在线观看免费中文在| 在线看三级毛片| 日韩欧美三级三区| 日韩国内少妇激情av| 午夜精品久久久久久毛片777| 午夜a级毛片| 亚洲国产欧美网| 午夜福利视频1000在线观看| 97人妻精品一区二区三区麻豆| 亚洲 欧美一区二区三区| 免费人成视频x8x8入口观看| 日韩人妻高清精品专区| 国产人伦9x9x在线观看| 免费观看人在逋| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产精品久久男人天堂| 欧美性猛交╳xxx乱大交人| 亚洲精品在线美女| 欧美极品一区二区三区四区| 综合色av麻豆| 一区二区三区激情视频| 老司机福利观看| 综合色av麻豆| 亚洲精品在线美女| 欧美极品一区二区三区四区| a级毛片a级免费在线| 全区人妻精品视频| 久久精品国产99精品国产亚洲性色| 哪里可以看免费的av片| 黄色成人免费大全| 国产69精品久久久久777片 | 国产免费男女视频| 欧美日韩福利视频一区二区| 色综合亚洲欧美另类图片| 一级毛片精品| 日韩欧美免费精品| 亚洲av五月六月丁香网| 18美女黄网站色大片免费观看| 国产精品久久久久久久电影 | 午夜免费激情av| 国产男靠女视频免费网站| 亚洲第一电影网av| 国产精品日韩av在线免费观看| 国产精品一区二区三区四区免费观看 | 国产亚洲av高清不卡| 18禁黄网站禁片午夜丰满|