頁(yè)巖氣藏多場(chǎng)耦合壓裂水平井產(chǎn)能預(yù)測(cè)

1、目錄1頁(yè)巖氣藏國(guó)內(nèi)外研究現(xiàn)狀21.1頁(yè)巖氣開(kāi)發(fā)利用現(xiàn)狀21.2頁(yè)巖氣數(shù)值模擬現(xiàn)狀32氣藏壓裂水平井穩(wěn)定產(chǎn)能預(yù)測(cè)模型72.1氣藏一裂縫滲流模型72.2水平井井筒壓降模型82.3穩(wěn)定產(chǎn)能預(yù)測(cè)藕合模型的建立及其求解103頁(yè)巖氣藏產(chǎn)能動(dòng)態(tài)分析9-11123.1吸附規(guī)律模擬及對(duì)產(chǎn)能影響研究123.1.1模型描述及基本參數(shù)123.1.2吸附氣對(duì)產(chǎn)能影響研究133.1.3等溫吸附規(guī)律對(duì)產(chǎn)能影響研究143.1.4吸附氣含量對(duì)產(chǎn)能影響研究163.2儲(chǔ)層滲透率對(duì)產(chǎn)能影響研究183.2.1模型描述及基本參數(shù)183.2.2基巖滲透率對(duì)產(chǎn)能影響研究203.2.3裂縫滲透率對(duì)產(chǎn)能影響研究21結(jié)論22參考文獻(xiàn)231頁(yè)巖氣藏

2、國(guó)內(nèi)外研究現(xiàn)狀頁(yè)巖氣主要存在于富有機(jī)質(zhì)的泥頁(yè)巖及其夾層中，包括吸附氣及游離氣，吸附氣以吸附狀態(tài)、游離氣以游離狀態(tài)存在。頁(yè)巖不僅是生成頁(yè)巖氣的源巖，又是匯集和儲(chǔ)存頁(yè)巖氣的儲(chǔ)層及蓋層。與常規(guī)天然氣不同的是頁(yè)巖氣田開(kāi)采壽命一般可達(dá)30年50年，甚至更長(zhǎng)，因此頁(yè)巖氣的優(yōu)點(diǎn)是開(kāi)采時(shí)間長(zhǎng)及產(chǎn)氣周期長(zhǎng)。并且多數(shù)的源巖分布廣厚度大，而且含頁(yè)巖氣量較高，導(dǎo)致多數(shù)頁(yè)巖氣井能夠保證穩(wěn)定的長(zhǎng)期生產(chǎn)。因此，能否有效合理的對(duì)頁(yè)巖氣進(jìn)行開(kāi)發(fā)生產(chǎn)，對(duì)保持中國(guó)天然氣供給的持續(xù)穩(wěn)定發(fā)展有著十分重要的意義。頁(yè)巖氣藏的儲(chǔ)層一般具有低孔隙度、低滲透率的特征，因而生產(chǎn)率和采收率較低，常規(guī)天然氣采收率60%以上，而頁(yè)巖氣僅為5%55%，不

3、能達(dá)到經(jīng)濟(jì)收益的要求。為提高產(chǎn)量使得開(kāi)采獲得的經(jīng)濟(jì)效益最大化，壓裂工藝則是現(xiàn)今開(kāi)采低滲透油氣藏的關(guān)鍵改造方法之一。評(píng)價(jià)是否為有效壓裂的主要看經(jīng)過(guò)壓裂的生產(chǎn)井滲流規(guī)律及產(chǎn)量變化規(guī)律，進(jìn)而為開(kāi)發(fā)頁(yè)巖氣藏提供了原理依據(jù)，所以對(duì)頁(yè)巖氣井壓裂模擬進(jìn)行深入的研究是十分有必要的。根據(jù)數(shù)據(jù)信息推測(cè)得知，中國(guó)頁(yè)巖氣藏的主要儲(chǔ)集區(qū)域是南部，主要為海相頁(yè)巖地層。除了海相頁(yè)巖地層，陸相沉積盆地的頁(yè)巖地層也具備頁(yè)巖氣儲(chǔ)集的基本條件，主要區(qū)域在松遼、鄂爾多斯、吐哈和準(zhǔn)噶爾等地。我國(guó)主要盆地和地區(qū)頁(yè)巖氣資源量為（1530）×1012m3與美國(guó)28.3 ×1012m3相當(dāng)，經(jīng)濟(jì)價(jià)值巨大，美國(guó)正直頁(yè)巖氣藏的

4、快速開(kāi)發(fā)階段，而中國(guó)只是剛剛起步，至今仍處于對(duì)成藏機(jī)理和勘探開(kāi)發(fā)的研究階段。因此對(duì)于頁(yè)巖氣的數(shù)值模擬研究顯得尤為重要，本文在有限元法基本原理的基礎(chǔ)上，考慮了頁(yè)巖氣吸附解吸特性，并對(duì)頁(yè)巖氣進(jìn)行三維壓裂數(shù)值模擬研究。1.1頁(yè)巖氣開(kāi)發(fā)利用現(xiàn)狀全球頁(yè)巖氣資源總量約456×1012m3時(shí)，相當(dāng)于其他兩種非常規(guī)天然氣一煤層氣和致密砂巖氣儲(chǔ)量的總和，世界主要地區(qū)頁(yè)巖氣儲(chǔ)量情況見(jiàn)表1-1。1821年的第一口工業(yè)性頁(yè)巖氣鉆井標(biāo)志著美國(guó)進(jìn)入頁(yè)巖氣的生產(chǎn)研究階段。隨著技術(shù)的成熟，美國(guó)頁(yè)巖氣年產(chǎn)量己從2007年的340×1012m3躍升至2010年的超過(guò)900×1012m3時(shí)，到2015

