Лабораторна робота №4



Виділення та промислова оцінка нафтоносних та газоносних колекторів (визначення ефективної товщини, коефіцієнта нафто-газонасичення та залишкового нафтогазонасичення) за даними комплексу методів ГДС

 

Мета роботи. Ознайомитись із методикою визначення коефіцієнта нафтогазонасичення, залишкового нафтогазонасичення та ефективної товщини за даними комплексу методів ГДС.

4.1 Теоретичні відомості

Визначення коефіцієнта нафтогазонасичення

Коефіцієнти нафтонасичення K н і газонасичення K г порід визначаються за параметром насичення Р н(Р г)= r нп/ r вп. Для розрахунку Р н необхідно знати питомий опір нафтоносної або газоносної породи r нп, що досліджується, та її опір r вп при 100 %-му насиченні пор пластовою водою.

Опір r нп нафтоносної породи визначається за допомогою діаграм уявного та ефективного опорів. Опір r вп розраховується за даними коефіцієнта пористості K п породи та опору пластових вод r в:

 

, (4.1)

 

або [якщо колекторські властивості нафтоносного (газоносного) об’єкта достатньо стабільні] приймається рівним його опору за межами контуру нафтоносності (газоносності) і, зокрема, на ділянках колектора, що знаходиться нижче початкового водонафтового контакту.

Після визначення r нп і r вп розраховується величина параметра насичення:

 

. (4.2)

 

Величина параметра насичення практично не залежить від коефіцієнта пористості колекторів даного типу і є обернено-степеневою функцією коефіцієнта водонасичення K в породи:

, (4.3)

 

де а і n – структурні елементи, які залежать від типу колектора, структури порового простору та типу заповнювача пор. При цьому показник степеня різко збільшується з переходом від гідрофільних до гідрофобних колекторів.

Коефіцієнт водонасичення пов'язаний із коефіцієнтом нафтогазонасичення наступним чином:

 

. (4.4)

 

За розрахованим значенням Р н визначаються коефіцієнти K н або K г (K нг=1- K в) (при атмосферному тиску) за номограмою (Рис. 4.1).

Зокрема, наприклад, для Р н=90 будуть отримані наступні значення K в і K н:

1) карбонатний колектор – K в=7.4 %, K н=92.6 %;

2) піщано-глинистий гідрофільний колектор – K в=11 %, K н=89 %;

3) піщано-глинистий слабо-гідрофобний колектор – K в=16.1 %, K н=83.9 %;

4) піщано-глинистий гідрофобний колектор – K в=26.5 %, K н=73.5 %.

Як видно з цього прикладу, неврахування гідрофобності колектора може призвести до значних похибок при визначенні коефіцієнта нафтонасичення колектора. Другим фактором, який понижує точність визначення коефіцієнта K н, є глинистість.

У тому випадку, коли тонкі глинисті прошарки чергуються з нафто- і газоносними піщаниками та свердловина розкриває пласти з кутом зустрічі близьким до 90°, за даними виміру уявного опору великими зондами визначається повздовжній питомий опір r t, н пачки. При цьому для нафтогазоносної глинистої породи:

 

1, 2, 3 – відповідно, піщано-глинисті гідрофільні, слабо-гідрофобні та гідрофобні породи; 4 – карбонатні породи

Рисунок 4.1 – Залежність параметра насичення P н

від коефіцієнта водонасичення К в

(нафтогазонасичення K нг) породи

, (4.5)

 

де r гл, r нп – відповідно питомий опір глин і нафтогазонасичених пісковиків потужністю h п; c гл – вміст глинистих прошарків сумарною потужністю h гл [ c гл= h гл/(h гл+ h п)] у сумарній потужності досліджуваного об’єкта, яка рівна h гл+ h п.

У тому випадку, коли пачка досліджуваних піщано-глинистих відкладів водоносна, за даними виміру уявного опору великими градієнт-зондами визначається питомий опір:

 

