董 浩，馬 輪，韋樂樂，王麗娜，王 芳，宋德彬，趙思琦，王 云

(1.西北大學(xué)大陸動力學(xué)國家重點實驗室/西北大學(xué)地質(zhì)學(xué)系，陜西西安 710069;2.中國石油勘探開發(fā)研究院西北分院/中國石油天然氣集團(tuán)公司油藏描述重點實驗室，甘肅蘭州 730020;3.中國石油天然氣管道局第四工程分公司，河北廊坊 065000;4.中國石油集團(tuán)長城鉆探工程有限公司測井公司遼河項目部，遼寧盤錦 124010;5.中石油管道公司西安輸油氣分公司，陜西西安 7100185)

由于頁巖氣儲層滲流條件非常差，目前對于頁巖氣藏的開發(fā)還局限于單井。頁巖氣藏單井開發(fā)效果的影響因素［1～8］很多，既有基質(zhì)滲透率、天然裂縫發(fā)育程度等決定頁巖氣藏開發(fā)效果的客觀因素，又有壓裂改造程度、生產(chǎn)制度等決定頁巖氣藏開發(fā)效果的主觀因素。另外，頁巖氣解吸作用和擴(kuò)散作用等也會對產(chǎn)能產(chǎn)生不可忽略的影響［9～15］。本文通過對以上諸多因素進(jìn)行全面分析，為頁巖氣藏體積壓裂的產(chǎn)能影響因素進(jìn)行權(quán)重分析，進(jìn)而為開發(fā)不同條件的頁巖氣藏提供理論指導(dǎo)。為了方便對比，本文采用一套模型參數(shù)(表1)。

表1 模型基礎(chǔ)參數(shù)表Table 1 basic parameters of model

1 基質(zhì)滲透率

頁巖氣藏基質(zhì)是氣體的主要儲存空間，氣體首先通過基質(zhì)運移到天然或人工裂縫中，然后由裂縫流入井筒，因此頁巖氣藏基質(zhì)的滲透率大小對產(chǎn)能的影響很大［16］。考慮到模型計算能力不能過大，假設(shè)水平段長為600m，壓裂人工裂縫6段，每段分為3簇，簇間距為25m，人工裂縫滲透率為200mD，縫長為150m，天然裂縫與人工裂縫正交分布，長度由內(nèi)向外逐漸減小，縫間距為25m，通過局部網(wǎng)格加密實現(xiàn)，水力壓裂后裂縫示意圖如圖1所示。

圖1 人工裂縫和天然裂縫分布示意圖Fig.1 Distribution of Artificial fractures and natural fractures

通過調(diào)整頁巖基質(zhì)滲透率來模擬其對單井產(chǎn)量的影響，基質(zhì)滲透率Km由0.0001mD增加到0.01mD，模擬得到日產(chǎn)氣量和10年累計產(chǎn)氣量(圖2)。從圖2可以看出，基質(zhì)滲透率越大，初始日產(chǎn)量越高，日產(chǎn)氣量下降越快。當(dāng)基質(zhì)滲透率Km小于0.001mD時，基質(zhì)滲透率的增加對日產(chǎn)氣量的影響不大，基質(zhì)滲透率Km從0.0001mD增加到0.0005mD時日產(chǎn)氣量增加不大，生產(chǎn)1000天后兩者的日產(chǎn)氣量基本沒有差別;而當(dāng)基質(zhì)滲透率Km大干0.001mD時，基質(zhì)滲透率的增加對日產(chǎn)氣量的增加影響較大。由圖2可見，當(dāng)基質(zhì)滲透率Km從0.001mD增加到0.005mD時，日產(chǎn)氣量增加幅度較大，特別是初期日產(chǎn)氣量的增幅更大。出現(xiàn)這種現(xiàn)象的原因是，當(dāng)基質(zhì)滲透率過小時 (小于0.001mD)，氣體在基質(zhì)中的流動很困難，基質(zhì)滲透率即使成倍的增加對氣體的流動改善也非常有限;而當(dāng)基質(zhì)滲透率相對較大時 (大于0.001mD)，稍微增加基質(zhì)滲透率就能較大幅度地提高日產(chǎn)氣量?？梢?，當(dāng)基質(zhì)滲透率大干0.001mD時，滲透率的增加對于頁巖氣產(chǎn)能的影響要大。

圖2 不同基質(zhì)滲透率條件下產(chǎn)氣量對比圖Fig.2 Comparison of gas production rate for different matrix permebilities