5、年預(yù)計(jì)年產(chǎn)量將達(dá)到1800×1012m3。對(duì)于美國(guó)天然氣生產(chǎn)預(yù)計(jì)增長(zhǎng)額來(lái)說(shuō)，頁(yè)巖氣做出了最大的貢獻(xiàn)，預(yù)計(jì)到2035年其頁(yè)巖氣生產(chǎn)量占天然氣生產(chǎn)總量的一半。加拿大占據(jù)了將近一半的北美天然氣市場(chǎng)份額，在北美發(fā)揮著重要的作用。其國(guó)內(nèi)的天然氣主要產(chǎn)區(qū)是不列顛哥倫比亞省，而在過(guò)去的10年里，該省天然氣生產(chǎn)的增長(zhǎng)額幾乎都是來(lái)自于致密氣和頁(yè)巖氣。墨西哥的頁(yè)巖氣生產(chǎn)具有潛力，但現(xiàn)處早期開(kāi)發(fā)階段。2011年2月，墨西哥北部的科阿韋拉州開(kāi)采出了第一批頁(yè)巖氣，并且這口頁(yè)巖氣井能夠口產(chǎn)量為29×105ft3，并計(jì)劃再鉆10口井，目標(biāo)是開(kāi)發(fā)天然氣和凝析油。續(xù)美國(guó)和加拿大以外，印度是第一個(gè)從頁(yè)巖中發(fā)現(xiàn)

6、天然氣的亞洲國(guó)家。同為亞洲國(guó)家的我國(guó)，具有的頁(yè)巖氣資源前景和開(kāi)發(fā)潛力與美國(guó)大致體相等。一些學(xué)者認(rèn)為：由于我國(guó)與美國(guó)的頁(yè)巖氣成藏的因素及地層條件大體相同，并且我國(guó)主要成藏地區(qū)的氣體儲(chǔ)量預(yù)算為26×1012m3時(shí)，與美國(guó)的28.3×1012m3可以相比，若以2010年美國(guó)1000×108m3時(shí)的頁(yè)巖氣產(chǎn)量為例進(jìn)行比較計(jì)算，然而我國(guó)目前的天然氣供需缺量?jī)H為300×108m3了時(shí)，也就是說(shuō)我國(guó)頁(yè)巖氣產(chǎn)量只需是美國(guó)的三分之一就可以填補(bǔ)國(guó)內(nèi)供需缺口。然而當(dāng)前我國(guó)頁(yè)巖氣發(fā)展仍面臨很多問(wèn)題及挑戰(zhàn)。雖然我國(guó)頁(yè)巖氣儲(chǔ)量豐富，但是由于我國(guó)頁(yè)巖氣開(kāi)發(fā)較晚，造成頁(yè)巖氣資源的勘能力

7、及水準(zhǔn)較低且技術(shù)水平不成體系。加強(qiáng)科研人員對(duì)頁(yè)巖氣開(kāi)采的進(jìn)一步研究及對(duì)外國(guó)先進(jìn)技術(shù)的學(xué)習(xí)，從而合理有效開(kāi)發(fā)頁(yè)巖氣新能源成為擺在我們面前的嚴(yán)峻課題。1.2頁(yè)巖氣數(shù)值模擬現(xiàn)狀國(guó)內(nèi)外以頁(yè)巖氣滲流理論為基礎(chǔ)的數(shù)值模擬方面的研究己建立的數(shù)學(xué)模型大致可以分為三類：經(jīng)驗(yàn)?zāi)Ｐ?、平衡吸附模型、非平衡吸附模型?)經(jīng)驗(yàn)?zāi)Ｐ蚅idine模型、Aireg模型及McFall模型等均為經(jīng)驗(yàn)?zāi)Ｐ?，其模型主要是?duì)可觀察的物理現(xiàn)象進(jìn)行簡(jiǎn)要的數(shù)學(xué)描述，在模型應(yīng)用時(shí)輸入的參數(shù)較少，但此模型理論基礎(chǔ)不足，精度較低。2)平衡吸附模型Kissell模型、Mckee模型及Bumb模型等為平衡吸附模型，其中最典型的為Bumb模型，此模型假設(shè)

8、儲(chǔ)層介質(zhì)為是單孔隙介質(zhì)，吸附氣在頁(yè)巖氣藏?cái)?shù)值模擬不中能被忽略，采用蘭氏等溫吸附方程對(duì)吸附頁(yè)巖氣進(jìn)行描述，并假設(shè)氣體擴(kuò)散是瞬間完成的，即吸附氣與游離氣壓力時(shí)刻保持平衡狀態(tài)。但模型由于忽略了吸附頁(yè)巖氣的解吸過(guò)程，所以不能反映客觀存在的解吸時(shí)間，其預(yù)測(cè)的產(chǎn)量高于實(shí)際產(chǎn)量。3)非平衡吸附模型頁(yè)巖氣的吸附解吸、擴(kuò)散和滲流為一個(gè)相互影響相互制約的整體過(guò)程，擴(kuò)散模型認(rèn)為擴(kuò)散不可忽略。非平衡模型又可分為基于菲克第一擴(kuò)散定律的擬穩(wěn)態(tài)模型(包括Psu-1模型、Psu-2模型、Psu-3模型及Comet模型等)以及基于菲克第二擴(kuò)散定律的非穩(wěn)態(tài)模型(包括Smith模型、Sugarwat模型、Chen模型等)。非平衡吸

9、附模型的假設(shè)條件為：1)頁(yè)巖氣藏由基巖一裂縫雙重介質(zhì)組成；2)頁(yè)巖氣藏中氣水均是微可壓縮的；3)還考慮了頁(yè)巖氣的吸附和游離兩種形態(tài)，其中裂縫中僅存在游離態(tài)的頁(yè)巖氣，而基巖中不僅存在游離態(tài)頁(yè)巖氣，還存在部分吸附態(tài)頁(yè)巖氣；4)開(kāi)采后地層壓力變化引起吸附氣解吸并擴(kuò)散到裂縫中的擴(kuò)散過(guò)程?；贔ick擴(kuò)散定律，建立了基質(zhì)中氣體的解吸擴(kuò)散的非平衡模型，再與裂縫中流體的滲流模型禍合。非平衡吸附擴(kuò)散模型較好的反映了頁(yè)巖氣在頁(yè)巖儲(chǔ)層中的賦存及運(yùn)移機(jī)理，較平衡吸附模型有所進(jìn)步。Bustin等建立了氣一水兩相雙重介質(zhì)模型，氣井的產(chǎn)能和氣體飽和度等很多因素有關(guān)。含氣飽和度、氣體吸附能力及地層壓力都與氣藏含氣量成正比，

