WWW.MASH.DOBROTA.BIZ
БЕСПЛАТНАЯ  ИНТЕРНЕТ  БИБЛИОТЕКА - онлайн публикации
 

«ШАНДРЫГОЛОВ ЗАХАР НИКОЛАЕВИЧ СИСТЕМА ПРОГРАММНО-ИНФОРМАЦИОННОГО ОБЕСПЕЧЕНИЯ ПРОЦЕССА АДАПТАЦИИ ГЕОЛОГОТЕХНОЛОГИЧЕСКИХ МОДЕЛЕЙ ГАЗОВЫХ ЗАЛЕЖЕЙ СЕВЕРА ЗАПАДНОЙ СИБИРИ ...»

На правах рукописи

ШАНДРЫГОЛОВ ЗАХАР НИКОЛАЕВИЧ

СИСТЕМА ПРОГРАММНО-ИНФОРМАЦИОННОГО

ОБЕСПЕЧЕНИЯ ПРОЦЕССА АДАПТАЦИИ ГЕОЛОГОТЕХНОЛОГИЧЕСКИХ МОДЕЛЕЙ ГАЗОВЫХ ЗАЛЕЖЕЙ СЕВЕРА

ЗАПАДНОЙ СИБИРИ

05.13.01 - Системный анализ, управление и обработка информации

(нефтегазовая отрасль)

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

5 ДЕК 2013 Тюмень-2013

Работа вьшолнена в ООО «ТюменНИИгипрогаз» и в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Тюменском государственном нефтегазовом университете» на кафедре «Моделирование и управление процессами нефтегазодобычи»

Научный Красовский Александр Викторович руководитель: кандидат технических наук Заместитель генерального директора по научным и проектным работам в области разработки и эксплуатации газовых месторождений ООО «ТюменНИИгипрогаз»

Тян Владимир Константинович Официальные доктор технических наук, профессор, оппоненты:

декан нефтетехнологического факультета самарского государственного университета, заведующий кафедрой «Трубопроводный транспорт»

Лазарева Валентина Георгиевна кандидат технических наук заместитель начальника отдела моделирования ОАО «НК «Роснефть» - ЗАО «РОСПАН ИНТЕРНЕШНЛ»



Ведущая организация: ООО «Геология резервуара» ГЕОТЕКС Холдинг, г. Тюмень

Защита состоится «28» декабря 2013 г. в «14°°» на заседании диссертационного совета Д 212.273.08 при ФГБОУ НПО «Тюменский государственный нефтегазовый университет» по адресу: 625000 г. Тюмень, ул. Володарского, д. 38, ауд. 219 (1-й корпус) .

С диссертацией можно ознакомится в библиотечно-издательском комплексе ФГБОУ НПО «Тюменский государственный нефтегазовый университет», по адресу 625039, г. Тюмень, ул. Мельникайте, 72 .

(www.tsogu.ru)

Автореферат разослан «27» ноября 2013 г .

Ученый се1фетарь Руднева Лариса диссертационного Николаевна совета

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. Согласно данным ведущих статистических агентств к 2030 году мировой спрос на природный газ вырастет в среднем на 30%, что является хорошей возможностью для закрепления позиции России как крупнейшей газодобываюшей державы в доле мировой экономики. В условиях роста спроса и снижения добывных возможностей в газовой отрасли приоритетными являются не только задачи освоения новых месторождений, но и эффективной эксплуатации уже находящихся в промышленной разработке запасов газа .

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

Однако стоит учесть, что технологические показатели, получаемые при расчёте геолого-технологической модели, имеют погрешность, величина которой зависит от качества модели, которая в свою очередь зависит от точности определенных исходных параметров .



Большое количество исходных параметров определено с большой долей ошибки, что вызвано либо отсутствием данных по конкретному объекту моделирования, либо в погрешности определения того или иного параметра при проведении исследований .

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

Описанный выше процесс определен во многих источниках как "адаптация модели па исторически данные" и требует большое количество времени и ресурсов как трудовых, так и вычислительных .

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

Объект исследования — Геолого-технологические модели газовых залежей .

Предмет исследования - Методы адаптации геологотехнологических моделей залежей УВС .

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

