基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法
【專利摘要】本發(fā)明公開了一種基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，包括：根據(jù)巖心實(shí)測數(shù)據(jù)與測井?dāng)?shù)據(jù)計(jì)算得到巖心孔隙度與滲透率的測井解釋模型，然后利用孔隙度、滲透率、地層厚度以及壓縮系數(shù)計(jì)算得到存儲(chǔ)系數(shù)與滲流系數(shù)的雙累積百分?jǐn)?shù)，再利用比值法，計(jì)算每個(gè)測井點(diǎn)的斜率δ值，作出δ曲線，根據(jù)預(yù)先確立的不同的臨界斜率δo值，將δ曲線上對(duì)應(yīng)的測井點(diǎn)劃分不同等級(jí)的“甜點(diǎn)”。本發(fā)明提出的“δ曲線”，是孔滲曲線的綜合響應(yīng)，提出的“δo值”也為動(dòng)態(tài)可變?cè)u(píng)價(jià)標(biāo)準(zhǔn)，簡單而非常高效地一次性多標(biāo)準(zhǔn)的實(shí)現(xiàn)了低滲致密砂巖儲(chǔ)層中單井垂向上的優(yōu)勢“甜點(diǎn)”的劃分與定量評(píng)價(jià)。
【專利說明】
基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明設(shè)及石油開采領(lǐng)域中儲(chǔ)存地質(zhì)學(xué)研究技術(shù)領(lǐng)域，具體設(shè)及一種基于存儲(chǔ)系 數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法。
【背景技術(shù)】
[0002] 我國低滲透儲(chǔ)層的分布區(qū)域廣、儲(chǔ)集巖類型多、資源潛力大，已成為現(xiàn)今油氣勘探 的重點(diǎn)。如吐哈盆地溫吉桑區(qū)塊=工河組與己喀地區(qū)八道灣組、鄂爾多斯盆地的陜北斜坡 中部延長組、準(zhǔn)鳴爾盆地瑪湖凹陷北斜坡區(qū)S疊系百口泉組、珠江口盆地文昌A凹陷古近系 珠海組等低孔-低滲砂巖儲(chǔ)層。
[0003] 已有的"甜點(diǎn)"評(píng)價(jià)技術(shù)有W下幾種，第一種:基于地震屬性切片，結(jié)合高分辨率層 序地層學(xué)與沉積相帶的分布，來預(yù)測平面上"甜點(diǎn)"的分布;第二種:基于巖石物理相的計(jì) 算，建立起巖石物理相"甜點(diǎn)"綜合評(píng)價(jià)指標(biāo)體系，運(yùn)只是一種評(píng)價(jià)指標(biāo);第=種:基于相控 多級(jí)重構(gòu)測井方法，是針對(duì)泥巖提出來的甜點(diǎn)系數(shù)。
[0004] 然而，對(duì)于低滲致密砂巖儲(chǔ)層，很多情況下，"甜點(diǎn)'儲(chǔ)層與"非甜點(diǎn)'儲(chǔ)層的地球物 理響應(yīng)差異小，利用地震資料識(shí)別難度很大，再者，基于地震屬性識(shí)別低滲致密砂巖儲(chǔ)層 "甜點(diǎn)"在區(qū)域上很有效果，但由于地震資料在垂向分辨率的限制，所W無法進(jìn)行單井的垂 向識(shí)別，所W再利用完地震屬性識(shí)別甜點(diǎn)區(qū)域之后，需要進(jìn)一步地利用測井進(jìn)行垂向識(shí)別。
[0005] 對(duì)于基于巖石物理相和成巖作用的識(shí)別方法，利用了多參數(shù)，所需要的資料必須 齊全，且步驟有點(diǎn)復(fù)雜，要求對(duì)地質(zhì)知識(shí)深厚掌握。另外，利用的評(píng)價(jià)參數(shù)是一個(gè)二維概念， 未考慮厚度，缺乏=維轉(zhuǎn)化，并且需要人為地劃分幾類，缺乏自動(dòng)化，且建立的標(biāo)準(zhǔn)，全區(qū)只 有一個(gè)，缺少靈活機(jī)動(dòng)性。
[0006] 因此，如何在垂向井段利用測井資料快速高效地建立可變標(biāo)準(zhǔn)的=維低滲致密儲(chǔ) 層"甜點(diǎn)"剖面的方法與技術(shù)，仍是本領(lǐng)域亟待解決的問題之一。

【發(fā)明內(nèi)容】