10、從而增加生產(chǎn)能力。頁(yè)巖基質(zhì)滲透率也是影響頁(yè)巖氣生產(chǎn)的因素之一，其主要與頁(yè)巖基質(zhì)結(jié)構(gòu)有關(guān)。Bustin對(duì)部分因素進(jìn)行對(duì)比，總結(jié)了各因素對(duì)頁(yè)巖氣藏產(chǎn)能的影響，最終發(fā)現(xiàn)擴(kuò)散系數(shù)及應(yīng)力敏感裂縫滲透率對(duì)頁(yè)巖氣藏產(chǎn)能影響較大。John Lee等應(yīng)用蘭氏等溫吸附定律在Warren&Root提出的經(jīng)典雙重介質(zhì)模型中，得到了考慮吸附的雙重介質(zhì)模型并對(duì)其求解，最終繪制了基于等溫吸附方程的雙重介質(zhì)模型理論圖版。頁(yè)巖氣主要吸附在頁(yè)巖基質(zhì)表面，而且頁(yè)巖氣儲(chǔ)層一般都具有天然裂縫，用經(jīng)典雙重介質(zhì)模型是合適的。由于吸附態(tài)頁(yè)巖氣主要吸附在基質(zhì)顆粒表面，且頁(yè)巖大多具有天然裂縫，所以用此經(jīng)典雙重介質(zhì)模型是較適合的。John

11、 Lee還提出了竄流系數(shù)，此系數(shù)是在串流系數(shù)中反映吸附效應(yīng)的一個(gè)變量。最后John Lee發(fā)現(xiàn)在不同的層中，吸附解吸是同時(shí)并且可逆發(fā)生的。李斌建立考慮孔隙中非平衡吸附氣解吸及擴(kuò)散過(guò)程的氣一水兩相裂縫滲流問(wèn)題的數(shù)學(xué)模型，對(duì)頁(yè)巖氣產(chǎn)量和產(chǎn)水率隨時(shí)間的變化給出了較好的預(yù)測(cè)。Watson1,2等建立了理想雙重介質(zhì)模型，并確定了此模型的適用性，提出一種新方法來(lái)分析天然裂縫氣藏的生產(chǎn)數(shù)據(jù)。Carlson等認(rèn)為由于頁(yè)巖的致密性和低滲透性，氣體分子擴(kuò)散流動(dòng)是頁(yè)巖氣在頁(yè)巖中的流動(dòng)主要方式，即便不考慮頁(yè)巖氣藏非常規(guī)的儲(chǔ)集方式和流動(dòng)性，僅考慮頁(yè)巖孔隙的致密性，泥盆紀(jì)頁(yè)巖也可視為特殊的雙重介質(zhì)模型，該模型的優(yōu)點(diǎn)是：對(duì)

12、于短期的預(yù)測(cè)頁(yè)巖氣產(chǎn)能結(jié)果很精準(zhǔn)，但該模型的缺點(diǎn)是：長(zhǎng)期的預(yù)測(cè)產(chǎn)能結(jié)果會(huì)偏小。然而對(duì)于Watson和Carlson提出的兩個(gè)模型，共同點(diǎn)是考慮的因素較理想，忽略了吸附氣的對(duì)產(chǎn)能及儲(chǔ)量的巨大影響，假設(shè)氣體只以游離態(tài)為存在狀態(tài)，進(jìn)而認(rèn)為氣體擴(kuò)散為頁(yè)巖氣藏的主要運(yùn)移方式，并且都采取雙重介質(zhì)模型對(duì)氣體進(jìn)行產(chǎn)能預(yù)測(cè)分析Ozkan等考慮到頁(yè)巖的致密性，認(rèn)為基質(zhì)中孔隙為納米級(jí)的，頁(yè)巖氣在基質(zhì)中的擴(kuò)散流動(dòng)發(fā)生在納米孔隙中，而在裂縫中流動(dòng)時(shí)則考慮裂縫滲透率的應(yīng)力敏感性，最終建立了非達(dá)西滲流的擴(kuò)散機(jī)制的基巖一裂縫雙重介質(zhì)模型。在前人研究的基礎(chǔ)上，一些研究者又建立了水平層狀雙重非均勻介質(zhì)的擬穩(wěn)態(tài)和非穩(wěn)態(tài)模型，對(duì)含氣

13、層與上、下層都有水交換的情形，采用擬三維模型，并考慮了地層壓裂的情況，發(fā)展了多維多組分頁(yè)巖氣生產(chǎn)模型，三孔隙度雙滲透率的三重介質(zhì)模型Zhang X3等用軟件Eclipse通過(guò)對(duì)蘭氏等溫吸附定律進(jìn)行擴(kuò)展，得到可模擬多組分吸附氣解吸的數(shù)學(xué)模型。并進(jìn)一步研究了一些數(shù)據(jù)對(duì)頁(yè)巖氣藏產(chǎn)能的影響，主要是氣藏地層參數(shù)和人工壓裂參數(shù)，模擬了地層參數(shù)及壓裂系數(shù)對(duì)頁(yè)巖氣產(chǎn)能的影響，并對(duì)這些參數(shù)進(jìn)行了敏感性分析，參數(shù)主要包括：基巖的孔隙度及滲透率、裂縫的孔隙度及滲透率、基巖一裂縫的禍合參數(shù)、基巖分塊及裂縫的半長(zhǎng)、間距及高度、巖石壓縮系數(shù)以及氣體體積含量。結(jié)合頁(yè)巖氣藏開(kāi)采現(xiàn)場(chǎng)情況，部分學(xué)者發(fā)現(xiàn)頁(yè)巖氣藏基巖的致密性使開(kāi)采

14、難度增大，壓裂技術(shù)成為頁(yè)巖氣藏有效開(kāi)采必要手段，然而天然形成裂縫與人工壓裂裂縫相比較，無(wú)論在裂縫形態(tài)還是裂縫的導(dǎo)流能力上都有很大差別，因此一些研究人員對(duì)頁(yè)巖氣藏產(chǎn)能進(jìn)行研究時(shí)采用多重介質(zhì)模型，多重介質(zhì)模型即基巖一小裂縫一大裂縫介質(zhì)模型。而對(duì)于多重介質(zhì)而言，依據(jù)擴(kuò)散氣體分子運(yùn)動(dòng)的平均自由程的長(zhǎng)短和基質(zhì)顆?？椎赖拇笮》譃椋悍瓶藬U(kuò)散、努森擴(kuò)散、表面擴(kuò)散和晶體擴(kuò)散。若孔道的直徑遠(yuǎn)大于擴(kuò)散分子的平均自由程時(shí)，以菲克擴(kuò)散為主要擴(kuò)散方式。當(dāng)孔道的直徑比擴(kuò)散分子的平均自由程小的時(shí)候，則以努森擴(kuò)散為主要擴(kuò)散方式，同時(shí)，隨著擴(kuò)散的氣體分子與孔道壁之間的碰撞，吸附的氣體分子還會(huì)沿著孔道壁表面移動(dòng)而形成表面擴(kuò)散。氣體