圖3所示為不同基質(zhì)滲透率頁巖氣藏10年末儲層壓力圖，從圖中可以看出，當(dāng)基質(zhì)滲透率Km為0.0001mD時，10年末的壓力降還圍繞在人工裂縫和天然裂縫之間，說明滲透率太小氣體流動困難，壓力外擴(kuò)的速度比較慢。隨著基質(zhì)滲透率的增加，10年末的壓力降范圍逐漸擴(kuò)大。低滲透率時人工裂縫間的氣體難以波及到，氣體無法動用，然而當(dāng)基質(zhì)滲透率為0.005mD時，人工裂縫間的壓力降逐漸連成一片，說明人工裂縫間的區(qū)域也得到了動用，期間的頁巖氣隨著壓力的下降而流入裂縫，進(jìn)入井筒被采出。當(dāng)基質(zhì)滲透率Km為0.01mD時，裂縫間的壓力降與裂縫內(nèi)的壓力降基本一致，說明裂縫間的頁巖氣得到了很好的動用，這也是不同基質(zhì)滲透率對產(chǎn)能影響的根本原因。

圖3 不同基質(zhì)滲透率條件下的10年末儲層壓力圖Fig.3 Formation pressure distribution after ten－year production for different matrix permeabilities

2 人工裂縫參數(shù)

目前頁巖氣藏的開發(fā)主要通過大規(guī)模水力壓裂方法，因此人工裂縫對于頁巖氣藏產(chǎn)能的影響很大。人工裂縫參數(shù)主要包括裂縫長度Lf(即裂縫半長)、裂縫導(dǎo)流能力Df及裂縫條數(shù)Nf。模擬過程中，假設(shè)水平段長1000m，天然裂縫仍然平行于水平井筒分布，具體參數(shù)如表1所示。

2.1 裂縫長度

通過調(diào)整裂縫半長Lf來模擬裂縫長度 (簡稱縫長)對頁巖氣藏單井日產(chǎn)氣量和累計產(chǎn)氣量的影響 (圖4)。圖4a所示為裂縫半長Lf從75m增加到200m時日產(chǎn)氣量的變化，從圖中可以看出，隨著裂縫半長的增加，日產(chǎn)氣量增加;且裂縫半長越大，初期日產(chǎn)氣量越大 (當(dāng)裂縫半長Lf為75m時，初期日產(chǎn)氣量為23×104m3，當(dāng)裂縫半長Lf為200m時，初期日產(chǎn)氣量為37×104m3。同時，隨著裂縫半長的增加，日產(chǎn)氣量的遞減率增大，到10年末時，裂縫半長為200m的日產(chǎn)氣量只比裂縫半長為75m的日產(chǎn)氣量多1000m3左右。從圖4b中可以看出，隨著裂縫半長的增加，10年累計產(chǎn)氣量增加，且產(chǎn)量的增加與裂縫半長的增加基本對應(yīng)，即裂縫半長每增加25m，累計產(chǎn)氣量的增幅基本相同。

圖5所示為不同裂縫半長下10年末的壓力分布圖，圖中僅顯示了裂縫周圍相同區(qū)域內(nèi)的壓力分布。從圖5中可以看出，裂縫半長越小，10年末的壓力波及區(qū)域越小，這充分說明丁人工裂縫長度對產(chǎn)量的影響程度。眾所周知，低滲透油氣藏中，人工裂縫的作用相當(dāng)于交通運輸中的高速公路，由于油氣藏的低滲透性，應(yīng)盡可能地修長高速公路，使更多的油氣流入其中，這樣才更有利于疏導(dǎo)交通。若裂縫半長較小，裂縫控制范圍就小，其外的油氣藏很難流進(jìn)裂縫這條高速公路里，也就很難被采出。顯然，對于頁巖氣藏的超低滲透率的儲層來說，裂縫半長的意義更加重要。

圖4 不同裂縫半長條件下的產(chǎn)氣量對比圖Fig.4 Comparison of gas production rate for different fracture lengths

2.2 裂縫導(dǎo)流能力和裂縫條數(shù)

裂縫導(dǎo)流能力為裂縫滲透率與裂縫寬度的乘積。裂縫導(dǎo)流能力的大小表征了裂縫運輸油氣的能力，裂縫導(dǎo)流能力越大，輸氣能力越強;反之，輸氣能力越弱［17］。通過調(diào)整人工裂縫的裂縫導(dǎo)流能力來模擬其對頁巖氣藏單井產(chǎn)能的影響，當(dāng)裂縫導(dǎo)流能力從10D·cm增加到40D·cm時，累計產(chǎn)氣量幾乎沒有增加，這與低滲透油氣藏中裂縫導(dǎo)流能力越大產(chǎn)量越高不吻合，這一現(xiàn)象主要是由頁巖氣藏的滲透率極低造成的。當(dāng)裂縫導(dǎo)流能力為10D·cm時，人工裂縫的輸氣能力已遠(yuǎn)超過頁巖基質(zhì)的輸氣能力，裂縫導(dǎo)流能力再增大到40D·cm時產(chǎn)量也不會有明顯增加。這一現(xiàn)象也可用高速公路來比喻，當(dāng)車流量很小時，雙車道和4車道的高速公路的運輸能力可能是相當(dāng)?shù)模黾榆嚨罃?shù)不會對運輸能力有大的提高。

