К оглавлению журнала

 

УДК 550.83:622.276.66

© В.И. Ибраев, 1993

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

В.И. ИБРАЕВ (ЗАПСИБНИИГЕОФИЗИКА)

Использование эволюционно-тектонического подхода применительно к изучаемым (разведываемым) объектам позволяет на основе данных сейсморазведки и ГИС выйти на надежное прогнозирование и оценку аномальных пластовых давлений.

В основу методики прогноза аномальных пластовых давлений положены:

Методика прогноза реализована в пакете прикладных программ "GEODIN-EC" интерпретационной системы ИНТЕРСЕЙС и в виде аналогичного пакета на ПЭВМ IBM PC/AT.

Методика и технология прогноза аномальных давлений опробована на различных площадях Западной Сибири. В качестве примера рассмотрим решение такой задачи на бованенковской группе месторождений. На территории общей площадью 110 км проведены сейсмические исследования масштаба 1:100 000. Число профилей-160, кратность наблюдений 12-24. На исследуемой территории семь куполовидных поднятий: 1 - Бованенковское, 2 - Северо-Бованенковское, 3 - Муртинское, 4 - Нейтинское, 5 - Остромысовское, 6 - Нерстинское, 7 - Седское (рис. 1).

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

HR=Hi- [H - (Fmin+ K)]

где HR и Hi - соответственно рассчитанная и исходная глубины реконструируемой поверхности; Н - глубина поверхности приведения; Fminминимум функции поверхности приведения; К - поправочный коэффициент, определяемый по стратиграфической привязке ОВ.

Такие расчеты по поверхности фундамента проведены во временном интервале от волжского яруса до настоящего времени.

Рельеф поверхности фундамента в волжское время незначительно всхолмлен, но структурные формы Бованенковского, Остромысовского и Седского куполовидных поднятий резко выражены, прослеживается каньон широтного простирания севернее Седского куполовидного поднятия. Наибольшие изменения рельефа и связанное с этим фактором перераспределение напряжений происходило в валанжин-барремское время. В барреме завершилось формирование всех перечисленных куполовидных поднятий по амплитуде и конфигурации, кроме Муртинского (см. рис. 1). Каньон широтного простирания и Муртинская структура оформились в альбское время. Дальнейшее развитие форм рельефа носит унаследованный характер.

Для определения напряженных зон реконструированную поверхность фундамента по кровле уватской свиты, горизонт Г (сеноман), продифференцировали по секторам в 10° и с заданным шагом по направлениям Х и Y. Карта, полученная в результате суммирования производных поверхностей, представляет собой распределение напряжений поверхности фундамента (рис. 2). На нее вынесены тектонические нарушения, выявленные и протрассированные на временных разрезах.

Зонам с отметкой 0 соответствует максимальное напряжение, в целом оно изменяется от -40 до +50 усл. ед. Приуроченность разломов, выделенных по временным разрезам, к зонам максимальных напряжений свидетельствует о высокой достоверности результатов, полученных площадным способом. На карте распределения напряжений фундамента отчетливо прослеживается зональность северо-западного простирания. В области этих зон, но уже в осадочном чехле, возможно наличие залежей углеводородов с проявлениями аномальных давлений.

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

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

Количественная оценка влияния факторов литологии и насыщения на Vинт устанавливается по данным акустического каротажа (АК), выполненного в скважине, исследованной методами ВСП и СК. Анализ качества акустического каротажа в поисково-разведочных скважинах показывает, что наиболее помехоустойчивой из регистрируемых величин является характеристика Т1, которая и использована для расчета пластовых скоростей.

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

После приведения пластовых скоростей, полученных по АК, к условиям измерения ВСП построена зависимость Vp = F(H) раздельно для чистых глин, газа и водонефтенасыщенных песчаников. По данным зависимостям установлены количественные соотношения скорости продольных волн в объектах различной литологии и насыщения относительно среднестатистических глин. Из полученных соотношений с учетом коэффициентов разнонасыщенной песчанистости достаточно просто рассчитать пластовую скорость в чистых глинах по измеряемой интервальной скорости.

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