15、擴(kuò)散的決定因素主要是：孔隙形狀、大小、連通性、多元?dú)怏w性質(zhì)和狀態(tài)，擴(kuò)散系數(shù)越大，擴(kuò)散速率越高。因?yàn)轫?yè)巖氣藏有機(jī)物基質(zhì)與無(wú)機(jī)物基質(zhì)孔徑分布范圍較大，菲克擴(kuò)散、努森擴(kuò)散、表面擴(kuò)散和晶體擴(kuò)散在頁(yè)巖氣基質(zhì)中可能同時(shí)存在，但以努森擴(kuò)散為主4。Wu Yushu等將裂縫分為人工裂縫和天然裂縫，考慮應(yīng)力敏感和克林肯伯格效應(yīng)，建立了裂縫氣藏多重致密介質(zhì)模型?；鶐r中考慮克林肯伯格效應(yīng)，裂縫中考慮慣性作用的高速非達(dá)西流。分析了克林肯伯格效應(yīng)對(duì)產(chǎn)能的影響，對(duì)比了雙重介質(zhì)模型和多重介質(zhì)模型的區(qū)別。通過(guò)油藏?cái)?shù)值模擬軟件模擬了一些壓裂參數(shù)對(duì)頁(yè)巖氣產(chǎn)能的影響，參數(shù)包括：裂縫導(dǎo)流能力、裂縫間距及解吸等。Cipolla等為如何改

16、進(jìn)壓裂設(shè)計(jì)、完井方式和油田開(kāi)發(fā)方案來(lái)提高最終采收率和產(chǎn)氣速度提供了理論指導(dǎo)。Kalantari-Dahaghi模擬了二氧化碳驅(qū)開(kāi)發(fā)頁(yè)巖氣的可行性。用DFN模型對(duì)不確定的參數(shù)進(jìn)行有效計(jì)算及校正。模型中考慮人工裂縫，同時(shí)對(duì)包含天然裂縫和人工裂縫模型的基巖和裂縫的關(guān)鍵參數(shù)進(jìn)行敏感性分析，對(duì)頁(yè)巖氣產(chǎn)能最敏感的因素指標(biāo)有：吸附類型選擇、天然裂縫滲透率、基質(zhì)一裂縫連接因子和水力壓裂裂縫參數(shù)(包括裂縫間距、裂縫高度、裂縫半長(zhǎng)和裂縫導(dǎo)流能力)。在參考基礎(chǔ)地質(zhì)資料的基礎(chǔ)上，可對(duì)這些參數(shù)在合理范圍內(nèi)進(jìn)行調(diào)整以用于歷史擬合研究。最終得出結(jié)論人工壓裂的處理是成功開(kāi)發(fā)頁(yè)巖氣藏的關(guān)鍵所在。不同于Wu Yushu的多重介質(zhì)

17、模型，Moridis等建立了考慮多組分吸附的頁(yè)巖氣等效介質(zhì)模型，模型考慮克林肯伯格效應(yīng)，并且認(rèn)為擴(kuò)散的影響不可忽略，模型假設(shè)在介質(zhì)中，氣體的流動(dòng)是達(dá)西流動(dòng)或高速非達(dá)西流動(dòng)同時(shí)存在的，氣體在頁(yè)巖中吸附滿足蘭氏等溫吸附定律，用氣水狀態(tài)方程的TOUGH+模擬器模擬頁(yè)巖氣流動(dòng)，比較等溫和非等溫、達(dá)西流和高速非達(dá)西流的差別，研究了吸附曲線、裂縫類型對(duì)產(chǎn)能的影響，并對(duì)比等效介質(zhì)、雙孔隙和雙滲透率3種模型模擬頁(yè)巖氣流動(dòng)時(shí)的差別，研究表明雙滲透率模型與實(shí)際擬合最好，等溫和非等溫情況下模型預(yù)測(cè)的結(jié)果相差不大。2氣藏壓裂水平井穩(wěn)定產(chǎn)能預(yù)測(cè)模型利用勢(shì)的疊加原理和流體力學(xué)理論，建立穩(wěn)定和不穩(wěn)定條件下氣藏一裂縫內(nèi)流體滲

18、流與水平井井筒內(nèi)流體流動(dòng)的禍合模型，氣藏模型考慮水平井的三維滲流特征，裂縫模型考慮裂縫內(nèi)的壓力損失，井筒壓降模型考慮井筒壁面摩擦和裂縫流體流入所引起的壓力降落等因素的影響，并利用數(shù)值解法對(duì)模型進(jìn)行求解，同時(shí)分析各個(gè)參數(shù)對(duì)壓裂水平井產(chǎn)能的影響，這對(duì)于壓裂水平井的生產(chǎn)具有重要的指導(dǎo)意義。2.1氣藏一裂縫滲流模型假設(shè)：（1）非均質(zhì)儲(chǔ)層上下邊界為不滲透邊界，水平方向半徑為Re，厚度為h，水平滲透率為Kh，地層中心有一口水平井，長(zhǎng)度為L(zhǎng)；（2）在水平段壓開(kāi)N條垂直裂縫，裂縫等距離分布并且貫穿整個(gè)儲(chǔ)層厚度，裂縫半長(zhǎng)為xf，裂縫滲透率為Kf，寬度為Wf。（3）水平段沒(méi)有進(jìn)行補(bǔ)孔，流體將先從地層流入裂縫，然后

19、沿裂縫進(jìn)入井筒，裂縫的產(chǎn)量為qgf，壓裂水平井產(chǎn)量為Qgf，壓裂水平井的產(chǎn)量即為各條裂縫產(chǎn)量之和。水平段井筒中的流動(dòng)方式包括井筒內(nèi)流體的主軸流動(dòng)和裂縫內(nèi)流體向井筒的徑向流動(dòng)。8假定氣體在水平井筒內(nèi)作等溫穩(wěn)定單相流動(dòng)，根據(jù)流體力學(xué)理論，由于壁面摩擦和流體匯流的影響，水平井筒內(nèi)存在一定的壓力損失。當(dāng)裂縫條數(shù)N為奇數(shù)時(shí)，裂縫中點(diǎn)橫坐標(biāo)分別為0，士d，士2d，士N0d，其中No= (N-1)/2，為y軸左右裂縫的條數(shù)，d=L/N，為裂縫與裂縫之間的間距。根據(jù)復(fù)位勢(shì)理論和勢(shì)的疊加原理，并根據(jù)真實(shí)氣體的狀態(tài)方程并將壓裂水平井產(chǎn)量換算為地面標(biāo)準(zhǔn)狀況下氣體產(chǎn)量，則第j條裂縫的勢(shì)近似為fj(md，0)：其中，m

