RU2165520C2 - Способ моделирования влияния взаимодействия скважин на водную фракцию, производимую подземным месторождением углеводородов (варианты) - Google Patents
Способ моделирования влияния взаимодействия скважин на водную фракцию, производимую подземным месторождением углеводородов (варианты) Download PDFInfo
- Publication number
- RU2165520C2 RU2165520C2 RU96123972/03A RU96123972A RU2165520C2 RU 2165520 C2 RU2165520 C2 RU 2165520C2 RU 96123972/03 A RU96123972/03 A RU 96123972/03A RU 96123972 A RU96123972 A RU 96123972A RU 2165520 C2 RU2165520 C2 RU 2165520C2
- Authority
- RU
- Russia
- Prior art keywords
- wells
- production
- well
- significant data
- water fraction
- Prior art date
Links
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/20—Displacing by water
Landscapes
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
- Automatic Analysis And Handling Materials Therefor (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Levels Of Liquids Or Fluent Solid Materials (AREA)
- Physical Or Chemical Processes And Apparatus (AREA)
Abstract
Изобретение относится к нефтедобывающей промышленности. Задачей изобретения является оптимизация добычи углеводородов за счет определения взаимодействия между различными скважинами, влияющих на величину водной фракции в вытекающих жидкостях, производимых одной или несколькими эксплуатационными скважинами из некоторой зоны подземного месторождения жидких углеводородов в процессе эксплуатации этого месторождения, подвергающегося промывке при помощи жидкости под давлением, нагнетаемой через одну или несколько нагнетательных скважин, или промывке водой из некоторой близлежащей водоносной зоны. Для этого осуществляют отбор совокупности значащих данных, исходя из измерений, полученных из истории нагнетания промывочной жидкости в данное месторождение и из истории добычи жидкого продукта, вытекающего из совокупности эксплуатационных скважин данной нефтеносной зоны. Затем итеративно устанавливают некоторую линейную оптимизированную модель, связывающую изменения во времени величин этих водных фракций с изменениями во времени упомянутых значащих данных, относящихся к параметрам других скважин, входящих в упомянутую систему скважин. 2 с. и 8 з.п. ф-лы, 14 ил.
Description
Настоящее изобретение касается способа моделирования влияний взаимодействия между различными скважинами на водную фракцию в вытекающих жидкостях, производимых подземным месторождением углеводородов в процессе его эксплуатации, подвергающимся промывке жидкостью под давлением, с целью оптимизации производительности данного месторождения.
Добыча воды представляет собой основную проблему при эксплуатации нефтяных месторождений. Бывает, что эксплуатанты нефтяного месторождения оказываются перед лицом ситуации, когда водяная фракция в дебите той или иной скважины является весьма значительной, тогда как степень рекуперации нефти на месте остается низкой, что наглядно демонстрирует неэффективность осуществляемой промывки. Эти обстоятельства могут привести к прекращению добычи на данной скважине со всеми вытекающими отсюда экономическими последствиями в результате не лучших решений по подавлению этих появлений воды. Именно в рамках эксплуатации стратифицированных источников, подвергающихся промывке водой, например, отмечаются иногда сложные изменения этой водной фракции.
Известен способ локальной обработки скважины, где происходит появление воды, путем закупоривания зон скважины, откуда поступает вода, при помощи закачивания цементного раствора, соответствующих полимеров, гелей и т.п. с использованием специального оборудования для разграничения обрабатываемых зон. Эта технология достаточно трудна в использовании, поскольку необходимо сначала хорошо разграничить критические зоны. При этом выполняемые операции являются тяжелыми и дорогостоящими и экономически не оправдываются на скважинах, которые часто эксплуатируются на пределе рентабельности.
Для ограничения слишком больших поступлений воды известно также применение к резервуарам глобальной обработки путем закачивания в них, например, полимеров. Однако эффективность этой технологии остается достаточно низкой и особенно трудно прогнозируемой.
Физика источника или месторождения обычно является весьма сложной. Рассматривается случай скважины, пересекающей некоторое число i уровней источников, рассматриваемых как гидравлически независимые в масштабе ближайшего окружения данной скважины (i=2 в случае фиг. 1). Под действием откачивания, например, с некоторым дебитом Q, определяемым производительностью насоса, давление в донной части скважины стабилизируется на уровне Pwf (динамическое давление). При этом функционирование скважины выражается следующими соотношениями:
Qi = IPi (Pi-Pwf),
где Pi представляет собой давление в слое i. Общий дебит скважины Q состоит из суммы вкладов Qi всех слоев i, причем каждый такой вклад зависит от показателя производительности IPi рассматриваемого слоя и от разности приложенных давлений Pi - Pwf. Водная фракция скважины fw представляет собой среднее значение водных фракций fwi каждого слоя, уравновешенное его вкладом в общий дебит скважины.
Qi = IPi (Pi-Pwf),
где Pi представляет собой давление в слое i. Общий дебит скважины Q состоит из суммы вкладов Qi всех слоев i, причем каждый такой вклад зависит от показателя производительности IPi рассматриваемого слоя и от разности приложенных давлений Pi - Pwf. Водная фракция скважины fw представляет собой среднее значение водных фракций fwi каждого слоя, уравновешенное его вкладом в общий дебит скважины.
Выражение Qi = IPi (Pi-Pwf) показывает, что любое изменение ÆPi давления Pi в некотором слое выражается в изменении ÆQi дебита Qi этого слоя, и, если водные фракции слоев различны, в изменении водной фракции скважины в зависимости от модификаций вкладов, относящихся к каждому слою, в общую производительность скважины. Изменение ÆPi давления в некотором слое может быть вызвано, в частности, изменением расходов закачивания или дебитов соседних скважин. В то же время в том случае, когда давления в различных слоях отличаются друг от друга, изменение давления добычи P выражается в изменении распределения дебитов (di).
Между тем, в том случае, когда давления Pi в различных слоях существенно отличаются друг от друга, любое изменение воздействия на скважину или ее нагрузки: расхода Q насоса или давления Pwf в скважине выражается изменением водной фракции, ее увеличением или уменьшением в соответствии с относительным распределением насыщенностей и давлений в каждом слое.
И хотя физические характеристики месторождения будут весьма сложными и давления часто будут из-за недостатка точных измерений переменными с неизвестными законами этих изменений, метод в соответствии с предлагаемым изобретением позволяет, тем не менее, в системе скважин, пронизывающих зону подземного месторождения углеводородов, подвергающегося в процессе эксплуатации промывке жидкостью под давлением (независимо от того, идет ли речь о закачиваемой под землю жидкости или о жидкости, поступающей из соседней водоносной зоны), провести удовлетворительное моделирование влияний взаимодействия между скважинами на водную фракцию в вытекающей жидкости, производимой по меньшей мере одной из этих скважин, с целью оптимизации добычи углеводородов на данном месторождении.
Способ в соответствии с предлагаемым изобретением отличается тем, что он включает следующее:
отбор совокупности значащих данных исходя из необработанных данных, извлеченных из истории закачивания промывочной жидкости в данное месторождение и истории добычи вытекающей жидкости из одной или нескольких эксплуатационных скважин;
итеративное установление линейной оптимизированной модели, связывающей изменения во времени значащих данных, относящихся к водной фракции в добыче упомянутой эксплуатационной скважины, с изменениями во времени значащих данных, относящихся к другим скважинам данной системы скважин.
отбор совокупности значащих данных исходя из необработанных данных, извлеченных из истории закачивания промывочной жидкости в данное месторождение и истории добычи вытекающей жидкости из одной или нескольких эксплуатационных скважин;
итеративное установление линейной оптимизированной модели, связывающей изменения во времени значащих данных, относящихся к водной фракции в добыче упомянутой эксплуатационной скважины, с изменениями во времени значащих данных, относящихся к другим скважинам данной системы скважин.
Факторы взаимодействия скважин, затрагивающие производство воды, наглядно выявляются реализованной таким образом моделью, поэтому инженеры данного месторождения оказываются в состоянии с лучшими шансами на успех воздействовать на различные параметры: выбор скважин закачивания жидкости под давлением, расходы закачивания, дебиты добычи и т.д. с тем, чтобы повысить эффективность промывания и увеличить степень извлечения нефти.
В соответствии с изобретением отбор значащих данных содержит частотную фильтрацию изменений необработанных данных, относящихся, например, к водной фракции этой эксплуатационной скважины, с одной стороны, и к водным фракциям других скважин данной системы скважин, с другой стороны.
В соответствии с вариантом выполнения предлагаемого изобретения отбор значащих данных содержит, например, выявление флюктуаций на низкой частоте, которая значительно ниже частотного диапазона, в котором были измерены или получены необработанные данные, затрагивающие водную фракцию.
В соответствии с вариантом выполнения предлагаемого изобретения сбор значащих данных содержит отбор среди производительных скважин и/или скважин закачивания некоторого сокращенного числа скважин, характеризующихся наиболее сильным взаимодействием с упомянутой производительной скважиной.
Отбор значащих данных может содержать, например, предварительную статистическую обработку необработанных данных и в случае необходимости отбор среди них некоторой совокупности данных, представляющих равномерное распределение во времени.
В соответствии со способом использования предлагаемого изобретения данный метод включает приложение к одной или к нескольким скважинам закачивания или добычи произвольных воздействий, изменяющих входные необработанные данные таким образом, чтобы лучше отобрать те скважины, которые представляют взаимодействия.
В соответствии со способом использования предлагаемого изобретения подходящим для случая, когда осуществляется моделирование эффектов взаимодействий, оказываемых различными скважинами данной системы скважин на водные фракции в вытекающих жидкостях, производимых соответственно различными эксплуатационными скважинами, подвергающимися промывке жидкостью под давлением, с целью оптимизации производительности данного месторождения, предпочтительно осуществляется дополнительная общая оптимизация различных построенных моделей, принимая во внимание перекрестные взаимодействия между значащими данными, эффективно действующими в каждой из этих моделей, с тем, чтобы добиться максимальной общей производительности данной нефтеносной зоны.
Реализация модели, точно предсказывающей поведение скважин, в результате осуществления метода в соответствии с предлагаемым изобретением позволяет надежно оценить эффективность обработки скважин и сделать это лучше, чем удается при помощи существующих на сегодняшний день методов, реализуемых исходя из некоторого среднего более или менее репрезентативного поведения или характеристик этих скважин. Такая модель, распространенная в масштабе некоторой системы скважин, представляет собой инструмент оптимизации добычи нефти на данном месторождении.
Моделирование, осуществленное в соответствии с предлагаемым методом, дает следующий эффект:
уточнение картины данного месторождения, поскольку качественная интерпретация выявленных взаимовлияний позволяет уточнить гидравлические связи между скважинами и корреляции пластов;
усовершенствование диагностики состояния промывки месторождения, поскольку изменения водной фракции непосредственно связаны с контрастами насыщения между различными слоями, то есть с состоянием их промывки. Анализ взаимных влияний позволяет лучше выполнить отбор скважин-кандидатов для превентивной обработки с целью сдерживания прихода воды и даже усовершенствовать оперативные условия такой обработки. В то же время, исходя из корреляционных связей между скважинами и сопоставления поведения нескольких эксплуатационных скважин, можно получить информацию, касающуюся состояния промывки относительно единицы поверхности каждого слоя данного месторождения.
уточнение картины данного месторождения, поскольку качественная интерпретация выявленных взаимовлияний позволяет уточнить гидравлические связи между скважинами и корреляции пластов;
усовершенствование диагностики состояния промывки месторождения, поскольку изменения водной фракции непосредственно связаны с контрастами насыщения между различными слоями, то есть с состоянием их промывки. Анализ взаимных влияний позволяет лучше выполнить отбор скважин-кандидатов для превентивной обработки с целью сдерживания прихода воды и даже усовершенствовать оперативные условия такой обработки. В то же время, исходя из корреляционных связей между скважинами и сопоставления поведения нескольких эксплуатационных скважин, можно получить информацию, касающуюся состояния промывки относительно единицы поверхности каждого слоя данного месторождения.
Другие характеристики и преимущества способа в соответствии с предлагаемым изобретением будут показаны в приведенном ниже описании не являющихся ограничительными примеров реализации этого изобретения со ссылками на приведенные чертежи, на которых:
фиг. 1 представляет собой схематический вид эксплуатационной скважины, производящей вытекающую жидкость с двух уровней месторождения, рассматриваемых в качестве гидравлически независимых друг от друга в масштабе некоторой зоны, окружающей данную эксплуатационную скважину;
фиг. 2 представляет собой схематическую иллюстрацию связи, которая существует между возмущениями, затрагивающими расход закачивания промывочной жидкости и/или добычи вытекающей жидкости для соседних скважин;
фиг. 3 представляет собой схематическую иллюстрацию способа зависимости, устанавливаемой выбранной линейной моделью;
фиг. 4 представляет собой возможную схему размещения друг относительно друга рассматриваемых скважин, на практических данных которых предлагаемый метод был опробован;
фиг. 5 схематически представляет изменения в функции времени необработанных измерений водной фракции скважины;
фиг. 6 схематически представляет изменения в функции времени среднемесячных значений водной фракции скважины;
фиг. 7 схематически представляет частотный спектр средних величин водной фракции скважины;
фиг. 8 схематически представляет изменения в функции времени среднемесячных значений водной фракции скважины, показанные в виде пунктирной кривой, скорректированные в результате фильтрации высоких частот спектра, показанного на фиг. 7, и представляющие собой выходные данные;
фиг. 9 схематически представляет частотный спектр величин месячного дебита эксплуатационной скважины, задействованного в данной модели;
фиг. 10 схематически представляет изменения в функции времени величин месячного дебита, производимого скважиной, скорректированных в результате фильтрации высоких частот спектра, показанного на фиг. 9;
фиг. 11 схематически представляет частотный спектр средних величин месячных объемов воды, закачиваемой в скважину;
фиг. 12 схематически представляет изменения в функции времени среднемесячных расходов скважины закачивания, скорректированные в результате фильтрации высоких частот спектра, показанного на фиг. 11;
фиг. 13 схематически представляет примеры функций взаимной корреляции между водной фракцией скважины и соответственно месячными дебитами производства скважин;
фиг. 14 схематически представляет сопоставление результатов модели, полученных для скважины с реальными измерениями, выполненными на этой скважине.
фиг. 1 представляет собой схематический вид эксплуатационной скважины, производящей вытекающую жидкость с двух уровней месторождения, рассматриваемых в качестве гидравлически независимых друг от друга в масштабе некоторой зоны, окружающей данную эксплуатационную скважину;
фиг. 2 представляет собой схематическую иллюстрацию связи, которая существует между возмущениями, затрагивающими расход закачивания промывочной жидкости и/или добычи вытекающей жидкости для соседних скважин;
фиг. 3 представляет собой схематическую иллюстрацию способа зависимости, устанавливаемой выбранной линейной моделью;
фиг. 4 представляет собой возможную схему размещения друг относительно друга рассматриваемых скважин, на практических данных которых предлагаемый метод был опробован;
фиг. 5 схематически представляет изменения в функции времени необработанных измерений водной фракции скважины;
фиг. 6 схематически представляет изменения в функции времени среднемесячных значений водной фракции скважины;
фиг. 7 схематически представляет частотный спектр средних величин водной фракции скважины;
фиг. 8 схематически представляет изменения в функции времени среднемесячных значений водной фракции скважины, показанные в виде пунктирной кривой, скорректированные в результате фильтрации высоких частот спектра, показанного на фиг. 7, и представляющие собой выходные данные;
фиг. 9 схематически представляет частотный спектр величин месячного дебита эксплуатационной скважины, задействованного в данной модели;
фиг. 10 схематически представляет изменения в функции времени величин месячного дебита, производимого скважиной, скорректированных в результате фильтрации высоких частот спектра, показанного на фиг. 9;
фиг. 11 схематически представляет частотный спектр средних величин месячных объемов воды, закачиваемой в скважину;
фиг. 12 схематически представляет изменения в функции времени среднемесячных расходов скважины закачивания, скорректированные в результате фильтрации высоких частот спектра, показанного на фиг. 11;
фиг. 13 схематически представляет примеры функций взаимной корреляции между водной фракцией скважины и соответственно месячными дебитами производства скважин;
фиг. 14 схематически представляет сопоставление результатов модели, полученных для скважины с реальными измерениями, выполненными на этой скважине.
Водная фракция скважины увеличивается с течением времени даже в том случае, если расходы закачиваемой воды и дебиты эксплуатационных скважин остаются постоянными. Это - дрейф, возникающий вследствие постоянной промывки слоев водой под давлением, а также постепенного замещения нефти водой в данном месторождении. В данном случае речь идет о достаточно медленно развивающемся явлении, проявляющемся с момента проникновения воды в эксплуатационные скважины и растягивающемся на несколько лет. Таким образом можно считать, что водная фракция той или иной скважины состоит из некоторого дрейфа и флюктуаций, возникающих вследствие возмущений на соседних скважинах:
fw = дрейф = Æfw (флюктуации)
Определение измененной водной фракции fw скважины получается при этом в результате учета упомянутого дрейфа, связанного с суммарным производством жидкости на этой скважине, и путем моделирования связи, которая существует между возмущениями, возникающими вследствие изменений расхода закачиваемой в недра жидкости и/или изменений дебита соседних скважин в соответствии со схемой, представленной на фиг. 2.
fw = дрейф = Æfw (флюктуации)
Определение измененной водной фракции fw скважины получается при этом в результате учета упомянутого дрейфа, связанного с суммарным производством жидкости на этой скважине, и путем моделирования связи, которая существует между возмущениями, возникающими вследствие изменений расхода закачиваемой в недра жидкости и/или изменений дебита соседних скважин в соответствии со схемой, представленной на фиг. 2.
Как уже было сказано, метод в соответствии с предлагаемым изобретением содержит определение некоторой линейной системы, которая связывает изменения водной фракции скважины с изменениями закачивания промывочной жидкости и добычи на соседних скважинах. Выбирается, например, саморегрессивная модель типа ARX, выбранная в библиотеке программного математического обеспечения типа MATLAB®, хорошо известного специалистам, которая позволяет установить передаточную функцию, которая может существовать между двумя сигналами. Эта передаточная функция характеризует соответствующую физическую систему.
Линейная модель ARX, связывающая некоторый входной сигнал x с некоторым сигналом y, как это схематически представлено на фиг. 3, характеризуется следующим уравнением:
A(g) y(t) = B(g) x(t-nk) + e(t),
при nk: задержка или запаздывание;
g: оператор запаздывания
A(g) = 1 + a1g-1 +...+ anag-na, на порядок A(g)
B(g) = b1 + b2g-1 +...+ bnbg-nb, nb порядок B(g)
Говоря более определенно:
y(t)+a1y(t-l)+. . . +anay(t-na)=b1x(t-nk)+ +b2x(t-nk-1) +...+ bnbx(t-nk-nb+1)+e(t)
В том случае, если na = 0, данная модель является трансверсной: выход зависит только от входов.
A(g) y(t) = B(g) x(t-nk) + e(t),
при nk: задержка или запаздывание;
g: оператор запаздывания
A(g) = 1 + a1g-1 +...+ anag-na, на порядок A(g)
B(g) = b1 + b2g-1 +...+ bnbg-nb, nb порядок B(g)
Говоря более определенно:
y(t)+a1y(t-l)+. . . +anay(t-na)=b1x(t-nk)+ +b2x(t-nk-1) +...+ bnbx(t-nk-nb+1)+e(t)
В том случае, если na = 0, данная модель является трансверсной: выход зависит только от входов.
В том случае, если na≠0, данная модель является рекурсивной: выход зависит не только от входов, но также и от предыдущих выходов.
Для упрощения записи определяется линейная модель с одним единственным входом x. Очевидно, однако, что такая модель легко может быть обобщена на множество входов.
Можно убедиться в том, что упомянутый выбор линейной модели был вполне легитимным, рассчитывая для этого индивидуальные изменения водной фракции некоторой скважины, соответствующие n различным возмущениям, и убеждаясь в том, что общее изменение водной фракции, являющееся результатом действия всех n возмущений, представленных одновременно, оказываются равными сумме рассчитанных индивидуальных изменений с точностью до влияний дрейфа.
Отбор значащих данных
Для моделирования взаимодействий, существующих между n скважинами нагнетания промывочной жидкости или производства продукта W1, W2, ... Wn, используются необработанные данные эксплуатации, полученные из истории добычи и нагнетания, и на основе этих данных формируются значащие данные.
Для моделирования взаимодействий, существующих между n скважинами нагнетания промывочной жидкости или производства продукта W1, W2, ... Wn, используются необработанные данные эксплуатации, полученные из истории добычи и нагнетания, и на основе этих данных формируются значащие данные.
Исторические справки добычи образованы измеренными данными: измерения нагнетаемых и добываемых расходов, измерения водной фракции и т.п. с более или менее регулярной периодичностью отбора образцов. Эти измерения часто бывают "зашумлены" и представляют обычно большую дисперсию. Таким образом, прежде всего важно сделать их значащими в большей степени следующим путем:
исключение аберрантных или отклоняющихся от некоторой нормы измерений, связанных с шумовыми эффектами, и устранение частей наиболее высокой частоты спектра изменений необработанных измерений, причем все это делается, в частности, при помощи статистических методов или методов обработки сигналов, известных в данной области техники;
переоценка в случае необходимости, исходя из необработанных измерений, полученных на слишком нерегулярных интервалах, набора данных с постоянным шагом квантования.
исключение аберрантных или отклоняющихся от некоторой нормы измерений, связанных с шумовыми эффектами, и устранение частей наиболее высокой частоты спектра изменений необработанных измерений, причем все это делается, в частности, при помощи статистических методов или методов обработки сигналов, известных в данной области техники;
переоценка в случае необходимости, исходя из необработанных измерений, полученных на слишком нерегулярных интервалах, набора данных с постоянным шагом квантования.
Выполняется также отбор скважин, данные по которым будут приняты во внимание, среди эксплуатационного поля W2, ..., Wn, в наибольшей степени способных к взаимодействию с данными по скважине W1, для которой в данном случае моделируется водная фракция. Для достижения этой цели осуществляется взаимная корреляция для каждой из пар скважин W1, W2, ..., W1, Wn, значащих данных, полученных ранее, и водной фракции скважины W1, и среди скважин W2, ..., Wn оставляются те, для которых коэффициент взаимной корреляции является наиболее высоким.
Выбирая значащие данные скважин, в наибольшей степени способных к взаимодействию, их применяют в качестве входных данных для выбранной линейной модели и определяют специальное уравнение, моделирующее взаимодействия между отобранными скважинами. Осуществляя затем анализ и интерпретацию полученных результатов репрезентативной модели, можно воздействовать на факторы, способные уменьшить водную фракцию моделируемых скважин, и увеличить таким образом производство нефти.
Описанная таким образом операция моделирования может быть повторена для моделирования водных фракций в добыче нескольких эксплуатационных скважин в данной зоне месторождения нефти, связывая их со значащими данными других скважин, располагающихся в данной зоне.
Может случиться, что будут наблюдаться перекрестные взаимодействия между моделируемыми водными фракциями вследствие того, что значащие данные добычи одной или нескольких эксплуатационных скважин, соответствующие водные фракции которых были промоделированы, сами участвуют в одной или в нескольких других моделях, установленных для других эксплуатационных скважин. В этом случае осуществляют общую оптимизацию различных моделей, полученных с учетом этих перекрестных взаимодействий, таким образом, чтобы обеспечить максимальную величину общей добычи в пределах данной зоны рассматриваемого месторождения.
Правильность систематического подхода, выбранного для определения метода моделирования водной фракции в дебите данной эксплуатационной скважины, была проверена исходя из реальных данных эксплуатационного нефтяного поля некоторого стратифицированного и гетерогенного месторождения, промываемого путем нагнетания воды под давлением. В частности, история водной фракции скважины этого поля была удовлетворительным образом смоделирована при помощи саморегрессивной модели ARX, содержащей на входе месячные значения добычи или нагнетания, более или менее задержанные, на нескольких соседних скважинах.
Пример моделирования изменения водной фракции скважины W1
В данном случае рассматривается группа из двенадцати скважин, проходящих сквозь данное нефтяное месторождение и схематически представленных на фиг. 4, четырех скважин нагнетания (W4, W5, W6 и W3) и восьми эксплуатационных скважин (W7, W8, W10, W9, W2, W1 и W11). Распределение различных нагнетательных и эксплуатационных скважин W1, ..., W12 в данном случае является относительно равномерным. Расстояние между различными скважинами в данном случае составляет примерно 500 м. Приведенные ниже примеры практической реализации метода в соответствии с предлагаемым изобретением основываются на моделировании изменений водной фракции в центральной эксплуатационной или производительной скважине W1.
В данном случае рассматривается группа из двенадцати скважин, проходящих сквозь данное нефтяное месторождение и схематически представленных на фиг. 4, четырех скважин нагнетания (W4, W5, W6 и W3) и восьми эксплуатационных скважин (W7, W8, W10, W9, W2, W1 и W11). Распределение различных нагнетательных и эксплуатационных скважин W1, ..., W12 в данном случае является относительно равномерным. Расстояние между различными скважинами в данном случае составляет примерно 500 м. Приведенные ниже примеры практической реализации метода в соответствии с предлагаемым изобретением основываются на моделировании изменений водной фракции в центральной эксплуатационной или производительной скважине W1.
В данном случае подлежит идентификации следующая система: выходные данные представляют собой величину водной фракции рассматриваемой эксплуатационной скважины W1, потенциальными входными данными являются объемы нагнетаемой в недра воды и объемы вытекающей жидкости, произведенной десятью соседними скважинами, от W2 до W12.
1. Отбор значащих данных.
а) Выходные параметры
В данном случае в распоряжении исследователя имеются необработанные измерения водной фракции, реализованные при помощи отбора испытательных образцов на выходе данной скважины через весьма нерегулярные интервалы времени (от нескольких дней до одного месяца), а также среднемесячные величины, полученные в результате нахождения средних величин необработанных измерений, выполненных на протяжении календарного месяца, и при любом числе располагаемых в данном случае измерений.
В данном случае в распоряжении исследователя имеются необработанные измерения водной фракции, реализованные при помощи отбора испытательных образцов на выходе данной скважины через весьма нерегулярные интервалы времени (от нескольких дней до одного месяца), а также среднемесячные величины, полученные в результате нахождения средних величин необработанных измерений, выполненных на протяжении календарного месяца, и при любом числе располагаемых в данном случае измерений.
На фиг. 5 схематически представлены изменения необработанных данных, относящихся к водной фракции эксплуатационной скважины W1, в некоторый момент времени, рассматриваемый как начальный. Отмечаются весьма резкие изменения "высокой частоты", характеризующиеся некоторой дисперсией, связанной с шумовыми возмущениями или с ошибками измерений, вокруг некоторого более медленного изменения (на относительно низкой частоте). Здесь речь идет о выявлении этих изменений, которые соответствуют так называемым "значащим" изменениям упомянутой водной фракции (связанным с взаимными влияниями).
Решение о фильтрации "высокочастотных" компонентов может состоять, например, в использовании среднемесячных величин водной фракции, которые имеются в распоряжении исследователя, с более регулярным шагом отбора образцов, составляющим примерно 30 суток. Эти средние величины в меньшей степени подвержены воздействию шумов или помех, чем упомянутые выше необработанные измерения (см. фиг. 6), причем получение среднего значения соответствует определенной фильтрации высоких частот. Здесь легче различаются медленные изменения водной фракции, исключая самую верхнюю часть спектра частот средних величин водной фракции, схематически представленного на фиг. 7, можно получить диаграмму значащих измерений, показанную схематически на фиг. 8.
б) Входные параметры модели
Данные о расходах нагнетания и дебитах производства всех двенадцати рассматриваемых в данном случае скважин представляют собой месячные величины, выраженные в кубических метрах в месяц. На фиг. 9 и 10 можно видеть, например, изменения расходов соответственно одной из эксплуатационных скважин W11 и одной из нагнетательных скважин W4 с ежемесячным отбором испытательных образцов. Их гистограммы (не представленные на приведенных в приложении фигурах) имеют форму распределения типа гауссова распределения.
Данные о расходах нагнетания и дебитах производства всех двенадцати рассматриваемых в данном случае скважин представляют собой месячные величины, выраженные в кубических метрах в месяц. На фиг. 9 и 10 можно видеть, например, изменения расходов соответственно одной из эксплуатационных скважин W11 и одной из нагнетательных скважин W4 с ежемесячным отбором испытательных образцов. Их гистограммы (не представленные на приведенных в приложении фигурах) имеют форму распределения типа гауссова распределения.
2. Обработка измерений
В первую очередь осуществляется отбор комплекта данных с равномерным временным распределением. Для учета возможных расхождений между периодами отбора испытуемых образцов производится оценка при помощи интерполяции комплекта значений, равномерно отстоящих друг от друга во времени с достаточно точно определенным шагом (например, раз в месяц).
В первую очередь осуществляется отбор комплекта данных с равномерным временным распределением. Для учета возможных расхождений между периодами отбора испытуемых образцов производится оценка при помощи интерполяции комплекта значений, равномерно отстоящих друг от друга во времени с достаточно точно определенным шагом (например, раз в месяц).
Затем осуществляется фильтрация данных, в частности, выходных данных.
На приведенной фиг. 7, на которой схематически представлен спектр усредненных измерений водной фракции для скважины W1, можно видеть, что относительно низкие частоты обладают достаточно большой спектральной энергией, которая выражается во временной области через медленные и более значимые изменения данной водной фракции. Для устранения наиболее высоких частот с относительно малой энергией, которые можно с наибольшей вероятностью приписать шумам или помехам, а также ошибкам измерений, выбирается применение фильтрации с пропусканием низкочастотной полосы. Частота среза такого полосного фильтра принимается равной 0,5 · 10-7 Гц, что соответствует периоду среза в 231,48 суток или 7,7 месяца. Однако можно модифицировать частоту среза этого полосного фильтра и сохранить, например, пиковое значение на уровне 1,1 · 10-7 Гц для того случая, когда ситуация будет соответствовать возможной интерференции, а также проверить, будет ли данная модель усовершенствованной или нет.
Действительная диаграмма изменения водной фракции скважины W1 после фильтрации показана схематически на фиг. 8.
Производится также фильтрация входных данных. Аналогичным образом уменьшается ширина спектров, связанных соответственно с необработанными измерениями входных параметров, полученными соответственно на эксплуатационной скважине W11 (фиг. 9) и на нагнетательной скважине W4 (фиг. 11) путем применения высокочастотных фильтров, следствием чего является сглаживание результирующих диаграмм изменений (фиг. 10 и 12). Здесь можно выбрать, например, ту же частоту среза или граничную частоту, что и для выходных данных.
3. Выбор наиболее значимых входов при помощи взаимной корреляции
Система с двенадцатью входами является весьма сложной. Чем больше входов в данной модели и, следовательно, коэффициентов, тем меньше будет ошибка сглаживания модели на интервале отладки, но модель при этом будет слишком специфической на этом интервале и не будет, таким образом, надежной при экстраполяции во времени. Следовательно, предпочтительный вариант состоит в том, чтобы сохранить в модели только те входы, которые наиболее существенно влияют на характеристики выхода.
Система с двенадцатью входами является весьма сложной. Чем больше входов в данной модели и, следовательно, коэффициентов, тем меньше будет ошибка сглаживания модели на интервале отладки, но модель при этом будет слишком специфической на этом интервале и не будет, таким образом, надежной при экстраполяции во времени. Следовательно, предпочтительный вариант состоит в том, чтобы сохранить в модели только те входы, которые наиболее существенно влияют на характеристики выхода.
Для осуществления выбора этих наиболее значимых входов реализуется взаимная корреляция между выходом данной модели (отфильтрованная и усредненная величина водной фракции скважины W1) и каждым из входов этой модели. Полученные в результате этого одиннадцать функций взаимной корреляции классифицируются в порядке возрастания их максимума. На фиг. 13 можно видеть схематический пример сопоставления двух функций взаимной корреляции. Этот пример показывает, что расход на скважине W8 оказывает большее влияние на водную фракцию скважины W1, чем расход на более удаленной скважине W12, которая не может продемонстрировать существенного влияния на упомянутую водную фракцию.
4. Полученная оптимальная модель
Выход данной модели представляет собой величину, отфильтрованную и усредненную, водной фракции fw(W1) для скважины W1. Оставленные в модели входы представляют собой отфильтрованные величины расходов на следующих скважинах:
gW8: дебит добычи на скважине W8 (м3 в месяц);
gW11: дебит добычи на скважине W11 (м3 в месяц);
gW4: расход нагнетания на скважине W4 (м3 в месяц).
Выход данной модели представляет собой величину, отфильтрованную и усредненную, водной фракции fw(W1) для скважины W1. Оставленные в модели входы представляют собой отфильтрованные величины расходов на следующих скважинах:
gW8: дебит добычи на скважине W8 (м3 в месяц);
gW11: дебит добычи на скважине W11 (м3 в месяц);
gW4: расход нагнетания на скважине W4 (м3 в месяц).
fw(W1) (t)=0,9132fw(W1) (t-1)-0,6465fw(W1)(t-2)
- 0,0028g центрир. (W8)(t-1)
+ 0,5546e - 3g центрир. (W11)(t-1)
- 0,0020g центрир. (W4)(t-2)+0,0011g центрир. (W4)(t-3)
+ 69,6992
Выполненное здесь так называемое "центрирование" состоит в извлечении постоянного компонента сигнала, который представляет его среднее значение:
x (центрированный) = x - среднее значение (х).
На фиг. 13 можно сопоставить выходную величину, рассчитанную с использованием описанной выше модели (сплошная линия) с истинным выходом (пунктирная линия).
Такое сравнение показывает, что данная модель является вполне удовлетворительной и надежной: обеспечивается хорошая экстраполяция на более чем 19 предыдущих месяцев по отношению к периоду идентификации и на 6 месяцев, следующих за этим периодом, причем сама идентификация осуществляется на протяжении 16 месяцев.
Большое значение имеет выбор числа коэффициентов и задержек. Для получения устойчивой оптимальной модели необходимо сохранить возможно меньшее число коэффициентов. Задержки могут быть выбраны в соответствии с удалением "входных" скважин по отношению к "выходной" скважине.
Claims (10)
1. Способ моделирования в системе скважин, проходящих сквозь зону подземного месторождения углеводородов, в процессе эксплуатации эффектов взаимодействия между несколькими скважинами этой системы скважин, влияющих на величину водной фракции в вытекающей жидкости, производимой по меньшей мере одной эксплуатационной скважиной этой же системы скважин, подвергающейся промывке при помощи жидкости под давлением, с целью оптимизации производительности данного месторождения, отличающийся тем, что осуществляют отбор совокупности значащих данных, исходя из необработанных данных, извлеченных из истории нагнетания промывочной жидкости в данное месторождение и из истории добычи жидкого продукта, втекающего из одной или нескольких эксплуатационных скважин, и итеративное установление некоторой оптимизированной линейной модели, связывающей изменения во времени значащих данных, относящихся к водной фракции в добыче эксплуатационной скважины, с изменениями во времени значащих данных, относящихся к параметрам других скважин, входящих в упомянутую систему скважин.
2. Способ моделирования по п.1, отличающийся тем, что отбор значащих данных сопровождается частотной фильтрацией изменений необработанных данных.
3. Способ моделирования по п.1, отличающийся тем, что отбор значащих данных содержит операцию частотной фильтрации изменений необработанных данных, относящихся к водной фракции эксплуатационной скважины, с одной стороны, и необработанных данных, относящихся к другим скважинам упомянутой системы скважин, с другой стороны.
4. Способ моделирования по п.3, отличающийся тем, что отбор значащих данных содержит, например, выявление флюктуаций на некоторой низкой частоте, располагающейся значительно ниже того диапазона частот, с которым производились измерения необработанных данных, затрагивающих водную фракцию.
5. Способ моделирования по одному из пп.1 - 4, отличающийся тем, что включает выбор среди других скважин системы скважин некоторого ограниченного числа скважин, представляющих наиболее сильные взаимодействия с эксплуатационной скважиной.
6. Способ моделирования по одному из пп.1 - 5, отличающийся тем, что выбор скважин, представляющих наиболее сильные взаимодействия с данной эксплуатационной скважиной, осуществляется путем установления попарной взаимной корреляции значащих данных, связанных с величиной водной фракции этой эксплуатационной скважины, соответственно со значащими данными, связанными с другими скважинами из упомянутой системы скважин.
7. Способ моделирования по одному из пп.1 - 6, отличающийся тем, что отбор значащих данных содержит статистическую обработку необработанных данных.
8. Способ моделирования по одному из пп.1 - 7, отличающийся тем, что отбор значащих данных выполняется среди необработанных данных, принадлежащих к некоторой совокупности данных с регулярным распределением во времени.
9. Способ моделирования по одному из пп.1 - 8, отличающийся тем, что прикладывают к одной или к нескольким скважинам произвольных воздействий, модифицирующих входные необработанные данные тким образом, чтобы обеспечить наилучший выбор скважин, представляющих наиболее сильные взаимодействия.
10. Способ моделирования в системе скважин, проходящих сквозь зону подземного месторождения углеводородов, в процессе эксплуатации эффектов взаимодействия, воздействующих со стороны различных скважин этой системы скважин на величины водных фракций вытекающих жидкостей, производимых соответственно несколькими эксплуатационными скважинами, подвергающимися промывке при помощи жидкости под давлением, с целью оптимизации производительности данного месторождения, отличающийся тем, что осуществляют отбор совокупности значащих данных, исходя из совокупности необработанных данных, извлеченных из истории нагнетания промывочной жидкости в данное месторождение и истории добычи вытекающей жидкости на одной или нескольких эксплуатационных скважинах данной системы скважин, итеративное установление линейных моделей, связывающих изменения во времени значащих данных, относящихся к величине водной фракции соответственно в добыче нескольких эксплуатационных скважин упомянутой системы скважин, с изменениями во времени значащих данных, относящихся к другим скважинам этой системы скважин и общую оптимизацию различных созданных моделей с учетом перекрестных взаимодействий между значащими данными, эффективно действующими в каждой из этих моделей, таким образом, чтобы добиться максимальной общей производительности данной зоны месторождения.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9515338A FR2742794B1 (fr) | 1995-12-22 | 1995-12-22 | Methode pour modeliser les effets des interactions entre puits sur la fraction aqueuse produite par un gisement souterrain d'hydrocarbures |
FR9515338 | 1995-12-22 |
Publications (2)
Publication Number | Publication Date |
---|---|
RU96123972A RU96123972A (ru) | 1999-02-20 |
RU2165520C2 true RU2165520C2 (ru) | 2001-04-20 |
Family
ID=9485846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU96123972/03A RU2165520C2 (ru) | 1995-12-22 | 1996-12-20 | Способ моделирования влияния взаимодействия скважин на водную фракцию, производимую подземным месторождением углеводородов (варианты) |
Country Status (6)
Country | Link |
---|---|
US (1) | US5881811A (ru) |
AR (1) | AR005253A1 (ru) |
FR (1) | FR2742794B1 (ru) |
GB (1) | GB2308691B (ru) |
NO (1) | NO308868B1 (ru) |
RU (1) | RU2165520C2 (ru) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8898018B2 (en) | 2007-03-06 | 2014-11-25 | Schlumberger Technology Corporation | Methods and systems for hydrocarbon production |
RU2573746C2 (ru) * | 2010-07-30 | 2016-01-27 | Эксонмобил Апстрим Рисерч Компани | Система и способы для прогнозирования поведения скважины |
Families Citing this family (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6101447A (en) * | 1998-02-12 | 2000-08-08 | Schlumberger Technology Corporation | Oil and gas reservoir production analysis apparatus and method |
US6853921B2 (en) | 1999-07-20 | 2005-02-08 | Halliburton Energy Services, Inc. | System and method for real time reservoir management |
US6980940B1 (en) * | 2000-02-22 | 2005-12-27 | Schlumberger Technology Corp. | Intergrated reservoir optimization |
MXPA03006977A (es) | 2001-02-05 | 2004-04-02 | Schlumberger Holdings | Optimizacion de sistemas de red de reservorio, pozo y superficie. |
DE60207549D1 (de) * | 2001-04-24 | 2005-12-29 | Exxonmobil Upstream Res Co | Verfahren zur verbesserung von produktionszuweisung in einem integrierten reservoir und oberflächenströmungssystem |
AU2003276456A1 (en) * | 2002-11-15 | 2004-06-15 | Schlumberger Technology B.V. | Optimizing well system models |
US7584165B2 (en) * | 2003-01-30 | 2009-09-01 | Landmark Graphics Corporation | Support apparatus, method and system for real time operations and maintenance |
CA2514516C (en) * | 2003-03-26 | 2012-11-20 | Exxonmobil Upstream Research Company | Performance prediction method for hydrocarbon recovery processes |
US7627461B2 (en) * | 2004-05-25 | 2009-12-01 | Chevron U.S.A. Inc. | Method for field scale production optimization by enhancing the allocation of well flow rates |
WO2005121840A1 (en) * | 2004-06-07 | 2005-12-22 | Exxonmobil Upstream Research Company | Method for solving implicit reservoir simulation matrix equation |
MX2008009308A (es) * | 2006-01-20 | 2008-10-03 | Landmark Graphics Corp | Administracion del sistema de produccion dinamica. |
EP2223126B1 (en) * | 2007-12-07 | 2018-08-01 | Landmark Graphics Corporation, A Halliburton Company | Systems and methods for utilizing cell based flow simulation results to calculate streamline trajectories |
EP2247820A4 (en) * | 2007-12-13 | 2016-02-24 | Exxonmobil Upstream Res Co | PARALLEL ADAPTIVE DATA DISTRIBUTION IN A STORAGE SIMULATION USING AN UN-STRUCTURED GRID |
WO2010039326A1 (en) * | 2008-09-30 | 2010-04-08 | Exxonmobil Upstream Research Company | Self-adapting iterative solver |
US20100082724A1 (en) * | 2008-09-30 | 2010-04-01 | Oleg Diyankov | Method For Solving Reservoir Simulation Matrix Equation Using Parallel Multi-Level Incomplete Factorizations |
US9574433B2 (en) * | 2011-08-05 | 2017-02-21 | Petrohawk Properties, Lp | System and method for quantifying stimulated rock quality in a wellbore |
CN107701172B (zh) * | 2017-09-22 | 2020-07-24 | 中石化石油工程技术服务有限公司 | 基于线性模型的页岩气水平井初期最高产能的预测方法 |
CN109667568B (zh) * | 2018-12-29 | 2021-05-11 | 中国石油大学(华东) | 一种用于分层注水工艺中层段组合的确定方法及装置 |
CN111502616B (zh) * | 2019-01-30 | 2022-03-29 | 中国石油天然气股份有限公司 | 注水参数的确定方法、装置及存储介质 |
CN110500083B (zh) * | 2019-08-05 | 2022-05-10 | 中国石油天然气股份有限公司 | 一种油水井动态连通性判别方法 |
RU2720718C1 (ru) * | 2020-02-04 | 2020-05-13 | Публичное акционерное общество «Татнефть» имени В.Д. Шашина | Способ эксплуатации нефтяного пласта |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU529348B2 (en) * | 1977-10-07 | 1983-06-02 | Schlumberger Technology B.V. | Means for determining characteristics of subsurface formations |
US5444619A (en) * | 1993-09-27 | 1995-08-22 | Schlumberger Technology Corporation | System and method of predicting reservoir properties |
-
1995
- 1995-12-22 FR FR9515338A patent/FR2742794B1/fr not_active Expired - Fee Related
-
1996
- 1996-12-18 NO NO965462A patent/NO308868B1/no not_active IP Right Cessation
- 1996-12-20 US US08/769,804 patent/US5881811A/en not_active Expired - Fee Related
- 1996-12-20 AR ARP960105855A patent/AR005253A1/es unknown
- 1996-12-20 RU RU96123972/03A patent/RU2165520C2/ru active
- 1996-12-20 GB GB9626510A patent/GB2308691B/en not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8898018B2 (en) | 2007-03-06 | 2014-11-25 | Schlumberger Technology Corporation | Methods and systems for hydrocarbon production |
RU2573746C2 (ru) * | 2010-07-30 | 2016-01-27 | Эксонмобил Апстрим Рисерч Компани | Система и способы для прогнозирования поведения скважины |
Also Published As
Publication number | Publication date |
---|---|
FR2742794B1 (fr) | 1998-01-30 |
NO308868B1 (no) | 2000-11-06 |
GB9626510D0 (en) | 1997-02-05 |
AR005253A1 (es) | 1999-04-28 |
FR2742794A1 (fr) | 1997-06-27 |
NO965462D0 (no) | 1996-12-18 |
NO965462L (no) | 1997-06-23 |
GB2308691A (en) | 1997-07-02 |
GB2308691B (en) | 2000-02-23 |
US5881811A (en) | 1999-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2165520C2 (ru) | Способ моделирования влияния взаимодействия скважин на водную фракцию, производимую подземным месторождением углеводородов (варианты) | |
US7966166B2 (en) | Method for determining a set of net present values to influence the drilling of a wellbore and increase production | |
US9140108B2 (en) | Statistical reservoir model based on detected flow events | |
US5798982A (en) | Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models | |
CN109709607B (zh) | 一种预测薄层砂岩储层厚度方法及装置 | |
US8504341B2 (en) | Methods, systems, and computer readable media for fast updating of oil and gas field production models with physical and proxy simulators | |
US20110191029A1 (en) | System and method for well test design, interpretation and test objectives verification | |
AU2007211291A1 (en) | Methods, systems, and computer-readable media for fast updating of oil and gas field production models with physical and proxy simulators | |
WO2007100827A2 (en) | Monte carlo simulation of well logging data | |
Wagner | Evaluating data worth for ground-water management under uncertainty | |
CA2808858A1 (en) | Wellbore real-time monitoring and analysis of fracture contribution | |
Konrad et al. | Hydraulic behavior of fault zones in pump tests of geothermal wells: a parametric analysis using numerical simulations for the Upper Jurassic aquifer of the North Alpine Foreland Basin | |
Almeida et al. | Reducing uncertainties of reservoir properties in an automatized process coupled with geological modeling considering scalar and spatial uncertain attributes | |
Zhou et al. | Analytical reservoir simulation and its applications to conventional and unconventional resources | |
Vrabie et al. | Digital twin for downhole pressure gauges: Model and field case study | |
Bahaa et al. | Soft computation application: Utilizing artificial neural network to predict the fluid rate and bottom hole flowing pressure for gas-lifted oil wells | |
Ivšinović et al. | Application of the bootstrap method on a large input data set-case study western part of the Sava Depression | |
CN112114357B (zh) | 缝洞型储层连通性预测方法及装置 | |
Costa et al. | Evaluation of an uncertainty reduction methodology based on Iterative Sensitivity Analysis (ISA) applied to naturally fractured reservoirs | |
Zalavadia et al. | Application of a sparse hybrid data-driven and physics model in unconventional reservoirs for production forecasting | |
CN119195693B (zh) | 结合油藏分布场景的油藏隔离阀自适应控制方法及装置 | |
CN115062552B (zh) | 一种缝洞油藏氮气吞吐效果预测方法及系统 | |
RU2808507C2 (ru) | Способ определения распределения объема закачанных в скважину жидкостей по зонам пласта вдоль ствола скважины | |
Sankaran et al. | Hybrid Data-Driven and Physics-Informed Reservoir Modeling for Unconventional Reservoirs | |
Mousavi Eshkelani et al. | A new multi-level local–global history matching procedure assisted by rate transient analysis and numerical simulation |