Прогноз напряженного состояния пород по данным сейсморазведки основывался на сопоставлении измеряемых интервальных скоростей, приведенных к условиям глинистого разреза, с интервальными скоростями продольных волн в нормально уплотняющихся с глубиной глинах (без аномально высоких поровых давлений - АВПоД). Следовательно, в каждом изучаемом разрезе необходимо иметь зависимость Vнорм.гл = F(H). Эту зависимость можно получить двумя способами: по данным АК в глинистых интервалах без АВПоД и расчетным способом, включая интервалы с аномальным поровым давлением по данным сейсморазведки.

Первый способ предполагает наличие глинистых отложений без АВПоД во всем диапазоне глубин изучаемого разреза. На практике это условие не всегда выполняется. В частности, на Бованенковской площади на глубине более 2300 м такие отложения отсутствуют, что не позволяет охарактеризовать основные аномально проявляющиеся отложения в интервалах органического вещества М, Во, Б и глубже. Поэтому в ППП "GEODIN-EC" был реализован второй подход, в котором решается прямая задача нахождения Уинт в нормально уплотняющихся глинах по известным его величинам, коэффициентам аномальности поровых давлений и газо- и водонасыщенной песчанистости.

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

По изложенной методике произведена оценка коэффициента аномальности в интервалах между сейсмическими границами Г, Г2, М1, М, Во, Б, Т2, А. Для этих интервалов были построены карты среднеинтервальных значений коэффициента аномальности пластовых и поровых давлений (рис. 3).

Для оценки достоверности полученных результатов выполнено сопоставление среднеинтервальных коэффициентов аномальности, определенных по результатам работы ППП "GEODIN-EC". и материалам ГИС в скважинах, не участвовавших в настройке интерпретационного алгоритма.

В пакете реализовано два варианта прогноза: с использованием единой зависимости изменения с глубиной скорости нормального уплотнения глин и поинтервальной оценкой положения "линии глин".

Имеющиеся статистические данные по материалам с аномально высоким пластовым и поровым давлением не позволяют однозначно отдать предпочтение первому или второму варианту. На рис. 4 приведены результаты по варианту 1 в сопоставлении с данными, полученными по материалам ГИС. Принимая результаты интерпретации материалов ГИС как достаточно апробированные и многократно подтвержденные испытанием скважины за истину, можно проведенное сопоставление считать доказательством достаточно высокой эффективности результатов площадного прогноза. Так, 92 % приведенных прогнозных значений удовлетворяют требованиям точности оценки технологических параметров бурения в зонах развития аномально высоких пластовых давлений.

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

ABSTRACT

The problem of areal forecast of anomaly pressures in terrigenous deposits is regarded on the base of geological-geophysical data. Effective forecasting is based on evolutional-tectonic approach for seismic and logging data analysis. The method of forecast is realized in the packet of applied programms "GEODIN-EC" of interpretation system INTERSELS and is used at one of Western Siberia fields.

РИС. 1. ПЛОЩАДНАЯ РЕКОНСТРУКЦИЯ ПОВЕРХНОСТИ ФУНДАМЕНТА БОВАНЕНКОВСКОЙ ГРУППЫ МЕСТОРОЖДЕНИЙ

РИС. 2. КАРТА РАСПРЕДЕЛЕНИЯ НАПРЯЖЕНИЙ ПОВЕРХНОСТИ ФУНДАМЕНТА.


Цифры в кружках
- куполовидные поднятия

РИС. 3. КАРТЫ СРЕДНЕИНТЕРВАЛЬНЫХ ЗНАЧЕНИЙ КОЭФФИЦИЕНТА АНОМАЛЬНОСТИ ПЛАСТОВЫХ И ПОРОВЫХ ДАВЛЕНИЙ В ИНТЕРВАЛАХ СЕЙСМИЧЕСКИХ ГРАНИЦ:

а - Во - Б, б - Б - T2, в - Т2 А

РИС. 4. СОПОСТАВЛЕНИЕ РЕЗУЛЬТАТОВ ПО ГИС И СЕЙСМОРАЗВЕДКЕ