20、=-No+j- l。在x軸上距原點(diǎn)較遠(yuǎn)處取點(diǎn)(Re，0)，R。為供給半徑，則供給邊界的勢(shì)e為由式(2-1)和式(2-2整理可得到當(dāng)裂縫條數(shù)N為偶數(shù)時(shí)，為了計(jì)算的方便，取d=L/2N，取N0=N-1，則每條裂縫對(duì)應(yīng)的橫坐標(biāo)為士d，士3d，士Nod。同樣可得到其中，i從一N0的速度遞增，j=1，2，N，m=- N0+2 (j-1)。由于裂縫的半長(zhǎng)遠(yuǎn)遠(yuǎn)大于水平井筒的半徑，所以裂縫內(nèi)的流體從裂縫邊緣向井筒周圍的聚集，如果忽略重力的影響，可近似看作地層厚度為Wf，流動(dòng)半徑h/2，邊界壓力為pfj的平面徑向流，不考慮裂縫的表皮因子造成的壓降，則有如下表達(dá)式成立將式（2-3）+式（2-5）或式（2-4）+（

21、2-5）得由，得進(jìn)行整理，得2.2水平井井筒壓降模型圖2.1為壓裂水平井井筒示意圖，取井筒上第j條裂縫和第j+1條裂縫之間的一小段進(jìn)行研究，第j條裂縫左端的出口壓力為p1j、出口速度為v1J，右端的進(jìn)口壓力為pzj、進(jìn)口速度為v2j。氣體從第j+1條裂縫的左端流到第j條裂縫左端的過(guò)程中，考慮井筒內(nèi)壓降，由動(dòng)量定理得其中，m為質(zhì)量流量，m=pAv，p為氣體的密度，A為井筒橫截面面積，為井筒壁面剪切應(yīng)力。因?yàn)榱芽p內(nèi)的氣體徑向流入水平井筒，所以有對(duì)式(2-9)進(jìn)行整理得其中，壁面剪切應(yīng)力造成的摩擦壓力降，記為將代入式（2-12）得式(2-11)右邊最后一項(xiàng)可認(rèn)為是由于動(dòng)量變化造成的加速壓力降，記為由

22、于忽略了除裂縫之外的其它井筒部分氣體的流入，此壓力降只有在裂縫處才會(huì)產(chǎn)生，所以上式又可寫為根據(jù)以上的分析可知，井筒內(nèi)的壓降可分為兩部分：摩擦壓力降和加速壓力降。下面對(duì)每條裂縫兩段的壓力p1j和p2j進(jìn)行具體求解。將式(2-13)寫成微分形式氣體的密度的表達(dá)式為利用氣體的狀態(tài)方程可得到速度v的表達(dá)式所以式(2-16)可以寫成對(duì)式(2-19)在第(j+1)條裂縫左端到第j條裂縫右端之間一段井筒距離內(nèi)進(jìn)行積分，可得其中，當(dāng)N為奇數(shù)時(shí)，；當(dāng)N為偶數(shù)時(shí)，。將和v的表達(dá)式代入式(2-15)得在末端裂縫處流體的流動(dòng)可看作是在一彎曲的管子中的流動(dòng)，不會(huì)產(chǎn)生因?yàn)橘|(zhì)量的增加而引起的速度變化，所以有第j段井筒壁面的

23、摩擦系數(shù)fj則可用下式計(jì)算得到其中，。當(dāng)Rej<2000時(shí)，屬于層流；當(dāng)Rej>4000時(shí)，屬于紊流；2000<Rej<4000時(shí)，屬于過(guò)渡流，f取上述兩者的加權(quán)平均值，即當(dāng)屬過(guò)渡流時(shí)其中，f1、f2分別為層流和紊流摩擦系數(shù)，為權(quán)值，一般取0.10.3即可。2.3穩(wěn)定產(chǎn)能預(yù)測(cè)藕合模型的建立及其求解 利用流量和壓力的禍合關(guān)系把氣藏模型和井筒模型連接起來(lái)得到禍合模型。壓裂水平井第j段井筒內(nèi)氣體的流量Qj的計(jì)算方法如下式(2-8)中的第j條裂縫底部的井筒壓力近似取為將計(jì)算摩擦壓力降的式(2-20)和加速壓力降的式(2-21)、式(2-22)與式(2-8)聯(lián)立，由于沒(méi)有引入新的

