RU2275660C1 - Method of selecting and forecasting areas with different types of geologic profile - Google Patents
Method of selecting and forecasting areas with different types of geologic profile Download PDFInfo
- Publication number
- RU2275660C1 RU2275660C1 RU2004124365/28A RU2004124365A RU2275660C1 RU 2275660 C1 RU2275660 C1 RU 2275660C1 RU 2004124365/28 A RU2004124365/28 A RU 2004124365/28A RU 2004124365 A RU2004124365 A RU 2004124365A RU 2275660 C1 RU2275660 C1 RU 2275660C1
- Authority
- RU
- Russia
- Prior art keywords
- data
- points
- spectral
- clusters
- different types
- Prior art date
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Изобретение относится к нефтяной геологии и может быть использовано для оптимизации размещения новых скважин на исследуемом объекте по комплексу данных сейсморазведки и исследования данных по скважинам.The invention relates to petroleum geology and can be used to optimize the placement of new wells at the studied object according to a set of seismic data and research data on wells.
Используя данные наземной сейсморазведки и результаты исследования скважин (данные по каротажу и по анализу керна), рассчитывают по сейсморазведочным данным спектрально-временные параметры (СВП), несущие в себе косвенную информацию о типах разреза в заданном временном (глубинном) интервале, а по скважинным данным составляют суждение о емкостно-фильтрационных свойствах соответствующих геологических пластов. Далее данные одновременно по всем СВП подвергают многомерному анализу на предмет выявления обособленных скоплений точек (кластеров) в пространстве СВП, отвечающих различным типам разреза. Прогноз выполняется путем сопоставления значений СВП в точках территории с значениями их для выделенных кластеров. Технический результат - повышение надежности обоснования заложения разведочных и эксплуатационных скважин и, тем самым, повышение эффективности геологоразведочных работ.Using the data of ground-based seismic exploration and the results of well research (data on logging and core analysis), spectral-temporal parameters (SVP) are calculated using seismic data, which carry indirect information about the types of sections in a given time (depth) interval, and according to borehole data make a judgment on the capacitance-filtration properties of the corresponding geological formations. Further, the data for all SVPs are simultaneously subjected to multivariate analysis to identify isolated clusters of points (clusters) in the SVP space corresponding to different types of sections. The forecast is carried out by comparing the values of the SVP at the points of the territory with their values for the selected clusters. The technical result is an increase in the reliability of the justification for the laying of exploration and production wells and, thereby, an increase in the efficiency of exploration work.
Геологической основой заявленного предложения является то обстоятельство, что литологические и вслед за ними физические (упругие) свойства пород сменяются в той или иной степени дискретно при переходе от одного литотипа к другому. Вследствие этого точки в пространстве СВП, соответствующие разным типам разреза, образуют обособленные скопления с центрами рассеяния, соответствующими типичным для этих типов сочетаниям признаков.The geological basis of the proposed proposal is the fact that the lithological and after them physical (elastic) properties of the rocks change to one degree or another discretely during the transition from one lithotype to another. As a result of this, points in the SVP space corresponding to different types of sections form separate clusters with scattering centers corresponding to combinations of features typical of these types of sections.
Известны способы решения задачи классификации типов геологической среды, основанные на распознавании образов (см. «Интерпретация данных сейсморазведки с использованием методов распознания образов» Галагон Е.А., Кузнецов О.А., Литвинов А.Я., Тальвирский Д.Б., Обзор ВИЭМС, Разведочная геофизика, 1984 г., 49 с.). В основе этих способов лежит анализ расстояний р между точками, соответствующими объектам, в m-мерном признаком пространствеKnown methods for solving the problem of classifying types of geological environment based on pattern recognition (see "Interpretation of seismic data using pattern recognition methods" Galagon EA, Kuznetsov OA, Litvinov A.Ya., Talvirsky DB, VIEMS Review, Exploration Geophysics, 1984, 49 pp.). The basis of these methods is the analysis of the distances p between points corresponding to objects in the m-dimensional sign of space
При этом предполагается, что все признаки, характеризующие объекты, равно информативны. В то же время, очевидно, что вариации слабо информативных признаков, обусловленные посторонними причинами, вносят в расчеты расстояний помеху, затрудняющую решение задачи.It is assumed that all the signs characterizing the objects are equally informative. At the same time, it is obvious that variations of poorly informative features, caused by extraneous reasons, introduce interference into the distance calculations, which complicates the solution of the problem.
На практике, однако, проверка справедливости этих гипотез и расчет значений параметров распределения, входящих в решающее правило, весьма затруднительны из-за малой выборки эталонных объектов отдельных типов. Попытки преодолеть эти трудности путем анализа облака точек в проекциях на координатные плоскости наталкиваются на громоздкость перебора проекций (большое число сочетаний исследуемых признаков по два) и на трудности разделения гетерогенной совокупности. Дело в том, что при взаимной корреляции физических параметров, часто имеющей место на практике, облака кластеров ориентированы под углом к осям этих параметров, и на проекциях распределения кластеров перекрывают друг друга. Сказанное является причиной невысокой эффективности и ненадежности способов, основанных на методах распознавания образов.In practice, however, verifying the validity of these hypotheses and calculating the distribution parameters included in the decision rule are very difficult due to the small sample of reference objects of individual types. Attempts to overcome these difficulties by analyzing a point cloud in projections onto coordinate planes encounter the cumbersome enumeration of projections (a large number of combinations of the studied features in two) and the difficulties of separating a heterogeneous population. The fact is that with the mutual correlation of physical parameters, which often takes place in practice, cluster clouds are oriented at an angle to the axes of these parameters, and overlap each other on the projections of the distribution of clusters. The foregoing is the reason for the low efficiency and unreliability of methods based on image recognition methods.
Наиболее близким к заявленному является способ, основанный на определении и сопоставлении спектрально-временных образов по сейсмическим профилям и в окрестности скважин (Методика картирования типов геологического разреза в межскважинном пространстве по данным сейсморазведки, Славкин B.C., Копелевич Е.А., Давыдова Е.А., Мушин И.А., Геофизика, том 4, М., ЕАГО, 1999, с.21-24). Недостатком этого способа является то, что решающее правило не формализовано, а основывается на экспертных оценках. Известны также попытки развития того же способа, предусматривающие построение карт полей СВП по всей исследуемой территории и анализ характера этих полей в окрестности эталонных скважин (Давыдова Е.А. Методика картирования нефтепродуктивных типов геологического разреза по спектрально-временным параметрам, с.14-16). Такой подход также несвободен от недостатков. Анализировать одновременно несколько (4-6) карт затруднительно, а если предварительно выбрать из них по экспертным оценкам одну-две, то в отобранных картах неизбежно присутствует помеха, а в отброшенных - неиспользуемая полезная информация.Closest to the claimed one is a method based on the determination and comparison of spectral-temporal images from seismic profiles and in the vicinity of wells (Methodology for mapping the types of geological sections in interwell space according to seismic data, Slavkin BC, Kopelevich EA, Davydova EA , Mushin I.A., Geophysics, Volume 4, M., EAGO, 1999, pp.21-24). The disadvantage of this method is that the decision rule is not formalized, but based on expert estimates. Attempts are also known to develop the same method, involving the construction of maps of SVP fields throughout the study area and an analysis of the nature of these fields in the vicinity of reference wells (Davydova EA, Methodology for mapping oil-producing types of a geological section by spectral-temporal parameters, p.14-16) . This approach is also not free from disadvantages. It is difficult to analyze several (4-6) cards at the same time, and if you first select one or two of them according to expert estimates, then in the selected cards there will inevitably be interference, and in discarded ones, unused useful information.
Технической задачей, на решение которой направлено данное изобретение, является повышение надежности и точности выявления и картирования областей с разными типами разреза, а значит более обоснованное заложение новых разведочных и эксплуатационных скважин.The technical problem to which this invention is directed is to increase the reliability and accuracy of identifying and mapping areas with different types of sections, and therefore more reasonable laying of new exploration and production wells.
Задача, которая решается предлагаемым способом, выражается в прогнозе участков с различными типами разреза одновременно по нескольким (m=4-6) спектрально-временным признакам с использованием количественных критериев.The problem, which is solved by the proposed method, is expressed in the prediction of areas with different types of sections simultaneously for several (m = 4-6) spectral-temporal signs using quantitative criteria.
Способ выделения и прогноза участков с различными типами геологического разреза включает проведение наземной сейсморазведки и исследований скважин (данные по каротажу и по анализу керна). По сейсморазведочным данным рассчитывают спектрально-временные параметры (СВП). Эти параметры характеризуют спектральную плотность энергетического спектра: - частотного (K1-К3) и временного (К4-K6) - в трех модификациях по каждой оси. Конкретно, K1 - это отношение энергии высоких и низких частот; К2 и К3 - удельная спектральная плотность энергетического частотного спектра СВАН-колонки (по оси частот f), умноженная на средневзвешенную (для К2) или на максимальную (К3) частоту спектра. Параметры К4, К5 и К6 по оси времен t симметричны параметрам K1, К2 и К3; формулы их отличаются заменой аргумента f на t. Как показали исследования, указанные СВП несут в себе косвенную информацию о типах разреза в заданном временном (глубинном) интервале. Одновременно составляют суждение о типах разреза соответствующих геологических пластов по скважинным данным.A method for identifying and predicting sections with different types of geological sections includes conducting surface seismic and well surveys (data on logging and core analysis). The seismic data calculate the spectral-temporal parameters (SVP). These parameters characterize the spectral density of the energy spectrum: - frequency (K 1 -K 3 ) and temporary (K 4 -K 6 ) - in three modifications on each axis. Specifically, K 1 is the ratio of the energy of high and low frequencies; K 2 and K 3 - specific spectral density of the energy frequency spectrum of the SWAN column (along the frequency axis f), multiplied by the weighted average (for K 2 ) or the maximum (K 3 ) frequency of the spectrum. The parameters K 4 , K 5 and K 6 along the time axis t are symmetric to the parameters K 1 , K 2 and K 3 ; their formulas differ by replacing the argument f with t. As studies have shown, these SVPs carry indirect information about the types of sections in a given time (deep) interval. At the same time, they make judgments about the types of sections of the corresponding geological formations according to well data.
Далее выполняется пересчет анализируемых СВП (с необходимым сглаживанием) на регулярную сетку в пределах исследуемой территории. Собственно анализу подлежат m-мерные векторы значений СВП (m=4-6) в точках фактических наблюдений и в дополнительных точках, квазиравномерно покрывающих территорию; значения векторов получаются путем интерполяции из узловых значений сетки. С использованием компьютерных программ анализируют структуру облака точек, соответствующих этим векторам, в проекциях на координатные плоскости признакового пространства. Подмножества точек из участков с разными типами разреза образуют в этом облаке обособленные скопления (кластеры).Then, the analyzed SVPs (with the necessary smoothing) are recalculated to a regular grid within the study area. The actual analysis is subject to m-dimensional vectors of SVP values (m = 4-6) at the points of actual observations and at additional points that quasi-uniformly cover the territory; vector values are obtained by interpolation from the nodal values of the grid. Using computer programs, the structure of the cloud of points corresponding to these vectors is analyzed in projections onto the coordinate planes of the feature space. Subsets of points from areas with different types of sections form isolated clusters (clusters) in this cloud.
Важной особенностью разработанного способа является использование формального аппарата факторного (компонентного) анализа - переход от пространства исходных, признаков (спектрально-временных параметров) к пространству соответствующих им факторов. Факторизация разворачивает оси в соответствии с упорядоченностью размещения точек (ориентировкой облаков кластеров) в m-мерном признаковом пространстве, что повышает контрастность разделения кластеров в проекциях на координатные плоскости. При этом факторы ранжируются по убыванию оценки изменчивости (разброса точек) вдоль их осей, поэтому в рассмотрении оставляются только mf<m главных факторов, суммарная доля изменчивости по которым не менее 90% от общей дисперсии. Остальные факторы исключаются из анализа как слабо информативные. Одновременно тем самым в несколько раз сокращается размерность задачи таксономии - без сокращения размерности анализируемых данных, поскольку каждый из оставленных для анализа факторов представляет собой линейную комбинацию всех m-признаков и косинусов углов поворота.An important feature of the developed method is the use of the formal apparatus of factor (component) analysis - the transition from the space of source features (spectral-temporal parameters) to the space of the factors corresponding to them. Factorization rotates the axes in accordance with the ordering of the distribution of points (orientation of the clouds of clusters) in the m-dimensional feature space, which increases the contrast of the separation of clusters in projections onto coordinate planes. In this case, the factors are ranked in descending order of variability (spread of points) along their axes, therefore, only mf <m main factors are left in consideration, the total share of variability along which is at least 90% of the total variance. Other factors are excluded from the analysis as poorly informative. At the same time, the dimension of the taxonomy problem is reduced several times - without reducing the dimension of the data being analyzed, since each of the factors left for analysis is a linear combination of all m-features and cosines of rotation angles.
Манипулируя различными выборками объектов анализа для факторизации, задают различные параметры поворота факторных осей и подбирают вариант, дающий наиболее контрастное разделение кластеров с точками, соответствующими скважинам с разными типами разреза, в проекции на факторные координатные плоскости.Manipulating different samples of analysis objects for factorization, they set various parameters for the rotation of factor axes and select the option that gives the most contrast separation of clusters with points corresponding to wells with different types of sections in a projection onto factor coordinate planes.
Собственно прогноз участков с разными типами разреза выполняется путем сравнения векторов СВП для тестируемых точек на территории (точек геофизических профилей и др.) с кластерами выделенных типов, преобразовав предварительно компоненты этих векторов в факторы. В итоге тестируемые точки получают пометки принадлежности к тому или иному типу.Actually, the prediction of sections with different types of sections is performed by comparing the SVP vectors for the tested points on the territory (points of geophysical profiles, etc.) with clusters of selected types, previously converting the components of these vectors into factors. As a result, the tested points receive marks of belonging to one or another type.
Результатом, достигнутым при использовании данного предложения, является повышение надежности картирования зон разных типов разреза. Это обеспечивает резкое снижение затрат на бурение новых разведочных и эксплуатационных скважин и, таким образом, повышает эффективность геологоразведочных работ.The result achieved by using this proposal is to increase the reliability of mapping zones of different types of sections. This provides a sharp reduction in the cost of drilling new exploration and production wells and, thus, increases the efficiency of exploration.
Claims (1)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2004124365/28A RU2275660C1 (en) | 2004-08-12 | 2004-08-12 | Method of selecting and forecasting areas with different types of geologic profile |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2004124365/28A RU2275660C1 (en) | 2004-08-12 | 2004-08-12 | Method of selecting and forecasting areas with different types of geologic profile |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2275660C1 true RU2275660C1 (en) | 2006-04-27 |
Family
ID=36655640
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2004124365/28A RU2275660C1 (en) | 2004-08-12 | 2004-08-12 | Method of selecting and forecasting areas with different types of geologic profile |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2275660C1 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014066504A1 (en) * | 2012-10-23 | 2014-05-01 | Schlumberger Canada Limited | Identifying formation, matrix and fluid related characteristics from subsurface data using factor analysis |
CN104991274A (en) * | 2015-07-03 | 2015-10-21 | 中国石油大学(华东) | Single-trap level favorable region optimal selection method under multi-geological factor quantitative constraints |
CN105221144A (en) * | 2015-09-21 | 2016-01-06 | 中国石油天然气股份有限公司 | Determine method and the device of reservoir reserve |
CN112180464A (en) * | 2020-11-03 | 2021-01-05 | 中国石油化工股份有限公司 | Reservoir physical property identification method |
CN112444868A (en) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | Seismic facies analysis method based on improved K-means algorithm |
-
2004
- 2004-08-12 RU RU2004124365/28A patent/RU2275660C1/en not_active IP Right Cessation
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014066504A1 (en) * | 2012-10-23 | 2014-05-01 | Schlumberger Canada Limited | Identifying formation, matrix and fluid related characteristics from subsurface data using factor analysis |
US9223055B2 (en) | 2012-10-23 | 2015-12-29 | Schlumberger Technology Corporation | Identifying formation, matrix and fluid related characteristics from subsurface data using factor analysis |
CN104991274A (en) * | 2015-07-03 | 2015-10-21 | 中国石油大学(华东) | Single-trap level favorable region optimal selection method under multi-geological factor quantitative constraints |
CN104991274B (en) * | 2015-07-03 | 2017-10-03 | 中国石油大学(华东) | The Favorable Areas method for optimizing of single trap level under a variety of geologic(al) factor quantity constraints |
CN105221144A (en) * | 2015-09-21 | 2016-01-06 | 中国石油天然气股份有限公司 | Determine method and the device of reservoir reserve |
CN105221144B (en) * | 2015-09-21 | 2018-04-06 | 中国石油天然气股份有限公司 | Determine the method and device of reservoir reserve |
CN112444868A (en) * | 2019-08-30 | 2021-03-05 | 中国石油化工股份有限公司 | Seismic facies analysis method based on improved K-means algorithm |
CN112444868B (en) * | 2019-08-30 | 2024-04-09 | 中国石油化工股份有限公司 | Seismic phase analysis method based on improved K-means algorithm |
CN112180464A (en) * | 2020-11-03 | 2021-01-05 | 中国石油化工股份有限公司 | Reservoir physical property identification method |
CN112180464B (en) * | 2020-11-03 | 2024-05-24 | 中国石油化工股份有限公司 | Reservoir physical property identification method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11693140B2 (en) | Identifying hydrocarbon reserves of a subterranean region using a reservoir earth model that models characteristics of the region | |
RU2474846C2 (en) | Method and apparatus for multidimensional data analysis to identify rock heterogeneity | |
West et al. | Interactive seismic facies classification using textural attributes and neural networks | |
AU2013337322B2 (en) | Seismic waveform classification system and method | |
US7519476B1 (en) | Method of seismic interpretation | |
CN111596978A (en) | Web page display method, module and system for lithofacies classification by artificial intelligence | |
US20210317726A1 (en) | Allocating resources for implementing a well-planning process | |
Venkataramanan et al. | An unsupervised learning algorithm to compute fluid volumes from NMR T1-T2 logs in unconventional reservoirs | |
US11815650B2 (en) | Optimization of well-planning process for identifying hydrocarbon reserves using an integrated multi-dimensional geological model | |
EP3047096A1 (en) | Identifying geological formation depth structure using well log data | |
CN108717211B (en) | A kind of prediction technique of the Effective source rocks abundance in few well area | |
CN111506861B (en) | Method for calculating crack strength of favorable region of target layer | |
Ding et al. | A novel neural network for seismic anisotropy and fracture porosity measurements in carbonate reservoirs | |
WO2018013004A1 (en) | Method for predicting effective collector capacity | |
US20220120933A1 (en) | Method of detection of hydrocarbon horizontal slippage passages | |
CN112946751B (en) | Method and system for predicting plane distribution of different reservoirs | |
RU2275660C1 (en) | Method of selecting and forecasting areas with different types of geologic profile | |
CN116168224A (en) | Machine learning lithology automatic identification method based on imaging gravel content | |
CN112444859A (en) | Shale reservoir fracture identification method and system for cooperative metamorphic ant body | |
US11434757B2 (en) | Direct hydrocarbon indicators analysis informed by machine learning processes | |
RU2718135C1 (en) | Method and system for prediction of effective thicknesses in inter-well space during construction of geological model based on spectral curves clustering method | |
RU22830U1 (en) | AUTOMATED SYSTEM FOR SEARCHING, EXPLORING AND EVALUATING THE OPERATIONAL PROPERTIES OF DEPOSITS AND DEPOSITS OF MINERAL RESOURCES AND FORECAST OF TECTONIC AND LITHOPHYSICAL PROPERTIES OF GEOLOGICAL MEDIA | |
Velis | Statistical segmentation of geophysical log data | |
Zhao et al. | Application of gradient structure tensor method in CBM fracture identification and sweet spot prediction | |
Mahmoud | Automatic characterization and quantitative analysis of seismic facies in naturally fractured reservoir: Case study of Amguid Messaoud field, Algeria. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20180813 |