劉海松 匡建超， 霍志磊 王 眾

（1.成都理工大學(xué)能源學(xué)院，四川 成都 610059；2.成都理工大學(xué)管理科學(xué)學(xué)院，四川 成都 610059）

0 引言

以往對(duì)儲(chǔ)層層內(nèi)非均質(zhì)性的研究，大多采用將各層的儲(chǔ)層參數(shù)簡(jiǎn)單相加進(jìn)行比較分析的方法，此法過(guò)于簡(jiǎn)單、系統(tǒng)性不強(qiáng)且對(duì)事物刻畫(huà)不準(zhǔn)確［1-2］。近年來(lái)，許多學(xué)者在儲(chǔ)層非均質(zhì)性定量研究方面做了大量工作，包括：常規(guī)反演和神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)相結(jié)合技術(shù)［3］、非均質(zhì)綜合指數(shù)法和模糊數(shù)學(xué)模型［4］、主因子分析法［5］、局部變差函數(shù)法［6］、儲(chǔ)層相控隨機(jī)建模［7］等。然而這些評(píng)價(jià)方法又都過(guò)于復(fù)雜，實(shí)用性不強(qiáng)。針對(duì)這一問(wèn)題，筆者構(gòu)建了一種較為簡(jiǎn)單、實(shí)用的評(píng)價(jià)方法——灰色集對(duì)分析評(píng)價(jià)法。

1 灰色集對(duì)分析評(píng)價(jià)方法

集對(duì)分析的主要思路是對(duì)不確定性系統(tǒng)中的集對(duì)作同異反分析，然后計(jì)算集對(duì)的同異反聯(lián)系度。集對(duì)是集對(duì)分析的基礎(chǔ)元素，而聯(lián)系度是集對(duì)分析的關(guān)鍵環(huán)節(jié)［8］。

設(shè)有聯(lián)系的集合A和B。A有n項(xiàng)表征其特性，即A=（a1，a2，…，an），B也有n項(xiàng)表征其特性，B=（b1，b2，…，bn）。A和B構(gòu)成集對(duì)H（A，B）。若要討論集對(duì)H（A，B）的關(guān)系，就需要計(jì)算H（A，B）的聯(lián)系度，其表達(dá)式為：

式中，s為同一性的個(gè)數(shù)；f為差異性的個(gè)數(shù)；i為差異不確定性系數(shù)；p為對(duì)立性的個(gè)數(shù)，在區(qū)間（-1，1）視不同情況取值；j為對(duì)立系數(shù)，且j≡-1。

1.1 評(píng)價(jià)方法

基于系統(tǒng)分析的觀點(diǎn)，選擇儲(chǔ)層定量評(píng)價(jià)指標(biāo)體系，并確定相應(yīng)等級(jí)的評(píng)價(jià)標(biāo)準(zhǔn)。將評(píng)價(jià)樣本的某一個(gè)指標(biāo)值xl（l=1，2，…，m；m為評(píng)價(jià)指標(biāo)數(shù)）看成一個(gè)集合，同時(shí)把相應(yīng)指標(biāo)的評(píng)價(jià)標(biāo)準(zhǔn)看成另一個(gè)集合Bk（k=1，2，…，K；K為評(píng)價(jià)等級(jí)數(shù)），與Al可構(gòu)成1個(gè)集對(duì)。根據(jù)集對(duì)分析原理，可用聯(lián)系度μl來(lái)描述集對(duì)的關(guān)系：

式中，al為指標(biāo)值與該指標(biāo)第k級(jí)標(biāo)準(zhǔn)的同一度；bl，1為指標(biāo)值與該指標(biāo)第k級(jí)標(biāo)準(zhǔn)相差1級(jí)的差異度；bl，2為指標(biāo)值與該指標(biāo)第k級(jí)標(biāo)準(zhǔn)相差2級(jí)的差異度；bl，K-2為指標(biāo)值與該指標(biāo)第k級(jí)標(biāo)準(zhǔn)相差K-2級(jí)的差異度；cl為指標(biāo)值與該指標(biāo)第k級(jí)標(biāo)準(zhǔn)的對(duì)立度（相差K-1級(jí)）；i1、i2、iK-2分別為第1差異度、第2差異度和第K-2差異度。將Bk特定為某1級(jí)指標(biāo)構(gòu)成的集合，則式（2）各參數(shù)可看成：al為xl隸屬于1級(jí)標(biāo)準(zhǔn)的可能程度；bl，1為指標(biāo)值xl隸屬于 2 級(jí)標(biāo)準(zhǔn)的可能程度；bl，2為指標(biāo)值xl隸屬于3級(jí)標(biāo)準(zhǔn)的可能程度；bl，K-2為指標(biāo)值xl隸屬于K-1級(jí)標(biāo)準(zhǔn)的可能程度；cl為指標(biāo)值xl隸屬于K級(jí)標(biāo)準(zhǔn)的可能程度［9］。