24、變量，這樣就會(huì)產(chǎn)生一個(gè)有3N個(gè)未知數(shù)、3N個(gè)方程的方程組，所以此方程組可以封閉求解。由于在上述的方程組中qgtj、Plj和P2j為未知數(shù)，共3N個(gè)，且它們之間為復(fù)雜的非線性關(guān)系，可采取迭代方法求解，即先假定一組裂縫的產(chǎn)量的初值代入式(2-20、式(2-21)和式(2-22)中，從第1條裂縫開(kāi)始進(jìn)行計(jì)算，交替利用式(2-20)、式(2-21)和式(2-22)，逐次計(jì)算出Plj和P2j的一組值：P11， P12，P1N，P1N。然后再將Plj和P2j的值代入到式(2-8)中計(jì)算qgtj，如此反復(fù)循環(huán)，直到滿足一定精度為止，最后可以求出氣井的總產(chǎn)量5,6Qgf： 3頁(yè)巖氣藏產(chǎn)能動(dòng)態(tài)分析9-11吸附氣

25、在頁(yè)巖氣藏的滲流模擬過(guò)程中是不可忽視的，本章通過(guò)對(duì)頁(yè)巖氣吸附規(guī)律的分析對(duì)比，總結(jié)出吸附氣對(duì)頁(yè)巖氣藏產(chǎn)能的影響規(guī)律。壓裂技術(shù)是頁(yè)巖氣藏開(kāi)采過(guò)程中經(jīng)常用到的手段，本章通過(guò)對(duì)基巖及裂縫的滲透率對(duì)比分析，總結(jié)出裂縫滲透率及基巖滲透率對(duì)產(chǎn)能影響規(guī)律。3.1吸附規(guī)律模擬及對(duì)產(chǎn)能影響研究3.1.1模型描述及基本參數(shù)通過(guò)借鑒國(guó)外頁(yè)巖儲(chǔ)層典型參數(shù)建立了頁(yè)巖氣藏三維模型。設(shè)頁(yè)巖氣藏地層模型長(zhǎng)1210m，寬810m，頁(yè)巖儲(chǔ)層深度2000m，有效厚度30m，模擬時(shí)間10a，壓裂采用體積壓裂方法，縫網(wǎng)尺寸為600m×200m，計(jì)算參數(shù)見(jiàn)表3.1。表3.1模擬計(jì)算的有關(guān)參數(shù)計(jì)算參數(shù)參數(shù)值計(jì)算參數(shù)參數(shù)值原始地層溫

26、度（K）361滲透率（10-3m2）0.001氣體摩爾質(zhì)（kg/mol）0.016孔隙度（%）5氣體常數(shù)（Pa·m3/(mol·k)）8.314標(biāo)況下氣體密度（kg/m3）0.717內(nèi)邊界壓力（MPa）4原始地層壓力（MPa）20蘭氏體積（m3/m3）13氣體粘度（mPa.s）0.067蘭氏壓力（MPa）3壓縮系數(shù)（MPa-1）1×10-4總地質(zhì)儲(chǔ)量（m3）5.43×108井筒半徑（m）0.1采用混合單元網(wǎng)格進(jìn)行剖分，模擬時(shí)將裂縫視為一個(gè)面，用三角形面單元進(jìn)行剖分，地層用任意四面體單元剖分。模型網(wǎng)格剖分圖如圖3.1所示，網(wǎng)格加密區(qū)為采用體積壓裂區(qū)域。(a

27、)頂部平面網(wǎng)格圖 (b)立體網(wǎng)格剖分圖圖3.1網(wǎng)格剖分圖圖3.2分別為模擬模型0.5年、2.5年、5年、8年的壓力分布圖，可以看出生產(chǎn)過(guò)程中各時(shí)間的壓力分布情況7。(a)0.5年時(shí)的壓力分布圖 (b)2.5年時(shí)的壓力分布圖(c)5年時(shí)的壓力分布圖 (d)8年時(shí)的壓力分布圖圖3.2壓力分布圖3.1.2吸附氣對(duì)產(chǎn)能影響研究圖3.3吸附氣對(duì)日產(chǎn)氣量影響曲線圖3.4吸附氣對(duì)累積產(chǎn)氣量影響曲線為了對(duì)考慮吸附氣與不考慮吸附氣進(jìn)行比較，建立了考慮與不考慮吸附氣2個(gè)方案，見(jiàn)圖3.3，圖3.4，從圖中可以看出：頁(yè)巖氣的產(chǎn)能主要來(lái)源于游離氣，但是解吸氣所占的比例也不可忽視，頁(yè)巖氣井的產(chǎn)量是游離氣生產(chǎn)和吸附氣解吸共

28、同作用的結(jié)果。由圖3.3可以看出：頁(yè)巖氣的穩(wěn)產(chǎn)時(shí)間較長(zhǎng)，后期生產(chǎn)速度逐漸減小。進(jìn)一步比較考慮吸附和不考慮吸附的兩個(gè)模型(即產(chǎn)氣量曲線及游離氣產(chǎn)量曲線)計(jì)算的口產(chǎn)氣量和累積產(chǎn)氣量，從圖3.3，3.4可以看出：生產(chǎn)初期考慮吸附和不考慮吸附的兩個(gè)模型計(jì)算的口產(chǎn)氣量和累積產(chǎn)氣量差別較??；生產(chǎn)中后期差別比較大，主要因?yàn)樯a(chǎn)初期游離氣對(duì)產(chǎn)能的影響比吸附氣大，隨著壓力的降低吸附氣對(duì)產(chǎn)能的影響變大。總體來(lái)看，在頁(yè)巖氣生產(chǎn)過(guò)程中，吸附對(duì)頁(yè)巖氣井累積產(chǎn)氣量有較大的影響，所以不能忽略吸附的影響。3.1.3等溫吸附規(guī)律對(duì)產(chǎn)能影響研究本文運(yùn)用蘭氏等溫吸附定律對(duì)吸附氣的解吸附過(guò)程進(jìn)行描述，為了研究相同吸附氣體積含量條件下

29、解吸附規(guī)律對(duì)頁(yè)巖氣藏產(chǎn)能的影響，在相同開(kāi)發(fā)條件的前提下，共設(shè)計(jì)5套模擬計(jì)算方案，方案數(shù)據(jù)見(jiàn)表3.2。設(shè)每套方案的初始吸附氣體積含量?jī)r(jià)相同，且吸附氣儲(chǔ)量占天然氣總地質(zhì)儲(chǔ)量的50%，對(duì)不同方案的產(chǎn)能進(jìn)行對(duì)比分析，見(jiàn)圖3.53.6，其中圖3.5為不同蘭氏體積VL和蘭氏壓力PL對(duì)應(yīng)的口產(chǎn)氣量，圖3.6為不同蘭氏體積VL和蘭氏壓力PL對(duì)應(yīng)的累積產(chǎn)氣量。表3.2模擬方案參數(shù)方案1方案2方案3方案4方案5VL(m3/m3)2015131313PL(MPa)333.532.5表3.3不同方案對(duì)應(yīng)氣井初始產(chǎn)氣量參數(shù)方案1方案2方案3方案4方案5初始產(chǎn)氣量(m3/d)3.4085.1985.5456.6808.7

30、13從圖3.5我們可以看出：在相同吸附氣體積含量的前提下，蘭氏體積VL和蘭氏壓力PL對(duì)氣井的初始產(chǎn)量和產(chǎn)量遞減速度有直接影響，5種方案所對(duì)應(yīng)的初始產(chǎn)氣量見(jiàn)表3.6，從中可以看出：相同吸附氣體積含量時(shí)，隨著蘭氏體積VL增大，氣井初始產(chǎn)氣量減少；隨著蘭氏壓力PL增大，氣井初始產(chǎn)氣量減少。這主要是由于蘭氏體積越大表示此地層吸附頁(yè)巖氣能力越強(qiáng)，越不易發(fā)生解吸，則開(kāi)采的初始產(chǎn)氣量越低；而蘭氏壓力越大表示地層解吸對(duì)壓力要求越高，越不易發(fā)生解吸，則開(kāi)采的初始產(chǎn)量越低。從圖3.5，中可以看出：隨蘭氏體積VL的增大，頁(yè)巖氣井產(chǎn)量遞減速度加快；隨蘭氏壓力PL的增大，頁(yè)巖氣井產(chǎn)量遞減速度加快。這主要是由于蘭氏體積V

