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 PDFInfo
- 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
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
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. Построение сетки.
Claims (1)
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)
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)
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 |
-
2013
- 2013-10-15 RU RU2013146133/08A patent/RU2542874C1/en not_active IP Right Cessation
Patent Citations (6)
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)
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 |