Для достижения цели работы поставлены следующие задачи:

1) На основе анализа существующих методов адаптации выявить их недостатки, определить наиболее корректный метод, применимый для геолого-технологических моделей газовых залежей;

2) Путем проведения численных экспериментов исследовать и количественно оценить влияние адаптируемых параметров геологотехнологической модели на динамику пластового давления и провести их ранжирование .

3) Выявить области корректировки адаптируемых параметров по кубам фильтрационно-емкостных свойств геолого-технологической модели путем анализа среднеквадратичного отклонения и величины градиентов пластового давления .

4) Разработать метод построения куба абсолютной проницаемости, основанный на обработке результатов газодинамических исследований, методах декомпозиции и композиции геологотехнологической модели .

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

6) Реализовать разработанные методы и алгоритмы в системе программно-информационного обеспечения .

Методы исследований. Решение поставленных задач осуществлялось на основе изучения литературных источников, научных публикаций, посвященных гидродинамическому моделированию залежей УВС, а также методов системного и сравнительного анализа, численных математических методов, теории разработки газовых месторождений. Численное решение уравнений фильтрации проведено с помощью гидродинамического симулятора «Eclipse» и методов подземной гидродинамики .





Научная новизна работы:

1) На основе анализа величин градиентов и погрешности пластового давления обосновано распределение корректирующих коэффициентов по кубам объемных параметров гидродинамической модели .

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

3) Предложены алгоритмы систематизации геолого-промысловой информации и ее обработки, обеспечивающие возможность расчета размерных, продуктивных и физических характеристик аналитического аквифера .

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

Практическая значимость работы:

Предложенные методы алгоритмизированы. Создано программное обеспечение, которое существенно ускоряет процесс создания геологотехнологических моделей и повьппает качество расчетов .

Прогнозные технологические показатели, полученные по результатам расчета геолого-технологических моделей, адаптированных при помощи описанных методов и алгоритмов, использовались на этапе проектировании разработки Юбилейного, Ямсовейского, Медвежьего, Южно-Русского, Западно-Таркосалинского, Комсомольского месторождений .

Основные положения, выносимые на защиту:

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

2) Метод определения абсолютной проницаемости по кубу гидродинамической модели с учетом промысловых исследований и приемов секторного моделирования .

3) Новый метод распределения корректирующих коэффициентов по кубам объемных параметров гидродинамической модели, основанный на анализе и обработке массивов значения градиентов и погрешности пластового давления .

Соответствие диссертации паспорту научной специальности .

Диссертация соответствует п.4 - Разработка методов и алгоритмов решения задач системного анализа, оптимизации, управления, принятия решений и обработки информации, п.5 - Разработка спехщального математического и программного обеспечения систем анализа, оптимизации, управления, пршеттия решений и обработки информации, п.6 - Методы идентификации систем управления на основе ретроспективной, текущей и экспертной информации, области исследований научной специальности 05.13.01 - Системный анализ, управление и обработка информации (нефтегазовая отрасль) .

Апробация работы .

Основные положения работы доложены и обсуждены на следующих научных и научно практических конференциях:

- региональная научно-техническая конференция «Компьютерное моделирование и системный анализ в нефтегазовой отрасли и образовании», Тюмень, 2010 г.;

- ХУ! научно-практическая конференции молодых ученных и специалистов ООО "ТюменНИИгипрогаз" «Проблемы развития газовой промышленности Сибири», Тюмень, 2010 г.;

— IV научно-практическая конференции молодых ученых и специалистов «Газовой отрасли - энергию молодых ученых», "СевКавНИПИгаз", Ставрополь 2011 г.;

— XVII научно-практическая конференции молодых ученных и специалистов ООО "ТюменНИИгипрогаз" «Проблемы развития газовой промышленности Сибири», Тюмень, 2012 г.;

— научно-технические совещания в ООО «Газпром добыча Ноябрьск», ООО «Газпром добыча Надым» и ОАО «Газпром» в 2011-2013 гг .

Публикации .

Основное содержание работы

отражено в 8 публикациях, из них 2 - в изданиях, рекомендованных ВАК Минобрнауки РФ .