31、L越大表示此地層吸附頁(yè)巖氣能力越強(qiáng)，越不易發(fā)生解吸，前期開(kāi)采主要依靠游離氣，隨著游離氣的減少后期產(chǎn)量沒(méi)有解吸氣的維持，急劇下降，則產(chǎn)量的遞減速度加快。圖3.5不同蘭氏體積和蘭氏壓力對(duì)應(yīng)的日產(chǎn)氣量從圖3.6我們可以看出：在相同吸附氣體積含量的前提下，蘭氏體積VL和蘭氏壓力PL對(duì)氣井的采氣速度及相同開(kāi)發(fā)時(shí)間累積產(chǎn)氣量有直接影響。5種方案所對(duì)應(yīng)的相同開(kāi)發(fā)時(shí)間累積產(chǎn)氣量見(jiàn)表3.4，從中可以看出：在相同的開(kāi)發(fā)時(shí)間前提下進(jìn)行對(duì)比，隨蘭氏體積VL的增加，頁(yè)巖氣井累積產(chǎn)氣量降低；隨蘭氏壓力PL的增加，頁(yè)巖氣井累積產(chǎn)氣量降低。這主要是由于蘭氏體積及蘭氏壓力越大說(shuō)明地層的吸附能力越強(qiáng)且越不易發(fā)生解吸，解吸氣在累積

32、產(chǎn)量中的貢獻(xiàn)越小，累積產(chǎn)量則越小。表3.4不同方案所對(duì)應(yīng)氣的累積產(chǎn)氣量(108m3)參數(shù)方案1方案2方案3方案4方案510.0490.1020.1220.1600.20730.120.2460.2760.3550.48550.1710.3410.3730.4840.66880.2230.4360.4850.6150.826100.2290.4490.510.6490.861采氣速度指一個(gè)氣田(或開(kāi)發(fā)區(qū))年采氣量與地質(zhì)儲(chǔ)量(或可采儲(chǔ)量)的比值，是衡量氣井生產(chǎn)速度快慢的指標(biāo)，也是氣田開(kāi)發(fā)中極為重要的開(kāi)發(fā)指標(biāo)。5種方案所對(duì)應(yīng)的采氣速度見(jiàn)表3.4，從中可以看出：隨蘭氏體積臉越大，頁(yè)巖氣井采氣速度越??；

33、隨蘭氏壓力PL越大，百巖氣采氣速度VL越小。表3.5不同方案對(duì)應(yīng)采氣速度方案方案1方案2方案3方案4方案5采氣速度(%)0.0089360.0187660.0224530.0293830.038141圖3.6不同蘭氏體積和蘭氏壓力對(duì)應(yīng)的累積產(chǎn)氣量3.1.4吸附氣含量對(duì)產(chǎn)能影響研究頁(yè)巖氣通常以游離態(tài)和吸附態(tài)賦存于儲(chǔ)層中，其中吸附氣的存在不僅增加了氣體的賦存方式，也對(duì)頁(yè)巖氣藏產(chǎn)能有所影響。為了研究吸附氣含量對(duì)頁(yè)巖氣藏產(chǎn)能的影響，分別模擬吸附氣含量為0， 10%， 30%， 50%， 70%， 90%和100%這七種方案，并對(duì)每個(gè)方案的頁(yè)巖氣井產(chǎn)量及地層壓力進(jìn)行對(duì)比分析，見(jiàn)圖3.63.8。表3.6不

34、同吸附氣含量對(duì)應(yīng)氣井初始產(chǎn)氣量吸附氣含量（%）01030507090100初始產(chǎn)氣量(104 m3/d)6.974836.810076.3876550.0293830.0381414.63.77圖3.7為不同吸附氣含量條件下口產(chǎn)氣量曲線圖，由圖3.8我們可以看出：吸附氣含量對(duì)氣井的初始產(chǎn)量及產(chǎn)量遞減速度有直接影響。7套方案所對(duì)應(yīng)的初始產(chǎn)氣量見(jiàn)表3.5，從中可以看出：氣井初始產(chǎn)氣量隨吸附氣含量增高而減少。這主要由于開(kāi)采初期產(chǎn)量主要來(lái)源于游離氣，吸附氣含量越高，則游離氣含量越低，初始產(chǎn)量初始產(chǎn)量則越低。從圖3.7中可以看出：氣井產(chǎn)量遞減速度隨吸附氣含量增高而減緩，這主要是由于吸附氣解吸受地層壓力影

35、響，吸附氣含量越高，開(kāi)采過(guò)程中地層壓力變化越小，解吸速度的越緩慢，產(chǎn)量變化越小，這也正是頁(yè)巖氣藏能長(zhǎng)期穩(wěn)定生產(chǎn)的原因。圖3.7不同吸附氣含量條件下日產(chǎn)氣量曲線圖表3.7不同吸附氣含量相同開(kāi)發(fā)時(shí)問(wèn)所對(duì)應(yīng)氣的累積產(chǎn)氣量(10g/m3)吸附器含量（%）0103050709010010.1650.1570.1430.1300.1140.0940.08430.3500.3380.3100.2760.2390.2010.18150.4690.4490.4090.3670.3210.2720.24780.5730.5560.5100.4530.3970.3370.304100.5840.5690.5310.

36、4790.4130.3420.305圖3.8為不同吸附氣含量條件下累積產(chǎn)氣量曲線圖，由圖3.6我們可以看出：吸附氣含量對(duì)采氣速度和相同開(kāi)發(fā)時(shí)間對(duì)應(yīng)的累積產(chǎn)氣量有直接影響。表3.8不同吸附氣含量對(duì)應(yīng)采氣速度吸附氣含量（%）01030507090100采氣速度(%)0.0303870.0289290.0263530.0239450.0209650.0173040.015412圖3.8吸附氣含量不同的條件下累積產(chǎn)氣量曲線圖7種方案分別在365天、500天、1000天、2000天、3000天所對(duì)應(yīng)的累積產(chǎn)氣量見(jiàn)表3.6，從中可以看出：氣井累積產(chǎn)氣量隨吸附氣含量增高而減少。這主要是由于吸附氣含量越高，氣