圖5 不同裂縫半長條件下的10年末儲層壓力分布圖Fig.5 Formation pressure distribution for different fracture lengths

由于頁巖氣藏儲層的致密性和人工裂縫的局限性，天然裂縫的復(fù)雜程度對頁巖氣單井產(chǎn)能具有決定性的意義。若天然裂縫非常發(fā)育，則很少的人工裂縫就能產(chǎn)生很好的溝通作用，形成很大的改造體積;若天然裂縫極其不發(fā)育，則在目前的工藝條件下，再多的人工裂縫也無法對儲層產(chǎn)生較好的改造作用?；趯θ斯ち芽p參數(shù)的優(yōu)化，選取人工裂縫條數(shù)為6條、裂縫導(dǎo)流能力為10D·cm、裂縫半長為200m，通過調(diào)整天然裂縫參數(shù)分析其對頁巖氣單井產(chǎn)能的影響。將裂縫條數(shù)從2條增加到12條來模擬其對頁巖氣單井產(chǎn)能的影響。為了模擬不同的裂縫條數(shù)，需要對模擬區(qū)域進(jìn)行不同形式的網(wǎng)格加密，工作量很大，為了計算簡便且不影響結(jié)果，以加密12條縫的情況對模擬區(qū)域進(jìn)行加密，而后根據(jù)不同的需要對表示裂縫的網(wǎng)格進(jìn)行不同的設(shè)置來模擬人工裂縫，當(dāng)裂縫條數(shù)達(dá)到10條或12條時，10年末的壓力已經(jīng)覆蓋了整個加密區(qū)域，對加密區(qū)域外僅有少部分波及，在這種情況下可以將加密區(qū)域視為一個小氣藏，外部的氣體很難流入裂縫中被采出來。值得注意的是，當(dāng)裂縫條數(shù)從2條增加到6條時的累計產(chǎn)氣量增加非常巨大，而當(dāng)裂縫條數(shù)大于6條之后，累計產(chǎn)氣量隨裂縫條數(shù)的增加的幅度不明顯。這是因為當(dāng)裂縫條數(shù)大于6條后，多余的裂縫條數(shù)對于改造體積的增加意義不大，因此10年末的累計產(chǎn)氣量增加不大。由此可知，對于本例，6條裂縫的效果最佳。

3 天然裂縫參數(shù)

3.1 裂縫導(dǎo)流能力

按照之前的高速公路理論，天然裂縫導(dǎo)流能力相當(dāng)于鏈接高速公路的輔路的運輸能力，導(dǎo)流能力大則允許進(jìn)入人工裂縫的油氣流大，在油氣量足夠大的情況下，裂縫導(dǎo)流能力越大，高速公路的運輸能力就越強。因此，調(diào)整天然裂縫導(dǎo)流能力從0.1D·cm至5.0D·cm，分析裂縫導(dǎo)流能力的增加對累計產(chǎn)氣量的影響。隨著天然裂縫導(dǎo)流能力的增加，10年累計產(chǎn)氣量逐漸增加。相比于天然裂縫導(dǎo)流能力從1.0D·cm增加到5.0D·cm時的累計產(chǎn)氣量增幅，天然裂縫導(dǎo)流能力從0.1D·cm增加到0.5D·cm時，累計產(chǎn)氣量的增幅更大，這主要是因為裂縫導(dǎo)流能力從0.1D·cm到0.5D·cm之間的天然裂縫的滲流性已經(jīng)比基質(zhì)系統(tǒng)有很大的提升，多余的裂縫導(dǎo)流能力增量也起不到大幅增加產(chǎn)氣量的作用［18，19］。

3.2 天然裂縫間距

天然裂縫間距，即兩條天然裂縫間的距離。假設(shè)天然裂縫與水平井筒平行，與人工裂縫垂直，因此天然裂縫的不同間距表示人工裂縫與儲層的連通性。天然裂縫間距小，說明其密度大，與人工裂縫的連通性好;天然裂縫間距大，說明其密度小，與人工裂縫的連通性差。由于前文假設(shè)人工裂縫半長為200m，所以對于間距大于200m的天然裂縫的模擬是沒有意義的。基于此，天然裂縫導(dǎo)流能力設(shè)為0.5D·cm，假設(shè)天然裂縫的間距分別為25m、50m、75m及100m，模擬不同間距條件下頁巖氣藏單井產(chǎn)能情況。為了考慮極端情況，增加間距為零的情形表征不存在天然裂縫的情況。

