RU2542874C1 - Method for adaptive breakdown of three-dimensional grid into layers when creating three-dimensional digital geologic models of formations - Google Patents

Method for adaptive breakdown of three-dimensional grid into layers when creating three-dimensional digital geologic models of formations Download PDF

Info

Publication number
RU2542874C1
RU2542874C1 RU2013146133/08A RU2013146133A RU2542874C1 RU 2542874 C1 RU2542874 C1 RU 2542874C1 RU 2013146133/08 A RU2013146133/08 A RU 2013146133/08A RU 2013146133 A RU2013146133 A RU 2013146133A RU 2542874 C1 RU2542874 C1 RU 2542874C1
Authority
RU
Russia
Prior art keywords
grid
layers
dimensional
layering
layer
Prior art date
Application number
RU2013146133/08A
Other languages
Russian (ru)
Inventor
Антон Александрович Колесниченко
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 RU2013146133/08A priority Critical patent/RU2542874C1/en
Application granted granted Critical
Publication of RU2542874C1 publication Critical patent/RU2542874C1/en

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

FIELD: physics.
SUBSTANCE: in the method, creation or correction of grid layers is carried out by successively searching source data which characterise layering; new grid layers are created at layering boundaries based on the source data under the condition of absence, in the allowable region vertically, of grid layers created in advance, otherwise the closest grid layer on the detected layering boundary is corrected vertically, the grid layer is corrected horizontally within a triangulation grid constructed on points of intersection of the source data with the corrected grid layer in the vicinity of said point.
EFFECT: providing accurate display of non-uniformity of the geologic structure of a formation in a three-dimensional grid with the least number of layers.
2 dwg

Description

Изобретение относится к горному делу, в частности к способам разбивки трехмерной сетки на слои, и может быть использовано в геологии при создании трехмерных цифровых геологических моделей пластов.The invention relates to mining, in particular to methods for dividing a three-dimensional grid into layers, and can be used in geology when creating three-dimensional digital geological models of formations.

Трехмерная сетка - это ячеистый каркас, в объеме которого выполняются все основные этапы геологического моделирования. Трехмерная сетка строится внутри так называемой «зоны» или пласта. Пласт представляет собой объем между двумя поверхностями, расположенными сверху и снизу. Для точного воспроизведения геологического строения моделируемых пластов весь объем делится на мелкие ячейки. Совокупность этих ячеек и есть трехмерная геологическая сетка. Процесс создания трехмерной геологической сетки заключается в разбиении объема пласта на мелкие ячейки в соответствии с заданными правилами, определяемыми непосредственно геологом, проводящим моделирование. Правила разбиения объема на ячейки необходимо задавать как по латерали («нарезка» сетки на «столбцы» и «ряды»), так и по вертикали («нарезка» сетки на слои) [1]. Известны два способа разбивки трехмерной сетки на слои:A three-dimensional grid is a cellular framework in the volume of which all the main stages of geological modeling are carried out. A three-dimensional grid is built inside the so-called “zone” or formation. The layer is the volume between two surfaces located above and below. To accurately reproduce the geological structure of the simulated formations, the entire volume is divided into small cells. The combination of these cells is a three-dimensional geological grid. The process of creating a three-dimensional geological grid consists in dividing the reservoir volume into small cells in accordance with predetermined rules, determined directly by the geologist conducting the simulation. The rules for dividing the volume into cells must be set both laterally (“slicing” the grid into “columns” and “rows”), and vertically (“slicing” the grid into layers) [1]. There are two known methods for breaking a three-dimensional grid into layers:

1. Пропорциональная разбивка. При данном типе разбивки вся толща пласта делится на заданное равное количество слоев, вне зависимости от его общей толщины. Таким образом, при пропорциональной разбивке любая область пласта содержит заданное количество слоев, однако толщина каждого слоя в разных частях пласта может меняться.1. Proportional breakdown. With this type of breakdown, the entire thickness of the reservoir is divided into a given equal number of layers, regardless of its total thickness. Thus, with a proportional breakdown, any area of the formation contains a predetermined number of layers, however, the thickness of each layer in different parts of the formation can vary.

