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

    柴達(dá)木盆地現(xiàn)今大地?zé)崃髋c晚古生代以來構(gòu)造-熱演化

    2015-03-07 03:57:58李宗星高俊鄭策劉成林馬寅生趙為永
    地球物理學(xué)報(bào) 2015年10期
    關(guān)鍵詞:徑跡柴達(dá)木盆地熱導(dǎo)率

    李宗星, 高俊, 鄭策,3, 劉成林, 馬寅生, 趙為永

    1 中國地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所頁巖油氣調(diào)查評價(jià)重點(diǎn)實(shí)驗(yàn)室, 北京 100081 2 中國地質(zhì)大學(xué)(北京) 水資源與環(huán)境工程北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083 3 中國石油大學(xué)(北京), 北京 102249 4 中石油青海油田分公司勘探開發(fā)研究院, 甘肅敦煌 736202

    ?

    柴達(dá)木盆地現(xiàn)今大地?zé)崃髋c晚古生代以來構(gòu)造-熱演化

    李宗星1, 高俊2*, 鄭策1,3, 劉成林1, 馬寅生1, 趙為永4

    1 中國地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所頁巖油氣調(diào)查評價(jià)重點(diǎn)實(shí)驗(yàn)室, 北京 100081 2 中國地質(zhì)大學(xué)(北京) 水資源與環(huán)境工程北京市重點(diǎn)實(shí)驗(yàn)室, 北京 100083 3 中國石油大學(xué)(北京), 北京 102249 4 中石油青海油田分公司勘探開發(fā)研究院, 甘肅敦煌 736202

    依據(jù)鉆孔系統(tǒng)穩(wěn)態(tài)測溫、靜井溫度資料與實(shí)測熱導(dǎo)率數(shù)據(jù)分析了柴達(dá)木盆地地溫場分布特征,建立了柴達(dá)木盆地?zé)釋?dǎo)率柱,新增了17個(gè)大地?zé)崃鲾?shù)據(jù). 柴達(dá)木盆地現(xiàn)今地溫梯度介于17.1~38.6 ℃·km-1,平均為28.6±4.6 ℃·km-1,大地?zé)崃鹘橛?2.9~70.4 mW·m-2,平均 55.1±7.9 mW·m-2. 盆地不同構(gòu)造單元地溫場存在差異,昆北逆沖帶、一里坪坳陷屬于“高溫區(qū)”,祁南逆沖帶屬于“中溫區(qū)”,三湖坳陷、德令哈坳陷及歐龍布魯克隆起屬于“低溫區(qū)”,盆地現(xiàn)今地溫場分布特征受控于地殼深部結(jié)構(gòu)、盆地構(gòu)造等因素. 以現(xiàn)今地溫場為基礎(chǔ),采用磷灰石、鋯石裂變徑跡年齡分布特征定性分析與徑跡長度分布數(shù)據(jù)定量模擬相結(jié)合,研究了柴達(dá)木盆地晚古生代以來的沉積埋藏、抬升剝蝕和熱演化史,并結(jié)合區(qū)域構(gòu)造背景,對柴達(dá)木盆地構(gòu)造演化過程進(jìn)行了探討,研究表明柴達(dá)木盆地晚古生代以來經(jīng)歷了六期(254.0—199 Ma,177—148.6 Ma,87—62 Ma,41.1—33.6 Ma,9.6—7.1 Ma,2.9—1.8 Ma)構(gòu)造運(yùn)動,六期構(gòu)造事件與研究區(qū)構(gòu)造演化的動力學(xué)背景相吻合. 其中白堊紀(jì)末期(87—62 Ma)的構(gòu)造事件導(dǎo)致了柴達(dá)木盆地東部隆升并遭受剝蝕,歐龍布魯克隆起形成雛形,柴達(dá)木盆地北緣在弱擠壓環(huán)境下形成坳陷盆地; 中新世末的兩期構(gòu)造事件(9.6—7.1 Ma和2.9—1.8 Ma)使柴達(dá)木盆地遭受強(qiáng)烈擠壓,盆地快速隆升,構(gòu)造變形強(qiáng)烈,基本形成現(xiàn)今的構(gòu)造面貌.

    地溫梯度; 大地?zé)崃? 裂變徑跡; 構(gòu)造-熱演化; 柴達(dá)木盆地

    1 引言

    柴達(dá)木盆地是青藏高原內(nèi)部最大的山間盆地,盆地的構(gòu)造-熱演化特征對于研究青藏高原隆升機(jī)制有非常重要的意義(Yin et al.,2008). 柴達(dá)木盆地平均海拔為~2800 m,整體呈現(xiàn)三角形狀,西北邊界為阿爾金斷裂帶,東北邊界為祁連山—南山?jīng)_斷帶,南部邊界為東昆侖轉(zhuǎn)換擠壓系統(tǒng)(祁漫塔格—東昆侖逆沖帶和昆侖斷裂帶),具有獨(dú)特的巖石圈深部動力學(xué)與盆地構(gòu)造-熱演化特征. 柴達(dá)木盆地也是我國西部大型的含油氣盆地,地層發(fā)育較齊全,泥巖和膏巖相當(dāng)發(fā)育,形成良好蓋層的同時(shí)也對盆地現(xiàn)今地溫分布產(chǎn)生大的影響(邱楠生,2001). 盆地的大地?zé)崃骱蜆?gòu)造-熱演化史為區(qū)域動力學(xué)與盆地構(gòu)造演化研究提供了重要參數(shù),同時(shí)也為盆地進(jìn)行油氣資源評價(jià)提供了基礎(chǔ)數(shù)據(jù),直接服務(wù)于生產(chǎn).

    柴達(dá)木盆地構(gòu)造-熱演化研究工作還較薄弱. 前人對柴達(dá)木盆地西部3 km以淺的溫度分布(張業(yè)成等,1990; 王鈞等,1990)進(jìn)行了研究,計(jì)算了盆地西部地區(qū)部分鉆孔的大地?zé)崃?沈顯杰等,1994; 李國樺,1992; 邱楠生,2001). 任戰(zhàn)利(1993)在早期應(yīng)用流體包裹體,鏡質(zhì)體反射率數(shù)據(jù)對柴達(dá)木盆地西部的地?zé)嵫莼愤M(jìn)行探索性分析. Qiu(2002,2003)利用磷灰石裂變徑跡(AFT)參數(shù)和鏡質(zhì)體反射率動力學(xué)模型對柴達(dá)木盆地西部地區(qū)熱歷史進(jìn)行了研究. 近年來,不少的研究者利用磷灰石 (鋯石) 的裂變徑跡對柴達(dá)木盆地西部新生代沉積源區(qū)、構(gòu)造熱事件進(jìn)行了研究(王世明等,2008; 孫國強(qiáng)等,2009; 高軍平等,2011). 但總體來看,(1) 研究區(qū)域上,已開展的研究工作主要集中在盆地西部,其他區(qū)域開展的工作少,盆地東部尚處于空白; (2)地質(zhì)時(shí)代上,現(xiàn)有的工作主要側(cè)重盆地新生代以來的演化史,中生代時(shí)期的研究很少,古生代時(shí)期的更是尚未涉及; (3)研究領(lǐng)域上,不同區(qū)域的差異性熱演化史、盆地構(gòu)造成因與熱演化關(guān)系研究尚比較缺乏; (4)數(shù)據(jù)資料上,前期現(xiàn)今地溫場研究的基礎(chǔ)數(shù)據(jù)多為20世紀(jì)90年代之前,特別是低溫?zé)崮甏鷮W(xué)數(shù)據(jù)也與現(xiàn)今不可同日而語.

    近年來隨著油氣勘探工作的深入,在柴達(dá)木盆地部署了更多的鉆井,新增了一大批地溫?cái)?shù)據(jù),同時(shí)先進(jìn)的鉆孔溫度連續(xù)采集系統(tǒng)使大規(guī)模、高精度深井穩(wěn)態(tài)測溫成為可能,光學(xué)掃描法(Optical Scanning Method)在巖石熱導(dǎo)率測試中的廣泛應(yīng)用(Popov et al.,1999),也使巖石熱導(dǎo)率測試工作實(shí)現(xiàn)了批量化、高精度化. 在柴達(dá)木盆地北緣魚卡—馬海尕秀斜坡帶新部署的柴頁1井、紅山構(gòu)造的庫1井等鉆孔,鉆遇了厚層的侏羅系地層,柴達(dá)木盆地東部的歐龍布魯克隆起的柴頁2井、石淺1井,德令哈坳陷帶的埃北1井、德參1井等均鉆遇了石炭系地層,大量來自古生代、中生代地層的巖芯和露頭樣品為研究盆地古生代以來的構(gòu)造-熱演化史提供了先決條件. 筆者對柴達(dá)木盆地新部署的柴頁1井、龍5井、柴頁2井進(jìn)行了系統(tǒng)穩(wěn)態(tài)測溫工作,分析整理了14口鉆孔的靜井溫度數(shù)據(jù),并對盆地鉆井和露頭的巖石熱物性進(jìn)行了測試和研究,建立了柴達(dá)木盆地巖石熱導(dǎo)率柱,揭示了柴達(dá)木盆地現(xiàn)今地溫梯度和大地?zé)崃鞣植继卣?分析了盆地地殼結(jié)構(gòu)對地溫場分布的影響; 以盆地現(xiàn)今地溫場特征為基礎(chǔ),采用裂變徑跡年齡分布特征定性分析與徑跡長度分布數(shù)據(jù)定量模擬相結(jié)合,研究了柴達(dá)木盆地晚古生代以來的沉積埋藏、抬升剝蝕和熱演化史,并結(jié)合區(qū)域構(gòu)造背景,對柴達(dá)木盆地構(gòu)造演化過程進(jìn)行了探討.

    2 地質(zhì)背景

    柴達(dá)木盆地屬于塔里木—中朝板塊分離出的、夾持于秦祁昆地槽褶皺帶之間的微型古陸(王鴻禎,1982,1985). 柴達(dá)木盆地基底總體表現(xiàn)為菱形,為古生代褶皺基底和元古代結(jié)晶基底構(gòu)成的雙重結(jié)構(gòu),盆地莫霍面深度變化介于55~63 km,地殼厚度在盆地中部最大,向盆地邊緣減薄(趙俊猛等,2006).柴達(dá)木盆地的成盆演化經(jīng)歷了晚三疊世末之前多島陸間洋、伸展裂谷洋、殘留海槽階段和晚三疊世末印支運(yùn)動之后的內(nèi)陸盆地形成階段(劉和甫,2001; 鄭孟林等,2004). 盆地南界發(fā)育東昆侖中央斷裂(簡稱昆中斷裂),北界為宗務(wù)隆山斷裂,西界為阿爾金斷裂,發(fā)育了南華系—古生界、中生界和新生界三套沉積地層,中新生界沉積巖最大連續(xù)厚度為17200 m,新生代柴達(dá)木盆地可劃分為6個(gè)一級構(gòu)造單元(圖1),分別為: 三湖坳陷、歐龍布魯克隆起、德令哈拗陷、昆北坳陷、一里坪坳陷、祁南逆沖帶(楊超等,2012).

    圖1 柴達(dá)木盆地構(gòu)造位置與測溫鉆孔分布圖(據(jù)Yin等(2008)修改)Fig.1 The tectonic setting of the Qaidam basin and distribution of the wells (Modified from Yin et al. (2008))

    3 柴達(dá)木盆地現(xiàn)今大地?zé)崃鞣植继卣?/h2>

    大地?zé)崃魇堑厍騼?nèi)部熱作用過程最直接的地表顯示(汪集旸和黃少鵬,1990),它是地溫場最重要的表征,比其他單項(xiàng)地?zé)豳Y料(溫度、地溫梯度)更能反映一個(gè)地區(qū)地?zé)釄龅幕咎攸c(diǎn)(陳墨香,1989).在一維穩(wěn)態(tài)條件下,熱流量(q)是巖石熱導(dǎo)率(k,W/(m·K))與垂直地溫梯度(dT/dZ, ℃·km-1)的乘積,即q=-k(dT/dZ),單位為mW·m-2,負(fù)號表示熱流方向與地溫梯度方向相反,其分布特征與現(xiàn)代地殼運(yùn)動和構(gòu)造活動之間呈相關(guān)性.

    3.1 測溫?cái)?shù)據(jù)

    地層溫度研究的主要依據(jù)是鉆井測溫,旨在借助于測量井液溫度顯示地下巖層的原始溫度,依據(jù)溫度資料類型可分系統(tǒng)穩(wěn)態(tài)測溫、靜井溫度、準(zhǔn)穩(wěn)態(tài)測溫和瞬態(tài)測溫四種(饒松等,2013).

    3.1.1 系統(tǒng)穩(wěn)態(tài)測溫

    鉆井過程中,由于地層打開、泥漿注入等因素的影響,地層原始地溫的分布遭到破壞. 井孔中溫度達(dá)到平衡的狀態(tài)后系統(tǒng)測得的溫度為系統(tǒng)穩(wěn)態(tài)測溫,該測溫資料能客觀地反映地層的真實(shí)溫度(Qiu,2003; 唐曉音等,2014). 2013年至2014年,筆者先后兩次在柴達(dá)木盆地北緣、東部進(jìn)行了鉆井測溫,得到了柴頁1井、柴頁2井、龍5井共3口鉆孔的穩(wěn)態(tài)測溫?cái)?shù)據(jù).這些鉆孔的靜井時(shí)間最短是1個(gè)月,最長的超過4個(gè)月,鉆孔內(nèi)的溫度基本恢復(fù)到了開孔前的初始狀態(tài). 工作使用鉆孔溫度連續(xù)采集系統(tǒng),配置鉑電阻探頭,測量分辨率為0.1 ℃,數(shù)據(jù)記錄間隔為0.01,從孔口開始記錄溫度,測量深度最淺的有875 m,最深的是2164 m(表1).

    表1 測溫井基礎(chǔ)資料與地溫、熱導(dǎo)率測試結(jié)果Table 1 The results of temperature,thermal conductivity and basic data of wells

    續(xù)表1

    續(xù)表1

    注:表中*代表原位熱導(dǎo)率; 編號1—17的鉆孔數(shù)據(jù)為本研究新增,編號18—56 的鉆孔數(shù)據(jù)是基于Qiu(2003)的成果增加經(jīng)緯度坐標(biāo)后匯編.

    圖2 柴達(dá)木盆地新增測溫鉆孔溫度-深度曲線1—柴頁1(Chaiye1); 2—龍5(long5); 3—柴頁2(Chaiye2); 4—建西1(Jianxi1); 5—臺南4(Tainan4); 6—建參1(Jiancan1); 7—澀19(Se19); 8—澀21 (Se21); 9—沙中151 (Shazhong151); 10—臺南3 (Tainan3); 11—臺南5 (Tainan5); 12—沙中30 (Shazhong30); 13—沙中100 (Shazhong100); 14—澀22 (Se22); 15—花14 (Hua14); 16—花20 (Hua20);17—七26 (Qi26) . 其中,編號1—3(紅色曲線)為穩(wěn)態(tài)測溫?cái)?shù)據(jù);編號4—17(黑色曲線)為整理分析后的靜井測溫曲線.Fig.2 The change of geo-temperature with depth in the 17 boreholes of Qaidam basin

    圖3 柴達(dá)木盆地現(xiàn)今地溫梯度分布圖Fig.3 Distribution pattern of present-day geothermal gradient in Qaidam basin

    3.1.2 靜井溫度

    靜井溫度,指關(guān)井?dāng)?shù)天或長期關(guān)井后對不同的含油層段進(jìn)行試油,有時(shí)數(shù)次關(guān)井并將儀器下放至不同深度進(jìn)行溫度測量,從而獲得鉆孔的系統(tǒng)試油溫度,它在一定程度上可以替代研究區(qū)系統(tǒng)穩(wěn)態(tài)測溫,這類資料是地溫場研究的主要依據(jù)之一.

    3.1.3 準(zhǔn)穩(wěn)態(tài)測溫和瞬態(tài)測溫

    準(zhǔn)穩(wěn)態(tài)測溫,是指完井后靜井1~3天內(nèi)所測的溫度數(shù)據(jù),與穩(wěn)定測溫相比,其誤差較大.而在完井后不到一天所進(jìn)行的測溫則叫瞬時(shí)測溫?cái)?shù)據(jù),由于靜井時(shí)間短,地溫與井溫尚未達(dá)到平衡,測得的井溫曲線不能反映真正的地溫情況,油田中鉆孔測溫?cái)?shù)據(jù)一般屬于這類溫度數(shù)據(jù).

    上述各類溫度數(shù)據(jù)的精度各異,系統(tǒng)穩(wěn)態(tài)測溫資料和靜井溫度數(shù)據(jù)比較可靠,是筆者研究柴達(dá)木盆地現(xiàn)今地溫場分布特征的主要依據(jù).

    3.2 柴達(dá)木盆地現(xiàn)今地溫梯度分布特征

    筆者對新增的鉆孔測溫?cái)?shù)據(jù)用最小二乘法進(jìn)行線性擬合,獲取了17個(gè)鉆孔的地溫梯度數(shù)據(jù),并對前人研究的原始測溫?cái)?shù)據(jù)整理分析,共得到了56口鉆孔的地溫梯度數(shù)據(jù)(表1),在此基礎(chǔ)上,繪制了柴達(dá)木盆地現(xiàn)今地溫梯度分布圖(圖3). 研究表明,柴達(dá)木盆地現(xiàn)今地溫梯度為17.1~38.6 ℃·km-1,平均為28.6±4.6 ℃·km-1,較塔里木盆地平均地溫梯度22.6±3.0 ℃·km-1高(馮昌格等,2009),也高于準(zhǔn)噶爾盆地21.3±3.7 ℃·km-1(饒松等,2013),但低于鄂爾多斯盆地29.3 ℃·km-1(饒松等,2013),更低于中國東部中新生代盆地(He et al.,2001,2002; Yang et al.,2004; Li et al.,2013; Yuan et al.,2009; Zuo et al.,2011).柴達(dá)木盆地西部的昆北逆沖帶、一里坪坳陷帶地溫梯度最高,值介于22.0~38.6 ℃·km-1,平均 29.7 ℃·km-1,尤其是東柴山、牛鼻梁及澀北一號、二號穹窿地帶地溫梯度平均可達(dá)34.2 ℃·km-1. 柴北緣的祁南逆沖帶地溫梯度相對低,高的地溫梯度值分布在魚卡—馬海一帶,平均為29.2 ℃·km-1,柴達(dá)木盆地東部的三湖坳陷、德令哈坳陷及歐龍布魯克隆起區(qū)地溫梯度最低,平均值小于24 ℃·km-1.

    3.3 巖石熱導(dǎo)率

    巖石熱導(dǎo)率表示巖石傳熱的特性,是研究一個(gè)區(qū)域大地?zé)崃鳌⑴璧責(zé)嵫莼仨毜臒嵛镄詤?shù),其物理意義為沿?zé)醾鲗?dǎo)方向在單位厚度巖石兩側(cè)的溫度差為1 ℃時(shí)單位時(shí)間內(nèi)所通過的比熱流量,單位為W/(m·K).

    圖4 柴達(dá)木盆地巖石熱導(dǎo)率與深度關(guān)系剖面圖Fig.4 The relationship between thermal conductivity and depth by different lithology in Qaidam basin

    圖5 柴達(dá)木盆地?zé)釋?dǎo)率統(tǒng)計(jì)直方圖Fig.5 Histogram of thermal conductivity in Qaidam basin

    3.3.1 樣品采集和分析測試

    筆者對柴達(dá)木盆地臺南3、柴頁1、柴頁2等10個(gè)鉆孔,石灰溝、阿勒格爾泰山、石拐子溝等6條露頭剖面共計(jì)151塊樣品(圖4)進(jìn)行了巖石熱導(dǎo)率測試. 測試采用的是光學(xué)掃描技術(shù)(Popov et al.,1999),測試儀器為德國生產(chǎn)的TCS(Thermal Conductivity Scanning)熱導(dǎo)率自動掃描儀,其測量范圍為0.2~25 W/(m·K),測量精度為±3%.所有樣品都進(jìn)行了飽水校正,以消除樣品成分、結(jié)構(gòu)等因素對巖石的熱導(dǎo)率影響(沈顯杰等,1994).

    3.3.2 測試結(jié)果和柴達(dá)木盆地?zé)釋?dǎo)率柱

    結(jié)果顯示柴達(dá)木盆地的巖石熱導(dǎo)率介于0.601~5.520 W/(m·K)之間,平均為2.208 W/(m·K),主要集中在1.820~3.020 W/(m·K)之間(圖4),占測試樣品總數(shù)的73.8%(圖5),與中國西部的塔里木盆地、準(zhǔn)格爾盆地相當(dāng)(馮昌格等,2009; 饒松等,2013). 柴達(dá)木盆地不同巖類熱導(dǎo)率統(tǒng)計(jì)分析(表2)表明,凝灰?guī)r、白云巖、灰?guī)r、礫巖的熱導(dǎo)率最大,粉沙巖、砂巖熱導(dǎo)率次之,泥巖、砂質(zhì)泥巖、油頁巖最小,且?guī)r石熱導(dǎo)率總體隨著埋深的增大,地層的變老而增大. 第四系和侏羅系的巖石熱導(dǎo)率均較小,歸因于樣品埋藏較淺,比較松散,且樣品中泥巖、油頁巖和碳質(zhì)泥巖占了相當(dāng)大的比例,這也證實(shí)了巖石的結(jié)構(gòu)、巖性對熱導(dǎo)率的控制作用.

    研究發(fā)現(xiàn)熱導(dǎo)率的值與砂泥含量比成正相關(guān). 統(tǒng)計(jì)數(shù)據(jù)顯示砂質(zhì)泥巖熱導(dǎo)率的平均值與隨深度變化的斜率均比泥巖大,砂巖熱導(dǎo)率的平均值與隨深度變化的斜率均比粉砂巖大(圖5),從泥巖、砂質(zhì)泥巖、泥質(zhì)粉砂巖、砂巖的過渡中,熱導(dǎo)率值隨砂泥比增加相應(yīng)變大.

    3.4 柴達(dá)木盆地現(xiàn)今大地?zé)崃鞣植继卣?/p>

    依據(jù)本研究所計(jì)算地溫梯度和實(shí)測的巖石熱導(dǎo)率數(shù)據(jù),得到了17個(gè)高質(zhì)量的大地?zé)崃髦?并整理分析前人的原始數(shù)據(jù),共得到了56個(gè)盆地的大地?zé)崃髦?表1),在此基礎(chǔ)上編制了柴達(dá)木盆地現(xiàn)今大地?zé)崃鞣植紙D(圖6).

    柴達(dá)木盆地現(xiàn)今大地?zé)崃髦到橛?2.9~70.4 mW·m-2, 平均 55.1±7.9 mW·m-2,高于我國西部的塔里木盆地、準(zhǔn)噶爾盆地(馮昌格等,2009; 饒松等,2013),低于中國東部及海域沉積盆地(He et al.,2001,2002; Yang et al.,2004; Li et al.,2013; Yuan et al.,2009; Zuo et al.,2011),也低于我國大陸地區(qū)大地?zé)崃髌骄?3.0±24.2 mW·m-2(Hu et al.,2000). 柴達(dá)木盆地現(xiàn)今大地?zé)崃鞣植继卣髋c地溫梯度分布相近,一里坪坳陷、昆北逆沖帶等地區(qū)的熱流值最高,介于56.0~70.4 mW·m-2,平均59.1 mW·m-2. 祁南逆沖帶內(nèi)熱流分布不均勻,熱流值沿SSE方向增大,冷湖附近熱流最低,平均46 mW·m-2,魚卡—馬海尕秀一帶熱流最高,平均值為55 mW·m-2. 三湖坳陷、德令哈坳陷及歐龍布魯克隆起帶熱流普遍偏低,平均值小于45 mW·m-2.柴達(dá)木盆地現(xiàn)今大地?zé)崃鞣植际遣痪鶆虻?盆地現(xiàn)今熱流分布狀況是構(gòu)造演化的結(jié)果和響應(yīng),從而反映了柴達(dá)木盆地在晚期發(fā)育階段經(jīng)歷了復(fù)雜的構(gòu)造運(yùn)動和深部熱活動.

    表2 柴達(dá)木盆地巖石熱導(dǎo)率柱Table 2 Thermal conductivity data of different formation in Qaidam basin

    柴達(dá)木盆地現(xiàn)今地溫場分布特征表明: 柴西昆北逆沖帶、一里坪坳陷無論是地溫梯度還是大地?zé)崃骶哂谂璧仄渌貐^(qū),屬于“高溫區(qū)”; 柴北緣祁南逆沖帶次之,屬于“中溫區(qū)”; 柴東三湖坳陷、德令哈坳陷及歐龍布魯克隆起區(qū)最低,屬于“低溫區(qū)”.

    4 柴達(dá)木盆地晚古生代以來構(gòu)造-熱演化史

    盆地構(gòu)造-熱演化史是盆地動力研究的主要內(nèi)容. 盆地的沉積埋藏、抬升剝蝕和熱演化過程是研究的核心內(nèi)容. 裂變徑跡定年是近年來地質(zhì)年代學(xué)研究的重要領(lǐng)域,廣泛地應(yīng)用于沉積盆地物源、構(gòu)造運(yùn)動時(shí)間、剝蝕量恢復(fù)、盆地古地溫演化及斷層活動性研究中(趙紅格等,2003; Armstrong,2005; 余心起等,2013; 高成等,2014; 陳漢林等,2014; 李振華等,2014; 王先美等,2008). 基于此,筆者利用磷灰石、鋯石裂變徑跡數(shù)據(jù)研究了柴達(dá)木盆地晚古生代以來的構(gòu)造-熱演化史.

    4.1 裂變徑跡定年原理

    裂變徑跡是指礦物內(nèi)238U原子核自發(fā)裂變碎片形成的輻射損傷痕跡,其密度和長度受時(shí)間和溫度的控制(Gleadow et al.,1983). 裂變徑跡定年遵從同位素定年的基本原理,但它測量的是放射性衰變對礦物晶體的物理損傷,而不是另一種同位素,該定年體系是確定巖石低溫?zé)釟v史的一項(xiàng)技術(shù),測試的目標(biāo)礦物主要是磷灰石和鋯石,其中,磷灰石封閉溫度100~120 ℃(Green et al.,1989),鋯石185~240 ℃(Brandon et al.,1998; Tagami and O′Sullivan,2005).

    4.2 樣品采集和實(shí)驗(yàn)方法

    7件磷灰石樣品和12件鋯石樣品采自鉆孔巖芯和剖面露頭. 巖芯為柴東上石炭統(tǒng)克魯克組,露頭樣品采自柴北緣聯(lián)合溝剖面的上奧陶統(tǒng)灘間山群、小賽什騰山剖面的上奧陶統(tǒng)灘間山群及柴西阿爾金南坡的上奧陶統(tǒng)灘間山群(圖1). 樣品巖性主要是玄武巖、砂巖、輝長巖,裂變徑跡樣品的制備與分析在中國地質(zhì)大學(xué)地質(zhì)過程與礦產(chǎn)資源重點(diǎn)實(shí)驗(yàn)室裂變徑跡實(shí)驗(yàn)室完成,測試方法詳見文獻(xiàn)(Gallagher et al.,1998).

    磷灰石樣品,是將采集的樣品利用傳統(tǒng)方法分離出磷灰石單礦物,用環(huán)氧樹脂粘在光薄片上,研磨拋光出礦物最大內(nèi)表面. 在室溫20 ℃左右下,用5 N的HNO3蝕刻18~22 s揭示自發(fā)徑跡. 將低鈾白云母外探測器與磷灰石一同封裝好送入反應(yīng)堆輻照,輻照在中國原子能科學(xué)研究院進(jìn)行,之后在20 ℃下用40%HF蝕刻18 min揭示誘發(fā)徑跡,中子注量利用CN5鈾玻璃標(biāo)定(Bellemans et al., 1995). 測定年齡時(shí)每個(gè)樣品任意選取20個(gè)左右質(zhì)量好、平行C軸的顆粒進(jìn)行測年(假設(shè)有足夠的顆粒). 根據(jù)IUGS(International Union of Geological Sciences,國際地質(zhì)科學(xué)聯(lián)合會)推薦的ζ常數(shù)法和標(biāo)準(zhǔn)裂變徑跡年齡方程(Hurford and Green,1982,1983)計(jì)算年齡值. 利用Durango和Fishcanyon磷灰石標(biāo)樣標(biāo)定的Zeta(ζ)參數(shù)值為353±18. 封閉徑跡的長度在偏光顯微鏡下進(jìn)行測量,徑跡長度分析中,對每個(gè)樣品測量100條水平(如果有100條的話)、狹窄的封閉徑跡(Gleadow et al.,1986).

    圖6 柴達(dá)木盆地現(xiàn)今熱流分布圖Fig.6 Distribution pattern of present-day heat flow in Qaidam basin

    鋯石樣品,是將采集的樣品利用傳統(tǒng)方法分離出鋯石單礦物,采用外探測器法將鋯石樣品顆粒用Teflon進(jìn)行固定,經(jīng)磨平和拋光后制成薄片,放入KOH-NaOH共溶混合物中進(jìn)行蝕刻,蝕刻溫度228 ℃,時(shí)間20~26 h; 與低鈾白云母外探測器貼緊后,送中國原子能研究所核反應(yīng)堆反射孔道照射. 反應(yīng)堆輻射通量為1×1015cm-2.低鈾白云母外探測器蝕刻條件為40%HF溶液,室溫,18 min. Zeta(ζ)標(biāo)樣選用美國國家標(biāo)準(zhǔn)局CN2鈾標(biāo)準(zhǔn)玻璃,Zeta(ζ)值為101.2±3.2. 在偏光顯微鏡下進(jìn)行裂變徑跡測量,運(yùn)用Zeta常數(shù)校準(zhǔn)法計(jì)算鋯石裂變徑跡年齡.

    4.3 測試結(jié)果

    磷灰石裂變徑跡年齡與長度測試結(jié)果見表3. 柴東上石炭統(tǒng)克魯克組三個(gè)樣品(Zk3-1-19,Zk3-1-74,Zk3-1-81),AFT年齡集中在61—72 Ma,小于其相應(yīng)的地層年齡;平均封閉徑跡長度為11.8~12.7 μm,小于初始徑跡長度(16.3±0.9 μm),呈單峰態(tài)分布,顯然這三件樣品經(jīng)歷了中度退火行為,最高溫度超過120 ℃. 柴北緣聯(lián)合溝上奧陶統(tǒng)灘間山群的4個(gè)AFT樣品(C1203-10,C1203-12,C1203-18,C1203-39),AFT年齡集中在39—58 Ma,遠(yuǎn)小于其相應(yīng)的地層年齡,平均封閉徑跡長度為10.3~11.4 μm,小于初始徑跡長度(16.3±0.9 μm),C1203-12,C1203-18呈似雙峰態(tài)分布. 磷灰石裂變徑跡年齡的2均大于或遠(yuǎn)遠(yuǎn)大于5%,表明所分析樣品的單顆粒年齡差異屬于統(tǒng)計(jì)誤差,并且樣品的磷灰石裂變徑跡年齡為單一的年齡,沒有多組年齡現(xiàn)象(Galbraith,1984).

    鋯石裂變徑跡年齡測試結(jié)果詳見表4. 測試結(jié)果表明,6個(gè)樣品P(2)>5%,其中,2個(gè)采自柴東上石炭統(tǒng)的樣品(zk3-1-74,zk3-1-81),ZFT年齡集中在155—195 Ma,小于其相應(yīng)的地層年齡,3個(gè)為柴北緣上奧陶統(tǒng)的樣品(C1204-01,C1203-10,DXS17001),ZFT年齡集中在143—162 Ma,遠(yuǎn)小于其相應(yīng)的地層年齡; 1個(gè)采自柴西阿爾金南坡上奧陶統(tǒng)的樣品(DMY15801),ZFT年齡集中在177 Ma,小于其相應(yīng)的地層年齡.未通過2驗(yàn)證(P(2)<5%)的樣品,用BinomFit軟件(Brandon,2002)進(jìn)行單顆粒裂變徑跡年齡數(shù)據(jù)的分解和頻度統(tǒng)計(jì)分析(圖7). 研究發(fā)現(xiàn),柴達(dá)木盆地東部上石炭統(tǒng)樣品Zk3-1-19存在254.0 Ma和144.0 Ma兩個(gè)年齡分組,樣品Zk3-2-84年齡集中在249.6 Ma和148.6 Ma兩個(gè)頻段; 柴北緣上奧陶統(tǒng)樣品C1203-02年齡集中在249.6 Ma和148.6 Ma兩個(gè)頻段,樣品C1203-18年齡集中在166.7 Ma和124.4 Ma兩個(gè)頻段,樣品C1203-58年齡集中在220.0 Ma和151.2 Ma兩個(gè)頻段; 柴西阿爾金南坡上奧陶統(tǒng)樣品DMY15201存在101.3 Ma和77.8 Ma兩個(gè)年齡分組. 綜合12個(gè)樣品測試結(jié)果,分析得出鋯石裂變徑跡年齡主要分布在254.0—199 Ma、177—148.6 Ma兩個(gè)時(shí)間段.

    表3 柴達(dá)木盆地磷灰石裂變徑跡分析數(shù)據(jù)Table 3 Apatite fission track data of the sample from Qaidam basin

    注:n為統(tǒng)計(jì)的磷灰石顆粒數(shù),ρd,ρs,ρi分別表示標(biāo)準(zhǔn)鈾玻璃跡密度、樣品自發(fā)裂變徑跡密度和誘發(fā)裂變徑跡密度;P(2)為2檢驗(yàn)的概率,當(dāng)P(2)>5%時(shí),通常認(rèn)為所測單顆粒年齡屬于同組年齡,否則, 屬于不同年齡組(Green, 1981);L為樣品平均封閉徑跡長度;t±1σ為樣品中值裂變徑跡年齡, 由Zeta (Hurford, 1990)法計(jì)算得出:ζ=353±18.

    表4 柴達(dá)木盆地鋯石裂變徑跡分析數(shù)據(jù)Table 4 Zircon fission track data of the sample from Qaidam basin

    注:n為統(tǒng)計(jì)的鋯石顆粒數(shù),ρd,ρs,ρi分別表示標(biāo)準(zhǔn)鈾玻璃跡密度、樣品自發(fā)裂變徑跡密度和誘發(fā)裂變徑跡密度;P(2)為2檢驗(yàn)的概率, 當(dāng)P(2)>5%時(shí),通常認(rèn)為所測單顆粒年齡屬于同組年齡, 否則, 屬于不同年齡組(Green, 1981);t±1σ為樣品中值裂變徑跡年齡, 由Zeta (Hurford, 1990)法計(jì)算得出:ζ=101.2±3.2.

    圖7 鋯石裂變徑跡年齡組分齡概率分布圖Fig.7 Probability density plots with best fit for the peak values of the zircon fission track age

    4.4 熱演化史模擬

    為了深入研究柴達(dá)木盆地構(gòu)造-熱演化過程,利用Ketcham等(1999)單組分退火模型和Monte Carlo法進(jìn)行磷灰石裂變徑跡樣品的時(shí)間-溫度歷史模擬,模擬軟件采用Hefty(Ketcham,2005). 磷灰石樣品封閉徑跡長度值和中值年齡值是模擬的基本參數(shù),動力學(xué)參數(shù)D值選取初始值1.65,初始徑跡長度采用16.3 μm.

    裂變徑跡熱模擬必須以區(qū)域地質(zhì)背景為基礎(chǔ),充分考慮區(qū)域及周緣的構(gòu)造背景,并以此建立模擬的邊界條件. 本區(qū)熱史模擬的邊界條件包括: (1)地表溫度: 柴達(dá)木盆地常年平均地表溫度為10 ℃左右,考慮到溫度的季節(jié)變化,本文取地表溫度范圍為5~30 ℃; (2) 區(qū)域構(gòu)造事件: ① 岡底斯—念青唐古拉地塊群與歐亞大陸拼合.早白堊世,岡底斯—念青唐古拉地塊群與歐亞大陸拼合,碰撞作用導(dǎo)致南祁連山和阿爾金山發(fā)生隆升(葛肖虹等,2001),柴達(dá)木盆地遭受擠壓. ② 印度板塊和歐亞板塊碰撞.自始新世以來,印度板塊與歐亞板塊一直處于擠壓與碰撞狀態(tài),導(dǎo)致了青藏高原及周緣包括柴達(dá)木盆地的隆升(滕吉文等,1999; 徐志琴等,2006). (3) 磷灰石裂變徑跡表觀年齡:測試結(jié)果表明,樣品表觀年齡集中在75—58 Ma和41—39 Ma兩個(gè)階段. (4) 模擬起算時(shí)間:筆者采用Corrigan(1991) 的估算公式t=[lo/(l-σ)]tFTA,其中l(wèi)o是新鮮徑跡長度,l是觀測到的平均徑跡長度,tFTA是觀測到的徑跡年齡,σ是觀測到的長度標(biāo)準(zhǔn)方差. 經(jīng)計(jì)算,本研究熱史模擬起算時(shí)間約在80 Ma. Ketcham(2005)根據(jù)實(shí)測與模擬的徑跡長度和年齡檢驗(yàn)值GOF(Goodness of Fit)將模擬結(jié)果分為3部分: 可以接受的熱史曲線范圍、較好的熱史曲線范圍和最佳熱歷史擬合曲線. 若年齡與長度GOF檢驗(yàn)都大于0.05,表明模擬結(jié)果是可以接受的,當(dāng)它們超過0.5時(shí),模擬結(jié)果是較好的. 本次熱歷史模擬樣品的長度分布檢驗(yàn)值為0.77~0.97,年齡檢驗(yàn)值(GOF)為0.94~0.99,均大于0.5(圖8),說明熱史模擬結(jié)果都是高質(zhì)量的.

    模擬結(jié)果(圖8)表明,所有溫度-時(shí)間曲線形態(tài)大體接近,反映了樣品白堊紀(jì)末以來的熱歷史. 柴達(dá)木盆地東部的樣品~80Ma冷卻到部分退火帶(60~120 ℃)之上,即在80 Ma之前,樣品發(fā)生了冷卻,此后樣品經(jīng)歷了約40 Ma緩慢的增溫過程; 樣品在~33.6 Ma,~7.1 Ma,~1.8 Ma三個(gè)階段經(jīng)歷了冷卻,冷卻速率分別為2.9 ℃/Ma、12.3 ℃/Ma、11.5 ℃/Ma. 柴達(dá)木盆地北緣的樣品約87 Ma冷卻到部分退火帶之上,冷卻速率~2.1 ℃/Ma,此后樣品處于緩慢增溫狀態(tài); ~41.1 Ma,~9.6 Ma,~2.9 Ma樣品發(fā)生了冷卻,冷卻速率分別為0.8 ℃/Ma、3.4 ℃/Ma、6.8 ℃/Ma. 總的來看,柴達(dá)木盆地磷灰石裂變徑跡樣品在87—62 Ma,41.1—33.6 Ma, 9.6—7.1 Ma,2.9—1.8 Ma四個(gè)時(shí)期發(fā)生了冷卻,同時(shí)還發(fā)現(xiàn),柴北緣相比柴東地區(qū)樣品發(fā)生冷卻的時(shí)間要早,推測青藏高原隆升對盆地構(gòu)造演化的作用是沿SW—NE向逐步推進(jìn)的過程.

    5 討論

    5.1 地殼結(jié)構(gòu)對柴達(dá)木盆地現(xiàn)今地溫場分布的控制作用

    前人已對柴達(dá)木盆地的地殼結(jié)構(gòu)和性質(zhì)做了系統(tǒng)研究(滕吉文,1974; Zhu and Helmberger,1998; 趙國澤等,2004; 李永華等,2006; 趙俊猛等,2006),得出柴達(dá)木盆地莫霍面深度變化介于55~63 km,地殼厚度在盆地中部(如達(dá)布遜湖)最大約63 km,向盆地邊緣減薄,盆地西部、西南緣地殼厚度~50 km,魚卡、大柴達(dá)木附近地殼亦較薄,厚度為~50 km(趙俊猛等,2006).筆者研究發(fā)現(xiàn)柴達(dá)木盆地相比我國中、東部其他盆地(如松遼盆地)莫霍面埋藏較深(崔軍平等,2007),地殼厚度大,消弱地幔的熱能通過地殼巖石向地面?zhèn)鲗?dǎo),導(dǎo)致整個(gè)盆地地溫梯度相對低,屬于溫盆類型.昆北逆沖帶、一里坪坳陷、柴北緣祁南逆沖帶地殼厚度較薄,區(qū)域莫霍面埋藏淺,蓋層受到地幔熱影響較其他區(qū)域大,屬于“中—高溫區(qū)”.總體講,柴達(dá)木盆地地溫梯度分布特征與地殼厚度成鏡像關(guān)系,地殼厚度小的地區(qū)地溫高,地殼厚的地區(qū)地溫低,盆地現(xiàn)今地溫場分布特征受控于深部地殼結(jié)構(gòu)和性質(zhì).

    5.2 柴達(dá)木盆地晚古生代以來構(gòu)造-熱演化

    研究表明,柴達(dá)木盆地晚古生代以來在254.0—199 Ma,177—148.6 Ma,87—62 Ma,41.1—33.6 Ma,9.6—7.1 Ma,2.9—1.8 Ma 六個(gè)時(shí)期經(jīng)歷了顯著的構(gòu)造事件.

    254.0—199 Ma,該裂變徑跡年齡組記錄了三疊紀(jì)中—晚期發(fā)生的構(gòu)造事件(印支運(yùn)動).該構(gòu)造事件揭開了柴達(dá)木盆地中、新生代構(gòu)造演化的序幕,也導(dǎo)致了阿爾金山的隆升(任收麥等,2009),盆地西部紅柳溝一帶、盆地北緣冷湖、結(jié)綠素一帶、盆地東部大煤溝普遍發(fā)現(xiàn)侏羅系底部不整合.177—148.6 Ma,是盆地侏羅紀(jì)中—晚期發(fā)生的構(gòu)造運(yùn)動(燕山運(yùn)動早期),這次構(gòu)造事件的發(fā)生與班公湖—怒江侏羅紀(jì)洋盆的完全閉合和岡底斯塊體與羌塘塊體的最終碰撞有關(guān)(滕吉文等,1999).在柴達(dá)木盆地中表現(xiàn)為中、下侏羅統(tǒng)與上侏羅統(tǒng)(J1+2與J3-K)或中、下侏羅統(tǒng)與白堊系(J1+2與K)之間的不整合現(xiàn)象,同時(shí)該構(gòu)造事件使區(qū)域的古氣候、古地理環(huán)境發(fā)生了根本性的變化,即由早期的溫暖-潮濕氣候轉(zhuǎn)變成為后期的干燥氧化氣候條件(湯良杰等,2000).87—62 Ma,該年齡組記錄了盆地白堊紀(jì)末期—古近紀(jì)早期的構(gòu)造事件,該構(gòu)造事件與早白堊世末岡底斯—念青唐古拉地塊群與歐亞大陸拼合有關(guān)(徐志琴等,2006).裂變徑跡熱模擬表明,該構(gòu)造運(yùn)動導(dǎo)致柴達(dá)木盆地東部隆升并遭受剝蝕,推測該構(gòu)造時(shí)期形成現(xiàn)今歐龍布魯克隆起帶的雛形.相比之下,柴北緣經(jīng)受該構(gòu)造運(yùn)動的影響相對強(qiáng),剝露速率~50 m/Ma,至始新世中期,柴北緣處于緩慢沉降狀態(tài),在弱擠壓環(huán)境下形成坳陷盆地,這與周建勛等(2003)通過平衡剖面所獲得的結(jié)論相同.

    圖8 柴達(dá)木盆地樣品AFT熱史反應(yīng)結(jié)果1 可接受的熱史曲線,即GOF在0.05~0.5之間;2 較好的熱史曲線,即GOF大于0.5;3 最好的熱史曲線.Fig.8 Inversed results of thermal history by AFT from Qadiam basin 1 Acceptable paths with GOF between 0.05 and 0.5; 2 Good paths with GOF more than 0.5; 3 Best fit.

    41.1—33.6 Ma,該期構(gòu)造事件(喜馬拉雅運(yùn)動早期)與印度板塊和歐亞板塊碰撞有關(guān)(徐志琴等,2006),在盆地中表現(xiàn)為下干柴溝組與下伏地層(路樂河組)之間的不整合,下油砂山組與上干柴溝組之間在盆地邊緣的不整合現(xiàn)象.王非等(2001)對南祁連南部中晚三疊世花崗質(zhì)巖體中的鉀長石、黑云母、白云母進(jìn)行40Ar/39Ar定年時(shí),分析認(rèn)為巖體約30 Ma經(jīng)歷一次快速的冷卻事件,也佐證了筆者研究結(jié)果.9.6—7.1 Ma,這組裂變徑跡年齡記錄了印度板塊與歐亞大陸碰撞作用在中新世末的進(jìn)一步增強(qiáng),促使了青藏高原的強(qiáng)烈隆升(潘裕生,1999),碰撞作用的遠(yuǎn)程傳導(dǎo)誘發(fā)了青藏高原東北緣構(gòu)造變形.裂變徑跡模擬結(jié)果表明,該時(shí)期盆地快速隆升,并遭受剝蝕,剝露速率~300 m/Ma,盆地周緣發(fā)生顯著構(gòu)造變形.張培震等(2006)以六盤山、積石山及其相鄰盆地為對象,開展了青藏高原東北緣晚新生代擴(kuò)展與隆升研究,發(fā)現(xiàn)這些地區(qū)也在該時(shí)期(5—10 Ma,或約8 Ma)發(fā)生了明顯的構(gòu)造變形.2.9—1.8 Ma,該組年齡反映了上新世與第四紀(jì)之間一次強(qiáng)烈的構(gòu)造運(yùn)動(喜馬拉雅運(yùn)動晚期),熱史模擬表明,該期構(gòu)造活動使盆地遭受強(qiáng)烈擠壓、隆升,從而導(dǎo)致了先存斷裂普遍遭受強(qiáng)烈改造,先期褶皺得到進(jìn)一步發(fā)展,同時(shí)也使得新地層卷入褶皺.新生代地層中自生礦物的O和C同位素測試結(jié)果,以及酒西盆地沉積與青藏高原隆升響應(yīng)關(guān)系的研究均顯示了該期強(qiáng)烈的構(gòu)造作用(湯濟(jì)廣,2007).

    6 結(jié)論

    (1) 柴達(dá)木盆地現(xiàn)今地溫梯度介于17.1~38.6 ℃·km-1,平均為28.6±4.6 ℃·km-1,現(xiàn)今大地?zé)崃髦到橛?2.9~70.4 mW·m-2,平均 55.1±7.9 mW·m-2,柴達(dá)木盆地總體屬于溫盆.不同構(gòu)造單元地溫場分布存在差異,昆北逆沖帶、一里坪坳陷屬于“高溫區(qū)”,祁南逆沖帶屬于“中溫區(qū)”,三湖坳陷、德令哈坳陷及歐龍布魯克隆起區(qū)屬于“低溫區(qū)”.現(xiàn)今熱流分布反映了柴達(dá)木盆地在晚期發(fā)育階段經(jīng)歷了復(fù)雜的構(gòu)造事件和深部熱背景,盆地現(xiàn)今地溫場分布特征受控于地殼深部結(jié)構(gòu)、盆地構(gòu)造等因素.

    (2) 柴達(dá)木盆地晚古生代以來經(jīng)歷了六期顯著的構(gòu)造事件.Ⅰ(254.0—199 Ma) 該構(gòu)造事件揭開了柴達(dá)木盆地中、新生代構(gòu)造演化的序幕; Ⅱ(177—148.6 Ma) 區(qū)域由早期的溫暖-潮濕氣候轉(zhuǎn)變成為后期的干燥氧化氣候; Ⅲ(87—62 Ma) 柴達(dá)木盆地東部緩慢隆升遭受剝蝕,歐龍布魯克隆起帶形成雛形,柴北緣在弱擠壓環(huán)境下形成坳陷盆地; Ⅳ(41.1—33.6 Ma) 是盆地對印度板塊和歐亞板塊碰撞早期的響應(yīng); Ⅴ(9.6—7.1 Ma) 青藏高原的強(qiáng)烈隆升誘發(fā)了柴達(dá)木盆地快速隆升,并遭受剝蝕,盆地周緣變形強(qiáng)烈; Ⅵ(2.9—1.8 Ma) 盆地遭受強(qiáng)烈擠壓、隆升與構(gòu)造變形,導(dǎo)致了先存斷裂普遍遭受強(qiáng)烈改造,先期褶皺得到進(jìn)一步發(fā)展.總的來說,柴達(dá)木盆地成盆早,但構(gòu)造定型晚,成盆過程相對寧靜,但晚期改造強(qiáng)烈,控制盆地構(gòu)造變形的主要喜山期構(gòu)造運(yùn)動中(9.6—7.1 Ma)、晚期(2.9—1.8 Ma),同時(shí)青藏高原隆升作用對其東北緣的傳遞、擴(kuò)展及地球動力學(xué)作用仍然有待于進(jìn)一步研究.

    致謝 兩位審稿專家提出了寶貴的修改意見,中國石油大學(xué)(北京)邱楠生教授無私提供了大量原始測溫、熱物性及古地溫?cái)?shù)據(jù),長江大學(xué)饒松老師在鉆井穩(wěn)態(tài)測溫中給予了很大幫助; 地溫資料收集得到中國石油青海油田勘探開發(fā)研究院的大力支持,在此一并致謝!

    Armstrong P A. 2005. Thermochronometers in sedimentary basins.ReviewsinMineralogy&Geochemistry,58(1): 499-525.

    Bellemans F,De Corte F,Van Den Haute P. 1995. Composition of SRM and CN U-doped glasses: Significance for their use as thermal neutron fluence monitors in fission track dating.RadiationMeasurements,24(2): 153-160.

    Brandon M T, Roden-Tice M K, Garver J I. 1998. Late Cenozoic exhumation of the Cascadia accretionary wedge in the Olympic Mountains, northwest Washington State.Geol.Soc.Am.Bull., 110(8): 985-1009.

    Brandon M T. 2002. Decomposition of mixed grain age distributions using binomfit.OnTrack, 24: 13-18.

    Chen H L, Chen S Q, Lin X B. 2014. A review of the Cenozoic tectonic evolution of Pamir syntax.AdvancesinEarthScience(in Chinese), 29(8): 890-902.

    Chen M X. 1989. Method for determining terrestrial heat flow in Meso-Cenozoic sedimentary basins.ChineseJournalofGeology(in Chinese), (2): 151-161. Corrigan J. 1991. Inversion of apatite fission track data for thermal history information.JournalofGeophysicalResearch, 96(B6): 10347-10360. Cui J P, Ren Z L, Xiao H, et al. 2007. Study on temperature distribution and controlling factors in the Hai Lar Basin.ChineseJournalofGeology(in Chinese), 42(4): 656-665.

    Feng C G, Liu S W, Wang L S, et al. 2009. Present-day geothermal regime in Tarim basin, northwest China.ChineseJournalofGeophysics(in Chinese), 52(11): 2752-2762, doi: 10.3969/j.issn.0001-5733.2009.11.010.

    Galbraith R F. 1984. On statistical estimation in fission track dating.JournaloftheInternationalAssociationforMathematicalGeology, 16(7): 653-669. Gallagher K, Brown R, Johnson C. 1998. Fission track analysis and its application to geological problems.AnnualReviewofEarthandPlanetarySciences, 26: 519-572.

    Gao C, Liu J, Li D W, et al. 2014. The constraints of the fission track thermochronology on active time of the Southern Tibetan Detachment System in Cho Oyu, Tibet.EarthScienceFrontiers(in Chinese), 21(6): 372-380.

    Gao J P, Fang X M, Song C H, et al. 2011. Tectonic-Thermo events of Northern Tibetan Plateau: Evidence from detrital apatite fission track data in Western Qaidam Basin.JournalofJilinUniversity(EarthScienceEdition) (in Chinese), 41(5): 1466-1475.

    Ge X H, Liu Y J, Ren S M, et al. 2001. Re-understanding on some academic problems of the Altun fault.ChineseJournalofGeology(in Chinese), 36(3): 319-325.

    Gleadow A J W, Duddy I R, Lovering J F. 1983. Fission track analysis: A new tool for the evaluation of thermal histories and hydrocarbon potential.APEAJ., 23: 93-102.

    Gleadow A J W, Duddy I R, Green P F, et al. 1986. Confined fission track lengths in apatite: A diagnostic tool for thermal history analysis.Contrib.Mineral.Petrol., 94(4): 405-415.

    Green P F. 1981. A new look at statistics in fission-track dating.Nucl.Tracks, 5(1-2): 77-86.

    Green P F, Duddy I R, Laslett G M, et al. 1989. Thermal annealing of fission tracks in apatite 4. Quantitative modelling techniques and extension to geological timescales.Chem.Geol., 79(2): 155-182.

    He L J, Wang K L, Xiong L P, et al. 2001. Heat flow and thermal history of the South China Sea.PhysicsoftheEarthandPlanetaryInteriors, 126(3-4): 211-220.

    He L J, Xiong L P, Wang K L. 2002. Heat flow and thermal modeling of the Yinggehai Basin, South China Sea.Tectonophysics, 351(3): 245-253. Hu S B, He L J, Wang J Y. 2000. Heat flow in the continental area of China: a new data set.EarthPlanet.Sci.Lett., 179(2): 407-419.

    Hurford A J, Green P F. 1982. A user′s guide to fission track dating calibration.EarthPlanet.Sci.Lett., 59(2): 343-354.

    Hurford A J, Green P F. 1983. The zeta age calibration of fission-track dating.ChemicalGeology, 41: 285-317.

    Hurford A J. 1990. Standardization of fission track dating calibration: Recommendation by the Fission Track Working Group of the I. U. G. S. Subcommission on Geochronology.ChemicalGeology:IsotopeGeoscienceSection, 80(2): 171-178.

    Ketcham R A, Donelick R A, Carlson W D. 1999. Variability of apatite fission track annealing Kinetics Ⅲ: Extrapolation to geological time scales.AmericanMineralogist, 84: 1235-1255. Ketcham R A. 2005. Forward and inverse modeling of low-temperature thermochronometry data.ReviewsinMineralogyandGeochemistry, 58(1): 275-314.

    Li G H. 1992. Characteristics of heat fluid and thermal structure of the crust in the Qaidam Basin [Ph. D. thesis] (in Chinese). Beijing: Institute of Geology, Chinese Academy of Sciences.

    Li Y H, Wu Q J, An Z H, et al. 2006. The Poisson ratio and crustal structure across the NE Tibetan Plateau determined from receiver functions.ChineseJournalofGeophysics(in Chinese), 49(5): 1359-1368.

    Li Z H, Chen G, Ding C, et al. 2014. Paleogeothermal Gradient recovery in middle of Yanshanian in northern Junggar Basin.GeologicalScienceandTechnologyInformation(in Chinese), 33(4): 31-36.

    Li Z X, Xu M, Zhao P, et al. 2013. Geothermal regime and hydrocarbon kitchen evolution in the Jianghan Basin.ScienceChinaEarthSciences, 56(2): 240-257.

    Liu H P. 2001. Geodynamic scenario of coupled basin and mountain system.EarthScience—JournalofChinaUniversityofGeosciences(in Chinese), 26(6): 581-596.

    Pan Y S. 1999. Formation and uplifting of the Qinghai-Tibet plateau.EarthScienceFrontiers(in Chinese), 6(3): 153-163.Popov Y A, Pevzner S L, Pimenov V P, et al. 1999. New geothermal data from the Kola superdeep well SG-3.Tectonophysics, 306(3-4): 345-366.

    Qiu N S. 2001. Research on heat flow and temperature distribution of the Qaidam Basin.JournalofChinaUniversityofMining&Technology(in Chinese), 30(4): 412-415.

    Qiu N S. 2002. Tectono-thermal evolution of the Qaidam Basin, China: evidence from Ro and apatite fission track data.PetroleumGeoscience, 8(3): 279-285.

    Qiu N S. 2003. Geothermal regime in the Qaidam basin, northeast Qinghai-Tibet Plateau.Geol.Mag., 140(6): 707-719.

    Rao S, Hu S B, Zhu C Q, et al. 2013. The characteristics of heat flow and lithospheric thermal structure in Junggar basin, northwest China.ChineseJournalofGeophysics(in Chinese), 56(8): 2760-2770, doi: 10.6038/cjg20130824.Ren S M, Ge X H, Liu Y J, et al. 2009. Multi-stage strike-slip motion and uplift along the Altyn Tagh fault since the Late Cretaceous.GeologicalBulletinofChina(in Chinese), 23(9-10): 926-932.

    Ren Z L. 1993. Geothermal evolution history of the Qaidam Basin: evidence from fluid inclusions and Ro. ∥ Zhao C Y, Liu C Y et al. eds. Developments in Oil and Gas Basin Geology (in Chinese). Xi′an: Northwestern University Press, 235-247.

    Shen X J, Li G H, Wang J A, et al. 1994. Terrestrial heat flow measurement and calculation of statistical heat flow in Qaidam Basin.ChineseJournalofGeophysics(ActaGeophysicaSinica) (in Chinese), 37(1): 56-65.

    Sun G Q, Su L, Wang X H, et al. 2009. Fission track evidences of tectonic evolution in west Qaidam basin.NaturalGasIndustry(in Chinese), 29(2): 27-31.

    Tagami T, O’Sullivan P B. 2005. Fundamentals of fission-track thermochronology.Rev.Mineral.Geochem., 58(1): 19-47.

    Tang J G. 2007. tectonic evolution and its control for hydrocarbon accumulation of Mesozoic-Cenozoic multicycle superimposed reformation basin in the west of northern Qaidam basin [Ph. D. thesis] (in Chinese). Beijing: China University of Geosciences. Tang L J, Jin Z J, Zhang M L, et al. 2000. Tectonic evolution and oil (gas) pool-forming stage in northern Qaidam Basin.PetroleumExplorationandDevelopment(in Chinese), 27(2): 36-39.

    Tang X Y, Hu S B, Zhang G C, et al. 2014. Characteristic of surface heat flow in the Pearl River Mouth Basin and its relationship with thermal lithosphere thickness.ChineseJournalofGeophysics(in Chinese), 57(6): 1857-1867, doi: 10.6038/cjg20140617.

    Teng J W. 1974. Deep reflected waves and the structure of the earth crust of the eastern part of Chaidam Basin.ChineseJournalofGeophysics(ActaGeophysicaSinica) (in Chinese), 17(2): 122-135.

    Teng J W, Zhang Z J, Wang G J, et al. 1999. The deep internal dynamical processes and new model of continental-continental collision in Himalayan collision orogenic zone.ChineseJournalofGeophysics(in Chinese), 42(4): 481-494. Wang F, Luo Q H, Li Q, et al. 2001. The cooling event around 30 Ma in northern edge of Qaidam Basin-Constrains from40Ar/39Ar and Fission Track Thermochronology.BulletinofMineralogyPetrologyandGeochemistry(in Chinese), 20(4): 228-230.

    Wang H Z. 1982. The main stages of crustal development of China.EarthScience—JournalofWuhanCollegeofGeology(in Chinese), (3): 155-177.

    Wang H Z. 1985. Atlas of the Palaeogeography of China (in Chinese). Beijing: China Map Publishing House.

    Wang J, Huang S Y, Huang G S, et al. 1990. Basic Characteristics of China Temperature Distribution (in Chinese). Beijing: Seismological Press, 147-155.

    Wang J Y, Huang S P. 1990. Compilation of heat flow data in the continental area of China (2th Edition).SeismologyandGeology(in Chinese), 12(4): 351-363, 366.

    Wang S M, Ma C Q, Se Z B, et al. 2008. Apatite fission track analyses of Cenozoic sedimentary source and basin thermal history in West Qaidam Basin.GeologicalScienceandTechnologyInformation(in Chinese), 27(5): 29-36.

    Wang X M, Zhong D L, Wang Y. 2008. A case of application using apatite fission track to restrict the time of brittle fault movement.ProgressinGeophysics(in Chinese), 23(5): 1444-1455.

    Xu Z Q, Yang J T, Li H B, et al. 2006. The Qinghai-Tibet Plateau and continental dynamics: A review on terrain tectonics, collisional orogenesis, and processes and mechanisms for the rise of the plateau.GeologyinChina(in Chinese), 33(2): 221-238.

    Yang C, Chen Q H, Ren L Y, et al. 2012. Tectonic Units of the Qaidam Basin.JournalofSouthwestPetroleumUniversity(Science&TechnologyEdition) (in Chinese), 34(1): 25-33.

    Yang S C, Hu S B, Cai D S, et al. 2004. Present-day heat flow, thermal history and tectonic subsidence of the East China Sea Basin.MarineandPetroleumGeology, 21(9): 1095-1105.

    Yin A, Dang Y Q, Wang L C, et al. 2008. Cenozoic tectonic evolution of Qaidam basin and its surrounding regions (Part 1): The southern Qilian Shan-Nan Shan thrust belt and northern Qaidam basin.GeologicalSocietyofAmericaBulletin, 120(7-8): 813-846.

    Yu X Q, Liu J L, Zhang D H, et al. 2013. Uprising period and elevation of the Wenyu granitic pluton in the Xiaoqinling district, Central China.ChineseScienceBulletin, 58(35): 4459-4471.

    Yuan Y S, Zhu W L, Mi L J, et al. 2009. “Uniform geothermal gradient” and heat flow in the Qiongdongnan and Pearl River Mouth Basins of the South China Sea.MarineandPetroleumGeology, 26(7): 1152-1162.

    Zhang P Z, Zheng D W, Yin G M, et al. 2006. Discussion on late Cenozoic growth and rise of northeastern margin of the Tibetan plateau.QuaternarySciences(in Chinese), 26(1): 5-13.

    Zhang Y C, Hu J J, Liu C F. 1990. The Basic Characteristic of Earth Temperature and its Relationship with Oil and Gas in Qaidam Basin (in Chinese). Beijing: China Building Industry Press, 14-46.

    Zhao G Z, Tang J, Zhan Y, et al. 2004. A study on the relationship between the electrical structure of the crust and the deformation of the block in the northeastern margin of the Qinghai Tibet Plateau.ScienceinChina:Ser.D(in Chinese), 34(10): 908-918.

    Zhao H G, Liu C Y. 2003. Detachment gliding structures of Langgu sag.JournalofNorthwestUniversity(NaturalScienceEdition)(in Chinese),33(3): 315-319.

    Zhao J M, Tang W, Li Y S, et al. 2006. Lithospheric density and geomagnetic intensity in northeastern margin of the Tibetan plateau and tectonic implications.EarthScienceFrontiers(in Chinese), 13(5): 391-400.

    Zheng M L, Li M J, Cao C C, et al. 2004. Characteristics of structures of various levels in the Qaidam Cenozoic Basin.ActaGeologicaSinica(in Chinese), 78(1): 26-35.

    Zhou J X, Xu F Y, Hu Y. 2003. Mesozoic and Cenozoic tectonism and its control on hydrocarbon accumulation in the northern Qaidam Basin of China.ActaPetroleiSinica(in Chinese), 24(1): 19-24.

    Zhu L P, Helmberger D V. 1998. Moho offset across the northern margin of the Tibetan Plateau.Science, 281(5380): 1170-1172. Zuo Y H, Qiu N S, Zhang Y, et al. 2011. Geothermal regime and hydrocarbon kitchen evolution of the offshore Bohai Bay Basin, North China.AAPGBulletin, 95(5): 749-769.

    附中文參考文獻(xiàn)

    陳漢林, 陳沈強(qiáng), 林秀斌. 2014. 帕米爾弧形構(gòu)造帶新生代構(gòu)造演化研究進(jìn)展. 地球科學(xué)進(jìn)展, 29(8): 890-902.

    陳墨香. 1989. 確定中、新生代沉積盆地大地?zé)崃鞯姆椒? 地質(zhì)科學(xué), (2): 151-161.

    崔軍平, 任戰(zhàn)利, 肖暉等. 2007. 海拉爾盆地地溫分布及控制因素研究. 地質(zhì)科學(xué), 42(4): 656-665.

    馮昌格, 劉紹文, 王良書等. 2009. 塔里木盆地現(xiàn)今地?zé)崽卣? 地球物理學(xué)報(bào), 52(11): 2752-2762, doi: 10.3969/j.issn.0001-5733.2009.11.010.

    高成, 劉嬌, 李德威等. 2014. 裂變徑跡熱年代學(xué)對西藏卓奧友峰藏南拆離系活動時(shí)間的約束. 地學(xué)前緣, 21(6): 372-380.

    高軍平, 方小敏, 宋春暉等. 2011. 青藏高原北部中-新生代構(gòu)造-熱事件: 來自柴西碎屑磷灰石裂變徑跡的制約. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 41(5): 1466-1475.

    葛肖虹, 劉永江, 任收麥等. 2001. 對阿爾金斷裂科學(xué)問題的再認(rèn)識. 地質(zhì)科學(xué), 36(3): 319-325.

    李國樺. 1992. 柴達(dá)木盆地大地?zé)崃魈卣骷暗貧峤Y(jié)構(gòu)分析[博士論文]. 北京: 中國科學(xué)院地質(zhì)研究所.

    李永華, 吳慶舉, 安張輝等. 2006. 青藏高原東北緣地殼S波速度結(jié)構(gòu)與泊松比及其意義. 地 球物理學(xué)報(bào), 49(5): 1359-1368.

    李振華, 陳剛, 丁超等. 2014. 準(zhǔn)噶爾盆地北部燕山中期古地溫梯度恢復(fù). 地質(zhì)科技情報(bào), 33(4): 31-36.

    劉和甫. 2001. 盆地-山嶺耦合體系與地球動力學(xué)機(jī)制. 地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào), 26(6): 581-596.

    潘裕生. 1999. 青藏高原的形成與隆升. 地學(xué)前緣, 6(3): 153-163.

    邱楠生. 2001. 柴達(dá)木盆地現(xiàn)代大地?zé)崃骱蜕畈康販靥卣? 中國礦業(yè)大學(xué)學(xué)報(bào), 30(4): 412-415.

    饒松, 胡圣標(biāo), 朱傳慶等. 2013. 準(zhǔn)噶爾盆地大地?zé)崃魈卣髋c巖石圈熱結(jié)構(gòu). 地球物理學(xué)報(bào), 56(8): 2760-2770, doi: 10.6038/cjg20130824.

    任收麥, 葛肖虹, 劉永江等. 2009. 晚白堊世以來沿阿爾金斷裂帶的階段性走滑隆升. 地質(zhì)通報(bào), 23(9-10): 926-932.

    任戰(zhàn)利. 1993. 柴達(dá)木盆地地?zé)嵫莼? 來自流體包裹體和Ro的證據(jù).∥趙重遠(yuǎn), 劉池陽等. 含油氣盆地地質(zhì)學(xué)研究進(jìn)展. 西安: 西北大學(xué)出版社, 235-247.

    沈顯杰, 李國樺, 汪輯安等. 1994. 青海柴達(dá)木盆地大地?zé)崃鳒y量與統(tǒng)計(jì)熱流計(jì)算. 地球物理學(xué)報(bào), 37(1): 56-65.

    孫國強(qiáng), 蘇龍, 王旭紅等. 2009. 柴達(dá)木盆地西部地區(qū)構(gòu)造演化的裂變徑跡揭示. 天然氣工業(yè), 29(2): 27-31.

    湯濟(jì)廣. 2007. 柴達(dá)木北緣西段中、新生代多旋回疊加改造型盆地構(gòu)造演化及對油氣成藏的控制作用[博士論文]. 北京: 中國地質(zhì)大學(xué).

    湯良杰, 金之鈞, 張明利等. 2000. 柴達(dá)木盆地北緣構(gòu)造演化與油氣成藏階段. 石油勘探與開發(fā), 27(2): 36-39.

    唐曉音, 胡圣標(biāo), 張功成等. 2014. 珠江口盆地大地?zé)崃魈卣骷捌渑c熱巖石圈厚度的關(guān)系. 地球物理學(xué)報(bào), 57(6): 1857-1867, doi: 10.6038/cjg20140617.

    滕吉文. 1974. 柴達(dá)木東盆地的深層地震反射波和地殼構(gòu)造. 地球物理學(xué)報(bào), 17(2): 122-135.

    滕吉文, 張中杰, 王光杰等. 1999. 喜馬拉雅碰撞造山帶的深層動力過程與陸-陸碰撞新模型. 地球物理學(xué)報(bào), 42(4): 481-494.

    王非, 羅清華, 李齊等. 2001. 柴達(dá)木盆地北緣30 Ma左右的去頂剝蝕作用——40Ar/39Ar熱年代學(xué)定量制約. 礦物巖石地球化學(xué)通報(bào), 20(4): 228-230.

    王鴻禎. 1982. 中國地殼構(gòu)造發(fā)展的主要階段. 地球科學(xué)——武漢地質(zhì)學(xué)院學(xué)報(bào), (3): 155-177.

    王鴻禎. 1985. 中國古地理圖集. 北京: 地圖出版社.

    王鈞, 黃尚瑤, 黃歌山等. 1990. 中國地溫分布的基本特征. 北京: 地震出版社, 147-155.

    汪集旸, 黃少鵬. 1990. 中國大陸地區(qū)大地?zé)崃鲾?shù)據(jù)匯編(第二版). 地震地質(zhì), 12(4): 351-363, 366.

    王世明, 馬昌前, 佘振兵等. 2008. 柴西新生代沉積源區(qū)及盆地?zé)釟v史的磷灰石裂變徑跡分析. 地質(zhì)科技情報(bào), 27(5): 29-36.

    王先美, 鐘大賚, 王毅. 2008. 利用磷灰石裂變徑跡約束脆性斷裂活動的時(shí)限. 地球物理學(xué)進(jìn)展, 23(5): 1444-1455.

    徐志琴, 楊經(jīng)綏, 李海兵等. 2006. 青藏高原與大陸動力學(xué)——地體拼合、碰撞造山及高原隆升的深部驅(qū)動力. 中國地質(zhì), 33(2): 221-238.

    楊超, 陳清華, 任來義等. 2012. 柴達(dá)木盆地構(gòu)造單元劃分. 西南石油大學(xué)學(xué)報(bào)(自然科學(xué)版), 34(1): 25-33.

    余心起, 劉俊來, 張德會等. 2013. 小秦嶺文峪花崗巖山體的隆升時(shí)代和幅度. 科學(xué)通報(bào), 58(33): 3416-3428.

    張培震, 鄭德文, 尹功明等. 2006. 有關(guān)青藏高原東北緣晚新生代擴(kuò)展與隆升的討論. 第四紀(jì)研究, 26(1): 5-13.

    張業(yè)成, 胡景江, 劉春鳳. 1990. 柴達(dá)木盆地地溫基本特征及其與油氣關(guān)系. 北京: 中國建筑工業(yè)出版社, 14-46.

    趙國澤, 湯吉, 詹艷等. 2004. 青藏高原東北緣地殼電性結(jié)構(gòu)和地塊變形關(guān)系的研究. 中國科學(xué)(D輯), 34(10): 908-918.

    趙紅格, 劉池陽. 2003. 廊固凹陷的拆離滑脫構(gòu)造. 西北大學(xué)學(xué)報(bào)(自然科學(xué)版), 33(3): 315-319.

    趙俊猛, 唐偉, 黎益仕等. 2006. 青藏高原東北緣巖石圈密度與磁化強(qiáng)度及動力學(xué)含義. 地學(xué)前緣, 13(5): 391-400.

    鄭孟林, 李明杰, 曹春潮等. 2004. 柴達(dá)木盆地新生代不同層次構(gòu)造特征. 地質(zhì)學(xué)報(bào), 78(1): 26-35.

    周建勛, 徐鳳銀, 胡勇. 2003. 柴達(dá)木盆地北緣中、新生代構(gòu)造變形及其對油氣成藏的控制. 石油學(xué)報(bào), 24(1): 19-24.

    (本文編輯 何燕)

    Present-day heat flow and tectonic-thermal evolution since the late Paleozoic time of the Qaidam basin

    LI Zong-Xing1, GAO Jun2*, ZHENG Ce1,3, LIU Cheng-Lin1, MA Yin-Sheng1, ZHAO Wei-Yong4

    1KeyLaboratoryofShaleOilandGasGeologicalSurvey,InstituteofGeomechanics,ChineseAcademyofGeologicalSciences,Beijing100081,China2BeijingKeyLaboratoryofWaterResourcesandEnvironmentandEngineering,ChinaUniversityofGeosciences,Beijing100083,China3ChinaUniversityofPetroleum,Beijing102249,China4ResearchInstituteofPetroleumExplorationandDevelopment,QinghaiOilfield,CNPC,GansuDunhuang736202,China

    The Qaidam basin is the largest intermountain basin on the northern margin of the Tibetan plateau, and is one of three major petroliferous basins in western China. The characteristics of tectonic-thermal evolution of this basin are important for understanding the uplift mechanism of the Tibetan plateau. This work studies the terrestrial heat flow and thermal history of the basin, which would help research the regional dynamics and tectonic evolution of the basin and can be applied to the petroleum resource assessment in this region.This work employed the advanced temperature continuous acquisition system which measures the deep well steady state temperature at a large scale and high accuracy. Then mass of rock thermal conductivity was tested by the optical scanning method. Borehole temperature measurement data were fitted linearly using the least squares method, yielding 17 borehole profiles of geothermal gradients. Furthermore, original temperature measurement data from previous studies were re-analyzed, resulting in 56 geothermal gradients and the measured heat flow data of boreholes. Based on the present geothermal field, the time-temperature history of the samples was simulated by the combination of analysis of apatite fission track ages with modeling of fission track length distribution. For this simulation, the one-component annealing model, Monte Carlo method, and the software Hefty were used. Apatite samples enclosed track length values and the median age were adopted to simulate the basic parameters, with kinetic parameter (Dpar) as the initial value of 1.65, and the initial track length 16.3 μm. To the zircon fission track samples which did not pass2test (P(2)<5%), decomposition on single grain fission track age data and frequency statistics analysis were carried out. Based on the analysis of thermochronological dating, we have studied the burial history, uplift and denudation and thermal history of the Qaidam basin. Then we analyzed the tectonic evolution of the Qaidam basin combining with regional tectonic setting.Geothermal gradients of the Qaidam basin vary from 17.1 to 38.6 ℃·km-1with an average of 28.6±4.6 ℃·km-1. The thermal conductivity of the basin varies from 0.601 W/(m·K) to 5.520 W/(m·K) with an average of 2.208 W/(m·K), mostly in the range 1.820 W/(m·K) to 3.020 W/(m·K). Present-day terrestrial heat flow in this area varies from 32.9 to 70.4 mW·m-2, with an average of 55.1±7.9 mW·m-2. The test shows that the apatite fission track age varies from 61 Ma to 72 Ma, less than the formation age. The mean closed track length varies from 11.8 to 12.7μm, less than the initial length of the track (16.3±0.9 μm), implying a single-mode distribution. The test results of zircon samples show that the zircon fission track ages concentrate in 155 Ma to 195 Ma, less than the formation age, which have two time periods, i.e. from 254 Ma to 199 Ma and from 177 Ma to 148.6 Ma.There exits the distribution difference of the geothermal field in the Qaidam basin. The Kunbei thrust belt and Yiliping depression are “hot zones”, the Qilian thrust belt is a “warm zone”, and the Sanhu depression, Oulongbuluke uplift, and Delingha depression are “cold zones”. Such distributions of the geothermal field might be controlled by crustal structure in the basin. Based on the present-day geothermal field, the burial history, uplift and denudation and thermal history of Qaidam basin are studied by the combination of analysis of apatite fission track and zircon fission track ages with modeling of fission track length distribution. In terms of the analysis of thermochronological dating, we further explore the tectonic evolution of the Qaidam basin combining with regional tectonic setting. The results show that the Qaidam basin had experienced tectonic movements of six phases (254.0—199 Ma, 177—148.6 Ma, 87—62 Ma, 41.1—33.6 Ma, 9.6—7.1 Ma, 2.9—1.8 Ma); these tectonic events just coincided with the regional tectonic settings of the study area. The eastern Qaidam basin was uplifted and suffered denudation. The prototype of the Oulongbuluke uplift was also formed in the end of Cretaceous (41.1—33.6 Ma), however, the depression basin developed under the weak compressive tectonic regime in the northern Qaidam basin. The present tectonic features of the Qaidam basin appeared after uplift and great structural deformation, which was related with the intense extrusion since the end of Miocene (9.6—7.1 Ma and 2.9—1.8 Ma).

    Geothermal gradient; Heat flow; Fission track; Tectono-thermal evolution; Qaidam basin

    10.6038/cjg20151021.

    Li Z X, Gao J, Zheng C, et al. 2015. Present-day heat flow and tectonic-thermal evolution since the late Paleozoic time of the Qaidam basin.ChineseJ.Geophys. (in Chinese),58(10):3687-3705,doi:10.6038/cjg20151021.

    國家自然科學(xué)基金項(xiàng)目(41302202),中國地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所基本科研項(xiàng)目(DZLXJK201305)資助.

    李宗星,男,1982年生,構(gòu)造地質(zhì)學(xué)博士,助理研究員,主要從事盆地地溫場、盆地構(gòu)造-熱演化研究.E-mail:lizongxing@cags.ac.cn

    *通訊作者 高俊,女,1982年生,水文地質(zhì)學(xué)在讀博士,主要從事裂縫分析和盆地模擬研究.E-mail:gjpursuedoc@163.com

    10.6038/cjg20151021

    P314

    2015-01-28,2015-06-06收修定稿

    李宗星, 高俊, 鄭策等. 2015. 柴達(dá)木盆地現(xiàn)今大地?zé)崃髋c晚古生代以來構(gòu)造-熱演化.地球物理學(xué)報(bào),58(10):3687-3705,

    猜你喜歡
    徑跡柴達(dá)木盆地熱導(dǎo)率
    基于蒙特卡羅模擬方法的圓筒形固體核徑跡氡探測器探測效率的研究
    空位缺陷對單層石墨烯導(dǎo)熱特性影響的分子動力學(xué)
    裂變徑跡LA-ICP-MS/FT法原理、實(shí)驗(yàn)流程和應(yīng)用
    柴達(dá)木盆地北緣鋰多金屬礦成礦條件及找礦潛力
    連續(xù)碳纖維鋁基復(fù)合材料橫向等效熱導(dǎo)率的模擬分析
    Si3N4/BN復(fù)合陶瓷熱導(dǎo)率及其有限元分析
    柴達(dá)木盆地西部古近系石膏及其硫同位素分布特征
    金屬熱導(dǎo)率的第一性原理計(jì)算方法在鋁中的應(yīng)用
    核乳膠中質(zhì)子徑跡計(jì)算機(jī)輔助識別方法研究
    柴達(dá)木盆地柴北緣侏羅系烴源巖特征分析
    成年美女黄网站色视频大全免费 | 午夜老司机福利剧场| 成人无遮挡网站| 亚洲人与动物交配视频| av黄色大香蕉| 插逼视频在线观看| 九九久久精品国产亚洲av麻豆| 久久影院123| av电影中文网址| 国产精品久久久久久精品古装| 国产一区二区三区av在线| 久久精品国产亚洲av天美| 丁香六月天网| 一本久久精品| 亚洲欧洲日产国产| 最新中文字幕久久久久| 婷婷色综合大香蕉| 视频区图区小说| 国产精品秋霞免费鲁丝片| 国产欧美日韩综合在线一区二区| 欧美精品一区二区免费开放| 五月伊人婷婷丁香| 成人国产麻豆网| 五月玫瑰六月丁香| 免费观看性生交大片5| 亚洲国产精品成人久久小说| 男人添女人高潮全过程视频| 中文天堂在线官网| 18禁裸乳无遮挡动漫免费视频| 啦啦啦视频在线资源免费观看| 欧美3d第一页| 国产精品99久久久久久久久| 亚洲精品久久久久久婷婷小说| 国产精品三级大全| 一级黄片播放器| h视频一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 色吧在线观看| 国产精品熟女久久久久浪| 少妇精品久久久久久久| 三级国产精品欧美在线观看| 97在线视频观看| 老司机影院成人| 亚洲国产精品国产精品| 精品国产露脸久久av麻豆| 精品一区二区免费观看| 亚洲欧美清纯卡通| 国产淫语在线视频| 汤姆久久久久久久影院中文字幕| 亚洲情色 制服丝袜| 97精品久久久久久久久久精品| 久久婷婷青草| 男女边摸边吃奶| 日韩av免费高清视频| av不卡在线播放| 在线观看一区二区三区激情| 高清黄色对白视频在线免费看| 女人精品久久久久毛片| 国模一区二区三区四区视频| 欧美一级a爱片免费观看看| 国产精品久久久久久久电影| 国产成人精品婷婷| 韩国av在线不卡| 国产精品国产av在线观看| 久久久久久久久久久久大奶| 久久午夜综合久久蜜桃| 纯流量卡能插随身wifi吗| 不卡视频在线观看欧美| 菩萨蛮人人尽说江南好唐韦庄| 卡戴珊不雅视频在线播放| 亚洲国产成人一精品久久久| 亚洲人成77777在线视频| 草草在线视频免费看| 91精品国产国语对白视频| 亚洲av不卡在线观看| 欧美日本中文国产一区发布| 精品久久久久久电影网| 久久国内精品自在自线图片| 欧美变态另类bdsm刘玥| 亚洲色图 男人天堂 中文字幕 | 各种免费的搞黄视频| 成人无遮挡网站| 精品亚洲成国产av| 国产成人精品无人区| 国产一区亚洲一区在线观看| 午夜激情福利司机影院| 亚洲久久久国产精品| 大香蕉久久成人网| 综合色丁香网| 亚洲精品一区蜜桃| 人妻夜夜爽99麻豆av| 国产精品免费大片| 在现免费观看毛片| 男人爽女人下面视频在线观看| 五月开心婷婷网| 日韩成人伦理影院| videosex国产| 一区二区日韩欧美中文字幕 | 制服人妻中文乱码| 在线观看免费高清a一片| 在线 av 中文字幕| 国产精品一二三区在线看| 丝瓜视频免费看黄片| 久久这里有精品视频免费| 在线精品无人区一区二区三| 国产精品久久久久久久电影| 亚洲精品亚洲一区二区| 精品国产一区二区久久| 91国产中文字幕| 亚洲一区二区三区欧美精品| 女性生殖器流出的白浆| 成人免费观看视频高清| 熟女电影av网| 午夜福利网站1000一区二区三区| 亚洲久久久国产精品| 新久久久久国产一级毛片| 成人午夜精彩视频在线观看| 看免费成人av毛片| 边亲边吃奶的免费视频| 日本与韩国留学比较| 亚洲国产精品专区欧美| 日本欧美国产在线视频| 少妇的逼水好多| 只有这里有精品99| 性色avwww在线观看| 国产一区有黄有色的免费视频| 大话2 男鬼变身卡| 日韩强制内射视频| 特大巨黑吊av在线直播| 精品国产露脸久久av麻豆| 777米奇影视久久| 熟女人妻精品中文字幕| 国产欧美亚洲国产| 亚洲综合色网址| 狂野欧美白嫩少妇大欣赏| 热99久久久久精品小说推荐| 国产极品粉嫩免费观看在线 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品视频人人做人人爽| 哪个播放器可以免费观看大片| 国产成人免费无遮挡视频| 中文字幕人妻熟人妻熟丝袜美| 国产免费一级a男人的天堂| 国产亚洲最大av| 欧美3d第一页| 亚洲av日韩在线播放| 欧美国产精品一级二级三级| 国产欧美日韩一区二区三区在线 | 少妇 在线观看| 人妻一区二区av| 精品国产乱码久久久久久小说| 亚洲国产最新在线播放| 18+在线观看网站| 五月伊人婷婷丁香| 美女国产高潮福利片在线看| 亚洲中文av在线| 中国三级夫妇交换| 一边摸一边做爽爽视频免费| 91精品一卡2卡3卡4卡| 久久鲁丝午夜福利片| 国产又色又爽无遮挡免| www.av在线官网国产| 国产精品一国产av| 成人国语在线视频| 青春草亚洲视频在线观看| 国产亚洲精品久久久com| 免费看av在线观看网站| 中文字幕久久专区| 伦理电影大哥的女人| 亚洲精品久久久久久婷婷小说| 成人国语在线视频| 18禁在线播放成人免费| 日韩亚洲欧美综合| 中文字幕制服av| 少妇被粗大的猛进出69影院 | 欧美一级a爱片免费观看看| 另类亚洲欧美激情| 亚洲国产av影院在线观看| 搡老乐熟女国产| 欧美3d第一页| 国产在线免费精品| 只有这里有精品99| 国产片内射在线| 久久人人爽人人片av| 99久久精品一区二区三区| 国产免费福利视频在线观看| 各种免费的搞黄视频| 搡老乐熟女国产| av卡一久久| 国产精品久久久久久久久免| 国产午夜精品一二区理论片| 高清黄色对白视频在线免费看| 免费不卡的大黄色大毛片视频在线观看| xxx大片免费视频| 亚洲经典国产精华液单| 成人亚洲欧美一区二区av| 人妻夜夜爽99麻豆av| 亚洲av中文av极速乱| 亚洲av成人精品一二三区| 国产精品欧美亚洲77777| 99热全是精品| 精品国产露脸久久av麻豆| 青春草亚洲视频在线观看| 久久精品国产亚洲av天美| 久久久久久久久久久久大奶| 亚洲,一卡二卡三卡| 亚洲国产毛片av蜜桃av| 久久精品久久久久久噜噜老黄| 亚洲国产精品成人久久小说| 色94色欧美一区二区| 国产国拍精品亚洲av在线观看| 飞空精品影院首页| 亚洲激情五月婷婷啪啪| 亚洲综合精品二区| 国产伦理片在线播放av一区| xxx大片免费视频| 成人二区视频| 日韩av在线免费看完整版不卡| 高清不卡的av网站| 久久毛片免费看一区二区三区| 亚洲伊人久久精品综合| 久久久精品区二区三区| 国产成人精品无人区| 啦啦啦视频在线资源免费观看| 日本91视频免费播放| 汤姆久久久久久久影院中文字幕| 久久久久精品性色| 女人久久www免费人成看片| 日韩三级伦理在线观看| 久久久久久久久久久丰满| 日韩一本色道免费dvd| 亚洲,一卡二卡三卡| 精品人妻在线不人妻| 搡女人真爽免费视频火全软件| 欧美精品亚洲一区二区| 黑丝袜美女国产一区| 国产高清国产精品国产三级| 国产日韩欧美在线精品| av福利片在线| 性色av一级| 大话2 男鬼变身卡| 午夜91福利影院| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费少妇av软件| av免费观看日本| 亚洲av不卡在线观看| 高清午夜精品一区二区三区| 99久久综合免费| 热99久久久久精品小说推荐| 国产精品免费大片| 午夜福利视频在线观看免费| 日韩精品有码人妻一区| 少妇人妻精品综合一区二区| 午夜福利在线观看免费完整高清在| 日日爽夜夜爽网站| 麻豆成人av视频| 国产欧美另类精品又又久久亚洲欧美| 亚洲少妇的诱惑av| 五月伊人婷婷丁香| 简卡轻食公司| 久久久久久久久久人人人人人人| 精品少妇久久久久久888优播| 国产免费又黄又爽又色| 亚洲综合色惰| 99热这里只有是精品在线观看| 人妻一区二区av| 婷婷色综合大香蕉| 午夜免费鲁丝| 99视频精品全部免费 在线| 一本—道久久a久久精品蜜桃钙片| 熟女人妻精品中文字幕| 亚洲精品乱码久久久v下载方式| 国产精品人妻久久久久久| 母亲3免费完整高清在线观看 | 下体分泌物呈黄色| 少妇被粗大的猛进出69影院 | 大又大粗又爽又黄少妇毛片口| 少妇高潮的动态图| 精品国产一区二区久久| 熟女人妻精品中文字幕| 国产伦理片在线播放av一区| 日产精品乱码卡一卡2卡三| av天堂久久9| 国产精品久久久久久精品电影小说| 欧美亚洲日本最大视频资源| 日韩一区二区视频免费看| 免费人妻精品一区二区三区视频| 自拍欧美九色日韩亚洲蝌蚪91| 丰满迷人的少妇在线观看| 丝瓜视频免费看黄片| 免费观看av网站的网址| 各种免费的搞黄视频| 亚洲欧美中文字幕日韩二区| 在线精品无人区一区二区三| 国产亚洲欧美精品永久| 99热这里只有是精品在线观看| 亚洲久久久国产精品| 26uuu在线亚洲综合色| 大香蕉97超碰在线| 久久毛片免费看一区二区三区| 91国产中文字幕| 亚洲,一卡二卡三卡| 亚洲人成网站在线播| 国产有黄有色有爽视频| 日本猛色少妇xxxxx猛交久久| 春色校园在线视频观看| 十分钟在线观看高清视频www| 日韩欧美精品免费久久| 日韩视频在线欧美| freevideosex欧美| 18禁裸乳无遮挡动漫免费视频| 亚洲av免费高清在线观看| 久热久热在线精品观看| 青春草视频在线免费观看| 国产熟女欧美一区二区| 国产亚洲午夜精品一区二区久久| 国产乱来视频区| 久久精品人人爽人人爽视色| 草草在线视频免费看| 日韩亚洲欧美综合| 日日摸夜夜添夜夜添av毛片| 又大又黄又爽视频免费| 午夜91福利影院| 两个人免费观看高清视频| 黄片无遮挡物在线观看| 男男h啪啪无遮挡| 久热这里只有精品99| 精品久久久久久久久亚洲| 夜夜骑夜夜射夜夜干| 99久久精品一区二区三区| 国产精品无大码| 插阴视频在线观看视频| 欧美日本中文国产一区发布| 色视频在线一区二区三区| 国产午夜精品一二区理论片| 国产高清国产精品国产三级| 国产精品国产av在线观看| 国产成人精品婷婷| h视频一区二区三区| 亚洲精品色激情综合| 国产乱人偷精品视频| 久久精品久久久久久噜噜老黄| 另类精品久久| 亚洲色图 男人天堂 中文字幕 | 精品久久久精品久久久| 成人免费观看视频高清| 观看美女的网站| 国产又色又爽无遮挡免| 亚洲欧洲精品一区二区精品久久久 | 国产精品99久久99久久久不卡 | 欧美日韩亚洲高清精品| 国产一区二区在线观看日韩| 免费黄网站久久成人精品| 2022亚洲国产成人精品| 欧美3d第一页| 亚洲精品日韩av片在线观看| 国产黄色视频一区二区在线观看| 男女边摸边吃奶| 熟女人妻精品中文字幕| 天天躁夜夜躁狠狠久久av| 热99国产精品久久久久久7| 99国产综合亚洲精品| 国产免费视频播放在线视频| 亚洲国产日韩一区二区| 亚洲精品av麻豆狂野| 欧美人与善性xxx| 欧美性感艳星| 亚洲在久久综合| 校园人妻丝袜中文字幕| 免费大片黄手机在线观看| 国精品久久久久久国模美| 欧美日本中文国产一区发布| 啦啦啦视频在线资源免费观看| 精品酒店卫生间| 久久精品久久精品一区二区三区| 国产精品.久久久| 黄色毛片三级朝国网站| 亚洲av二区三区四区| 丝袜美足系列| 国产av精品麻豆| 18禁在线无遮挡免费观看视频| 久久国产精品大桥未久av| 美女内射精品一级片tv| 精品国产一区二区三区久久久樱花| 久久午夜福利片| 欧美人与善性xxx| 亚州av有码| 久久久久精品性色| 国产精品久久久久久av不卡| 亚洲精品国产av蜜桃| 乱码一卡2卡4卡精品| 最新的欧美精品一区二区| 中文字幕最新亚洲高清| 丰满饥渴人妻一区二区三| 国产精品一区二区在线观看99| 日日撸夜夜添| 蜜桃久久精品国产亚洲av| 我的老师免费观看完整版| 久久影院123| 国产日韩欧美亚洲二区| 国产精品 国内视频| 王馨瑶露胸无遮挡在线观看| 国产午夜精品一二区理论片| av免费在线看不卡| 欧美丝袜亚洲另类| 欧美三级亚洲精品| 天美传媒精品一区二区| 日韩精品免费视频一区二区三区 | 亚洲国产日韩一区二区| 日本猛色少妇xxxxx猛交久久| 女的被弄到高潮叫床怎么办| 国产免费一级a男人的天堂| 色婷婷久久久亚洲欧美| 亚洲av中文av极速乱| 国产乱来视频区| 一区在线观看完整版| 女的被弄到高潮叫床怎么办| 亚洲伊人久久精品综合| 男人添女人高潮全过程视频| 男女国产视频网站| 国产探花极品一区二区| 久久韩国三级中文字幕| 人体艺术视频欧美日本| 国产精品久久久久久精品电影小说| 久久久久久久久久久久大奶| 另类精品久久| 少妇人妻 视频| 99久久人妻综合| 亚洲精品色激情综合| 男女无遮挡免费网站观看| 91久久精品国产一区二区三区| 97超碰精品成人国产| 欧美性感艳星| 久久久久久久久久久免费av| av电影中文网址| 人成视频在线观看免费观看| 人妻系列 视频| 男人操女人黄网站| 亚洲少妇的诱惑av| 国产亚洲欧美精品永久| 国产成人一区二区在线| 91久久精品国产一区二区成人| 免费高清在线观看日韩| 女性被躁到高潮视频| 日韩av不卡免费在线播放| 国产亚洲一区二区精品| 在线观看人妻少妇| 亚洲伊人久久精品综合| 精品一区二区三卡| 人人澡人人妻人| 九草在线视频观看| 99久国产av精品国产电影| 成人毛片a级毛片在线播放| 亚洲五月色婷婷综合| 国产 一区精品| 欧美成人午夜免费资源| 亚洲av福利一区| 一区二区日韩欧美中文字幕 | 大片免费播放器 马上看| 你懂的网址亚洲精品在线观看| 美女国产高潮福利片在线看| 亚洲精品,欧美精品| 免费黄网站久久成人精品| 成人国产av品久久久| 亚洲精品日韩在线中文字幕| 国产高清国产精品国产三级| 男人添女人高潮全过程视频| 亚洲欧美成人综合另类久久久| 啦啦啦啦在线视频资源| 亚洲精品aⅴ在线观看| 在现免费观看毛片| 大陆偷拍与自拍| av国产精品久久久久影院| 丝袜在线中文字幕| 日产精品乱码卡一卡2卡三| 久久韩国三级中文字幕| 一个人看视频在线观看www免费| av又黄又爽大尺度在线免费看| 日日啪夜夜爽| 日韩精品有码人妻一区| 亚洲精品中文字幕在线视频| 婷婷色麻豆天堂久久| 精品一区二区三区视频在线| 国产av码专区亚洲av| 高清不卡的av网站| 91精品三级在线观看| 亚洲情色 制服丝袜| 国产又色又爽无遮挡免| 日本猛色少妇xxxxx猛交久久| 亚洲国产毛片av蜜桃av| 一本色道久久久久久精品综合| 中文字幕免费在线视频6| 五月伊人婷婷丁香| 国产高清三级在线| 国产片内射在线| 看十八女毛片水多多多| 老司机亚洲免费影院| av播播在线观看一区| 九九爱精品视频在线观看| 亚洲av福利一区| 亚洲久久久国产精品| 亚洲国产最新在线播放| 最新的欧美精品一区二区| 18禁观看日本| 久久久国产一区二区| 一个人看视频在线观看www免费| 午夜老司机福利剧场| 久久久亚洲精品成人影院| 女的被弄到高潮叫床怎么办| 爱豆传媒免费全集在线观看| 久久久久久久久久成人| 午夜福利影视在线免费观看| freevideosex欧美| 97在线人人人人妻| 黑人猛操日本美女一级片| 少妇高潮的动态图| 国产片特级美女逼逼视频| 亚洲熟女精品中文字幕| 免费黄频网站在线观看国产| 制服人妻中文乱码| 99精国产麻豆久久婷婷| 少妇人妻精品综合一区二区| 极品人妻少妇av视频| 精品人妻在线不人妻| 国产av国产精品国产| 久久99一区二区三区| 亚洲成人一二三区av| 黑人高潮一二区| 尾随美女入室| 国语对白做爰xxxⅹ性视频网站| 亚洲精品乱码久久久久久按摩| 久久久久久久大尺度免费视频| 精品午夜福利在线看| 日韩一区二区视频免费看| 亚洲色图综合在线观看| 91久久精品电影网| 亚洲精品国产av蜜桃| 国产精品三级大全| 女性生殖器流出的白浆| 51国产日韩欧美| 精品国产乱码久久久久久小说| 久久婷婷青草| 中文字幕av电影在线播放| 我的老师免费观看完整版| 中文字幕亚洲精品专区| 国国产精品蜜臀av免费| 国产男女超爽视频在线观看| 三级国产精品片| 最近最新中文字幕免费大全7| 人妻一区二区av| 久久国内精品自在自线图片| 欧美人与性动交α欧美精品济南到 | 亚洲怡红院男人天堂| 少妇被粗大的猛进出69影院 | 丝袜脚勾引网站| 精品国产乱码久久久久久小说| 久久久a久久爽久久v久久| 日韩中文字幕视频在线看片| 秋霞伦理黄片| xxx大片免费视频| √禁漫天堂资源中文www| 国产日韩欧美在线精品| 日韩一区二区视频免费看| 大码成人一级视频| 国产色婷婷99| 精品久久蜜臀av无| 中文字幕久久专区| 丰满乱子伦码专区| a级片在线免费高清观看视频| 黑人欧美特级aaaaaa片| 伦理电影大哥的女人| 午夜91福利影院| 欧美三级亚洲精品| 简卡轻食公司| 中文欧美无线码| 久久久久久久大尺度免费视频| 搡女人真爽免费视频火全软件| videos熟女内射| 国产亚洲欧美精品永久| 少妇的逼水好多| 综合色丁香网| 亚洲人成77777在线视频| 一本大道久久a久久精品| 亚洲精品aⅴ在线观看| 亚洲精品成人av观看孕妇| 欧美日韩精品成人综合77777| av线在线观看网站| 99久久中文字幕三级久久日本| 久久这里有精品视频免费| 18禁在线播放成人免费| 国产在线免费精品| 欧美 亚洲 国产 日韩一| 中文天堂在线官网| 久久久精品区二区三区| 少妇丰满av| 黑人巨大精品欧美一区二区蜜桃 | 人人妻人人澡人人看| 狠狠精品人妻久久久久久综合| 久久久久久伊人网av| 久久精品国产鲁丝片午夜精品| 观看av在线不卡| 自线自在国产av| 日韩三级伦理在线观看| 精品一品国产午夜福利视频| 欧美日韩视频高清一区二区三区二| 久久久久久久久久久免费av| 热99国产精品久久久久久7| freevideosex欧美| 999精品在线视频| 亚洲怡红院男人天堂| 国产成人精品一,二区| 亚洲人成77777在线视频| 亚洲国产欧美在线一区| 考比视频在线观看| 黑人欧美特级aaaaaa片|