[0007] 為了解決上述技術(shù)問題，本發(fā)明提供了一種基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層 甜點(diǎn)評(píng)價(jià)方法，能夠?qū)ξ镄圆町惡艽蟆?甜點(diǎn)儲(chǔ)層"很薄的低滲致密砂巖儲(chǔ)層簡單而有效的 評(píng)價(jià)。
[000引本發(fā)明解決上述技術(shù)問題的技術(shù)方案如下：
[0009]本發(fā)明提供了 一種基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，包括： [0010] Sl、根據(jù)采集取樣的巖屯、實(shí)測數(shù)據(jù)和現(xiàn)阱數(shù)據(jù)計(jì)算得到巖屯、孔隙度和滲透率的測 井解釋模型；
[0011] S2、根據(jù)所述的巖屯、孔隙度和滲透率的測井解釋模型W及地層厚度、壓縮系數(shù)計(jì) 算得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流系數(shù)；
[0012] S3、根據(jù)每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流系數(shù)計(jì)算得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù) 百分?jǐn)?shù)和滲流系數(shù)百分?jǐn)?shù)，進(jìn)而得到每一個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比和累積 滲流系數(shù)百分?jǐn)?shù)的百分比；
[001引S4、利用比值法，計(jì)算每一個(gè)現(xiàn)阱點(diǎn)的斜率S值，作出S曲線；
[0014] S5、根據(jù)預(yù)先確定的臨界斜率So值，將S曲線上對(duì)應(yīng)的測井點(diǎn)劃分為不同等級(jí)的 "甜點(diǎn)'。
[001引本發(fā)明的有益效果為:提出的"S曲線"，是孔滲曲線的綜合響應(yīng)，提出的"50值"也 為動(dòng)態(tài)可變?cè)u(píng)價(jià)標(biāo)準(zhǔn)，簡單而非常高效地一次性多標(biāo)準(zhǔn)的實(shí)現(xiàn)了低滲致密砂巖儲(chǔ)層中單井 垂向上的優(yōu)勢"甜點(diǎn)"的劃分與定量評(píng)價(jià)。
[0016] 在上述技術(shù)方案的基礎(chǔ)上，本發(fā)明還可W做出如下改進(jìn)。
[0017] 進(jìn)一步的，所述步驟Sl中對(duì)研究區(qū)域的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)利用多元線性回 歸方法計(jì)算得到巖屯、孔隙度和滲透率的測井解釋模型：
[0018] O 二a*AC+腳gr+ 丫 *CNL+e*RT+〇;
[0019] L巧=a*AC+腳gr+ 丫 *CNL+e*RT+〇;
[0020] 或者，L巧=a+腳巫；
[00別]式中：O為孔隙度，％ ;K為滲透率，mD;AC為測井聲波時(shí)差值，us/s;C化為中子測井 值，％;GR為測井自然伽馬測井，API;RT為電阻率，歐米;a、e、丫、e、o為系數(shù)，常數(shù)，其數(shù)值由 多元線性回歸得到。
[0022] 進(jìn)一步的，所述步驟S2通過W下方式得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流系數(shù)：
[0023] 存儲(chǔ)系數(shù)：S= O XCtXh;
[0024] 滲流系數(shù):P=K XhAi;
[00巧]厚度h:取測井序列的間隔值，0.1m、0.125m、0.25m、0.5m。
[0026] 其中，Ct為巖石壓縮系數(shù)，1/M化;h為測井序列的間隔值，h為常數(shù)，按照不同的測 井深度，h值不同，y為地層原油粘度，為常數(shù)。
[0027] 進(jìn)一步的，所述步驟S3中通過如下方式得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)百分?jǐn)?shù)和滲 流系數(shù)百分?jǐn)?shù)：
[0028] 垂向井段所有的測井點(diǎn)的孔隙度O、滲透率K、地層厚度h、存儲(chǔ)系數(shù)P和滲流系數(shù)S 組成如下數(shù)據(jù)矩陣：
[0029]
[0030] 按照如下公式計(jì)算垂向井段每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)百分?jǐn)?shù)和滲流系數(shù)百分?jǐn)?shù)：
[0031]
[0032] 其中，Ki--第i個(gè)測井點(diǎn)的滲透率；
[00削巫i--第i個(gè)測井點(diǎn)的孔隙度；
[0034] hi一一第i個(gè)測井點(diǎn)對(duì)應(yīng)的地層厚度；
[0035] Si---第i個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)；
[0036] Pi---第i個(gè)測井點(diǎn)的滲流系數(shù)；
[0037] SSi---第i個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)的百分?jǐn)?shù)；
[0038] PPi---第i個(gè)測井點(diǎn)的滲流系數(shù)的百分?jǐn)?shù)；
[0039] n---測井點(diǎn)的總個(gè)數(shù)。
[0040] 進(jìn)一步的，通過如下方式計(jì)算每一個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比和累 積滲流系數(shù)百分?jǐn)?shù)的百分比：
[0041]
[0042] 其中，Sri---第i個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比；
[0043] Pri---第i個(gè)測井點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)的百分比。
[0044] 進(jìn)一步的，所述臨界斜率So值根據(jù)采油工藝與技術(shù)手段動(dòng)態(tài)可調(diào)。
【附圖說明】
[0045] 圖1為本發(fā)明實(shí)施例一的一種基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法 流程圖；
[0046] 圖2為實(shí)施例一中提供的自動(dòng)化運(yùn)算程序運(yùn)行時(shí)的界面圖；
[0047] 圖3為實(shí)施例二提供的垂向井段的低滲儲(chǔ)層"甜點(diǎn)"綜合評(píng)價(jià)效果圖。
【具體實(shí)施方式】
[0048] W下結(jié)合附圖對(duì)本發(fā)明的原理和特征進(jìn)行描述，所舉實(shí)例只用于解釋本發(fā)明，并 非用于限定本發(fā)明的范圍。
[0049] 實(shí)施例一、一種基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法。
[0050] 參見圖1，Sl、根據(jù)采集取樣的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)計(jì)算得到巖屯、孔隙度和滲 透率的測井解釋模型。
[0051] 具體的，首先簡單介紹一下低滲透率"甜點(diǎn)"的概念，"甜點(diǎn)"是多種地質(zhì)作用的綜 合反映，是沉積作用、成巖作用、后期構(gòu)造作用和流體改造作用形成的有利巖石物理相流動(dòng) 單元，它們具有相對(duì)優(yōu)質(zhì)儲(chǔ)層的滲流、儲(chǔ)集及含油氣與非均質(zhì)特征，集中地反映該區(qū)特低滲 透儲(chǔ)層含油有利區(qū)形成的地質(zhì)特點(diǎn)。是專口針對(duì)低滲儲(chǔ)層、致密儲(chǔ)層的專業(yè)"術(shù)語"。
[0052] 在整個(gè)測井過程中，可W先將單井的一部分?jǐn)?shù)據(jù)作為采樣數(shù)據(jù)，根據(jù)采樣的巖屯、 實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)計(jì)算擬合得到巖屯、孔隙度和滲透率的測井解釋模型。在本實(shí)施例中， 所說的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)主要包括測井聲波時(shí)差值A(chǔ)C,單位為us/s，中子測井值 CNL，測井自然伽瑪測井測井值GR，單位為API,電阻率，單位為歐米。垂向井段的每一個(gè)測井 點(diǎn)均有對(duì)應(yīng)的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)。
[0053] 其中，可W采用多種方法根據(jù)采樣測井點(diǎn)的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)計(jì)算得出巖 屯、孔隙度和滲透率的測井解釋模型。在本實(shí)施例中，可W采樣測井點(diǎn)的巖屯、實(shí)測數(shù)據(jù)和測 井?dāng)?shù)據(jù)利用多元線性回歸方法計(jì)算得到巖屯、孔隙度和滲透率的測井解釋模型：
[0054] O =a*AC+腳gr+丫 *CNL+e*RT+〇; (1)
[0055] L巧=a*AC+腳gr+丫 *CNL+e*RT+〇; (2)
[0化6] Lg = Q+腳巫； (3)
[0化7] 式中：巫為孔隙度，％ ;K為滲透率，mD;AC為測井聲波時(shí)差值，us/s;C化為中子測井 值，％;GR為測井自然伽馬測井，API;RT為電阻率，歐米;a、e、丫、e、o為系數(shù)，常數(shù)，其數(shù)值由 多元線性回歸得到;K1、K2為滲透率的兩種計(jì)算方法，選用更為準(zhǔn)確的計(jì)算方法。
[0058] 需要說明的是，在采樣的數(shù)據(jù)中，每一個(gè)采樣測井點(diǎn)的孔隙度?、滲透率K，W及每 一個(gè)采樣測井點(diǎn)的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)（即4(：八化、61?、1^)都是已知的，因此，根據(jù)上 述公式（1)、（2)和(3)可得出擬合系數(shù)a、e、y、e、o。對(duì)于計(jì)算滲透率的兩個(gè)公式(2)和(3)， 將計(jì)算出來的擬合系數(shù)a、e、丫、e、〇代入公式(2)和(3)中，得到計(jì)算滲透率的兩個(gè)公式，根 據(jù)兩個(gè)公式的多元線性回歸擬合度，來確定公式(2)和公式(3)的擬合準(zhǔn)確性，即分別計(jì)算 AC、GR、CN和LRT運(yùn)四個(gè)參數(shù)在不同的擬合系數(shù)的條件下與滲透率K的相關(guān)性，相關(guān)性更高的 計(jì)算公式為更為準(zhǔn)確的滲透率計(jì)算公式，因此，本實(shí)施例中提出了兩種滲透率的計(jì)算公式， 選擇相關(guān)性更高的計(jì)算公式為滲透率的計(jì)算公式。
[0059] S2、根據(jù)所述的巖屯、孔隙度和滲透率的測井解釋模型W及地層厚度、壓縮系數(shù)計(jì) 算得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流系數(shù)；
[0060] 具體的，本步驟根據(jù)上述步驟1得出的巖屯、孔隙度和滲透率的測井解釋模型，并根 據(jù)整個(gè)垂向井段的每一個(gè)測井點(diǎn)的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)，計(jì)算出每一個(gè)測井點(diǎn)的巖屯、 孔隙度和滲透率。
[0061] 然后，根據(jù)整個(gè)垂向井段的每一個(gè)測井點(diǎn)的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)的巖屯、孔隙 度和滲透率，W及研究井段的地層厚度和預(yù)設(shè)好的巖石壓縮系數(shù)計(jì)算出每一個(gè)測井點(diǎn)的存 儲(chǔ)系數(shù)和滲流系數(shù)。本實(shí)施例中，可通過W下方公式計(jì)算每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流 系數(shù)：
[0062] 存儲(chǔ)系數(shù)：S=巫 XCtXh; (4)
[0063] 滲流系數(shù):P=KXhAi; (5)
[0064] 其中，地層厚度h:取測井序列的間隔值，0.1m、0.125m、0.25m、0.5m。
[0065] Ct為巖石壓縮系數(shù)，1/M化;h為測井序列的間隔值，h為常數(shù)，按照不同的測井深 度，h值不同，y為地層原油粘度，mPa. S，或者厘泊。
[0066] S3、根據(jù)每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流系數(shù)計(jì)算得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù) 百分?jǐn)?shù)和滲流系數(shù)百分?jǐn)?shù)，進(jìn)而得到每一個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比和累積 滲流系數(shù)百分?jǐn)?shù)的百分比。
[0067] 具體的，根據(jù)上述步驟2的公式(4)和公式(5)可W計(jì)算出每一個(gè)測井點(diǎn)的存儲(chǔ)系 數(shù)和滲流系數(shù)，并可通過如下的公式計(jì)算每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)百分?jǐn)?shù)和滲流系數(shù)百分 數(shù)：
[0068] 根據(jù)每一個(gè)測井點(diǎn)的巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)，W及結(jié)合前述公式（1)、公式(2)、 公式(3)和公式(4)，可W計(jì)算出每一個(gè)測井點(diǎn)的巖屯、孔隙度?、滲透率K、存儲(chǔ)系數(shù)S和滲流 系數(shù)P，因此，對(duì)于研究井段有如下數(shù)據(jù)矩陣：
[0069]
'、 "X 、 '* / 、 "' / ± . ) \ " " ? " y X <* X K " ?'t "X V n /
[0070] 然后可按照如下公式計(jì)算垂向井段每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)百分?jǐn)?shù)和滲流系數(shù) 百分?jǐn)?shù)：
[0071]
[0072] 其中，Si---第i個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)；
[0073] Pi---第i個(gè)測井點(diǎn)的滲流系數(shù)；
[0074] SSi---第i個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)的百分?jǐn)?shù)；
[0075] PPi---第i個(gè)測井點(diǎn)的滲流系數(shù)的百分?jǐn)?shù)；
[0076] n---測井點(diǎn)的總個(gè)數(shù)。
[0077] 最后，計(jì)算出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)百分?jǐn)?shù)和滲流系數(shù)百分?jǐn)?shù)之后，可按照如 下計(jì)算公式計(jì)算出垂向井段每一個(gè)測井點(diǎn)的累計(jì)存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比和累計(jì)滲流系 數(shù)百分?jǐn)?shù)的百分比：
[007引
[0079] 其中，Sri---第i個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比；
[0080] Pri---第i個(gè)測井點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)的百分比。
[0081] S4、利用比值法，計(jì)算每一個(gè)測井點(diǎn)的斜率S值，作出S曲線。
[0082] 具體的，計(jì)算出了垂向井段每一個(gè)測井點(diǎn)的累計(jì)存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比和累計(jì) 滲流系數(shù)百分?jǐn)?shù)的百分比之后，可通過如下公式計(jì)算每一個(gè)測井點(diǎn)的斜率S值：
[0083]
[0084] 其中，Si---第i點(diǎn)的斜率；
[00化]Hi---第i個(gè)測井點(diǎn)；
[0086] Sri---第i個(gè)測井深度點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比；
[0087] Pri---第i個(gè)測井深度點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)的百分比；
[0088] Sri+1---第i+1個(gè)測井深度點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比；
[0089] Pri+1---第i+1個(gè)測井深度點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)的百分比。
[0090] 根據(jù)上述公式可求出每一個(gè)測井點(diǎn)的斜率S值，即S1，S2，的，M……Sn，作出"S曲 線"。
[0091] S5、根據(jù)預(yù)先確定的臨界斜率So值，將S曲線上對(duì)應(yīng)的測井點(diǎn)劃分為不同等級(jí)的 "甜點(diǎn)'。
[0092] 具體的，上述步驟S4繪制出了每一個(gè)測井點(diǎn)的"S曲線"，然后根據(jù)預(yù)先確定的臨界 斜率So值，將整個(gè)垂向井段的每一個(gè)測井點(diǎn)為劃分不同等級(jí)的"甜點(diǎn)"。
[0093] 其中，臨界斜率So值根據(jù)采油工藝與技術(shù)手段動(dòng)態(tài)可調(diào)，比如，采油工藝復(fù)雜或者 所使用的技術(shù)手段比較高，貝柯W將臨界斜率M直設(shè)定得大一些，運(yùn)樣對(duì)于劃分出來的"甜 點(diǎn)"等級(jí)會(huì)更為準(zhǔn)確一些。
[0094] 需要說明的是，根據(jù)現(xiàn)有的很多實(shí)踐經(jīng)驗(yàn)，在低滲儲(chǔ)層中，通常有如下表1的規(guī)律：
[0095] 表 1
[0096]
[0097] 實(shí)施例二、
[009引為了加深對(duì)實(shí)施例一提供的低滲儲(chǔ)層"甜點(diǎn)"評(píng)價(jià)方法的進(jìn)一步理解，本實(shí)施例提 供了一個(gè)具體的例子來說明。
[0099] 涵洲11-7油田位于北部灣盆地涵西南凹陷南部斜坡、涵西南低凸起向東北傾末 端，為3號(hào)扭張斷裂帶東帶上的斷塊構(gòu)造。涵洲11-7油田包括涵洲11-4N區(qū)塊、涵洲11-7區(qū)塊 及涵洲11-8區(qū)塊，油田面積約200km2。目的層段主要沉積體系為緩坡型扇S角洲沉積體系 (主要為流=段)和近源型辮狀河=角洲沉積體系（主要為流一段），伴生沉積體系為重力流 沉積體系，湖泊沉積體系為沉積背景相，另外在流=段底部發(fā)育粗粒陸上沖積扇沉積體系。 根據(jù)對(duì) WZl 1-7-1，WZl 1-7-2, WZl l-7-2Sa，WZl 1-7-3, WZl 1-7-4, WZl 1-8-2及WZl 1-4N-6運(yùn)7 口 井近千個(gè)樣品的常規(guī)物性統(tǒng)計(jì)與分析，研究區(qū)流一段與流=段的孔隙度、滲透率數(shù)值差異 較大。
[0100] 流一段孔隙度在10~25 %之間的樣品占80 % W上，滲透率有兩個(gè)分布區(qū)間：一個(gè) 介于（0.1~10) X 10-3皿2之間，運(yùn)部分樣品占43 %左右，另外一個(gè)100 X IO-Smi2W上，占 41 %左右。從整個(gè)流一段物性資料來看，其孔隙度在1.87~27.46%之間，平均為15.87%， 滲透率在(0.001~23327) X 10-3皿2之間，平均為1061.4X 10-3皿2。也就是說，流一段W中 低孔、特低滲儲(chǔ)層為主，同時(shí)中滲、高滲甚至特高滲儲(chǔ)層均有發(fā)育，流=段80%的樣品孔隙 度小于15% ,90% W上的樣品滲透率小于IOX 10-3皿2,因此流S段發(fā)育的主要是低孔特低 滲儲(chǔ)層。通過對(duì)流=段物性資料的統(tǒng)計(jì)，其孔隙度主要分布在0.13~21.91%之間，平均為 11.77%，滲透率分布在(0.01~34.6) X 10-3皿2之間，平均為5.56 X 10-3皿2。不同的油層 組，滲透率變化不大，均W特低滲儲(chǔ)層占絕對(duì)優(yōu)勢。因此，涵洲11-7油田流沙港組=段砂巖 的儲(chǔ)層品質(zhì)較差，總體屬于低孔低滲砂巖?，F(xiàn)在就對(duì)WZ11-7-1井的流=段進(jìn)行低滲儲(chǔ)層"甜 點(diǎn)"識(shí)別工作。
[0101] 步驟SI,巖屯、孔隙度和滲透率的測井解釋模型建立，具體方法為：
[0102] 本實(shí)施例根據(jù)巖屯、實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)，利用多元回歸方法得到如下巖屯、孔隙度 和滲透率的測井解釋模型：巫=0.04188+0.5709CNL+0.00163AC-0.0 Ol 12GR+0.0002397RT; R = O.7366
[0103] LnK = -IO. 5649巧.6254CNL+0.1849AC-0.0476GR+0.07：MRT ;R = 0.7621
[0104] 步驟S2:累積存儲(chǔ)系數(shù)與滲流系數(shù)的計(jì)算。
[0105] 1)每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)、滲流系數(shù)計(jì)算。
[0106] 依據(jù)公式:滲流系數(shù)= P = KXhAi;存儲(chǔ)系數(shù):S= O XCt Xh;分別計(jì)算每個(gè)測井點(diǎn) 的P、s。
[0107] 依照實(shí)驗(yàn)測試，該層位未發(fā)育超壓，巖石壓縮系數(shù)Ct變化不大，取0.23 X 10-4MPa- 1。原油的粘度y為0.8厘泊。厚度h取測井系列間隔0.1m，如下表2給出了計(jì)算結(jié)果與數(shù)據(jù)形 式：
[010引 表2
[0109]
[0110]
[0111] 2)累積存儲(chǔ)系數(shù)與滲流系數(shù)的計(jì)算。
[0112] 先計(jì)算全井段的累積滲流系數(shù)、存儲(chǔ)系數(shù)之和，再逐個(gè)計(jì)算每個(gè)測井點(diǎn)的累積滲 流系數(shù)百分?jǐn)?shù)與累積存儲(chǔ)系數(shù)百分?jǐn)?shù)，按如下公式：
[0113]
[0114] 巧甲，
[0115] Si---第i個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)；
[0116] Pi---第i個(gè)測井點(diǎn)的滲流系數(shù)；
[0117] SSi---第i個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)的百分?jǐn)?shù)；
[0118] PPi---第i個(gè)測井點(diǎn)的滲流系數(shù)的百分?jǐn)?shù)；
[0119] n---研究層位測井點(diǎn)的總個(gè)數(shù)。
[0120] 最后，計(jì)算垂向井段每一個(gè)測井點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)的百分比、累積存儲(chǔ)系 數(shù)百分撕的巧.格々n下A管/A井.
[0121]
[0122] 式中，
[0123] Sri---第i個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比；
[0124] Pri---第i個(gè)測井點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)的百分比。
[0125] 表3給出了本實(shí)施例的各項(xiàng)計(jì)算結(jié)果：
[0126] 表3
[0127]
[012 引
[0129] 步驟S3: "S曲線'的繪制。
[0130] if當(dāng)每個(gè)娜I#點(diǎn)對(duì)欣的魁莖S值，利用如下公式計(jì)算；
[0131]
[0132] 式中，
[0133] Si---第i點(diǎn)的斜率；
[0134] 化---第i個(gè)測井點(diǎn)深度；
[0135] Sri---第i個(gè)測井深度點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)；
[0136] Pri---第i個(gè)測井深度點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)；
[0137] Sri+1---第i+1個(gè)測井深度點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)；
[0138] Pri+1---第i+1個(gè)測井深度點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)；
[0139] 根據(jù)上述公式求出每一個(gè)測井點(diǎn)的S值，即Sl，52,的，M……細(xì)，作出"S曲線"。
[0140] 表4,給出了數(shù)據(jù)計(jì)算結(jié)果與形式。
[0141] 表4