Структура и объем работы. Диссертация состоит из введения, пяти разделов, заключения и списка литературы. Объем диссертации составляет 125 страниц, содержит 52 рисунка, 2 приложения. В библиографии представлено 102 наименований работ российских и зарубежных авторов .

СОДЕРЖАНИЕ РАБОТЫ

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

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

В целом принцип гидродинамического моделирования газовых залежей можно описать уравнением двухфазной фильтрации (1) .

Ы Ы """" '

–  –  –

где: д - скорость фильтрации, - коэффициенты абсолютных проницаемостей по осям координат, tj - динамическая вязкость, рплотность флюида, crcos(9 - произведение коэффициента поверхностного натяжения на косинус краевого угла смачивания на границе газ-водапорода, J - функции, вязкости, т - коэффициент пористости .

Система решается с использованием численных методов. Существует большое количество коммерческих программ, реализующих численное решение уравнений фильтращш. К наиболее распространенным можно отнести: Eclipse (ScMimiberger), Tempest (Roxar), TimeZYX и tNavigator (RFD) .

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

Этап адаптации пластового давления в геолого-технологической модели залежи является одним из основополагающих, так как именно качество адаптации модели по пластовому давлению определяет правильность расчета технологических показателей на прогнозный период .

Большинство процессов автоматической адаптации можно разделить на два этапа:

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

Так, в работах Яанса для адаптации пластового давления используется метод нелинейной регрессии. Пласт разделяется на зоны с различными значениями проводимостей и удельных объемов. Давление в модели является функцией этих двух параметров, которые изменяются формальным способом.

С помощью регрессионного анализа выбираются значения этих независимых параметров, при которых минимизируется следующее соотношение (2):

min: Е = Wi (pi^^n " Р1вычислУ (2) где: Wi-весовой коэффициент (коэффициент доверия) наблюдаемого значения; р„аб,-наблюдаемое значение (факт); Р1еычисл- значение, вычисленное по модели .

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

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

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

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

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

Так, например, для получения диапазона варыфования абсолютной проницаемосга пласта сопоставлены значения, полученные по результатам геофизических и керновых исследований. Охфеделена относительная погрешность, которая варьируется от 8 до 15% в зависимости от месгорояодения .

По результатам анализа установлено, »гго процесс адаптации пластового давления имеет наибольший эффект при модификации таких параметров, как:

1) объемные параметры залежи (изменение данного параметра на ±5% позволяет скорректировать пластовое давление от 7 до 17 %);

2) фильтрационные характеристики пород (проницаемость) (изменение данного параметра на ±20% позволяет скорректировать пластовое давление от 3 до 10 %);

3) параметры водонапорных горизонтов (изменение данного параметра на ±15% позволяет скорректировать пластовое давление от 3 до 7 %);

4) изменение других параметров оказывает влияние на пластовое давление менее 3%, что незначительно .

Порядок корректировки параметров определен исходя из степени его влияния на пластовое давление .

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

Процесс адаптации осуществляется по алгоритму, приведенному в виде блок-схемы на рисунке 1 .

0,93 1,05, Рисунок 1 - Алгоритм процесса адаптации пластового давления в гидродинамической модели газовой залежи .

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

В разделе рассмотрены решения задачи адаптации пластового давления в таких программных продуктах как «Н1з1огу+», «ЕпАВЬЕ»

(Кохаг) и «8тОрЬ (Schlumbeгger) .

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

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

Для адаптации исходными инструментами выступает два набора точек: модельные и реальные. При этом каждой модельной точке должна соответствовать только одна реальная точка, относящаяся к тому же типу и объекту .

Разницу между замеренным и фактическим показателем обозначим за Ri.

Тогда:

Ri = Wi* wa(0i - Ci), (4) где: Wi.Wci- коэффициенты доверия; Oj- наблюдаемое значение; q значение, вычисленное по модели .

Каждой паре точек будет соответствовать один остаток .

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

Mean Square - средний корень квадрата ошибки, далее - СКО или а):

(5) а

–  –  –

Коэффициент 1/2 не влияет на ход и результат решения и необходим для более простого нахождения производной .

Для соблюдения условия минимума функции нескольких

–  –  –

оценки погрешности пластового давления, который представлен на рисунке 2 .

Рисунок 2 - Блок-схема алгоритма распределения корректирующих коэффициентов по кубам объемных параметров гидродинамической модели на основе анализа величин градиентов и оценки погрешности пластового давления В третьем разделе описан метод построения куба абсолютной проницаемости, основанный на обработке результатов газодинамических исследований, методах декомпозшщи и композиции геологотехнологической модели. Под секторной моделью подразумевается гидродинамическая модель высокой точности, ограниченная определенной площадью вокруг исследуемого объекта. В качестве исследуемого объекта выступает сквагаша (либо группа скважин). Целью построения секторных моделей в рамках поставленной задачи является воспроизведение кривой восстановления давления (далее КВД) по скважинам во время их остановки и уточнение абсолютной проницаемости по вскрытой скважиной части пласта по результатам интерпретации гидродинамических исследований .

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

Интерпретация КВД происходит при помощи специализированного программного обеспечения Saphir (KAPPA) и PanSystem. В результате интерпретации определяется большое количество параметров

–  –  –

где: К — корреляционный коэффициент; К1 „р - проницаемость г-ой ячейки, полученная по зависимости; к ^ и к п - соответственно первая и последняя ячейка вскрываемая скважиной; проницаемость, полученная по результатам интерпретации КВД. Щ эфф - эффективная толщина г-ой ячейки, вычисляемая как:

где: Ог^ - толщина /-ой ячейки; ЫТС1 - песчанистость /-ой ячейки .

Далее проницаемость, изначально заданная во вскрытых скважинами ячейках, умножается на соответствующей корреляционный коэффициент .

Таким образом, если в формуле (10) К1„р принять как уже скорректированную, то корреляционный коэффициент будет равен 1 .

Учитывая, что коррекция проницаемости по формуле (10) касается только ячеек, вскрытых скважинами, необходимо распределение (интерполяция) ее значений по всей модели залежи. Для этого, исходя из глубины середины вскрытой ячейки, получено соответствие абсолютной глубины к пересчитанной проницаемости с привязкой к инклинометрии скважин. Полученные данные загружаются в пакет геологического моделирования, где значения пересчитанной проницаемости распределяются в межскважинном пространстве, а также в ячейках, не вскрытых перфорацией, использовался алгоритм интерполяции Кригинга .

Алгоритм разработанного метода в виде блок-схемы представлен на рисунке 3 .

В четвертом разделе рассмотрен метод обработки и систематизации геолого-промысловой информации для расчета размерных, продуктивных и физических хараьстеристик аналитических моделей водонапорных пластов. За основу взяты аналитические модели водонапорных пластов Картера-Трейси и Фетковича .

Предполагается, что давление на внешней границе водоносного пласта не меняется.

Тогда интенсивность притока вычисляется как:

где: р^^- давление на внешней границе пласта, предполагаемое ПОСТОЯННЫМ, Па; р - среднее давление на границе водоносного и продуктивного пластов (предполагается равньв! среднему давлению рд в водоносном пласте). Па .

Рисунок 3 - Блок-схема алгоритма распределения абсолютной проницаемости по площади и разрезу газовой залежи с помощью секторного моделирования .

В конечно-разностном виде по аналогии определяется С^у = С^^у .

где alJ - коэффициент распределения притока по блокам. Положим V СР"^^ + Р") дая интервала от С" до Тогда — См,..( р ^ - рГ;) - \Сац^Р''^' \г -_ 1 (12) = .

Переменный водоносный горизонт представлен так называемой моделью в виде аоршка» .

Если водоносный пласт относительно мал и имеет непроницаемые границы, то во все моменты времени данный пласт будет находиться в приблизительном равновесном состоянии с продуктивным пластом .

Поэтому фильтрация будет наблюдаться только при изменении давления на границе продукпшного и водоносного пластов, т. е. при Я А - - с - УрА где — ингенсивность притока в единицах объема в пластовых условиях, сут.; с - общая сжимаемость (воды и породы) в водоносном пласте (с = с^ -I- Сц), Па''; Ур^- поровый объем водоносного пласта, м' .



Конечное разностное уравнение (1.17) принимает вид:

–  –  –

где Ац- площадь сечения, через которую происходит приток из водоносного пласта в блоке у, м^ .

Суммирование площадей осуществляется по всем блокам, где имеется приток из водоносного пласта. В неоднородном пласте приток лучше распределять в соответствии с проводимостями .

Выше описаны два крайних случая. Для водоносного пласта в виде «горшка» приток не зависит от давления (среднее давление в водоносном пласте меняется вслед за давлением р). Для стащюнарного водоносного пласта не изменяется: pJ^ = В действительности интенсивность притока зависит от промежуточных условий, ее можно получить решением уравнения нестационарной фильтрации в водоносном пласте. Данный подход был развит Ван Эвердингеном и Херстом (1949) для аналитических расчетов, что позволяет определить для некоторых простых случаев .

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

Обычно эти результаты получают с помощью «функции влияния» QI(t') и по формулам:

для постоянного давления:

' (16) = -р);

о для постоянного расхода:

р Т - р а ) = Р1Шл- (17) Функции Q/(t) и P/(t) зависят от геометрии и свойств водоносного пласта. Для многих случаев функцж влияния табулированы. Так как для реального водоносного пласта р и переменны, точное решение можно получить по принципу суперпозиции в сочетании с расчетом материального баланса по пласту. Упрощенный, пригодный для применения на ЭВМ метод предложен Картером и Трейси (1960). По этому методу для расхода на границе используется функция влияния

PI(t)H в результате получено выражение:

Qm/ = + - Р")у], (18) где a(t) и b(t) - функции P/(t) и полного притока W(t) .

Простой подход, пригодный для произвольного конечного водоносного пласта, предложен Фетковичем (1971).

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

M-Q^tj' = - рГ'). (19) где Ра - среднее давление в водоносном пласте в момент времени Па; - коэффициент продуктивности для водоносного пласта .

После завершения шага по времени определяют полный приток и с помощью уравнения материального баланса подсчитывают новое значение для водоносного пласта. Феткович показал, что для определения /а на временном шаге даже для нестационарного режима водоносного пласта можно использовать псевдостационарные или стационарные уравнения. Таким образом, его метод совершенно аналогичен методу учета отдельных скважин, и при соответствующем выборе функции расхода учитьшает все три типа поведения водоносного пласта .

На основе более подробного описания математических моделей Фетковича и Картера-Трейси составлены алгоритмы расчета их параметров .

В пятом разделе описывается проведение численных экспериментов на гидродинамических моделях месторождений западной Сибири .

Для адаптации модели последовательно применены все три метода .

Оценочным параметром настройки является среднеквадратичное отклонение модельно пластового давления от значений фактических замеров .

На первом этапе применен метод, основанный на анализе среднеквадратичного отклонения и величины градиентов пластового давления, позволяющий с помощью методов интерполяции распределить области корректировки адаптируемых параметров по кубам фильтрационноемкосгных свойств геолого-технологической модели (далее - метод 1) .

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

На втором этапе, согласно предложенному методу построения куба абсолютной проницаемости, основаьшый на обработке результатов газодинамических исследований, методах декомпозиции и композиции геолого-технологической модели (далее - метод 2), основьшаясь на результатах ГДИС и интерпретации 1СВД в специализированньк программных продуктах, выполнено построение секторных моделей в районе кустов скважин. Пример модели куста скважин Юбилейного месторояодения представлен на рисунке 4 .

Исходя из результатов интерпретации газодинамических исследований по скважинам, во вскрываемых скважинами частях пласта уточнена абсолютная проницаемость породы. Далее уточненная по скважинам проницаемость распределена по всему моделируемому объекту .

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

Ниже в таблице 1 представлена динамика СКО при применении всех трех описанных методов .

Таблица 1 - Динамика СКО

–  –  –

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

Это позволяет сделать вывод о высокой точности полученных по этой модели технологических показателей разработки на прогнозный период .

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1) Анализ существующих методов адаптации геолого-технологических моделей позволил выявить, что к наиболее существенным недостаткам, среди рассмотренных методов, можно отнести отсутствие единственного рещения и получение "нефизичных" значений адаптируемых параметров. К наиболее подходящим методам для адаптации геолого-технологических моделей можно отнести градиентные методы, которые исключают ряд недостатков обнаруженных в рассмотренных методах и имеют хорошею сходимость .

2) В результате численных экспериментов установлено, что на процесс адаптации наиболее существенное влияние оказывает:

- объемные параметры залежи (изменение данного параметра на ±5% позволяет скорректировать пластовое давление от 7 до 17 %);

- фильтрационные характеристики пород (проницаемость) (изменение данного параметра на ±20% позволяет скорректировать пластовое давление от 2 до 10 %);

- параметры водонапорных горизонтов (изменение данного параметра на ±15% позволяет скорректировать пластовое давление от 1 до 7 %);

3) Установлено, что оценка градиентов изменяемых параметров геологотехнологической модели газовой залежи и на ее основе корректировка порового объема позволяет повысить точность адаптации на 20-40% .

4) С помощью разработанного метода распределения абсолютной проницаемости по кубу гидродинамической модели с учетом гидродинамических исследований и приемов секторного моделирования возможно снижения среднеквадратичного отклонения до 45% .

5) Доказано, что оценка и расчет параметров водоносного аналитического аквифера с помощью разработанного автором метода позволяет повысить точность адаптации пластового давления до 30%

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

7) На основе разработанных алгоритмов создана система программного обеспечения, позволяющая автоматизировать процесс адаптации пластового давления геолого-технологических моделей сеноманских газовых залежей, что приводит к снижению временных затрат на 50-60% .

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

Ведущие периодические издания:

1. АВ. Красовский. Создание математической модели уникального газового месторождения / АЗ. Красовский, СЮ. Свенгский, З.Н. ПЬвдрыголов, А.О. Лысов // Научно-технический журнал. Бурение и нефхъ - 2013-№11.-С.26-28 .

2. А.В. Красовский. Новые математические методы адаптации геолого-технологических моделей сеноманских газовых залежей / А.В .

Красовский, С.Ю. Свенгский, З.Н. Шандрыголов, М.А. Казанцев // Журнал «Нефть и газ» - 2013-№11.-С.110-112 .

Другие издания:

3. А.О. Лысов. Методика выбора п^аметров адаптации пластового давления гидродинамических моделей газовых месторождений / А.О. Лысов, З.Н .

Шащфыголов, ело. Светский // Новые информационные технологии в нефтегазовой отрасли и офазовании. Материалы V всероссийской научнотехнической конференции с междун^одным учасшем - 2012 - С. 3 7 ^ .

4. АВ. К]расовский. Адашация пластового давления модели сеноманской газовой залежи Западно-т^косалинского НГКМ / АВ. Красовский, СЮ .

Свенгский, АО. Лысов, 3 Л. Шацдрыголов // Материалы восьмой Всероссийской научно-технической конференции «Геология и нефтегазоносность Западносибирского мегабассейна (опьп-, инновации)» - 2012 - С. 50-53 .

5. АВ. Красовский. Адаптация пластового давления гидродинамических моделей газовых залежей на историю разработки с использованием функции уплотнения пород пласта / АВ. Красовский, СЮ. Свенгский, АО. Лысов, З Л Шшщзыголов // Сборник научных трудов ООО «ТюменНИИгипрогаз»-2013-С. 27-30 .

6. АВ. 1фасовский. Моделирование разработки нефтяной оторочки пласта АПу Западно-Т^косалинского месторождения / АВ. Щзасовский, МГ. Мавлевдинов, З.Н. Шандрыголов // Сборник тезисов и докладов ООО «ТюменНИИгипрогаз»

'Проблемы развитая газовой промышленносш Сибири"- 2010 - С. 34-37 .

7. А.В. Красовский. Алгоритм расчета п^аметров модели водоносного горизонта на примере сеноманской залежи Медвежьего месторождения / АВ. Ьфасовский, Е.В. Рауданен, З.Н. Шандрыголов С.Ю. Свенгский // Сборник тезисов и докладов ООО «Газпром добьиа Надым» "УП Научно-практическая конференция молодых специалистов и ученых"- 2013 - С. 24-25 .

8. В.В. Воробьев. Моделирование режимов работы экспериментальной двухзабойной скважины на туронской газовой залежи / В.В. Воробьев, М.Г .

Мавлетдинов, А.В. Красовский, АС. Широкова З.Н. Шандрыголов // Научнотехнический журнал. Наука и ТЭК. - 2011.-№ 1.-С. 19-20 Подписано к пеяаш 26.11.2013 г .

Формат бумага 60x841/16 .

Усл. печ. л. 1,00. Заказ № 191. Тираж 150 экз .

ООО «ТюменНИИгапрогаз», ООВ




Похожие работы:

«Проект на 27.07.2015г. ПРОГРАММА Всероссийский День Строителя 2015 Конгресс-парк гостиницы "Radisson Royal", г. Москва, Кутузовский проспект, 2/1, стр. 1. Организатор: Министерство строительства и жилищно-коммунального хозяйства Российской Федерации, при участии Правительства Москвы, Правительст...»

«БК-20400М Паспорт БК-20400М Паспорт Блок конденсаторов БК-20400М ПАСПОРТ г. Киев. 2013 г. © НПП "ВЭЛ" 2013 © НПП "ВЭЛ" 2013 БК-20400М Паспорт БК-20400М Паспорт 1. Введение Настоящий паспорт является документом, удостоверяющим гарантированные предприятием-изготовителем основные параметры и технические характеристики блока конденсат...»

«Общество с ограниченной ответственностью "ПРОСАМ" ОКП 40 1760 УТВЕРЖДАЮ Генеральный директор ООО “ПРОСАМ” С.И. Михель “” 2007 г. МАШИНА ЭЛЕКТРОННАЯ КОНТРОЛЬНО–РЕГИСТРИРУЮЩАЯ “ОКА–102К” Инструкция по эксплуатации ПРАУ.466137.013–20 ИЭ Подпись и дата СОГЛАСО...»

«С.А. Колюбин ДИНАМИКА РОБОТОТЕХНИЧЕСКИХ СИСТЕМ Санкт-Петербург МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ УНИВЕРСИТЕТ ИТМО С.А. Колюбин ДИНАМИКА РОБОТОТЕХНИЧЕСКИХ СИСТЕМ РЕКОМЕНДОВАНО К ИСПОЛЬЗОВАНИЮ В УНИВЕРСИТЕТЕ ИТМО по направлениям 15.04.06 Мехатроника и роб...»

«отзыв Официального оппонента на диссертационную работу Исламова Рустэма Рильевича на тему "Совершенствование системы мониторинга технического состояния протяженных участков магистральных нефтегазопроводов применением волоконно-оптических сенсоров деформаций", представленную на соискание ученой степени...»

«База нормативной документации: www.complexdoc.ru ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР НЕФТЕПРОДУКТЫ МЕТОД ОПРЕДЕЛЕНИЯ ВЯЗКОСТИ АВТОМАТИЧЕСКИМ КАПИЛЛЯРНЫМ ВИСКОЗИМЕТРОМ ГОСТ 7163-84 ГОСУДАРСТВЕННЫЙ КОМИТЕТ СССР ПО СТАНДАРТАМ Москва ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР НЕФТЕПРОДУКТЫ Метод определения вяз...»

«В последнее десятилетие во многих развитых странах, в том числе и в Бетребителей – это молодые люди в возрасте от 15 до 25 лет. Наркомания среларуси, наблюдается тенденция к увеличению среди де...»

«Магфурова Светлана Олеговна, Маклакова Евгения Михайловна К ВОПРОСУ О ЛЕКСИКО-ГРАММАТИЧЕСКОЙ ХАРАКТЕРИСТИКЕ БЕЗЛИЧНЫХ ГЛАГОЛЬНЫХ ЛЕКСЕМ СО ЗНАЧЕНИЕМ ФИЗИОЛОГИЧЕСКОГО СОСТОЯНИЯ В рамках теории лексической грамматики при определении лексико-грамматического статуса безличных глагольных лексем в ст...»

«Лекция 1. Разреженные газы, вакуум, вакуумные устройства в криогенной технике: сосуды Дьюара, вакуумные насосы. Вакуум (от лат. vacuum — пустота) — состояние материи в отсутствии вещества. Также его иногда называют безвоздушным пространством, х...»




 
2019 www.mash.dobrota.biz - «Бесплатная электронная библиотека - онлайн публикации»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.