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

    頁(yè)巖氣高效開(kāi)采的力學(xué)問(wèn)題與挑戰(zhàn)1)

    2017-03-10 02:58:51柳占立孟慶國(guó)詹世革黃克智
    力學(xué)學(xué)報(bào) 2017年3期

    柳占立 莊 茁 孟慶國(guó) 詹世革 黃克智

    ?(清華大學(xué)航天航空學(xué)院應(yīng)用力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,北京100084)?(國(guó)家自然科學(xué)基金委員會(huì)數(shù)理學(xué)部,北京100085)

    頁(yè)巖氣高效開(kāi)采的力學(xué)問(wèn)題與挑戰(zhàn)1)

    柳占立?,2)莊 茁?孟慶國(guó)?詹世革?黃克智?

    ?(清華大學(xué)航天航空學(xué)院應(yīng)用力學(xué)教育部重點(diǎn)實(shí)驗(yàn)室,北京100084)?(國(guó)家自然科學(xué)基金委員會(huì)數(shù)理學(xué)部,北京100085)

    頁(yè)巖氣是指賦存于富含有機(jī)質(zhì)泥頁(yè)巖中以吸附和游離狀態(tài)為主要存在方式的天然氣,中國(guó)資源量豐富,地域分布廣泛.頁(yè)巖氣開(kāi)采能緩解我國(guó)常規(guī)油氣產(chǎn)量不足、煤化石燃料引起環(huán)境污染等問(wèn)題,已成為中國(guó)綠色能源開(kāi)發(fā)的重要領(lǐng)域.盡管北美頁(yè)巖氣“革命”取得了成功,目前也僅有預(yù)期產(chǎn)量5%~15%的采收率.與北美地區(qū)相比,中國(guó)頁(yè)巖氣埋藏深,賦存條件差,自然豐度低,因此,高效開(kāi)采面臨更多的困難和挑戰(zhàn).近年來(lái),圍繞國(guó)家重大能源戰(zhàn)略需求,瞄準(zhǔn)技術(shù)發(fā)展前沿,學(xué)術(shù)界和工業(yè)界聯(lián)合對(duì)頁(yè)巖氣高效開(kāi)采的關(guān)鍵科學(xué)和技術(shù)問(wèn)題展開(kāi)研究.本文結(jié)合近三年四川、重慶地區(qū)的頁(yè)巖氣試驗(yàn)區(qū)塊遇到的新問(wèn)題,針對(duì)中國(guó)未來(lái)3500m以下深部開(kāi)采的新挑戰(zhàn),如地質(zhì)沉積、裂縫發(fā)育構(gòu)造不同、上覆壓力增加、水平應(yīng)力場(chǎng)變化等新問(wèn)題,介紹和總結(jié)了目前中國(guó)頁(yè)巖氣高效開(kāi)采面臨的力學(xué)科學(xué)問(wèn)題,主要包括多重耦合下的安全優(yōu)質(zhì)鉆完井力學(xué)理論和方法、水力壓裂體積改造和多尺度縫網(wǎng)形成機(jī)制、多尺度滲流力學(xué)特性與解吸附機(jī)理等.“深部頁(yè)巖氣高效開(kāi)采”的研究面向國(guó)家重大能源需求,科學(xué)意義重大,工程背景明確,需要工程力學(xué)、石油工程、地球物理、化學(xué)工程和環(huán)境工程等多學(xué)科專家合作,開(kāi)展理論研究、物理模擬、數(shù)值模擬及現(xiàn)場(chǎng)試驗(yàn)等綜合應(yīng)用基礎(chǔ)研究,取得高效開(kāi)采頁(yè)巖油氣理論與技術(shù)的突破.學(xué)科交叉是研究頁(yè)巖氣高效開(kāi)采問(wèn)題、突破技術(shù)瓶頸的橋梁,只有力學(xué)與石油工程、地球科學(xué)等學(xué)科實(shí)現(xiàn)深度交叉融合,才能更加有效地推動(dòng)頁(yè)巖油氣等非常規(guī)油氣資源的開(kāi)發(fā).

    頁(yè)巖氣,工程地質(zhì)力學(xué),鉆井完井理論,水力壓裂,多尺度滲流

    引言

    2015年美國(guó)頁(yè)巖氣年產(chǎn)量為4307.82億m3,約占美國(guó)天然氣總產(chǎn)量的56.22%,超過(guò)我國(guó)2015年天然氣總年產(chǎn)量(1270.97億m3).據(jù)美國(guó)能源部能源信息管理局2013年預(yù)測(cè),中國(guó)的頁(yè)巖氣儲(chǔ)量排名世界第一,占比全球儲(chǔ)量的36%,是美國(guó)的1.5倍,達(dá)到1115萬(wàn)億m3[1].加快頁(yè)巖氣資源勘探與開(kāi)發(fā)已成為頁(yè)巖氣資源大國(guó)的共同目標(biāo).北美頁(yè)巖氣的商業(yè)化開(kāi)發(fā)給世界各國(guó)的能源結(jié)構(gòu)調(diào)整帶來(lái)巨大影響,特別是在當(dāng)前我國(guó)不合理的能源消費(fèi)結(jié)構(gòu)背景下,頁(yè)巖油氣開(kāi)采將成為我國(guó)綠色能源開(kāi)發(fā)的新領(lǐng)域,從而成為緩解原油產(chǎn)量不足,降低減少煤化石燃料環(huán)境污染的有效途徑.“十三五”期間,我國(guó)將推動(dòng)能源結(jié)構(gòu)不斷優(yōu)化調(diào)整,天然氣等清潔能源需求持續(xù)加大,為頁(yè)巖油氣大規(guī)模開(kāi)發(fā)提供了寶貴機(jī)遇[23].

    盡管北美頁(yè)巖氣 “革命”取得了成功,但是美國(guó)頁(yè)巖油氣田的產(chǎn)量數(shù)據(jù)表明,目前只有預(yù)期產(chǎn)量5%~15%的油氣開(kāi)采出來(lái)[4],其原因是多方面的,如水平井的井段縮短和井筒完整性沒(méi)有達(dá)到預(yù)期,頁(yè)巖層改造的裂縫網(wǎng)絡(luò)沒(méi)有完好形成,高比例的裂縫簇沒(méi)有產(chǎn)出(據(jù)6個(gè)區(qū)塊統(tǒng)計(jì),30%的裂縫簇沒(méi)有產(chǎn)出),油氣滲流和解吸附尚未形成流動(dòng)等.相比北美以海相沉積為主的頁(yè)巖氣層,其產(chǎn)區(qū)主要分布在比較穩(wěn)定的大地構(gòu)造巖層內(nèi),而中國(guó)南方海相頁(yè)巖氣儲(chǔ)層存在于復(fù)雜山地區(qū)域和經(jīng)歷了后期地震的強(qiáng)烈改造,頁(yè)巖儲(chǔ)層以裂隙和微孔為主,層理、節(jié)理及夾層發(fā)育,儲(chǔ)層埋藏深,頁(yè)巖氣賦存環(huán)境差[5].以四川和塔里木盆地為例,頁(yè)巖氣埋藏大多在3000至6000m,且地表環(huán)境較差,多山區(qū)或沙漠,經(jīng)歷地震重構(gòu)、嚴(yán)重缺水和道路崎嶇,因此,施工成本和開(kāi)發(fā)難度很高.盡管條件困難,在“十二五”期間,頁(yè)巖氣勘探開(kāi)發(fā)取得了重大突破,2015年形成以重慶涪陵為代表的年產(chǎn)50億m3的氣田(可供3200萬(wàn)人口的年用氣量),成為北美之外第一個(gè)實(shí)現(xiàn)規(guī)?;虡I(yè)開(kāi)發(fā)的國(guó)家.但從現(xiàn)有頁(yè)巖油氣開(kāi)采技術(shù)層面看,中國(guó)的頁(yè)巖油氣高效開(kāi)采仍面臨諸多理論、方法和工程技術(shù)的挑戰(zhàn).

    因此,頁(yè)巖油氣高效開(kāi)采中的基礎(chǔ)科學(xué)理論和工程技術(shù)問(wèn)題亟需解決,主要分為4個(gè)方面:(1)頁(yè)巖油氣開(kāi)采中的鉆井完井力學(xué)理論與方法;(2)水力壓裂體積改造和多尺度縫網(wǎng)形成機(jī)制;(3)多尺度多物理場(chǎng)滲流力學(xué)特性與解吸附機(jī)理;(4)頁(yè)巖油氣開(kāi)采新概念與新技術(shù).學(xué)科交叉是研究頁(yè)巖油氣高效開(kāi)采問(wèn)題、突破技術(shù)瓶頸的橋梁,只有力學(xué)與石油工程、地球科學(xué)等學(xué)科實(shí)現(xiàn)深度交叉融合,才能更加有效地推動(dòng)頁(yè)巖油氣等非常規(guī)油氣資源的開(kāi)發(fā),實(shí)現(xiàn)“十三五”規(guī)劃中“力爭(zhēng)實(shí)現(xiàn)頁(yè)巖氣產(chǎn)量300億立方米”的目標(biāo).

    1 頁(yè)巖氣高效開(kāi)采的力學(xué)研究現(xiàn)狀和發(fā)展趨勢(shì)

    頁(yè)巖儲(chǔ)層具有低孔、低滲透率的物性特征(孔隙度一般在4%~6%,滲透率小于1.0nD),油氣溢出阻力遠(yuǎn)大于常規(guī)天然油氣,基本無(wú)自然產(chǎn)能,流體的滲透通道主要是裂縫網(wǎng)絡(luò)系統(tǒng),因此對(duì)頁(yè)巖儲(chǔ)層進(jìn)行大規(guī)模壓裂改造是形成頁(yè)巖氣工業(yè)產(chǎn)能的主要手段[67].現(xiàn)有頁(yè)巖氣增產(chǎn)方法普遍采用水壓致裂與水平井技術(shù)來(lái)增加儲(chǔ)層的宏觀透氣性[810].整體來(lái)看,我國(guó)頁(yè)巖氣開(kāi)發(fā)基礎(chǔ)理論和技術(shù)尚處于初步研究階段,例如,(1)還沒(méi)有完全掌握頁(yè)巖特征,造成前期勘探開(kāi)發(fā)遇到了諸多問(wèn)題,亟需建立全面客觀的頁(yè)巖可壓性評(píng)價(jià)指標(biāo)以指導(dǎo)我國(guó)頁(yè)巖氣的高效開(kāi)發(fā);(2)在頁(yè)巖水平井鉆井過(guò)程中,出現(xiàn)井壁垮塌嚴(yán)重、垮塌周期明顯等問(wèn)題,直接導(dǎo)致水平井段縮短、體積改造規(guī)模降低,水平井井筒完整性因頁(yè)巖氣井特有的開(kāi)發(fā)方式和井下工況受到嚴(yán)重挑戰(zhàn);(3)水力壓裂僅憑井口流量控制,裂縫簇間距沒(méi)有優(yōu)化,縫網(wǎng)形成機(jī)制尚不完全清楚,彌散裂縫密度不能定量評(píng)估;(4)頁(yè)巖基質(zhì)中氣體的微流動(dòng)(吸附、解吸、擴(kuò)散與滲流)是影響產(chǎn)氣量的決定因素,這也是頁(yè)巖氣不同于常規(guī)天然氣藏的主要區(qū)別,由于對(duì)頁(yè)巖油氣滲流機(jī)理、開(kāi)采理論認(rèn)識(shí)不清,導(dǎo)致現(xiàn)有的水平井壓裂技術(shù)基本停留在工程尺度和經(jīng)驗(yàn)?zāi)7聦用妫荒茉u(píng)估采收率,無(wú)法在更廣泛的區(qū)域形成工業(yè)產(chǎn)能.下面分4個(gè)方面敘述力學(xué)研究現(xiàn)狀和發(fā)展趨勢(shì).

    1.1 鉆井完井力學(xué)理論與方法

    此頁(yè)巖非彼頁(yè)巖,已經(jīng)成為工程界和學(xué)術(shù)界對(duì)頁(yè)巖非均質(zhì)性的共識(shí).由于我國(guó)南方海相頁(yè)巖氣儲(chǔ)層位于復(fù)雜山地,具有后期強(qiáng)改造特征,頁(yè)巖水平井鉆井完井工程面臨諸多地質(zhì)力學(xué)問(wèn)題的挑戰(zhàn),包括巖石物理和地質(zhì)力學(xué)參數(shù)的非均質(zhì)性、初始地應(yīng)力場(chǎng)特征的系統(tǒng)認(rèn)知,巖石非線性變形、強(qiáng)度各向異性與剪切摩擦等力學(xué)行為,礦物組分、有機(jī)質(zhì)的微觀力學(xué)特性以及礦物顆粒支撐和黏土膠結(jié)控制的滑移機(jī)理.因此,高效開(kāi)采頁(yè)巖氣,必須能夠評(píng)價(jià)頁(yè)巖氣儲(chǔ)層、鉆井和完井的質(zhì)量[1112].

    (1)頁(yè)巖氣儲(chǔ)層、鉆井和完井的質(zhì)量評(píng)價(jià)體系

    儲(chǔ)層質(zhì)量主要包括頁(yè)巖有機(jī)質(zhì)含量、儲(chǔ)層物性以及地層孔隙壓力,特別是區(qū)塊各異的儲(chǔ)層傷害機(jī)制等;鉆井質(zhì)量主要包括水平井靶體深度,水平段長(zhǎng)度、軌跡方位的最優(yōu)化以及全壽命周期的井筒完整性等;完井質(zhì)量主要包括頁(yè)巖脆性、微裂縫發(fā)育以及水平地應(yīng)力的非均勻性等.三者之間不是相互排斥或孤立,而是相互依存、滲透和轉(zhuǎn)化.例如,2009年—2016年3月,在四川長(zhǎng)寧—威遠(yuǎn)頁(yè)巖示范區(qū)的開(kāi)發(fā)過(guò)程中,套管變形達(dá)到30%.在壓裂的101口井中(水平井90口),32口井出現(xiàn)了不同程度的套管變形(變形點(diǎn)達(dá)47個(gè)),導(dǎo)致橋塞無(wú)法坐封到位,壓裂段數(shù)減少,單井產(chǎn)量降低,加上井筒完整性差,井生命周期降低,從而嚴(yán)重影響了頁(yè)巖氣的采收率.套管變形成為制約長(zhǎng)寧—威遠(yuǎn)區(qū)塊頁(yè)巖氣開(kāi)發(fā)的瓶頸之一.

    (2)頁(yè)巖油氣儲(chǔ)層巖石物理、地質(zhì)力學(xué)各向異性模型

    工程地質(zhì)部門基于地球物理技術(shù)和地質(zhì)力學(xué)理論,通過(guò)反演屬性橫向約束和井震結(jié)合,搭建層面、斷層和裂縫組成的構(gòu)造模型,揭示山地頁(yè)巖強(qiáng)構(gòu)造活動(dòng)特征的異常地層壓力和初始地應(yīng)力場(chǎng)分布規(guī)律,提供了“甜點(diǎn)”區(qū)位和建井鉆井指南.為了指導(dǎo)現(xiàn)場(chǎng)施工,需要優(yōu)化設(shè)計(jì)井工廠的水平井靶體深度、水平段長(zhǎng)度和軌跡方位,建立頁(yè)巖氣單井產(chǎn)量敏感性評(píng)價(jià)方法及分析預(yù)測(cè)模型.

    (3)力學(xué)與化學(xué)耦合作用下的井壁失穩(wěn)機(jī)理

    在被應(yīng)力損傷的頁(yè)巖與鉆井液發(fā)生化學(xué)耦合作用下,井壁失穩(wěn)導(dǎo)致套管變形的事件經(jīng)常發(fā)生,導(dǎo)致水平井段縮短和產(chǎn)量降低.針對(duì)長(zhǎng)裸眼水平井和工廠化鉆井完井過(guò)程中普遍存在的油基鉆井液成本高、環(huán)保壓力大,水基鉆井液井壁易垮塌,建井周期長(zhǎng)等主要技術(shù)難題,需要研發(fā)減小損害、降低濾失、平衡活度的環(huán)保型鉆井液體系.

    (4)深部地層復(fù)雜應(yīng)力環(huán)境高精度宏細(xì)觀力學(xué)測(cè)試解釋新方法

    隨著我國(guó)對(duì)深部(大于3500m)頁(yè)巖氣開(kāi)采的需求增大,通過(guò)建立高精度的宏細(xì)觀力學(xué)測(cè)試解釋新方法,研究復(fù)雜應(yīng)力條件下的頁(yè)巖宏細(xì)觀變形破壞機(jī)理和強(qiáng)度特性,揭示層理面控制的頁(yè)巖強(qiáng)度各向異性,黏土礦物控制的頁(yè)巖穩(wěn)態(tài)、非穩(wěn)態(tài)與脆性開(kāi)裂特性;基于Biot本構(gòu)關(guān)系和各向異性強(qiáng)度準(zhǔn)則對(duì)應(yīng)力模型進(jìn)行三維井壁穩(wěn)定性校核,指導(dǎo)鉆井液流量和密度設(shè)計(jì);這些研究都是深部開(kāi)采必須解決的基礎(chǔ)力學(xué)問(wèn)題.

    1.2 水力壓裂縫網(wǎng)形成機(jī)制

    水力壓裂是用帶有相位角的射孔槍沿井筒分段微爆起裂,在巖層孔壁形成初始裂縫,在水壓驅(qū)動(dòng)下沿著地層最大主應(yīng)力(壓應(yīng)力)方向形成擴(kuò)展裂縫簇,裂縫速度為幾米每秒的量級(jí)[1318].含有支撐劑的壓裂液作用在裂縫面上阻止裂縫閉合,形成擴(kuò)展的裂縫網(wǎng)絡(luò).這種縫網(wǎng)在各向異性的非均質(zhì)頁(yè)巖中是隨機(jī)分布的,而在各向同性的均質(zhì)砂巖中則形成幣狀裂縫.

    我國(guó)在頁(yè)巖水力壓裂工程方面已經(jīng)有所突破,比較成功的區(qū)塊是重慶涪陵焦石壩區(qū)塊,但是對(duì)于水力壓裂縫網(wǎng)形成機(jī)理的認(rèn)識(shí)尚不完全清楚,這也嚴(yán)重制約了頁(yè)巖油氣的高效開(kāi)發(fā).我國(guó)深層頁(yè)巖比重高,高應(yīng)力效應(yīng)引發(fā)的力學(xué)行為變化突出;多段多簇射孔、多種工作液壓裂施工等綜合因素控制著頁(yè)巖裂縫形態(tài).針對(duì)以上特征,要大幅改善體積壓裂效果,必須對(duì)水力壓裂縫網(wǎng)形成機(jī)制的關(guān)鍵問(wèn)題開(kāi)展研究.

    (1)頁(yè)巖復(fù)雜縫網(wǎng)的可壓性評(píng)價(jià)

    評(píng)價(jià)頁(yè)巖儲(chǔ)層壓裂后形成復(fù)雜裂縫網(wǎng)絡(luò)的能力是縫網(wǎng)形成機(jī)制分析和儲(chǔ)層質(zhì)量評(píng)價(jià)的重要方面.狹義上可壓性僅指脆性,廣義上則包括所有控制縫網(wǎng)復(fù)雜程度的因素[1922].礦場(chǎng)上使用最多的Jarvie礦物評(píng)價(jià)模型和Rickman力學(xué)評(píng)價(jià)模型本質(zhì)上為經(jīng)驗(yàn)公式.從地質(zhì)學(xué)角度分析,脆性礦物含量、孔隙度、沉積作用、膠結(jié)礦物等均對(duì)脆性有影響;從力學(xué)角度分析,可以從應(yīng)力與應(yīng)變關(guān)系、力學(xué)強(qiáng)度參數(shù)、巖樣破碎特征等方面提取不同的參數(shù)來(lái)表征脆性.除脆性外,天然弱面和地應(yīng)力也是決定裂縫復(fù)雜程度的重要因素.從真三軸大物模實(shí)驗(yàn)以及弱面交互模型出發(fā)提取地應(yīng)力差、弱面逼近角等參數(shù)作為評(píng)價(jià)參數(shù);實(shí)驗(yàn)發(fā)現(xiàn)天然弱面膠結(jié)面特征、注入流體粘度和速率也有重要影響.目前關(guān)于可壓性評(píng)價(jià)的研究雖然有所進(jìn)展,但仍無(wú)統(tǒng)一認(rèn)識(shí),在應(yīng)力特征、弱面力學(xué)參數(shù)、射孔與壓裂工程參數(shù)等綜合可壓性評(píng)價(jià)方面,尚無(wú)成熟的模型.

    (2)裂縫起裂機(jī)制

    隨著頁(yè)巖氣開(kāi)發(fā)深度增加,閉合應(yīng)力亦大幅增加.在不同壓裂液的擾動(dòng)下,頁(yè)巖具有不同的起裂模式.從考慮不同角度下井眼或單個(gè)射孔孔眼的圍壓應(yīng)力,計(jì)算起裂方位與壓力,逐漸發(fā)展到考慮頁(yè)巖橫觀各向同性、應(yīng)力--滲流--Biot本構(gòu)的耦合模型,以及井壁單個(gè)天然弱面、水化應(yīng)力、彈塑性變形等滲流--應(yīng)力--化學(xué)耦合起裂模型.大量的實(shí)驗(yàn)及現(xiàn)場(chǎng)觀測(cè)表明井筒工作液的擾動(dòng)存在滲吸水化行為與前置酸液的力學(xué)損傷行為,但關(guān)于力學(xué)本構(gòu)及滲流規(guī)律耦合的模型尚未建立.北美地區(qū)6個(gè)區(qū)塊的現(xiàn)場(chǎng)監(jiān)測(cè)表明裂縫簇失效率超過(guò)30%,因此,在應(yīng)力作用下的多孔眼起裂競(jìng)爭(zhēng)與裂縫簇形成機(jī)制受到關(guān)注,其研究工作剛剛開(kāi)始.

    (3)多物理場(chǎng)耦合的裂縫擴(kuò)展模式

    這是頁(yè)巖水力壓裂研究的熱點(diǎn),可以歸結(jié)為如下問(wèn)題:天然弱面的非連續(xù)、高應(yīng)力彈塑性變形、分段壓裂的多裂縫簇、圍壓作用下多場(chǎng)耦合.在0.3~1.0m尺度范圍的物理模型真三軸實(shí)驗(yàn)中,通過(guò)聲發(fā)射或計(jì)算機(jī)斷層掃描技術(shù),初步揭示了裂縫在圍壓影響下界面處的兩類行為:一類是裂縫穿越界面伴隨雙擴(kuò)、轉(zhuǎn)向行為;另一類是進(jìn)入界面受限后止裂、分叉和起裂行為.在理論模型方面,通過(guò)數(shù)值模擬引入位移不連續(xù)建立多裂縫擴(kuò)展模型,分析了裂縫間的應(yīng)力陰影效應(yīng)及競(jìng)爭(zhēng)機(jī)制;裂縫交互及多裂縫隨機(jī)轉(zhuǎn)向擴(kuò)展規(guī)律;模擬縫內(nèi)支撐劑沉降和濃度分布;基于自相似原理模擬正交復(fù)雜裂縫三維形態(tài),能夠考慮裂縫干擾、濾失和支撐劑輸送,但無(wú)法反映裂縫隨機(jī)擴(kuò)展與剪切斷裂過(guò)程.對(duì)于巖層在多相流體作用下裂縫尖端的張剪復(fù)合開(kāi)裂機(jī)制,需要綜合考慮弱面結(jié)構(gòu)及各向異性斷裂準(zhǔn)則、高圍壓下流固耦合效應(yīng)、縫內(nèi)支撐劑力學(xué)暫堵與流體濾失耦合效應(yīng)等.

    (4)支撐劑輸運(yùn)機(jī)制

    在頁(yè)巖壓裂中,支撐劑輸運(yùn)是在復(fù)雜裂縫通道中的固液兩相流動(dòng),也是支撐劑離散流、裂縫連續(xù)流與基質(zhì)濾失流的多尺度流動(dòng)過(guò)程.現(xiàn)場(chǎng)監(jiān)測(cè)反饋支撐劑的輸運(yùn)距離遠(yuǎn)遠(yuǎn)低于壓裂液的響應(yīng)距離.目前對(duì)單一裂縫中的支撐劑運(yùn)移、沉降實(shí)驗(yàn)和理論模擬較多,認(rèn)識(shí)到在清水壓裂中平板裂縫內(nèi)的支撐劑輸運(yùn)是依靠翻越砂堤維持平衡高度推進(jìn),而對(duì)復(fù)雜裂縫中的支撐劑輸運(yùn)研究較少,僅有少量實(shí)驗(yàn)研究了垂直正交裂縫中支撐劑的輸運(yùn),提出了支撐劑從主縫進(jìn)入次縫的兩種控制機(jī)制:砂堤表面顆粒向次縫滾落的重力機(jī)制和主縫流速大于臨界流速時(shí)才能攜帶支撐劑進(jìn)入次縫的動(dòng)能機(jī)制.在數(shù)值模擬方面,采用計(jì)算流體力學(xué)方法模擬了正交平板裂縫中清水輸送支撐劑,得到了與物模實(shí)驗(yàn)類似的結(jié)論.目前實(shí)驗(yàn)和模擬均未解決粗糙裂縫壁面的動(dòng)邊界、壁面高溫下壓裂液的流變性、壁面多孔多裂縫下濾失行為等問(wèn)題,反映出對(duì)支撐劑的輸運(yùn)機(jī)制尚未完全明晰.

    (5)縫網(wǎng)形成機(jī)制和支撐劑作用原理

    頁(yè)巖壓裂由近及遠(yuǎn)形成主裂縫簇+分支裂縫+次生裂縫組成的縫網(wǎng).研究發(fā)現(xiàn)主裂縫及分支裂縫以張性I型裂縫為主,次生裂縫以剪切滑移II型裂縫為主.支撐劑以多層及單層鋪置形式進(jìn)入張性裂縫,卻難以進(jìn)入剪切錯(cuò)位裂縫,后者是依靠巖石粗糙面形成的自支撐裂縫,這可能是現(xiàn)場(chǎng)監(jiān)測(cè)到支撐劑輸運(yùn)長(zhǎng)度遠(yuǎn)低于壓裂液響應(yīng)距離的原因之一.針對(duì)支撐劑多層鋪置短期導(dǎo)流能力的問(wèn)題,從支撐劑尺度、強(qiáng)度、組合、層數(shù)等方面已開(kāi)展較多研究.對(duì)于單層鋪置,閉合應(yīng)力的影響比較明顯.對(duì)于自支撐裂縫,基于實(shí)際巖板的導(dǎo)流實(shí)驗(yàn)研究表明,閉合壓力、壁面粗糙度、楊氏模量、剪切滑移量、非達(dá)西效應(yīng)對(duì)導(dǎo)流能力的影響突出,不同的頁(yè)巖導(dǎo)流能力差異最高可達(dá)2個(gè)數(shù)量級(jí).對(duì)于滑移量計(jì)算,主要基于直剪實(shí)驗(yàn)擬合或者根據(jù)無(wú)限大地層單裂縫應(yīng)變計(jì)算獲得.目前需要解決的問(wèn)題是剪切錯(cuò)位量計(jì)算、粗糙裂縫壁面的數(shù)學(xué)描述及重構(gòu)、高閉合應(yīng)力下裂縫的接觸力學(xué)、粗糙壁面條件下裂縫導(dǎo)流能力預(yù)測(cè)等.

    (6)裂縫監(jiān)測(cè)評(píng)價(jià)方法

    目前頁(yè)巖裂縫監(jiān)測(cè)評(píng)價(jià)主要針對(duì)改造體積(SRV)開(kāi)展,直接法包括微地震和測(cè)斜儀監(jiān)測(cè)等,間接法包括裂縫破裂體模擬、裂縫擴(kuò)展擬合、產(chǎn)量擬合等.井下微地震監(jiān)測(cè)是應(yīng)用最廣泛的監(jiān)測(cè)方法,通過(guò)微地震事件監(jiān)測(cè)可以確定人工裂縫的整體趨勢(shì),但仍需要對(duì)震源機(jī)制深入研究以確定裂縫的產(chǎn)生性質(zhì)、流體導(dǎo)流能力和維持時(shí)間等參數(shù).目前震源機(jī)制的反演方法主要基于縱、橫波振幅值進(jìn)行,初步認(rèn)識(shí)到不同井況水力壓裂過(guò)程中誘導(dǎo)的微地震震源機(jī)制是有差異的.進(jìn)一步需要對(duì)震源機(jī)制三維顯示方法、優(yōu)化觀測(cè)方式、結(jié)合巖石物理和施工參數(shù)進(jìn)行綜合研究.基于改造體積評(píng)價(jià)技術(shù),需要開(kāi)展縫網(wǎng)跨尺度表征、裂縫延伸特征、儲(chǔ)層參數(shù)反演等相關(guān)研究.

    1.3 多尺度滲流力學(xué)特性與解吸附機(jī)理

    頁(yè)巖油氣開(kāi)采中的滲流與解吸附是一個(gè)多尺度、多場(chǎng)、多相耦合的力學(xué)問(wèn)題[2325].立足我國(guó)頁(yè)巖油氣開(kāi)發(fā)實(shí)際,研究頁(yè)巖油氣藏滲流力學(xué)特性與解吸附機(jī)理,形成頁(yè)巖油氣多尺度滲流機(jī)理和開(kāi)采理論,對(duì)于頁(yè)巖油氣藏的高效開(kāi)發(fā)具有重要的科學(xué)價(jià)值和實(shí)際工程意義.

    (1)頁(yè)巖油氣藏孔隙結(jié)構(gòu)特征和儲(chǔ)集方式.

    頁(yè)巖孔隙類型多、結(jié)構(gòu)復(fù)雜,具有多尺度性.主要由四類不同的孔隙介質(zhì)構(gòu)成:有機(jī)質(zhì)中分布的納米級(jí)粒內(nèi)孔隙、無(wú)機(jī)礦物中納米至微米級(jí)粒間孔隙、發(fā)育豐富的微米至毫米級(jí)天然裂縫以及更大尺度的水壓裂縫.基巖孔隙是頁(yè)巖油氣藏的主要儲(chǔ)集空間,基巖非常致密,主要為納米孔隙,滲透率極低.納米級(jí)孔隙直徑在5~800nm之間,大部分在100nm左右,孔喉直徑一般在10~20nm.目前一般采用掃描電鏡、聚焦離子束掃描電鏡、原子力顯微鏡、納米計(jì)算機(jī)斷層掃描等手段觀測(cè)頁(yè)巖的孔隙結(jié)構(gòu)分布.吸附氣和游離氣共存,儲(chǔ)集方式多樣.游離氣儲(chǔ)集在孔隙空間和微裂縫中,吸附氣儲(chǔ)集在有機(jī)質(zhì)中,吸附氣可占總儲(chǔ)集量的20%~85%.對(duì)于頁(yè)巖氣儲(chǔ)層的吸附/解吸機(jī)理研究,目前沿用煤層氣固、氣界面吸附理論,即認(rèn)為滿足固、氣界面吸附的Langmuir等溫吸附/解吸關(guān)系式.實(shí)際上頁(yè)巖通常以泥巖形式在潮汐沼澤或者深水盆地的水環(huán)境下沉積,頁(yè)巖有機(jī)質(zhì)(干酪根)的降解及甲烷氣的產(chǎn)生是甲烷菌在水環(huán)境下無(wú)氧作用的結(jié)果.頁(yè)巖氣藏均存在一定的含水飽和度,成藏過(guò)程在水環(huán)境下產(chǎn)氣與聚集.因此,頁(yè)巖氣成藏與生產(chǎn)存在氣、液、固三相界面的平衡吸附與解吸.而固、液界面吸附與解吸主要與孔隙水中的氣體溶解度密切相關(guān),對(duì)壓力不敏感.目前國(guó)內(nèi)外對(duì)于頁(yè)巖氣儲(chǔ)層氣、液、固三相系統(tǒng)下的吸附及解吸機(jī)理研究匱乏,在氣藏開(kāi)發(fā)中存在吸附量與產(chǎn)量之間的矛盾,因此,亟需開(kāi)展頁(yè)巖氣實(shí)際成藏過(guò)程、生產(chǎn)過(guò)程和多尺度運(yùn)移規(guī)律的研究.

    (2)頁(yè)巖油氣在微納米孔隙內(nèi)的流動(dòng)規(guī)律

    微納米孔隙為頁(yè)巖儲(chǔ)層的主要儲(chǔ)集空間,流體在微納米孔隙內(nèi)流動(dòng),其毛細(xì)管壓力及相滲曲線與常規(guī)油氣藏存在差異,亟需開(kāi)展頁(yè)巖微納尺度多相流運(yùn)移機(jī)制的研究.國(guó)內(nèi)外在微納米級(jí)致密孔隙介質(zhì)物理模擬實(shí)驗(yàn)方面存在著微孔隙油水分布規(guī)律監(jiān)測(cè)難度大、出口微量油水計(jì)量不準(zhǔn)確等技術(shù)瓶頸,因此難以通過(guò)實(shí)驗(yàn)手段揭示致密微孔隙內(nèi)流體滲流規(guī)律.目前國(guó)內(nèi)對(duì)于微納米多孔介質(zhì)內(nèi)流動(dòng)機(jī)理的研究主要是采用數(shù)字巖心的微觀流動(dòng)模擬方法,包括圖像處理、基礎(chǔ)資料獲取及數(shù)字巖心構(gòu)建等.

    (3)頁(yè)巖油氣儲(chǔ)層介質(zhì)的多尺度流動(dòng)特性

    頁(yè)巖油氣藏儲(chǔ)層介質(zhì)具有明顯的多尺度性,基于連續(xù)介質(zhì)假設(shè)和達(dá)西方程的傳統(tǒng)滲流力學(xué)無(wú)法準(zhǔn)確刻畫(huà)頁(yè)巖氣藏的流動(dòng)機(jī)制,需要從分子尺度和介觀尺度、達(dá)西尺度和油藏尺度及尺度升級(jí)方法進(jìn)行研究和論證.在分子尺度上,形成了基于分子動(dòng)力學(xué)模擬的頁(yè)巖油氣吸附解吸及運(yùn)移機(jī)理的研究方法,研究甲烷、二氧化碳及水在黏土、有機(jī)質(zhì)及干酪根中的吸附規(guī)律,闡明頁(yè)巖油氣吸附解吸行為.在介觀尺度上,建立了頁(yè)巖多尺度數(shù)字巖心的構(gòu)建方法,以及基于數(shù)字巖心的單相和氣水兩相流動(dòng)模擬方法.根據(jù)分子尺度上流體流動(dòng)機(jī)制,分別建立了考慮微尺度效應(yīng)、氣體高壓影響和有機(jī)質(zhì)吸附影響的格子Boltzmann及孔隙網(wǎng)絡(luò)模型流動(dòng)模擬方法,能夠進(jìn)行大密度比氣水兩相流模擬,揭示了頁(yè)巖氣單相和氣水兩相孔隙尺度滲流規(guī)律.在達(dá)西尺度上,形成了基于均勻化理論的頁(yè)巖基巖尺度升級(jí)方法,考慮基巖不同孔隙介質(zhì)的賦存方式和相應(yīng)運(yùn)移機(jī)理,建立了考慮有機(jī)質(zhì)分布特征和相應(yīng)運(yùn)移機(jī)制的尺度升級(jí)數(shù)學(xué)模型,研究了微觀孔隙結(jié)構(gòu)特征對(duì)頁(yè)巖油氣宏觀流動(dòng)的影響.在油藏尺度上,建立了應(yīng)力場(chǎng)和滲流場(chǎng)耦合的頁(yè)巖裂縫擴(kuò)展及油氣開(kāi)采數(shù)值模擬方法.

    (4)頁(yè)巖油氣藏應(yīng)力、滲流、溫度三場(chǎng)耦合數(shù)值模型

    目前頁(yè)巖油氣藏的數(shù)值模擬大多未考慮應(yīng)力場(chǎng)與滲流場(chǎng)的耦合或者在滲流場(chǎng)中未考慮儲(chǔ)層耦合運(yùn)移機(jī)制.頁(yè)巖油氣藏通常發(fā)育天然微裂縫,人工壓裂后會(huì)形成多尺度離散裂縫網(wǎng)絡(luò),隨著儲(chǔ)層地應(yīng)力場(chǎng)的改變,裂縫會(huì)隨之開(kāi)啟或閉合,對(duì)滲流場(chǎng)產(chǎn)生重要影響,而滲流場(chǎng)的改變反之影響地應(yīng)力場(chǎng).目前相關(guān)的商業(yè)軟件對(duì)于此類強(qiáng)耦合問(wèn)題尚未形成模擬功能.因此,亟需建立一套適用于頁(yè)巖氣藏應(yīng)力、滲流、溫度三場(chǎng)耦合的全隱式裂縫性介質(zhì)數(shù)值模擬方法和技術(shù).

    1.4 頁(yè)巖氣開(kāi)采新概念與新技術(shù)

    深部頁(yè)巖油氣開(kāi)采是未來(lái)我國(guó)頁(yè)巖油氣開(kāi)采的重點(diǎn),然而現(xiàn)有鉆完井理論、壓裂技術(shù)和滲流解吸附方式難以滿足其高效開(kāi)發(fā)需求[7,2627].截止2016年9月,全國(guó)完鉆頁(yè)巖氣井960口,累產(chǎn)氣約114億m3,初步形成3500m淺開(kāi)發(fā)工程技術(shù)鏈.然而,四川盆地3500~4000m頁(yè)巖氣資源量為3.5萬(wàn)億m3,是 3500m淺資源量的兩倍多,可采儲(chǔ)量 0.9萬(wàn)億m3,其開(kāi)采已經(jīng)納入“十三五”規(guī)劃.不僅是我國(guó),美國(guó)2016年也成功試采開(kāi)發(fā)了超4000m深層超高壓頁(yè)巖氣藏.

    北美頁(yè)巖氣低成本開(kāi)發(fā)得益于關(guān)鍵技術(shù)的進(jìn)步.在鉆井方面,北美正研發(fā)大直徑連續(xù)管欠平衡鉆井技術(shù),并取得成功,水平井鉆速達(dá)30m/h.在完井壓裂方面,美國(guó)2016年嘗試超長(zhǎng)水平井低成本開(kāi)發(fā)頁(yè)巖氣,水平段長(zhǎng)度超過(guò)5600m,作業(yè)成本降低40%.相比北美頁(yè)巖氣開(kāi)采條件,我國(guó)開(kāi)采頁(yè)巖氣從地質(zhì)地貌條件、鉆井和壓裂工程、鉆完井成本等方面都更加復(fù)雜;部分地區(qū)水資源匱乏,生態(tài)環(huán)境脆弱;采用常規(guī)頁(yè)巖氣開(kāi)采方法和概念在深部頁(yè)巖和一些復(fù)雜地質(zhì)條件下的儲(chǔ)層中開(kāi)采成本高、效果差,甚至在壓裂后完全沒(méi)有效果.

    我國(guó)頁(yè)巖氣開(kāi)采的形勢(shì)嚴(yán)峻,面臨著3個(gè)方面的挑戰(zhàn):(1)鉆完井工程技術(shù)難題:建井成本居高,占總成本50%以上;水平井鉆速低,縫網(wǎng)壓裂改造有效性難以評(píng)定.(2)資源條件復(fù)雜:儲(chǔ)層埋深大,地質(zhì)構(gòu)造復(fù)雜、地表多丘陵,井場(chǎng)布置困難;深層地應(yīng)力高,井眼穩(wěn)定和巖石破碎難.(3)生態(tài)環(huán)境脆弱:所有入井液均帶來(lái)水資源和土壤污染風(fēng)險(xiǎn).我國(guó)水資源匱乏,但同時(shí)壓裂耗水量巨大.我國(guó)人均水資源量?jī)H2100m3,單井壓裂用水量達(dá)1.8萬(wàn)m3.針對(duì)這3個(gè)方面的挑戰(zhàn),亟需提高改造效率,增大泄氣溝通面積;創(chuàng)新技術(shù)方法,降低單井作業(yè)成本;亟需發(fā)展環(huán)境友好型技術(shù),減少耗水和廢物排放.最終實(shí)現(xiàn)優(yōu)快、低成本、環(huán)保的技術(shù)突破.發(fā)展頁(yè)巖氣高效開(kāi)采的新理念、新技術(shù),為實(shí)際施工提供新的思路.

    (1)連續(xù)管鉆完井一體化、新型壓裂液及無(wú)水壓裂技術(shù)

    連續(xù)管鉆完井一體化技術(shù),采用連續(xù)管微小井眼鉆井方法,在頁(yè)巖氣儲(chǔ)層中鉆成不同空間方位、數(shù)量、長(zhǎng)度的分支井眼,成井后直接采用連續(xù)管一體化實(shí)施水力壓裂,利用井間應(yīng)力干擾作用控制形成三維復(fù)雜裂縫網(wǎng)絡(luò).需要重點(diǎn)研究在頁(yè)巖層理弱面與地應(yīng)力競(jìng)爭(zhēng)下,徑向井壓裂裂縫起裂的模式和機(jī)理,從而控制裂縫起裂的位置及擴(kuò)展模式.壓裂液潤(rùn)滑裂縫壁面,導(dǎo)致表面能降低,從而影響裂縫擴(kuò)展.探索在徑向井多裂縫競(jìng)爭(zhēng)擴(kuò)展中的流體壓力分配和人工、天然裂縫作用等機(jī)制.使用連續(xù)管鉆完井壓裂聯(lián)作,效率高和成本低.壓裂流體選用超臨界二氧化碳或液氮等無(wú)水壓裂流體替代水壓裂,具有黏度低、表面張力低、擴(kuò)散能量強(qiáng)、壓裂更易形成復(fù)雜縫網(wǎng)等特性,從而達(dá)到安全高效、無(wú)水、環(huán)保開(kāi)采頁(yè)巖氣的目標(biāo).

    (2)頁(yè)巖氣多級(jí)壓裂井產(chǎn)能評(píng)價(jià)

    只有做到掌握真實(shí)氣體、限域效應(yīng)等對(duì)于微納尺度頁(yè)巖氣的傳輸特征與機(jī)理的影響,并傳遞到宏觀傳輸模型;定量表征基質(zhì)孔隙、天然孔隙以及巖心尺度微納米孔隙流動(dòng)與天然裂縫滲流的相互耦合,闡明其非均質(zhì)滲流特性,建立巖心尺度的滲透率模型,才能建立顯式表征非平面復(fù)雜裂縫網(wǎng)中多級(jí)壓裂井產(chǎn)能的預(yù)測(cè)模型.需要考慮自發(fā)滲析和基質(zhì)、裂縫中的多相流動(dòng),建立壓裂液返排預(yù)測(cè)模型,探索關(guān)井時(shí)期壓裂液在基質(zhì)中的自吸,以及復(fù)雜裂縫網(wǎng)、天然裂縫發(fā)育程度對(duì)于壓裂液返排的影響,從而提高壓裂液的返排率.在此方面,光纖生產(chǎn)剖面測(cè)試技術(shù)和壓裂示蹤劑動(dòng)態(tài)檢測(cè)技術(shù)是產(chǎn)能評(píng)價(jià)的重要手段,可準(zhǔn)確評(píng)價(jià)壓裂井段產(chǎn)量貢獻(xiàn)率,為優(yōu)化產(chǎn)能預(yù)測(cè)模型和完井設(shè)計(jì)提供基礎(chǔ)數(shù)據(jù).

    (3)基質(zhì)與裂隙相互作用的原位科學(xué)管理

    我國(guó)頁(yè)巖基質(zhì)的滲透率極低,在納達(dá)西水平;基質(zhì)與裂隙之間傳導(dǎo)性能的差異高達(dá)幾個(gè)數(shù)量級(jí),導(dǎo)致頁(yè)巖氣滲透率的演化成為一個(gè)漫長(zhǎng)的過(guò)程,而且基質(zhì)與裂隙之間一直處于非平衡態(tài),準(zhǔn)確定義和描述這一動(dòng)態(tài)變化是實(shí)現(xiàn)頁(yè)巖氣有效開(kāi)采的前提條件.為了應(yīng)對(duì)這些挑戰(zhàn),需要實(shí)現(xiàn)水力壓裂創(chuàng)造的頁(yè)巖傳氣能力與頁(yè)巖基質(zhì)氣體供應(yīng)能力相匹配.因此,理解頁(yè)巖基質(zhì)傳導(dǎo)性能和基質(zhì)、裂隙之間相互作用,在原位應(yīng)力、壓力以及溫度條件下的演化過(guò)程,并在此基礎(chǔ)上實(shí)現(xiàn)基質(zhì)與裂隙相互作用的原位科學(xué)管理,實(shí)現(xiàn)頁(yè)巖氣可持續(xù)開(kāi)采,成為目前頁(yè)巖油氣高效開(kāi)采的發(fā)展趨勢(shì).

    2 頁(yè)巖氣高效開(kāi)采的力學(xué)研究面臨的挑戰(zhàn)

    2.1 鉆井完井力學(xué)理論與方法

    頁(yè)巖水平井鉆井完井工程面臨諸多地質(zhì)力學(xué)問(wèn)題的挑戰(zhàn),包括巖石物理和地質(zhì)力學(xué)參數(shù)的非均質(zhì)性、初始地應(yīng)力場(chǎng)特征的系統(tǒng)認(rèn)知,巖石非線性變形、強(qiáng)度各向異性與剪切摩擦等力學(xué)行為,礦物組分、有機(jī)質(zhì)的微觀力學(xué)特性以及礦物顆粒支撐和黏土膠結(jié)控制的滑移機(jī)理.因此,需要重點(diǎn)突破如下方面的力學(xué)研究.

    基于地球物理、鉆錄測(cè)井和巖石力學(xué)一體化的地質(zhì)力學(xué)環(huán)境非均質(zhì)性、各向異性表征方法,實(shí)現(xiàn)多尺度裂縫結(jié)構(gòu)識(shí)別、異常地層壓力預(yù)測(cè)和地應(yīng)力場(chǎng)數(shù)值模擬.開(kāi)發(fā)深部地層復(fù)雜應(yīng)力環(huán)境高精度巖石力學(xué)測(cè)試解釋新方法,研究復(fù)雜應(yīng)力條件下頁(yè)巖的宏細(xì)觀變形破壞機(jī)理和強(qiáng)度特性,揭示層理面、礦物組分、膠結(jié)結(jié)構(gòu)等控制的頁(yè)巖強(qiáng)度與脆性特性,以及基于Biot本構(gòu)關(guān)系的三維井壁穩(wěn)定性校核方法.鉆井液與超低滲頁(yè)巖壓力傳遞效應(yīng)的高精度測(cè)試技術(shù),研究鉆井液與井壁圍巖應(yīng)力損傷與化學(xué)腐蝕耦合作用下的壓力、活度和損傷演化效應(yīng),周期性失穩(wěn)機(jī)理與規(guī)律.研究多維度、跨尺度的地質(zhì)力學(xué)屬性特征及其演化規(guī)律,壓裂過(guò)程中礦物組分、層面遮擋、地應(yīng)力約束的裂縫擴(kuò)展特征,地層壓力的各向異性變化,地層變形位移風(fēng)險(xiǎn)的量化表征等.

    2.2 水力壓裂縫網(wǎng)形成機(jī)制研究

    頁(yè)巖水力壓裂的目的是形成匹配地層滲流的有效網(wǎng)絡(luò)裂縫,已有研究尚在初步認(rèn)識(shí)階段,在高應(yīng)力和物理化學(xué)作用下起裂機(jī)理、復(fù)雜介質(zhì)下大型物模實(shí)驗(yàn)與數(shù)值模擬、復(fù)雜動(dòng)邊界支撐劑輸運(yùn)規(guī)律、頁(yè)巖不同支撐模式下的長(zhǎng)效導(dǎo)流能力、高精度壓裂隨機(jī)裂縫實(shí)時(shí)動(dòng)態(tài)監(jiān)測(cè)評(píng)價(jià)等方面仍存在諸多問(wèn)題和挑戰(zhàn).

    在裂縫起裂方面,工作液自發(fā)滲吸及誘導(dǎo)裂縫機(jī)理、頁(yè)巖水化及酸損傷本構(gòu)、彈塑性頁(yè)巖多簇多孔眼起裂機(jī)制目前仍認(rèn)識(shí)不清.在裂縫擴(kuò)展方面,多相流體作用下裂縫尖端應(yīng)力場(chǎng)、裂縫剪切滑移機(jī)制、高應(yīng)力力學(xué)非線性與流固耦合非線性理論、支撐劑--壓裂液--裂縫多相作用與動(dòng)態(tài)擴(kuò)展耦合機(jī)制、非連續(xù)的多裂縫隨機(jī)擴(kuò)展理論需要尋求突破.目前亟需開(kāi)展如下研究工作.

    (1)建立裂縫簇?cái)U(kuò)展歸一化條件.在裂縫簇形成過(guò)程中,水平井筒中的流量從微爆孔口分配給各壓裂段內(nèi)的裂縫簇,簇內(nèi)水壓與分配流量和射孔摩阻的乘積成正比,而裂縫尖端處地應(yīng)力主導(dǎo)的應(yīng)力強(qiáng)度因子與裂縫面上水壓力構(gòu)成驅(qū)動(dòng)裂縫擴(kuò)展的競(jìng)爭(zhēng)機(jī)制.工程上希望形成平行擴(kuò)展,行穩(wěn)至遠(yuǎn)的裂縫簇,而不希望出現(xiàn)裂縫之間的屏蔽效應(yīng),后者是北美6個(gè)頁(yè)巖氣區(qū)塊中30%裂縫簇沒(méi)有產(chǎn)出的原因之一.裂縫擴(kuò)展和止裂既要滿足流量的增加和保持,又要滿足應(yīng)力強(qiáng)度因子或者能量釋放率表征的斷裂準(zhǔn)則.由此流量增加率和斷裂準(zhǔn)則構(gòu)成了裂縫簇?cái)U(kuò)展的歸一化條件.在該條件中,應(yīng)力強(qiáng)度因子和斷裂韌性是地球物理勘測(cè)數(shù)據(jù)的已知條件,井口流量是工程控制條件,射孔摩阻是可以計(jì)算的參數(shù).這樣,由此建立裂縫簇?cái)U(kuò)展歸一化條件可以評(píng)估和優(yōu)化裂縫簇間距和擴(kuò)展長(zhǎng)度.

    (2)建立基于裂縫密度的采收率評(píng)估模型.在水平井的壓裂過(guò)程中,不但需要形成主裂縫簇,更重要的是形成成百上千彌散裂縫的縫網(wǎng)結(jié)構(gòu),這些裂縫包含水力裂縫本身的隨機(jī)擴(kuò)展和分叉,以及它們與天然裂縫的交匯.如何對(duì)裂縫簇和縫網(wǎng)進(jìn)行力學(xué)表征和數(shù)值模擬,對(duì)流固耦合的斷裂力學(xué)理論和模擬計(jì)算提出全新的挑戰(zhàn).根據(jù)北美5個(gè)區(qū)塊26個(gè)月油氣產(chǎn)量調(diào)查統(tǒng)計(jì),保證頁(yè)巖氣充分解吸附的合理裂縫間距是0.1m,該數(shù)據(jù)提供了一種逆向推算思路,目前國(guó)際上由此形成了兩種分析模型,基于本構(gòu)層次的損傷模型和基于多尺度自相似裂縫密度模型.國(guó)外這些統(tǒng)計(jì)數(shù)據(jù)如何對(duì)比我國(guó)頁(yè)巖氣的現(xiàn)場(chǎng)裂縫間距數(shù)據(jù),結(jié)合井口控制流量和孔口摩阻系數(shù),可以推算裂縫密度和評(píng)估采收率.研究的目標(biāo)是能夠開(kāi)發(fā)頁(yè)巖氣工程模擬軟件,指導(dǎo)現(xiàn)場(chǎng)工程.

    (3)發(fā)展大物模實(shí)驗(yàn)技術(shù)和提高數(shù)據(jù)處理水平.水壓過(guò)程對(duì)裂縫擴(kuò)展的影響主要體現(xiàn)在流體對(duì)裂縫尖端亞臨界擴(kuò)展的影響,和水力的加載方式,即地應(yīng)力圍壓下裂縫面內(nèi)分布力加載的影響.水力裂縫面臨的另一個(gè)問(wèn)題是后期的檢測(cè)和表征.在室內(nèi)的大物模實(shí)驗(yàn)中,如中石油勘探院廊坊分院的 1.0立方米三維水壓實(shí)驗(yàn)設(shè)備,是目前亞洲最大尺度的露頭巖石的物理實(shí)驗(yàn)平臺(tái),樣品尺寸為762mm×762mm×914mm,最大載荷為69MPa.以聲發(fā)射和截面觀測(cè)為主的表征方式仍然與現(xiàn)場(chǎng)幾千米深的壓裂施工差別很大,無(wú)法給出實(shí)際工程中流量與裂縫速度、擴(kuò)展長(zhǎng)度和多裂縫相互作用影響間距等真實(shí)數(shù)據(jù).因此,水力裂縫的準(zhǔn)確監(jiān)測(cè)與表征仍然是亟待解決的問(wèn)題.

    (4)揭示支撐劑運(yùn)輸和裂縫導(dǎo)流機(jī)制.在支撐劑輸運(yùn)方面,動(dòng)邊界下高溫、高壓、高應(yīng)力下固液多相流理論、多重介質(zhì)下流體競(jìng)爭(zhēng)濾失理論、多裂縫壁面效應(yīng)與支撐劑分布預(yù)測(cè)模型是需要解決的難題.在裂縫導(dǎo)流能力方面,頁(yè)巖--支撐劑接觸與頁(yè)巖粗糙錯(cuò)位縫接觸動(dòng)力學(xué)理論、高溫高應(yīng)力下低濃度支撐裂縫與自支撐裂縫的長(zhǎng)效導(dǎo)流實(shí)驗(yàn)與理論需要進(jìn)行研究.在裂縫評(píng)價(jià)方面,跨尺度三維網(wǎng)絡(luò)裂縫表征方法、高精度三維震源機(jī)制與力學(xué)反演理論、基于施工參數(shù)的裂縫延伸反演理論需要重點(diǎn)突破.

    (5)發(fā)展頁(yè)巖壓裂數(shù)值模擬軟件.目前國(guó)際上水力壓裂數(shù)值模擬形成了一些軟件體系,從整體區(qū)塊的壓裂設(shè)計(jì),到單井壓裂設(shè)計(jì)和實(shí)時(shí)監(jiān)測(cè)分析等方面,都有相應(yīng)的軟件,主要有 FracpropPT,Meyer,TerraFrac,E-StimPlan等,這些軟件的開(kāi)發(fā)和使用大大促進(jìn)了水力壓裂技術(shù)的發(fā)展.但是這些軟件大多采用非常簡(jiǎn)化的力學(xué)模型,認(rèn)為儲(chǔ)層是各向同性的均勻材料,采用解析或半解析的裂縫模型,把裂縫簡(jiǎn)化成對(duì)稱平面裂紋,忽略了裂縫復(fù)雜的幾何形態(tài)及多個(gè)裂縫之間的相互作用,難以應(yīng)對(duì)非均質(zhì)各向異性頁(yè)巖中復(fù)雜縫網(wǎng)系統(tǒng)的模擬.針對(duì)頁(yè)巖壓裂的復(fù)雜特點(diǎn),在頁(yè)巖壓裂領(lǐng)域急需發(fā)展符合真實(shí)物理狀況的水力壓裂數(shù)值模擬工具.

    2.3 多尺度滲流力學(xué)特性與解吸附機(jī)理研究

    頁(yè)巖儲(chǔ)層因其多尺度的儲(chǔ)集空間和多樣性的油氣賦存方式,在開(kāi)采過(guò)程中的運(yùn)移是應(yīng)力、多相滲流的耦合過(guò)程,常規(guī)的油氣滲流理論不再適用,亟需開(kāi)展頁(yè)巖油氣多尺度滲流特征及開(kāi)采理論的研究,面臨以下挑戰(zhàn)性問(wèn)題:

    (1)頁(yè)巖氣、液、固三相系統(tǒng)吸附/解吸機(jī)理及規(guī)律.前期關(guān)于頁(yè)巖氣藏的吸附/解吸機(jī)理的研究只是局限于氣、固兩相系統(tǒng),其吸附/解吸機(jī)理與頁(yè)巖氣藏礦場(chǎng)生產(chǎn)特征差別較大,大規(guī)模水力壓裂造成很多水滯留在多孔介質(zhì)中,使得有機(jī)質(zhì)的氣體很難解吸釋放.工程實(shí)踐和實(shí)驗(yàn)證明:水相的存在對(duì)吸附/解吸影響較大,因此有必要研究氣、液、固三相系統(tǒng)頁(yè)巖氣的吸附/解吸機(jī)制.

    (2)頁(yè)巖微納尺度多相流運(yùn)移機(jī)制及其尺度升級(jí).頁(yè)巖不同則孔隙介質(zhì)運(yùn)移機(jī)制不同,亟需構(gòu)建考慮不同孔隙介質(zhì)類型的頁(yè)巖數(shù)字巖心,全面考慮耦合運(yùn)移機(jī)制,開(kāi)展微納尺度氣水、油水兩相流動(dòng)模擬,揭示頁(yè)巖儲(chǔ)層微納尺度滲流規(guī)律,開(kāi)展尺度升級(jí)研究,得到考慮微觀運(yùn)移機(jī)制的宏觀滲流模型.

    (3)應(yīng)力場(chǎng)與滲流場(chǎng)全耦合的頁(yè)巖油氣藏?cái)?shù)值模擬.在實(shí)際工程中,應(yīng)力場(chǎng)與滲流場(chǎng)發(fā)生相互影響,因此,需要建立頁(yè)巖油氣藏力學(xué)模型,開(kāi)展頁(yè)巖油氣藏的應(yīng)力場(chǎng)與滲流場(chǎng)完全耦合的數(shù)值模擬研究,形成高效求解方法.

    (4)分段壓裂水平井非穩(wěn)態(tài)壓力產(chǎn)量分析及開(kāi)發(fā)潛力預(yù)測(cè)模型.分段壓裂水平井是頁(yè)巖油氣藏開(kāi)發(fā)的主要手段,目前非穩(wěn)態(tài)壓力(試井)解釋及產(chǎn)能預(yù)測(cè)方法未能全面考慮頁(yè)巖儲(chǔ)層內(nèi)的耦合運(yùn)移機(jī)制,因此亟需建立基于頁(yè)巖儲(chǔ)層耦合運(yùn)移機(jī)制的滲流數(shù)學(xué)模型,開(kāi)展頁(yè)巖油氣藏非穩(wěn)態(tài)壓力(試井)解釋及產(chǎn)能預(yù)測(cè)方法研究.另外,在井工廠模式下,需要考慮同步壓裂應(yīng)力干擾水平井和人工裂縫,揭示井工廠開(kāi)采模式下的水平井位置以及分段壓裂縫的空間配置,形成井工廠優(yōu)化設(shè)計(jì)方法.

    2.4 頁(yè)巖氣開(kāi)采新概念與新技術(shù)

    頁(yè)巖基質(zhì)的超低滲和基質(zhì)與裂縫傳導(dǎo)性能之間的高差異是頁(yè)巖氣藏的兩個(gè)最大特征.頁(yè)巖的滲透性可以低到納達(dá)西水平,基質(zhì)與裂縫之間傳導(dǎo)性能的差異可以高到幾個(gè)數(shù)量級(jí).這兩個(gè)特征決定了頁(yè)巖滲透率的演化將是一個(gè)漫長(zhǎng)的過(guò)程,而且基質(zhì)與裂縫之間一直處于非穩(wěn)定非平衡的狀態(tài),準(zhǔn)確定義和描述這個(gè)動(dòng)態(tài)系統(tǒng)是頁(yè)巖氣高效開(kāi)采的核心理論模型.目前的研究不但要關(guān)注水力壓裂創(chuàng)造頁(yè)巖氣通道的能力,也要關(guān)注頁(yè)巖氣能否從基質(zhì)中流出.美國(guó)頁(yè)巖氣開(kāi)采面臨可持續(xù)性的挑戰(zhàn),中國(guó)頁(yè)巖氣開(kāi)采面臨技術(shù)和方法上的突破.為了應(yīng)對(duì)這些挑戰(zhàn),提出了實(shí)現(xiàn)水力壓裂創(chuàng)造的頁(yè)巖氣傳氣能力與基質(zhì)氣體供應(yīng)能力相匹配,避免“斷氣”現(xiàn)象發(fā)生,最終達(dá)到頁(yè)巖氣可持續(xù)性開(kāi)采的新理論概念.研究重點(diǎn)是理解頁(yè)巖基質(zhì)傳導(dǎo)性能和基質(zhì)--裂縫之間相互作用在原位應(yīng)力、壓力及溫度條件下是如何演化的,并在此基礎(chǔ)上實(shí)現(xiàn)從制造裂縫到基質(zhì)與裂縫相互作用的原位科學(xué)管控.

    面對(duì)頁(yè)巖氣深部開(kāi)采的挑戰(zhàn),發(fā)展連續(xù)管鉆完井一體化技術(shù).頁(yè)巖堅(jiān)硬致密,鉆井破巖效率低,徑向井眼成孔難,如何鉆成不同空間方位、數(shù)量和長(zhǎng)度分支井眼?如何通過(guò)徑向孔眼間應(yīng)力干擾,設(shè)計(jì)、控制并形成復(fù)雜裂縫網(wǎng)絡(luò)?如何準(zhǔn)確評(píng)價(jià)頁(yè)巖氣井產(chǎn)能?回答這些問(wèn)題,需要發(fā)展連續(xù)管鉆完井一體化技術(shù):采用連續(xù)管鉆多層、多孔徑向井,加速解吸、擴(kuò)散、滲流速度,徑向井壓裂一體化形成復(fù)雜裂縫網(wǎng)絡(luò),增大泄氣面積,從而實(shí)現(xiàn)“一井多層、一層多孔、一孔多縫”的復(fù)雜三維裂縫網(wǎng)絡(luò),達(dá)到連續(xù)管鉆完井壓裂聯(lián)作,高效率、低成本開(kāi)采頁(yè)巖油氣的目的.同時(shí)建立壓裂縫網(wǎng)形成機(jī)制和基于微納尺度流動(dòng)的產(chǎn)能評(píng)價(jià)方法,

    總的來(lái)說(shuō),連續(xù)管鉆完井一體化壓裂縫網(wǎng)形成機(jī)制研究;新型壓裂液及無(wú)水壓裂機(jī)理研究;頁(yè)巖油氣井的開(kāi)采壽命優(yōu)化研究;安全環(huán)保頁(yè)巖油氣開(kāi)采新概念和新方法研究等問(wèn)題是目前頁(yè)巖油氣高效開(kāi)采的力學(xué)研究面臨的挑戰(zhàn).

    3 結(jié)論

    根據(jù)近三年四川、重慶地區(qū)的頁(yè)巖氣試驗(yàn)區(qū)塊遇到的新問(wèn)題,中國(guó)未來(lái)3500m以下深部開(kāi)采的新挑戰(zhàn),如地質(zhì)沉積、裂縫發(fā)育構(gòu)造不同,上覆壓力增加,水平應(yīng)力場(chǎng)變化等新問(wèn)題,目前我國(guó)頁(yè)巖油氣高效開(kāi)采面臨的主要力學(xué)科學(xué)問(wèn)題如下:

    (1)多重耦合下的安全優(yōu)質(zhì)鉆完井力學(xué)理論和方法

    井壁圍巖穩(wěn)定性的多尺度、多物理場(chǎng)耦合力學(xué)理論研究;深部頁(yè)巖(3500m以下)水平井穩(wěn)定性控制理論研究;鉆采過(guò)程中頁(yè)巖儲(chǔ)層傷害機(jī)理和地質(zhì)力學(xué)特征演化研究;地層--水泥環(huán)--套管系統(tǒng)在復(fù)雜工況條件的破壞規(guī)律研究;水基鉆井液適應(yīng)性研究;非均質(zhì)各向異性含孔隙頁(yè)巖本構(gòu)理論研究.

    (2)多尺度水力壓裂縫網(wǎng)形成機(jī)理

    深部頁(yè)巖壓裂致裂機(jī)理研究;多簇裂縫起裂和擴(kuò)展與流場(chǎng)分布機(jī)理研究;頁(yè)巖裂縫網(wǎng)的跨尺度表征與采收率評(píng)價(jià)體系研究;多物理場(chǎng)耦合和多尺度水力壓裂數(shù)值模擬和裂縫管控技術(shù)研究;大物模實(shí)驗(yàn)和現(xiàn)場(chǎng)數(shù)據(jù)對(duì)數(shù)模的驗(yàn)證研究;裂縫縫隙的支撐劑輸運(yùn)機(jī)制與接觸力學(xué)理論模型研究.

    (3)多尺度滲流力學(xué)特性與解吸附機(jī)理

    深部頁(yè)巖基質(zhì)滲流的微觀機(jī)制研究;數(shù)字巖心再現(xiàn)技術(shù)研究;多尺度、多相、多物理場(chǎng)耦合的滲流機(jī)理研究;尺度升級(jí)方法和滲透率計(jì)算研究;裂縫擴(kuò)展和縫網(wǎng)演化過(guò)程中油氣體解吸、驅(qū)替和滲流特性研究.井工廠開(kāi)采模式下同步壓裂裂縫擴(kuò)展和多相多場(chǎng)耦合流動(dòng)的三維數(shù)值模擬的大規(guī)模計(jì)算.

    “深部頁(yè)巖油氣高效開(kāi)采”的研究面向國(guó)家重大能源需求,科學(xué)意義重大,工程背景明確,需要工程力學(xué)、石油工程、地球物理、化學(xué)工程和環(huán)境工程等多學(xué)科交叉合作,開(kāi)展理論研究、物理模擬、數(shù)值模擬及現(xiàn)場(chǎng)試驗(yàn)等綜合應(yīng)用基礎(chǔ)研究,取得高效開(kāi)采頁(yè)巖油氣理論與技術(shù)的突破.

    致謝 感謝上海大學(xué)周哲瑋教授,中國(guó)石油大學(xué)(華東)姚軍教授、薛世峰教授,中國(guó)石油大學(xué)(北京)李根生院士、陳勉教授,北京大學(xué)張東曉教授,中石油鉆井院楊恒林高工、田中蘭高工、項(xiàng)德貴高工、陳朝偉高工,西南石油大學(xué)郭建春教授,中國(guó)科學(xué)院力學(xué)研究所趙亞溥研究員、袁泉子副研究員、劉曰武研究員、蘇業(yè)旺研究員,中國(guó)科學(xué)院巖土力學(xué)研究所劉繼山教授等對(duì)本文的貢獻(xiàn).

    1 U.S.Energy Information Administration.Shale gas production drives world natural gas production growth.2006.http://www.eia.gov/todayinenergy/detail.php?id=27512

    2王南,裴玲,雷丹鳳等.中國(guó)非常規(guī)天然氣資源分布及開(kāi)發(fā)現(xiàn)狀.油氣地質(zhì)與采收率,2015,22(1):26-31(Wang Nan,Pei Ling,Lei Danfeng,et al.Analysis of unconventional gas resources distribution and development status in China.Petroleum Geology and Recovery E ffi ciency,2015,22(1):26-31(in Chinese))

    3董大忠,鄒才能,楊樺等.中國(guó)頁(yè)巖氣勘探開(kāi)發(fā)進(jìn)展與發(fā)展前景.石油學(xué)報(bào),2012,33(s1):107-114(Dong Dazhong,Zou Caineng,Yang Hua,et al.Progress and prospects of shale gas exploration and development in China.Acta Petrolei Sinica,2012,33(s1):107-114(in Chinese))

    4 Ba?zant ZP,Salviato M,Chau VT,et al.Why fracking works.Journal of Applied Mechanics,2014,81(10):101010

    5李玉喜,張金川.我國(guó)非常規(guī)油氣資源類型和潛力.國(guó)際石油經(jīng)濟(jì),2011,3:61-67(Li Yuxi,Zhang Jinchuan.Unconventional gas&oil resource type and potential of China.International Petroleum Economy,2011,3:61-67(in Chinese))

    6 Nelson PH.Pore-throat sizes in sandstones,tight sandstones,and shales.AAPG Bulletin,2009,93(3):329-340

    7鄒才能,朱如凱,吳松濤等.常規(guī)與非常規(guī)油氣聚集類型,特征,機(jī)理及展望——以中國(guó)致密油和致密氣為例.石油學(xué)報(bào),2012,33(2):173-187(Zou Caineng,Zhu Rukai,Wu Songtao,et al.Types,Characteristics,genesis and prospects of conventional and unconventional hydrocarbon accumulations:taking tight oil and tight gas in China as an instance.Acta Petrolei Sinica,2012,33(2):173-187(in Chinese))

    8 Nobakht M,Clarkson CR,Kaviani D.New type curves for analyzing horizontal well with multiple fractures in shale gas reservoirs//Proceedingsof the SPE UnconventionalResources Conference,Calgary,Alberta,Canada,15-17 November 2011.Society of Petroleum Engineers,2011

    9 Lewis AM,Hughes RG.Production data analysis of shale gas reservoirs//Proceedings of the SPE Annual Technical Conference and Exhibition,Denver,Colorado,USA,21-24 September 2008.Society of Petroleum Engineers,2008

    10朱光普,姚軍,樊冬艷等.頁(yè)巖氣藏壓裂水平井試井分析.力學(xué)學(xué)報(bào),2015,47(6):945-954(Zhu Guangpu,Yao Jun,Fan Dongyan,et al.Pressure transient analysis of fractured horizontal well in shale gas reservoir.Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):945-954(in Chinese))

    11 Arthur JD,Bohm B,Layne M.Evaluating implications of hydraulic fracturing in shale gas reservoirs.SPE 121038,2009

    12 Akbarnejad-Nesheli B,Valko PP,Lee WJ.Relating fracture network characteristics to shale gas reserve estimation.SPE 154841,2012 Warpinski NR,Branagm PT.Altered-stress fracturing.SPE 17533,1989

    13 Soliman MY,East L,Adams D.Geomechnics aspects of multiple fracturing of horizontal and vertical wells.SPE 86992,2004

    14 Cheng Y.Impacts of the number of perforation clusters and cluster spacing on production performance of horizontal shale gas wells.SPE 138843,2010

    15 Wikel K.Geomechanics:bridging the gap from geophysics to engineering in unconventional reservoirs.First Break,2011,29:71-81

    16 Wu K,Olson JE.Investigation of critical in situ and injection factors in multi-frac treatments:guidelines for controlling fracture complexity.SPE 163821,2013

    17陳作,薛承瑾,蔣廷學(xué)等.頁(yè)巖氣井體積壓裂技術(shù)在我國(guó)的應(yīng)用建議.天然氣工業(yè),2010,30(10):30-32(ChenZuo,XueChengjin,Jiang Tingxue,et al.Proposals for the application of fracturing by stimulated reservoir volume(SRV)in shale gas wells in China.Natural Gas Industry,2010,30(10):30-32(in Chinese))

    18張廣明,劉勇,劉建東等.頁(yè)巖儲(chǔ)層體積壓裂的地應(yīng)力變化研究.力學(xué)學(xué)報(bào),2015,47(6):965-972(Zhang Guangming,Liu Yong,Liu Jiandong,et al.Research on the geostress change of shale reservoir volume fracturing.Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):965-972(in Chinese))

    19陳昀,金衍,陳勉.基于能量耗散的巖石脆性評(píng)價(jià)方法.力學(xué)學(xué)報(bào),2015,47(6):984-993(Chen Yun,Jin Yan,Chen Mian.A rock brittleness evaluation method based on energy dissipation.Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):984-993(in Chinese))

    20隋麗麗,楊永明,鞠楊等.巖石可壓裂性分形描述方法初探.力學(xué)與實(shí)踐,2014,36(6):753-756(Sui Lili,Yang Yongming,Ju Yang,et al.Fractal description of rock fracture behavior.Mechanics in Engineering,2014,36(6):753-756(in Chinese))

    21 Chong KK,Grieser WV,Passman A.A completions guide book to shale-play development:A review of successful approaches toward shale-play stimulation in the last two decades//Proc.of Canadian Unconventional Resources and International Petroleum Conference,Calgary,Alberta,Canada,SPE133874,2010

    22曾青冬,姚軍,孫致學(xué).頁(yè)巖氣藏壓裂縫網(wǎng)擴(kuò)展數(shù)值模擬.力學(xué)學(xué)報(bào),2015,47(5):994-999(Zeng Qingdong,Yao Jun,Sun Zhixue.Numerical modeling of fracture network propagation in shale reservoirs.Chinese Journal of Theoretical and Applied Mechanics,2015,47(5):994-999(in Chinese))

    23魏明強(qiáng),段永剛,方全堂等.頁(yè)巖氣藏孔滲結(jié)構(gòu)特征和滲流機(jī)理研究現(xiàn)狀.油氣藏評(píng)價(jià)與開(kāi)發(fā),2011,1(4):73-77(Wei Mingqiang,Duan Yonggang,Fang Quantang,et al.Current research situation of porosity&permeability characteristics and seepage mechanism of shale gas reservoir.Reservoir Evaluation and Development,2011,1(4):73-77(in Chinese))

    24 Loucks RG,Reed RM,Ruppel SC,et al.Morphology,genesis,and distribution of nanometer-scale pores in siliceous.Mudstones of the Mississippian Barnett Shale.Journal of Sedimentary Research,2009,79(11-12):848-861

    25樊冬艷,姚軍,孫海等.考慮多重運(yùn)移機(jī)制耦合頁(yè)巖氣藏壓裂水平井?dāng)?shù)值模擬.力學(xué)學(xué)報(bào),2015,47(6):906-915(Fan Dongyan,Yao Jun,Sun Hai,et al.Numerical simulation of multi-fractured horizontal well in shale gas reservoir considering multiple gas transport mechanisms.Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):906-915(in Chinese))

    26鄒才能,朱如凱,吳松濤等.常規(guī)與非常規(guī)油氣聚集類型,特征,機(jī)理及展望——以中國(guó)致密油和致密氣為例.石油學(xué)報(bào),2012,33(2):173-187(Zou Caineng,Zhu Rukai,Wu Songtao,et al.Types,Characteristics,genesis and prospects of conventional and unconventional hydrocarbon accumulations:taking tight oil and tight gas in China as an instance.Acta Petrolei Sinica,2012,33(2):173-187(in Chinese))

    27李道倫,鄭德溫,方朝合等.頁(yè)巖氣井組分比例變化規(guī)律研究.力學(xué)學(xué)報(bào),2015,47(6):899-905(Li Daolun,Zheng Dewen,Fang Chaohe,et al.Sensitivity study of fl wing gas composition for wells in shale gas reservoir.Chinese Journal of Theoretical and Applied Mechanics,2015,47(6):899-905(in Chinese))

    28張東曉,楊婷云.頁(yè)巖氣開(kāi)發(fā)綜述.石油學(xué)報(bào),2013,34(4):792-801(Zhang Dongxiao,Yang Tingyun.An overview of shale-gas production.Acta Petrolei Sinica,2013,34(4):792-801(in Chinese))

    PROBLEMS AND CHALLENGES OF MECHANICS IN SHALE GAS EFFICIENT EXPLOITATION1)

    Liu Zhanli?,2)Zhuang Zhuo?Meng Qingguo?Zhan Shige?Huang Keh-Chih??(Applied Mechanics Lab.,School of Aerospace Engineering,Tsinghua University,Beijing 100084,China)?(Department of Mathematical&Physical Science,National Natural Science Foundation of China,Beijing 100085,China)

    Shale gas is unconventional natural gas stored in shale in free or absorbed forms and sometimes in free flui phase.The exploitation of shale gas has become a promising fiel of green energy development in China.Although great success has been achieved in shale gas revolution in North America with the technique of hydraulic fracturing,there is only 5%~15%of the stored oil and gas could be exploited,which is still a puzzle for petroleum engineers.Compared with the North America,China’s shale gas reservoirs are deep burial,the geologic construction conditions are complicated and natural quality is low,therefore,efficient exploitation is facing more difficulties and challenges.In recent years,aimingat the national major energy strategy and the frontier of technological development,China’s academia and industry have carried out the preliminary study on some of the key scientifi and technical issues.Around the new issues encountered in the shale gas extraction in Sichuan and Chongqing areas in recent three years,this paper introduces and summarizes the key mechanics problems and challenges that the high efficient shale gas extraction is facing,mainly includes the multifiel coupling safe and high quality drilling mechanics,hydraulic fracturing and multi-scale fracture network formation mechanism and multi-scale seepage and desorption mechanism of shale gas,to solve the challenges in deep exploitation below 3500 meters in China,such as geologic sedimentation,di ff erent fracture development,increasing overburden pressure,the change of horizontal stress,etc.The deep shale gas exploitation is not only to adapt to the national energy demand,but also has scientifi and engineering significance To realize the efficient exploitation of shale oil and gas,it needs the interdisciplinary collaboration of mechanical engineering,petroleum engineering,geophysics,chemical engineering and environmental engineering to carry out basic theoretical research,physical simulation,numerical simulation and fiel experiment.It has been recognized that interdisciplinary research is the bridge and the key to breakthrough the technology bottleneck and realize the efficient exploitation of shale gas.It is necessary of the deep collaboration between mechanics,petroleum engineering,earth science and other disciplines to promote the development of shale gas and other unconventional oil and gas resources.

    shale gas,engineering geology mechanics,well drilling and completion theory,hydraulic fracture,multiscale seepage

    TE37

    :A

    10.6052/0459-1879-16-399

    2016–12–27 收稿,2017–04–26 錄用,2017–05–03 網(wǎng)絡(luò)版發(fā)表.

    1)國(guó)家自然科學(xué)基金資助項(xiàng)目(11532008,11372157).

    2)柳占立,副教授,主要研究方向:塑性力學(xué),斷裂力學(xué)和復(fù)合材料力學(xué).E-mail:liuzhanli@tsinghua.edu.cn

    柳占立,莊茁,孟慶國(guó),詹世革,黃克智.頁(yè)巖氣高效開(kāi)采的力學(xué)問(wèn)題與挑戰(zhàn).力學(xué)學(xué)報(bào),2017,49(3):507-516

    Liu Zhanli,Zhuang Zhuo,Meng Qingguo,Zhan Shige,Huang Keh-Chih.Problems and challenges of mechanics in shale gas efficient exploitation.Chinese Journal of Theoretical and Applied Mechanics,2017,49(3):507-516

    色吧在线观看| 日韩强制内射视频| 嫩草影院入口| 欧美成人午夜免费资源| 欧美激情极品国产一区二区三区 | 国产男女超爽视频在线观看| 欧美日韩视频精品一区| av天堂久久9| 水蜜桃什么品种好| 国产一区有黄有色的免费视频| 免费看不卡的av| 精品一品国产午夜福利视频| 亚洲高清免费不卡视频| 九色成人免费人妻av| 久久久精品94久久精品| 99九九线精品视频在线观看视频| 日本av手机在线免费观看| 色视频在线一区二区三区| 精品国产国语对白av| 丝袜脚勾引网站| 免费看av在线观看网站| 99久久人妻综合| 一级爰片在线观看| 欧美+日韩+精品| 麻豆成人av视频| 成人亚洲欧美一区二区av| 免费大片黄手机在线观看| 成年人免费黄色播放视频 | 亚洲av免费高清在线观看| 亚洲高清免费不卡视频| 人妻人人澡人人爽人人| 亚洲精品国产av成人精品| 男人和女人高潮做爰伦理| 99久久精品国产国产毛片| 下体分泌物呈黄色| 国产精品一二三区在线看| 蜜桃久久精品国产亚洲av| 中文在线观看免费www的网站| 欧美人与善性xxx| 只有这里有精品99| av在线播放精品| 亚洲av免费高清在线观看| 韩国av在线不卡| 亚洲精品中文字幕在线视频 | 欧美人与善性xxx| 日本黄色日本黄色录像| 少妇被粗大的猛进出69影院 | 日本猛色少妇xxxxx猛交久久| 91久久精品国产一区二区三区| 久久综合国产亚洲精品| 久久韩国三级中文字幕| 一级毛片我不卡| 在线 av 中文字幕| 欧美丝袜亚洲另类| 免费黄网站久久成人精品| 成年女人在线观看亚洲视频| 久久av网站| 老司机影院毛片| 街头女战士在线观看网站| 亚洲国产精品一区三区| 欧美日韩亚洲高清精品| 国产高清不卡午夜福利| 少妇人妻久久综合中文| 春色校园在线视频观看| 久久久精品免费免费高清| 久久久久久久久久久久大奶| 亚洲情色 制服丝袜| 免费观看的影片在线观看| 久久久久国产精品人妻一区二区| 国产精品99久久久久久久久| 一区二区av电影网| 日韩欧美精品免费久久| 久久狼人影院| 一区二区三区精品91| 偷拍熟女少妇极品色| 久久精品国产亚洲网站| 极品少妇高潮喷水抽搐| 国产欧美日韩精品一区二区| 男人和女人高潮做爰伦理| 国产精品一区www在线观看| 又黄又爽又刺激的免费视频.| 不卡视频在线观看欧美| 大话2 男鬼变身卡| 又大又黄又爽视频免费| 有码 亚洲区| av卡一久久| 国产精品一二三区在线看| 国产色婷婷99| 你懂的网址亚洲精品在线观看| 精品久久久久久电影网| 一二三四中文在线观看免费高清| a级毛色黄片| 久久久欧美国产精品| 肉色欧美久久久久久久蜜桃| 视频区图区小说| av天堂久久9| 老司机亚洲免费影院| 九九爱精品视频在线观看| 在线精品无人区一区二区三| 在线观看av片永久免费下载| 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 亚洲精品日本国产第一区| 亚洲图色成人| 亚洲国产日韩一区二区| av播播在线观看一区| 一级片'在线观看视频| 久久6这里有精品| 日韩 亚洲 欧美在线| 欧美一级a爱片免费观看看| 黄色日韩在线| 免费观看a级毛片全部| 爱豆传媒免费全集在线观看| 日本爱情动作片www.在线观看| 女性被躁到高潮视频| 国产精品蜜桃在线观看| 少妇被粗大的猛进出69影院 | av免费观看日本| 九色成人免费人妻av| 亚洲va在线va天堂va国产| 人妻少妇偷人精品九色| 国产老妇伦熟女老妇高清| 成人二区视频| 特大巨黑吊av在线直播| 97在线人人人人妻| 国产av一区二区精品久久| 亚洲,欧美,日韩| 日日啪夜夜撸| 最新中文字幕久久久久| 国产一区二区在线观看av| 欧美精品国产亚洲| 婷婷色综合www| 久久这里有精品视频免费| 自拍偷自拍亚洲精品老妇| 国产av国产精品国产| 又爽又黄a免费视频| 国产免费福利视频在线观看| 色哟哟·www| 99热这里只有是精品在线观看| www.av在线官网国产| 成人国产av品久久久| 26uuu在线亚洲综合色| 嫩草影院入口| 在线播放无遮挡| 五月伊人婷婷丁香| 成年人午夜在线观看视频| 久久毛片免费看一区二区三区| 国产真实伦视频高清在线观看| av卡一久久| 夜夜爽夜夜爽视频| 自拍偷自拍亚洲精品老妇| av卡一久久| 中文字幕久久专区| 婷婷色综合大香蕉| 下体分泌物呈黄色| 十八禁网站网址无遮挡 | 日韩欧美一区视频在线观看 | 又黄又爽又刺激的免费视频.| 国内少妇人妻偷人精品xxx网站| 99九九线精品视频在线观看视频| 国产高清三级在线| 少妇人妻一区二区三区视频| 夜夜骑夜夜射夜夜干| 国产精品国产av在线观看| 精品久久久久久久久亚洲| 一本久久精品| 99国产精品免费福利视频| 精品久久久久久电影网| 天堂俺去俺来也www色官网| 大香蕉久久网| √禁漫天堂资源中文www| 国产一级毛片在线| 99视频精品全部免费 在线| 国产亚洲午夜精品一区二区久久| 亚洲图色成人| 街头女战士在线观看网站| 精品国产国语对白av| 午夜视频国产福利| 久久99蜜桃精品久久| av国产久精品久网站免费入址| h视频一区二区三区| 亚洲熟女精品中文字幕| 简卡轻食公司| 夫妻性生交免费视频一级片| 成人漫画全彩无遮挡| 亚洲美女搞黄在线观看| 寂寞人妻少妇视频99o| 在线观看国产h片| 亚洲图色成人| 国产真实伦视频高清在线观看| 亚洲精品国产av成人精品| 国产精品一二三区在线看| 久久女婷五月综合色啪小说| 午夜久久久在线观看| 在线亚洲精品国产二区图片欧美 | 在线观看三级黄色| 99热这里只有精品一区| 国产极品天堂在线| 欧美 亚洲 国产 日韩一| 亚洲无线观看免费| 欧美精品国产亚洲| 久久精品国产鲁丝片午夜精品| 欧美精品一区二区大全| 亚洲国产最新在线播放| 内射极品少妇av片p| 爱豆传媒免费全集在线观看| 熟妇人妻不卡中文字幕| 日本与韩国留学比较| 欧美日韩国产mv在线观看视频| 国产av一区二区精品久久| 亚洲av成人精品一二三区| 精品人妻熟女av久视频| 大陆偷拍与自拍| 日韩强制内射视频| 99热这里只有是精品50| 亚洲怡红院男人天堂| 亚洲人与动物交配视频| 免费观看在线日韩| 国产精品久久久久久av不卡| 免费看不卡的av| 七月丁香在线播放| 18禁动态无遮挡网站| 日韩欧美 国产精品| 国产成人免费无遮挡视频| 久久婷婷青草| 各种免费的搞黄视频| 免费看av在线观看网站| 男女边摸边吃奶| 亚洲成人手机| 国产在线视频一区二区| 噜噜噜噜噜久久久久久91| 国产乱人偷精品视频| 这个男人来自地球电影免费观看 | 老女人水多毛片| 国产午夜精品久久久久久一区二区三区| 最近手机中文字幕大全| 久久久国产精品麻豆| 免费av中文字幕在线| 亚洲图色成人| 美女主播在线视频| 99热这里只有是精品在线观看| 日本免费在线观看一区| 卡戴珊不雅视频在线播放| 人妻 亚洲 视频| 精品一区二区免费观看| 丝袜在线中文字幕| 汤姆久久久久久久影院中文字幕| 丝袜脚勾引网站| 欧美最新免费一区二区三区| 最近中文字幕高清免费大全6| 青春草视频在线免费观看| 久久精品国产自在天天线| 狂野欧美白嫩少妇大欣赏| 最近中文字幕2019免费版| 美女xxoo啪啪120秒动态图| 国产亚洲av片在线观看秒播厂| 日本av免费视频播放| 曰老女人黄片| 国产成人午夜福利电影在线观看| 日本-黄色视频高清免费观看| 午夜视频国产福利| 精品熟女少妇av免费看| 免费观看a级毛片全部| 观看美女的网站| 大片免费播放器 马上看| 亚洲国产欧美日韩在线播放 | 亚洲美女黄色视频免费看| 校园人妻丝袜中文字幕| 少妇 在线观看| 日韩在线高清观看一区二区三区| 边亲边吃奶的免费视频| 精品午夜福利在线看| 美女cb高潮喷水在线观看| 国产日韩欧美在线精品| 欧美精品人与动牲交sv欧美| 我要看日韩黄色一级片| 热re99久久精品国产66热6| 国产成人精品久久久久久| 久久99热这里只频精品6学生| h视频一区二区三区| 中国三级夫妇交换| 国产欧美日韩一区二区三区在线 | 偷拍熟女少妇极品色| 日韩精品免费视频一区二区三区 | 女性生殖器流出的白浆| 久久97久久精品| 欧美高清成人免费视频www| 五月玫瑰六月丁香| 国产一区有黄有色的免费视频| 日韩不卡一区二区三区视频在线| 欧美日韩亚洲高清精品| av在线观看视频网站免费| 女性被躁到高潮视频| 久久久久久久久久久免费av| 久久影院123| 美女中出高潮动态图| 极品人妻少妇av视频| 欧美成人精品欧美一级黄| 最新的欧美精品一区二区| 我要看黄色一级片免费的| 精品酒店卫生间| 在线观看国产h片| 精品一区二区三卡| 国产成人91sexporn| 国产精品久久久久久久久免| 国产精品福利在线免费观看| 国产69精品久久久久777片| 日韩欧美 国产精品| 一区二区三区乱码不卡18| 在线精品无人区一区二区三| 久久久精品94久久精品| 少妇被粗大的猛进出69影院 | 国产精品一区二区性色av| 亚洲精品久久午夜乱码| av线在线观看网站| 水蜜桃什么品种好| 免费黄频网站在线观看国产| 高清黄色对白视频在线免费看 | 亚洲人与动物交配视频| 精品少妇黑人巨大在线播放| 精品国产一区二区久久| 久久人人爽av亚洲精品天堂| 热re99久久精品国产66热6| 国产精品.久久久| 国产伦精品一区二区三区视频9| 一级毛片久久久久久久久女| 国产精品伦人一区二区| 99国产精品免费福利视频| 一级毛片 在线播放| 日韩欧美精品免费久久| 自线自在国产av| av福利片在线观看| 亚洲美女视频黄频| av国产久精品久网站免费入址| 久久久国产欧美日韩av| 人妻夜夜爽99麻豆av| 久久午夜综合久久蜜桃| av线在线观看网站| 国产高清不卡午夜福利| 九色成人免费人妻av| 亚洲不卡免费看| 久久精品国产自在天天线| 亚洲精品456在线播放app| 在线观看免费日韩欧美大片 | 日日摸夜夜添夜夜添av毛片| 国产精品一区二区在线不卡| 人人妻人人爽人人添夜夜欢视频 | 亚洲丝袜综合中文字幕| 久久久国产一区二区| 欧美一级a爱片免费观看看| 一本色道久久久久久精品综合| 99久久人妻综合| 亚洲美女视频黄频| 中文字幕人妻熟人妻熟丝袜美| 亚洲国产av新网站| 精品亚洲乱码少妇综合久久| 亚洲综合色惰| 99热这里只有是精品50| 亚洲内射少妇av| 大又大粗又爽又黄少妇毛片口| 精品午夜福利在线看| 国产一区有黄有色的免费视频| 亚洲精品乱久久久久久| 日韩电影二区| 男女边吃奶边做爰视频| 在线播放无遮挡| 亚洲国产精品一区三区| 不卡视频在线观看欧美| 91精品伊人久久大香线蕉| 日韩熟女老妇一区二区性免费视频| 久久 成人 亚洲| 女人久久www免费人成看片| 91精品国产国语对白视频| 欧美激情极品国产一区二区三区 | 美女xxoo啪啪120秒动态图| 91精品一卡2卡3卡4卡| 国产日韩一区二区三区精品不卡 | 久久国产精品大桥未久av | 日韩av在线免费看完整版不卡| 国产精品久久久久久精品电影小说| av.在线天堂| 日本欧美国产在线视频| 免费观看av网站的网址| 亚洲欧美中文字幕日韩二区| 亚洲激情五月婷婷啪啪| 另类精品久久| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久精品电影小说| 欧美精品人与动牲交sv欧美| 亚洲精品第二区| 纵有疾风起免费观看全集完整版| 亚洲精品第二区| 国产黄频视频在线观看| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻一区二区三区麻豆| 在线观看三级黄色| 国产午夜精品一二区理论片| 青青草视频在线视频观看| 国产成人freesex在线| 校园人妻丝袜中文字幕| 国产伦精品一区二区三区四那| 亚洲av男天堂| 卡戴珊不雅视频在线播放| 精品一区二区三区视频在线| 亚洲,欧美,日韩| 国产一区二区在线观看av| a级毛色黄片| 国产永久视频网站| 六月丁香七月| a级一级毛片免费在线观看| 精品国产乱码久久久久久小说| 国产乱来视频区| 免费高清在线观看视频在线观看| 一本一本综合久久| 久久久久国产网址| 99九九线精品视频在线观看视频| 国产在线男女| 一级毛片 在线播放| 亚洲精品乱码久久久v下载方式| 日本午夜av视频| 国产中年淑女户外野战色| 国产永久视频网站| 最后的刺客免费高清国语| 大香蕉久久网| av在线播放精品| 国产精品人妻久久久久久| 亚洲av福利一区| 熟妇人妻不卡中文字幕| freevideosex欧美| 国产美女午夜福利| 一级毛片黄色毛片免费观看视频| 久久99热这里只频精品6学生| 在线观看一区二区三区激情| 国产精品福利在线免费观看| 国内精品宾馆在线| 国产伦在线观看视频一区| √禁漫天堂资源中文www| 91成人精品电影| 日本wwww免费看| 欧美成人精品欧美一级黄| 国产男女超爽视频在线观看| 又粗又硬又长又爽又黄的视频| 又爽又黄a免费视频| 欧美日韩精品成人综合77777| 最近2019中文字幕mv第一页| 亚洲电影在线观看av| 国内精品宾馆在线| 日日啪夜夜撸| 熟妇人妻不卡中文字幕| 国产亚洲91精品色在线| 欧美老熟妇乱子伦牲交| 午夜福利视频精品| 国产精品一区二区性色av| 久久人人爽人人片av| 日本av免费视频播放| 国产深夜福利视频在线观看| 久久人人爽人人爽人人片va| 国产伦在线观看视频一区| 毛片一级片免费看久久久久| 在线看a的网站| 国产乱来视频区| 日韩精品有码人妻一区| 亚洲第一av免费看| 噜噜噜噜噜久久久久久91| 九草在线视频观看| 日本欧美国产在线视频| 国产精品偷伦视频观看了| 国产爽快片一区二区三区| 18禁裸乳无遮挡动漫免费视频| 人妻 亚洲 视频| 精品少妇久久久久久888优播| 国产男女超爽视频在线观看| 国产成人aa在线观看| 国产午夜精品久久久久久一区二区三区| 欧美成人午夜免费资源| 免费久久久久久久精品成人欧美视频 | 亚洲av中文av极速乱| 嘟嘟电影网在线观看| 久久久久久久久大av| 在线观看免费视频网站a站| 一级毛片久久久久久久久女| 最近中文字幕2019免费版| 国产成人aa在线观看| 免费观看a级毛片全部| 精品亚洲成国产av| 成人影院久久| 夫妻性生交免费视频一级片| 亚洲人与动物交配视频| 视频区图区小说| 青春草视频在线免费观看| 一级毛片久久久久久久久女| 熟女电影av网| 国产精品国产三级国产专区5o| 亚洲欧美成人综合另类久久久| 国产日韩欧美在线精品| 国产综合精华液| 插逼视频在线观看| 麻豆精品久久久久久蜜桃| 亚洲怡红院男人天堂| 午夜免费观看性视频| 久久97久久精品| 大话2 男鬼变身卡| 国产一级毛片在线| 大码成人一级视频| 一个人看视频在线观看www免费| 一本久久精品| 国产精品偷伦视频观看了| 高清在线视频一区二区三区| 国产视频首页在线观看| 国产精品一区二区在线观看99| 国产精品久久久久久精品电影小说| 亚洲熟女精品中文字幕| 成人美女网站在线观看视频| 美女xxoo啪啪120秒动态图| 99精国产麻豆久久婷婷| 国产日韩欧美视频二区| 七月丁香在线播放| 中文资源天堂在线| av卡一久久| 久久免费观看电影| 久久精品国产a三级三级三级| 国产精品欧美亚洲77777| 五月开心婷婷网| 国产 精品1| 最黄视频免费看| 免费人成在线观看视频色| 久久久久久久精品精品| 一级毛片黄色毛片免费观看视频| 丰满饥渴人妻一区二区三| 日韩大片免费观看网站| 国产精品99久久久久久久久| 日本欧美国产在线视频| 69精品国产乱码久久久| 男人添女人高潮全过程视频| 日日摸夜夜添夜夜爱| 永久免费av网站大全| 亚洲国产欧美日韩在线播放 | 在线精品无人区一区二区三| 亚洲四区av| 街头女战士在线观看网站| 97在线人人人人妻| 国产高清国产精品国产三级| 欧美区成人在线视频| 亚洲精品成人av观看孕妇| 色婷婷av一区二区三区视频| 国产精品不卡视频一区二区| 春色校园在线视频观看| 久久国产乱子免费精品| 青春草视频在线免费观看| 高清不卡的av网站| 久久久久久久久久久久大奶| 国产深夜福利视频在线观看| 久久午夜综合久久蜜桃| 欧美日韩一区二区视频在线观看视频在线| 建设人人有责人人尽责人人享有的| 丰满人妻一区二区三区视频av| 尾随美女入室| 老司机亚洲免费影院| 久久精品国产自在天天线| 秋霞伦理黄片| 嫩草影院入口| 自拍偷自拍亚洲精品老妇| 99热这里只有是精品50| 久久精品久久精品一区二区三区| 欧美精品人与动牲交sv欧美| 免费黄网站久久成人精品| 亚洲精品国产成人久久av| 亚洲成色77777| 国产乱来视频区| 日韩中文字幕视频在线看片| 久久久久国产网址| 欧美xxxx性猛交bbbb| 少妇人妻 视频| a级片在线免费高清观看视频| 亚洲精品国产av成人精品| 亚洲精品日韩av片在线观看| 日韩一区二区视频免费看| 久久午夜福利片| 日产精品乱码卡一卡2卡三| 国产精品麻豆人妻色哟哟久久| 久久精品久久精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 久久精品国产鲁丝片午夜精品| 美女xxoo啪啪120秒动态图| av.在线天堂| 三上悠亚av全集在线观看 | 国产黄片美女视频| a级毛色黄片| 欧美精品国产亚洲| 亚洲av二区三区四区| 日本黄色片子视频| 2018国产大陆天天弄谢| 国产男女超爽视频在线观看| 超碰97精品在线观看| 一级a做视频免费观看| 亚洲成色77777| 国产69精品久久久久777片| 国产国拍精品亚洲av在线观看| 欧美最新免费一区二区三区| 亚洲av综合色区一区| av又黄又爽大尺度在线免费看| 亚洲欧美日韩东京热| 午夜激情福利司机影院| 日韩中文字幕视频在线看片| 久久久久久久久久久免费av| 久久99热6这里只有精品| 日韩av免费高清视频| 亚洲av二区三区四区| 亚洲av综合色区一区| 2018国产大陆天天弄谢| 99精国产麻豆久久婷婷| 欧美 日韩 精品 国产| 狂野欧美激情性xxxx在线观看| 成年美女黄网站色视频大全免费 |