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

    二連盆地白音查干凹陷和烏里雅斯太凹陷巖石熱導(dǎo)率測試與分析*

    2020-03-01 13:58:36余如洋黃少鵬張炯許威柯婷婷左銀輝周勇水
    巖石學(xué)報 2020年1期
    關(guān)鍵詞:查干白音熱導(dǎo)率

    余如洋 黃少鵬, 2, 3 張炯 許威 柯婷婷 左銀輝 周勇水

    1. 西安交通大學(xué)人居環(huán)境與建筑工程學(xué)院,西安 7100492. 深圳大學(xué)深地科學(xué)與綠色能源研究院,深圳 5180603. 美國密歇根大學(xué)地球與環(huán)境科學(xué)系,安娜堡 MI 48109-10054. 成都理工大學(xué)能源學(xué)院,成都 6100595. 中國石化中原油田,濮陽 4570011.

    巖石的熱物性參數(shù)主要包括巖石熱導(dǎo)率、生熱率、熱擴(kuò)散率、熱容、比熱、密度等,其中巖石熱導(dǎo)率是研究一個地區(qū)大地?zé)崃鳌⑸畈繜釥顟B(tài)和巖石圈熱結(jié)構(gòu)及地?zé)豳Y源評價中不可或缺的參數(shù)。由于熱流值是巖石熱導(dǎo)率與地溫梯度的乘積,如果熱導(dǎo)率不準(zhǔn)確,計(jì)算得到的熱流值必然偏離研究區(qū)域的真實(shí)情況,依據(jù)地表熱流、熱導(dǎo)率和生熱率計(jì)算的深部地溫場也會出現(xiàn)偏差。

    巖石熱導(dǎo)率的大小主要取決于巖石的礦物組成,但是也受到巖石結(jié)構(gòu)和構(gòu)造、孔隙度、飽水狀況、溫度和壓力等影響(Clauser and Huenges, 1995; Alishaevetal., 2012; Pasqualeetal., 2015; Ghanbarian and Daigle, 2016)。巖石樣品熱導(dǎo)率通常是在室溫常壓下對不含水的干巖樣進(jìn)行測試得到的,而干巖樣的熱導(dǎo)率值一般都比原位(自然狀態(tài))持有水分巖石的值偏小,其原因在于不含水巖樣孔隙中充填著空氣,空氣熱導(dǎo)率為0.025W/(m·K),該值小于水的熱導(dǎo)率0.607W/(m·K)(Nagaraju and Roy, 2014)。此外,巖石在地下原位的溫度和壓力狀況明顯有別于室溫常壓條件。有研究(Abdulagatovaetal., 2010)表明,在溫度為15~250℃之間,壓力在100MPa范圍內(nèi),壓力每增加10MPa,巖石的熱導(dǎo)率升高1.3%~2.0%;溫度每升高100K,熱導(dǎo)率較常溫下降4.8%~6%。因此,要得到切合實(shí)際的熱導(dǎo)率數(shù)據(jù),需要對有孔隙或裂隙巖石樣品的熱導(dǎo)率數(shù)據(jù)進(jìn)行飽水校正和溫壓校正。目前已經(jīng)有很多基于理論研究的模型和經(jīng)驗(yàn)公式用于計(jì)算巖石的原位熱導(dǎo)率,但是這些模型和經(jīng)驗(yàn)公式都存在一定的局限性(Abdulagatovaetal., 2009)。不同的研究區(qū)需要應(yīng)用不同的方法進(jìn)行對比分析,篩選合適的模型或者經(jīng)驗(yàn)公式進(jìn)行校正。

    二連盆地位于內(nèi)蒙-大興安嶺古生代碰撞造山帶,是一個中-新生代沉積斷陷盆地,已經(jīng)探明含有豐富的油氣資源,但迄今地?zé)嵫芯砍潭热匀缓艿?,熱物性?shù)據(jù)十分稀少,這嚴(yán)重制約了二連盆地大地?zé)崃?、巖石圈熱結(jié)構(gòu)分析及地?zé)豳Y源潛力評估等研究工作的展開。比如根據(jù)幾個零星的低熱流估算值(肖偉等, 2004),前人把二連盆地定性為冷盆(藺文靜等, 2013)。可是盆地周邊新生代火山活動強(qiáng)烈,而且與二連盆地相鄰的銀根-額濟(jì)納旗盆地屬于高熱流區(qū)域(左銀輝等, 2013),二連盆地內(nèi)平均地溫梯度大于35℃/km(Xuetal., 2019; Zhangetal., 2020),顯然與“冷盆”不相稱。但是目前未有實(shí)測熱導(dǎo)率數(shù)據(jù)的報道,這更加突顯了豐富二連盆地巖石熱導(dǎo)率數(shù)據(jù)的重要性。本文在室溫常壓下對白音查干凹陷和烏里雅斯太凹陷98塊鉆孔巖芯樣品進(jìn)行了巖石熱導(dǎo)率測試,并搜集了中原油田在烏里雅斯太凹陷已有的31塊巖芯熱導(dǎo)率數(shù)據(jù),分析了這兩個地區(qū)的熱導(dǎo)率特征,并采用簡化的有效熱導(dǎo)率理論公式對實(shí)測熱導(dǎo)率進(jìn)行飽水校正。在此基礎(chǔ)上,利用Abdulagatovaetal. (2010)建立的溫壓校正公式對其進(jìn)行井下原位溫度和壓力校正,由此得到巖石原位熱導(dǎo)率,結(jié)合地層巖性組成,建立了白音查干凹陷和烏里雅斯太凹陷的地層熱導(dǎo)率,為這兩個地區(qū)的后續(xù)大地?zé)崃饔?jì)算和相關(guān)研究提供可靠的巖石熱導(dǎo)率基礎(chǔ)數(shù)據(jù)。

    1 地質(zhì)背景

    二連盆地位于華北板塊與西伯利亞板塊所夾持的中亞造山帶東段,是在晚古生界褶皺基底和侏羅系殘留盆地基礎(chǔ)上,經(jīng)由中生代后期強(qiáng)烈拉伸、裂陷而形成的早白堊世中、小型裂谷盆地群(程三友等, 2011)。

    圖1 研究區(qū)構(gòu)造單元劃分與地層剖面 (a)二連盆地構(gòu)造單元(劉震等, 2007);(b)白音查干凹陷構(gòu)造帶劃分及熱物性測試取樣井位置;(c)烏里雅斯太凹陷洼槽分布及熱物性測試取樣井位置;(d)白音查干凹陷A-A′剖面地層(Zuo et al., 2016);(e)烏里雅斯太凹陷B-B′剖面地層(易士威等, 2006)Fig.1 Schematic tectonic map and formation profiles of the study areas (a) tectonic division of the Erlian Basin (modified after Liu et al., 2007); (b) and (c) divisions of the structural zones and the locations of the core sampling boreholes of the Baiyinchagan Sag and Uliastai Sag, respectively; (d) A-A′ formation profile of the Baiyinchagan Sag (modified after Zuo et al., 2016); (e) B-B′ formation profile of the Uliastai Sag (modified after Yi et al., 2006)

    白音查干凹陷位于二連盆地西緣(圖1a),為受北部塔拉斷層和南部白音翁特?cái)鄬庸餐刂频幕畎枷?,呈東西向展布,是發(fā)育在早白堊世時期的陸相斷陷型盆地(Bonnettietal., 2014)。盆地基底是古生界地層Pz,為淺海相和海陸交互相碎屑巖夾碳酸鹽巖、巖漿巖、火山碎屑巖建造,剝蝕比較嚴(yán)重(崔永謙等, 2011)。古生界以上即進(jìn)入白堊系,缺失三疊系和侏羅系,自上而下沉積層是二連達(dá)布蘇組(K2er)、賽漢塔拉組(K1bs)、都紅木組(K1bd)、騰格爾組(K1bt)及阿爾善組(K1ba)(圖1d)(趙澄林等, 1996)。二連達(dá)布蘇組(K2er)為河流相沉積,巖性為含礫泥巖、泥巖和砂(礫)巖不等厚互層。賽罕塔拉組(K1bs)以河流相沉積為主,巖性以含礫砂巖、粗砂巖為主,夾薄層泥巖或含礫泥巖。都紅木組(K1bd)可分為一段和二段,都紅木組一段為半深湖、深湖亞相沉積,主要以泥巖為主,夾粉砂巖和泥質(zhì)粉砂巖;都二段為湖泊相沉積,以泥巖、白云質(zhì)泥巖為主夾粉砂巖和泥質(zhì)粉砂巖等。騰格爾組(K1bt)為淺湖-深湖亞相沉積,以泥巖、白云質(zhì)泥巖為主,夾砂、礫巖。阿爾善組(K1ba)早期為河流濱湖相沉積,主要以礫狀砂巖、砂、礫巖為主,夾有薄層泥巖、粉砂巖等,中晚期則為湖相沉積,巖性主體為暗色泥巖夾砂巖。

    烏里雅斯太凹陷位于二連盆地東北角,為北東走向的西北斷東南超的單斷箕狀凹陷,自南向北可以分為南洼槽、中洼槽和北洼槽(圖1c)。其地層自上而下劃分為賽罕塔拉組(K1bs)、騰格爾組二段(K1bt2)騰格爾組一段(K1bt1)和阿爾善組(K1ba)(圖1e)(趙澄林等, 1996)。賽罕塔拉組(K1bs)為河流、沼澤相沉積,巖性由砂質(zhì)礫巖、含礫砂巖和砂巖與泥巖、含礫泥巖組成,構(gòu)成了兩個下粗上細(xì)的次級正旋回層,該組巖性比較突出的特征是砂巖疏松、泥巖含礫。騰格爾組二段(K1bt2)以廣湖盆淺水體的高水位體系域沉積為主。巖性由砂質(zhì)礫巖、含礫砂巖、粉砂巖與泥巖組成下粗上細(xì)的正旋回層。騰格爾組一段(K1bt1)總體上為較深湖泊相環(huán)境下沉積的泥巖,其中騰一(K1bt1)下部以深-淺湖相泥巖沉積為主,含砂礫巖、細(xì)砂巖;騰一(K1bt1)中下部為水進(jìn)退積型沉積,泥巖發(fā)育;中部泥巖與砂礫巖互層沉積;上部為水退進(jìn)積型沉積過程,以砂巖、砂礫巖與泥巖間互沉積為主。阿爾善組(K1ba)早期快速充填了一套以河流相為主的粗碎屑。中期發(fā)生湖侵,后期整體抬升,遭受剝蝕。其巖性以粉砂巖、細(xì)砂巖、砂礫巖與泥巖互層沉積為主。該凹陷中洼槽構(gòu)造發(fā)育的典型特征是賽漢塔拉組(K1bs)缺失。此外,烏里雅斯太凹陷騰格爾組一段相當(dāng)于白音查干凹陷的騰格爾組,騰格爾組二段相當(dāng)于白音查干凹陷的都紅木組。

    2 熱導(dǎo)率測試結(jié)果

    2.1 樣品選取與測試

    本次巖芯測試樣品源自中國石油化工集團(tuán)公司中原油田巖芯庫,采樣原則是樣品盡可能覆蓋盆地內(nèi)出露的地層,巖芯類型上盡可能包括盆地內(nèi)部各地層代表巖性。在烏里雅斯太凹陷地區(qū),我們選取了7口鉆井的55塊巖芯樣品,其巖性主要是砂巖、泥巖和礫巖。我們還收集了中原油田已有的15口鉆井的31塊巖芯的熱導(dǎo)率數(shù)據(jù),其巖性主要是砂巖,地層包括白堊系騰格爾組和阿爾善組。在白音查干凹陷地區(qū),選取了6口鉆井的43塊巖芯樣品,其巖性主要是砂巖、泥巖,含少量的礫巖、片巖和玄武巖,其地層包括白堊系賽漢塔拉組、都紅木組、騰格爾組和阿爾善組,具體的井位信息和各鉆井巖芯選取數(shù)量見表1。

    本文熱導(dǎo)率測試儀器為德國Lippmann and Rauen GbR公司生產(chǎn)的光學(xué)掃描熱導(dǎo)儀TCS(Thermal Conductivity Scanning)。該款儀器熱導(dǎo)率測試范圍為0.2~25.0W/(m·K),測試分辨率為0.001W/(m·K),測量精度為±3%。這一類型儀器已成功應(yīng)用于國內(nèi)外巖石熱導(dǎo)率測試和研究工作當(dāng)中,如俄羅斯科拉半島的超深鉆(Popovetal., 1999)及中國蘇魯-大別變質(zhì)帶的大陸超深鉆孔(Heetal., 2008)、新疆塔里木盆地(Liuetal., 2011)、青海共和盆地(Zhangetal., 2018)等。本文巖芯測試樣品形狀均為圓柱形,直徑為6~9cm,厚度為3~6cm。為了使每個巖芯樣品表面粗糙度在1mm以內(nèi),對每個樣品均進(jìn)行了拋光處理。此外,為了減小巖樣表面不同光反射系數(shù)的影響,在拋光面均勻涂抹厚約25~40μm的水溶性黑色油漆,晾干后,所有巖芯樣品均在室溫常壓條件下進(jìn)行三次熱導(dǎo)率測試,然后取三次測量結(jié)果的平均值作為熱導(dǎo)率測試值。

    表1研究區(qū)熱導(dǎo)率測試巖芯取樣井位

    Table 1 The core sampling boreholes of this study

    研究區(qū)鉆孔名稱北緯東經(jīng)取樣數(shù)(個)烏里雅斯太凹陷中康246° 27'31.3″118° 02'17.2″2中康1046° 24'44.0″117° 52'10.7″2中康2446° 21'18.4″117° 46'1.2″2中康17-1046° 25'15.6″117° 50'16.2″18中康2546° 20'27.7″117° 46'53.6″5中康2846° 24'55.6″117° 48'43.4″9中康1846° 25'17.9″117° 49'14.4″17白音查干凹陷查3042°04'14.6″108° 09'45.3″6錫4142° 06'45.9″108° 05'35.7″3錫242° 02'12.6″107° 46'47.0″5翁142° 05'58.9″108° 08'48.8″10達(dá)642° 05'51.5″108° 03'57.7″8烏142° 11'23.2″108° 30'14.6″11

    2.2 測試結(jié)果

    實(shí)驗(yàn)室測試結(jié)果(見電子版附表1)統(tǒng)計(jì)分析(圖2)表明,白音查干凹陷泥巖熱導(dǎo)率范圍為1.10~2.70W/(m·K),平均值為1.86±0.45W/(m·K);砂巖熱導(dǎo)率范圍為0.89~3.14W/(m·K),平均值為1.81±0.54W/(m·K);片巖熱導(dǎo)率范圍為3.11~4.91W/(m·K),平均值為4.07±0.78W/(m·K);玄武巖熱導(dǎo)率范圍1.23~1.40W/(m·K),平均值為1.29±0.09W/(m·K)。烏里雅斯太凹陷泥巖熱導(dǎo)率范圍為1.69~3.15W/(m·K),平均熱導(dǎo)率2.06±0.26W/(m·K);砂巖熱導(dǎo)率范圍為1.23~3.24W/(m·K),平均值為2.12±0.41W/(m·K);礫巖的熱導(dǎo)率范圍為1.32~2.86W/(m·K),平均值為2.31±0.42W/(m·K)。從熱導(dǎo)率數(shù)據(jù)來看,研究區(qū)內(nèi)片巖平均熱導(dǎo)率最高,玄武巖最低。在沉積巖中,泥巖、砂巖和礫巖三種巖性熱導(dǎo)率比較接近,礫巖最高,砂巖變化范圍最大,這可能與砂巖的自身結(jié)構(gòu)有關(guān)系。白音查干凹陷測試樣品主要分布于100~3500m之間,沉積巖熱導(dǎo)率隨深度增大有增加的趨勢;烏里雅斯太凹陷地區(qū)鉆孔測試樣品則主要分布于800~2600m,其沉積巖熱導(dǎo)率隨深度同樣有增大的趨勢(圖3)。白音查干凹陷巖石熱導(dǎo)率主要分布在1.00~2.50W/(m·K),烏里雅斯太凹陷則主要分布在1.50~2.50W/(m·K),將兩個地區(qū)巖石熱導(dǎo)率進(jìn)行疊加,其熱導(dǎo)率數(shù)據(jù)主要介于1.50~2.50W/(m·K)之間(圖2)。

    圖2 實(shí)測熱導(dǎo)率數(shù)據(jù)統(tǒng)計(jì)直方圖 (a)白音查干凹陷;(b)烏里雅斯太凹陷;(c)兩個地區(qū)疊加數(shù)據(jù)Fig.2 Histograms of the measured thermal conductivities (a) Baiyinchagan Sag; (b) Uliastai Sag; (c) Baiyinchagn Sag and Uliastai Sag combined

    圖3 白音查干凹陷(a)和烏里雅斯太凹陷(b)實(shí)測熱導(dǎo)率與深度的關(guān)系Fig.3 Measured thermal conductivity versus depth from Baiyinchagan Sag (a) and Uliastai Sag (b)

    3 熱導(dǎo)率校正

    本文熱導(dǎo)率數(shù)據(jù)均是在室溫常壓下對巖石干樣進(jìn)行測試得到,但是巖石在地下原位的溫度和壓力狀況明顯有別于實(shí)驗(yàn)室的條件,而且在原位高壓下巖石孔隙處于飽水狀態(tài)。因此,需要對干樣熱導(dǎo)率進(jìn)行飽水、溫度和壓力校正。我們分別對白音查干凹陷的3口油井和烏里雅斯太凹陷的4口油井進(jìn)行了測溫工作(Xuetal., 2019; Zhangetal., 2020),獲得了這兩個凹陷地溫隨深度變化的實(shí)際數(shù)據(jù),為這兩個地區(qū)的實(shí)測巖石熱導(dǎo)率進(jìn)行溫度和壓力校正提供了基礎(chǔ)數(shù)據(jù)支撐。

    3.1 校正方法

    3.1.1 飽水校正

    Hashin-Shtrikman模型中定義兩相介質(zhì)的各向同性均勻巖石的有效熱導(dǎo)率λeff有如下形式的上下限:

    (1)

    式中λf為介質(zhì)孔隙中流體的熱導(dǎo)率(W/(m·K));λm為基質(zhì)熱導(dǎo)率(W/(m·K));φ為孔隙度(%),取樣品實(shí)測值或者補(bǔ)償中子測井值計(jì)算得到的孔隙度值。對式(1)進(jìn)行一定的化簡后可以得到:

    (2)

    在室溫常壓測試下,由于干樣巖芯孔隙中飽含空氣,空氣熱導(dǎo)率λf為0.026W/(m·K),可視為隔熱相忽略不計(jì),那么巖石基質(zhì)熱導(dǎo)率與干樣熱導(dǎo)率具有以下關(guān)系:

    (3)

    將式(3)化簡得到基質(zhì)熱導(dǎo)率表達(dá)式:

    (4)

    將式(2)、(3)、(4)進(jìn)行替代并化簡,可得:

    (5)

    (6)

    式中,λdry為干巖樣實(shí)測熱導(dǎo)率值(W/(m·K));當(dāng)計(jì)算巖石飽水情況下的有效熱導(dǎo)率數(shù)值時,λf取值為0.607W/(m·K)。

    3.1.2 溫壓校正

    多孔巖石熱導(dǎo)率不僅受孔隙中所充填的流體影響,也隨所在地層溫度和壓力的不同而變化。Abdulagatovaetal. (2010)對油田中孔隙度為13%的飽和砂巖進(jìn)行了溫度和壓強(qiáng)實(shí)驗(yàn),發(fā)現(xiàn)巖石的熱導(dǎo)率隨溫度增加而減小,而隨著壓強(qiáng)的增加而增大,并推導(dǎo)出了砂巖巖石熱導(dǎo)率與溫度T和壓強(qiáng)P的經(jīng)驗(yàn)公式:

    (7)

    式中,λsf(T)為當(dāng)溫度為T時高壓下狀態(tài)下熱導(dǎo)率與常壓狀態(tài)下熱導(dǎo)率的差值(W/(m·K));P0為常數(shù),取值32MPa(Abdulagatovaetal., 2009);λ0, T為常壓下,溫度為T時的熱導(dǎo)率(W/(m·K));P為壓強(qiáng)(Pa)。

    (8)

    其中,ρ為地層密度(kg/m3);g為重力加速度(m/s2);h為巖石樣品原位深度(m)。

    在鉆孔溫度為15~250℃之間,壓力小于100MPa時,壓力每增加10MPa,巖石的熱導(dǎo)率增加1.3%~2.0%,在本研究中取值1.5%;溫度每升高100K,熱導(dǎo)率較常溫下降4.8%~6%,取平均值5.4%,那么λsf(T)和λT可簡化為:

    圖4 白音查干凹陷(a)和烏里雅斯太凹陷(b)飽水校正值與干樣熱導(dǎo)率對比Fig.4 Comparison of water saturation corrected and the measured thermal conductivities from Baiyinchagan Sag (a) and Uliastai Sag (b)

    (9)

    (10)

    式中,λm為室溫常壓下的熱導(dǎo)率實(shí)測值(W/(m·K));T0為常溫(K)。

    為了估算巖芯樣品的原位溫度,本文采用7口油井的實(shí)際測溫?cái)?shù)據(jù)(Xuetal., 2019; Zhangetal., 2020),利用線性回歸分析方法分別計(jì)算得到白音查干凹陷和烏里雅斯太凹陷的地溫梯度,其中白音查干凹陷為37.1℃/km,烏里雅斯太凹陷為39.4℃/km。在公式(8)中我們?nèi)〉貙悠骄芏葹?.5×103kg/m3,常溫取值288K,結(jié)合公式(8)、(9)和(10),把溫度壓強(qiáng)校正公式(7)化簡為式(11),校正后的熱導(dǎo)率λc為:

    (11)

    其中巖石原位深度h根據(jù)實(shí)際情況取值,(T-288)為巖芯樣品原位溫度相較于常溫288K的差值。

    3.2 校正結(jié)果

    根據(jù)補(bǔ)償中子測井值計(jì)算得到白音查干凹陷巖石孔隙度在2%~20%之間,對于烏里雅斯太凹陷,中原油田地質(zhì)資料顯示,烏里雅斯太凹陷北洼槽阿爾善組平均孔隙度為11.7%,騰格爾組平均孔隙度為14.4%,考慮到泥巖結(jié)構(gòu)相對致密,研究區(qū)的泥巖孔隙度均取值2.0%,烏里雅斯太凹陷阿爾善組和騰格爾組的砂巖和礫巖孔隙度取值分別為11.7%和14.4%。

    圖4和表2表明,經(jīng)飽水校正后的巖石熱導(dǎo)率均高于干樣實(shí)測熱導(dǎo)率,但是變化不是很大。烏里雅斯太凹陷泥巖熱導(dǎo)率校正后比校正前增大不到1.5%,而礫巖和砂巖熱導(dǎo)率的數(shù)值校正前后變化在10%以內(nèi)。白音查干凹陷泥巖熱導(dǎo)率校正幅度平均在2.5%左右, 砂巖則在15%以內(nèi)。泥巖熱

    表2巖石樣品熱導(dǎo)率飽水校正結(jié)果

    Table 2 Summary of thermal conductivities after water saturation correction

    研究區(qū)巖性熱導(dǎo)率范圍(W/(m·K))實(shí)測值校正值孔隙度(%)校正量(%)范圍平均值白音查干凹陷泥巖1.10~2.701.16~2.7421.3~5.32.4砂巖0.89~3.141.12~3.266~203.9~25.812.2烏里雅斯太凹陷泥巖1.69~3.151.71~3.1720~1.41.1砂巖1.23~3.241.38~3.4011.7、14.44.4~11.97.4礫巖1.32~2.861.47~3.0211.7、14.45.5~11.37

    圖5 白音查干凹陷(a)和烏里雅斯太凹陷(b)巖石樣品熱導(dǎo)率溫度和壓力綜合校正結(jié)果Fig.5 Measured versus corrected thermal conductivities of the core samples from Baiyinchagan Sag (a) and Uliastai Sag (b)

    導(dǎo)率校正前后變化不大(<5%),這是由于泥巖結(jié)構(gòu)致密且孔隙度較小(均值為2%)的緣故。

    本文中巖石熱導(dǎo)率的溫度和壓力校正是假設(shè)巖石已處于飽水情況下進(jìn)行的。白音查干凹陷泥巖和砂巖熱導(dǎo)率溫度和壓力校正值相對于干樣實(shí)測值平均增大4.5%和14.8%(圖5、表3)。烏里雅斯太凹陷泥巖綜合校正后熱導(dǎo)率平均增加4.3%,對于砂巖和礫巖而言,綜合校正后平均分別增加10.7%和10.2%(圖5、表3)。熱導(dǎo)率溫度和壓力校正值相對于飽水校正值校正前后變化不大,其中白音查干凹陷泥巖和砂巖平均增加2%,烏里雅斯太凹陷泥巖、砂巖和礫巖平均增加3%。校正幅度較小(<5%)是因?yàn)闇囟刃U蛪毫πU谝欢ǔ潭壬舷嗷サ窒?,且?guī)r芯埋藏深度較淺,地層壓力較小,溫度較低,故而校正幅度較小。

    經(jīng)過飽水、溫度和壓力校正后,兩個地區(qū)校正后的巖石熱導(dǎo)率相較于實(shí)驗(yàn)室條件下實(shí)測熱導(dǎo)率均有所增加。白音查干凹陷的干樣實(shí)測熱導(dǎo)率范圍為0.89~4.91W/(m·K),平均值為2.06±0.86W/(m·K),校正后熱導(dǎo)率范圍為1.12~4.93W/(m·K),平均值為2.18±0.85W/(m·K),烏里雅斯太凹陷實(shí)測熱導(dǎo)率范圍為1.23~3.24W/(m·K),平均值為2.12±0.37W/(m·K),校正后熱導(dǎo)率范圍為1.42~3.51W/(m·K),平均值為2.29±0.39W/(m·K)。其中,砂巖和礫巖校正前后變化值較大,各巖性熱導(dǎo)率校正情況見表3。白音查干凹陷校正后熱導(dǎo)率主要分布在1.00~2.50W/(m·K)之間,烏里雅斯太凹陷則主要分布在1.50~2.50W/(m·K)之間,將兩個地區(qū)熱導(dǎo)率數(shù)據(jù)疊加,校正后的巖石熱導(dǎo)率主要分布在1.50~2.50W/(m·K)之間(圖6)。此外,校正后巖石熱導(dǎo)率和實(shí)測熱導(dǎo)率整體上均呈現(xiàn)隨深度增大的趨勢(圖7)。

    表3巖石樣品熱導(dǎo)率溫度和壓力校正結(jié)果

    Table 3 Summary of core sample thermal conductivity corrections for the combined temperature-pressure effects

    研究區(qū)巖性熱導(dǎo)率范圍(W/(m·K))平均校正量(%)干樣實(shí)測值飽水校正值溫壓校正值溫壓校正綜合校正白音查干凹陷泥巖1.10~2.701.16~2.741.15~2.8124.5砂巖0.89~3.141.12~3.261.12~3.27214.8烏里雅斯太凹陷泥巖1.69~3.151.71~3.171.77~3.2734.3砂巖1.23~3.241.38~3.401.42~3.51310.7礫巖1.32~2.861.47~3.021.51~3.11310.2

    4 地層熱導(dǎo)率

    以地層為單位進(jìn)行統(tǒng)計(jì),利用校正后的熱導(dǎo)率數(shù)據(jù)計(jì)算各組地層各巖性熱導(dǎo)率的算術(shù)平均值,加之統(tǒng)計(jì)得到的不同地層巖性厚度比,利用加權(quán)平均公式可以計(jì)算不同地層的熱導(dǎo)率λ,計(jì)算公式如下:

    λ=λ1P1+λ2P2+…+λnPn

    (12)

    式中λ1、λ2、λn為各巖性的平均熱導(dǎo)率(W/(m·K));P1、P2、…、Pn為各巖性在地層中的厚度比(%)。

    白音查干凹陷主要巖性是砂巖和泥巖,還存在少量礫巖、變質(zhì)巖和火成巖,由于后三種巖性占比小(小于1%),在計(jì)算地層熱導(dǎo)率時可忽略不計(jì)。本文中選取泥巖和砂巖樣品,共計(jì)36塊,建立白音查干凹陷地區(qū)的地層熱導(dǎo)率(表4)。其中賽罕塔拉組、都紅木組、騰格爾組和阿爾善組地層熱導(dǎo)率分別為1.24W/(m·K)、2.22W/(m·K)、2.06W/(m·K)和2.10W/(m·K),整個白堊系地層熱導(dǎo)率為2.00W/(m·K)。同樣利用干樣實(shí)測值計(jì)算了各個地層的熱導(dǎo)率,其結(jié)果顯示白堊系地層熱導(dǎo)率為1.84W/(m·K),利用校正后的巖石熱導(dǎo)率計(jì)算得到的白堊系地層熱導(dǎo)率相較于利用實(shí)測巖石熱導(dǎo)率計(jì)算得到的白堊系地層熱導(dǎo)率高8.7%。

    同樣利用校正后的巖石熱導(dǎo)率,建立了烏里雅斯太凹陷的地層熱導(dǎo)率(表5)。烏里雅斯太凹陷各個地層熱導(dǎo)率均大于2.00W/(m·K),白堊系地層熱導(dǎo)率為2.17W/(m·K),而侏羅-二疊系的地層熱導(dǎo)率為2.84W/(m·K),這可能與該地層埋藏較深,巖石結(jié)構(gòu)相對致密有關(guān)。此外,利用干樣實(shí)測值計(jì)算了該地區(qū)的地層熱導(dǎo)率,其結(jié)果相較于利用校正后的熱導(dǎo)率計(jì)算得到的白堊系地層熱導(dǎo)率低6.8%。利用校正后的地層熱導(dǎo)率值、地溫梯度、地層生熱率等熱物性參數(shù)(Xuetal., 2019; Zhangetal., 2020)計(jì)算得到白音查干凹陷和烏里雅斯太凹陷平均大地?zé)崃髦禐?4.2mW/m2和85.5mW/m2,明顯大于中國大陸地區(qū)平均熱流值61.5mW/m2(姜光政等, 2016)。

    圖6 熱導(dǎo)率校正值統(tǒng)計(jì)直方圖 (a)白音查干凹陷;(b)烏里雅斯太凹陷;(c)白音查干和烏里雅斯太凹陷熱導(dǎo)率數(shù)據(jù)疊加Fig.6 Histograms of the corrected thermal conductivity values (a) Baiyinchagan Sag; (b) Uliastai Sag; (c) Baiyinchagan Sag and Uliastai Sag combined

    表4白音查干凹陷地層巖石熱導(dǎo)率

    Table 4 Thermal conductivities of rock formations in the Baiyinchagan Sag

    地層巖性樣品數(shù)熱導(dǎo)率范圍(W/(m·K))平均熱導(dǎo)率(W/(m·K))地層厚度占比(%)加權(quán)平均值(W/(m·K))地層熱導(dǎo)率(W/(m·K))白堊系賽罕塔拉都紅木騰格爾阿爾善泥巖21.15~1.291.22±0.1067.2砂巖41.12~1.751.29±0.3132.8泥巖51.51~2.201.85±0.2753.3砂巖22.04~3.272.65±0.8746.7泥巖41.45~2.791.98±0.5862.8砂巖111.80~2.592.19±0.2037.2泥巖61.94~2.812.21±0.3145.2砂巖21.52~2.451.99±0.6654.81.242.222.062.102.00

    表5烏里雅斯太凹陷地層巖石熱導(dǎo)率

    Table 5 Thermal conductivities of rock formations in the Uliastai Sag

    地層巖性樣品數(shù)熱導(dǎo)率范圍(W/(m·K))平均熱導(dǎo)率(W/(m·K))巖性含量(%)加權(quán)熱導(dǎo)率(W/(m·K))地層熱導(dǎo)率(W/(m·K))白堊系騰一上騰一中騰一下阿爾善礫巖41.51~3.112.53±0.701.1泥巖111.94~2.332.13±0.1373.3砂巖111.42~2.562.11±0.6525.6礫巖62.26~2.852.61±0.223.0泥巖151.81~2.292.09±0.1666.4砂巖161.56~3.152.41±0.4130.6礫巖12.132.131.0泥巖22.03~2.262.14±0.1664.0砂巖82.08~3.092.33±0.3235.0泥巖31.77~2.061.94±0.1557.2砂巖62.13~3.262.47±0.4242.82.132.202.212.242.17侏羅系、二疊系泥巖32.31~3.272.72±0.5057.1砂巖22.52~3.513.01±0.7042.92.842.84

    表6研究區(qū)及周邊盆地白堊系巖石熱導(dǎo)率

    Table 6 Thermal conductivities of the Cretaceous rocks in the study area and surrounding basins

    地區(qū)巖性樣品數(shù)熱導(dǎo)率范圍(W/(m·K))均值±標(biāo)準(zhǔn)差(W/(m·K))數(shù)據(jù)來源白音查干凹陷烏里雅斯太凹陷查干凹陷塔里木盆地準(zhǔn)噶爾盆地泥巖151.10~2.701.86±0.45砂巖200.89~3.141.81±0.54泥巖311.69~3.152.06±0.26砂巖311.23~3.242.12±0.41礫巖111.32~2.862.31±0.42泥巖351.46~2.562.07±0.28砂巖261.14~2.922.13±0.41礫巖191.75~3.452.53±0.52泥巖181.519~2.4241.908±0.264砂巖450.523~2.9581.565±0.662礫巖111.336~3.9292.521±0.748泥巖131.665~2.4071.996±0.234砂巖111.302~2.1691.700±0.283實(shí)測實(shí)測+收集左銀輝等, 2013馮昌格等, 2009饒松等, 2013

    5 討論

    5.1 研究區(qū)與周邊地區(qū)熱導(dǎo)率差異

    在二連盆地白音查干凹陷和烏里雅斯太凹陷采集測量的各類巖石樣品中均是砂巖熱導(dǎo)率變化范圍最大,其中白音查干凹陷砂巖熱導(dǎo)率范圍為0.89~3.14W/(m·K),烏里雅斯太凹陷為1.24~3.24W/(m·K),與典型砂巖熱導(dǎo)率范圍(2.60~4.80W/(m·K))(Norden and F?rster, 2006; Liuetal., 2011)有所差別。兩個研究區(qū)砂巖具有變化范圍較廣的熱導(dǎo)率值,在其他盆地研究中也有所體現(xiàn)。馮昌格等(2009)對塔里木盆地940多塊實(shí)測巖芯熱導(dǎo)率進(jìn)行了整理,其結(jié)果顯示塔里木盆地從寒武紀(jì)地層到新近紀(jì)地層的巖石熱導(dǎo)率中,砂巖的變化范圍最大(表6)。對比白音查干凹陷都紅木組和烏里雅斯太凹陷騰格爾組二段,白音查干凹陷騰格爾組和烏里雅斯太凹陷騰格爾中段,可以發(fā)現(xiàn)其砂、泥巖熱導(dǎo)率數(shù)值范圍相差不大,這是由于兩個地區(qū)在同一地質(zhì)時期沉積環(huán)境類似。白音查干凹陷和烏里雅斯太凹陷地區(qū)砂、泥巖間的巖石熱導(dǎo)率平均值相差不大,因此不能簡單地以巖石的熱物性數(shù)據(jù)來區(qū)別巖石類型(李香蘭等, 2015),必須在研究區(qū)內(nèi)采集相當(dāng)數(shù)量的有代表性的巖石進(jìn)行實(shí)測。

    表7我國西北部幾個主要沉積盆地實(shí)測巖石熱導(dǎo)率

    Table 7 The measured thermal conductivities of several major sedimentary basins in Northwest China

    盆地地層 柴達(dá)木盆地準(zhǔn)噶爾盆地塔里木盆地海拉爾盆地銀額盆地二連盆地白堊系侏羅系2.111.461.671.46~1.702.17~2.352.00~2.172.191.79——2.84數(shù)據(jù)來源邱楠生, 2002崔軍平等, 2007左銀輝等, 2013本文工作

    圖7 白音查干凹陷(a)和烏里雅斯太凹陷(b)熱導(dǎo)率校正值與深度的關(guān)系Fig.7 Corrected thermal conductivities versus depths from Baiyinchagan Sag (a) and Uliastai Sag (b)

    圖8 白音查干凹陷砂巖熱導(dǎo)率與孔隙度關(guān)系Fig.8 The relationship between the thermal conductivity of sandstone and porosity in the Baiyinchagan Sag

    本文基于百余塊巖芯樣品的實(shí)測結(jié)果得到白音查干凹陷和烏里雅斯太凹陷白堊系地層熱導(dǎo)率分別為2.00W/(m·K)和2.17W/(m·K),明顯大于肖偉等(2004)報道的二連盆地?zé)釋?dǎo)率值1.80W/(m·K),但與位于二連盆地西側(cè)的銀額盆地查干凹陷地層熱導(dǎo)率范圍2.16~2.35W/(m·K)(左銀輝等, 2013)較為接近,此外,也與柴達(dá)木盆地白堊系地層熱導(dǎo)率2.11W/(m·K)(邱楠生, 2002)非常接近。但二連盆地白堊系熱導(dǎo)率高于海拉爾盆地、準(zhǔn)噶爾盆地和塔里木盆地白堊系地層熱導(dǎo)率(表7),可能是由于研究區(qū)和這幾個盆地經(jīng)歷了不同的沉積成巖作用以及物理—地球化學(xué)演化過程。

    5.2 熱導(dǎo)率影響因素

    巖石熱導(dǎo)率大小受巖石組成礦物成分、孔隙度、含水飽和度、溫度、壓力等因素影響(Abdulagatovaetal., 2009; Alishaevetal., 2012),但最主要的影響因素是礦物成分、孔隙度和含水飽和度(Clauser and Huenges, 1995; Alishaevetal., 2012; Nagaraju and Roy, 2014; Chenetal., 2017)。本文實(shí)測數(shù)據(jù)表明,研究區(qū)內(nèi)砂巖熱導(dǎo)率變化范圍較大,這應(yīng)該是因?yàn)檠芯繀^(qū)砂巖孔隙度變化范圍大且砂巖礦物組成不一樣的緣故。白音查干凹陷砂巖孔隙度變化范圍為5%~20%。烏里雅斯太凹陷因?yàn)槲镌磪^(qū)母巖性質(zhì)的變化,巖石類型以長石巖屑砂巖、巖屑長石砂巖、砂礫巖為主,少量長石砂巖,巖屑砂巖則更少,南、北次洼也有所差別,北次洼石英、長石含量比南次洼高,南次洼無長石砂巖類,北次洼缺少巖屑砂巖類(趙澄林等, 1996)。在上述巖石組成礦物中石英熱導(dǎo)率為7.8W/(m·K),長石為2.3W/(m·K)(Horai, 1971)。在孔隙度較小的巖石中,巖石熱導(dǎo)率的大小與石英含量存在正相關(guān)關(guān)系(Vosteen and Schellschmidt, 2003)。針對各個巖性礦物組成對巖石熱導(dǎo)率的影響,Jorandetal. (2015)進(jìn)行過具體分析。

    在地質(zhì)體原位情況下,大多數(shù)的巖石都具有孔隙,孔隙中充填有液體或氣體,這實(shí)際上是一種多相介質(zhì)系統(tǒng)。干燥和飽水巖石為二相介質(zhì),含水而非飽和的巖石為三相介質(zhì)。因充填于巖石孔隙中的空氣(0.026W/(m·K))和水(0.607W/(m·K))的熱導(dǎo)率值比巖石的固體成分低,故相同巖性巖石熱導(dǎo)率一般會隨著孔隙度的增大而降低。如圖8所示,白音查干凹陷砂巖熱導(dǎo)率整體上是隨著孔隙度的增加而降低,而圖中相同孔隙度的砂巖熱導(dǎo)率出現(xiàn)了不同的值,這應(yīng)該與巖石的礦物組成有關(guān)。此外,空氣的熱導(dǎo)率比水低(水的熱導(dǎo)率約為空氣熱導(dǎo)率的23倍),因而干燥巖樣熱導(dǎo)率值小于飽水巖樣熱導(dǎo)率值。當(dāng)巖樣孔隙度較大,含水而非飽和時,其熱導(dǎo)率則隨著含水率的增加而增加(Nagaraju and Roy, 2014)。目前,比較常用的飽水校正方法是幾何均值模型(Geometric mean),但是近來研究發(fā)現(xiàn)當(dāng)巖石孔隙度大于6%時,該方法校正后的熱導(dǎo)率數(shù)值明顯偏大(Nagaraju and Roy, 2014; Guoetal., 2017),故本文選取有效熱導(dǎo)率理論對研究區(qū)的實(shí)測熱導(dǎo)率進(jìn)行飽水校正。

    此外,巖石熱導(dǎo)率會隨著地層壓力增大而增加,這與巖石壓實(shí)和和膠結(jié)作用有關(guān)。隨著巖石埋藏深度加大,地層年代的增加,在上覆地層壓力的作用下,碎屑顆粒緊密堆積,沉積物體積變小,密度變大,同時隨著膠結(jié)作用的進(jìn)行,孔隙度隨著深度增加而逐漸減小(劉偉等, 2002)。但是隨著壓力的進(jìn)一步增加,巖石中的裂縫和孔隙會逐漸閉合,以致熱導(dǎo)率趨于一個恒定值(Walsh and Decker, 1966; G?rgülüetal., 2008)。對于溫度而言,溫度升高會使礦物顆粒間產(chǎn)生熱阻抗,導(dǎo)致熱傳導(dǎo)率降低。溫度較低時,傳導(dǎo)效應(yīng)大于輻射效應(yīng);但隨溫度升高,輻射效應(yīng)會逐漸增大,甚至超過傳導(dǎo)效應(yīng),但是增大至一定程度后就不再增加(Clauser and Huenges, 1995)。一般規(guī)律是隨著溫度的升高,巖石熱導(dǎo)率會降低(Jhaetal., 2016)。本文中,巖芯樣品所在深度范圍為117.8~3159.5m之間,對應(yīng)的壓強(qiáng)和溫度范圍分別為2.89~77.41MPa和13.5~118.6℃,其泥巖熱導(dǎo)率溫度壓力綜合校正幅度在2%~3.1%之間,砂巖和礫巖校正幅度在1.7%~3.0%之間。校正后巖石熱導(dǎo)率增加較小,原因在于在采集的巖芯樣品采樣深度范圍內(nèi),巖石熱導(dǎo)率隨壓力增大有所增加而隨溫度升高略有降低,溫度效應(yīng)部分抵消了壓力的影響。

    6 結(jié)論

    基于二連盆地白音查干凹陷和烏里雅斯太凹陷巖芯樣品的巖石熱導(dǎo)率測量與統(tǒng)計(jì)分析,得到以下主要結(jié)論:

    (1)二連盆地129塊鉆孔巖芯樣品實(shí)驗(yàn)室條件下實(shí)測熱導(dǎo)率值范圍0.89~4.91W/(m·K),樣品熱導(dǎo)率隨深度增加有所增加,整體呈現(xiàn)正態(tài)分布的特點(diǎn),主要分布在1.50~2.50W/(m·K)之間,平均值為2.06±0.86W/(m·K);

    (2)本文測試的樣品包括泥巖、砂巖、礫巖、頁巖和玄武巖等五類巖石,孔隙度在2%~20%之間,對于孔隙度較小的泥巖,飽水對巖石熱導(dǎo)率的影響較小,可以忽略,但對于孔隙度較大的砂、礫巖,熱導(dǎo)率飽水校正量可高達(dá)20%;

    (3)本文采集樣品深度在117.8~3159.5m范圍內(nèi),對應(yīng)的溫度和壓力范圍分別為13.5~118.6℃和2.89~77.41MPa,溫度和壓力對巖石熱導(dǎo)率的影響在一定程度上相互抵消,因而綜合影響不明顯(<5%);

    (4)經(jīng)過對實(shí)驗(yàn)室條件下干樣測試結(jié)果進(jìn)行了原位條件下的飽水和溫度壓力綜合校正,研究區(qū)不同巖性的巖石樣品平均熱導(dǎo)率值分別為泥巖2.08±0.36W/(m·K)、砂巖2.28±0.50W/(m·K)、礫巖2.53±0.44W/(m·K)、頁巖4.16±0.76W/(m·K)、玄武巖1.33±0.09W/(m·K);

    (5)根據(jù)地層厚度加權(quán)平均得到的白音查干凹陷和烏里雅斯太凹陷白堊系地層熱導(dǎo)率分別為2.00W/(m·K)和2.17W/(m·K),與銀額盆地和柴達(dá)木盆地的地層熱導(dǎo)率相近,比海拉爾盆地、準(zhǔn)噶爾盆地和塔里木盆地的地層熱導(dǎo)率高。結(jié)合鉆孔測溫?cái)?shù)據(jù),計(jì)算得到二連盆地大地?zé)崃髦到橛?4~85mW/m2,明顯高于此前文獻(xiàn)估算的43~66mW/m2和中國大陸地區(qū)平均熱流值61.5mW/m2。

    致謝感謝翟明國院士和匿名審稿人對本文初稿的關(guān)注和修改意見。西安交通大學(xué)地?zé)崤c環(huán)境研究實(shí)驗(yàn)室劉植和李毅在本文寫作過程中提出了建設(shè)性意見,在此表示感謝!

    猜你喜歡
    查干白音熱導(dǎo)率
    白音和新中國劇社的往事
    公民與法治(2022年3期)2022-07-29 00:57:30
    喚醒芬芳
    空位缺陷對單層石墨烯導(dǎo)熱特性影響的分子動力學(xué)
    連續(xù)碳纖維鋁基復(fù)合材料橫向等效熱導(dǎo)率的模擬分析
    Si3N4/BN復(fù)合陶瓷熱導(dǎo)率及其有限元分析
    中國女?dāng)z影家聚焦一代天驕查干蘇魯克大典采風(fēng)活動花絮及作品選登(下)
    旅游縱覽(2019年8期)2019-08-14 18:29:26
    查干凹陷蘇紅圖組火山巖儲集特征及主控因素
    查干凹陷火山巖裂縫發(fā)育特征及其影響因素
    錄井工程(2017年3期)2018-01-22 08:40:23
    金屬熱導(dǎo)率的第一性原理計(jì)算方法在鋁中的應(yīng)用
    星星
    中國詩歌(2015年12期)2015-11-17 19:30:54
    人妻夜夜爽99麻豆av| 日韩一本色道免费dvd| 成年免费大片在线观看| 亚洲av电影在线观看一区二区三区 | 极品少妇高潮喷水抽搐| 97热精品久久久久久| 久久99蜜桃精品久久| 午夜日本视频在线| 2021天堂中文幕一二区在线观| 亚洲电影在线观看av| 在线看a的网站| 久久久亚洲精品成人影院| 欧美日韩在线观看h| 亚洲av成人精品一二三区| 欧美精品国产亚洲| 国产免费一区二区三区四区乱码| 人妻系列 视频| 成人二区视频| 黄色欧美视频在线观看| 国产成人freesex在线| 午夜福利在线观看免费完整高清在| av福利片在线观看| 熟女电影av网| 内射极品少妇av片p| 男人爽女人下面视频在线观看| 一级毛片我不卡| 热99国产精品久久久久久7| 99热全是精品| 亚洲va在线va天堂va国产| 国产精品三级大全| 亚洲婷婷狠狠爱综合网| 卡戴珊不雅视频在线播放| 亚洲精品日本国产第一区| 啦啦啦啦在线视频资源| 国产成人精品一,二区| 人人妻人人爽人人添夜夜欢视频 | 丝袜脚勾引网站| 久久精品国产亚洲网站| 99热这里只有是精品50| 伦理电影大哥的女人| 色网站视频免费| 女人被狂操c到高潮| 中文字幕制服av| 久久久久久久国产电影| 国产欧美日韩精品一区二区| 十八禁网站网址无遮挡 | 高清av免费在线| 人妻 亚洲 视频| 熟女人妻精品中文字幕| 日本三级黄在线观看| 亚洲欧美成人精品一区二区| 制服丝袜香蕉在线| 精品一区二区免费观看| 欧美成人a在线观看| 最近2019中文字幕mv第一页| 日韩视频在线欧美| 一级毛片黄色毛片免费观看视频| 精华霜和精华液先用哪个| 少妇的逼好多水| av国产久精品久网站免费入址| 极品教师在线视频| 黄色配什么色好看| av网站免费在线观看视频| 国产欧美另类精品又又久久亚洲欧美| 国产精品嫩草影院av在线观看| 老师上课跳d突然被开到最大视频| 成人免费观看视频高清| 亚洲精品,欧美精品| 亚洲自偷自拍三级| 插阴视频在线观看视频| 免费黄网站久久成人精品| 午夜视频国产福利| 日韩人妻高清精品专区| 国产av不卡久久| 午夜福利视频精品| 内地一区二区视频在线| 国产成人精品福利久久| 国产真实伦视频高清在线观看| 精品久久久久久电影网| 国产精品女同一区二区软件| 亚洲av免费高清在线观看| 亚洲国产成人一精品久久久| 欧美人与善性xxx| 全区人妻精品视频| 26uuu在线亚洲综合色| 欧美日韩精品成人综合77777| 男女边吃奶边做爰视频| 在线 av 中文字幕| 最近2019中文字幕mv第一页| 婷婷色综合www| 69av精品久久久久久| 精品久久国产蜜桃| 亚洲精品久久午夜乱码| 插阴视频在线观看视频| 在线免费十八禁| 最近中文字幕高清免费大全6| 精品久久国产蜜桃| 久久精品久久精品一区二区三区| 18禁在线无遮挡免费观看视频| 亚洲精品视频女| 精品熟女少妇av免费看| 男的添女的下面高潮视频| 99热网站在线观看| 少妇丰满av| 哪个播放器可以免费观看大片| 中文欧美无线码| 国产精品久久久久久久电影| 欧美精品一区二区大全| 日韩av在线免费看完整版不卡| 美女脱内裤让男人舔精品视频| 十八禁网站网址无遮挡 | 日本av手机在线免费观看| 国产亚洲精品久久久com| 91在线精品国自产拍蜜月| 日韩大片免费观看网站| 嫩草影院新地址| 最近2019中文字幕mv第一页| 亚洲人成网站高清观看| 岛国毛片在线播放| 午夜福利在线观看免费完整高清在| 99re6热这里在线精品视频| 成人欧美大片| 国产熟女欧美一区二区| 亚洲精品亚洲一区二区| 免费观看a级毛片全部| 亚洲精品aⅴ在线观看| 国产精品偷伦视频观看了| 少妇人妻久久综合中文| 美女内射精品一级片tv| 亚洲国产精品成人久久小说| 在线观看一区二区三区激情| 一级毛片aaaaaa免费看小| 国产午夜精品一二区理论片| 免费看日本二区| 麻豆精品久久久久久蜜桃| 中文字幕久久专区| 国产老妇女一区| 亚洲精品视频女| av在线观看视频网站免费| 国产精品不卡视频一区二区| 午夜免费男女啪啪视频观看| 免费观看av网站的网址| 成人免费观看视频高清| 国产精品国产av在线观看| 国产免费福利视频在线观看| 丰满人妻一区二区三区视频av| 夫妻午夜视频| 国产淫片久久久久久久久| 蜜桃亚洲精品一区二区三区| 干丝袜人妻中文字幕| 精品人妻熟女av久视频| 国内精品宾馆在线| 国产高清国产精品国产三级 | 亚洲伊人久久精品综合| 身体一侧抽搐| 亚洲精品视频女| 亚洲av一区综合| 2021少妇久久久久久久久久久| 青春草视频在线免费观看| 欧美亚洲 丝袜 人妻 在线| 精品少妇久久久久久888优播| 国产精品国产三级国产专区5o| 91在线精品国自产拍蜜月| 69人妻影院| 国产高潮美女av| 国产成人免费观看mmmm| 亚洲国产av新网站| 久久久久久久午夜电影| 99久国产av精品国产电影| 国产色爽女视频免费观看| 久久女婷五月综合色啪小说 | 伦精品一区二区三区| av.在线天堂| 精品人妻偷拍中文字幕| 亚洲色图av天堂| 伊人久久精品亚洲午夜| av福利片在线观看| 爱豆传媒免费全集在线观看| 欧美日本视频| 亚洲av中文av极速乱| av播播在线观看一区| 99re6热这里在线精品视频| 免费看av在线观看网站| 国产成年人精品一区二区| 岛国毛片在线播放| 嫩草影院入口| 一区二区三区精品91| 日韩在线高清观看一区二区三区| 一本久久精品| 深爱激情五月婷婷| 国产精品爽爽va在线观看网站| 日韩人妻高清精品专区| 亚洲电影在线观看av| 国产探花极品一区二区| 亚洲人成网站在线观看播放| 色网站视频免费| 国产精品av视频在线免费观看| 黄色怎么调成土黄色| 国产精品蜜桃在线观看| 久久热精品热| 美女内射精品一级片tv| 99久久精品一区二区三区| 亚洲性久久影院| 成人午夜精彩视频在线观看| 别揉我奶头 嗯啊视频| 亚洲国产色片| 最近2019中文字幕mv第一页| 久久久久精品性色| 校园人妻丝袜中文字幕| 七月丁香在线播放| 欧美精品国产亚洲| 精品午夜福利在线看| 免费电影在线观看免费观看| 亚洲精品日韩在线中文字幕| 在线a可以看的网站| 午夜免费观看性视频| 久久精品国产鲁丝片午夜精品| 狂野欧美激情性xxxx在线观看| 国产精品国产av在线观看| 国产女主播在线喷水免费视频网站| 久久久久久国产a免费观看| 免费黄色在线免费观看| 久久久久久久久久久免费av| 国国产精品蜜臀av免费| 如何舔出高潮| 男女边摸边吃奶| 国国产精品蜜臀av免费| 亚洲天堂国产精品一区在线| 中文欧美无线码| 性色avwww在线观看| 成人欧美大片| 久久久成人免费电影| 国产精品久久久久久精品电影小说 | 免费大片18禁| 久久久a久久爽久久v久久| 亚洲精品亚洲一区二区| 亚洲综合精品二区| 国产成人一区二区在线| 亚洲自偷自拍三级| 成人亚洲精品av一区二区| 天堂中文最新版在线下载 | 又黄又爽又刺激的免费视频.| 精品人妻偷拍中文字幕| 日韩免费高清中文字幕av| 中国美白少妇内射xxxbb| 九九久久精品国产亚洲av麻豆| av女优亚洲男人天堂| 亚洲,欧美,日韩| 亚洲综合色惰| 欧美成人精品欧美一级黄| 肉色欧美久久久久久久蜜桃 | www.av在线官网国产| 午夜激情久久久久久久| 丝瓜视频免费看黄片| 日韩人妻高清精品专区| 久久久久九九精品影院| 九草在线视频观看| 午夜激情久久久久久久| 极品少妇高潮喷水抽搐| 久久国产乱子免费精品| 狂野欧美激情性bbbbbb| 一级毛片黄色毛片免费观看视频| 女人十人毛片免费观看3o分钟| 精品久久久噜噜| 97人妻精品一区二区三区麻豆| 免费观看性生交大片5| 特大巨黑吊av在线直播| 亚洲精品视频女| 好男人视频免费观看在线| 亚洲精华国产精华液的使用体验| 在线观看三级黄色| 国产亚洲av片在线观看秒播厂| 国内揄拍国产精品人妻在线| 亚洲国产精品成人综合色| av黄色大香蕉| 亚洲天堂av无毛| 成人无遮挡网站| 欧美潮喷喷水| 免费电影在线观看免费观看| 欧美日韩亚洲高清精品| 日本黄大片高清| 精品久久国产蜜桃| 精品人妻一区二区三区麻豆| 亚洲欧美精品专区久久| 亚洲国产精品成人久久小说| 国产高清不卡午夜福利| 综合色丁香网| 人妻 亚洲 视频| 少妇被粗大猛烈的视频| 少妇熟女欧美另类| 国产午夜福利久久久久久| 久久综合国产亚洲精品| 一级毛片aaaaaa免费看小| 91aial.com中文字幕在线观看| 性色avwww在线观看| 80岁老熟妇乱子伦牲交| 国产精品国产三级国产av玫瑰| 麻豆精品久久久久久蜜桃| 插逼视频在线观看| av国产精品久久久久影院| 夜夜爽夜夜爽视频| 免费电影在线观看免费观看| 国产久久久一区二区三区| 国产极品天堂在线| 男女边吃奶边做爰视频| 成人美女网站在线观看视频| 亚洲美女视频黄频| 亚洲精品中文字幕在线视频 | 欧美潮喷喷水| 国产一区亚洲一区在线观看| av黄色大香蕉| 久久99热这里只有精品18| 女人被狂操c到高潮| 国内少妇人妻偷人精品xxx网站| 一区二区三区免费毛片| 成人亚洲精品一区在线观看 | 久久午夜福利片| 十八禁网站网址无遮挡 | 永久网站在线| 久久久久久国产a免费观看| 欧美xxxx性猛交bbbb| 女人被狂操c到高潮| 国产欧美日韩精品一区二区| 三级国产精品欧美在线观看| 亚洲av福利一区| 毛片一级片免费看久久久久| 91精品国产九色| 色网站视频免费| 一级a做视频免费观看| 人妻系列 视频| 在线精品无人区一区二区三 | 亚洲国产av新网站| 国产精品嫩草影院av在线观看| 看黄色毛片网站| eeuss影院久久| 97热精品久久久久久| 夜夜爽夜夜爽视频| av在线观看视频网站免费| 久久6这里有精品| 黄色欧美视频在线观看| 日日啪夜夜撸| 毛片一级片免费看久久久久| 免费少妇av软件| 色网站视频免费| 亚洲国产av新网站| 97热精品久久久久久| 日日啪夜夜爽| 国产熟女欧美一区二区| 国产成人福利小说| 干丝袜人妻中文字幕| 五月开心婷婷网| 精品午夜福利在线看| 青春草国产在线视频| 久久鲁丝午夜福利片| 亚洲精品色激情综合| 最后的刺客免费高清国语| 日韩 亚洲 欧美在线| 日韩一区二区视频免费看| 人妻夜夜爽99麻豆av| 亚洲av一区综合| 亚洲精品乱码久久久v下载方式| 特级一级黄色大片| 国产精品国产三级国产av玫瑰| 久久97久久精品| 精品人妻一区二区三区麻豆| 免费观看av网站的网址| 日韩大片免费观看网站| 狂野欧美激情性bbbbbb| 岛国毛片在线播放| 777米奇影视久久| 国产乱来视频区| 色网站视频免费| 一级毛片电影观看| a级毛色黄片| 国产精品.久久久| 2018国产大陆天天弄谢| 天堂俺去俺来也www色官网| 欧美最新免费一区二区三区| 精品久久久久久久久av| 欧美丝袜亚洲另类| 两个人的视频大全免费| 国产精品久久久久久精品电影小说 | 欧美xxⅹ黑人| 中文乱码字字幕精品一区二区三区| 亚洲国产高清在线一区二区三| 成人无遮挡网站| 人人妻人人澡人人爽人人夜夜| 亚洲欧美日韩卡通动漫| 免费av不卡在线播放| 97热精品久久久久久| 欧美区成人在线视频| 欧美成人一区二区免费高清观看| 国产成人精品福利久久| 精品少妇久久久久久888优播| 日本熟妇午夜| 91精品伊人久久大香线蕉| 国产精品一区二区三区四区免费观看| 日韩成人av中文字幕在线观看| 国产成人午夜福利电影在线观看| 青青草视频在线视频观看| 久热这里只有精品99| 久久精品综合一区二区三区| 免费av观看视频| 听说在线观看完整版免费高清| 日日摸夜夜添夜夜添av毛片| 少妇裸体淫交视频免费看高清| 天天一区二区日本电影三级| 欧美zozozo另类| 伦理电影大哥的女人| 在线观看免费高清a一片| 欧美激情久久久久久爽电影| 成年女人在线观看亚洲视频 | 欧美变态另类bdsm刘玥| 精品一区在线观看国产| 色5月婷婷丁香| 青青草视频在线视频观看| 国产成人一区二区在线| 国产免费视频播放在线视频| 久久久久久久久久人人人人人人| 好男人在线观看高清免费视频| 国产免费福利视频在线观看| 最近2019中文字幕mv第一页| av播播在线观看一区| 嘟嘟电影网在线观看| 国内揄拍国产精品人妻在线| 亚洲国产成人一精品久久久| 国内少妇人妻偷人精品xxx网站| 日韩欧美精品免费久久| 人妻一区二区av| 国产一区亚洲一区在线观看| 麻豆成人午夜福利视频| 亚洲性久久影院| 亚洲国产精品专区欧美| 一级黄片播放器| 天堂网av新在线| 国产黄片美女视频| 精品人妻视频免费看| 一级毛片 在线播放| 熟女人妻精品中文字幕| 国产伦精品一区二区三区四那| 日本色播在线视频| 久久久久网色| 国产精品.久久久| 黑人高潮一二区| 99久久人妻综合| 91午夜精品亚洲一区二区三区| 日本猛色少妇xxxxx猛交久久| 特大巨黑吊av在线直播| 一区二区av电影网| 热99国产精品久久久久久7| 国产精品不卡视频一区二区| xxx大片免费视频| 天堂俺去俺来也www色官网| 女人久久www免费人成看片| 九色成人免费人妻av| 99热这里只有是精品在线观看| 夜夜看夜夜爽夜夜摸| 在线观看国产h片| 老女人水多毛片| 亚洲激情五月婷婷啪啪| 久久久色成人| 久久女婷五月综合色啪小说 | 一个人观看的视频www高清免费观看| 免费少妇av软件| 亚洲成色77777| 丝袜美腿在线中文| 亚洲欧美一区二区三区黑人 | 观看美女的网站| www.av在线官网国产| 一级毛片我不卡| 大话2 男鬼变身卡| 亚洲精品成人av观看孕妇| 亚洲精品,欧美精品| 一级av片app| 免费播放大片免费观看视频在线观看| 免费少妇av软件| 九色成人免费人妻av| 人人妻人人爽人人添夜夜欢视频 | 黄片无遮挡物在线观看| 亚洲国产精品国产精品| 日本猛色少妇xxxxx猛交久久| 亚洲精华国产精华液的使用体验| 69av精品久久久久久| 高清午夜精品一区二区三区| 久久精品夜色国产| 精品久久久久久久久亚洲| 极品少妇高潮喷水抽搐| av一本久久久久| 精品酒店卫生间| 好男人视频免费观看在线| 我的老师免费观看完整版| 韩国av在线不卡| 97人妻精品一区二区三区麻豆| 久久精品国产a三级三级三级| 亚洲精品色激情综合| 小蜜桃在线观看免费完整版高清| 久久久久久久精品精品| 人体艺术视频欧美日本| 国产精品99久久99久久久不卡 | 亚洲精品一区蜜桃| 国产成人精品婷婷| 欧美老熟妇乱子伦牲交| av国产久精品久网站免费入址| 边亲边吃奶的免费视频| 啦啦啦啦在线视频资源| 99久久精品热视频| 内地一区二区视频在线| 最近中文字幕2019免费版| 国产片特级美女逼逼视频| 一区二区三区精品91| 国产在线一区二区三区精| 有码 亚洲区| 亚洲最大成人中文| 99久久精品国产国产毛片| 婷婷色麻豆天堂久久| 午夜精品国产一区二区电影 | 直男gayav资源| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩一区二区视频免费看| 少妇裸体淫交视频免费看高清| 免费av不卡在线播放| 国产成人aa在线观看| 亚洲精品日韩av片在线观看| 国产日韩欧美亚洲二区| 国产精品精品国产色婷婷| 岛国毛片在线播放| 男人爽女人下面视频在线观看| 国产精品国产三级专区第一集| 熟妇人妻不卡中文字幕| 国产一区亚洲一区在线观看| 欧美+日韩+精品| av国产久精品久网站免费入址| 男女边摸边吃奶| 99久久人妻综合| 麻豆国产97在线/欧美| 免费看光身美女| 免费观看性生交大片5| 伊人久久精品亚洲午夜| 国产精品久久久久久av不卡| 精品酒店卫生间| 亚洲国产精品成人综合色| 成年版毛片免费区| 久久人人爽av亚洲精品天堂 | 91狼人影院| 国产高清国产精品国产三级 | 免费观看av网站的网址| 日韩亚洲欧美综合| 国内精品美女久久久久久| 国产免费一区二区三区四区乱码| 亚洲经典国产精华液单| 久久6这里有精品| 成人亚洲精品一区在线观看 | 亚洲综合色惰| 国产亚洲av片在线观看秒播厂| 色吧在线观看| 涩涩av久久男人的天堂| 国产成人91sexporn| 免费黄频网站在线观看国产| 日韩三级伦理在线观看| 久热这里只有精品99| 亚洲aⅴ乱码一区二区在线播放| 国产女主播在线喷水免费视频网站| 欧美极品一区二区三区四区| 亚洲内射少妇av| 黄色一级大片看看| 水蜜桃什么品种好| 国产精品秋霞免费鲁丝片| 看黄色毛片网站| 亚洲国产最新在线播放| 神马国产精品三级电影在线观看| 免费黄网站久久成人精品| 国产精品三级大全| 少妇高潮的动态图| 狂野欧美白嫩少妇大欣赏| 精品久久久久久电影网| 国产欧美日韩精品一区二区| 亚洲高清免费不卡视频| 2021天堂中文幕一二区在线观| av线在线观看网站| av在线天堂中文字幕| 亚洲成色77777| 国产精品国产三级国产专区5o| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品乱久久久久久| 亚洲人成网站在线观看播放| 六月丁香七月| 午夜日本视频在线| tube8黄色片| 久久精品久久精品一区二区三区| 两个人的视频大全免费| 欧美高清性xxxxhd video| 噜噜噜噜噜久久久久久91| 老司机影院毛片| 啦啦啦啦在线视频资源| 高清av免费在线| 成人亚洲精品一区在线观看 | 国产精品蜜桃在线观看| 日韩av免费高清视频| 五月开心婷婷网| 一级黄片播放器| 亚洲自偷自拍三级| 黄色配什么色好看| 51国产日韩欧美| www.色视频.com| 特级一级黄色大片| 午夜爱爱视频在线播放| 一级毛片 在线播放| 久久亚洲国产成人精品v| 亚洲av免费高清在线观看| 99热国产这里只有精品6| 亚洲欧美成人精品一区二区| 国产亚洲av嫩草精品影院| 久久国产乱子免费精品| 欧美潮喷喷水| 亚洲国产最新在线播放|