2. Параллельная разбивка. При параллельной разбивке вся толща пласта разделяется на слои заданной толщины, а границы слоев проводятся параллельно либо кровле, либо подошве этого пласта, либо какой-нибудь другой поверхности.2. Parallel breakdown. With a parallel breakdown, the entire thickness of the formation is divided into layers of a given thickness, and the boundaries of the layers are parallel to either the roof, or the sole of this formation, or some other surface.

Количество слоев, на которые делится весь объем пласта, зависит от необходимой степени детальности создаваемой модели, т.е. задается минимально необходимая толщина слоев, а следовательно, чем детальнее необходима модель, тем большее количество слоев необходимо задавать. Затем на ячейки сетки вдоль траекторий скважин, пробуренных через пласт, выполняется перенос (осреднение) результатов интерпретации ГИС (геофизических исследований скважин) - кривых фаций, литологии, пористости, нефтенасыщенности и др. Все эти данные напрямую характеризуют слоистое строение пласта, его неоднородность. При их переносе на сеточную область, границы слоев которой не соответствуют границам слоистости исходных скважинных данных, происходит искажение исходных данных и исчезновение особо тонких прослоев. Таким образом, недостатком всех способов является то, что в процессе формирования слоев не учитывается исходная информация, по которой будет проводиться дальнейшее создание геологической модели. Вторым существенным недостатком является то, что в целях достижения высокой точности моделей необходимо увеличивать количество слоев в сетке модели, что, в свою очередь, приводит к увеличению размерности модели и усложняет возможность ее дальнейшего использования.The number of layers into which the entire volume of the reservoir is divided depends on the necessary degree of detail of the created model, i.e. the minimum required layer thickness is set, and therefore, the more detailed the model is needed, the more layers must be specified. Then, meshing along the trajectories of wells drilled through the reservoir, the results are transferred (averaged) of the interpretation of well logs (geophysical surveys of the wells) - facies curves, lithology, porosity, oil saturation, etc. All these data directly characterize the layered structure of the formation, its heterogeneity. When they are transferred to the grid region, the layer boundaries of which do not correspond to the layered boundaries of the original well data, the initial data are distorted and especially thin layers disappear. Thus, the disadvantage of all the methods is that in the process of forming the layers, the initial information on which the further creation of the geological model will be carried out is not taken into account. The second significant drawback is that in order to achieve high accuracy of the models, it is necessary to increase the number of layers in the model grid, which, in turn, leads to an increase in the dimension of the model and complicates the possibility of its further use.

Наиболее близким способом, взятым за прототип, является способ пропорциональной разбивки модели на слои [1, прил. 3, стр.247]. Недостатками данного способа являются ошибки при переносе исходных данных на сетку модели и увеличение размерности сетки модели при увеличении количества слоев в целях минимизации ошибок переноса исходных данных.The closest way, taken as a prototype, is the method of proportional dividing the model into layers [1, adj. 3, p. 247]. The disadvantages of this method are errors when transferring the source data to the model grid and increasing the dimension of the model grid with an increase in the number of layers in order to minimize the errors of transferring the source data.

Задачей изобретения является создание способа адекватного отображения неоднородности геологического строения пласта в трехмерной сетке модели при наименьшем количестве слоев.The objective of the invention is to provide a method for adequately displaying heterogeneity of the geological structure of the formation in a three-dimensional grid of the model with the least number of layers.

Поставленная задача решается тем, что в способе адаптивной разбивки трехмерной сетки на слои используется исходная информация по скважинам, характеризующая слоистую неоднородность пласта, и производится корректировка слоев сетки находящихся в определенной близости от границ слоистости по исходным данным.The problem is solved in that in the method of adaptive dividing of a three-dimensional grid into layers, the initial information on the wells, which characterizes the layered heterogeneity of the formation, is used, and the grid layers located in a certain vicinity of the layering boundaries are corrected according to the initial data.