設(shè)評(píng)價(jià)樣本為集合A，1級(jí)評(píng)價(jià)標(biāo)準(zhǔn)為集合B，則集對(duì)的K元聯(lián)系度可定義為：

式中wl為指標(biāo)l的權(quán)重。

儲(chǔ)層定量評(píng)價(jià)一般有K-1個(gè)門限值s1，s2，…，sK-1，對(duì)越小越優(yōu)指標(biāo)，某樣本值xl與該指標(biāo)1級(jí)標(biāo)準(zhǔn)的聯(lián)系度μl為（K=3為例）：

對(duì)越大越優(yōu)指標(biāo)，某樣本值xl與該指標(biāo)1級(jí)標(biāo)準(zhǔn)的聯(lián)系度μl為（K=3為例）：

1.2 確定權(quán)重

采用灰色關(guān)聯(lián)分析法確定權(quán)重的步驟為：① 根據(jù)評(píng)價(jià)目的確定評(píng)價(jià)指標(biāo)體系，收集評(píng)價(jià)數(shù)據(jù)；②原始數(shù)據(jù)均值化處理；③ 確定母序列和子序列；④ 求出關(guān)聯(lián)度系數(shù)和關(guān)聯(lián)度；⑤ 歸一化處理，得出權(quán)重［11-12］。選取某一因素列作為母序列，記為｛xt（0）｝，t=1，2，…，n；因素集合中剩余的列作為子序列，記為｛xt（d）｝，d=1，2，…，m，其中n表示樣本數(shù)，m表示屬性數(shù)目。先對(duì)原始數(shù)據(jù)進(jìn)行均值化處理，得到處理后的無(wú)量綱矩陣，確定母序列和子序列，再運(yùn)用下面的公式計(jì)算關(guān)聯(lián)度系數(shù)：

式中，Lt（d，0）為第d個(gè)因素與母序列的關(guān)聯(lián)度系數(shù)；ρ為分辨系數(shù)，取值范圍為（0，1），通常取0.5；rd，0為第d個(gè)序列與母序列的關(guān)聯(lián)度；wd為各指標(biāo)的權(quán)重。

2 應(yīng)用實(shí)例

2.1 研究區(qū)的地質(zhì)特征

XC氣田構(gòu)造處于四川盆地川西坳陷中段孝泉—豐谷北東東向隆起帶的西段，為孝泉—新場(chǎng)復(fù)式背斜的新場(chǎng)局部圈閉。該氣藏儲(chǔ)層屬低孔、低滲的典型致密儲(chǔ)層，同時(shí)受沉積、成巖以及構(gòu)造等因素的影響，尤其是宏觀、微觀裂縫的差異發(fā)育，使得儲(chǔ)層橫、縱向非均質(zhì)性強(qiáng)，儲(chǔ)集性能差異大［13］。XC氣田須二氣藏儲(chǔ)層無(wú)論層間、層內(nèi)還是平面的非均質(zhì)性均很強(qiáng)。以TX22、TX23、TX24和TX274個(gè)砂層組為例，對(duì)儲(chǔ)層層內(nèi)的非均質(zhì)性進(jìn)行定量評(píng)價(jià)。

2.2 非均質(zhì)性評(píng)價(jià)

取樣TX22、TX23、TX24和TX274個(gè)砂層組，選取與其層內(nèi)非均質(zhì)性相關(guān)的滲透率變異系數(shù)、滲透率突進(jìn)系數(shù)及滲透率級(jí)差3個(gè)評(píng)價(jià)指標(biāo)組成指標(biāo)體系。滲透率變異系數(shù)度量?jī)?chǔ)層滲透率變化的程度，滲透率級(jí)差表征儲(chǔ)層滲透率兩極分化的程度，滲透率突進(jìn)系數(shù)表征滲透率最大值超過(guò)一般水平的程度，計(jì)算公式為：

式中，Vk為滲透率變異系數(shù)；Ki為層內(nèi)某樣品的滲透率（i=1，2，…，n），mD；為層內(nèi)所有樣品滲透率的平均值，mD；Tk為滲透率突進(jìn)系數(shù)；Kmax為層內(nèi)最大滲透率，mD；Jk為滲透率級(jí)差；Kmin為層內(nèi)最小滲透率，mD。

選取滲透率變異系數(shù)列為母序列，運(yùn)用式（7）～（9）求出各因素的權(quán)重w1=（0.43，0.41，0.16），w2=（0.44，0.40，0.16）。研究區(qū)屬于低孔、低滲致密砂巖氣藏，儲(chǔ)層、基質(zhì)的非均質(zhì)性均影響氣藏的合理開(kāi)發(fā)。因此，為定量評(píng)價(jià)儲(chǔ)層內(nèi)非均質(zhì)性的強(qiáng)弱，從儲(chǔ)層非均質(zhì)性參數(shù)（表1）和基質(zhì)非均質(zhì)性參數(shù)（表2）兩方面評(píng)價(jià)儲(chǔ)層的非均質(zhì)性強(qiáng)弱程度。