, (4.6)

 

де r t, в – питомий повздовжній опір пачки.

Таким чином, у глинистих пісках і піщаниках величина відношення:

 

(4.7)

 

не дорівнює параметру насичення Р н.

Розв’язуючи рівняння (4.7) відносно Р н= r нп/ r вп, одержимо:

 

. (4.8)

З формули (4.8) випливає, що дійсне значення Р н завжди більше величини Р н, г розрахованої за даними відношень питомих опорів, які виміряні у свердловині в нафтогазоносному і водоносному колекторах.

У тому випадку, коли точне значення r вп невідоме, при c гл<0.5 і 1- r вп/ r гл < 0.2, а також при c гл<0.25 і 0.2<1- r вп/ r гл<0.5 або c гл<0,1 і 0.5<1- r вп/ r гл<1.5 розрахунок Р н може бути виконаний за наближеною формулою:

 

. (4.9)

 

За рубежем коефіцієнт водонасичення нафтоносних колекторів, у тому числі глинистих, обчислюють за різними формулами. У більшості випадків у них враховується шарувате залягання глинистих мінералів і обернено-квадратична залежність між Р н і K в. До таких формул відносяться:

для частого чергування піщаних і глинистих прошарків:

 

, (4.10)

 

а також наступні рівняння для розсіяного та шаруватого включення глинистих мінералів:

Хоссена:

 

, (4.11)

 

Долля:

 

, (4.12)

 

Пупона і Лоу:

 

, (4.13)

 

Симанду:

 

(4.14)

 

a гл – коефіцієнт, що враховує особливості впливу провідності глин і лінійно змінюється від 0 до 1;

Де Вітте:

 

(4.15)

 

Р. Десбрант на підставі аналізу даних по одній зі свердловин, які пробурені у Сахарі, показує, що величини K в, які визначені різними методами, мають близькі значення. За методом Хоссена в більшості випадків одержують завищені K в, по методу Сіманду – найбільш близькі до дійсного.

Визначення коефіцієнтів нафтогазонасичення за даними нейтронних методів

У основу визначення коефіцієнта нафтонасичення за даними нейтронних методів покладені залежності інтенсивності випромінювання від змісту в пластових водах елементів з аномальними нейтронними властивостями (наприклад хлору). Однак істотне ускладнення виникає через. невелику глибинність дослідження цими методами. Ця обставина при наявності зони проникнення фільтрату промивної рідини, концентрація солей у якому істотно відрізняється від їх концентрації в пластових водах, створює труднощі для внесення відповідних виправлень.

При достатній мінералізації пластових вод за хлором найбільш інформативні імпульсні методи. При зазначеній умові нейтронна поглинаюча активність перебуває в прямій залежності від коефіцієнта водонасичення K в порового простору.

Коефіцієнт K н за даними імпульсного методу одержують при рішенні наступного рівняння:

 

, (4.16)

 

де величини A n п нп, A n п ск, A n п в і An п ве – відповідно, нейтронно-поглинаюча активність нафтогазоносного пласта, скелету породи, води та вуглеводневих елементів.

З рівняння (4.16) знаходимо:

 

. (4.17)

 

Якщо в розрізі свердловини є водонасичений об’єкт із таким ж мінеральним складом і пористістю, як і нафтонасичений пласт, то рішення завдання істотно полегшується. Для цього об’єкта будемо мати:

 

. (4.18)

Різниця нейтронних поглинаючих активностей у водонасиченому та нафтонасиченому об’єктах становить:

 

, (4.19)

 

отже,

 

. (4.20)

 

Визначення коефіцієнта об’ємного газонасичення

Більші можливості нейтронні методи мають при визначенні коефіцієнта газонасичення порового простору колекторів. Рішення цього завдання ґрунтується на наступних положеннях.

У газоносному колекторі у зв’язку з меншим вмістом водню, чим у водоносному заданої пористості, і зменшенням його густини коефіцієнт нейтронної пористості визначається співвідношенням:

 