При формировании слоев сетки производится перебор всех исходных данных и проводится поиск границ слоистости. В случае обнаружения такой границы проводится проверка на наличие уже существующего слоя сетки в пределах допустимой области (ДО). ДО - это максимальное расстояние, на которое возможна корректировка слоя сетки по вертикали, определяется пользователем. При отсутствии существующего слоя в ДО проводится создание нового слоя сетки. Вновь созданный слой распространяется на всю область модели по пропорциональной или параллельной схеме. При наличии в ДО одного или нескольких ранее созданных слоев сетки производится корректировка, перенос по вертикали, ближайшего слоя сетки на соответствующую границу слоистости. Для сохранения плавности поверхности слоя с помощью аппроксимации проводится корректировка слоя по горизонтали в пределах триангуляционной сетки, построенной по точкам пересечения исходных данных с корректируемым слоем сетки и определенными величинами поправок в этих точках. Таким образом, в результате перебора всех исходных данных формируется нарезка трехмерной сетки на слои, полностью характеризующая слоистость в исходной информации при использовании минимального числа слоев сетки.When forming the grid layers, all source data are enumerated and a search is made for the boundaries of layering. If such a boundary is detected, a check is carried out for the presence of an existing mesh layer within the allowable area (DO). DO - this is the maximum distance at which it is possible to adjust the grid layer vertically, is determined by the user. In the absence of an existing layer in DO, a new mesh layer is created. The newly created layer extends to the entire area of the model in a proportional or parallel pattern. If there are one or several previously created mesh layers in the DO, the vertical adjustment of the nearest mesh layer to the corresponding layering boundary is carried out. To preserve the smoothness of the layer surface by means of approximation, the layer is horizontally adjusted within the triangulation grid, constructed from the points of intersection of the initial data with the corrected grid layer and certain correction values at these points. Thus, as a result of enumerating all the initial data, a three-dimensional grid is cut into layers, which completely characterizes the layering in the initial information using the minimum number of grid layers.

Сущность осуществления предложенного способа поясняется чертежом, на котором:The essence of the proposed method is illustrated in the drawing, on which:

Фиг.1. Вид на трехмерную сетку (13) сбоку. По скважинам (1, 2, 3) слоистость характеризуется наличием (4) и отсутствием (5) коллектора. Слои сетки (6) создаются на границах слоев, определенных в скважине (1), прослеживаются по всей площади и корректируются (8) в скважине (2), так как граница слоистости находится в пределах допустимой области (7) возможной корректировки. По скважине (2) новые слои не проводятся, так как имеющиеся границы слоистости уже охарактеризованы имеющимися скорректированными слоями сетки. В скважине (3) имеется граница слоистости (9), не имеющая в допустимой области слоев сетки, поэтому проводится создание нового слоя сетки (10).Figure 1. Side view of a three-dimensional grid (13). In wells (1, 2, 3), layering is characterized by the presence of (4) and the absence of (5) reservoir. Layers of the grid (6) are created at the boundaries of the layers defined in the well (1), are traced over the entire area and are adjusted (8) in the well (2), since the layering boundary is within the allowable region (7) of possible adjustment. No new layers are drawn along well (2), since the existing layering boundaries are already characterized by the existing adjusted mesh layers. In the well (3), there is a layering boundary (9), which does not have a grid layer in the allowable region; therefore, a new grid layer is created (10).

Фиг.2. Вид на трехмерную сетку (13) сверху. Область влияния (11) определенной подвижки (14) слоя в скважине (2) по площади трехмерной сетки. Триангуляционная сеть по исходным данным (12).Figure 2. View of the three-dimensional grid (13) from above. The influence area (11) of a certain displacement (14) of the layer in the well (2) over the area of the three-dimensional grid. Triangulation network according to the source data (12).