由表2可知:當(dāng)無天然裂縫時，單井初期日產(chǎn)氣量較低，存在天然裂縫的井初期日產(chǎn)氣量較高。隨著天然裂縫間距減小，初期日產(chǎn)量快速增加。天然裂縫間距由25m增加到50m時，初期日產(chǎn)氣量大幅度下降，而從50m增加到100m過程中，初期日產(chǎn)氣量下降的幅度減緩。這說明當(dāng)天然裂縫間距為25m時，能夠相對較好地對裂縫泄氣區(qū)域進(jìn)行覆蓋，而一旦天然裂縫密度降低，達(dá)不到覆蓋裂縫泄氣區(qū)域所需要的密度時，初期日產(chǎn)氣量就會大幅下降。因此，對于頁巖氣藏的開發(fā)，天然裂縫的發(fā)育程度越強越有利。

表2 不同天然裂縫間距時的初期日產(chǎn)氣量表Table 2 initial daily rate for different fracture interval

3.3 天然裂縫分布

實際地層中天然裂縫分布非常復(fù)雜，很難通過數(shù)值模擬方法準(zhǔn)確表征，更無法進(jìn)行因素分析。為了對天然裂縫的分布狀況進(jìn)行模擬，考慮了以下3種情況，即不存在天然裂縫的情況、天然裂縫不完全連通的情況及天然裂縫完全連通的情況。不存在天然裂縫的情況可以理解為天然裂縫完全不連通，頁巖氣的開發(fā)依靠基質(zhì)和人工裂縫滲流;天然裂縫不完全連通，是指存在相互不溝通的天然裂縫;天然裂縫完全連通，是指天然裂縫完全貫通整個人工裂縫系統(tǒng)，人工裂縫與天然裂縫形成裂縫網(wǎng)絡(luò)［20～23］。

對于天然裂縫完全連通的情形，改造區(qū)域內(nèi)的頁巖氣初期動用比較充分;而對于天然裂縫不完全連通或無天然裂縫的情形，改造區(qū)域內(nèi)的頁巖氣在初期無法被充分動用，因此初期日產(chǎn)氣量比較小。但是，隨著改造區(qū)域內(nèi)氣體的采出，完全連通情況下的日產(chǎn)氣量的下降速度比另外兩種情況快得多。

4 敏感因素分析及改造體積

由前述分析可以看出，不同參數(shù)對產(chǎn)能的影響最終都?xì)w結(jié)為改造體積的大小對產(chǎn)能的影響，這一點與低滲透油藏非常相類似。在低滲透油藏的開發(fā)中，單井泄油面積對產(chǎn)能有直接影響，而泄油面積的大小又與地層滲透率、人工裂縫參數(shù)等因素有關(guān)［24，25］。因此，對頁巖氣藏產(chǎn)能的影響因素進(jìn)行深入研究非常必要，即找到影響改造體積的幾個關(guān)鍵因素，這對頁巖氣藏初期開發(fā)方案制訂有著非常重要的意義。本文首先采用正交試驗方法對以上各因素進(jìn)行綜合評價，找出3個對產(chǎn)能影響最大的因素;然后根據(jù)其影響程度的大小設(shè)計方法繪制圖版;最后通過一種BP人工神經(jīng)網(wǎng)絡(luò)的方法確定任意條件下的改造體積，為頁巖氣藏快速產(chǎn)能預(yù)測提供有效方法。

4.1 正交試驗方法

單因素或雙因素試驗，實施起來較為簡單。但是當(dāng)考察的因素在3個或3個以上時，進(jìn)行全方位的試驗規(guī)模將異常龐大。通常情況下，當(dāng)試驗因素過多時，試驗方案太多而導(dǎo)致試驗不可行。正交方法設(shè)計可以解決這一難題，它是一種高效的針對多因素、多水平的優(yōu)化方法。它通過從總體實驗方法中挑選出具有代表性的試驗，對它們進(jìn)行評估來找出最優(yōu)的組合，即可以從一部分試驗了解到全部試驗。當(dāng)試驗因素比較多時，正交法可以大幅度降低試驗的次數(shù)，從而使開展試驗成為可能。

正交試驗包括兩個要素:因素和水平。因素是指試驗的項目，如前述的基質(zhì)滲透率、人工裂縫長度等;水平是指每個因素的取值個數(shù)，如基質(zhì)滲透率從0.0001mD增加到0.0005mD為5個水平。正交試驗具有水平均勻性和搭配均勻性的特點。水平均勻性是指所選的方法對每個因素及其對應(yīng)的水平都是均勻分配的;搭配均勻性是指每個因素的每個水平在所有方案中出現(xiàn)的次數(shù)相同，且任意兩個因素的搭配以相同的次數(shù)出現(xiàn)。

