RU2432450C2 - Procedure for development of non-uniform massive or multipay gas-oil or oil-gas-condensate field - Google Patents

Procedure for development of non-uniform massive or multipay gas-oil or oil-gas-condensate field Download PDF

Info

Publication number
RU2432450C2
RU2432450C2 RU2009143585/03A RU2009143585A RU2432450C2 RU 2432450 C2 RU2432450 C2 RU 2432450C2 RU 2009143585/03 A RU2009143585/03 A RU 2009143585/03A RU 2009143585 A RU2009143585 A RU 2009143585A RU 2432450 C2 RU2432450 C2 RU 2432450C2
Authority
RU
Russia
Prior art keywords
classes
oil
parameters
gas
class
Prior art date
Application number
RU2009143585/03A
Other languages
Russian (ru)
Other versions
RU2009143585A (en
Inventor
Расим Нагимович Дияшев (RU)
Расим Нагимович Дияшев
Ринат Гатинович Харисов (RU)
Ринат Гатинович Харисов
Виктор Николаевич Рябченко (RU)
Виктор Николаевич Рябченко
Анатолий Александрович Савельев (RU)
Анатолий Александрович Савельев
Николай Александрович Зощенко (RU)
Николай Александрович Зощенко
Original Assignee
Общество с Ограниченной Ответственностью "ТНГ-Групп"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Общество с Ограниченной Ответственностью "ТНГ-Групп" filed Critical Общество с Ограниченной Ответственностью "ТНГ-Групп"
Priority to RU2009143585/03A priority Critical patent/RU2432450C2/en
Publication of RU2009143585A publication Critical patent/RU2009143585A/en
Application granted granted Critical
Publication of RU2432450C2 publication Critical patent/RU2432450C2/en

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

FIELD: oil and gas production.
SUBSTANCE: procedure consists in drilling boreholes, in recovery of core, in geophysical survey, in testing and in hydro-dynamic survey of separate intervals by thickness of cross-section. According to the inventions database is created on results of measurements of various geo-physical parametres accepted for each region for each 0.05-0.15 m of pay thickness. Parameters with asymmetrical logging normal distribution are converted to normal distribution; all parametres are standardised by subtraction of average and by division to mean square deviation. There is made hierarchic procedure of classes arrangement consisting in successive interlayers uniting in such a way that obtained classes have minimal parametres spread around centre - dispersion. Ranges of changes and average values of these parametres for each singled out class of pays are calculated. There is determined nature of lithology, filtration-capacitive properties and fluid saturation of each class on qualitative level with consideration of fundamental links of these parametres and mineral element composition of rock. There are additionally used results of core analysis, of tests and hydro-dynamic survey in intervals of separate classes-pays; the said characteristics are determined at quantitative level. Singled out classes between boreholes are correlated. There is built geological-hydrodynamic model of producing deposit on base of which there is evaluated efficiency of system of development for various versions of inclusion of singled out classes into operational objects. There are made calculations of process indices of development with determination of levels of hydrocarbon extractions under various modes and there are evaluated prognosis coefficients of oil extraction for separate operational objects and field in whole. The best version is realised in practice.
EFFECT: differentiation of non-uniform massive or multi-pay producing cross section represented with non-collectors and collectors with various lithology and fluid saturation into typical intervals - classes-pays in each borehole using data of all performed complex of survey.
2 cl, 1 ex, 8 tbl, 10 dwg

Description

Изобретение относится к нефтегазовой промышленности и может найти применение при освоении и разработке неоднородных массивных или многопластовых газонефтяных или нефтегазоконденсатных месторождений.The invention relates to the oil and gas industry and may find application in the development and development of heterogeneous massive or multi-layer gas-oil or oil-gas condensate fields.

Известен способ разработки неоднородных массивных или многопластовых газонефтяных или нефтегазоконденсатных месторождений путем создания самостоятельных эксплуатационных объектов, включающих пласты-коллекторы с «близкими» фильтрационно-емкостными свойствами (ФЕС) и насыщенные углеводородами с «близкими» характеристиками, используя данные геофизических исследований (ГИС), геолого-технологических исследований (ГТИ), анализа керна, пластовых флюидов, результатов испытаний пластов и их гидродинамических исследований (ГДИ) (Патент РФ №2346148, опублик. 10.02.2009).There is a method of developing heterogeneous massive or multi-layer gas-oil or oil and gas condensate fields by creating independent production facilities, including reservoir reservoirs with "close" reservoir properties and saturated with hydrocarbons with "close" characteristics, using data from geophysical surveys (GIS), geological -technological studies (GTI), core analysis, formation fluids, reservoir test results and their hydrodynamic studies (GDI) (RF Patent No. 2346148, published. 02.10.2009).

Известный способ не позволяет формализовать процесс расчленения разреза и объективно выделять в разрезе месторождения пласты с «близкими» ФЕС и насыщенными углеводородами с «близкими» характеристиками, ибо:The known method does not allow to formalize the process of partitioning the section and objectively distinguish in the section of the deposit formations with "close" FES and saturated hydrocarbons with "close" characteristics, because:

1. Определение с «близкими» свойствами носит качественный характер;1. A definition with “close” properties is of a qualitative nature;

2. Часто те или иные геофизические методы исследования оказываются неэффективными по геологическим и технологическим причинам (плотные породы с высоким сопротивлением, трещиноватые и кавернозные породы, сильное поглощение промывочной жидкости и т.д.);2. Often, certain geophysical research methods are ineffective for geological and technological reasons (dense rocks with high resistance, fractured and cavernous rocks, strong absorption of flushing fluid, etc.);

3. Керн отбирается не из каждой скважины, часто недостаточен его вынос;3. The core is not selected from each well; its removal is often insufficient;

4. Испытаниями и гидродинамическими исследованиями охватываются в разных случаях разные интервалы и толщины пласта и оцениваются средние фильтрационные характеристики всего интервала испытаний (исследований);4. Tests and hydrodynamic studies cover in different cases different intervals and thicknesses of the reservoir and average filtration characteristics of the entire test (research) interval are evaluated;

5. Результаты решения задачи расчленения разреза зависят от профессионализма исполнителей.5. The results of solving the problem of partitioning a section depend on the professionalism of the performers.

Поэтому на начальном этапе освоения месторождения, когда принимаются базовые решения, определяющие принципы разработки на весь срок, такие как выделение эксплуатационных объектов и создание систем их разработки, информационное обеспечение оказывается неадекватным сложности неоднородных массивных или многопластовых месторождений. Следствием являются неполный охват процессом вытеснения по толщине и объему эксплуатационных объектов и всего месторождения, низкие коэффициенты извлечения углеводородов.Therefore, at the initial stage of field development, when basic decisions are made that determine the principles of development for the entire period, such as the allocation of production facilities and the creation of systems for their development, information support is inadequate to the complexity of heterogeneous massive or multi-layer deposits. The consequence is incomplete coverage of the displacement process in terms of thickness and volume of production facilities and the entire field, low hydrocarbon recovery rates.

Наиболее близким к предложенному изобретению по технической сущности является способ, который включает бурение скважин с отбором керна и проведение стандартного для данного геологического разреза комплекса ГИС, анализ керна и создание общей по месторождению базы данных по пористости и проницаемости по керну, выделение геологических объектов по индексу перколяции с высокой степенью корреляции для групп коллекторов, имеющих некоторые особенности строения структуры порового пространства, и определение проницаемости коллекторов, используя удельные сопротивления по диаграммам каротажа сопротивлений (Патент РФ №2211329, опублик. 27.08.2003 г. - прототип).Closest to the proposed invention in technical essence is a method that includes drilling wells with coring and conducting a standard GIS complex for a given geological section, analyzing the core and creating a database for the porosity and permeability for the core of the field, identifying geological objects by the percolation index with a high degree of correlation for groups of reservoirs having some structural features of the pore space structure, and determination of reservoir permeability, using zuya resistivity diagrams resistivity logging (RF patent №2211329, publish 27.08.2003, the -. the prototype).

Известный способ не позволяет дифференцировать разрез каждой скважины, пробуренной на месторождении с неоднородными массивными или многопластовыми продуктивными отложениями, а служит лишь для повышения точности определения одного параметра - коэффициента проницаемости по результатам каротажа, используя петрофизические зависимости, построенные с предварительным выделением групп пластов, отличающихся по структуре порового пространства. Известный способ не может быть применен для коллекторов трещиноватых и кавернозных, т.к. для них по керну параметры пористости и проницаемости адекватно не характеризуют емкостные и фильтрационные характеристики всего коллектора. Кроме того, известный способ не ориентирован на использование всего потенциала имеющихся в настоящее время методов ГИС.The known method does not allow to differentiate the section of each well drilled in the field with heterogeneous massive or multilayer productive deposits, but serves only to increase the accuracy of determining one parameter - permeability coefficient from the logging results using petrophysical dependencies constructed with preliminary identification of reservoir groups that differ in structure pore space. The known method cannot be applied to fractured and cavernous reservoirs, because for them, core parameters of porosity and permeability do not adequately characterize the capacitive and filtration characteristics of the entire reservoir. In addition, the known method is not focused on using the full potential of the currently available GIS methods.

В предложенном способе решается задача дифференциации неоднородного массивного или многопластового продуктивного разреза, представленного неколлекторами и коллекторами с различной литологией и флюидонасыщенностью на характерные интервалы-классы пластов в каждой скважине, используя данные всего выполненного комплекса ГИС.The proposed method solves the problem of differentiating a heterogeneous massive or multi-layer productive section, represented by non-reservoirs and reservoirs with different lithology and fluid saturation into characteristic intervals-classes of formations in each well, using data from the entire well logging system performed.

Задача решается тем, что в способе разработки неоднородных массивных и многопластовых газонефтяных или нефтегазоконденсатных месторождений, включающем бурение скважин, проведение комплекса геофизических исследований, согласно изобретению создают базу данных по результатам измерений различных, принятых для каждого региона, геофизических параметров, например электрических, радиоактивных, волновых, газового каротажа для каждой 0,05-0,15 м толщины пласта, параметры с ассиметричным (логнормальным) распределением преобразуют к нормальному распределению, стандартизуют все параметры вычитанием среднего и делением на среднеквадратичное отклонение, после чего выполняют иерархическую процедуру построения классов, состоящую в последовательном объединении прослоев таким образом, что полученные классы имеют минимальный разброс параметров вокруг центра (дисперсию), выделяя по резкому росту дисперсии естественное для данного разреза число классов, которое в зависимости от неоднородности, общей толщины и многопластовости разреза и количества фаз флюида может составить до 8-10 с отличающимися для них геофизическими параметрами, вычисляют диапазоны изменения и средние значения этих параметров для каждого выделенного класса пластов и, учитывая фундаментальные связи этих параметров с минеральным и элементным составом пород, определяют характер литологии, фильтрационно-емкостные свойства и флюидонасыщение на качественном уровне, используя дополнительные данные по керну, испытаниям и гидродинамическим исследованиям, выполненным в интервалах тех или иных выделенных классов - на количественном уровне, проводят корреляцию выделенных классов между скважинами, строят геолого-гидродинамическую модель продуктивной залежи (месторождения), на которой оценивают эффективность системы разработки при различных вариантах включения выделенных классов в эксплуатационные объекты расчетами технологических показателей разработки с определением уровней добычи углеводородов на режимах истощения и закачки вытесняющих агентов с оценкой прогнозных коэффициентов нефтеизвлечения отдельных эксплуатационных объектов и месторождения в целом и лучший вариант реализуют на практике.The problem is solved in that in a method for developing heterogeneous massive and multi-layer gas-oil or oil-gas condensate fields, including drilling wells, conducting a complex of geophysical surveys, according to the invention, a database is created by measuring various geophysical parameters adopted for each region, for example, electrical, radioactive, wave , gas logging for each 0.05-0.15 m of formation thickness, parameters with asymmetric (lognormal) distribution are converted to normal Allocation of the parameters is standardized by subtracting the mean and dividing by the standard deviation, after which they perform a hierarchical class construction procedure consisting in sequentially combining interlayers in such a way that the resulting classes have a minimum dispersion of parameters around the center (variance), highlighting the sharpness of dispersion that is natural for a given section number of classes, which, depending on the heterogeneity, total thickness and multi-layer nature of the section and the number of fluid phases, can be up to 8-10 s geophysical parameters differing for them, calculate the ranges of variation and average values of these parameters for each selected class of formations and, taking into account the fundamental relationships of these parameters with the mineral and elemental composition of rocks, determine the nature of lithology, filtration-capacitive properties and fluid saturation at a qualitative level using additional data according to the core, tests and hydrodynamic studies performed in the intervals of certain distinguished classes - at a quantitative level, conduct to correlation of distinguished classes between wells, a geological and hydrodynamic model of a productive reservoir (field) is built, on which the effectiveness of the development system is assessed for various options for incorporating the distinguished classes into production facilities by calculating technological development indicators with determining hydrocarbon production levels at depletion and injection modes of displacing agents with an estimate forecast oil recovery coefficients of individual production facilities and the field as a whole and the best option p Alize in practice.

Сущность изобретенияSUMMARY OF THE INVENTION

В известных технических решениях по расчленению разреза месторождения в качестве базовой информации используются результаты анализа керна, как правило, отбираемого при бурении первых разведочных скважин. Объем отбираемого керна практически всегда ограничен не только по числу скважин, но и по разрезу из-за плохого его выноса. Результаты определения ФЕС пластов по керну адекватно не отражают характеристики пород в пластовых условиях, особенно при наличии трещиноватости, т.к. происходит изменение его свойств при отборе, выносе на поверхность в связи с изменением термобарических условий, а также подготовке образцов к анализам. Однако полученные на таком искаженном материале данные далее используются для построения петрофизических зависимостей «керн-керн» и «керн-ГИС» и последующей интерпретации по ним данных геофизических исследований скважин (ГИС) всего месторождения или даже группы месторождений. Построенные по такой информационной базе зависимости «пористость-проницаемость» для отдельных групп коллекторов на месторождении по известному решению несут в себе указанные выше погрешности и не позволяют получать адекватные величины параметров для пластов конкретных скважин и обоснованно объединять их в эксплуатационные объекты для самостоятельной разработки, что приводит к уменьшению охвата пластов процессом вытеснения углеводородов по толщине и площади и снижению коэффициента извлечения углеводородов.In well-known technical solutions for the exploration of a section of a field, the core analysis results, as a rule, taken during drilling of the first exploratory wells, are used as basic information. The volume of core samples is almost always limited not only by the number of wells, but also by section because of its poor removal. The results of the determination of reservoir properties by core do not adequately reflect the characteristics of the rocks in the reservoir conditions, especially in the presence of fracturing, because there is a change in its properties during sampling, removal to the surface in connection with a change in thermobaric conditions, as well as the preparation of samples for analysis. However, the data obtained on such distorted material are then used to construct the petrophysical dependences of “core-core” and “core-GIS” and subsequent interpretation of the data of geophysical studies of wells (GIS) of the entire field or even a group of fields. According to a well-known solution, the porosity-permeability dependences constructed for such individual information bases for individual groups of reservoirs contain the above errors and do not allow obtaining adequate parameter values for specific wells and reasonably combine them into production facilities for independent development, which leads to reduce the coverage of formations by the process of displacement of hydrocarbons in thickness and area and to reduce the coefficient of hydrocarbon recovery.

В известных технических решениях отсутствует и формализованное комплексное использование всех данных измерений геофизическими методами в скважине. Между тем прямые отсчеты показаний электрических, радиоактивных, ядерно-физических, геолого-технологических измерений несут в себе многопараметрическую информацию об особенностях строения коллекторов, их фильтрационно-емкостных свойствах и характере флюидонасыщения.In the known technical solutions, there is no formalized integrated use of all measurement data by geophysical methods in the well. Meanwhile, direct readings of the readings of electrical, radioactive, nuclear-physical, geological and technological measurements carry multiparameter information about the structural features of the reservoirs, their filtration-capacitive properties and the nature of fluid saturation.

В предложенном способе комплексно анализируются все параметры геофизических измерений в скважине по всему продуктивному разрезу месторождения с шагом 0,05-0,15 м с выделением интервалов-классов, отличающихся по их ФЕС и характеру флюидонасыщения. Последующее использование этих результатов при выделении самостоятельных эксплуатационных объектов и их разработке наиболее оптимальными для характеристик каждого из этих эксплуатационных объектов методами воздействия позволяет увеличить коэффициент извлечения углеводородов по сравнению с известными техническими решениями. Эта цель достигается тем, что в отличие от известного способа разработки неоднородных массивных или многопластовых газонефтяных и нефтегазоконденсатных месторождений путем бурения скважин, проведения комплекса геофизических исследований, выборочного отбора керна, испытания и гидродинамических исследований отдельных интервалов по толщине разреза скважин и выделения эксплуатационных объектов с «близкими» характеристиками, создают базу данных параметров всех выполненных геофизических измерений в скважине с их количественными значениями с шагом 0,05-0,15 м по разрезу. Этими параметрами могут быть результаты электрических измерений (кажущееся сопротивление - КС, естественный потенциал - ПС, индукционное сопротивление - ИК и др.), радиоактивных (гамма-каротаж - ГК, нейтронный гамма-каротаж - НГК, спектральный гамма-каротаж - ГКС, содержание радиоактивных элементов - калий, торий, уран и др.), волновых (интервальные времена пробега продольной и поперечной волн - DTP и DTS и волны Лэмба-Стоунли - DTST, коэффициенты затухания продольной, поперечной и волны Лэмба-Стоунли - ALPHAP, ALPHAS, ALPHAST и др.), газового каротажа (содержание СH4, C2H6, С3Н8, С4Н10, C5H12, С6Н14 и суммарное газосодержание), а также современных ядерно-физических (импульсный нейтронный гамма-каротаж - ИНГК, углеродный - кислородный каротаж - с/о и др.). Распределение каждого параметра проверяют на нормальность и те из них, которые имеют ассиметричное распределение, преобразуют к нормальному, например, путем логарифмирования, стандартизуют все параметры вычитанием среднего и делением на среднеквадратичное отклонение, после чего выполняют иерархическую процедуру построения классов, состоящую в последовательном объединении прослоев, начиная с минимального их числа, таким образом, что полученные классы имеют минимальный разброс параметров вокруг центра (дисперсию), выделяя по резкому росту дисперсии естественное для данного разреза число классов, которое в зависимости от неоднородности, общей толщины и многопластовости разреза и количества фаз флюида может составить до 8-10 с отличающимися для них геофизическими параметрами. Точность отнесения отдельных слоев с толщиной 0,05-0,15 м к тому или иному классу составляет до 95% и более. Вычисляют диапазоны изменения и средние значения этих параметров для каждого выделенного класса пластов. Путем сопоставительного анализа средних значений геофизических параметров для выделенных классов и, учитывая фундаментальные связи этих параметров с минеральным и элементным составом пород, характером флюидонасыщенности, определяют литологию, строение и ФЕС коллекторов, типы насыщающих пласты флюидов на качественном уровне, используя ограниченные данные по керну, испытаниям и гидродинамическим исследованиям, выполненным в интервалах тех или иных выделенных классов пластов - на количественном уровне. Выполняют корреляцию выделенных и охарактеризованных качественно и количественно пластов в виде классов между скважинами, что представляет надежную основу для построения 3-мерной геолого-гидродинамической модели месторождения, наполняя ее физико-химическими свойствами флюидов по данным PVT-исследований, на которой выполняют повариантные расчеты технологических показателей разработки месторождения при различных схемах формирования эксплуатационных объектов, т.е. при разных составах эксплуатационных объектов, включающих в разных сочетаниях те или иные классы пластов в разрезе месторождения. Число выделяемых при многомерной классификации фильтрационных интервалов - классов зависит от неоднородности, общей толщины и многопластовости разреза и фазового состояния пластового флюида, а количество эксплуатационных объектов - от ожидаемых коэффициентов нефтеизвлечения. Лучшие сочетания пластов в эксплуатационных объектах по их количеству, характеристикам, а также технологическим показателям - коэффициентам нефтеизвлечения - внедряют на практике.In the proposed method, all parameters of geophysical measurements in the well are comprehensively analyzed over the entire productive section of the field with a step of 0.05-0.15 m with the allocation of interval classes that differ in their physical properties and the nature of fluid saturation. The subsequent use of these results in the identification of independent production facilities and their development by the most optimal impact methods for the characteristics of each of these production facilities allows increasing the hydrocarbon recovery coefficient in comparison with the known technical solutions. This goal is achieved by the fact that, in contrast to the known method of developing heterogeneous massive or multi-layer gas and oil and oil and gas condensate fields by drilling wells, conducting a complex of geophysical surveys, selective coring, testing and hydrodynamic studies of individual intervals along the thickness of the section of the wells and identifying production facilities with "close »Characteristics, create a database of parameters of all performed geophysical measurements in the well with their quantitative values with increments of 0.05-0.15 m along the section. These parameters can be the results of electrical measurements (apparent resistance - CS, natural potential - PS, induction resistance - IR, etc.), radioactive (gamma-ray log - GC, neutron gamma-ray log - NGC, spectral gamma-ray log - GC, content radioactive elements - potassium, thorium, uranium, etc.), wave (interval travel times of longitudinal and transverse waves - DTP and DTS and Lamb-Stoneley waves - DTST, attenuation coefficients of longitudinal, transverse and Lamb-Stoneley waves - ALPHAP, ALPHAS, ALPHAST etc.), gas logging (content H 4, C 2 H 6, C 3 H 8, C 4 H 10, C 5 H 12, C 6 H 14 and the total gas content), as well as modern nuclear physics (pulsed neutron-gamma logging - INGK, the carbon - oxygen logging - s / o, etc.). The distribution of each parameter is checked for normality, and those that have an asymmetric distribution are converted to normal, for example, by logarithm, standardize all parameters by subtracting the mean and dividing by the standard deviation, and then perform a hierarchical class construction procedure consisting in sequentially combining layers, starting with their minimum number, so that the classes obtained have a minimum spread of parameters around the center (variance), highlighting in sharp Ost natural dispersion for a given number of classes of the section, which, depending on the inhomogeneity and the total thickness of the multilayered and the number section of the fluid phases may amount to 8-10 with differing for their geophysical parameters. The accuracy of classifying individual layers with a thickness of 0.05-0.15 m to one or another class is up to 95% or more. The ranges of variation and the average values of these parameters for each selected class of formations are calculated. A comparative analysis of the average values of geophysical parameters for the selected classes and, taking into account the fundamental relationships of these parameters with the mineral and elemental composition of rocks, the nature of fluid saturation, determine the lithology, structure and reservoir properties of reservoirs, types of fluid-saturated reservoirs at a qualitative level using limited core data, tests and hydrodynamic studies carried out in the intervals of certain distinguished classes of formations - at a quantitative level. Correlation of the identified and qualitatively and quantitatively stratified formations in the form of classes between the wells is performed, which provides a reliable basis for constructing a 3D geological and hydrodynamic model of the field, filling it with the physicochemical properties of the fluids according to PVT studies, on which variant calculations of technological parameters are performed field development under various schemes for the formation of production facilities, i.e. with different compositions of production facilities, including in different combinations these or those classes of formations in the context of the field. The number of filtration intervals - classes allocated during multidimensional classification - depends on the heterogeneity, total thickness and multilayerness of the section and the phase state of the formation fluid, and the number of production facilities depends on the expected oil recovery coefficients. The best combinations of formations in production facilities by their number, characteristics, as well as technological indicators - oil recovery coefficients - are put into practice.

Пример конкретного исполненияConcrete example

Нефтегазоконденсатное месторождение N находится на стадии разведки и освоения. Продуктивный разрез характеризуется сложным геологическим строением. Основные продуктивные рифейские отложения с точки зрения получения достоверных результатов при комплексной интерпретации данных ГИС также относятся к категории сложных, что обусловлено многокомпонентным составом скелета пород, сложной структурой емкостного пространства, пониженными значениями фильтрационно-емкостных свойств (ФЕС): предельно низкий диапазон изменения коэффициента общей пористости (0,1-3%) соизмерим с абсолютной погрешностью определения этого параметра радиоактивными и акустическими методами (2-2,5%).Oil and gas condensate field N is at the exploration and development stage. The productive section is characterized by a complex geological structure. The main productive Riphean deposits from the point of view of obtaining reliable results with a comprehensive interpretation of well logging data also belong to the category of complex ones, which is due to the multicomponent composition of the rock skeleton, the complex structure of the capacitive space, and the low values of the filtration-capacitive properties (FES): extremely low range of variation of the coefficient of total porosity (0.1-3%) is commensurate with the absolute error in determining this parameter by radioactive and acoustic methods (2-2.5%).

Литологический состав рифейских коллекторов неоднородный (доломит-кварц-глина). Структура пустот сложная (поры, каверны, трещины). Коллекторы гидрофобные, что обусловливает их высокие электрические сопротивления - до 5000 Ом·м по данным бокового каротажа. В этом разрезе практически не эффективны электрические методы ГИС - ПС, МК, БМК, БКЗ, ИК, ДК. В интервалах коллекторов не образуется глинистой корки в процессе бурения.The lithological composition of Riphean reservoirs is heterogeneous (dolomite-quartz-clay). The structure of voids is complex (pores, cavities, cracks). The reservoirs are hydrophobic, which leads to their high electrical resistance - up to 5000 Ohm · m according to lateral logging. In this context, the electric GIS methods are practically ineffective - PS, MK, BMK, BKZ, IK, DK. In the intervals of the reservoirs, a clay crust does not form during drilling.

Имеет место сложная фильтрационная система коллекторов. Макротрещины, обеспечивающие наибольшие дебиты скважин, не могут быть учтены при построении модели коллектора по керновым данным. В таких условиях известные способы ГИС и анализа керна практически не позволяют дифференцировать разрез скважины. Геофизическая характеристика пород нередко выглядит достаточно однородной и не отражает изменения реальных их свойств.There is a complex filtration system of collectors. Macrocracks that provide the highest flow rates of wells cannot be taken into account when constructing a reservoir model from core data. In such conditions, well-known methods of well logging and core analysis practically do not allow differentiating a well section. The geophysical characteristics of the rocks often look quite homogeneous and do not reflect changes in their real properties.

Некоторые характеристики одной из залежей месторождения приведены в табл.1.Some characteristics of one of the deposits of the field are given in Table 1.

По предлагаемому способу, используя отсчеты параметров по комплексу геофизических методов исследований, результаты некоторых их преобразований, а также данные ГТИ, полученные по скважине ТК-505, была создана база данных для выполнения многопараметрической классификации разреза скважины. Использованная база данных включала более 2000 строк по глубине по 29 переменным (табл.2).According to the proposed method, using parameter readings for a complex of geophysical research methods, the results of some of their transformations, as well as GTI data obtained from the TK-505 well, a database was created to perform multi-parameter classification of the well section. The database used included more than 2000 rows in depth over 29 variables (Table 2).

Основные характеристики распределения параметров приведены в табл.3.The main characteristics of the distribution of parameters are given in table.3.

Ряд параметров, распределение которых не подчинялось нормальному закону, были преобразованы, чтобы привести их к логнормальному виду, и стандартизованы вычитанием среднего и делением на среднеквадратичное отклонение (табл.4).A number of parameters, the distribution of which did not obey the normal law, were converted to bring them to a lognormal form, and were standardized by subtracting the mean and dividing by the standard deviation (Table 4).

К полученным стандартизованным значениям был применен один из статистических методов - метод иерархической классификации. Дерево классификации скважины показано на фиг.1.One of the statistical methods, the hierarchical classification method, was applied to the obtained standardized values. The well classification tree is shown in FIG.

Метод состоит в иерархическом попарном слиянии классов (кластеров), начиная с классов, состоящих из единственного элемента (в нашем случае - набора параметров, отнесенного к одной глубине). Классы для слияния выбираются так, чтобы дисперсия (разброс вокруг среднего, полученного слиянием класса) была минимальной. Метод предполагает одинаковое распределение всех параметров и чувствителен к наличию больших значений (асимметричное распределение), поэтому все параметры сначала преобразуются к симметричному распределению, наиболее близкому к нормальному, например, логарифмированием. Естественное число классов для данного разреза, соответствующее выделению примерно одинаковых по характеристикам и отличных друг от друга по совокупности параметров, выбирается на основании анализа дерева слияния, на котором вертикальная ось соответствует дисперсии: число классов, при котором начинается резкий рост дисперсии, означает слияние несходных по совокупности характеристик классов и хорошо видно на фиг.1: при разделении на 5 и меньшее число классов дисперсия увеличивается.The method consists in hierarchical pairwise merging of classes (clusters), starting with classes consisting of a single element (in our case, a set of parameters related to the same depth). Classes for merging are chosen so that the variance (scatter around the average obtained by merging the class) is minimal. The method assumes the same distribution of all parameters and is sensitive to the presence of large values (asymmetric distribution), so all parameters are first converted to a symmetric distribution that is closest to normal, for example, logarithm. The natural number of classes for this section, corresponding to the selection of parameters that are approximately the same in characteristics and differ from each other in the set of parameters, is selected based on the analysis of the merger tree, on which the vertical axis corresponds to the variance: the number of classes at which a sharp increase in dispersion begins, means a merger of dissimilar the aggregate characteristics of the classes and can be clearly seen in figure 1: when divided into 5 or fewer classes, the variance increases.

Были рассмотрены варианты классификации с выделением от 2 до 8 классов. Ниже применение способа проиллюстрируем на примере с выделением 8 классов. Проведено описание каждого класса по использованным для классификации параметрам. Пример такого описания приводится для варианта с выделением 8 классов и ограниченного числа (9 из 29) параметров (фиг.2).Classification options were considered with a separation from 2 to 8 classes. Below, we illustrate the application of the method using an example with the allocation of 8 classes. A description of each class is carried out according to the parameters used for classification. An example of such a description is given for the option with the allocation of 8 classes and a limited number (9 of 29) of the parameters (figure 2).

На чертеже нижняя и верхняя границы прямоугольника соответствуют 25% и 75% квантилям распределения (в прямоугольнике находится 50% выборки), точка в прямоугольнике показывает медиану, как наиболее надежную оценку истинного среднего, горизонтальные линии (к которым идут пунктиры) показывают доверительные границы для выборки, соответствующие среднеквадратическому отклонению; все, что выходит за доверительные границы, считается выбросами и показывается круглыми значками. Снизу приведены номера классов, для которых строилась оценка по соответствующим переменным. Анализ этой информации позволяет дать оценку степени отличия тех или иных параметров между классами, разброс этих параметров и т.д.In the drawing, the lower and upper boundaries of the rectangle correspond to 25% and 75% of the quantiles of the distribution (the rectangle contains 50% of the sample), the point in the rectangle shows the median as the most reliable estimate of the true average, horizontal lines (to which the dashed lines go) show confidence boundaries for the sample corresponding to standard deviation; anything that goes beyond confidence limits is considered outliers and is shown with round icons. Below are the numbers of classes for which an estimate was constructed for the corresponding variables. The analysis of this information allows us to assess the degree of difference of various parameters between classes, the scatter of these parameters, etc.

На фиг.3 приведены результаты расчленения разреза скважины ТК-505, в табл.5 - средние значения использованных при классификации параметров для каждого класса при различных вариантах их выделения (от 2 до 8).Figure 3 shows the results of the breakdown of the section of the well TK-505, in table 5 - the average values used in the classification of parameters for each class with different options for their allocation (from 2 to 8).

С учетом фундаментальных связей этих параметров с минеральным и элементным составом пород выполнена смысловая интерпретация выделенных классов (интервалов) по разрезу скважины ТК-505. Приведем в качестве примера вариант выделения 8 классов (фиг.4).Taking into account the fundamental relationships of these parameters with the mineral and elemental composition of the rocks, a semantic interpretation of the selected classes (intervals) along the TK-505 well section is performed. We give as an example the option of allocating 8 classes (figure 4).

Вендские отложения (до 2340 м), представленные катангской (интервал 2256-2316 м) и оскобинской свитами (интервал глубин 2316-2340 м), относятся к классам 1,5 и немного 2. По вещественному составу (соотношению основных породосоставляющих компонент: кальцит, доломит, неорганический остаток - глина, кварц) породы, относящиеся к классу 1, представляют глины с карбонатным цементом (7:17:51,8); класс 5 аналогичен 1 и представлен соотношением (1,6:1,2:20,4). Высокой карбонатностью обладают породы класса 2 (4,9:44,2:19,1). Высокая глинистость классов 1 и 5 подтверждается и максимальным (среди классов) содержанием калия, тория и урана (К=1,195÷1,895, Th=2,026÷3,365, U=2,39÷2,814).Vendian deposits (up to 2340 m), represented by the Katanga (interval 2256-2316 m) and Oskobin formations (depth interval 2316-2340 m), belong to classes 1.5 and a little 2. In terms of material composition (the ratio of the main rock-making components: calcite, dolomite, inorganic residue - clay, quartz) rocks belonging to class 1, are clay with carbonate cement (7: 17: 51.8); class 5 is similar to 1 and is represented by the ratio (1.6: 1.2: 20.4). The rocks of class 2 (4.9: 44.2: 19.1) have high carbonate content. The high clay content of classes 1 and 5 is also confirmed by the maximum (among the classes) content of potassium, thorium and uranium (K = 1.195 ÷ 1.895, Th = 2.026 ÷ 3.365, U = 2.39 ÷ 2.814).

Породы, относящиеся ко 2 классу, менее радиоактивны (К=0,1917, Th=0,5868, U=0,5843) (см. табл.5).The rocks belonging to class 2 are less radioactive (K = 0.1917, Th = 0.5868, U = 0.5843) (see Table 5).

Значения общего водородосодержания (W) в классах 1, 5 максимальны и коррелируются с концентрациями Тh и К, минимальные величины ИК (электрометрия) свидетельствуют о наличии глинистости в рассматриваемых классах.The values of total hydrogen content (W) in classes 1, 5 are maximum and correlate with the concentrations of Th and K, the minimum values of IR (electrometry) indicate the presence of clay in the classes under consideration.

Характерно, что углы падения пластов глинистых отложений (классы 1, 5 имеют значения 10,73÷12,9°, доломиты класса 2 - 30,76° (они имеют угловое несогласие, унаследованное от рифейских отложений).It is characteristic that the angles of incidence of clay deposits (classes 1, 5 have values of 10.73 ÷ 12.9 °, dolomites of class 2 - 30.76 ° (they have angular disagreement inherited from Riphean deposits).

По суммарному газопоказанию худшими являются породы, относящиеся к классу 5. Породы классов 1 и 2 имеют примерно равные показатели, на порядок превышающие класс 5. Последние могут интерпретироваться как газосодержащие отложения.According to the total gas indication, the rocks belonging to class 5 are the worst. The rocks of classes 1 and 2 have approximately equal indicators, an order of magnitude higher than class 5. The latter can be interpreted as gas-bearing deposits.

Рифейские отложения (ниже 2340 м) представлены классами 2, 3, 4, 6, 7, 8. Отложения класса 6 (интервал 2467-2500 м) по содержанию основных минеральных компонент (6,5:84,6:8,3) относятся к наиболее чистым доломитам. Породы класса 4 (интервалы 2405-2435 и 2445-2451) могут быть отнесены к известковистым окварцованным доломитам (33,2:60,3:6,4). Интервалы пород класса 8 (интервалы 2433-2446 и 2450-2457 м) имеют соотношение (9,1:59,2:31,7) - битумные (или кварцесодержащие) доломитовые разности. Отложения классов 3 (интервалы 2340-2360 м с включением пород класса 2, и 2385-2396 м) и 7 (интервал 2358-2405 с включением класса 3) имеют сопоставимые количества кальцита и доломита, но различаются по неорганическому остатку - возможно отношения представлены доломитами с высоким содержанием кварца. По описанию образцов каменного материала отложения классов 3, 7 - это доломиты, разбитые сетью трещин и каверн. Встречаются и брекчии, залеченные гидроокислами железа и вторичным доломитом, а также черным углеродистым веществом.Riphean deposits (below 2340 m) are represented by classes 2, 3, 4, 6, 7, 8. Deposits of class 6 (interval 2467-2500 m) in terms of the content of the main mineral components (6.5: 84.6: 8.3) are to the purest dolomites. Class 4 rocks (intervals 2405-2435 and 2445-2451) can be attributed to calcareous silicified dolomites (33.2: 60.3: 6.4). The intervals of the rocks of class 8 (intervals 2433-2446 and 2450-2457 m) have the ratio (9.1: 59.2: 31.7) - bitumen (or quartz-containing) dolomite differences. Deposits of classes 3 (intervals 2340-2360 m with the inclusion of rocks of class 2, and 2385-2396 m) and 7 (interval 2358-2405 with the inclusion of class 3) have comparable amounts of calcite and dolomite, but differ in inorganic residue - possibly the ratios are represented by dolomites high in quartz. According to the description of the samples of stone material, deposits of classes 3, 7 are dolomites, broken by a network of cracks and caverns. There are breccias healed with iron hydroxides and secondary dolomite, as well as black carbonaceous matter.

В рассматриваемом случае отложения классов 3 и 7 приурочены к кровельной части рифея, где они практически доминируют до глубины 2405 м. Отложения классов 4 и 8 преобладают в интервале глубин 2405-2458 м. Ниже до глубины 2468 м - класс 2. В интервале глубин 2467-2500 - отложения класса 6.In this case, deposits of classes 3 and 7 are confined to the roofing part of the Riphean, where they practically dominate to a depth of 2405 m. Deposits of classes 4 and 8 prevail in the depth range of 2405-2458 m. Below to a depth of 2468 m - class 2. In the depth range of 2467 -2500 - Class 6 sediments.

Классы 7 и 3 находятся над нефтенасыщенной частью и характеризуются как газонасыщенные. Наиболее низкими величинами УЭС обладают породы класса 6, из которых была получена минерализованная вода со значимым дебитом.Classes 7 and 3 are located above the oil-saturated part and are characterized as gas-saturated. Class 6 rocks have the lowest resistivity values, from which mineralized water with a significant flow rate was obtained.

Максимальные значения Est, характерные для уплотненных карбонатов (9,18÷9,22*105 неп.м.), приурочены к классам 4 и 8 и именно к ним относятся минимальные значения Δtp и Δts (138,4÷138,3 мкс/м и 258,3÷258,4 мкс/м), что также свидетельствует о весьма низких коллекторских свойствах данных классов; судя по величине W (2,11-2,14%) породы классов 4 и 8 соответствуют кавернозно-трещинному типу, т.к. резкое азимутальное несогласие (практически на 180°) нижележащих пород класса 8 с классом 4 (DA меняется с 63,4 до 215,1°) при близких углах падения слоев 25,86÷26,89° является предпосылкой для деструктизации монолита пород; по значениям суммарного газопоказания носителем УВ и коллекторских свойств наиболее вероятно являются отложения класса 4 нежели 8, в таком случае вода, занимающая трещинно-кавернозное пространство пород класса 6, благодаря угловому несогласию отложений классов 8-4 и присутствию плотных доломитов класса 2 блокируется на глубинах 2456-2468 м; однозначно являются нефтегазосодержащими пластами породы классов 3, 4, 7. Породы класса 6 водонасыщены.The maximum Est values characteristic of compacted carbonates (9.18 ÷ 9.22 * 10 5 n.m.) are confined to classes 4 and 8 and the minimum values Δtp and Δts (138.4 ÷ 138.3 μs) belong to them. / m and 258.3 ÷ 258.4 μs / m), which also indicates the very low reservoir properties of these classes; judging by the value of W (2.11-2.14%), rocks of classes 4 and 8 correspond to the cavernous-fractured type, because a sharp azimuthal disagreement (almost 180 °) of the underlying rocks of class 8 with class 4 (DA varies from 63.4 to 215.1 °) at close angles of incidence of 25.86 ÷ 26.89 ° is a prerequisite for the destruction of the rock monolith; according to the values of the total gas indication, the carrier of hydrocarbons and reservoir properties are most likely to be deposits of class 4 rather than 8, in which case the water occupying the fissure-cavernous space of rocks of class 6 is blocked at the depths of 2456 due to the angular disagreement of deposits of classes 8-4 and the presence of dense dolomites of class 2 -2468 m; rocks of classes 3, 4, 7 are clearly oil and gas containing strata. Class 6 rocks are water saturated.

Керном была охарактеризована относительно небольшая часть разреза скв. ТК-505 (см. фиг.4). Как видно, пористость по керну дифференцируется очень слабо. Лишь на границе вендских и рифейскх отложений выделяется прослой сильно неоднородных пород, отнесенных к классу 1, - глины с карбонатным цементом (кальцит: доломит: неорганический остаток - глина, кварц - в соотношении 7:17:51,8). Те или иные классы пластов описываются по содержанию минеральных компонент, используя данные анализа керна, где он отобран.Kern characterized a relatively small part of the SLE section. TK-505 (see figure 4). As can be seen, core porosity is very poorly differentiated. Only on the border of the Vendian and Riphean deposits is a layer of highly heterogeneous rocks classified as class 1 distinguished - clay with carbonate cement (calcite: dolomite: inorganic residue - clay, quartz - in the ratio of 7: 17: 51.8). These or those classes of formations are described by the content of mineral components, using core analysis data, where it is selected.

Проведены испытания и промыслово-гидродинамические исследования 5 интервалов, которые совпадают с различными выделенными классами пластов (см. фиг.4). По их результатам получают количественную характеристику этих классов пластов: класс 6 - коллекторы водонасыщенные; класс 1 - неколлекторы без притока; классы 3 и 7 - нефтегазонасыщенные с различным газосодержанием и классы 4 и 8 - в основном нефтенасыщенные. Определяют также фильтрационные характеристики (табл.6).Tests and field-hydrodynamic studies of 5 intervals were carried out, which coincide with different distinguished classes of formations (see figure 4). According to their results, a quantitative characteristic of these classes of formations is obtained: class 6 — water-saturated reservoirs; class 1 - non collectors without inflow; classes 3 and 7 are oil and gas saturated with different gas contents and classes 4 and 8 are mainly oil saturated. Filtration characteristics are also determined (Table 6).

Наглядно демонстрируется разрешающая способность предложенного способа классификации при последовательном увеличении количества выделяемых классов (см. фиг.3). Так, при выделении 2-х классов четко отделились отложения рифея и венда (раздел на глубине 2340 м); 3-х классов - рифей разделился на нефтегазоконденсатную часть (2340-2454) и водонасыщенную (2454-2500); 4-х классов - разделились газонасыщенные и нефтенасыщенные интервалы (2440-2407 и 2407-2454). Отметим, что между частями разреза с различной флюидонасыщенностью имеется некоторая переходная зона, где смежные классы перекрываются. Это еще одно подтверждение о надежности расчленения разреза скважины. С увеличением числа классов начинают проявляться особенности строения коллекторов и т.д. Ясно, что каждый вариант может быть проанализирован, как это было изложено выше при варианте выделения 8 классов.The resolution of the proposed classification method is clearly demonstrated with a sequential increase in the number of distinguished classes (see figure 3). So, when distinguishing 2 classes, the Riphean and Vendian deposits were clearly separated (section at a depth of 2340 m); 3 classes - the Riphean was divided into the oil and gas condensate part (2340-2454) and water-saturated (2454-2500); 4 classes - gas-saturated and oil-saturated intervals (2440-2407 and 2407-2454) were divided. Note that there is a transition zone between parts of the section with different fluid saturation, where adjacent classes overlap. This is another confirmation of the reliability of the section of the well. With an increase in the number of classes, structural features of collectors, etc., begin to appear. It is clear that each option can be analyzed, as described above with the option of distinguishing 8 classes.

На сводной модели разреза скважины ТК-505 по литологии и флюидонасыщению (см. фиг.4) колонки 1-4 представляют результаты стандартных исследований по ГИС и анализа керна. Видно, что расчленение разреза и по параметрам, и по флюидонасыщению весьма неопределенное. Дифференциация разреза по предлагаемому способу более четкая как по флюидонасыщенности, так и по характеристикам коллекторов (колонки 5-7) (подробное описание было выше). Такая информация, полученная по группе скважин, позволяет провести корреляцию между ними по классам (пластам) и построить геолого-гидродинамическую модель месторождения для выполнения многовариантных технологических расчетов показателей разработки на прогноз.In the composite model of the TK-505 well section for lithology and fluid saturation (see Fig. 4), columns 1-4 present the results of standard well logging and core analysis. It can be seen that the dismemberment of the section in terms of parameters and fluid saturation is very uncertain. The section differentiation according to the proposed method is clearer both in terms of fluid saturation and reservoir characteristics (columns 5-7) (a detailed description was above). Such information obtained from a group of wells allows us to correlate between them by classes (formations) and build a geological and hydrodynamic model of the field to perform multivariate technological calculations of development indicators for the forecast.

Важнейшее достоинство предлагаемого способа заключается еще в том, что подтвержденные испытаниями и гидродинамическими исследованиями результаты о характере флюидонасыщения и продуктивных характеристиках в небольших интервалах разреза могут быть распространены на большие толщины, отнесенные к тому или иному классу пластов, включая пропущенные при стандартной интерпретации продуктивные пласты и пропластки.The most important advantage of the proposed method lies in the fact that the results confirmed by tests and hydrodynamic studies on the nature of fluid saturation and productive characteristics in small intervals of a section can be extended to large thicknesses assigned to one or another class of formations, including productive formations and interlayers that were missing in the standard interpretation .

В приводимом примере исполнения, имея результаты классификации только по одной скважине, построенная геолого-гидродинамическая модель является условной.In the given execution example, having the classification results for only one well, the constructed geological and hydrodynamic model is conditional.

Другие исходные данные для моделирования, в том числе о свойствах флюидов, приняты поданным недропользователя.Other initial data for modeling, including the properties of the fluids, are accepted by the subsoil user.

Гидродинамическая модель построена с использованием симулятора Eclipse 2003А на основе геологической модели, созданной в пакете Petrel 2005. Сетка ячеек имеет размерность 17×10×30. Общее число ячеек 5100. Сетка моделирования описана с использованием геометрии угловых точек (corner point geometry). Размеры ячеек в горизонтальной плоскости составляют 100 м. Сетка разделена на 30 слоев, размер ячеек по вертикали составляет 5,4 м (фиг.5). При создании модели использован набор функций псевдокомпозиционного моделирования:The hydrodynamic model was built using the Eclipse 2003A simulator based on the geological model created in the Petrel 2005 package. The mesh of cells has a dimension of 17 × 10 × 30. The total number of cells is 5100. The simulation grid is described using corner point geometry. The mesh sizes in the horizontal plane are 100 m. The grid is divided into 30 layers, the vertical mesh size is 5.4 m (FIG. 5). When creating the model, a set of functions of pseudo-compositional modeling was used:

- PVTG (PVT свойства конденсата),- PVTG (PVT condensate properties),

- PVTO (PVT свойства растворенного в нефти газа),- PVTO (PVT properties of gas dissolved in oil),

- BOGI (зависимость объемного коэффициента насыщенной газом нефти от давления),- BOGI (dependence of the volumetric coefficient of gas saturated oil on pressure),

- RSGI (зависимость газового фактора от давления),- RSGI (dependence of the gas factor on pressure),

- BGGI (зависимость объемного коэффициента конденсата от давления),- BGGI (dependence of the volumetric coefficient of the condensate on pressure),

- RVGI (зависимость количества конденсата от давления).- RVGI (dependence of the amount of condensate on pressure).

Начальное распределение свойств модели рассчитывалось по методу последовательного моделирования по Гауссу (Sequential Gaussian Simulation).The initial distribution of model properties was calculated by the method of sequential Gaussian modeling (Sequential Gaussian Simulation).

ВНК принят на глубине - 2458 м. ГНК - на глубине - 2404 м. Начальное пластовое давление - 21,9 МПа.VNK was adopted at a depth of 2458 m. GOC at a depth of 2404 m. The initial reservoir pressure was 21.9 MPa.

Рассмотрены 3 варианта моделирования:3 modeling options are considered:

- вариант 1 - отбор нефти из нефтенасыщенной части коллектора с одновременной закачкой газа в газовую шапку выше отметки ГНК;- option 1 - the selection of oil from the oil-saturated part of the reservoir with the simultaneous injection of gas into the gas cap above the mark of the GOC;

- вариант 2 - отбор нефти из нефтенасыщенной части коллектора с одновременной закачкой воды в водоносную область ниже отметки ВНК;- option 2 - the selection of oil from the oil-saturated part of the reservoir with the simultaneous injection of water into the aquifer below the VNK mark;

- вариант 3 - отбор нефти из нефтенасыщенной и газа из газонасыщенной части с закачкой воды ниже ВНК и сухого газа выше ГНК (сайклинг-процесс).- option 3 - selection of oil from the oil-saturated and gas from the gas-saturated part with the injection of water below the BHC and dry gas above the BHC (cycling process).

Вовлекаемые при этих вариантах классы пластов и интервалы их залегания приводятся в табл.7. На фиг.6 приводится результирующий график, показывающий динамику коэффициента нефтеизвлечения при разных вариантах разработки: на конец расчетного периода значения коэффициента нефтеизвлечения по вариантам 1, 2 и 3 составляют, соответственно, 0,360; 0,237 и 0,414. Худший результат разработки по варианту 2 с закачкой воды; лучший - по варианту 3 путем выделения в самостоятельные объекты эксплуатации нефтенасыщенной и газо (конденсато)-насыщенной частей разреза. Объемы добычи конденсата не приведены. Они еще улучшают показатели 3 варианта. Принимается, что после выработки запасов нефти начнется отбор газа на истощение.The classes of formations involved in these options and the intervals of their occurrence are given in Table 7. 6 is a resulting graph showing the dynamics of the oil recovery coefficient for different development options: at the end of the billing period, the oil recovery coefficient values for options 1, 2 and 3 are, respectively, 0.360; 0.237 and 0.414. Worst development result for option 2 with water injection; the best one, according to option 3, by separating the oil-saturated and gas (condensate) -saturated sections of the section into separate operating facilities. Condensate production volumes are not shown. They still improve the performance of 3 options. It is assumed that after the depletion of oil reserves, gas will begin to exhaust gas.

Однако наиболее эффективным будет выполнение 3-го варианта путем проводки горизонтальных скважин (ГС) в интервалах пластов-классов (табл.8). Поскольку важнейшее значение при проводке ГС имеет расстояние ее ствола от водо-, газонефтяных контактов и коллекторских свойств пластов, по которым проходит горизонтальный ствол, то использование результатов многопараметрической классификации разреза играет огромную роль для оптимизации профиля скважины. Сопоставление результатов моделирования процесса разработки вертикальными или горизонтальными скважинами в условиях анализируемого месторождения N показало, что коэффициент нефтеизвлечения может быть до 20 пунктов выше, чем вертикальными. При этих расчетах было принято, что длина ГС составляет 700 м, ствол скважины ориентирован перпендикулярно к основному направлению трещиноватости коллекторов. Коэффициент нефтеизвлечения, рассчитанный в проектном документе на разработку месторождения, составляет 38,9%, что существенно ниже достигаемой по предложенному способу величины (41,4+20=61,4%).However, the most effective will be the implementation of the 3rd option by drilling horizontal wells (horizontal wells) in the intervals of the reservoir classes (Table 8). Since the location of the wellbore from the water and gas-oil contacts and reservoir properties of the formations along which the horizontal wellbore is of prime importance when conducting a wellbore, the use of the results of multi-parameter classification of the section plays a huge role in optimizing the well profile. A comparison of the simulation results of the development process with vertical or horizontal wells in the conditions of the analyzed field N showed that the oil recovery coefficient can be up to 20 points higher than vertical ones. In these calculations, it was assumed that the length of the well is 700 m, the wellbore is oriented perpendicular to the main direction of fracture of the reservoirs. The oil recovery coefficient calculated in the project document for the development of the field is 38.9%, which is significantly lower than the value achieved by the proposed method (41.4 + 20 = 61.4%).

Таблица 1Table 1 ПоказателиIndicators ЗначениеValue Средняя глубина залегания, мAverage depth, m 20082008 Средняя эффективная газонефтенасыщенная толщина, мThe average effective gas-saturated thickness, m 97,897.8 Пористость, %Porosity,% 0,97-1,370.97-1.37 Средняя нефтенасыщенность пласта, доли ед.The average oil saturation of the reservoir, the share of units 0,512-0,8900.512-0.890 Проницаемость, мдарсиPermeability, mdarsi 1,7-8,41.7-8.4 Коэффициент песчанистости, доли ед.Sandiness coefficient, fractions of units 0,570.57 Коэффициент расчлененности, доли ед.The coefficient of dissection, the share of units 42,042.0 Начальная пластовая температура, °СInitial reservoir temperature, ° C 30,030,0 Начальное пластовое давление, МПаInitial reservoir pressure, MPa 21,121.1 Вязкость нефти в пластовых условиях, мПа·сOil viscosity in reservoir conditions, MPa · s 1,341.34 Плотность нефти в пластовых условиях, кг/м3 The density of oil in reservoir conditions, kg / m 3 0,7130.713 Давление насыщения нефти газом, МПаSaturation pressure of oil with gas, MPa 21,121.1 Газосодержание нефти, м3Gas content of oil, m 3 / t 174174

Таблица 2table 2 № п/пNo. p / p НазваниеTitle Ед. изм.Units rev. ОписаниеDescription 1one DTPDtp мкс/сμs / s Интервальное время пробега продольной волныLongitudinal Wavelength 22 ALPHADIP1ALPHADIP1 1/м1m Коэффициент затухания продольной волны дипольный 1Dipole longitudinal wave attenuation coefficient 1 33 ALPHADIP2ALPHADIP2 1/м1m Коэффициент затухания продольной волны дипольный 2Dipole longitudinal wave attenuation coefficient 2 4four DTSDTS мкс/мμs / m Интервальное время пробега поперечной волныShear Wave Interval 55 ALPHAPALPHAP 1/м1m Коэффициент затухания продольной волныLongitudinal wave attenuation coefficient 66 ALPHASALPHAS 1/м1m Коэффициент затухания поперечной волныShear wave attenuation coefficient 77 DTSTDTST мкс/мμs / m Интервальное время пробега волны Лэмба-СтоунлиLamb-Stoneley Wave Interval Run Time 88 ESTEst усл. ед.conv. units Энергия волны Лэмба-СтоунлиLamb-Stoneley Wave Energy 99 ALPHASTALPHAST 1/м1m Коэффициент затухания волны Лэмба-СтоунлиLamb-Stoneley wave attenuation coefficient 1010 POTAPota %% Концентрация калияPotassium concentration 11eleven THORThor PPMPPM Концентрация торияThorium concentration 1212 URANURAN PPMPPM Концентрация уранаUranium concentration 1313 DTSfDtsf мкс/мμs / m Интервальное время быстрой поперечной волныShear Wave Interval Time 14fourteen DTLMDTLM мкс/мμs / m Интервальное время пробега модельной волныInterval travel time of the model wave 15fifteen WW %% Общее водородосодержание, %Total hydrogen content,% 1616 DTDIP1DTDIP1 мкс/мμs / m Интервальное время пробега поперечной волны DIP1Shear Wave Interval DIP1 1717 DTDIP2DTDIP2 мкс/мμs / m Интервальное время пробега поперечной волны DIP2Shear Wave Interval DIP2 18eighteen DIPDip GRADGRAD Угол падения, градусыAngle of incidence, degrees 1919 DADA GRADGRAD Азимут падения, градусыAzimuth of a fall, degrees 20twenty ADIPADIP GRADGRAD Кажущийся угол падения, градусыApparent angle of incidence, degrees 2121 AZAPAZAP GRADGRAD Кажущийся азимут падения, градусыApparent azimuth of fall, degrees 2222 MethanMethan %, отн% rel МетанMethane 2323 EthanEthan %, отн% rel ЭтанEthane 2424 PropanPropan %, отн% rel ПропанPropane 2525 ButanButan %, отн% rel БутанButane 2626 PenthanPenthan %, отн% rel ПентанPentane 2727 GexanGexan %, отн% rel ГексанHexane 2828 GasSumGasum %% Суммарный газTotal gas 2929th IKIK Ом·мOhm m Индукционный каротажInduction Logging

Figure 00000001
Figure 00000001

Таблица 4Table 4 ПараметрParameter СреднееAverage СКОSKO ПараметрParameter СреднееAverage СКОSKO DTPDtp 0.8810.881 0.0060.006 DTDIP1DTDIP1 1.9691.969 0.0570.057 ALPHADIP1ALPHADIP1 2.4892.489 1.641.64 DTDIP2DTDIP2 2.1092.109 0.0550.055 ALPHADIP2ALPHADIP2 2.6372.637 1.8841.884 DIPDip 3.1373.137 0.7090.709 DTSDTS 1.8021.802 0.0650.065 DADA 159.868159.868 95.74795.747 ALPHAPALPHAP 1.9951.995 0.9090.909 ADIPADIP 3.1553.155 0.6850.685 ALPHASALPHAS 2.3122.312 1.9331.933 AZAPAZAP 181.839181.839 105.873105.873 DTSTDTST 704.607704.607 11.12911.129 MethanMethan -4.552-4.552 1.5291.529 ESTEst 450.326450.326 71.80871.808 EthanEthan -6.653-6.653 1.1121.112 ALPHASTALPHAST 3.3623.362 0.580.58 PropanPropan -8.7-8.7 1.5631.563 POTAPota -1.27-1.27 0.8150.815 ButanButan -7.771-7.771 0.9430.943 THORThor -0.263-0.263 0.8810.881 PenthanPenthan -5.863-5.863 0.6650.665 URANURAN -0.401-0.401 0.8270.827 GexanGexan -4.236-4.236 0.5150.515 DTSfDtsf 1.9241.924 0.070.07 GasSumGasum -3.32-3.32 0.7420.742 DTLMDTLM 706.711706.711 9.3029.302 IKIK 82.88582.885 32.39832.398 WW 1.0911.091 0.7270.727

Figure 00000002
Figure 00000003
Figure 00000002
Figure 00000003

Figure 00000004
Figure 00000004

Таблица 7Table 7 ВариантOption ТехнологияTechnology КлассыClasses Интервалы разрезаCut Intervals 1one отбор нефтиoil extraction 4,84.8 2404-2434; 2446-2451; 2434-2446; 2451-24582404-2434; 2446-2451; 2434-2446; 2451-2458 закачка газаgas injection 33 2341-23592341-2359 22 отбор нефтиoil extraction 4,84.8 2404-2434; 2446-2451; 2434-2446; 2451-24582404-2434; 2446-2451; 2434-2446; 2451-2458 закачка водыwater injection 66 2468-25002468-2500 33 отбор газаgas extraction 3,73,7 2386-2396; 2359-23862386-2396; 2359-2386 отбор нефтиoil extraction 4,84.8 2404-2434; 2434-24462404-2434; 2434-2446 закачка газаgas injection 33 2341-23592341-2359 закачка водыwater injection 66 2468-25002468-2500

Таблица 8Table 8 Назначение скважиныWell designation Класс пластаFormation class Интервал перфорации, мPerforation interval, m ГазодобывающаяGas producing 77 2359-23862359-2386 ГазонагнетательнаяGas injection 33 2341-23592341-2359 НефтедобывающаяOil producing 4four 2404-24342404-2434 ВодонагнетательнаяWater injection 66 2468-25002468-2500

Claims (2)

1. Способ разработки неоднородного массивного или многопластового газонефтяного или нефтегазоконденсатного месторождения, включающий бурение скважин, отбор керна, проведение комплекса геофизических исследований, испытания и гидродинамические исследования отдельных интервалов по толщине разреза, отличающийся тем, что создают базу данных по результатам измерений различных, принятых для каждого региона геофизических параметров для каждой 0,05-0,15 м толщины пласта, параметры с асимметричным логнормальным распределением преобразуют к нормальному распределению, стандартизуют все параметры вычитанием среднего и делением на среднеквадратичное отклонение, после чего выполняют иерархическую процедуру построения классов, состоящую в последовательном объединении прослоев таким образом, что полученные классы имеют минимальный разброс параметров вокруг центра - дисперсию, выделяя по резкому росту дисперсии естественное для данного разреза число классов с отличающимися для них геофизическими параметрами, вычисляют диапазоны изменения и средние значения этих параметров для каждого выделенного класса пластов, с учетом фундаментальных связей этих параметров с минеральным и элементным составами пород определяют характер литологии, фильтрационно-емкостные свойства и флюидонасыщение каждого класса на качественном уровне, дополнительно используют результаты анализа керна, испытаний и гидродинамических исследований в интервалах отдельных классов-пластов и указанные характеристики определяют на количественном уровне, проводят корреляцию выделенных классов между скважинами, строят геолого-гидродинамическую модель продуктивной залежи, на которой оценивают эффективность системы разработки при различных вариантах включения выделенных классов в эксплуатационные объекты расчетами технологических показателей разработки с определением уровней добычи углеводородов на различных режимах с оценкой прогнозных коэффициентов нефтеизвлечения отдельных эксплуатационных объектов и месторождения в целом и лучший вариант реализуют на практике.1. A method of developing a heterogeneous massive or multi-layer gas-oil or oil-gas condensate field, including drilling wells, coring, conducting a complex of geophysical studies, tests and hydrodynamic studies of individual intervals along the thickness of the section, characterized in that they create a database of different measurement results accepted for each region of geophysical parameters for each 0.05-0.15 m of formation thickness, parameters with asymmetric lognormal distribution transform to distribution, standardize all parameters by subtracting the mean and dividing by the standard deviation, and then perform a hierarchical class construction procedure consisting in sequentially combining interlayers in such a way that the resulting classes have a minimum dispersion of parameters around the center - variance, highlighting a sharp increase in variance that is natural for a given section number of classes with different geophysical parameters, calculate the ranges of variation and average values of these parameters for I of each selected class of formations, taking into account the fundamental relationships of these parameters with the mineral and elemental composition of rocks, determine the nature of lithology, filtration-capacity properties and fluid saturation of each class at a qualitative level, additionally use the results of core analysis, tests and hydrodynamic studies in the intervals of individual classes of formations and the indicated characteristics are determined at a quantitative level, the selected classes are correlated between wells, geological and hydrodynamic are built a model of a productive reservoir, on which the effectiveness of the development system is assessed for various options for incorporating the allocated classes into production facilities by calculating technological development indicators with determining hydrocarbon production levels in various modes with estimating the predicted oil recovery coefficients of individual production facilities and the field as a whole and the best option is implemented in practice . 2. Способ по п.1, отличающийся тем, что при разработке месторождения горизонтальными скважинами проводят оптимизацию профиля их стволов с использованием данных о продуктивности выделенных классов-пластов и расстояния до газонефтяного или водонефтяного контактов. 2. The method according to claim 1, characterized in that when developing a field by horizontal wells, they optimize the profile of their shafts using data on the productivity of the selected reservoir classes and the distance to the gas-oil or water-oil contacts.
RU2009143585/03A 2009-11-25 2009-11-25 Procedure for development of non-uniform massive or multipay gas-oil or oil-gas-condensate field RU2432450C2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2009143585/03A RU2432450C2 (en) 2009-11-25 2009-11-25 Procedure for development of non-uniform massive or multipay gas-oil or oil-gas-condensate field

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2009143585/03A RU2432450C2 (en) 2009-11-25 2009-11-25 Procedure for development of non-uniform massive or multipay gas-oil or oil-gas-condensate field

Publications (2)

Publication Number Publication Date
RU2009143585A RU2009143585A (en) 2011-05-27
RU2432450C2 true RU2432450C2 (en) 2011-10-27

Family

ID=44734653

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2009143585/03A RU2432450C2 (en) 2009-11-25 2009-11-25 Procedure for development of non-uniform massive or multipay gas-oil or oil-gas-condensate field

Country Status (1)

Country Link
RU (1) RU2432450C2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2493362C1 (en) * 2012-09-19 2013-09-20 Открытое акционерное общество "Татнефть" им. В.Д. Шашина Method of oil filed development
RU2513787C1 (en) * 2012-10-17 2014-04-20 Открытое акционерное общество "Всероссийский нефтегазовый научно-исследовательский институт имени академика А.П. Крылова" (ОАО "ВНИИнефть") Method for oil deposit development based on system address action
RU2567581C1 (en) * 2015-02-05 2015-11-10 Открытое акционерное общество "Татнефть" им. В.Д. Шашина Method for determination of formation intervals with viscous or highly viscous oil
RU2602249C1 (en) * 2015-10-20 2016-11-10 Шлюмберже Текнолоджи Б.В. Method for determining characteristics of gas-oil transition zone in cased well
RU2695418C1 (en) * 2019-01-22 2019-07-23 Общество с ограниченной ответственностью "ЛУКОЙЛ - Западная Сибирь" Method for development of non-uniform oil deposit

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115653580A (en) * 2022-11-07 2023-01-31 成都理工大学 Method for depicting internal structure of complex lithologic reservoir

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2493362C1 (en) * 2012-09-19 2013-09-20 Открытое акционерное общество "Татнефть" им. В.Д. Шашина Method of oil filed development
RU2513787C1 (en) * 2012-10-17 2014-04-20 Открытое акционерное общество "Всероссийский нефтегазовый научно-исследовательский институт имени академика А.П. Крылова" (ОАО "ВНИИнефть") Method for oil deposit development based on system address action
RU2567581C1 (en) * 2015-02-05 2015-11-10 Открытое акционерное общество "Татнефть" им. В.Д. Шашина Method for determination of formation intervals with viscous or highly viscous oil
RU2602249C1 (en) * 2015-10-20 2016-11-10 Шлюмберже Текнолоджи Б.В. Method for determining characteristics of gas-oil transition zone in cased well
US10309219B2 (en) 2015-10-20 2019-06-04 Schlumberger Technology Corporation Method for determining characteristics of a gas-oil transition zone in an uncased well
RU2695418C1 (en) * 2019-01-22 2019-07-23 Общество с ограниченной ответственностью "ЛУКОЙЛ - Западная Сибирь" Method for development of non-uniform oil deposit

Also Published As

Publication number Publication date
RU2009143585A (en) 2011-05-27

Similar Documents

Publication Publication Date Title
Sullivan et al. Architectural analysis of deep-water outcrops: Implications for exploration and development of the Diana sub-basin, western Gulf of Mexico
RU2432450C2 (en) Procedure for development of non-uniform massive or multipay gas-oil or oil-gas-condensate field
KR101853334B1 (en) Method for predicting exploitation site of unconventional resources applied to sequence stratigraphic analysis
Wang et al. Integrated reservoir characterization study of a carbonate ramp reservoir: Seminole San Andres Unit, Gaines County, Texas
AlRassas et al. CO2 storage capacity estimation under geological uncertainty using 3-D geological modeling of unconventional reservoir rocks in Shahejie Formation, block Nv32, China
Euzen et al. Well log cluster analysis: an innovative tool for unconventional exploration
Bhattacharya et al. Flow unit modeling and fine-scale predicted permeability validation in Atokan sandstones: Norcan East field, Kansas
Mahmud et al. Petrophysical interpretations of subsurface stratigraphic correlations, Baram Delta, Sarawak, Malaysia
Sadeq et al. Porosity and permeability analysis from well logs and core in fracture, vugy and intercrystalline carbonate reservoirs
Hughes et al. The unconventional Carboniferous reservoirs of the Greater Kirby Misperton gas field and their potential: North Yorkshire's sleeping giant
Shankar et al. Waterflood Performance Analyses for the Bhagyam Viscous Oil Reservoir
Ghadami et al. Enhanced History Matching and Prediction Using Integrated Analytical and Numerical Modeling Approach
Eikeland et al. Dry gas reinjection in a strong waterdrive gas/condensate field increases condensate recovery—case study: The Sleipner Øst Ty field, South Viking Graben, Norwegian North Sea
Elshahawi et al. Combining continuous fluid typing, wireline formation testers, and geochemical measurements for an improved understanding of reservoir architecture
Ozkaya Modeling finite-fracture networks in a partially fractured reservoir in the Middle East
Archer PD 6 (l) Reservoir Definition and Characterisation for Analysis and Simulation
Punanova Oil and gas possibility of crystalline basement taking into account development in it of non-structural traps of combined type
Hovorka et al. Comparing carbon sequestration in an oil reservoir to sequestration in a brine formation—field study
Ohno et al. Analysis of an interwell tracer test in a depleted heavy-oil reservoir
Gaviria Garcia Reservoir simulation of co2 sequestration and enhanced oil recovery in Tensleep Formation, Teapot Dome field
Walker et al. Three case studies of progress in quantitative seismic-engineering integration
Gutierrez Oseguera et al. Evaluation of Flow Units and Capillary Pressures of the Giant Chicontepec Tight Oil Paleochannel in Mexico and a Fresh Look at Drilling and Completions
Walters Organic geochemistry at varying scales: from kilometres to ångstroms
Kim et al. Post-Waterflooding Oil Reservoir Surveillance Using Pulsed Neutron Well Logging in Upper Burgan Formation, Northern Kuwait
Martyns-Yellowe Integrated Lithostratigraphic, Formation and Structural Evaluation of the Migrant Structure, Sable Subbasin, Offshore Nova Scotia.

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20131126