, (4.21)

 

де H в і H г – водневі еквіваленти водонасиченої і газонасиченої часток порового простору при термобаричних умовах залягання колектора; D K п n d – виправлення за зменшення густини газоносного пласта в порівнянні з водоносним.

Зіставляючи формулу (4.21) з формулою, що визначає Kп n для повністю водонасиченого колектора, приходимо до висновку, що газонасичення породи викличе зменшення Kп n на величину

 

, (4.22)

 

відповідно:

 

, (4.23)

 

де nHг с – середнє число атомів водню, що припадають на молекулу газу; М г с середня відносна молекулярна маса газу; рпл – пластовий тиск; Тпл – температура пласта; d г с – середня густина газу при атмосферному тиску і t=0; d в p – густина води при пластовому тиску; Z – виправлення за стискування реальних газів.

Таким чином, знаючи Kп колектора і визначивши Kп n на його обводненій та газонасиченій ділянках, можна при відомих рпл, Тпл, Z, nHг с, М г с і DKп n за формулою (4.23) розрахувати коефіцієнт газонасичення.

Інший спосіб визначення коефіцієнта газонасичення колектора полягає у використанні залежностей DIn g=f(Kг) (Рис. 4.2), In= F (Kг) і Itn=F(Kг) при Kп= const або номограм, які складені за даними залежностями. Останні будуються для заданих свердловинних умов, типу і розміру зонда, на підставі вимірів In g, In й Itn на моделях газоносних колекторів або в свердловинах у колекторах, коефіцієнт газонасичення яких надійно встановлений.

При розрахунку подвійного різницевого параметра:

 

(4.24)

 

за D I /n g оп приймається інтенсивність нейтронного g -випромінювання у водонасиченому колекторі заданого мінерального складу, а за D I //n g оп – те ж, у колекторі із встановленим високим (опорним) коефіцієнтом газонасичення.

зонд НГМ-60

Рисунок 4.2 – Палетка залежностей D I n g= f (K г)

Встановлений за залежностями, які приведеним на рисунку 4.2, коефіцієнт Kг відповідає дійсному значенню, якщо пластовий тиск, температура і склад газу в досліджуваному й опорному колекторах ідентичні. В іншому випадку дійсне значення Kг розраховується за формулою:

 

, (4.25)

 

де індексом “оп” відзначені параметри, що відносяться до опорного колектора, а без цього індексу – до досліджуваного колектора.

При оцінці коефіцієнта газонасичення необхідно знати коефіцієнт пористості. У зв’язку із цим рішення завдання полегшується при використанні комплексу геофізичних методів і у першу чергу гамма-гамма-каротажу, наприклад, шляхом розв’язку системи рівнянь відносно Kп, Kг, Kв і Kн:

 

. (4.26)

 

На рисунку 4.3 приведена палетка фірми «Шлюмберже», що дозволяє за відомими коефіцієнтом нейтронної пористості Kп n (з виправленням за літологію), густиною d п породи і d ск скелетної фракції, параметром Pw вологості визначити Kп, Kг, Kв і Kн=1-(Kг+Kв).

Перевагою нейтронних методів оцінки коефіцієнта газонасичення колектора є можливість рішення завдання у свердловинах, які обсаджені сталевими колонами. Це дозволяє широко застосовувати нейтронні методи для спостереження за виробітком газових покладів.

Визначення коефіцієнтів залишкового нафтогазонасичення

Можливість визначення коефіцієнта нафтогазонасичення колектора за геофізичним даними дозволяє оцінити коефіцієнт залишкового нафтогазонасичення Kн пп (нафтогазонасичення промитої породи) і розрахувати коефіцієнт витіснення.