由于需要對比的因素依然非常多，首先需要對各單因素分析定性篩選出比較敏感的因素重點研究。如前所述，人工裂縫長度對產(chǎn)能的影響比較大，為非常敏感的參數(shù);裂縫數(shù)量在大于6條后，其對產(chǎn)能的影響非常小，可以忽略不計;而裂縫導(dǎo)流能力相對于基質(zhì)滲透率非常大，因此它的敏感性很弱。對于天然裂縫參數(shù)，當(dāng)天然裂縫的導(dǎo)流能力大于0.5D·cm后，產(chǎn)能對導(dǎo)流能力不敏感;天然裂縫的間距對產(chǎn)能的影響很大;由于頁巖氣的開發(fā)效果受天然裂縫的連通性影響太大，為了避免天然裂縫的分布對結(jié)果的影響太大造成偏差，這里不把天然裂縫分布作為正交試驗的因素。另外，頁巖氣開發(fā)中氣體解吸對后期產(chǎn)能的影響不容忽略，基質(zhì)滲透率同樣對產(chǎn)能影響很大，因此認(rèn)為基質(zhì)滲透率、人工裂縫長度、天然裂縫間距及氣體吸附比為正交試驗的4個合理因素。

若要得到各因素對指標(biāo)影響的主次關(guān)系，需要對正交試驗各個方案的結(jié)果進(jìn)行極差分析。極差反應(yīng)了各因素的水平變化時，指標(biāo)的變化幅度，即當(dāng)基質(zhì)滲透率、人工裂縫長度、天然裂縫間距及氣體吸附比例變化時，頁巖氣藏單井10年累計產(chǎn)氣量的變化幅度。

4.2 影響因素的排序

為了模擬基質(zhì)滲透率、人工裂縫長度、天然裂縫間距和氣體吸附比等4個因素對產(chǎn)能影響的主次關(guān)系，假設(shè)水平段長1000m，人工壓裂6段，每段3簇，人工裂縫導(dǎo)流能力為10D·cm，天然裂縫導(dǎo)流能力為0.5D·cm，天然裂縫不完全連通，每個因素分別取5個水平 (表3)。

表3 正交試驗中各因素的水平值表Table 3 Factor values in orthogonal test

通過正交設(shè)計，得到了25種方案，對模擬結(jié)果進(jìn)行極差分析。首先將各因素、各水平對應(yīng)的指標(biāo)數(shù)值相加，如對應(yīng)基質(zhì)滲透率的Kl，值即為0.0001mD的指標(biāo)值之和，對應(yīng)基質(zhì)滲透率的K3值即為0.005mD的指標(biāo)值之和，對應(yīng)人工裂縫長度的K3值即為縫長為250m的指標(biāo)值之和，kl至k5為Kl至K5的平均值。R表示極差，即Max(Ki) －Min(Ki)，i=l，2，3，4，5。

由表4可以看出，從極差R值大小來看，基質(zhì)滲透率最大，人工裂縫長度次之，天然裂縫間距再次，氣體吸附比最小。也就是說，以上4個參數(shù)值改變相同幅度時，10年累計產(chǎn)氣量對基質(zhì)滲透率的變化最為敏感。可以認(rèn)為，以10年累計產(chǎn)氣量為指標(biāo)時，影響因素的主次順序為基質(zhì)滲透率、人工裂縫長度、天然裂縫間距和氣體吸附比。

表4 極差分析表Table 4 Range analysis table

4.3 層次分析法

層次分析法是將與決策有關(guān)的元素分解成多個層次，在此基礎(chǔ)上進(jìn)行定性和定量分析。層次分析法建模基本上可以分為以下4個步驟。

4.3.1 建立遞階層次結(jié)構(gòu)模型

分為最高層、中間層和最底層。最高層，即是預(yù)定的指標(biāo)，本例中就是10年累計產(chǎn)氣量;中間層，包含了實現(xiàn)目標(biāo)所涉及的中間環(huán)節(jié)，即基質(zhì)滲透率、人工裂縫長度、天然裂縫間距和氣體吸附比;最底層，包括了實現(xiàn)目標(biāo)的具體措施和方案。

4.3.2 構(gòu)造各層次中的所有判斷矩陣

由于進(jìn)行比較的因素多，所以提出一種因子兩兩比較的方案，即每次取兩個因子xi和xj，以aij表示xi和xj對z的影響大小之比，結(jié)果用矩陣A=(aij)n×n表示，稱A為z—x之間的成對比較判斷矩陣?？梢钥闯觯魓i拜口xj對z的影響之比為aij，則xi墨和xj對z的影響之比應(yīng)為aij=1/aij。

關(guān)于如何確定aij的值，建議引用數(shù)字1～9及其倒數(shù)作為標(biāo)度。表5列出了1～9標(biāo)度的含義。

