Решение для притока на неустановившемся режиме с постоянным забойным давлением

Последнее сообщение
AlNikS 851 15
Сен 12

Решение для притока на неустановившемся режиме с постоянным забойным давлением. Существует? Ткните носом?

 

Задача такая, на производстве при обосновании бурения мало кого интересует, какой дебит будет у скважин, когда воронка действительно разойдется и всё устаканится (этот дебит смутит инвесторов, т.к. окупаться скважина с ним будет долго). Интересуют запускные дебиты (не позднее месяца после запуска) и процент падения. Соответственно, есть некие производственные критерии, типа, скважина должна запуститься с дебитом по нефти Х и падать с процентом Y, тогда ОК. Если дебит меньше планки Х, то не ОК.

И вот есть у меня карты kh, есть фазухи, PVT, и я сижу и думаю, как бы мне построить карту запускных дебитов? Может, решение где-то на поверхности, я просто не вижу?

RomanK. 2139 16
Сен 12 #1

На фактические цифры ты маловероятно, что выйдешь. А что тебе это решение даст? Теоретически дебит в t=0 бесконечность.

AlNikS 851 15
Сен 12 #2

В t=0 мне хватит ума не брать дебит. Найду точки, где есть скважины, их запускной дебит, вычислю t при котором формула дает факт, буду это t использовать для расчета по всей площади, потом обрежу по пороговому дебиту и получу области куда выгодно бурить, а куда не выгодно.

RomanK. 2139 16
Сен 12 #3

Попробуй так, чтоли. В первом приближении
q(t)=q * EXP(-q/Q * t);
где q - стартовый дебит, Q - упругий запас
q=PI * (P - Pc)
PI-продуктивность, P-начальное пластовое; Pc - забойное давление
Q=ShB (P-Pc)
B - упругоемкость

RomanK. 2139 16
Сен 12 #4

Для неустановившегося режима, PI в начале будет выше, но я думаю иметь запас прочности это норм.

helgibh 68 12
Сен 12 #5

Глава 2 мануала, transient productivity index

sander 294 12
Сен 12 #6

по формуле Лейбензона попробуй, примерно как описал Роман, хотя есть ряд ограничений для ее применения но для первого приближения можно брать

kochichiro 924 16
Сен 12 #7

Если сапфир есть под рукой сделай test design с параметром q(p). Задай забойное и времени сколько нужно, потом на графике щелкаешь правой кнопкой dump data и снимешь точки для экселя или куда там тебе нужно. Формулировка уравнения на рисунке внизу

Stroncz 1116 17
Сен 12 #8

В уравнении Грингартена (1979) для неустановившегося режима если время устремить к нулю, то приток стремится к бесконечности, следовательно должно быть на малых временах ограничение по VLP. То есть при каком-то времени эти два решения сшиваются вместе. При каком именно - надо наверное в каждом конкретном случае смотреть (по аналогии со временем наступления псведоустойчивого режима). 

AlNikS 851 15
Сен 12 #9

RomanK. пишет:
Попробуй так, чтоли. В первом приближении q(t)=q * EXP(-q/Q * t); где q - стартовый дебит, Q - упругий запас q=PI * (P - Pc) PI-продуктивность, P-начальное пластовое; Pc - забойное давление Q=ShB (P-Pc) B - упругоемкость

Смущает сама концепция индекса продуктивности, он рассчитывается для установившегося режима.

AlNikS 851 15
Сен 12 #10

helgibh пишет:

Глава 2 мануала, transient productivity index

Какого мануала? :)

AlNikS 851 15
Сен 12 #11

sander пишет:

по формуле Лейбензона попробуй, примерно как описал Роман, хотя есть ряд ограничений для ее применения но для первого приближения можно брать

Погуглил, это ж вроде для газа?

AlNikS 851 15
Сен 12 #12

kochichiro пишет:

Это похоже на то, что нужно! 162.6 это для каких единиц?

RomanK. 2139 16
Сен 12 #13

Wasteland Rat пишет:

RomanK. пишет:
Попробуй так, чтоли. В первом приближении q(t)=q * EXP(-q/Q * t); где q - стартовый дебит, Q - упругий запас q=PI * (P - Pc) PI-продуктивность, P-начальное пластовое; Pc - забойное давление Q=ShB (P-Pc) B - упругоемкость