Коефіцієнт Kн пп залишкового нафтогазонасичення для заводнених ділянок покладу встановлюється за параметром Рнпп способами, які описані вище. Суттєвим є наступна обставина. На обводнених ділянках колектора коефіцієнт Kн п звичайно не перевершує 40–50 %. У цьому діапазоні значно знижується точність оцінки Kн. Зниження точності повинне компенсуватися більш прецизійним визначенням r п, що має в більшості випадків невисоке значення.

Рисунок 4.3 – Палетка для визначення коефіцієнтів K п, K г (K н) за d п, K п n і Р w (за даними фірми «Шлюмберже»)

У зв’язку із цим діаграми r у(r е) на ділянках нафтонасичених колекторів, які промиті пластовими водами, повинні реєструватися у великих масштабах. Реєстрацію діаграм мікрозондами та екранованим зондом для оцінки r пп із підвищеною точністю варто проводити багаторазово для набору необхідної статистики.

Найбільш точно параметр нафтонасичення встановлюється для ділянок колектора, які розташовані поблизу первинного водонафтового контакту, на яких:

 

. (4.27)

 

При використанні зондів оптимального розміру остання формула може бути замінена наступної:

 

, (4.28)

 

де r у пп і r у вп – оптимальні значення уявного опору, відповідно, в промитій і водоносній частинах колектора біля первинних водонафтового або водогазового контактів.

В іншому способі коефіцієнт залишкового нафтогазонасичення розраховується за формулою:

 

, (4.29)

 

де r нп і r пп, r в і r в пп, Пп п і П п пп – відповідно, електричні питомі опори колектора, його порових вод і поверхневих провідностей у первинному та промитому станах. В окремому випадку при витісненні нафти пластовими водами і r в= r в пп та Пп п= П п пп:

 

. (4.30)

 

За даними K н пп (отриманим на підставі вимірів Р н пп) складається карта змін K н пп по площі родовища та встановлюються залежності між K н пп та іншими параметрами колектора (глинистістю, проникністю, коефіцієнтом динамічної пористості, первинним нафтонасиченням та ін.). Такі карти, можуть використовуватись для визначення середніх значень β в і K н пп по площі родовища та для площ, що знову розбурюються, із колекторами, які характеризуються аналогічними петрофізичними властивостями.

За рубежем для оцінки коефіцієнта залишкового нафтонасичення використовуються залежності Тіксьє, Пупона і Лоу:

 

, (4.31)

 

де показник n змінюється від 2 до 5; а також лінійні зв’язки:

 

, (4.32)

. (4.33)

 

Згідно Дахнова встановлено зв’язок між K в з та K в зв:

 

. (4.34)

 

На рисунку 4.4 наведені залежності K в пп=f(K в), K н пп = 1- K в пп=F(K н), K н вил= K н- K н пп= F (K н) і коефіцієнта β в витіснення нафти від K в для різних значень показника а. Тут же для зіставлення наведені залежності K в пп=(K в)1/2 і K в пп=0,57× K в+0,43. З рисунку 4.4 випливає, що зі зменшенням а, тобто зі спрощенням конфігурації порових каналів, за інших рівних умов збільшуються K н вил, β в і K в пп, а з ростом K н підвищуються K н вил і β в та зменшується K в пп. Вибір параметра а повинен проводитися на підставі статистичної обробки даних вимірів K в пп (K н вил, β в) і K в(K н).

У газоносних колекторах коефіцієнт залишкового газонасичення може бути встановлений за даними вимірів I n, D I n g s I tn поблизу газоводяного контакту і наступній оцінці K г відповідно до залежностей цих параметрів від K г. Однак через невелику величину K г пп та суттєвого впливу свердловинних умов для розрахунку середнього значення K г пп із необхідною точністю потрібна велика кількість таких визначень.

1 – K в пп= f (K в) [ K н пп= f (K н)]; 2 – те ж, за Тіксьє K в пп=(K в)1/2; 3 – те ж, за Пупоном і Лоу K в пп=0,57× K в+0,43; 4 – K н вил= F (K н); 5 – b п= f (K в); 6 – β = f (K в) за Тіксьє