表5 9標(biāo)度法表Table 5 Nine scale analysis method

4.3.3 層次單排序及一致性檢驗

通過以上成對比較之后，雖然能一定程度上減少其他因素的干擾，但是仍然會客觀存在一定的非一致性，因此需要進(jìn)行一致性檢驗，檢驗過程可查閱相關(guān)參考文獻(xiàn)［24～26］。

4.3.4 層次總排序及一致性檢驗

以上得到了一組元素對上一層元素的權(quán)重向量，最終還需要得到最底層各元素對目標(biāo)層的排序權(quán)重，從而進(jìn)行方案選擇，層次總排序就是自上而下對單排序的權(quán)重進(jìn)行合成［27～31］。

以10年累計產(chǎn)氣量為目標(biāo)層，以基質(zhì)滲透率、人工裂縫長度、天然裂縫間距和氣體吸附比為準(zhǔn)備層，判斷準(zhǔn)備層對目標(biāo)層的排序。實驗總結(jié)得出，基質(zhì)滲透率的權(quán)重占到43%，對產(chǎn)能的影響最大;人工裂縫長度次之，權(quán)重為25%;天然裂縫間距再次，權(quán)重為13%;氣體吸附比最小，權(quán)重為13%。

4.4 BP試驗法

PlackettBurman試驗 (BP試驗)，實際上就是一種篩選試驗，它是對于試驗次數(shù)多，且沒有明確了解各因素對指標(biāo)的影響強弱次序的情況下采用的一種試驗設(shè)計方法。此方法對每個因子取兩水平來進(jìn)行分析，通過比較各個因子兩水平的差異與整體的差異來確定因子的顯著性。圖6a“+”表示最大值，“－”表示最小值，字母A至G表示不同的因素。本例中，只用A至D范圍內(nèi)的表，E至G的表刪除即可。此方法為定性分析方法，所選指標(biāo)為10年累計產(chǎn)氣量。圖6b為4個參數(shù)對指標(biāo)的敏感程度，從圖中可以看出，基質(zhì)滲透率和人工裂縫長度對頁巖氣藏單井產(chǎn)能的影響較為顯著，而天然裂縫間距和氣體吸附比對產(chǎn)能的影響不是很顯著。

圖6 BP試驗設(shè)計表和各參數(shù)對指標(biāo)的敏感程度圖Fig.6 BP test design table and sensibility degree for different parameters

綜合以上分析可以看出，頁巖儲層的基質(zhì)滲透率對產(chǎn)能的影響最為顯著，然而對于指定的頁巖儲層，基質(zhì)滲透率變化幅度比較小，因此若要實現(xiàn)較好的開發(fā)效果，人工改造規(guī)模就顯得非常重要。另外，天然裂縫的發(fā)育程度也是影響產(chǎn)能的一個很重要的因素?？傊?，對于指定的頁巖儲層，在基質(zhì)滲透率和天然裂縫發(fā)育程度固定的情況下，應(yīng)盡可能地增加人工改造的力度，盡最大可能增加改造體積，實現(xiàn)對頁巖氣藏的較好開發(fā)。

為了能夠得到任意參數(shù)組合下的改造體積形狀和大小，本文通過建立BP人工神經(jīng)網(wǎng)絡(luò)模型對這一問題進(jìn)行了研究。BP(Back Propagation)網(wǎng)絡(luò)是1986年由Rumelhart和McCelland為首的科學(xué)家小組提出的，是一種誤差逆?zhèn)鞑ニ惴ㄓ?xùn)練的多層前饋網(wǎng)絡(luò)，它是對人腦或自然神經(jīng)網(wǎng)絡(luò)若干基本特性的抽象和模擬，是目前應(yīng)用最廣泛的神經(jīng)網(wǎng)絡(luò)模型之一。BP網(wǎng)絡(luò)模型由3個層組成:輸入層、輸出層和隱含層，其中隱含層可以是一層或多層，每層由多個神經(jīng)元組成 (圖7)。BP網(wǎng)絡(luò)模型多層前饋網(wǎng)絡(luò)的算法可簡述如下。

圖7 BP神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)圖Fig.7 BP artificial neural network model

(1)假設(shè):

①輸入層有m0個神經(jīng)元;

②網(wǎng)絡(luò)有n層，輸出層為第n層，第k層有mk個神經(jīng)元;

③uk(i)表示第k層第i個神經(jīng)元所接受的信息，wk(i，j)表示第k－1層第j個神經(jīng)元到第k層第i個神經(jīng)元的權(quán)重，ak(i)表示第k層第i個神經(jīng)元的輸出，ao(j)表示輸入的第j個分量;

④層與層之間的神經(jīng)元存在信息交換，但同一層之間的神經(jīng)元無信息交換;