Смущает сама концепция индекса продуктивности, он рассчитывается для установившегося режима.

И это правильно, что сомневаешься. Эти формулы не подходят к твоему случаю с точки зрения точного описания.

Но с другой стороны, практически хорошо должны (я уверен и могут) описать вывод скважины на режим и решают твою задачу, исходя из исходных данных.

Я бы упростил более:

1. q (начальный, первые сутки) т/сут = (kh/mu) * (P(пластовое, ориентировочное) - P(забойное))

2. Построить корреляцию между расчетным/фактическим дебитом

3. q (последующие сутки) т/сут = q(начальный) * EXP(-I (интенсивность падения давления) * t)

4. Корреляция между интенсивностью падения давления и параметрами пласта

 

Подход идеальный.

valer 441 11
Сен 12 #14

Фактический запускной дебит сильно зависит от технологии бурения. Одну разведочную пробурили на инвертном растворе да ещё получили открытый фонтан. Дебит запускной 360 м3/сут. Средний по залежи на тот период 10 м3/сут.

helgibh 68 12
Сен 12 #15

Того самого, который изучали в HW, well test. Формула от kochichiro очень похожа на уравнение линейного стока, но коэффициентами может отличаться.

RomanK. 2139 16
Сен 12 #16

От технологии бурения можно получить меньший дебит, никак не больший :) Хорошее пластовое и продуктивность. Бывает.

AlNikS 851 15
Сен 12 #17

valer пишет:

Фактический запускной дебит сильно зависит от технологии бурения. Одну разведочную пробурили на инвертном растворе да ещё получили открытый фонтан. Дебит запускной 360 м3/сут. Средний по залежи на тот период 10 м3/сут.

В данный момент меня интересует потенциал пласта, а не качество будущего вскрытия. А когда разница в 30 раз, то это уже геологические неопределенности а не качество вскрытия, скорее всего.

AlNikS 851 15
Сен 12 #18

RomanK. пишет:

Wasteland Rat пишет:

RomanK. пишет:
Попробуй так, чтоли. В первом приближении q(t)=q * EXP(-q/Q * t); где q - стартовый дебит, Q - упругий запас q=PI * (P - Pc) PI-продуктивность, P-начальное пластовое; Pc - забойное давление Q=ShB (P-Pc) B - упругоемкость

Смущает сама концепция индекса продуктивности, он рассчитывается для установившегося режима.

И это правильно, что сомневаешься. Эти формулы не подходят к твоему случаю с точки зрения точного описания.

Но с другой стороны, практически хорошо должны (я уверен и могут) описать вывод скважины на режим и решают твою задачу, исходя из исходных данных.

Я бы упростил более:

1. q (начальный, первые сутки) т/сут = (kh/mu) * (P(пластовое, ориентировочное) - P(забойное))

2. Построить корреляцию между расчетным/фактическим дебитом

3. q (последующие сутки) т/сут = q(начальный) * EXP(-I (интенсивность падения давления) * t)

4. Корреляция между интенсивностью падения давления и параметрами пласта

 

Подход идеальный.

Тогда это можно еще более упростить:

1. Построить статистическую корреляцию между kh*(Pпласт-Pзаб) и входным дебитом существующих скважин.

2. Пересчитать карту потенциальных дебитов из карты kh*(Pпласт-Pзаб) по корреляции.

:)

Это не совсем то, что я хотел, но тоже решение.

Что-то все сложнее, чем я думал по началу, я думал формулы готовые есть по этому поводу :)

AlNikS 851 15
Сен 12 #19

helgibh пишет:

Того самого, который изучали в HW, well test. Формула от kochichiro очень похожа на уравнение линейного стока, но коэффициентами может отличаться.

Я по-моему вполне четко сказал, что мне нужно решение для постоянного забойного давления :)

helgibh 68 12
Сен 12 #20

Да, я заметил. Как раз в рассуждениях про transient productivity index об этом речь идет. И kochichiro привел формулу для Pwf=const

valer 441 11
Сен 12 #21

а не качество будущего вскрытияНет, я не о будущем. Вы же хотите использовать историю по существующему фонду.

RomanK. 2139 16
Сен 12 #22

Wasteland Rat пишет:

RomanK. пишет:

Wasteland Rat пишет:

RomanK. пишет:
Попробуй так, чтоли. В первом приближении q(t)=q * EXP(-q/Q * t); где q - стартовый дебит, Q - упругий запас q=PI * (P - Pc) PI-продуктивность, P-начальное пластовое; Pc - забойное давление Q=ShB (P-Pc) B - упругоемкость

Смущает сама концепция индекса продуктивности, он рассчитывается для установившегося режима.

И это правильно, что сомневаешься. Эти формулы не подходят к твоему случаю с точки зрения точного описания.

Но с другой стороны, практически хорошо должны (я уверен и могут) описать вывод скважины на режим и решают твою задачу, исходя из исходных данных.

Я бы упростил более:

1. q (начальный, первые сутки) т/сут = (kh/mu) * (P(пластовое, ориентировочное) - P(забойное))

2. Построить корреляцию между расчетным/фактическим дебитом

3. q (последующие сутки) т/сут = q(начальный) * EXP(-I (интенсивность падения давления) * t)

4. Корреляция между интенсивностью падения давления и параметрами пласта

 

Подход идеальный.

Тогда это можно еще более упростить:

1. Построить статистическую корреляцию между kh*(Pпласт-Pзаб) и входным дебитом существующих скважин.

2. Пересчитать карту потенциальных дебитов из карты kh*(Pпласт-Pзаб) по корреляции.

:)

Это не совсем то, что я хотел, но тоже решение.

Что-то все сложнее, чем я думал по началу, я думал формулы готовые есть по этому поводу :)

Это именно то, что тебе нужно. Можешь сравнить с формулой, приводимой коширо и ты заметишь, что они очень похожи. Сначала определись со стартовым дебитом, потом определись с динамикой падения. Экспоненциальное падение - это случай постоянного забойного давления.

sander 294 12
Сен 12 #23

Wasteland Rat пишет:

Погуглил, это ж вроде для газа?

возможно ты нашел для газа, там и для нефти есть...

AlNikS 851 15
Сен 12 #24

RomanK. пишет:

Это именно то, что тебе нужно. Можешь сравнить с формулой, приводимой коширо и ты заметишь, что они очень похожи. Сначала определись со стартовым дебитом, потом определись с динамикой падения. Экспоненциальное падение - это случай постоянного забойного давления.

Никакой корреляции не наблюдается, чтоб её. Видимо, технологический фактор сильнее геологического. Либо карты все неправильные :)

Dorzhi 970 17
Сен 12 #25

 Petroleum production systems, Michael Economides, страница 17, уравнение 2-7. решение притока для неустановившегося режима с постоянным забойным давлением - аналитическое решение уравнения диффузии, Earlougher (1977), точно такое же какое показал уже konichiro , коэффициенты для field units, время в часах.

Krichevsky 731 14
Сен 12 #26

>> В данный момент меня интересует потенциал пласта, а не качество будущего вскрытия.

Вот эту фразу не понял. Вы хотите посчитать запускной дебит и темп его падения, не задаваясь несовершенством скважины?

valer 441 11
Сен 12 #27

Это не технологический, это человеческий фактор!

Кто-же покажет дебит из бурения ниже проектного? 

RomanK. 2139 16
Сен 12 #29

Да бог с ним с теорией. Результат какой?

sander 294 12
Сен 12 #30

так это для вертикально-трещиноватых?

AlNikS 851 15
Сен 12 #31

RomanK. пишет:

Да бог с ним с теорией. Результат какой?

Результат узнаем через год где-то. :)

AlNikS 851 15
Сен 12 #32

sander пишет:

так это для вертикально-трещиноватых?

для вертикальной скважины с ГРП

RomanK. 2139 16
Сен 12 #33

А как же опробация по имеющимся данным? Или вот это "никакой корреляции не наблюдается" и был окончательным ответом?

AlNikS 851 15
Сен 12 #34

История добычи восстановлена каким-то макаром из накопленной добычи 90-х годов. По забойным давлениям вообще труба. Скважины работали фонтаном. Так что да :)

AlNikS 851 15
Сен 12 #35

Если кому-нибудь это что-нибудь скажет, речь идет о восстановлении упоротой в хлам Талинской площади. Ну чтоб ни у кого вопросов по верификации на исторические данные больше не возникало :)