Источники информацииInformation sources

1. Закревский К.Е. «Геологическое 3D моделирование». М.: ООО «ИПЦ „Маска"», 2009.1. Zakrevsky K.E. "Geological 3D modeling." M .: LLC "CPI" Mask "", 2009.

Приложение 3. Построение сетки.Appendix 3. Mesh construction.

Claims (1)

Способ адаптивной разбивки трехмерной сетки на слои при создании трехмерных цифровых геологических моделей пластов состоит в том, что создание или корректировка слоев сетки проводится при последовательном переборе исходных данных (ИД), характеризующих слоистость, при этом новые слои сетки создаются на границах слоистости по ИД при условии отсутствия в допустимой области по вертикали ранее созданных слоев сетки, иначе проводится корректировка по вертикали ближайшего слоя сетки на обнаруженную границу слоистости, а также корректируется слой сетки по горизонтали в пределах триангуляционной сетки, построенной по точкам пересечения исходных данных с корректируемым слоем сетки, в окрестностях данной точки. A method for adaptively dividing a three-dimensional grid into layers when creating three-dimensional digital geological models of formations is that the creation or adjustment of the grid layers is carried out by sequentially sorting the initial data (ID) characterizing the layering, while new grid layers are created at the boundaries of the layered by ID, provided the absence of previously created mesh layers in the permissible vertical region; otherwise, the vertical adjustment of the nearest mesh layer is carried out to the detected layering boundary, and is also adjusted I am the horizontal grid layer within the triangulation grid, constructed at the points of intersection of the source data with the adjusted grid layer, in the vicinity of this point.
RU2013146133/08A 2013-10-15 2013-10-15 Method for adaptive breakdown of three-dimensional grid into layers when creating three-dimensional digital geologic models of formations RU2542874C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2013146133/08A RU2542874C1 (en) 2013-10-15 2013-10-15 Method for adaptive breakdown of three-dimensional grid into layers when creating three-dimensional digital geologic models of formations

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2013146133/08A RU2542874C1 (en) 2013-10-15 2013-10-15 Method for adaptive breakdown of three-dimensional grid into layers when creating three-dimensional digital geologic models of formations

Publications (1)

Publication Number Publication Date
RU2542874C1 true RU2542874C1 (en) 2015-02-27

Family

ID=53289995

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2013146133/08A RU2542874C1 (en) 2013-10-15 2013-10-15 Method for adaptive breakdown of three-dimensional grid into layers when creating three-dimensional digital geologic models of formations

Country Status (1)

