DO: 10.24411/2619-0761-2020-10025 УДК 550.834.048.05
УТОЧНЕНИЕ СТРОЕНИЯ МЕСТОРОЖДЕНИЯ ЗЫХ ПРИМЕНЕНИЕМ АТРИБУТНОГО АНАЛИЗА ДАННЫХ СЕЙСМОРАЗВЕДКИ 30
Ахмедов Т.Р.
Азербайджанский государственный университет нефти и промышленности, г. Баку, Азербайджан
*E-mail: akhmedov. tofik@bk.ru
Аннотация. Статья посвящена кинематической интерпретации данных сейсморазведки 3D на площади Зых. Даются краткие сведения о площади исследования, об ее географическом положении, об отложениях, вскрытых глубоким бурением, о строении складки Зых, об истории геолого-геофизического изучения. Отмечается, что основными целями настоящих исследований является построение структурного каркаса по пластам балаханской свиты (VIII и IX) и пластам ПК и КаС, который будут служить основой построения 3D-геологических моделей месторождения. Предметом объемного моделирования геологического строения являются: форма залежей (структурный план, морфология поверхности флюидных контактов), внутренняя структура резервуара и залежей (пространственное взаимоотношение коллекторов и неколлекторов, распределение в разрезе различных по типу флюидов), фильтрационно-емкостные свойства пород-коллекторов и их нефте-, газонасыщенность.
Настоящие исследования охватывают вопросы, касающиеся структурного плана целевых горизонтов, т.е. форм залежей. Очевидно, что только пликативным вариантом строения невозможно составлять адекватную объемный модель геологического строения. Построены структурные карты по пластам балаханской свиты (VIII и IX), а также по кровле подкирмакинской свиты (пласт ПК).
Очевидно, что только пликативным вариантом строения невозможно составлять адекватную объемный модель геологического строения. Чтобы получить наиболее точную модель строения необходимо построить модели разломов, для чего использовался программный пакет PETREL.10.2. Для построения структурных карт использованы карты, построенные на производстве в пликативном варианте и трехмерные модели разрывных нарушений, построенные по результатам атрибутного анализа сейсмических данных; были использованы в основном два куба атрибутов сейсмического волнового поля - кубы когерентности и погружения. Для каждого из моделируемых пластов были построены карты атрибутов когерентности и погружений.
Отмечается, что данные, полученные на основе атрибутного анализа (DEP и когерентности), использованы для получения модели разломов и уточнения геологического строения месторождения Зых и для построения структурных карт в дизъюнктивном варианте целевых горизонтов. На основе построенной модели разломов и пликативных структурных карт построены структурные карты в дизъюнктивном варианте и объемные модели тектонических блоков (сегменты 3D-модели).
Наконец, полученные по атрибутному анализу модели сегментов для отдельных отделов продуктивной толщи нижнего плиоцена получили свое отражения в структурных картах в дизъюнктивном варианте целевых пластов, которые являлись конечными результатами настоящих исследований.
Ключевые слова: сейсморазведки 3D, месторождение Зых, балаханская свита (VIII и IX), подкирмакинская свита, атрибутный анализ.
О) ®
Содержимое этой работы может использоваться в соответствии с условиями лицензии Creative Commons Attribution 4.0. Любое дальнейшее распространение этой работы должно содержать указание на автора (ов) и название работы, цитирование в журнале и DOI.
Введение.
Месторождение Зых расположено в юго-восточной части Абшеронского полуострова, на восточной окраине гор. Баку, на южном продолжении платообразной возвышенности, простирающейся от сел. Сураханы на севере до Каспийского моря на юге.
На месторождении обнажены отложения верхнего и среднего абшерона (рис. 1).
Скважинами вскрыт разрез всего плиоценового комплекса [1]. Глубоким бурением вскрыт разрез всего плиоценового комплекса отложений до понтического яруса. Лито-логически он ничем не отличается от синхронного разреза других площадей Восточного Абшерона. Месторождение Зых входит в Восточно-Апшеронский синклинорий и охватывает Зыхский участок Карачухур-Зыхской антиклинали.
а)
Рассматриваемое месторождение приурочено к асимметричной брахиантиклинальной складке меридионального простирания с крутым (45...500) восточным и сравнительно более пологим (22...250) западными крыльями (рис. 2). Складка резко асимметрична в продольном сечении [1]. Короткая северная периклиналь ограничивается хорошо выраженной седловиной, а южная, сильно вытянутая периклиналь складки прослеживается до о. Песчаный (рис. 2). Складка осложнена рядом продольных нарушений, которые затухают в пределах верхнего отдела ПТ (за исключением двух основных сбросов). Нарушения имеют небольшие амплитуды и амплитуда их редко превышает 40...50 м. Поперечные нарушения незначительны, амплитуда их 10...20 м. Они приурочены к сводовой части складки. Южное погружение складки в районе Зыхского озера осложнено ископаемым грязевым вулканом [1, 2].
б)
а)
Рис. 1. Расположение месторождения Зых (а) и его космический снимок (б)
Рис. 2. Геологический разрез (а) месторождения Зых и временной разрез (б)
Площадь Зых неоднократно изучалась различными модификациями сейсмической разведки [1, 4]. В 1937...1941, 1945...1949, 1949...1951 и 1955 годах выполнены сейсмо-разведочные работы методом отраженных волн (МОВ), в 1977...1978 гг. - методом общей глубинной точки (МОГТ). Отдельные профили отработаны в 1993 и 1999 годах. В 1996 и 2003...2004 годах проводились детальные сейсморазведочные работы методом ОГТ. С целью изучения скоростной модели среды выполнены сейсмокаротаж-ные исследования СК в скв. 60 площади Карачухур, в 1949...1950 годах - в скв. 156, 191 площади Зых, в 2012 году в скв. 1856 и в 2014 году в скв. 2014 площади Говсаны проводилось ВСП (вертикальное сейсмическое профилирование). В 2011...12 гг. проводилась сейсморазведка 3D, охватившая еще и соседнее с востока месторождение Говсан, материалы которых легли в основу данных исследований [7].
Площадь Зых введена в разработку с 1935 г., когда был получен фонтан нефти из IX горизонта продуктивной толщи (ПТ). В 1938...1940 гг. выявлены залежи нефти в свитах нижнего отдела ПТ (НКП-надкирмакинская песчанистая, КС-кирмакинская свита, ПК-подкирмакинская и КаС-калинская свита). Отложения нижнего отдела продуктивной толщи: КаС (калинская свита), характеризуется чередованием глинистых и песчаных прослоев [2]; ПК (подкирмакинская свита) представлена толстослоистыми песчаными породами с тонкими линзовидными прослоями глин.
Отложения верхнего отдела продуктивной толщи: пласты VIII и IX (балаханская свита), представлены чередованием песчаных, алевритовых, глинистых пород.
Основными целями настоящих исследований является: построение структурного каркаса по пластам балахан-ской свиты (VIII и IX) и пластам ПК и КаС, который будут служить основой построения 3D-геологических моделей месторождения. Исследования включают: а -построение структурных карт на основе карт изохрон и разбивок в скважинах; б -построение модели разломов; в - составление структурных карт в дизъюнктивном варианте.
Основная часть. Построение структурных карт осуществлялся на основе карт изохрон и разбивок в скважинах. Для построения структурных карт были использованы карты изохрон, полученных ООО «ГЕОЛОГИЯ РЕЗЕРВУАРА», «Отчет по проведению обработки и комплексной интерпретации сейсморазведки МОГТ-3D и скважинных данных с целью до разведки контрактной площади месторождений Зых-Говсаны в Абшеронском нефтегазоносном районе», 2012 г. [1, 4].
Для построения структурных карт по пластам балаханской свиты (VIII и IX) были использованы карты изохрон по сейсмическим горизонтам БС-VIII и БС-IX. По результатам выполненной нами корреляции данных в скважинах были построены карты средних скоростей (интерполяция по скважинам), после чего были получены структурные карты по кровле пластов VIII и IX. В качестве подошвы моделирования пласта IX принята поверхность по кровле пласта SP (свита «перерыва»), которая получена через карту толщин, построенную по результатам интерполяции значений толщин пласта IX в скважинах [5].
Структурная карта по кровле подкирма-кинской свиты (пласт ПК) построена от сейсмического горизонта II через карту средних скоростей, построенную на основе раз б и вок в скважинах [3]. Структурная карта по кровле калинской свиты (пласт КаС) пол учена через карту общих толщин, построенную по результатам интерполяции значений толщин пласта ПК в скважинах. Подошва моделирования пласта КаС (ввиду недостаточной информации по скважинным данным) принята условно, при постоянной мощности отложений калинской свиты рав ной 100 м.
Имея ввиду, что карты изохрон были построены на производстве в пликативном варианте, то, на этом этапе, и структурные карты получаем без учета тектонических нарушений. Разломы будут встроены в структурные карты на следующем этапе, п о сле п о строения модели разломов [10, 11]. На ри с. 3 приведены структурные карты по СГ-II, БС-VIII и БС-IX.
Рис. 3. Структурные карты: а - по кровле пласта ПК; б - по кровле пласта IX; в - по кровле пласта VIII
Предметом объемного моделирования геологического строения являются: форма залежей (структурный план, морфология поверхности флюидных контактов), внутренняя структура резервуара и залежей (пространственное взаимоотношение
коллекторов и неколлекторов, распределение в разрезе различных по типу флюидов), фильтрационно-емкостные свойства пород-коллекторов и их нефте-, газонасыщенность.
Настоящие исследования охватывают вопросы, касающиеся структурного плана целевых горизонтов, т.е. форм залежей. Очевидно, что только пликативным вариантом строения невозможно составлять адекватную объемный модель геологического строения. Чтобы получить наиболее точную модель строения необходимо построить модели разломов, для чего использовался
а)
программный пакет PETREL.10.2 [8]. Объектами исследования являлись пластово-св одовые, литологически и тектонически экранированные залежи продуктивных пластов нижнего плиоцена.
Для выделения тектонических нарушений в волновом поле применялись кубы когерентности и погружения [9, 12]. Куб когерентности рассчитывается на основе исходного куба. Куб когерентности является кубом коэффициентов корреляции между трассами сейсмической записи. Та часть се й смического поля, где нет осложнений и разломов будет равна ~1, а места разломов -меньше единицы [13,14, 15]. Критерием наличия дизъюнктивных нарушений служит резкое изменение параметра когерентности (рис. 4).
Рис. 4. Карты атрибута когерентности в интервале пластов VIII (а) и IX (б)
Кубы погружений и азимутов рассчитываются на основе мгновенной фазы и мгновенной частоты. И также является атрибутом для выделения тектонических нарушений. Критерием достоверности картирования дизъюнктивных нарушений служит резкое изменение градиента в кубе погружений (рис. 5). Качество сейсмической записи не позволяет детально подойти к вопросу выделения разломов, поэтому, в рамках настоящей работы, модель разломов для пластов балаханской свиты (VIII и IX) принята единой, ровно как и для пластов нижнего отдела продуктивной толщи (подкирмакинской и калинской свит). Для каждого из моделируемых пластов были
построены карты атрибутов когерентности и погружений. Расчет производился по вре менному кубу в интервале соответствующего пласта. На основе анализа полученных карт атрибутов, а также принимая во внимание геофизические и промысловые характеристики скважин, попадающих в один тектонический блок, были откартиро-ваны тектонические нарушения. На рис. 6 приведена карта атрибутов когерентности моделируемого пласта, т.е. для интервала пластов ПК и КаС. А на рис. 7 приведена принятая модель разломов для пластов балаханской свиты. На рис. 8 приведена принятая модель разломов для пластов подкирмакинской и калинской свит.
- '' .,'Г 1 ■ >
ч.ЧЛ-
Ю
ш
Рис. 5. Карты атрибута DEP (погружений) в интервале пластов VIII (а) и IX (б)
а)
<УШ)1Л
Рис. 6. Карты атрибута когерентности в интервале пластов ПК (а) и КаС (б)
Принятая нами модель разломов
Для пластов VIII и IX поверхности разломов приняты
вертикальными
щшш
JS^umSS. ш »
Рис. 7. Принятая модель разломов для пластов балаханской свиты (VIII и IX) на основе атрибута DEP
Принятая нами модель разломов
Для пластов ПК и КаС поверхности разломов
приняты вертикальным»
жк
Жштл
Рис. 8. Принятая модель разломов для пластов подкирмакинской и калинской свит (ПК и КаС)
на основе когерентности
Данные, полученные на основе атрибутного анализа (БЕР и когерентности), использованы для получения модели разломов и уточнения геологического строения месторождения Зых и для построения структурных карт в дизъюнктивном варианте целевых горизонтов. На основе построенной модели разломов и пликативных структурных карт построены структурные карты в дизъюнктивном варианте и объемные
модели тектонических блоков (сегменты 3Б-модели). Модели сегментов для нижнего и верхнего отделов продуктивной толщи приведены на рис. 9.
Полученные по атрибутному анализу мо де л и сегментов для отдельных отделов продуктивной толщи нижнего плиоцена получили свое отражения в структурных картах в дизъюнктивном варианте целевых пластов (рис. 10 и 11).
а)
а)
а)
Рис. 9. Выделенные тектонические блоки для пластов VIII и IX (а), ПК и КаС (б)
б)
Рис. 10. Структурные карты по кровле пласта VIII (а) и по кровле пласта IX (б)
б)
Рис. 11. Структурные карты по кровле пласта ПК (а) и по кровле пласта КаС (б)
Заключение. 1. Структурные карты, полученные нами, по данным сейсморазведки 3D на данный момент наиболее достоверно и точно отображают строение месторождения Зых.
2. На основе построенных структурных 2D-поверхностей по стратиграфическим отметкам по кровле и подошве пластов, а также построенных моделей разломов, будут построены дискретные трехмерные сетки, т.е. они будут использованы для построения геологической (литологической, петрофизи-ческой) модели месторождения Зых.
Литература
1. Ализаде А.А., Ахмедов Г.А., Ахмедов А.М. и др. Геология нефтяных и газовых месторождений Азербайджана. M.: Недра, 1966. С. 311-313.
2. Ализаде, А.А. и др. Продуктивная толща Азербайджана. М.: Недра, 2018. Том I, II.
3. Ампилов Ю.П. Сейсмическая интерпретация: опыт и проблемы. М.: Геоинформ-марк, 2004. 277 с.
4. Ахмедов Т.Р. О геологической эффективности сейсморазведки при изучении не антиклинальных ловушек Азербайджана разного типа // Известия УГГУ. 2016. С. 41-45.
5. Ахмедов, Т.Р. и др. Роль боковых волн в формировании волнового поля скважинной сейсморазведки на примере площади Говсан // Азербайджанское нефтяное хозяйство. 2013. № 2. С. 9-14.
6. Барышев Ю.А. Прогноз продуктивности терригенных коллекторов по динамическим параметрам отраженных волн на Верхнечонской площади // Геофизика. 2001. №2. С. 27-32.
7. Воскресенский Ю.Н. Состояние и перспективы развития методов анализа амплитуд сейсмических отражений для прогнозирования залежей углеводородов // Геология, методы поисков, разведки и оценки месторождений топливно-энергетического сырья. М . : Геоинформцентр, 2002. Вып. 4-5. 77 с.
8. Кириллов А.С., Закревский K.E. Практикум по сейсмической интерпретации в Петрел. M.: Издательство МАИ-Принт, 2014. 288 с.
9. Кондратьев И.К., Бондаренко М.Т., Каменев С.П. Динамическая интерпретация данных сейсморазведки при решении задач нефтегазовой геологии // Геофизика. 1996. №5-6. С. 41-47.
10. Кондратьев О.К. Автоматизированные системы оценки качества сейсмограмм и временных сейсмических разрезов ОГТ // Геофизика. Специальный выпуск. 2002. С. 3-12.
11. Овчаренко А.В., Сафонов А.С., Шлезингер, А.Е. и др. Методические приемы интерпретации геофизических материалов при поисках, разведке и освоении месторождений углеводородов. М.: Научный мир, 2002. С. 43-57.
12. Урупов А.К. Основы трехмерной сейсморазведки. M.: Нефть и газ, РГУНГ. 2004. 584 с.
13. Chopra S., Marfurt K.J. Seismic attributes-A historical perspective // Geophysics. 2005. Т. 70. № 5. С. 3-28.
14. Fozao K.F., Fotso K., Djieto-Lordon A., Mbeleg M. Hydrocarbon inventory of the eastern part of the Rio Del Rey Basin using seismic attributes // Journal of petroleum exploration and production technology. 2018. Т.8. №3. С. 655-665.
15. Teyler M.N. Application of 3 D seismic attributes analysis workflow: A case study from NESS country KANSAS, USA, B.S APPL 3D SEISMIC ATTR. 2013. 30 с.
Контактные данные:
Ахмедов ТофикРашид оглы, эл. почта: akhmedov.tofik@bk.ru
© Ахмедов Т.Р., 2020
DETAILED STUDY OFZYKH FIELD SETTING BY USE OF ATTRIBUTE ANALYSIS OF 3D SEISMIC SURVEY DATA
T.R. Ahmеdov*
Azerbaijan State Oil and Industry University, Baku, Azerbaijan
* E-mail: akhmedov. tofik@bk.ru
Abstract. The paper is devoted to kinematic interpretation of 3D data acquired in Zykh area. The paper expounds in brief the data on study area, its geographical position, deposits recovered by deep drilling, setting of Zykh fold and the history of geological and geophysical studies covering the area. It has been noted that the major aim of these studies consists in outlining the structure for Balakhany suite layers (VIII and IX) and Qirmakyalty and Qala suite layers, which will be taken as a base for 3D geological models design. The subject of volumetric modeling of geological setting are the following: the form of deposit (structural plane, morphology of fluid contacts surface), inner structure of reservoir and deposits (spatial interrelation of reservoirs and non-reservoirs, distribution of various fluid types across the section), filtration and capacitance parameters of reservoir-rocks and their oil- and gas saturation.
The present study covers the questions regarding the structural plane of target horizons, i.e. the deposit forms. It is obvious that it is not possible to design the adequate volumetric model of geological setting by use of plicatives model only. The structural maps have been drawn for layers of Balakhany suite (VIII and IX), as well as for the top of Qirmakyalty suite (Qirmakyalty layer).
Evidently, it is not possible to design an adequate volumetric model of geological setting by use of plicatives version only. To gain more accurate model it is required to design models of faults and for this PETREL 10.2 software package has been applied. To design structural maps we have used maps drawn in plicative version and 3D models of faults drawn on the basis of attribute analysis of seismic data. In general, two cubes of seismic wavefield attributes have been applied -coherency and subsidence cubes.
It has been noted that the data acquired on the basis of attribute analysis (DEP and coherency) were used for design of faulting model and for acquisition of more accurate geological setting of Zykh field and for design of structural maps in disjunctive version of target horizons. Designed model of faulting and structural maps for plicatives were used to draw structural maps in disjunctive version and design of volumetric models of tectonic blocks (3D model segments).
Finally, models of segments acquired by attribute analysis for various parts of Productive Series of the Lower Pliocene have been reflected in structural maps in disjunctive version of the target layers, which are the final results of present study.
Keywords: 3D seismic survey, Zykh field, Balakhany suite (VIII and IX), Qyrmakyalty suite, attribute analysis.
References
1. Alizade, A.A., Ahmedov, G.A., Ahmedov, A.M. i dr. Geologiya neftyanyh i gazovyh mestorozhdenij Azerbajdzhana [Geology of oil and gas fields in Azerbaijan]. M.: Nedra, 1966. Pp. 311-313. (rus)
2. Alizade, A.A. i dr. Produktivnaya tolshcha Azerbajdzhana [Productive strata of Azerbaijan]. M.: Nedra, 2018. Tom I, II. (rus)
3. Ampilov, Yu.P. Sejsmicheskaya interpretaciya: opyt i problem [Seismic interpretation: experience and problems]. M.: Geoinform-mark, 2004. 277 p. (rus)
4. Ahmedov, T.R. O geologicheskoj effek-tivnosti sejsmorazvedki pri izuchenii ne anti-klinal'nyh lovushek Azerbajdzhana raznogo tipa [On the geological efficiency of seismic prospecting in the study of different types of non-anticlinal traps in Azerbaijan] // Izvestiya UGGU. 2016. Pp. 41-45. (rus)
5. Ahmedov, T.R. i dr. Rol' bokovyh voln v formirovanii volnovogo polya skvazhinnoj sejsmorazvedki na primere ploshchadi Govsan [The role of lateral waves in the formation of the wave field of borehole seismic survey on the example of the Hovsan area] // Azerbajdzhanskoe neftyanoe hozyajstvo. 2013. № 2. Pp. 9-14. (rus)
6. Baryshev, Yu.A. Prognoz produktivnosti terrigennyh kollektorov po dinamicheskim parametram otrazhennyh voln na Verhnechonskoj ploshchadi [Forecast of the productivity of terrigenous reservoirs by the dynamic parameters of reflected waves in the Verkhnechonskaya area] // Geofizika. 2001. №2. Pp. 27-32. (rus)
7. Voskresenskij, Yu.N. Sostoyanie i perspektivy razvitiya metodov analiza amplitud sejsmicheskih otrazhenij dlya prognozirovaniya zalezhej uglevodorodov [State and prospects for the development of methods for analyzing the amplitudes of seismic reflections for predicting hydrocarbon deposits] // Geologiya, metody poiskov, razvedki i ocenki mestorozhdenij toplivno-energeticheskogo syr'ya. M.: Geoinformcentr, 2002. Vyp. 4-5. 77 p. (rus)
8. Kirillov, A.S., Zakrevskij, K.E. Praktikum po sejsmicheskoj interpretacii v Petrel [State and prospects for the development of methods for analyzing the amplitudes of seismic reflections for predicting hydrocarbon deposits]. M.: Izdatel'stvo MAI-Print, 2014. 288 p. (rus)
9. Kondrat'ev, I.K., Bondarenko, M.T., Kamenev, S.P. Dinamicheskaya interpretaciya dannyh sejsmorazvedki pri reshenii zadach neftegazovoj geologii [Dynamic interpretation of seismic data in solving problems of oil and gas geology] // Geofizika. 1996. №5-6. Pp. 41 -47. (rus)
10. Kondrat'ev, O.K. Avtomatizirovannye si stemy ocenki kachestva sejsmogramm i vremennyh sejsmicheskih razrezov OGT [utomated systems for assessing the quality of seismograms and temporary seismic sections of CDP] // Geofizika. Special'nyj vypusk. 2002. Pp. 3-12. (rus)
11. Ovcharenko, A.V., Safonov, A.S., Shl ezinger, A.E. i dr. Metodicheskie priemy interpretacii geofizicheskih materialov pri poiskah, razvedke i osvoenii mestorozhdenij uglevodorodov [Methodological methods of interpretation of geophysical materials in the search, exploration and development of hydrocarbon deposits]. M.: Nauchnyj mir, 2002. Pp. 43-57. (rus)
12. Urupov, A.K. Osnovy trekhmernoj sejsmorazvedki [Fundamentals of 3D seismic exp l oration]. M.: Neft' i gaz, RGUNG. 2004. 584 p. (rus)
13. Chopra, S., Marfurt, K.J. Seismic attrib-utes-A historical perspective // Geophysics. 2005. Vol. 70. Issue 5. 30 p.
14. Fozao, K.F., Fotso, L., Djieto-Lordon A., Mbeleg M. Hydrocarbon inventory of the eastern part of the Rio Del Rey Basin using seismic attributes // Journal of petroleum exploration and production technology. 2018. Vol. 8. Issue 3. Pp. 655-665.
15. Teyler, M.N. Application of 3 D seismic attributes analysis workflow: A case study from NESS country KANSAS, USA, B.S APPL 3D SEISMIC ATTR. 2013. 30 p. Published 2013.
Contacts:
Tofiq R. Ahmedov, akhmedov.tofik@bk.ru
© Ahmedov, T.R., 2020
Ахмедов Т.Р. Уточнение строения месторождения Зых применением атрибутного анализа данных сейсморазведки 3D // Вектор ГеоНаук. 2020. Т.3. №3. С. 5-14. DOI: 10.24411/2619-0761-2020-10025.
Ahmеdov, T.R., 2020. Detailed study of Zykh field setting by use of attribute analysis of 3D seismic survey data. Vector of Geosciences. 3(3). Pp. 5-14. DOI: 10.24411/2619-0761-2020-10025.