Рисунок 4.4 – Залежності K в пп = f ( K в ) [ K н пп = f ( K н )]; K н вил = Ф ( K н ) і β в= F ( K в ) для різних значень показника степені а в рівнянні

Визначення коефіцієнта витіснення нафти

Коефіцієнт витіснення нафти

 

(4.35)

 

розраховують за даними вимірів K н і K н пп або K в зв, якщо залежність bп = f (K в зв) або K в з= f (K в зв) визначена за результатами лабораторних досліджень керна.

Для багатьох теригенних колекторів

 

, (4.36)

 

де показник а змінюється від 0,3 до 0,7. У цих умовах

 

(4.37)

 

припускаючи 0,01/ K в близьке до 0,1.

Кількість нафти, що вилучається з одиниці об’єму колектора – об’ємний коефіцієнт вилучення нафтонасичення c вил і коефіцієнт витіснення β в, можна знайти за величиною параметра вологості. Вирішуючи обернену задачу, використовуючи залежність Рw=f(w), яка отримана за лабораторними або розрахунковими даними, а також номограму на рисунку 4.5, знаходимо

 

, (4.38)

 

 

Рисунок 4.5 – номограма для розрахунку параметра P w, а також визначення об’ємної вологості w та коефіцієнта K в(K н) неглинистих колекторів за параметром P w

, (4.39)

 

і, відповідно,

 

, (4.40)

. (4.41)

 

Для гідрофільних колекторів при w <0,2¸0,3 можна прийняти Pw=aв/ w q, де aв»aпaн і q»(m+n)/2.

У цих умовах

 

, (4.42)

 

і

 

. (4.43)

 

Точність визначення c н вил і β в зростає з підвищенням мінералізації промивної рідини.

Визначення сумарної потужності колектора

Під сумарною потужністю S h нафтогазонасиченого колектори розуміється різниця між висотними оцінками залягання покрівлі та підошви колектора в досліджуваній свердловині (при повнім насиченні колектора нафтою і газом) або різниця між висотними оцінками залягання покрівлі колектора і водонафтового та газонафтового контактів (при плаваючому покладі).

Покрівля і підошва колектора встановлюються відомими способами інтерпретації діаграм геофізичних методів дослідження свердловин. При цьому переваги віддаються діаграмам найбільш інформативних геофізичних параметрів. Для теригенного розрізу це діаграми r у (r е) (особливо мікрометодів), U ПС, I g, для карбонатного – діаграми I ng, I gg, D T.

Визначення ефективної потужності колектора

Під ефективною потужністю колектора розуміють сумарну потужність (по вертикалі) нафто- і газонасичених прошарків, що володіють динамічною пористістю і залягають у межах досліджуваної нафтогазонасиченой частини колектора.

Ефективна потужність встановлюється такими способами.

1. За кривими мікрозондів (ділянкам позитивних збільшень D r) і відносно низьким значенням r у і r е (на кривих екранованого мікрозонда) з вирахуванням потужності щільних, які виділяються різко підвищеними значеннями r у (r е), і глинистих прошарків, на ділянках яких D r у близьке до нуля.

При оцінці ефективної потужності тонкошаруватого колектора за даними мікрозондів і екранованого мікрозонда діаграми повинні записуватися в масштабі глибин 1:50 та у масштабах r у і r е, при яких мінімальні відхилення діаграм від нульової лінії перевищують 10 мм.

2. За кривими екранованого зонда, градієнт-зонда, якщо потужності ефективно пористих прошарків точно визначаються за даними цих зондів.

3. За кривими мікрокаверноміра і кіркоміра, що чітко встановлює зміни товщини глинистої кірки на ділянках свердловини, яка перетинає динамічно пористі прошарки. Мікрокавернограма повинна реєструватися в масштабі діаметру 1:1 і у масштабах глибин 1:200 або 1:50. Останній використається в тих випадках, коли ефективна потужність h еф прошарків визначається з комплексу кривих мікрокаверноміра, кіркоміра та кривих мікрозондів, які зареєстровані у масштабі глибин 1:50.