Country Link
RU (1) RU2542874C1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113219529A (en) * 2021-04-30 2021-08-06 电子科技大学 Boundary-constrained three-dimensional structure geological model construction method
RU2787499C1 (en) * 2022-10-20 2023-01-09 Федеральное государственное бюджетное образовательное учреждение высшего образования "Удмуртский государственный университет" Method for mapping a non-anticlinal oil trap

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7248259B2 (en) * 2001-12-12 2007-07-24 Technoguide As Three dimensional geological model construction
EA200970543A1 (en) * 2006-12-05 2009-10-30 Конокофиллипс Компани METHOD AND DEVICE FOR STRENGTHENING LAYERS IN GEOLOGICAL MODELING
RU2382193C2 (en) * 2004-11-04 2010-02-20 Бейкер Хьюз Инкорпорейтед Multidimensional inversion of logging data using different scales and method of constructing images of deep-lying formations
WO2011059535A1 (en) * 2009-11-12 2011-05-19 Exxonmobil Upstream Research Company Method and apparatus for reservoir modeling and simulation
US20120026167A1 (en) * 2010-07-16 2012-02-02 Ran Longmin Method for generating a hex-dominant mesh of a geometrically complex basin
RU2481599C2 (en) * 2006-06-21 2013-05-10 ТЕРРАСПАРК ДЖИОСАЙЕНСИЗ, ЭлЭлСи Sedimentation system detection

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7248259B2 (en) * 2001-12-12 2007-07-24 Technoguide As Three dimensional geological model construction
RU2382193C2 (en) * 2004-11-04 2010-02-20 Бейкер Хьюз Инкорпорейтед Multidimensional inversion of logging data using different scales and method of constructing images of deep-lying formations
RU2481599C2 (en) * 2006-06-21 2013-05-10 ТЕРРАСПАРК ДЖИОСАЙЕНСИЗ, ЭлЭлСи Sedimentation system detection
EA200970543A1 (en) * 2006-12-05 2009-10-30 Конокофиллипс Компани METHOD AND DEVICE FOR STRENGTHENING LAYERS IN GEOLOGICAL MODELING
WO2011059535A1 (en) * 2009-11-12 2011-05-19 Exxonmobil Upstream Research Company Method and apparatus for reservoir modeling and simulation
US20120026167A1 (en) * 2010-07-16 2012-02-02 Ran Longmin Method for generating a hex-dominant mesh of a geometrically complex basin

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113219529A (en) * 2021-04-30 2021-08-06 电子科技大学 Boundary-constrained three-dimensional structure geological model construction method
RU2787499C1 (en) * 2022-10-20 2023-01-09 Федеральное государственное бюджетное образовательное учреждение высшего образования "Удмуртский государственный университет" Method for mapping a non-anticlinal oil trap

Similar Documents

Publication Publication Date Title
CN104879103B (en) A kind of seperated layer water injection effect analysis method
CN104809277B (en) A kind of Oil in Super-low Permeability compact reservoir horizontal well Geological Modeling
AU2013397958B2 (en) A simulation-to-seismic workflow construed from core based rock typing and enhanced by rock replacement modeling
CN105651676B (en) A kind of reservoir heterogeneity characterizing method under horizontal well rule well pattern
Lin et al. A California statewide three-dimensional seismic velocity model from both absolute and differential times
CN104992468A (en) Fracture-cavern type carbonate hydrocarbon reservoir three-dimensional geological modeling method
CN104297787B (en) The three-dimensional petrofacies data processing method and processing device of fluvial facies Low permeability and competent sand reservoir
US10309216B2 (en) Method of upscaling a discrete fracture network model
CA2997608C (en) History matching of hydrocarbon production from heterogenous reservoirs
CN102313790A (en) Submarine geologic body carbon dioxide sequestration potential assessment method
CN103114840A (en) Calculating method and device of organic carbon content of high-over high mature shale
CN106154322B (en) Log curve correction method and apparatus
CN110056346B (en) Oil reservoir three-dimensional original water saturation simulation method based on trend change function
Hryb et al. Unlocking the true potential of the Vaca Muerta Shale via an integrated completion optimization approach
CN104063626A (en) Predication method of outer-source entrapment fullness degree
CN104914470A (en) Carbonate rock fracture-cavity reservoir reserve correction method
CN106354893A (en) Method for geologic modeling with horizontal well data during drilling
Du et al. Generalization of dual-porosity-system representation and reservoir simulation of hydraulic fracturing-stimulated shale gas reservoirs
CN107765307A (en) Reef beach type Reservoir Body modeling method and system
RU2542874C1 (en) Method for adaptive breakdown of three-dimensional grid into layers when creating three-dimensional digital geologic models of formations
AU2013398344B2 (en) Local updating of 3D geocellular model
CN107368688B (en) Method and device for determining marine single sand body
CN105089659A (en) Conglomerate oil reservoir seepage unit recognizing method
CN115705519A (en) Horizontal well track optimization design method based on three-dimensional geological model
Xinghe et al. Constraining method of stochastic modeling for fluvial petroleum reservoir controlled by depositional facies using wells and seismic data

Legal Events

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

Effective date: 20171016