• <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
    久久国内精品自在自线图片| 一个人免费看片子| 人人妻人人添人人爽欧美一区卜| 欧美另类一区| 久久久国产一区二区| 一本色道久久久久久精品综合| 亚洲怡红院男人天堂| 国产精品国产三级专区第一集| 三级国产精品欧美在线观看| 我的女老师完整版在线观看| 91在线精品国自产拍蜜月| 简卡轻食公司| 丁香六月天网| 伦理电影大哥的女人| 亚洲精品亚洲一区二区| 国产男女内射视频| 亚洲精品国产av成人精品| 欧美丝袜亚洲另类| 午夜福利,免费看| 国产免费一级a男人的天堂| 亚洲av成人精品一二三区| 久久ye,这里只有精品| 久久国内精品自在自线图片| tube8黄色片| 中国美白少妇内射xxxbb| 男女免费视频国产| 精品午夜福利在线看| 在线观看免费高清a一片| 麻豆乱淫一区二区| h视频一区二区三区| 99精国产麻豆久久婷婷| 欧美日韩精品成人综合77777| √禁漫天堂资源中文www| 亚洲国产最新在线播放| 狠狠精品人妻久久久久久综合| 少妇的逼水好多| 免费观看av网站的网址| 免费播放大片免费观看视频在线观看| 日日摸夜夜添夜夜爱| 99热6这里只有精品| 精品一区二区三卡| 国产女主播在线喷水免费视频网站| 一区二区av电影网| 啦啦啦中文免费视频观看日本| 久久精品夜色国产| av专区在线播放| 一级毛片黄色毛片免费观看视频| 精品国产乱码久久久久久小说| 成人影院久久| 亚洲五月色婷婷综合| 老女人水多毛片| 看免费成人av毛片| 精品久久久久久久久亚洲| 国产视频首页在线观看| a级毛色黄片| 大又大粗又爽又黄少妇毛片口| .国产精品久久| 哪个播放器可以免费观看大片| 免费不卡的大黄色大毛片视频在线观看| 一区二区三区四区激情视频| 亚洲av不卡在线观看| 欧美老熟妇乱子伦牲交| 大香蕉久久成人网| 综合色丁香网| 久久精品国产亚洲网站| 欧美变态另类bdsm刘玥| freevideosex欧美| 看非洲黑人一级黄片| 国产精品国产三级专区第一集| 亚洲天堂av无毛| 欧美人与性动交α欧美精品济南到 | 亚洲情色 制服丝袜| 下体分泌物呈黄色| 亚洲人成网站在线观看播放| 精品国产露脸久久av麻豆| 日韩av不卡免费在线播放| 国产一区亚洲一区在线观看| 亚洲国产色片| av专区在线播放| 交换朋友夫妻互换小说| 黑人猛操日本美女一级片| a级毛片黄视频| 亚洲第一av免费看| 国产免费一级a男人的天堂| 国产欧美另类精品又又久久亚洲欧美| 精品一区二区三区视频在线| 国产精品国产av在线观看| 人成视频在线观看免费观看| 欧美国产精品一级二级三级| 三级国产精品片| 国产午夜精品一二区理论片| 超色免费av| 亚洲欧洲精品一区二区精品久久久 | 成人二区视频| 边亲边吃奶的免费视频| 又粗又硬又长又爽又黄的视频| 一本大道久久a久久精品| 久久久国产欧美日韩av| 永久免费av网站大全| 在线观看三级黄色| 九色亚洲精品在线播放| 亚洲性久久影院| 久久99精品国语久久久| 一区二区三区四区激情视频| 内地一区二区视频在线| 一级毛片黄色毛片免费观看视频| 免费久久久久久久精品成人欧美视频 | 免费av中文字幕在线| av女优亚洲男人天堂| 午夜影院在线不卡| 日韩人妻高清精品专区| 美女xxoo啪啪120秒动态图| 亚洲精品成人av观看孕妇| 国产午夜精品一二区理论片| 日本午夜av视频| 亚洲精品久久午夜乱码| 亚洲国产欧美日韩在线播放| 老司机影院毛片| 国产一区二区在线观看av| 美女内射精品一级片tv| 日日啪夜夜爽| 久久人妻熟女aⅴ| 亚州av有码| 2021少妇久久久久久久久久久| 日本91视频免费播放| 九九在线视频观看精品| 精品久久蜜臀av无| 一级黄片播放器| 久久亚洲国产成人精品v| 亚洲精品,欧美精品| 欧美性感艳星| 九九爱精品视频在线观看| 女人久久www免费人成看片| 亚洲欧洲国产日韩| 在现免费观看毛片| 亚洲三级黄色毛片| 在线观看免费日韩欧美大片 | 久久久国产精品麻豆| 亚洲欧洲国产日韩| 国产片特级美女逼逼视频| 日韩一区二区视频免费看| a级毛片在线看网站| 国产片内射在线| 天堂俺去俺来也www色官网| 久久久久久久久久成人| 色婷婷av一区二区三区视频| 精品一区二区三区视频在线| 国产又色又爽无遮挡免| 久久99精品国语久久久| 大片免费播放器 马上看| tube8黄色片| 交换朋友夫妻互换小说| 久久青草综合色| 精品99又大又爽又粗少妇毛片| 亚洲国产日韩一区二区| 中国三级夫妇交换| 青春草国产在线视频| 日韩,欧美,国产一区二区三区| 人人妻人人澡人人看| 午夜91福利影院| 日本91视频免费播放| 国产女主播在线喷水免费视频网站| 丰满乱子伦码专区| 亚洲精品乱码久久久v下载方式| 国产av码专区亚洲av| 我的女老师完整版在线观看| 视频中文字幕在线观看| 91在线精品国自产拍蜜月| 国产极品天堂在线| 亚洲综合色惰| av卡一久久| 99久久人妻综合| 国产成人精品福利久久| av播播在线观看一区| 观看美女的网站| 亚洲国产精品成人久久小说| 最近的中文字幕免费完整| 飞空精品影院首页| 亚洲第一av免费看| videosex国产| 国产精品偷伦视频观看了| 又粗又硬又长又爽又黄的视频| 99久久综合免费| 啦啦啦啦在线视频资源| 成年女人在线观看亚洲视频| 人人妻人人澡人人爽人人夜夜| av播播在线观看一区| 欧美另类一区| 简卡轻食公司| 涩涩av久久男人的天堂| 亚洲av国产av综合av卡| 校园人妻丝袜中文字幕| 亚洲第一av免费看| 精品国产乱码久久久久久小说| 亚洲欧美清纯卡通| 国产成人精品一,二区| 波野结衣二区三区在线| 最近手机中文字幕大全| 亚洲精品日本国产第一区| 女性被躁到高潮视频| 一区二区三区乱码不卡18| 人妻少妇偷人精品九色| 韩国高清视频一区二区三区| 99热这里只有精品一区| 国产毛片在线视频| 日韩人妻高清精品专区| 在线观看三级黄色| 日韩欧美精品免费久久| 中文字幕亚洲精品专区| 日韩成人伦理影院| 日韩人妻高清精品专区| 大码成人一级视频| 婷婷色综合www| 欧美成人午夜免费资源| 国产探花极品一区二区| 汤姆久久久久久久影院中文字幕| xxx大片免费视频| 春色校园在线视频观看| 在线精品无人区一区二区三| 菩萨蛮人人尽说江南好唐韦庄| 久久久久网色| 日韩熟女老妇一区二区性免费视频| av免费在线看不卡| 久久精品久久久久久久性| 一级毛片 在线播放| 99九九在线精品视频| 免费人成在线观看视频色| 欧美日韩在线观看h| 亚洲国产最新在线播放| 精品少妇黑人巨大在线播放| 高清av免费在线| 国产无遮挡羞羞视频在线观看| 夫妻性生交免费视频一级片| 观看美女的网站| tube8黄色片| 亚洲精品,欧美精品| 性色avwww在线观看| 考比视频在线观看| 成人国产av品久久久| 久久婷婷青草| 精品久久久久久电影网| 性色av一级| 99视频精品全部免费 在线| 777米奇影视久久| 建设人人有责人人尽责人人享有的| 亚洲成色77777| 国国产精品蜜臀av免费| 99久久人妻综合| 亚洲精品自拍成人| 夜夜骑夜夜射夜夜干| 久久久欧美国产精品| 日韩制服骚丝袜av| 国产无遮挡羞羞视频在线观看| 免费观看av网站的网址| 午夜福利视频精品| 久久久久久久久久成人| 丰满乱子伦码专区| 精品亚洲成a人片在线观看| 只有这里有精品99| 人妻少妇偷人精品九色| av.在线天堂| 日本av手机在线免费观看| 国产精品麻豆人妻色哟哟久久| 最新的欧美精品一区二区| 丰满少妇做爰视频| 久久久久国产网址| 9色porny在线观看| 久久久久久久亚洲中文字幕| 日本wwww免费看| 麻豆乱淫一区二区| 亚洲av欧美aⅴ国产| 91精品伊人久久大香线蕉| 精品一品国产午夜福利视频| 我的女老师完整版在线观看| 99热这里只有精品一区| av在线观看视频网站免费| 高清欧美精品videossex| 不卡视频在线观看欧美| a级毛片免费高清观看在线播放| 伊人久久国产一区二区| 日韩伦理黄色片| 久久精品夜色国产| 成人毛片60女人毛片免费| 日日摸夜夜添夜夜爱| 少妇被粗大的猛进出69影院 | 女人精品久久久久毛片| 午夜精品国产一区二区电影| 欧美精品一区二区免费开放| 亚洲人与动物交配视频| 亚洲国产精品一区二区三区在线| 王馨瑶露胸无遮挡在线观看| 亚洲精品美女久久av网站| 大片电影免费在线观看免费| 国产亚洲精品第一综合不卡 | 99九九线精品视频在线观看视频| 一区二区日韩欧美中文字幕 | 精品午夜福利在线看| 亚洲精品国产av蜜桃| 一级,二级,三级黄色视频| 日韩欧美精品免费久久| 成人免费观看视频高清| 免费观看性生交大片5| 国产av一区二区精品久久| 我的女老师完整版在线观看| 亚洲图色成人| 国产一区二区在线观看日韩| 99热这里只有是精品在线观看| 久久人妻熟女aⅴ| 久久国产亚洲av麻豆专区| 少妇的逼水好多| 各种免费的搞黄视频| 特大巨黑吊av在线直播| 亚洲熟女精品中文字幕| 国产69精品久久久久777片| 中文天堂在线官网| 国内精品宾馆在线| 美女脱内裤让男人舔精品视频| 日本爱情动作片www.在线观看| 欧美变态另类bdsm刘玥| 亚洲欧洲日产国产| 久久精品熟女亚洲av麻豆精品| 亚洲国产av影院在线观看| 丰满饥渴人妻一区二区三| 交换朋友夫妻互换小说| 五月天丁香电影| 久久国产精品大桥未久av| 国产成人av激情在线播放 | 亚洲伊人久久精品综合| 亚洲不卡免费看| 国产精品偷伦视频观看了| 丝袜脚勾引网站| 精品国产一区二区三区久久久樱花| 色视频在线一区二区三区| 大又大粗又爽又黄少妇毛片口| 国产成人91sexporn| 欧美亚洲日本最大视频资源| 欧美精品一区二区免费开放| 国产亚洲一区二区精品| 只有这里有精品99| 99热这里只有是精品在线观看| av网站免费在线观看视频| 一本大道久久a久久精品| 人人妻人人澡人人爽人人夜夜| 日本av手机在线免费观看| 另类亚洲欧美激情| 多毛熟女@视频| 久久久久久久亚洲中文字幕| 亚州av有码| 亚洲国产精品国产精品| www.色视频.com| 18禁观看日本| 精品国产乱码久久久久久小说| 丰满饥渴人妻一区二区三| 亚洲精品亚洲一区二区| 少妇 在线观看| 天天操日日干夜夜撸| 一级毛片我不卡| 久久鲁丝午夜福利片| 青春草亚洲视频在线观看| 欧美变态另类bdsm刘玥| av一本久久久久| 午夜免费观看性视频| 自拍欧美九色日韩亚洲蝌蚪91| 丁香六月天网| 免费不卡的大黄色大毛片视频在线观看| 99久久人妻综合| 久久久久精品久久久久真实原创| 久久久久久久国产电影| 日韩精品有码人妻一区| 久久精品国产亚洲av天美| 一区在线观看完整版| 另类精品久久| 国产精品一国产av| 久久久久久久大尺度免费视频| 欧美日韩在线观看h| a级毛色黄片| 亚洲精品久久久久久婷婷小说| 亚洲人与动物交配视频| av网站免费在线观看视频| 中文乱码字字幕精品一区二区三区| 高清午夜精品一区二区三区| 午夜福利影视在线免费观看| 欧美三级亚洲精品| av不卡在线播放| 夜夜看夜夜爽夜夜摸| 人妻人人澡人人爽人人| 精品熟女少妇av免费看| 国产精品久久久久成人av| 久久国产精品大桥未久av| 大片免费播放器 马上看| 女性生殖器流出的白浆| 极品人妻少妇av视频| 成年av动漫网址| 美女主播在线视频| 中文字幕制服av| 69精品国产乱码久久久| 欧美日韩亚洲高清精品| 久久精品国产自在天天线| 久久99热6这里只有精品| 最新的欧美精品一区二区| 秋霞在线观看毛片| 日本av免费视频播放| 青春草国产在线视频| av网站免费在线观看视频| 国产精品久久久久久精品电影小说| 日本黄大片高清| 五月天丁香电影| 国产成人精品无人区| 亚洲国产色片| 国产黄色免费在线视频| 免费看av在线观看网站| 亚洲天堂av无毛| 特大巨黑吊av在线直播| 九九久久精品国产亚洲av麻豆| 乱人伦中国视频| 美女国产视频在线观看| 天天影视国产精品| 多毛熟女@视频| 久久毛片免费看一区二区三区| 欧美另类一区| 国产成人一区二区在线| 天美传媒精品一区二区| 久久精品人人爽人人爽视色| 亚洲精品国产av成人精品| 日本-黄色视频高清免费观看| 一本一本综合久久| 男女啪啪激烈高潮av片| 热re99久久精品国产66热6| 欧美日韩视频精品一区| 纯流量卡能插随身wifi吗| 亚洲欧美成人综合另类久久久| 老熟女久久久| 国产精品人妻久久久影院| 成人综合一区亚洲| 男女高潮啪啪啪动态图| 高清视频免费观看一区二区| 亚洲色图综合在线观看| 国产男女超爽视频在线观看| 99久久综合免费| 国产熟女欧美一区二区| 女性被躁到高潮视频| 美女大奶头黄色视频| 少妇人妻精品综合一区二区| 一区二区三区精品91| 美女视频免费永久观看网站| 国产国语露脸激情在线看| 日日啪夜夜爽| 久久人人爽人人爽人人片va| 亚洲欧美成人精品一区二区| 一区二区三区四区激情视频| 99精国产麻豆久久婷婷| 久久久久精品性色| 日本与韩国留学比较| 如何舔出高潮| 国产成人免费观看mmmm| 在线观看免费视频网站a站| 男女高潮啪啪啪动态图| 免费看光身美女| 亚洲高清免费不卡视频| 久久ye,这里只有精品| 最近2019中文字幕mv第一页| 黑丝袜美女国产一区| 黄片播放在线免费| 久热这里只有精品99| 99久久中文字幕三级久久日本| 国产亚洲一区二区精品| 国产成人精品无人区| 亚洲高清免费不卡视频| 伦精品一区二区三区| 一级毛片电影观看| 一个人免费看片子| 18禁裸乳无遮挡动漫免费视频| 国产色婷婷99| 国产日韩一区二区三区精品不卡 | 五月玫瑰六月丁香| 日韩在线高清观看一区二区三区| 亚洲精品av麻豆狂野| 最近中文字幕2019免费版| 一本一本综合久久| 一区二区三区精品91| 久久久久久久国产电影| 久久久久视频综合| 免费看不卡的av| 观看美女的网站| .国产精品久久| 在现免费观看毛片| 大香蕉久久成人网| 亚洲av在线观看美女高潮| 免费不卡的大黄色大毛片视频在线观看| 精品一品国产午夜福利视频| 人人妻人人添人人爽欧美一区卜| 成人亚洲欧美一区二区av| 一本一本综合久久| 91国产中文字幕| 蜜桃久久精品国产亚洲av| 国产一区亚洲一区在线观看| 日韩成人伦理影院| 好男人视频免费观看在线| 满18在线观看网站| 日本黄色片子视频| 久久久久久伊人网av| 秋霞伦理黄片| 91在线精品国自产拍蜜月| 精品国产露脸久久av麻豆| 国产精品99久久99久久久不卡 | 亚洲一级一片aⅴ在线观看| 亚洲精品国产av蜜桃| 看非洲黑人一级黄片| 一个人免费看片子| 成人18禁高潮啪啪吃奶动态图 | 国产高清国产精品国产三级| 成人免费观看视频高清| 男人爽女人下面视频在线观看| 超色免费av| 久久精品国产a三级三级三级| 久久久国产一区二区| 久久久久久久久久成人| 亚洲精品国产av成人精品| 亚洲精品视频女| 国产爽快片一区二区三区| 亚洲综合色网址| 久久亚洲国产成人精品v| 亚洲高清免费不卡视频| 日韩人妻高清精品专区| 精品99又大又爽又粗少妇毛片| 亚洲国产欧美在线一区| 欧美激情国产日韩精品一区| 春色校园在线视频观看| 蜜臀久久99精品久久宅男| 自线自在国产av| 午夜激情久久久久久久| 精品少妇黑人巨大在线播放| 免费观看的影片在线观看| 2022亚洲国产成人精品| 欧美成人午夜免费资源| 亚洲伊人久久精品综合| 精品人妻一区二区三区麻豆| 桃花免费在线播放| 一个人免费看片子| 黑人高潮一二区| 夜夜爽夜夜爽视频| 在线观看免费日韩欧美大片 | 久久久久久久久久久久大奶| 五月天丁香电影| 人妻夜夜爽99麻豆av| 国产片特级美女逼逼视频| 色网站视频免费| av电影中文网址| 你懂的网址亚洲精品在线观看| 午夜91福利影院| 久久久久久久久久人人人人人人| 国产免费视频播放在线视频| 女人久久www免费人成看片| 日日摸夜夜添夜夜爱| 看十八女毛片水多多多| 国产成人a∨麻豆精品| 亚洲人成77777在线视频| 好男人视频免费观看在线| 欧美日韩综合久久久久久| av电影中文网址| 国产成人精品久久久久久| 丝袜脚勾引网站| 亚洲av中文av极速乱| 韩国高清视频一区二区三区| 街头女战士在线观看网站| 国产亚洲一区二区精品| 国产 一区精品| 女的被弄到高潮叫床怎么办| 国产熟女午夜一区二区三区 | 日本av免费视频播放| 久久久久久久久久久丰满| 伦理电影免费视频| 最近中文字幕2019免费版| 国产精品一二三区在线看| 七月丁香在线播放| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品久久久久久婷婷小说| 97在线人人人人妻| 9色porny在线观看| 亚洲欧美精品自产自拍| 美女视频免费永久观看网站| 一级a做视频免费观看| 性高湖久久久久久久久免费观看| 日韩 亚洲 欧美在线| 亚洲怡红院男人天堂| 日韩在线高清观看一区二区三区| 人妻人人澡人人爽人人| 中文欧美无线码| 欧美一级a爱片免费观看看| 国精品久久久久久国模美| 久久热精品热| 国产乱人偷精品视频| 久久久久人妻精品一区果冻| 国产精品久久久久久久电影| 日本黄色片子视频| 精品午夜福利在线看| 免费大片黄手机在线观看| 免费观看的影片在线观看| 亚洲av欧美aⅴ国产| 日韩一区二区三区影片| 免费高清在线观看视频在线观看| 欧美成人精品欧美一级黄| 亚洲国产精品专区欧美| tube8黄色片| 夫妻午夜视频| 日日摸夜夜添夜夜添av毛片| 日本av手机在线免费观看| 看十八女毛片水多多多| 最近手机中文字幕大全| 91久久精品电影网| 少妇 在线观看| 日本午夜av视频| 熟妇人妻不卡中文字幕| 国产精品秋霞免费鲁丝片| 午夜影院在线不卡|