37、體越不易被采出，累積產(chǎn)氣量越低。7種方案所對(duì)應(yīng)的采氣速度見(jiàn)表3.7，從中可以看出：采氣速度隨吸附氣含量增高而逐漸降低。這主要是由于吸附氣需解吸后才能開(kāi)采，吸附氣越高，氣體越不易被采出，采氣速度越低。圖3.9是吸附氣含量不同的條件下地層平均壓力曲線，從圖中我們可以看出：頁(yè)巖氣地質(zhì)儲(chǔ)量一定時(shí)，隨吸附氣含量增高，地層平均壓力下降速度緩慢。由于吸附氣在解吸過(guò)程中有壓力補(bǔ)給，壓力下降緩慢。圖3.9不同吸附氣含量條件下地層平均壓力曲線圖3.2儲(chǔ)層滲透率對(duì)產(chǎn)能影響研究3.2.1模型描述及基本參數(shù)表3.9模擬計(jì)算的有關(guān)參數(shù)計(jì)算參數(shù)參數(shù)值計(jì)算參數(shù)參數(shù)值原始地層溫度（K）361基質(zhì)滲透率（10-3m2）0.001

38、氣體摩爾質(zhì)（kg/mol）0.016孔隙度（%）5氣體常數(shù)（Pa·m3/(mol·k)）8.314標(biāo)況下氣體密度（kg/m3）0.717內(nèi)邊界壓力（MPa）4原始地層壓力（MPa）20蘭氏體積（m3/m3）13氣體粘度（mPa.s）0.067蘭氏壓力（MPa）3壓縮系數(shù)（MPa-1）1×10-4總地質(zhì)儲(chǔ)量（m3）5.43×108井筒半徑（m）0.1基質(zhì)滲透率（10-3m2）5通過(guò)借鑒國(guó)外頁(yè)巖儲(chǔ)層典型參數(shù)建立了頁(yè)巖氣藏三維模型，設(shè)頁(yè)巖氣藏地層模型為長(zhǎng)1210m，寬810m，有效厚度30m，模擬時(shí)間為10a，壓裂采用體積壓裂方法，縫網(wǎng)尺寸為600m×

39、;200m，模型計(jì)算參數(shù)見(jiàn)表3.8。采用混合單元網(wǎng)格進(jìn)行剖分，模擬時(shí)將裂縫視為一個(gè)面，用三角形面單元進(jìn)行剖分，地層用任意四面體單元剖分。模型網(wǎng)格剖分圖如圖3.10所示，網(wǎng)格加密區(qū)為采用體積壓裂區(qū)域。(a)頂部平面網(wǎng)格圖 (b)立體網(wǎng)格剖分圖圖3.10網(wǎng)格剖分圖圖3.11分別為模擬模型0.5年、2.5年、5年、8年的壓力分布圖，可以看出生產(chǎn)過(guò)程中各時(shí)間的壓力分布情況。(a)0.5年時(shí)的壓力分布圖 (b)2.5年時(shí)的壓力分布圖(c)5年時(shí)的壓力分布圖 (d)8年時(shí)的壓力分布圖圖3.11壓力分布圖3.2.2基巖滲透率對(duì)產(chǎn)能影響研究為比較不同基巖滲透率對(duì)頁(yè)巖氣藏產(chǎn)能的影響，共設(shè)計(jì)了5套方案，如表3.1

40、0所示。表3.10模擬方案參數(shù)方案1方案2方案3方案4方案5基巖滲透率(10-3m2)0.010.0050.0010.00050.0001圖3.12圖3.13為分別不同基巖滲透率下頁(yè)巖氣口產(chǎn)氣量和累積產(chǎn)氣量變化曲線。從圖中我們可以看出：基巖滲透率越大，頁(yè)巖氣口產(chǎn)量越大，累積產(chǎn)氣量越大；基巖滲透率越小，基巖滲透率對(duì)頁(yè)巖氣藏產(chǎn)能的影響越小，主要原因是：基巖滲透率越小導(dǎo)致基巖滲透率與裂縫滲透率的差距越大，這時(shí)裂縫成為主要的滲流通道，所以基巖對(duì)生產(chǎn)速度的影響越小。圖3.12不同基巖滲透率下頁(yè)巖氣日產(chǎn)氣量變化曲線圖3.13不同基巖滲透率下頁(yè)巖氣累積產(chǎn)氣量變化曲線3.2.3裂縫滲透率對(duì)產(chǎn)能影響研究為比較不

41、同裂縫滲透率對(duì)頁(yè)巖氣藏產(chǎn)能的影響，共設(shè)計(jì)了6套方案，如表3.11所示。表3.11模擬方案參數(shù)方案1方案2方案3方案4方案5方案6裂縫滲透率(10-3m2)2010510.50.1圖3.14圖3.15為不同裂縫滲透率下頁(yè)巖氣井產(chǎn)能變化的曲線。由圖3.14可知，裂縫滲透率越大，口產(chǎn)氣量越大，穩(wěn)產(chǎn)期越短，口產(chǎn)量下降的越早；裂縫滲透率越小，生產(chǎn)初期口產(chǎn)量較小，穩(wěn)產(chǎn)期越長(zhǎng)；從圖3.15我們可以看出：裂縫滲透率越大，相同生產(chǎn)時(shí)間時(shí)累積產(chǎn)氣量越大，但是裂縫滲透率小于10m2，不同裂縫滲透率累積產(chǎn)量的差別較大，當(dāng)裂縫滲透率大于10 m2時(shí)，裂縫滲透率對(duì)累積產(chǎn)量的影響變小，尤其生產(chǎn)后期累積產(chǎn)量己相差不大。圖3.14不同裂縫滲透率下頁(yè)巖氣日產(chǎn)氣量變化曲線圖3.15不同裂縫滲透率下頁(yè)巖氣累積產(chǎn)氣量變化曲線對(duì)比基巖與裂縫滲透率對(duì)頁(yè)巖氣產(chǎn)能的影響，裂縫滲透率對(duì)產(chǎn)能影響較大，基巖滲透率對(duì)產(chǎn)能影