⑤信息交換的方向是從輸入層到輸出層方向，沒有反向信息傳播。

(2)根據(jù)以上假設(shè)，網(wǎng)絡(luò)的輸入、輸出關(guān)系可表示為:

(3)對于具有多個隱含層的前饋神經(jīng)網(wǎng)絡(luò)，若激發(fā)函數(shù)為 S，即:f(x)=1/［1+exp(－x)］，且指標(biāo)函數(shù)取:

其中:

(4)權(quán)重迭代公式為:

然后選定學(xué)習(xí)的數(shù)據(jù)，隨機(jī)確定初始權(quán)矩陣，用學(xué)習(xí)數(shù)據(jù)計算網(wǎng)格輸出，再根據(jù)上述公式進(jìn)行反向修正，直到用完所有學(xué)習(xí)數(shù)據(jù)。BP網(wǎng)絡(luò)建模具有非線性映照能力、并行分布處理方式、自學(xué)習(xí)和自適應(yīng)能力、數(shù)據(jù)融合的能力和多變量系統(tǒng)等特點，是一種成功的模擬方法。

結(jié)合本文的目的，確定了采用BP網(wǎng)絡(luò)模型來解決改造面積形狀和大小為目標(biāo)函數(shù)的問題。用上述模擬實驗的結(jié)果為例，采用基質(zhì)滲透率、人工裂縫長度和天然裂縫間距與改造面積形狀和大小一一對應(yīng)的數(shù)據(jù)輸入模式，應(yīng)用MATLAB語言進(jìn)行了建模和計算。

在BP網(wǎng)絡(luò)模型對改造面積長寬比和改造面積長邊長的訓(xùn)練中對改造面積長寬比的訓(xùn)練次數(shù)達(dá)到261次，最終誤差為0.00095;對改造面積長邊長的訓(xùn)練次數(shù)為239次，最終誤差為0.00099;由此可見，BP網(wǎng)絡(luò)模型通過大量的訓(xùn)練對數(shù)值模擬結(jié)果實現(xiàn)很好的逼近 (圖8)。綜上所述，可以通過大量的數(shù)值模擬結(jié)果和BP網(wǎng)絡(luò)模型訓(xùn)練相結(jié)合的方法預(yù)測出不同基質(zhì)滲透率、人工裂縫長度和天然裂縫間距條件下改造面積的形狀和大小。

圖8 BP網(wǎng)絡(luò)模型對改造面積長寬比和長邊長的訓(xùn)練結(jié)果圖Fig.8 Exercise results for stimulated area length－wide ratio and area side length with BP network model

5 結(jié) 論

通過大量數(shù)值模擬對影響頁巖氣藏單井產(chǎn)能的因素進(jìn)行分析，包括基質(zhì)滲透率、人工裂縫參數(shù)、天然裂縫參數(shù)及氣體吸附比等，結(jié)果表明:

(1)隨著基質(zhì)滲透率、人工裂縫長度、人工裂縫條數(shù)、氣體吸附比增加及天然裂縫間距減小，單井產(chǎn)能增加，反之亦然。

(2)通過正交試驗方法得到了各因素對產(chǎn)能影響的主次關(guān)系，并采用層次分析法和BP試驗法對該關(guān)系進(jìn)行了驗證，得出對產(chǎn)能影響最大的3個因素為基質(zhì)滲透率、人工裂縫長度和天然裂縫間距。

(3)以上述3個因素為變量，通過大量的模擬計算繪制了圖版，建立了BP網(wǎng)絡(luò)模型對圖版數(shù)據(jù)進(jìn)行了訓(xùn)練，用以快速確定任意條件下的改造體積形狀和大小。

［1］鄒才能，董大忠，楊樺，等.中國頁巖氣形成條件及勘探實踐［J］.天然氣工業(yè)，2011，31(12):26－39．

［2］張林曄，李政，朱日房.頁巖氣的形成與開發(fā) ［J］．天然氣工業(yè)，2009，29(1):124－128.

［3］聶海寬，何發(fā)岐，包書景.中國頁巖氣地質(zhì)特殊性及其探對策 ［J］．天然氣工業(yè)，2011，31(11):111－116.

［4］孟慶峰，侯貴廷.阿巴拉契業(yè)盆地Marcellus頁巖氣藏地質(zhì)特征及啟示 ［J］.中國石油勘探，2012，17(1):67－73.

［5］劉成均，郭平，蒙春，等.頁巖氣開發(fā)利用及其前景分析 ［J］.石油科技論壇，2014，33(2):40－46．

［6］劉存忠.頁巖氣開發(fā)現(xiàn)狀及前景思考 ［J］．石油科技論壇，2013，32(3):32－37．

