К оглавлению журнала | |
УДК 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)]
где H
R и 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".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. КАРТА РАСПРЕДЕЛЕНИЯ НАПРЯЖЕНИЙ ПОВЕРХНОСТИ ФУНДАМЕНТА.
а
- Во - Б, б - Б - T2, в - Т2 – А