[0143] 步驟S4:依據(jù)So值大小評(píng)價(jià)儲(chǔ)層"甜點(diǎn)"等級(jí)。
[0144] 根據(jù)生產(chǎn)實(shí)踐與經(jīng)濟(jì)技術(shù)標(biāo)準(zhǔn)，確定儲(chǔ)層"甜點(diǎn)'的臨界值"MT,依據(jù)M直劃分 "甜點(diǎn)"儲(chǔ)層的等級(jí)。此處臨界值"So值標(biāo)準(zhǔn)"，每口井可隨時(shí)變動(dòng)，所W可W依據(jù)經(jīng)濟(jì)技術(shù)與 采油工藝的發(fā)展，隨時(shí)更換"M直標(biāo)準(zhǔn)"的大小，程序運(yùn)行圖，如圖2。
[0145] 此程序運(yùn)用VB代碼編譯（該算法思想，也可W通過C/C++、matlab編譯均可），程序 設(shè)有四個(gè)按鈕，物性計(jì)算按鈕對(duì)應(yīng)步驟SI,累計(jì)存儲(chǔ)系數(shù)滲流系數(shù)計(jì)算按鈕對(duì)應(yīng)步驟S2和 S3，"S值"的計(jì)算按鈕對(duì)應(yīng)步驟S4，"甜點(diǎn)"儲(chǔ)層等級(jí)評(píng)價(jià)按鈕對(duì)應(yīng)步驟S5。程序運(yùn)行時(shí)，先將 深度與測井?dāng)?shù)據(jù)集放到Excel表格內(nèi)，然后逐個(gè)點(diǎn)擊按鈕，計(jì)算產(chǎn)生的結(jié)果會(huì)在原始數(shù)據(jù)后 面的相應(yīng)的"列"一一顯示，形式參見表2、表3、表4,最后利用此結(jié)果進(jìn)行作圖分析。
[0146] -般性的，根據(jù)現(xiàn)有很多實(shí)踐經(jīng)驗(yàn)，在低滲儲(chǔ)層中，有如下規(guī)律：
[0147]
[0148] 綜合評(píng)價(jià)結(jié)果，如圖3。其中，評(píng)價(jià)的儲(chǔ)層"甜點(diǎn)"與測井?dāng)?shù)據(jù)對(duì)應(yīng)，并且深度已達(dá)到 0.1米的數(shù)量級(jí)，有利于儲(chǔ)層單井垂向上的精細(xì)刻畫。
[0149] 從的解釋結(jié)果來看，通過引進(jìn)"S曲線"與"So值"標(biāo)準(zhǔn)，將儲(chǔ)層的"甜點(diǎn)"在單井上細(xì) 分了成了I、ii、m類，并且考慮了地層厚度的大小。如果每口井定義不同的So值，可W獲得 更多的不同標(biāo)準(zhǔn)的"甜點(diǎn)"評(píng)價(jià)指標(biāo)，避免了全區(qū)使用一個(gè)分類評(píng)價(jià)標(biāo)準(zhǔn)的局限。
[0150] 本發(fā)明提供的一種基于存儲(chǔ)系數(shù)與滲流系數(shù)的低滲儲(chǔ)層"甜點(diǎn)"評(píng)價(jià)方法，與現(xiàn)有 技術(shù)相比，帶來的有益效果主要包括：
[0151] (1)引進(jìn)了 "S曲線"，將低滲儲(chǔ)層的存儲(chǔ)能力與滲流能力融為一條曲線，避免了傳 統(tǒng)的依據(jù)傾斜程度進(jìn)行分段，減少人為因素誤差，更能直觀地判別"甜點(diǎn)"儲(chǔ)層的等級(jí)；引進(jìn) 值，"甜點(diǎn)"只是一個(gè)相對(duì)的概念，而剛好用斜率5值，即變化梯度來很好的表征，賦予了 "甜點(diǎn)"數(shù)學(xué)的意義，便于量化計(jì)算；
[0152] (2)引進(jìn)"M直"，為可變的臨界值，可W為不同的井、不同的開發(fā)時(shí)期制定不同的 "甜點(diǎn)"識(shí)別標(biāo)準(zhǔn)，運(yùn)樣即克服了整個(gè)油田區(qū)塊單一的劃分方案；
[0153] (3)利用測井?dāng)?shù)據(jù)直接識(shí)別低滲致密儲(chǔ)層"甜點(diǎn)"，在地震屬性識(shí)別完W后進(jìn)一步 識(shí)別，提高了其垂向上的識(shí)別率；
[0154] (4)利用加入?yún)?shù)h(厚度)之后，使得"甜點(diǎn)"評(píng)價(jià)更顯得=維化，克服了只利用孔 隙度和滲透率二維平面概念識(shí)別"甜點(diǎn)"的不足；
[0155] (5)該方法步驟簡單，可W對(duì)全井解釋實(shí)現(xiàn)全面自動(dòng)化，為進(jìn)一步推動(dòng)儲(chǔ)層地質(zhì)評(píng) 價(jià)自動(dòng)化、一體化具有一定的指導(dǎo)意義。
[0156] 在本說明書的描述中，參考術(shù)語"實(shí)施例一"、"示例"、"具體示例"、或"一些示例" 等的描述意指結(jié)合該實(shí)施例或示例描述的具體方法、裝置或者特點(diǎn)包含于本發(fā)明的至少一 個(gè)實(shí)施例或示例中。在本說明書中，對(duì)上述術(shù)語的示意性表述不必須針對(duì)的是相同的實(shí)施 例或示例。而且，描述的具體特征、方法、裝置或者特點(diǎn)可W在任一個(gè)或多個(gè)實(shí)施例或示例 中W合適的方式結(jié)合。此外，在不相互矛盾的情況下，本領(lǐng)域的技術(shù)人員可W將本說明書中 描述的不同實(shí)施例或示例W及不同實(shí)施例或示例的特征進(jìn)行結(jié)合和組合。
[0157] W上所述僅為本發(fā)明的較佳實(shí)施例，并不用W限制本發(fā)明，凡在本發(fā)明的精神和 原則之內(nèi)，所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本發(fā)明的保護(hù)范圍之內(nèi)。
【主權(quán)項(xiàng)】
1. 一種基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，其特征在于，包括： 51、 根據(jù)采集取樣的巖心實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)計(jì)算得到巖心孔隙度和滲透率的測井解 釋模型； 52、 根據(jù)所述的巖心孔隙度和滲透率的測井解釋模型以及地層厚度、壓縮系數(shù)計(jì)算得 出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流系數(shù)； 53、 根據(jù)每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流系數(shù)計(jì)算得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)百分 數(shù)和滲流系數(shù)百分?jǐn)?shù)，進(jìn)而得到每一個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比和累積滲流 系數(shù)百分?jǐn)?shù)的百分比； 54、 利用比值法，計(jì)算每一個(gè)測井點(diǎn)的斜率δ值，作出δ曲線； 55、 根據(jù)預(yù)先確定的臨界斜率δ〇值，將δ曲線上對(duì)應(yīng)的測井點(diǎn)劃分為不同等級(jí)的"甜點(diǎn)"。2. 如權(quán)利要求1所述的基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，其特征在 于，所述步驟Sl中對(duì)采集取樣的巖心實(shí)測數(shù)據(jù)和測井?dāng)?shù)據(jù)利用多元線性回歸方法計(jì)算得到 巖心孔隙度和滲透率的測井解釋模型： Φ =a*AC+P*gr+ γ *CNL+e*RT+〇; LgK = a*AC+P*gr+ γ *CNL+e>!<RT+〇; 或者，LgK = α+β*Φ ; 式中：Φ為孔隙度，％ ;Κ為滲透率，mD; AC為測井聲波時(shí)差值，us/s; CNL為中子測井 值，％;GR為測井自然伽馬測井，API;RT為電阻率，歐米;α、β、γ、ε、σ為系數(shù)，常數(shù)，其數(shù)值由 多元線性回歸得到。3. 如權(quán)利要求2所述的基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，其特征在 于，所述步驟S2通過以下方式得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)和滲流系數(shù)： 存儲(chǔ)系數(shù):S=O XCtXh; 滲流系數(shù):P = K X h/μ; 厚度h:取測井序列的間隔值，0.1m、0.125m、0.25m、0.5m; 其中，Ct為巖石壓縮系數(shù)，l/MPa;h為測井序列的間隔值，h為常數(shù)，按照不同的測井深 度，h值不同，μ為地層原油粘度，為常數(shù)。4. 如權(quán)利要求3所述的基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，其特征在 于，所述步驟S3中通過如下方式得出每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)百分?jǐn)?shù)和滲流系數(shù)百分?jǐn)?shù)： 垂向井段所有的測井點(diǎn)的孔隙度Φ、滲透率Κ、地層厚度h、存儲(chǔ)系數(shù)P和滲流系數(shù)S組成 如下數(shù)據(jù)矩陣：按照如下公式計(jì)算垂向井段每一個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)百分?jǐn)?shù)和滲流系數(shù)百分?jǐn)?shù)： 其中，Ki--第i個(gè)測井點(diǎn)的滲透卒；Φ?一一第i個(gè)測井點(diǎn)的孔隙度； hi-一第i個(gè)測井點(diǎn)對(duì)應(yīng)的地層厚度； Si---第i個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)； Pi---第i個(gè)測井點(diǎn)的滲流系數(shù)； SSi---第i個(gè)測井點(diǎn)的存儲(chǔ)系數(shù)的百分?jǐn)?shù)； PPi---第i個(gè)測井點(diǎn)的滲流系數(shù)的百分?jǐn)?shù)； η 測井點(diǎn)的總個(gè)數(shù)。5. 如權(quán)利要求4所述的基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，其特征在 于，通過如下方式計(jì)算每一個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比和累積滲流系數(shù)百分 數(shù)的百分比：其中，Sri---第i個(gè)測井點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)的百分比； Pri---第i個(gè)測井點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)的百分比。6. 如權(quán)利要求5所述的基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，其特征在 于，所述步驟S4中所述的利用比值法，計(jì)算每一個(gè)測井點(diǎn)的斜率δ值具體包括： 其中，δ?---第i點(diǎn)的_ ,,Hi---第i個(gè)測井點(diǎn)； Sri---第i個(gè)測井深度點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)； Pri---第i個(gè)測井深度點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)； Sri+1---第i+1個(gè)測井深度點(diǎn)的累積存儲(chǔ)系數(shù)百分?jǐn)?shù)； Pri+1---第i+1個(gè)測井深度點(diǎn)的累積滲流系數(shù)百分?jǐn)?shù)。7.如權(quán)利要求1所述的基于存儲(chǔ)系數(shù)和滲流系數(shù)的低滲儲(chǔ)層甜點(diǎn)評(píng)價(jià)方法，其特征在 于，所述臨界斜率So值根據(jù)采油工藝與技術(shù)手段動(dòng)態(tài)可調(diào)。
【文檔編號(hào)】G06F17/50GK106021793SQ201610383330
【公開日】2016年10月12日
【申請(qǐng)日】2016年6月1日
【發(fā)明人】姚光慶, 陳孝君, 谷明輝, 龍海清, 黃銀濤, 毛千, 樊曉依
【申請(qǐng)人】中國地質(zhì)大學(xué)（武漢）