表1 儲(chǔ)層非均質(zhì)性參數(shù)表

表2 基質(zhì)非均質(zhì)性參數(shù)表

運(yùn)用式（4）計(jì)算可知，砂層組TX22、TX23、TX24層內(nèi)滲透率非均質(zhì)性程度屬嚴(yán)重非均質(zhì)型，而TX27層內(nèi)滲透率非均質(zhì)性程度屬非均質(zhì)型。而從基質(zhì)非均質(zhì)性方面看，TX22和TX24屬于嚴(yán)重非均質(zhì)型，TX23屬于相對(duì)均質(zhì)型，TX27屬于非均質(zhì)型。經(jīng)灰色集對(duì)分析與主成分分析的評(píng)價(jià)結(jié)果對(duì)比，得出兩者的評(píng)價(jià)結(jié)果是一致的。

3 結(jié)束語(yǔ)

結(jié)合灰色集對(duì)分析理論的特點(diǎn)，建立了灰色集對(duì)分析模型，并將其運(yùn)用到儲(chǔ)層非均質(zhì)性的定量評(píng)價(jià)中，并采用灰色關(guān)聯(lián)度確定權(quán)重?；疑瘜?duì)分析涉及的概念清晰，結(jié)構(gòu)簡(jiǎn)單，計(jì)算方便，且計(jì)算結(jié)果可靠。XC氣田須二氣藏的實(shí)際應(yīng)用表明，該方法是可行而有效的，今后可在儲(chǔ)層定量評(píng)價(jià)中加以應(yīng)用。

［1］張興平，衣英杰，夏冰，等.利用多種參數(shù)定量評(píng)價(jià)儲(chǔ)層層間非均質(zhì)性：以尚店油田為例［J］.油氣地質(zhì)與采收率，2004，11（1）：56-57.

［2］何琰.基于模糊綜合評(píng)判與層次分析的儲(chǔ)層定量評(píng)價(jià)：以包界地區(qū)須家河組為例［J］.油氣地質(zhì)與采收率，2011，18（1）：23-25.

［3］張輝，李茂，蔣麗萍，等.儲(chǔ)層非均質(zhì)性評(píng)價(jià)技術(shù)研究及應(yīng)用［J］. 中國(guó)海上油氣，2011，23（3）：175-178.

［4］楊少春.儲(chǔ)層非均質(zhì)性定量研究的新方法［J］.石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2000，24（1）：53-56.

［5］李祖兵，顏其彬，羅明高.非均質(zhì)綜合指數(shù)法在砂礫巖儲(chǔ)層非均質(zhì)性研究中的應(yīng)用：以雙河油田V下油組為例［J］.地質(zhì)科技情報(bào)，2007，26（6）：83-87.

［6］黃延忠，喬勇.利用局部變差函數(shù)描述儲(chǔ)層非均質(zhì)性［J］.石油天然氣學(xué)報(bào)，2010，32（2）：39-42.

［7］徐德英，胡望水，熊平，等.濮城油田沙一段儲(chǔ)層非均質(zhì)性模型研究［J］. 石油天然氣學(xué)報(bào)，2010，32（1）：47-50.

［8］趙克勤.集對(duì)分析及其初步應(yīng)用［M］.杭州：浙江科技出版社，2000.

［9］王文圣，金菊良，丁晶，等.水資源系統(tǒng)評(píng)價(jià)新方法：集對(duì)評(píng)價(jià)法［J］.中國(guó)科學(xué)E輯：技術(shù)科學(xué)，2009，39（9）：1 529-1 534.

［10］程乾生.屬性識(shí)別理論模型及其應(yīng)用［J］.北京大學(xué)學(xué)報(bào)：自然科學(xué)版，1997，33（1）：12-20.

［11］丁熊，譚秀成，羅冰，等.基于灰色模糊理論的多參數(shù)碳酸鹽儲(chǔ)層評(píng)價(jià)［J］.西南石油大學(xué)：自然科學(xué)版，2008，30（5）：88-92.

［12］劉吉余，彭志春，郭曉博.灰色關(guān)聯(lián)分析法在儲(chǔ)層評(píng)價(jià)中的應(yīng)用：以大慶薩爾圖油田北二區(qū)為例［J］.油氣地質(zhì)與采收率，2005，12（2）：13-15.

［13］崔衛(wèi)東，蔣志斌，曾德銘.四川盆地中部上三疊統(tǒng)須家河組儲(chǔ)層控制因素研究［J］.石油天然氣學(xué)報(bào)，2008，30（6）：214-216.