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

    裂谷盆地構造-熱演化模擬中幾個問題的討論

    2015-12-12 07:11:52劉瓊穎何麗娟
    地球物理學報 2015年2期
    關鍵詞:側翼巖石圈熱導率

    劉瓊穎,何麗娟

    1中國科學院地質(zhì)與地球物理研究所巖石圈演化國家重點實驗室,北京 100029

    2中國科學院大學,北京 100049

    1 引言

    沉積盆地構造-熱演化模擬是盆地模擬的主要內(nèi)容之一.構造沉降史與熱流演化史構成其兩大核心研究內(nèi)容.拉張模型在描述拉張盆地沉降和熱流演化方面取得了極大的成功,實現(xiàn)了構造和熱的完美結合(何麗娟和汪集暘,2007).

    McKenzie(1978)提出了一維瞬時均勻拉張模型.該模型認為,巖石圈的拉張在瞬時完成,由于巖石圈減薄,引起熱的軟流圈物質(zhì)被動上涌,導致巖石圈被加熱,在Airy均衡作用下形成初始沉降.拉張結束后,巖石圈進入熱冷卻階段,地表在初始沉降的基礎上繼續(xù)發(fā)生熱沉降.該模型是運動學模型的代表模型,能很好地預測盆地主要觀測數(shù)據(jù),為以后研究拉張盆地構造熱演化奠定了基礎.考慮到巖石圈拉張不可能瞬時完成,學者們進一步提出有限勻速拉張模型(Jarvis and McKenzie,1980).何麗娟(1999)提出了勻減速拉張速率模式,更接近實際情況.一些學者在均勻拉張模型的基礎上提出了縱向非均勻拉張模型,如雙層非均勻拉張模型(Royden and Keen,1980;Hellinger and Sclater,1983;Dehler et al.,1997)、與深度相關的拉張模型(Rowley and Sahagian,1986;Davis and Kusznir,2004;Kusznir and Karner,2007;Rosenbaum et al.,2008).上述拉張模型多假設巖石圈拉張變形為純剪變形,但實際上變形并非總是對稱的.Wernicke(1985)提出了簡單剪切拉張模型,但迄今研究者并未觀測到切穿整個巖石圈的斷裂.為此,學者們用耦合的簡單剪切/純剪切模型來描述不對稱拉張的情形,具有一定的實用性(Kusznir and Ziegler,1992;Clift et al.,2002;Baur et al.,2010;Chen et al.,2013).一些盆地在地質(zhì)歷史時期曾經(jīng)歷過多期拉張,且具有繼承性,多期拉張模型由此被提出(Wang et al.,1996;Dehler et al.,1997;He et al.,2002;Lin et al.,2002).White(1994)以及 White和Bellingham(2002)提出了利用沉降數(shù)據(jù)反演巖石圈拉張期間一維和二維應變速率的方法,該方法不需要知道巖石圈拉張持續(xù)時間、拉張期次等先驗信息,在一些地區(qū)得到了較好的應用(Xie et al.,2006;Shillington et al.,2008;Song et al.2010;Chen,2014).考慮到盆地拉張過程中可能的構造反轉,研究者們進一步建立了適用于多期伸展和擠壓交替的復雜過程的數(shù)值計算模型(Negredo et al.,1999;He and Wang,2004;趙長煜等,2010).以上模型均為運動學模型,巖石圈變形通過定義的速度場來實現(xiàn),它們在預測盆地熱流演化方面具有優(yōu)勢,但巖石圈流變性隱含在其均衡補償模式中,無法反映巖石圈流變性及其變化對盆地構造熱演化的影響,在預測盆地幾何形態(tài)方面存在一定的缺陷(何麗娟,2000).

    而在拉張模型的另一類模型,動力學模型中,巖石圈流變性則是通過本構方程明確體現(xiàn),能夠?qū)崿F(xiàn)熱學與流變學的耦合,深入理解盆地動力學的基本過程(Fernàndez and Ranalli,1997).根據(jù)其本構方程,動力學模型可分為平面應變模型(Braun and Beaumont,1989;Bassi,1991;Chéry et al.,1992;Huismans and Beaumont,2008;Van Avendonk et al.,2009;Huismans and Beaumont,2011)和平面應力模型(Sonder and England,1989;Bassi and Sabadini,1994;Polyansky,2002).由于目前對盆地成因的動力學根本機制還不甚了解,以及模型初始擾動的給定、巖石圈流變學參數(shù)的不確定性等都在一定程度上阻礙了動力學模型的發(fā)展.

    現(xiàn)有的拉張模型多采用固定的初始巖石圈、地殼厚度,軟流圈多被賦予和巖石圈地幔相同的熱參數(shù),模型上邊界的給定存在爭議,且模型中較少地考慮沉積蓋層的熱效應.本文擬重點探討拉張模型中初始地殼、巖石圈厚度、軟流圈對流、模型上邊界對構造熱演化的影響,以及沉積蓋層存在的情況下拉張盆地側翼抬升的演化.文中的模擬均使用運動學模型,假定巖石圈為純剪變形,利用有限元方法,在拉格朗日坐標系下,通過求解二維瞬態(tài)熱傳導方程來模擬盆地構造位移發(fā)生變化的同時,溫度場和熱流在時間、空間上的演化歷史.

    2 初始巖石圈、地殼厚度對構造熱演化的影響

    多數(shù)拉張模型都假定拉張開始前模型的初始巖石圈厚度為125km (McKenzie,1978;Jarvis and McKenzie,1980;He et al.,2002;Lin et al.,2002;He and Wang,2004;陳林等,2008;Song et al.,2010;Chen,2014),這也是大洋巖石圈厚度的漸進值(Parsons and Sclater,1977),該值在拉張開始前巖石圈未經(jīng)歷拉張或破壞的情況下往往是適用的.然而,在一些地區(qū)拉張開始前巖石圈厚度已經(jīng)發(fā)生了較大的變化,比如中生代的華北克拉通.這時,需要根據(jù)實際需要選取合適的初始巖石圈厚度.

    對于初始地殼厚度的選取,往往依賴于現(xiàn)今的Moho面深度.Dehler等(1997)在計算加拿大西緣夏洛特皇后盆地的構造熱演化時使用的初始地殼厚度為34km,He和Wang(2004)在模擬濟陽坳陷的熱史時給定的初始地殼厚度是36km,Chen(2014)模擬南海白云凹陷時采用的初始地殼厚度為30km,這些模型預測的剖面現(xiàn)今Moho面深度都和實際觀測值接近.

    圖1 初始地殼、巖石圈厚度對構造-熱演化的影響Fig.1 Effects of initial crustal and lithospheric thicknesses on tectono-thermal modeling

    表1 模型中使用的參數(shù)Table 1 Parameters and values used in modeling

    為了了解初始地殼、巖石圈厚度對構造熱演化的影響,本文建立了兩組模型(圖1).圖1a和b固定初始地殼厚度為35km,模型拉張10Ma,拉張系數(shù)為1.5,熱沉降20Ma,計算不同初始巖石圈厚度下的構造沉降史及熱流演化史.圖1c和d是固定初始巖石圈厚度為110km,改變初始地殼厚度,計算構造沉降量及熱流值.模擬中使用的參數(shù)見表1.結果顯示,當初始地殼厚度恒定,初始巖石圈厚度從70km增加到120km時,計算得到的構造沉降量會隨之減小,最大差異達660m以上;熱流也隨著初始巖石圈厚度的增加而降低,最大差異近25mW·m-2.當初始巖石圈厚度恒定,初始地殼厚度從31km增加到41km時,構造沉降量隨之增大,最大差異超過900m;熱流也與初始地殼厚度成正比,但最大差異不到6mW·m-2.由此看出,初始巖石圈、地殼厚度的選擇,對構造熱演化的影響非常大,尤其是對構造沉降量的影響尤為顯著.

    實際計算中,當拉張開始前的初始地殼、巖石圈厚度未知時,可給出多組初始厚度值.每給出一組初始地殼、巖石圈厚度值,通過模擬就可以計算得到一組現(xiàn)今熱流、Moho面深度及熱巖石圈厚度值,將計算值與觀測值相對比,通過不斷調(diào)整初始地殼、巖石圈厚度值,即可反推得到模擬剖面的最佳初始地殼、巖石圈厚度值.

    3 軟流圈對流對構造熱演化的影響

    大部分拉張模型是將軟流圈和巖石圈地幔賦予相同的熱導率,在計算中僅考慮熱傳導作用(McKenzie,1978;Keen and Dehler,1993;He et al.,2002;Lin et al.,2009;Chen,2014).而 Lachenbruch(1978)首次指出了拉張背景下對流在熱傳遞中的重要性.Steckler(1985)、Buck(1986)、van Wijk等(2008)都進一步指出被動裂谷作用會引起軟流圈的小范圍對流,其作為額外的熱源會對沉降、抬升史及溫度場產(chǎn)生重要影響.

    為了討論軟流圈對流對構造熱演化的影響,本文采用等效熱傳導率綜合考慮軟流圈熱傳導和熱對流因素(范桃園和安美建,2009).事實上,通過在熱傳導方程中提高熱導率以等效模擬對流的熱效應的方法已廣泛被學者們采用(Morgan and Chen,1993;Chen,2001;Spinelli and Harris,2011;楊少華和石耀霖,2013;Schmeling and Marquart,2014).在純對流的情況下,相對于正常的熱導率,提高的倍數(shù)往往為表征對流傳熱效率與單純的熱傳導效率之比的努塞爾特數(shù),Nu(Morgan and Chen,1993;Chen,2001;Cochran and Buck,2001;楊少華和石耀霖,2013).在研究Nu隨系統(tǒng)總體熱狀態(tài)變化特征中,實驗和數(shù)值模擬都給出如下規(guī)律(Schubert et al.,2001;Turcotte and Schubert,2002):

    式中,c,β為常數(shù),β值一般在0.3左右,Ra為瑞利數(shù).研究表明,當Ra接近5×105~106時,c=0.27,β=0.3185(Schubert et al.,2001).而在上地幔對流的情況下,Ra通常為106左右(Schubert et al.,2001;Jaupart and Mareschal,2011).利用上式可以求得此時Nu接近22,也就是說,當認為軟流圈地幔以純對流的形式熱傳遞的話,其等效熱導率應近似為22倍的正常地幔熱導率,約為70W·m-1K-1.本文計算了軟流圈等效熱導率從3.2W·m-1K-1(純熱傳導,和巖石圈地幔熱導率相同)到70W·m-1K-1(等效于純熱對流)變化時的構造沉降量、熱流及溫度場.結果顯示,軟流圈等效熱導率對熱沉降期的構造沉降量及熱流影響很大,而在拉張期其影響相對微弱.考慮了軟流圈的對流效應之后,相比于單純的熱傳導模型,構造沉降量變小,而熱流則增大.對流越強,這種影響越明顯(圖2a和b).在溫度場方面,整個過程中考慮軟流圈對流的模型其下地殼至模型底部的溫度明顯比僅考慮熱傳導的模型溫度高.在熱沉降階段,軟流圈純對流模型中溫度場變化甚微,而在軟流圈純熱傳導模型中則能明顯的觀察到溫度場的冷卻(圖2c).

    圖2 軟流圈等效熱導率對構造熱演化的影響(拉張期為10Ma,拉張系數(shù)為1.5,熱沉降期為20Ma)(a)構造沉降史;(b)熱流隨時間的變化;(c)溫度-深度曲線.Fig.2 Effects of asthenosphere equivalent thermal conductivity on tectono-thermal modeling(The extension lasts 10Ma with a stretching factor of 1.5,and then enters thermal subsidence until 30Ma)(a)Tectonic subsidence history;(b)Heat flow history;(c)Temperature-depth curve.

    軟流圈純傳導模型的熱流隨時間的降低速度要比實際情況快得多,而軟流圈純對流模型計算的沉降量要比實際小得多,二者都與實際偏差較大.事實上,純傳導的固體巖石圈與純對流的流體軟流圈之間存在一過渡層,即流變邊界層(Sleep,2006;Artemieva,2011;圖3),其間傳導與對流共同作用來傳遞熱量.在拉張模型中,同時考慮軟流圈的熱傳導與對流更接近真實的情況,純傳導與純對流應作為其下限與上限.在模擬中,得知地幔對流的努塞爾特數(shù)Nu或瑞利數(shù)Ra時,可相應的判斷軟流圈等效熱導率的上限.流變邊界層的厚度對巖石圈本身的結構特征并不敏感,而主要受軟流圈黏性系數(shù)η控制,其厚度與lg(η)成正比(He,2014).考慮到軟流圈對流的瑞利數(shù)Ra與黏性系數(shù)η成反比,而由式(1)可知努塞爾特數(shù)Nu正比于Raβ,由此可推斷地幔對流的努塞爾特數(shù)Nu越大,則流變邊界層的厚度越小.這也意味著在拉張模型中,流變邊界層厚度小的地方,流變邊界層中等效熱傳導率的上限變大,而其下純對流軟流圈的等效熱導率為Nu倍的正常地幔熱導率,也相應大于流變邊界層厚的地方的軟流圈等效熱導率.其對構造熱演化的影響由上述模擬結果可以得知,流變邊界層厚度小的地方,構造沉降量速率變慢,總構造沉降量變小,而熱沉降過程中熱衰減變緩,熱流降低的總幅度變小.

    4 模型上邊界對構造熱演化的影響

    圖3 溫度(T)、地溫梯度(G)隨深度的變化(純傳導的巖石圈與純對流的軟流圈之間為流變邊界層RBL,改自He,2014)Fig.3 Variations in temperature (T)and temperature gradient(G)with depth(The transition layer between the conductive lithosphere and the convective asthenosphere is the rheological boundary layer RBL,modified from He,2014).

    現(xiàn)有的拉張模型對其上邊界采用兩種處理方法.一種是整個過程一直固定上邊界,其并不隨計算出的沉降量而作調(diào)整(McKenzie,1978;Rowley and Sahagian,1986;Keen and Dehler,1993),是對模型的一種簡化;另一種則是根據(jù)均衡補償?shù)玫降臉嬙斐两盗繉δP蜕线吔缭诖瓜蛏喜粩噙M行調(diào)整,更真實的還原盆地的演化過程(Wang et al.,1996;He and Wang,2004;Wangen et al.,2008).為了探討固定與移動邊界對構造熱演化的影響,我們建立了以下兩組模型.一組模型的拉張時間為10Ma,之后裂后熱松弛持續(xù)100Ma(圖4a和b).另一組模型的拉張時間為40Ma,裂后熱沉降時間同為100Ma(圖4c和d).兩組模型都在拉張系數(shù)分別為2、4和6時,對比固定與移動上邊界模型的計算結果.結果顯示,當拉張時間為10Ma時,固定與移動上邊界對拉張期的熱流及沉降量影響甚微;但當拉張時間增加為40Ma時,模型上邊界的處理方式對拉張期的熱流及初始沉降產(chǎn)生了明顯的影響,尤其是在拉張系數(shù)較大的情況下,但這種影響小于熱沉降期兩種模型的構造-熱演化差異.當拉張系數(shù)為6時,拉張結束時兩種模型的熱流差異為0.3mW·m-2,初始沉降量差異為30m(圖4c和d).在熱沉降階段,無論拉張時間的長短,固定與移動上邊界模型的熱流及構造沉降的差異隨熱沉降的持續(xù)均逐漸變大,移動邊界模型的熱流大于固定邊界的,而構造沉降量小于固定邊界模型求得的值.但上述兩組拉張時間不同、熱沉降時間相同的模型在熱沉降結束時,由于上邊界處理方式不同造成的熱流、總沉降量及溫度場的差異相差并不大,以下結果以拉張時間為40Ma的一組模擬為例.拉張系數(shù)為2時,固定與移動上邊界模型的熱流及構造沉降的最大差值分別為0.7mW·m-2及85m;當拉張系數(shù)為4時,二者的熱流及構造沉降最大差異增大為1.7mW·m-2和162m;當拉張系數(shù)進一步增大到6時,兩種模型的熱流最大差異變?yōu)?.4mW·m-2,構造沉降量的最大差異為198m(圖4c和d).在溫度場方面,從圖4e—g可以看出,固定與移動上邊界模型的溫度結構有著明顯的差別,尤其是在巖石圈部分,相同深度上,固定邊界模型的溫度比移動邊界模型的溫度高幾十到上百攝氏度.當拉張系數(shù)分別為2、4和6時,兩種模型在同一深度上溫度的最大差異分別為95、161及185℃.

    總體上看,模型上邊界對構造熱演化的影響隨著拉張系數(shù)的增大而增強.相比于模型上邊界的給定對熱流及構造沉降量的影響,其對模型溫度場有著更為顯著的影響.從上述討論也可以看到,當計算主要針對熱流及構造沉降量時,在拉張系數(shù)小于2或是拉張時間小于10Ma且只關注于拉張期的熱流及構造沉降變化時,兩種模型的計算結果差異并不大,使用固定上邊界模型的計算結果也是可接受的.但在拉張強度較大的地區(qū),如中國南海的鶯歌海盆地等(拉張系數(shù)大于4,He et al.,2002),使用簡化的固定上邊界模型將對計算帶來較大的誤差.

    5 拉張盆地的側翼隆起

    裂谷的側翼隆起常作為拉張盆地及張裂大陸邊緣的邊界,其是拉張盆地模擬的重要組成.學者們提出了多種機制來解釋拉張盆地的側翼隆起:(1)由拉張巖石圈向未拉張區(qū)域的側向熱傳導(Alvarez et al.,1984;Buck et al.,1988;Leroy et al.,2008);(2)由拉張引起的小規(guī)模地幔對流造成的裂谷之下向側翼的熱傳遞(Keen,1985;Buck,1986);(3)裂谷側翼之下相比于地殼減薄更強的巖石圈地幔減?。℉ellinger and Sclater,1983;White and McKenzie,1988;Dehler et al.,1997);(4)拉張過程中通過向上撓曲對裂谷側翼的力學支持(Braun and Beaumont,1989;Weissel and Karner,1989;Chéry et al.,1992;Kooi et al.,1992;van der Beek et al.,1994).一般認為裂谷的側翼隆起會在拉張結束后的幾百萬年內(nèi)隨著熱衰減而消失(Weissel and Karner,1989;ten Brink and Stern,1992).而另有觀測數(shù)據(jù)顯示一些裂谷側翼的隆升可以長期維持,如非洲南部的大陡崖,有學者用具較大強度的巖石圈的撓曲均衡回彈來解釋這一現(xiàn)象(Weissel and Karner,1989;ten Brink and Stern,1992).但這些模型是在空盆載水的前提下討論的,忽略了沉積物的熱效應.本文重點討論在載水及載沉積物兩種情況下由側向熱傳導造成的盆地邊緣抬升的差異.

    圖4 固定與移動上邊界模型的構造熱演化結果對比圖a和b是拉張時間為10Ma、熱沉降時間為100Ma時的計算結果,圖c—g是拉張時間為40Ma、熱沉降時間為100Ma的結果,圖a和c為兩種模型計算的熱流差值的絕對值隨時間的變化,圖b和d為構造沉降量差值的絕對值,圖e—g為不同拉張系數(shù)下熱沉降結束時的溫度-深度曲線.Fig.4 Comparison of tectono-thermal evolution between fixed and moving upper boundary models(a)and(b)The extension lasts 10Ma.(c)—(g)The stretching period is 40Ma.In all figures,the thermal subsidence period lasts 100 Ma.(a)and(c)are absolute values of differences of the calculated heat flow between the two models.(b)and(d)are absolute values of differences of the calculated tectonic subsidence.(e),(f)and(g)are temperature-depth curves at the end of post-rift thermal subsidence stage with different stretching factors.

    圖5 載水和載沉積物模型的盆地側翼隆起效應拉張時間為10Ma,拉張系數(shù)為a圖所示,圖b—d為不同模型的構造沉降演化圖,圖e和f分別為拉張前x=-40km及x=-30km處兩點的抬升量隨時間的變化.Fig.5 Rift flank uplift of water and sediment loading modelsThe extension lasts 10Ma.(a)Distribution of stretching factor.(b),(c)and(d)Graphs of tectonic subsidence history for different models.(e)and(f)Variations of uplift with time at two points x=-40km and x=-30km,respectively.

    模擬采用二維非均勻拉張模型,拉張系數(shù)分布如圖5a所示,拉張10Ma,之后為熱沉降期.本文計算了載水、載熱導率分別為2.2W·m-1K-1及1.6W·m-1K-1的沉積物三組模型的構造沉降/抬升史.結果顯示,三組模型的盆地側翼隨著拉張逐漸抬升,但最大抬升量都并非出現(xiàn)在拉張剛結束時,而是在拉張結束后的6—10Ma內(nèi)側翼抬升到最高,之后慢慢消減(圖5b—f).上述計算中,載水模型的盆地側翼最大抬升量為273m,小于載沉積物模型的側翼最大抬升量;而熱導率為1.6W·m-1K-1的沉積物模型比熱導率為2.2W·m-1K-1的沉積物模型的側翼最大抬升量略大,二者分別為497m及489m(圖5e).模擬結果同時顯示,離拉張區(qū)域越遠,盆地側翼的抬升量越小(圖5e和f).值得注意的是,載水模型的盆地側翼隆升會在200Ma左右消失,而考慮了沉積物的模型其側翼隆升則會在200Ma左右穩(wěn)定于某個數(shù)值,也就是說側翼隆升并沒有隨時間消失.沉積物熱導率為2.2W·m-1K-1的模型的盆地側翼隆升最終維持在180m左右,而沉積物熱導率為1.6W·m-1K-1的模型的側翼隆升量則會穩(wěn)定于約245m(圖5e).這可能是由于沉積物與盆地基底熱導率差異產(chǎn)生的熱折射效應造成的.由于沉積物的熱導率較小,盆地基底的熱導率大于沉積蓋層的,隆起處的溫度低于凹陷區(qū)的溫度,凹陷處的熱會向相鄰隆起處轉移,使得盆地凹陷處的熱流降低,而隆起處的熱流增加(熊亮萍和高維安,1982;熊亮萍和張菊明,1984;Wang et al.,1985).在沉積蓋層與基底的交界處,具有最大的水平熱流值.經(jīng)過長時間的熱松弛后,溫度場達到穩(wěn)定.載水模型的溫度場可以在一定時間的熱冷卻后恢復到初始值.而由于沉積層的熱效應,載沉積物的模型盆地邊緣的溫度明顯比載水模型的溫度高,在不考慮剝蝕的情況下,前者的盆地側翼抬升能維持在某一水平.沉積蓋層與基底熱導率差異越大,這種效應越明顯,也因此產(chǎn)生更大的盆地側翼抬升.

    6 結論

    (1)拉張模型中改變初始地殼、巖石圈厚度會對計算的構造沉降量及熱流值產(chǎn)生很大影響.模擬中可根據(jù)需要選取多組初始地殼、巖石圈厚度值,將模擬結果與觀測值進行對比以確定最佳初始值.

    (2)軟流圈對流對構造熱演化的影響很大.軟流圈純對流及純傳導模型是拉張模型的兩種極限,在流變邊界層中應同時考慮熱傳導與對流.得知地幔對流的努塞爾特數(shù)Nu時,可相應的判斷軟流圈等效熱導率的上限.

    (3)模型上邊界對構造熱演化的影響隨著拉張系數(shù)的增大而增強,固定與移動上邊界對熱沉降期的熱流及構造沉降量的影響比拉張期的大.兩種模型得到的溫度場有著較為明顯的差別.

    (4)沉積蓋層的熱效應對裂谷盆地側翼隆起影響顯著.載沉積物模型的盆地邊緣抬升可以長時間維持,而載水模型的側翼抬升則會在拉張結束后逐漸消失.

    Alvarez F,Virieux J,Le Pichon X.1984.Thermal consequences of lithosphere extension over continental margins:the initial stretching phase.Geophys.J.Int.,78(2):389-411,doi:10.1111/j.1365-246X.1984.tb01956.x.

    Artemieva I.2011.The Lithosphere:an Interdisciplinary Approach.New York:Cambridge Univ.Press.

    Bassi G.1991.Factors controlling the style of continental rifting:insights from numerical modelling.EarthPlanet.Sci.Lett.,105(4):430-452,doi:10.1016/0012-821X(91)90183-I.

    Bassi G,Sabadini R.1994.The importance of subduction for the modern stress field in the Tyrrhenian area.Geophys.Res.Lett.,21(5):329-332,doi:10.1029/93GL03543.

    Baur F,Littke R,Wielens H,et al.2010.Basin modeling meets rift analysis—A numerical modeling study from the Jeanne d′Arc basin,offshore Newfoundland,Canada.Mar.Pet.Geol.,27(3):585-599,doi:10.1016/j.marpetgeo.2009.06.003.

    Braun J,Beaumont C.1989.A physical explanation of the relation between flank uplifts and the breakup unconformity at rifted continental margins.Geology,17(8):760-764,doi:10.1130/0091-7613(1989)017〈0760:APEOTR〉2.3.CO;2.

    Buck W R.1986.Small-scale convection induced by passive rifting:the cause for uplift of rift shoulders.EarthPlanet.Sci.Lett.,77(3-4):362-372,doi:10.1016/0012-821X(86)90146-9.

    Buck W R,Martinez F,Steckler M S,et al.1988.Thermal consequences of lithospheric extension:pure and simple.Tectonics,7(2):213-234,doi:10.1029/TC007i002p00213.

    Chen L.2014.Stretching factor estimation for the long-duration and multi-stage continental extensional tectonics:Application to the Baiyun Sag in the northern margin of the South China Sea.Tectonophysics,611(2014):167-180,doi:10.1016/j.tecto.2013.11.026.

    Chen L,Song H B,Dong C Z,et al.2008.2-D strain rate forward modeling.ChineseJ.Geophys.(in Chinese),51(6):1819-1825,doi:10.3321/j.issn:0001-5733.2008.06.023.

    Chen L,Zhang Z J,Song H B,et al.2013.Numerical modeling of extensional sedimentary basin formation with MATLAB:Application to the northern margin of the South China Sea.Comput.Geosci.,51:153-165,doi:10.1016/j.cageo.2012.07.014.

    Chen Y J.2001.Thermal effects of gabbro accretion from a deeper second melt lens at the fast spreading East Pacific Rise.J.Geophys.Res.:SolidEarth,106(B5):8581-8588,doi:10.1029/2000JB900420.

    Chéry J,Lucazeau F,Daignières M,et al.1992.Large uplift of rift flanks:A genetic link with lithospheric rigidity?EarthPlanet.Sci.Lett.,112(1-4):195-211,doi:10.1016/0012-821X(92)90016-O.

    Clift P,Lin J,Barckhausen U.2002.Evidence of low flexural rigidity and low viscosity lower continental crust during continental break-up in the South China Sea.Mar.Pet.Geol.,19(8):951-970,doi:10.1016/S0264-8172(02)00108-3.

    Cochran J R,Buck W R.2001.Near-axis subsidence rates,hydrothermal circulation,and thermal structure of mid-ocean ridge crests.J.Geophys.Res.:SolidEarth,106(B9):19233-19258,doi:10.1029/2001JB000379.

    Davis M,Kusznir N.2004.Depth-dependent lithospheric stretching at rifted continental margins.// Proc.NSF Rifted Margins Theoretical Institute,1:192-137.

    Dehler S A,Keen C E,Rohr K M M.1997.Tectonic and thermal evolution of Queen Charlotte Basin:lithospheric deformation and subsidence models.BasinRes.,9(3):243-261,doi:10.1046/j.1365-2117.1997.00043.x.

    Fan T Y,An M J.2009.Discussion on thickness of the Tarim Craton lithosphere and its mechanism.GeologicalReview(in Chinese),55(3):329-334.

    Fernàndez M,Ranalli G.1997.The role of rheology in extensional basin formation modelling.Tectonophysics,282(1-4):129-145,doi:10.1016/S0040-1951(97)00216-3.

    He L J.1999.Mutiple tectono-thermal modeling of Liaohe Bain in the Cenozoic.ChineseJ.Geophys.(in Chinese),42(1):62-68.

    He L J.2000.Advance in tectono-thermal modelling of sedimentary basins.Adv.EarthSci.(in Chinese),15(6):661-665.

    He L J.2014.Numerical modeling of convective erosion and peridotite-melt interaction in big mantle wedge:Implications for the destruction of the North China Craton.J.Geophys.Res.:SolidEarth,119(4):3662-3677,doi:10.1002/2013JB010657.

    He L J,Wang J Y.2004.Tectono-thermal modelling of sedimentary basins with episodic extension and inversion,a case history of the Jiyang Basin,North China.BasinRes.,16(4):587-599,doi:10.1111/j.1365-2117.2004.00245.x.

    He L J,Wang J Y.2007.Tectono-thermal modeling of sedimentary basins:review and outlook.ProgressinGeophys.(in Chinese),22(4):1215-1219.

    He L J,Xiong L P,Wang J Y.2002.Heat flow and thermal modeling of the Yinggehai Basin,South China Sea.Tectonophysics,351(3):245-253,doi:10.1016/s0040-1951(02)00160-9.

    Hellinger S J,Sclater J G.1983.Some comments on two-layer extensional models for the evolution of sedimentary basins.J.Geophys.Res.:SolidEarth,88(B10):8251-8269,doi:10.1029/JB088iB10p08251.

    Huismans R,Beaumont C.2011.Depth-dependent extension,twostage breakup and cratonic underplating at rifted margins.Nature,473(7345):74-78,doi:10.1038/nature09988.

    Huismans R S,Beaumont C.2008.Complex rifted continental margins explained by dynamical models of depth-dependent lithospheric extension.Geology,36(2):163-166,doi:10.1130/G24231A.1.

    Jarvis G T,McKenzie D P.1980.Sedimentary basin formation with finite extension rates.EarthPlanet.Sci.Lett.,48(1):42-52,doi:10.1016/0012-821X(80)90168-5.

    Jaupart C,Mareschal J C.2011.Heat Generation and Transport in the Earth.New York:Cambridge Univ.Press.

    Keen C.1985.The dynamics of rifting:deformation of the lithosphere by active and passive driving forces.Geophys.J.Int.,80(1):95-120,doi:10.1111/j.1365-246X.1985.tb05080.x.

    Keen C E,Dehler S A.1993.Stretching and subsidence:Rifting of conjugate margins in the North Atlantic Region.Tectonics,12(5):1209-1229,doi:10.1029/93tc00915.

    Kooi H,Cloetingh S,Burrus J.1992.Lithospheric necking and regional isostasy at extensional basins 1.Subsidence and gravity modeling with an application to the Gulf of Lions Margin(SE France).J.Geophys.Res.:SolidEarth,97(B12):17553-17571,doi:10.1029/92jb01377.

    Kusznir N,Ziegler P.1992.The mechanics of continental extension and sedimentary basin formation:a simple-shear/pure-shear flexural cantilever model.Tectonophysics,215(1):117-131,doi:10.1016/0040-1951(92)90077-J.

    Kusznir N,Karner G.2007.Continental lithospheric thinning and breakup in response to upwelling divergent mantle flow:application to the Woodlark,Newfoundland and Iberia margins.Geol.Soc.Spec.Publ.,282(1):389-419,doi:10.1144/SP282.16.Lachenbruch A H.1978.Heat flow in the Basin and Range province and thermal effects of tectonic extension.PureAppl.Geophys.,117(1-2):34-50,doi:10.1007/BF00879732.

    Leroy M,Gueydan F,Dauteuil O.2008.Uplift and strength evolution of passive margins inferred from 2-D conductive modelling.Geophys.J.Int.,172(1):464-476,doi:10.1111/j.1365-246X.2007.03566.x.

    Lin C S,Zhang Y M,Li S T,et al.2002.Quantitative modelling of multiphase lithospheric stretching and deep thermal history of some Tertiary rift basins in eastern China.ActaGeol.Sin.-Engl.,76(3):324-330,doi:10.1111/j.1755-6724.2002.tb00548.x.

    Lin G,Peng M L,Zhao C B,et al.2009.Numerical analysis and simulation experiment of lithospheric thermal structures in the South China Sea and the Western Pacific.J.EarthSci.,20(1):85-94,doi:10.1007/s12583-009-0009-0.

    McKenzie D.1978.Some remarks on the development of sedimentary basins.EarthPlanet.Sci.Lett.,40(1):25-32,doi:10.1016/0012-821X(78)90071-7.

    Morgan J P,Chen Y J.1993.The genesis of oceanic crust:Magma injection,hydrothermal circulation,and crustal flow.J.Geophys.Res.:SolidEarth,98(B4):6283-6297,doi:10.1029/92JB02650.

    Negredo A,F(xiàn)ernandez M,Torné M,et al.1999.Numerical modeling of simultaneous extension and compression:The Valencia Trough(western Mediterranean).Tectonics,18(2):361-374,doi:10.1029/1998TC900026.

    Parsons B,Sclater J G.1977.An analysis of the variation of ocean floor bathymetry and heat flow with age.J.Geophys.Res.,82(5):803-827,doi:10.1029/JB082i005p00803.

    Polyansky O P.2002.Dynamic causes for the opening of the Baikal Rift Zone:a numerical modelling approach.Tectonophysics,351(1-2):91-117,doi:10.1016/S0040-1951(02)00127-0.

    Rosenbaum G,Weinberg R F,Regenauer-Lieb K.2008.The geodynamics of lithospheric extension.Tectonophysics,458(1-4):1-8,doi:10.1016/j.tecto.2008.07.016.

    Rowley D B,Sahagian D.1986.Depth-dependent stretching:a different approach.Geology,14(1):32-35,doi:10.1130/0091-7613(1986)14〈32:DSADA〉2.0.CO;2.

    Royden L,Keen C.1980.Rifting process and thermal evolution of the continental margin of eastern Canada determined from subsidence curves.EarthPlanet.Sci.Lett.,51(2):343-361,doi:10.1016/0012-821X(80)90216-2.

    Schmeling H,Marquart G.2014.A scaling law for approximating porous hydrothermal convection by an equivalent thermal conductivity:theory and application to the cooling oceanic lithosphere.Geophys.J.Int.,197(2):645-664,doi:10.1093/gji/ggu022.

    Schubert G,Turcotte D L,Olson P.2001.Mantle Convection in the Earth and Planets.Cambridge:Cambridge Univ.Press.

    Shillington D J,White N,Minshull T A,et al.2008.Cenozoic evolution of the eastern Black Sea:A test of depth-dependent stretching models.EarthPlanet.Sci.Lett.,265(3-4):360-378,doi:10.1016/j.epsl.2007.10.033.

    Sleep N H.2006.Mantle plumes from top to bottom.EarthSci.Rev.,77(4):231-271,doi:10.1016/j.earscirev.2006.03.007.

    Sonder L J,England P C.1989.Effects of a temperature-dependent rheology on large-scale continental extension.J.Geophys.Res.:SolidEarth,94(B6):7603-7619,doi:10.1029/JB094iB06p07603.

    Song H B,Chen L,Zhang J,et al.2010.A MATLAB program for 1Dstrain rate inversion.Comput.Geosci.,36(1):16-23,doi:10.1016/j.cageo.2009.08.002.

    Spinelli G A,Harris R N.2011.Effects of the legacy of axial cooling on partitioning of hydrothermal heat extraction from oceanic lithosphere.J.Geophys.Res.:SolidEarth,116(B9):B09102,doi:10.1029/2011JB008248.

    Steckler M S.1985.Uplift and extension at the Gulf of Suez:indications of induced mantle convection.Nature,317(6033):135-139,doi:10.1038/317135a0.

    ten Brink U,Stern T.1992.Rift flank uplifts and hinterland basins:comparison of the Transantarctic Mountains with the Great Escarpment of southern Africa.J.Geophys.Res.:SolidEarth,97(B1):569-585,doi:10.1029/91JB02231.

    Turcotte D L,Schubert G.2002.Geodynamics.New York:Cambridge Univ.Press.

    Van Avendonk H J A,Lavier L L,Shillington D J,et al.2009.Extension of continental crust at the margin of the eastern Grand Banks,Newfoundland.Tectonophysics,468(1-4):131-148,doi:10.1016/j.tecto.2008.05.030.

    van der Beek P,Cloetingh S,Andriessen P.1994.Mechanisms of extensional basin formation and vertical motions at rift flanks:Constraints from tectonic modelling and fission-track thermochronology.EarthPlanet.Sci.Lett.,121(3-4):417-433,doi:10.1016/0012-821X(94)90081-7.

    van Wijk J,van Hunen J,Goes S.2008.Small-scale convection during continental rifting:Evidence from the Rio Grande rift.Geology,36(7):575-578,doi:10.1130/G24691A.1.

    Wang J Y,Wang J A,Xiong L P,et al.1985.Analysis of factors affecting heat flow density determination in the Liaohe Basin,North China.Tectonophysics,121(1):63-78,doi:10.1016/0040-1951(85)90268-9.

    Wang K L,He L J,Xiong L P,et al.1996.Thermal history model of a passive continental margin:Northern South China Sea.//Wang J Y,Li N S,eds.Geothermics in China.Beijing:Seismological Press,275-292.

    Wangen M,F(xiàn)jeldskaar W,F(xiàn)aleide J I,et al.2008.Forward modeling of stretching episodes and paleo heat flow of the Voring margin,NE Atlantic.J.Geodyn.,45(2-3):83-98,doi:10.1016/j.jog.2007.07.002.

    Weissel J K,Karner G D.1989.Flexural uplift of rift flanks due to mechanical unloading of the lithosphere during extension.J.Geophys.Res.:SolidEarth,94(B10):13919-13950,doi:10.1029/JB094iB10p13919.

    Wernicke B.1985.Uniform-sense normal simple shear of the continental lithosphere.Can.J.EarthSci.,22(1):108-125,doi:10.1139/e85-009.

    White N,McKenzie D.1988.Formation of the“steer′s head”geometry of sedimentary basins by differential stretching of the crust and mantle.Geology,16(3):250-253,doi:10.1130/0091-7613(1988)016〈0250:FOTSSH〉2.3.CO;2.

    White N.1994.An inverse method for determining lithospheric strain rate variation on geological timescales.EarthPlanet.Sci.Lett.,122(3):351-371,doi:10.1016/0012-821X(94)90008-6.

    White N,Bellingham P.2002.A two-dimensional inverse model for extensional sedimentary basins 1.Theory.J.Geophys.Res.:SolidEarth,107(B10):ETG 17-1-ETG 17-20,doi:10.1029/2001JB000173.

    Xie X N,Müller R D,Li S T,et al.2006.Origin of anomalous subsidence along the Northern South China Sea margin and its relationship to dynamic topography.Mar.Pet.Geol.,23(7):745-765,doi:10.1016/j.marpetgeo.2006.03.004.

    Xiong L P,Gao W A.1982.Characteristics of geotherm in uplift and depression.ChineseJ.Geophys.(in Chinese),25(5):448-456.

    Xiong L P,Zhang J M.1984.Mathematical simulation of refraction and redistribution of heat flow.ChineseJ.Geol.(in Chinese),19(4):445-454.

    Yang S H,Shi Y L.2013.Numerical simulation of Ningwu Ice Cave preserved for thousands of years in Shanxi Province.Sci.Technol.Rev.(in Chinese),31(13):23-28,doi:10.3981/j.issn.1000-7857.2013.13.002.

    Zhao C Y,Song H B,Qian R Y,et al.2010.Research on tectonothermal evolution modeling method for superimoposed basin with the Jianghan Basin as an example.ChineseJ.Geophys.(in Chinese),53(1):128-137,doi:10.3969/j.issn.0001-5733.2010.01.014.

    附中文參考文獻

    陳林,宋海斌,董崇志等.2008.二維應變速率正演模擬.地球物理學報,51(6):1819-1825,doi:10.3321/j.issn:0001-5733.2008.06.023.

    范桃園,安美建.2009.現(xiàn)今塔里木克拉通巖石圈厚度分析及機制探討.地質(zhì)論評,55(3):329-334.

    何麗娟.1999.遼河盆地新生代多期構造熱演化模擬.地球物理學報,42(1):62-68.

    何麗娟.2000.沉積盆地構造熱演化模擬的研究進展.地球科學進展,15(6):661-665.

    何麗娟,汪集暘.2007.沉積盆地構造熱演化研究進展:回顧與展望.地球物理學進展,22(4):1215-1219.

    熊亮萍,高維安.1982.隆起與拗陷地區(qū)地溫場的特點.地球物理學報,25(5):448-456.

    熊亮萍,張菊明.1984.熱流的折射和再分配的數(shù)學模擬.地質(zhì)科學,19(4):445-454.

    楊少華,石耀霖.2013.山西省寧武縣萬年冰洞持續(xù)存在機制數(shù)值模擬.科技導報,31(13):23-28,doi:10.3981/j.issn.1000-7857.2013.13.002.

    趙長煜,宋海斌,錢榮毅等.2010.疊合盆地構造熱演化模擬方法研究——以江漢盆地為例.地球物理學報,53(1):128-137,doi:10.3969/j.issn.0001-5733.2010.01.014.

    猜你喜歡
    側翼巖石圈熱導率
    空位缺陷對單層石墨烯導熱特性影響的分子動力學
    第四章 堅硬的巖石圈
    連續(xù)碳纖維鋁基復合材料橫向等效熱導率的模擬分析
    Si3N4/BN復合陶瓷熱導率及其有限元分析
    陶瓷學報(2020年5期)2020-11-09 09:23:04
    適用于張口受限患者的可調(diào)式印模托盤設計
    巖石圈磁場異常變化與巖石圈結構的關系
    地震研究(2017年3期)2017-11-06 21:54:14
    2014年魯?shù)?—5級地震相關斷裂的巖石圈磁異常分析
    地震研究(2017年3期)2017-11-06 01:58:51
    YB45型硬盒包裝機商標紙上膠機構的改進
    金屬熱導率的第一性原理計算方法在鋁中的應用
    蘆山7.0級地震前后巖石圈磁場異常變化研究
    地震研究(2014年1期)2014-02-27 09:29:41
    国产成+人综合+亚洲专区| 男女做爰动态图高潮gif福利片| АⅤ资源中文在线天堂| 最近在线观看免费完整版| 精品99又大又爽又粗少妇毛片 | 一区二区三区高清视频在线| 97超视频在线观看视频| 免费看十八禁软件| av黄色大香蕉| 亚洲专区国产一区二区| 午夜福利欧美成人| 俄罗斯特黄特色一大片| 老司机午夜福利在线观看视频| 欧美高清成人免费视频www| 每晚都被弄得嗷嗷叫到高潮| 久久精品影院6| 亚洲五月婷婷丁香| 高潮久久久久久久久久久不卡| 小说图片视频综合网站| 黄片小视频在线播放| 亚洲aⅴ乱码一区二区在线播放| 成人鲁丝片一二三区免费| 色综合婷婷激情| 一进一出抽搐gif免费好疼| 宅男免费午夜| 久久午夜亚洲精品久久| 国产高清三级在线| 亚洲电影在线观看av| 色综合亚洲欧美另类图片| 99国产精品一区二区蜜桃av| 欧美国产日韩亚洲一区| 全区人妻精品视频| 美女黄网站色视频| 欧美色视频一区免费| 男女那种视频在线观看| 麻豆久久精品国产亚洲av| 久久中文看片网| 一边摸一边抽搐一进一小说| 国产高清三级在线| 欧美一区二区精品小视频在线| 十八禁网站免费在线| 黄片大片在线免费观看| 国产 一区 欧美 日韩| 成人永久免费在线观看视频| 黄色成人免费大全| 超碰成人久久| 亚洲电影在线观看av| 亚洲av成人不卡在线观看播放网| 日本与韩国留学比较| 国产主播在线观看一区二区| 波多野结衣高清作品| 黄频高清免费视频| tocl精华| 国产成人欧美在线观看| 亚洲精品乱码久久久v下载方式 | 亚洲国产高清在线一区二区三| 精品国产美女av久久久久小说| av视频在线观看入口| 久久草成人影院| 少妇熟女aⅴ在线视频| 无遮挡黄片免费观看| 国产精品女同一区二区软件 | 日本 av在线| 免费人成视频x8x8入口观看| 亚洲avbb在线观看| 国语自产精品视频在线第100页| 男女午夜视频在线观看| 91在线精品国自产拍蜜月 | 国产毛片a区久久久久| 一级a爱片免费观看的视频| 亚洲av第一区精品v没综合| 精品久久蜜臀av无| 99视频精品全部免费 在线 | 人人妻人人看人人澡| 大型黄色视频在线免费观看| 我要搜黄色片| 狂野欧美激情性xxxx| 少妇熟女aⅴ在线视频| 少妇丰满av| 欧美黑人巨大hd| 久久伊人香网站| 精品不卡国产一区二区三区| 国产黄a三级三级三级人| 在线观看免费午夜福利视频| 国产亚洲精品一区二区www| 久久国产乱子伦精品免费另类| 午夜亚洲福利在线播放| 亚洲欧美精品综合久久99| 国产精品一区二区三区四区免费观看 | 亚洲中文字幕一区二区三区有码在线看 | 欧美3d第一页| 不卡av一区二区三区| 国产又色又爽无遮挡免费看| 性色avwww在线观看| 男女那种视频在线观看| av中文乱码字幕在线| 中亚洲国语对白在线视频| 日本在线视频免费播放| 两个人的视频大全免费| 午夜福利免费观看在线| 色综合欧美亚洲国产小说| 黄色日韩在线| 国产午夜精品久久久久久| 在线免费观看不下载黄p国产 | 日本 欧美在线| 中文在线观看免费www的网站| 欧美三级亚洲精品| 舔av片在线| 成人鲁丝片一二三区免费| 成人国产综合亚洲| 国产精品一区二区三区四区久久| 欧美色欧美亚洲另类二区| 成人永久免费在线观看视频| 国产精品久久久久久精品电影| 老鸭窝网址在线观看| 伦理电影免费视频| 欧美+亚洲+日韩+国产| 免费看美女性在线毛片视频| 日韩欧美在线二视频| 在线观看免费视频日本深夜| 久久久久久久午夜电影| 精品久久久久久成人av| 我的老师免费观看完整版| 久久精品综合一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 国语自产精品视频在线第100页| 黄片小视频在线播放| 成人av在线播放网站| 精品久久久久久久人妻蜜臀av| 亚洲真实伦在线观看| 国产黄色小视频在线观看| 成年女人毛片免费观看观看9| 99国产精品一区二区三区| 精华霜和精华液先用哪个| 日日干狠狠操夜夜爽| 亚洲精品乱码久久久v下载方式 | АⅤ资源中文在线天堂| 99久久精品国产亚洲精品| 欧美日韩国产亚洲二区| 午夜免费观看网址| 国产日本99.免费观看| 日韩国内少妇激情av| 男女午夜视频在线观看| 嫩草影院入口| 国产伦精品一区二区三区四那| 国产精品影院久久| 99国产精品一区二区蜜桃av| 丁香六月欧美| 色在线成人网| 国产亚洲av嫩草精品影院| 一级毛片女人18水好多| 免费看a级黄色片| 欧美在线一区亚洲| 国产1区2区3区精品| 国产av在哪里看| 亚洲aⅴ乱码一区二区在线播放| 午夜久久久久精精品| 亚洲国产高清在线一区二区三| 亚洲熟妇熟女久久| 一卡2卡三卡四卡精品乱码亚洲| 成年女人看的毛片在线观看| 亚洲中文av在线| 性欧美人与动物交配| av福利片在线观看| 女生性感内裤真人,穿戴方法视频| 色视频www国产| 国产乱人伦免费视频| 一二三四社区在线视频社区8| 18禁黄网站禁片午夜丰满| 日本熟妇午夜| 欧美黄色片欧美黄色片| 欧美黄色淫秽网站| 99视频精品全部免费 在线 | 欧美日韩精品网址| 亚洲中文av在线| 免费看日本二区| 欧美日韩综合久久久久久 | 国产亚洲欧美98| 在线观看66精品国产| 成年女人永久免费观看视频| 日韩欧美在线二视频| 两个人的视频大全免费| 精品久久久久久久末码| 香蕉久久夜色| 日韩三级视频一区二区三区| 欧美成人免费av一区二区三区| 精品国产亚洲在线| 听说在线观看完整版免费高清| 丰满人妻熟妇乱又伦精品不卡| 日韩免费av在线播放| av视频在线观看入口| xxx96com| 在线观看一区二区三区| 欧美色视频一区免费| 超碰成人久久| 午夜福利免费观看在线| 亚洲国产精品久久男人天堂| 999久久久国产精品视频| 男女之事视频高清在线观看| 天堂av国产一区二区熟女人妻| 无遮挡黄片免费观看| 这个男人来自地球电影免费观看| 免费在线观看日本一区| av国产免费在线观看| 一本一本综合久久| 亚洲色图av天堂| 麻豆国产97在线/欧美| 成熟少妇高潮喷水视频| 999久久久国产精品视频| 色播亚洲综合网| 日日夜夜操网爽| 国产一区二区在线观看日韩 | 国产单亲对白刺激| 国产爱豆传媒在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产激情欧美一区二区| 久久中文看片网| 日本黄色片子视频| 亚洲电影在线观看av| 少妇裸体淫交视频免费看高清| xxx96com| 日韩欧美一区二区三区在线观看| 99久国产av精品| 国产高潮美女av| 51午夜福利影视在线观看| 欧美一级a爱片免费观看看| 久久国产乱子伦精品免费另类| 给我免费播放毛片高清在线观看| 岛国在线观看网站| 欧美日韩中文字幕国产精品一区二区三区| av黄色大香蕉| 亚洲成人精品中文字幕电影| bbb黄色大片| av天堂中文字幕网| 91在线精品国自产拍蜜月 | 97碰自拍视频| 毛片女人毛片| 日本成人三级电影网站| 日日干狠狠操夜夜爽| 国产伦精品一区二区三区四那| 九九久久精品国产亚洲av麻豆 | 麻豆av在线久日| 中文字幕久久专区| 免费电影在线观看免费观看| 狂野欧美激情性xxxx| 欧美日韩综合久久久久久 | 少妇的逼水好多| 一进一出好大好爽视频| 色综合欧美亚洲国产小说| 国产乱人视频| 亚洲成av人片免费观看| 国产亚洲av高清不卡| 国产一区二区三区在线臀色熟女| 国产精品电影一区二区三区| 99热只有精品国产| 亚洲欧美精品综合久久99| 亚洲专区中文字幕在线| 最近在线观看免费完整版| 国产黄片美女视频| 最近最新中文字幕大全免费视频| 国产成+人综合+亚洲专区| 一二三四社区在线视频社区8| www.精华液| 人人妻,人人澡人人爽秒播| 日本一二三区视频观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲熟妇熟女久久| 国产人伦9x9x在线观看| 国产伦精品一区二区三区视频9 | 国产成人av激情在线播放| 国内揄拍国产精品人妻在线| 中出人妻视频一区二区| 久久欧美精品欧美久久欧美| 成人av在线播放网站| 女同久久另类99精品国产91| 母亲3免费完整高清在线观看| netflix在线观看网站| 88av欧美| 国产在线精品亚洲第一网站| 一进一出好大好爽视频| 色尼玛亚洲综合影院| 亚洲自拍偷在线| 在线观看免费午夜福利视频| 亚洲中文字幕日韩| 国产精品久久视频播放| 精品福利观看| 99久久无色码亚洲精品果冻| 婷婷亚洲欧美| 日本a在线网址| 中文字幕av在线有码专区| 黄片小视频在线播放| 亚洲男人的天堂狠狠| 无遮挡黄片免费观看| 国产成人aa在线观看| 日本黄色视频三级网站网址| www.www免费av| 国产精品一区二区三区四区免费观看 | 两个人视频免费观看高清| 午夜日韩欧美国产| a级毛片a级免费在线| 一级毛片女人18水好多| 国产精品一区二区三区四区久久| 久久久久久久久免费视频了| 18禁国产床啪视频网站| svipshipincom国产片| 成人三级黄色视频| 欧美黄色淫秽网站| 午夜福利18| bbb黄色大片| 亚洲成人久久性| 俺也久久电影网| 亚洲va日本ⅴa欧美va伊人久久| a级毛片a级免费在线| 啦啦啦韩国在线观看视频| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 亚洲成人中文字幕在线播放| 久久久国产成人精品二区| 日韩人妻高清精品专区| 欧美激情在线99| 精品国产乱子伦一区二区三区| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 欧美高清成人免费视频www| 国产黄a三级三级三级人| 亚洲av成人精品一区久久| 久久精品国产亚洲av香蕉五月| 特大巨黑吊av在线直播| 亚洲成人中文字幕在线播放| 老司机福利观看| 少妇人妻一区二区三区视频| 又爽又黄无遮挡网站| 国产男靠女视频免费网站| 热99re8久久精品国产| 一进一出抽搐gif免费好疼| 香蕉国产在线看| 一区二区三区激情视频| 黄片大片在线免费观看| 变态另类成人亚洲欧美熟女| 成人无遮挡网站| 日韩国内少妇激情av| 精品国产乱码久久久久久男人| 国产亚洲精品综合一区在线观看| 国产一区二区激情短视频| 免费在线观看日本一区| 亚洲中文av在线| 啦啦啦韩国在线观看视频| www.999成人在线观看| 国产又黄又爽又无遮挡在线| 免费电影在线观看免费观看| 91老司机精品| 国产欧美日韩精品亚洲av| 夜夜躁狠狠躁天天躁| 国产91精品成人一区二区三区| 伦理电影免费视频| 97超视频在线观看视频| 中文字幕久久专区| 久久久久久国产a免费观看| 99热精品在线国产| 国产精品日韩av在线免费观看| 国产av麻豆久久久久久久| 九色成人免费人妻av| 看免费av毛片| 琪琪午夜伦伦电影理论片6080| 美女被艹到高潮喷水动态| 99在线人妻在线中文字幕| 国产成人av激情在线播放| 男女床上黄色一级片免费看| 免费在线观看日本一区| 久久久国产精品麻豆| 波多野结衣巨乳人妻| 成年女人永久免费观看视频| netflix在线观看网站| 最近视频中文字幕2019在线8| 99国产精品99久久久久| 国产精品自产拍在线观看55亚洲| 免费在线观看视频国产中文字幕亚洲| 日韩欧美免费精品| 日本熟妇午夜| 亚洲五月天丁香| 国产精品一区二区精品视频观看| 性色avwww在线观看| 国产精品亚洲一级av第二区| 久久午夜亚洲精品久久| 18禁观看日本| 日本撒尿小便嘘嘘汇集6| 人人妻,人人澡人人爽秒播| xxxwww97欧美| 美女午夜性视频免费| 人妻久久中文字幕网| 男人舔女人下体高潮全视频| 伦理电影免费视频| 亚洲国产欧洲综合997久久,| 免费一级毛片在线播放高清视频| 国产精品精品国产色婷婷| 好看av亚洲va欧美ⅴa在| 国产伦人伦偷精品视频| 国产精品国产高清国产av| 99久久无色码亚洲精品果冻| 午夜两性在线视频| 久久久色成人| 波多野结衣高清无吗| 午夜激情欧美在线| 亚洲成人精品中文字幕电影| 精品99又大又爽又粗少妇毛片 | 麻豆一二三区av精品| 90打野战视频偷拍视频| 两性午夜刺激爽爽歪歪视频在线观看| 日韩av在线大香蕉| 又紧又爽又黄一区二区| 99在线视频只有这里精品首页| 精品电影一区二区在线| 国产视频内射| www日本黄色视频网| 色av中文字幕| 久久中文看片网| 中文字幕久久专区| 嫩草影院精品99| 国产精品亚洲美女久久久| 日本a在线网址| 国产黄片美女视频| 一本一本综合久久| 午夜福利在线在线| 欧美成人免费av一区二区三区| 婷婷精品国产亚洲av| 欧美xxxx黑人xx丫x性爽| 国产亚洲av高清不卡| 香蕉国产在线看| 日韩人妻高清精品专区| 国产单亲对白刺激| 国产精品久久久久久精品电影| 久久热在线av| 色综合站精品国产| 搞女人的毛片| 一边摸一边抽搐一进一小说| 成人av一区二区三区在线看| 最新在线观看一区二区三区| 一a级毛片在线观看| 日韩欧美国产在线观看| 国产一区二区在线av高清观看| 午夜日韩欧美国产| 欧美zozozo另类| 国产亚洲精品久久久久久毛片| 久久午夜亚洲精品久久| 精品人妻1区二区| 中出人妻视频一区二区| 国产精品一及| 国产伦精品一区二区三区四那| 久久中文字幕人妻熟女| 网址你懂的国产日韩在线| 禁无遮挡网站| av天堂在线播放| 99久久综合精品五月天人人| 免费电影在线观看免费观看| 波多野结衣巨乳人妻| 黄频高清免费视频| 九九在线视频观看精品| 狂野欧美白嫩少妇大欣赏| 村上凉子中文字幕在线| 亚洲第一电影网av| 国产成人福利小说| 一个人看视频在线观看www免费 | 精品一区二区三区av网在线观看| 综合色av麻豆| 国产乱人伦免费视频| 日韩有码中文字幕| 国产av在哪里看| 校园春色视频在线观看| 18禁国产床啪视频网站| 看黄色毛片网站| 最近最新中文字幕大全免费视频| 精品国产亚洲在线| 国产视频内射| 在线永久观看黄色视频| 精品国内亚洲2022精品成人| 性色av乱码一区二区三区2| 国产乱人伦免费视频| 后天国语完整版免费观看| 亚洲av成人不卡在线观看播放网| 特大巨黑吊av在线直播| 亚洲18禁久久av| 亚洲av成人一区二区三| 国产精品 国内视频| 黄色 视频免费看| 亚洲成a人片在线一区二区| 最新中文字幕久久久久 | 草草在线视频免费看| 啦啦啦观看免费观看视频高清| 黄色视频,在线免费观看| 国产精品 国内视频| 亚洲人成网站在线播放欧美日韩| 女警被强在线播放| 午夜a级毛片| 怎么达到女性高潮| 性色av乱码一区二区三区2| 国产亚洲av高清不卡| 黄色片一级片一级黄色片| 成人国产一区最新在线观看| 国产麻豆成人av免费视频| 久久香蕉精品热| av天堂在线播放| 男人的好看免费观看在线视频| 亚洲欧美精品综合一区二区三区| 亚洲欧美一区二区三区黑人| 国产精品电影一区二区三区| 亚洲av成人av| 日本一本二区三区精品| 丰满人妻一区二区三区视频av | 99久久成人亚洲精品观看| 国产三级黄色录像| 久久久久国产精品人妻aⅴ院| 俺也久久电影网| av天堂在线播放| 久久久久精品国产欧美久久久| 琪琪午夜伦伦电影理论片6080| 色老头精品视频在线观看| 男人和女人高潮做爰伦理| av中文乱码字幕在线| 伦理电影免费视频| 天天添夜夜摸| 麻豆成人av在线观看| 一进一出抽搐动态| 婷婷六月久久综合丁香| 人人妻人人看人人澡| 亚洲五月天丁香| 久久精品人妻少妇| 成人18禁在线播放| 国产一区二区三区在线臀色熟女| 在线a可以看的网站| 国产美女午夜福利| 久久婷婷人人爽人人干人人爱| 天堂av国产一区二区熟女人妻| 免费人成视频x8x8入口观看| 999久久久国产精品视频| 99热这里只有精品一区 | 国产在线精品亚洲第一网站| 国产伦人伦偷精品视频| 身体一侧抽搐| 成熟少妇高潮喷水视频| 成年女人永久免费观看视频| 国产精品1区2区在线观看.| 午夜激情福利司机影院| 九九久久精品国产亚洲av麻豆 | 老汉色av国产亚洲站长工具| 色老头精品视频在线观看| 在线观看日韩欧美| 午夜a级毛片| 中文字幕高清在线视频| 这个男人来自地球电影免费观看| 国产野战对白在线观看| 国产成人av激情在线播放| 亚洲七黄色美女视频| 国产高潮美女av| 欧美中文综合在线视频| 精品久久久久久久毛片微露脸| 亚洲人成网站在线播放欧美日韩| 欧美3d第一页| 天堂√8在线中文| 日韩精品中文字幕看吧| 国产麻豆成人av免费视频| 在线永久观看黄色视频| 好看av亚洲va欧美ⅴa在| 国产亚洲av嫩草精品影院| 国产精品久久久av美女十八| 日日夜夜操网爽| 中亚洲国语对白在线视频| 国产男靠女视频免费网站| 亚洲人成网站高清观看| 两个人视频免费观看高清| 精品电影一区二区在线| 亚洲欧洲精品一区二区精品久久久| 国产精品av久久久久免费| 亚洲电影在线观看av| 婷婷丁香在线五月| 色综合站精品国产| 午夜两性在线视频| 久久久久国产一级毛片高清牌| 国产午夜精品久久久久久| av天堂中文字幕网| 首页视频小说图片口味搜索| 亚洲精品美女久久av网站| 亚洲午夜精品一区,二区,三区| 国产成人aa在线观看| 久久天堂一区二区三区四区| 国产亚洲精品一区二区www| 好男人在线观看高清免费视频| 一区福利在线观看| 国产亚洲精品综合一区在线观看| 青草久久国产| 国产一区在线观看成人免费| 中文字幕av在线有码专区| 日本一二三区视频观看| 99re在线观看精品视频| 人人妻人人看人人澡| 免费观看的影片在线观看| 国产一区二区在线观看日韩 | 免费av不卡在线播放| 97超级碰碰碰精品色视频在线观看| 亚洲18禁久久av| 国产黄片美女视频| 国产69精品久久久久777片 | 色老头精品视频在线观看| x7x7x7水蜜桃| 国产精品一区二区三区四区久久| 国产乱人伦免费视频| 最新在线观看一区二区三区| 日韩三级视频一区二区三区| 黑人巨大精品欧美一区二区mp4| 欧美中文日本在线观看视频| 国产三级黄色录像| 亚洲成人免费电影在线观看| 99久国产av精品| 国产精品爽爽va在线观看网站| 动漫黄色视频在线观看| 蜜桃久久精品国产亚洲av|