［7］張煥芝，何艷青.全球頁巖氣資源潛力及開發(fā)現(xiàn)狀［J］.石油科技論壇，2010，29(6):53－57．

［8］張抗，頁巖氣革命帶來油氣地質(zhì)學(xué)和勘探學(xué)的重大創(chuàng)新 ［J］．石油科技論壇，2012，31(6):37－41．

［9］楊金華，田洪亮，郭曉霞，等，美國頁巖氣水平井鉆井提速提效案例與啟示 ［J］．石油科技論壇，2013，32(6):44－48．

［10］陳桂華，祝彥賀，徐強.頁巖氣成藏的四性特征及對下?lián)P子地區(qū)頁巖氣勘探的啟示 ［J］．中國石油勘探，2012，17(5):63－70．

［11］張金川，金之鈞，袁明生.頁巖氣成藏機(jī)理和分布［J］．大然氣工業(yè)，2004，24(7):15－18．

［12］蔣裕強，董大忠，漆麟，等.頁巖氣儲層的藎本特征及其評價 ［J］．天然氣工業(yè)，2010，30(10):7－12．

［13］張金川，姜生玲，唐玄，等.我國頁巖氣富集類型及資源特點 ［J］．天然氣工業(yè)，2009，29(12):109－114．

［14］張金川，汪宗余，聶海寬，等.頁巖氣及其勘探研究意義 ［J］．現(xiàn)代地質(zhì)，2008，(4):640－646.

［15］李欣，段勝楷，孫揚，等.美國頁巖氣勘探開發(fā)最新進(jìn)展 ［J］．天然氣工業(yè)，2011，(8):124－126，142．

［16］于榮澤，張曉偉，卞亞南，等.頁巖氣藏流動機(jī)理與產(chǎn)能影響因素分析 ［J］．天然氣工業(yè)，2012，(9):10－15，126．

［17］程遠(yuǎn)方，董丙響，時賢，等，頁巖氣藏三孔雙滲模型的滲流機(jī)理 ［J］.天然氣工業(yè)，2012，32(9):44 －47，130．

［18］謝維揚，李曉平.水力壓裂縫導(dǎo)流的頁巖氣藏水平井穩(wěn)產(chǎn)能力研究［J］．天然氣地球科學(xué)，2012，23(2):387－392．

［19］王永輝，盧擁軍，李永平，等.非常規(guī)儲層壓裂改造技術(shù)進(jìn)展及應(yīng)用 ［J］．石油學(xué)報，2012，33(增刊):149－158．

［20］Durbin P A，Iaccarino G．An approach to local refinement of structured grids［J］．Journal of Computational Physics，2002，181:639－653．

［21］Mehl S，Hill M C，Leake S A．Comparison of local grid refinement methods for Modflow ［J］．Ground Water，2006，44(6):792－796．

［22］Kilic A，Ertekin T．Application or a local grid refinement protocol in highly faultcd rescrvoir architectures ［J］．Journal of Canadian petroleum technology，2003，42(4):58－69．

［23］Shaoul J R，Behr A，Mtshedlishvili G．Developing a tool for 3D reservoir simulation of hydraulically fractured wells．IPTC10182，2005．

［24］景龍華，基于AHP的企業(yè)現(xiàn)金流量分析 ［J］．現(xiàn)代商貿(mào)工業(yè)，2010，22(9):182－183.

［25］郭金玉，張忠彬，孫慶云.層次分析法的研究與應(yīng)用［J］.中國安全科學(xué)學(xué)報，2008，18(5):148－153．

［26］鄧雪，李家銘，曾浩健，等.層次分析法權(quán)重計算方法分析及其應(yīng)用研究 ［J］．?dāng)?shù)學(xué)的實踐與認(rèn)識，2012，42(7):93－100．

［27］郭金玉，張忠彬，孫慶云.層次分析法的研究與應(yīng)用［J］.中國安全科學(xué)學(xué)報，2008，18(5):148－153．

［28］鄧雪，李家銘，曾浩健，等.層次分析法權(quán)重計算方法分析及其應(yīng)用研究 ［J］．?dāng)?shù)學(xué)的實踐與認(rèn)識，2012，42(7):93－100．

［29］楊莉娜，楊斌，魯洪江，等.五百梯氣田長興組儲層裂縫發(fā)育程度的常規(guī)測井神經(jīng)網(wǎng)絡(luò)預(yù)測 ［J］．中國石油勘探，2011，16(1):63－69．

［30］甘云雁，張士誠，陳利，等，復(fù)雜斷塊油藏不規(guī)則井網(wǎng)整體壓裂優(yōu)化設(shè)訃 ［J］.石油學(xué)報，2006，27(4):81－84.

［31］Rumelhart D E，McCelland J L．Parallel distribution processing(Vols I and II)［M］.Cambridge:MIT Press，1986．