RomanK. 2139 16
Сен 12 #36

Так ты участник операции ТНК-BP, по восстановлению разработки Талинской площади, там где намечены порядка 15 участков с разными видами воздействий и бубнами под предводительством Бриллианта?

voron4m 384 14
Сен 12 #37
Прошу прощения, но я не рекламирую ПО...
Контора не бедная, Ecrin скорее всего есть. Kochichiro уже предложил воспользоваться Saphir-ом.
Но в Ecrin есть достаточно мощный пакет по анализу и прогнозу добычи - Topaze.
Довольно просто построить модель скважины (любой), граничные условия, начальное и рабочее давления. Можешь построить различные варианты по скину, проницаемости, длина трещины или её проводимости - на любой изощрённый вкус. Хочешь поиграться фазовыми - есть такая опция есть в NL версии. Есть даже немного примитивный расчёт VLP. Великолепно рассчитывает дебиты на всех стадиях: начальный, установившийся и прочее. Нужны деклайны или PI - есть дополнительные графики.
Если надо по нескольким скважинам - пожалуйста.
Надо по месторождению - пожалуйста.
Неудобно сделан экспорт данных, но где это сделано хорошо? Можно Copy-Paste или dump data. Можно поковырять Word report в Program Files\Kappa\Ecrin и подписать интерфейс. К Saphir-у прикручено, работает без проблем.
voron4m 384 14
Сен 12 #38

В тему (по горизонталкам правда).

http://www.petroleumengineers.ru/node/5810

AlNikS 851 15
Сен 12 #39

voron4m пишет:

Контора не бедная, Ecrin скорее всего есть.

Был да сплыл... Я не в проектном институте.

kochichiro 924 16
Сен 12 #40

Еще лучше вообще тогда воспользоваться Rubis там тебе и граничные условия будут и можно как посоветовал Stroncz задать ограничение по пропускной способности лифта. А сетка Вороного позволит не заморачиваться с LGR для трещины гидроразрыва (поскольку генерируется автоматом) и сделать мелкие временные шаги вплоть до 0,1 часа и наблюдать изменение дебита при постоянном давлении вплоть от пускового неустановившегося режима до псевдо-установившегося с влиянием границ.

Короче первоначальная идея эволюционировала и пришла в итоге к гидродинамическому моделированиюSmile

 

voron4m 384 14
Сен 12 #41

Я тоже не в проектном институте работаю, в простом отделе по разработке. Wink

Все прогнозы по новым скважинам, переводам, фракам или прочим работам делаем как раз в Topaze.

RomanK. 2139 16
Сен 12 #42

Запускные дебиты сильно с простым Дюпюи и Топазом сходятся?

voron4m 384 14
Сен 12 #43

Вертикальных скважин нет. Горизонтальные раньше считали по Joshi. Расчётный дебит (начальный) из Topaze в 1,5 -2 раза выше, чем по Joshi. Месяца через 2-3 ( иногда до полугода в зависимости от проницаемости ) практически один в один. Сходимость проверена на практике по большинству новых скважин.Давай живой пример вертикальной скважины - прогоню в Topaze, но только в понедельник получится.

RomanK. 2139 16
Сен 12 #44

Завышение в первые месяцы это по факту?

voron4m 384 14
Окт 12 #45

RomanK. пишет:

Завышение в первые месяцы это по факту?

Да, по факту. Вот "живой" пример, дебиты с тест-сепаратора.

 

RomanK. 2139 16
Окт 12 #46

Расчётного значения не хватает (?)

AlNikS 851 15
Май 13 #47

RomanK. пишет:

Так ты участник операции ТНК-BP, по восстановлению разработки Талинской площади, там где намечены порядка 15 участков с разными видами воздействий и бубнами под предводительством Бриллианта?

Скажем так, я "участник" не со стороны Бриллианта :)
Сижу до сих пор ковыряю материалы и думаю, что же тут можно сделать. Возможно, было бы рентабельно вложиться в изобретение машины времени и отмотать обратно к началу разработки...

RomanK. 2139 16
Сен 13 #48

Известно следующее простое решение для пуска скважины с постоянным забойным давлением (точность 1%)

Поставлять в формулу Дарси за место радиуса контура

RomanK. 2139 16
Сен 13 #50

точность 99%

Страницы

Go to top