4. Методом анізотропії, якщо питомі електричні опори нафтонасичених і щільних прошарків значно відрізняються від опору глинистих. За даними ймовірного відношення r п гл/ r гл ш для заданого відношення r нп/ r гл за палеткою, яка приведена на рисунку 4.6, знаходять відношення (S h еф+S h пл)/S h, де S h еф і S h пл – сумарні потужності нафтонасичених і щільних прошарків; S h – сумарна потужність колектора.

Із суми S h еф+S h пл віднімають S h пл за діаграмами екранованого мікрозонда або знаходять середньостатистичне значення S h пл за законтурними свердловинами, де метод анізотропії дає величину S h пл/S h. Чим більше розходження між електричними питомими опорами r нп, r пл і r гл ш, а також між r пл, r вп і r гл ш, тим точніше встановлюється ефективна потужність. Метод рекомендується, в основному, для оцінки ефективної потужності піщано-глинистих колекторів.

5. За площею W ПС аномалії кривій U ПС Сумарна потужність нафтогазонасичених і щільних (неглинистих) прошарків дорівнює відношенню площі W ПС аномалії в межах нафтонасиченої частини колектора до величини статичного потенціалу:

 

, (4.44)

 

де поправочний коефіцієнт x ПС<1 знаходиться шляхом зіставлення ефективних потужностей, які отримані за даними методу U ПС та інших способів, або приблизно приймається рівним D U ¥ПС пл/D U ¥ПС (D U ¥ПС пл. – амплітуда аномалії D U ¥ПС у щільних прошарках великої потужності).

Рисунок 4.6 – Залежність r п гл/ r гл ш= f [ c п(c гл)]

Сумарну потужність щільних прошарків знаходять за даними мікрозондів та екранованого мікрозонда, потім розраховують:

 

. (4.45)

 

Цей спосіб рекомендується для теригенного розрізу при незначній різниці між електричними опорами колектора і порід, що вміщують (теригенний глинистий колектор), а також за умови насичення газом і нафтою всіх ефективних прошарків у межах всієї потужності колектора, що відмічається єдиною аномалією U ПС.

6. За площею W g аномалії кривої I g. Сумарна потужність неглинистих прошарків визначається відношенням W g кривої I g до амплітуди аномалії D I g= I g гл-1/ I g гл-0 в чистому неглинистому колекторі:

 

, (4.46)

 

де I g гл-1 і I g гл-0 – інтенсивності природного g -випромінювання в глині та колекторі, що не містить глинистих включень; x g=(I g гл-1- I g пл)/(I g гл-1- I g кол) (I g пл, I g кол – інтенсивності g -випромінювання, відповідно, в щільних породах і колекторі при d c= const).

Величину S h пл знаходять способами, які зазначені вище, і далі розраховують:

 

. (4.47)

 

Спосіб рекомендується для оцінки ефективної потужності піщано-глинистих колекторів у тих випадках, коли g -активність q g гл глинистих мінералів істотно перевищує g -активність мінералів q g ск, які входять до складу скелетної фракції.

Якщо перераховані способи визначення ефективної потужності використовуються при підрахунках запасів нафти і газу, то в розраховані значення S h еф повинна бути внесена поправка, яка враховує геометрію зустрічі пласта і свердловини. Якщо колектор має кут падіння a з азимутом q і розкривається свердловиною, яка викривлена під кутом d з азимутом j, шукана потужність S h еф розраховується за формулою:

 

, (4.48)

 

де S h еф д – ефективна потужність колектора, яка встановлена за діаграмами геофізичних параметрів у даній свердловині зазначеними способами, які зазначені вище.

 

 


Дата добавления: 2016-01-04; просмотров: 17; Мы поможем в написании вашей работы!

Поделиться с друзьями:






Мы поможем в написании ваших работ!