CN105184867B - The three dimensional contour line method of rule-based grid and Corner-point Grids technology - Google Patents
The three dimensional contour line method of rule-based grid and Corner-point Grids technology Download PDFInfo
- Publication number
- CN105184867B CN105184867B CN201510631247.8A CN201510631247A CN105184867B CN 105184867 B CN105184867 B CN 105184867B CN 201510631247 A CN201510631247 A CN 201510631247A CN 105184867 B CN105184867 B CN 105184867B
- Authority
- CN
- China
- Prior art keywords
- grid
- thickness
- corner
- stratum
- top bottom
- Prior art date
- Legal status (The legal status 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 status listed.)
- Expired - Fee Related
Links
- 238000005516 engineering process Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 title claims description 52
- 238000003325 tomography Methods 0.000 claims abstract description 46
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 22
- 238000011160 research Methods 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims description 13
- 230000035699 permeability Effects 0.000 claims description 5
- 230000000717 retained effect Effects 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 9
- 239000010410 layer Substances 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 6
- 239000003245 coal Substances 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000003491 array Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000012800 visualization Methods 0.000 description 3
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical group [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 2
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- RXZBMPWDPOLZGW-HEWSMUCTSA-N (Z)-roxithromycin Chemical compound O([C@@H]1[C@@H](C)C(=O)O[C@@H]([C@@]([C@H](O)[C@@H](C)C(=N\OCOCCOC)/[C@H](C)C[C@@](C)(O)[C@H](O[C@H]2[C@@H]([C@H](C[C@@H](C)O2)N(C)C)O)[C@H]1C)(C)O)CC)[C@H]1C[C@@](C)(OC)[C@@H](O)[C@H](C)O1 RXZBMPWDPOLZGW-HEWSMUCTSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 241001415846 Procellariidae Species 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011439 discrete element method Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
Landscapes
- Processing Or Creating Images (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The elevation data of stratum top bottom interface of the present invention first in borehole data, the regular grid of a certain survey region stratum top bottom surface is generated by interpolation algorithm, simultaneously, enter row interpolation using the formation thickness data in borehole data, the regional stratum thickness regular grid is generated, keeps the consistent of grid quantity；Then according to the boundary coordinate and stratal pinch-outs line number evidence in research area, border file is established, rejects the grid beyond border in grid file；Secondly, tomography binding is carried out to stratum top bottom surface regular grid according to layer data；It is finally based on the three-dimensional modeling that two kinds of regular grids realize stratum using Corner-point Grids technology.
Description
Technical field
Three dimensional contour line method of the present invention, more particularly to a kind of rule-based grid and Corner-point Grids technology
Three dimensional contour line method.
Background technology
For a long time, scientific research personnel took much count of Geologic modeling and visualization, early in international exploration geophysics in 1992
Scholar association and European Association of Exploration Geophysicists have just set up the SEG/EARG three-dimensional modelings committee, come within 1993 and add
The scholar Simon W.Houlding to put on airs propose the concept of three-dimensional geological modeling, and illustrate and realize geology three-dimensional visible
Some basic skills of change technology, including triangle network generating method, triangulation network surface model construction method, THREE DIMENSIONAL TRIANGULATION NET solidification side
Method, the delimitation of geological objects boundary and connection etc., substantially reflect the core achievement of 3 d visualization of geologic volume technology at that time.Pin
To the particularity and complexity of Geologic modeling, French professor Mallet proposes discrete smooth interpolation technology, it has also become well-known
Three-dimensional geological modeling business software GOCAD core technology, separately there is more famous commercial applications software, such as ROXAR companies
RMS, the Petrel softwares etc. of Schlumberger, substantially represent the highest of current geological prospecting three-dimensional visualization application
It is horizontal.
Domestic aspect, Beijing Titan Technology Co., Ltd. is in internationally famous Canadian Apollo science and technology group GIS
On basis of software, the thinking exploitation based on framework modeling has formed TITAN T3D, utilizes parallel or substantially parallel cross-sectional data
Set up the true three-dimensional entity model of the arbitrarily complicated shaped objects of three dimensions；The reason of Beijing Li Zheng Information Technology Co., Ltds is just
GIS Development platform software can utilize isopleth or the landform altitude point of distribution at random, establish dimensionally prime number
Word model.
Share more than 20 kinds of spatial modeling method in current achievement in research to be suggested, from the angle of three-dimensional geological model substantially
Four classes, i.e. surface model, body Model, mixed model and object model can be divided into.
1) the structure mould method based on surface model：The expression to three-dimensional geologic surface is laid particular emphasis on, such as topographical surface, geological stratification
Face etc., said three-dimensional body is expressed by the surface of body.The conventional structure mould method based on surface model has TIN structure mould
(TIN), grid structure mould (Grid), wire frame structure mould (Wire Frame), in addition, also including boundary representation structure mould (B-rep), section
Structure mould (Section), multi-DEMs structure mould etc., its shortcoming of the model for making to establish in this way are the absence of belonging to geologic body inside
Property expression, it is difficult to geologic body carry out spatial analysis.It is only used for Two-dimensional numerical simulation.
2) the structure mould based on body Model：The border of three-dimensional geologic and internal overall expression are laid particular emphasis on, using volume elements as base
This unit expresses 3D solid.Because the attribute of volume elements can be with independent description and storage, thus entity can be carried out three-dimensional
Spatial operation and analysis, currently used structure mould method have tetrahedron mesh (TEN), entity geometry (CSG), Octree structure
Mould (Octree), triangular prism structure mould (TP) and generalized tri-prism (GTP) etc., shortcoming are that the model storage that this method is established is empty
Between big, model unit irregularities, be not easy to algorithm traversal, calculating speed is slow.
3) structure mould is mixed：Because geological phenomenon is extremely complex, simple structure mould based on surface model and based on body Model
Structure mould all can not accurately take into account the border of expression three-dimensional geologic and the change of built-in attribute, and a variety of models is integrated and mixed
Conjunction can then utilize each single model possessed advantage when expressing different spaces entity, can realize to three-dimensional geologic more
Completely, effectively express, TIN-Grid data describe to can be regarded as a kind of mixed model in GIS.In addition, CSG-
Octree and TEN-Octree structures mould are also all typical mixed models, are memory space and amount of calculation the shortcomings that this model
It is bigger, belong to fine modeling method.
Corner-point Grids are a kind of structured grids, propose in the seventies of last century six, are ground in recent years in reservoir numerical simulation
Study carefully middle widely used, it can easily realize the edge fitting in region, be concentrated etc. suitable for fluid and surface force
Calculate, major advantage is that principle is simple, mess generation speed is fast, quality is good, data structure is simple.
The content of the invention
The technical problem to be solved in the present invention is that existing numerical simulation for oil-gas reservoir method operates when reservoir grid is established
Method is mostly complex, it is necessary to which simulating personnel has very sturdy professional knowledge, and software application is operated very
Skillfully, otherwise easily malfunctioning, cause final numerical simulation to fail, this undoubtedly improves the threshold that uses of reservoir simulation software, and
Add the difficulty of numerical simulation.
In order to solve this technical problem, the invention provides a kind of rule-based grid and the three-dimensional of Corner-point Grids technology
Stratum Modeling, including：
S1：For the research area of three dimension strata model to be established, extraction stratigraphic boundary, pinching line coordinates, tomography coordinate and
Borehole data；
S2：Stratum top bottom surface elevation data and formation thickness data are obtained according to the borehole data；Respectively enter step
S30 and step S40；
S30：Row interpolation is entered to the stratum top bottom surface elevation data, obtains stratum top bottom surface regular grid；Then utilize
Stratum top bottom surface regular grid is modified described in stratigraphic boundary and pinch-out boundary coordinate pair；
S40：Row interpolation is entered to the formation thickness data, obtains thickness regular grid；
In the step S30 and S40, in addition to：
Tomography binding is carried out at least one the stratum top bottom surface regular grid and thickness regular grid, while causes institute
The quantity of stratum top bottom surface regular grid and thickness regular grid, coordinate, arrangement mode is stated to be consistent；
S5：Depth assignment is carried out to the plane Corner-point Grids using the thickness regular grid, so as to be modeled, obtained
To three-dimensional Corner-point Grids；
S6：Threedimensional model is obtained using the three-dimensional Corner-point Grids.
Optionally, when being modified in the step S30 to the stratum top bottom surface regular grid, the stratum is removed
In top bottom surface regular grid outside pinch-out boundary coordinate connecting line, and the grid disjoint with the pinching line coordinates connecting line
Unit.
Optionally, in the step S30, also using the revised stratum top bottom surface rule of tomography coordinate pair
Grid carries out the process of tomography binding；In the step S40, in addition to by the way that thickness regular grid is retained and picked
Remove, at least so that grid quantity, coordinate and the distribution mode of the thickness regular grid and the revised stratum top bottom surface
Regular grid is consistent.
Optionally, the process of tomography binding includes：
S321：Fault line is obtained according to the tomography coordinate,
S322：Two nearest network computations of end points of detachment layer line are moved to endpoint location；
S323：Judge other network computations of near fault lines and the vertical range of the fault line, select wherein vertical
Closest some network computations, by it along being moved to perpendicular to the direction of fault line on the fault line, form tomography and tie up
The stratum top bottom surface regular grid after fixed.
Optionally, the step S323 further comprises：According to the linear equation of fault line, tomography is linearly searched for
The grid unit passed through, then each node of these grid units and fault line are carried out to the judgement of vertical range one by one, so
Afterwards, the minimum some nodes of vertical range are selected, for wherein each node, utilized the vertical line of the fault line of the node
The former coordinate of the node is replaced in the coordinate of the intersection point point of the fault line, the coordinate for completing all some nodes is replaced, handled
Afterwards, tomography binding is completed.
Optionally, the depth assignment in the step S5 is specially：Using described in the thickness determination in thickness regular grid
The corresponding modeling depth of three-dimensional Corner-point Grids, and then according to the depth modelling.
Optionally, in the S5：
If the thickness of grid unit is more than the minimum thickness unit of modeling demand, first by the thickness regular grid edge
Thickness direction splits into some cell cubes in units of the minimum thickness unit, then in modeling, in net corresponding to top surface
Modeling is formed and the cell cube quantity, some three-dimensional grids of thickness identical under lattice unit or corresponding to bottom surface on grid cell
Unit；
Optionally, in the S5：
If the thickness of grid unit is not more than the minimum thickness unit of modeling demand, in plane angle point corresponding to top surface
Modeling forms thickness and minimum thickness list under the grid cell of grid or corresponding to bottom surface on the grid cell of plane Corner-point Grids
First identical three-dimensional grid unit.
Optionally, in the step S5, in addition to attribute tax is carried out to the three-dimensional Corner-point Grids after depth assignment
The process of value.
Optionally, the attribute assignment comprises at least assigns to each grid cell porosity in the three-dimensional Corner-point Grids
The process of value and permeability assignment.
The present invention seek to based on existing data format and grid in field, propose a kind of simplification, easily
In realization and the oil-gas reservoir grid method for building up of numerical simulation can be quickly applied to, while general stratum can also be extended to
Or the quick three-dimensional modeling of rock stratum.
The elevation data of stratum top bottom interface of the present invention first in borehole data, is generated a certain by interpolation algorithm
The regular grid of survey region stratum top bottom surface, meanwhile, enter row interpolation using the formation thickness data in borehole data, generate
The regional stratum thickness regular grid, keep the consistent of grid quantity；Then according to the boundary coordinate and stratal pinch-outs in research area
Line number evidence, border file is established, reject the grid beyond border in grid file；Secondly, bottom is pushed up to stratum according to layer data
Face regular grid carries out tomography binding；It is finally based on two kinds of regular grids and realizes that the three-dimensional on stratum is built using Corner-point Grids technology
Mould.
Data format of the present invention has advantage.Stratum top bottom surface, formation thickness regular grid are turned to three
The data pre-processing of stratum modeling is tieed up, can be carried out by existing interpolation method and software, use the public affairs of ground popular software
Data format is opened, existing mature technology can be used, quickly realizes that data prepare；The regularization Grid square of stratum top bottom surface
In itself contained the absolute altitude (Z coordinate) and geodetic coordinates (X, Y-coordinate) on stratum, can be with to the three-dimensional geological model of foundation
Reach with true stratum carry out corresponding to effect.
The present invention is skillfully constructed using regular grid and the advantage of Corner-point Grids technology, modeling method.
Rule-based grid carries out Corner-point Grids in the present invention, can be carried out, subtracted using computer program rapid batch
Few manual intervention, improves modeling efficiency, and the workload of Man Graphics is greatly reduced.
Brief description of the drawings
Fig. 1 is the three dimensional contour line method of rule-based grid and Corner-point Grids technology in many embodiments of the present invention
Schematic flow sheet；
Fig. 2 is the schematic top plan view of top bottom surface regular grid in stratum in one embodiment of the invention；
Fig. 3 is the 3-D view of top bottom surface regular grid in stratum in one embodiment of the invention；
Fig. 4 is the amendment schematic diagram of top bottom surface regular grid in stratum in one embodiment of the invention；
Fig. 5 is the schematic diagram of one embodiment of the invention interrupting layer binding；
Fig. 6 is the schematic three dimensional views of bottom surface regular grid in one embodiment of the invention；
Fig. 7 is the schematic three dimensional views of three-dimensional Corner-point Grids in one embodiment of the invention；
Fig. 8 is the schematic diagram of threedimensional model in one embodiment of the invention；
Fig. 9 is the research area border schematic diagram in one embodiment of the invention；
Figure 10 is threedimensional model schematic diagram corresponding with Fig. 9.
Embodiment
Three-dimensional formation below with reference to Fig. 1 to Figure 10 to rule-based grid provided by the invention and Corner-point Grids technology
Modeling method is described in detail, wherein being in optional embodiment of the present invention, it is believed that those skilled in the art exist
Do not change spirit of the invention and it is modified and polished in the range of content.
It is pointed out that following examples content has write implementation process, the actually realization of this partial content exactly still
Need to carry out programming realization by OpenGL or other 3 d graphic libraries, mode is exactly to write to realize program in advance, program fortune
During row input corresponding to data, after the computing of backstage can display model, be desirable to clearly come out step herein.Other OpenGL
Etc simply some class libraries or function library, it is necessary to write the display that routine call these functions carry out Controlling model.Meanwhile this hair
The realization of bright scheme can be done based on OpenGL, and certainly other shape libraries can also be realized certainly, so being not limited to
OpenGL。
Fig. 1 is refer to, the invention provides the three dimensional contour line side of a kind of rule-based grid and Corner-point Grids technology
Method, including：
S1：For the research area of three dimension strata model to be established, extraction stratigraphic boundary, pinching line coordinates, tomography coordinate and
Borehole data；
S2：Stratum top bottom surface elevation data and formation thickness data are obtained according to the borehole data；Respectively enter step
S30 and step S40；
Two above step can be regarded as the Data Preparation Process of the present invention, and in the present invention, data prepare to be basis.
It is selected to establish after the research area of three dimension strata model, it is necessary first to extract three class master datas：Stratigraphic boundary/stratal pinch-outs
Line coordinates, tomography coordinate and borehole data.Also need to extract from borehole data purpose stratum top bottom surface elevation data, thickness
Degrees of data.If modeled for coal seam, the top bottom surface elevation data on purpose stratum can use seat earth contour line data.
S30：Row interpolation is entered to the stratum top bottom surface elevation data, obtains stratum top bottom surface regular grid；
Fig. 2 and Fig. 3 are refer to, row interpolation, regular grid such as Fig. 2 and Fig. 3 of generation are entered to stratum top bottom surface elevation data
It is shown.For example using the seat earth elevation data in certain research area's borehole data, using Kriging interpolation or other are conventional
The regular grid (20 × 20 regular grid as shown in Figure 2) of interpolation method generation seat earth, grid spacing can be according to number
Value simulation needs are set, four angle points of each grid unit respectively to should unit geodetic coordinates X, Y, in unit
It is worth for the seat earth absolute altitude value Z in the coordinates regional.If being modeled for coal seam, also seat earth contour line data can be used to enter
Row gridding.
S31：It is modified using stratum top bottom surface regular grid described in stratigraphic boundary and pinch-out boundary coordinate pair；
Fig. 4 is refer to, when being modified in the step S31 to the stratum top bottom surface regular grid, described in removal
In the top bottom surface regular grid of stratum outside pinch-out boundary coordinate connecting line, and it is disjoint with the pinching line coordinates connecting line
Grid unit.It is in practical operation, specifically can be described as reading border/stratal pinch-outs coordinate file, to stratum top bottom surface rule
Grid is modified.Repairing timing, the grid unit intersected with stratigraphic boundary/stratal pinch-outs coordinate line are retained.
S32：Tomography binding is carried out using the revised stratum top bottom surface regular grid of tomography coordinate pair, obtains plane
Corner-point Grids；
Fig. 5 is refer to, the process of tomography binding can include：
S321：Fault line is obtained according to the tomography coordinate,
S322：Two nearest network computations of end points of detachment layer line are moved to endpoint location；
S323：Judge other network computations of near fault lines and the vertical range of the fault line, select wherein vertical
Closest some network computations, by it along being moved to perpendicular to the direction of fault line on the fault line, form tomography and tie up
The stratum top bottom surface regular grid after fixed.
The step S323 further comprises：According to the linear equation of fault line, linearly search for tomography and passed through
Grid unit, then each node of these grid units and fault line are carried out to the judgement of vertical range, then, selection one by one
The minimum some nodes of vertical range, for wherein each node, utilized the vertical line of the fault line of the node described
The coordinate of the intersection point point of fault line replaces the former coordinate of the node, and the coordinate for completing all some nodes is replaced, after processing, is completed
Tomography is bound.
Specifically, tomography binding (when studying in area without tomography, skipping this step).Regular grid structure list in the plane
First, manifestation mode is simple, it is impossible to expresses fault information, it is impossible to embodies the excellent of Corner-point Grids expression complicated geological object constructional features
Gesture.Tomography is expressed, it is necessary to which tomography is tied in quadrilateral mesh by force in Corner-point Grids model (plane quadrilateral), make
Some Grid Edges therein turn into tomography side.The mesh point of near Faults, is exactly moved on fault line by so-called tomography binding,
To reach the purpose in fault mosaic to grid.The principle of binding is：It is consistent nearby, i.e., first the end points of every tomography is tied to
In grid, then the point on fault line is tied in nearest grid successively along fault line, until completing.Fig. 5 is rule
Grid binding fault line forms the schematic diagram of Corner-point Grids, and p1, p2 are two end points of fault line in Fig. 5 (a), binding method tool
Body is as follows：
The first step：The mesh point nearest apart from tomography head and the tail node is accommodated to first and last node.It is as full in searched in Fig. 5 (a)
The mesh node of sufficient condition is 0,7, then the two nodes is respectively moved into p1, p2 position；
Second step：According to fault line p1p2 linear equation, the grid unit passed through along p1p2 directions search tomography simultaneously adds
With mark, each node of these grid units and fault line are subjected to Distance Judgment one by one, the minimum node coordinate of distance is used
Coordinate on fault line is replaced, and changes the node coordinate of grid unit.As in grid unit A four nodes in Fig. 5 (a) except
0 has moved to outside end points p1 that remaining three nodes and tomography are closest for point 1, and the vertical line intersection point of point 1 will be crossed on straight line
Coordinate of the coordinate of point as new mesh node, four grid being attached thereto can deform upon.
3rd step：Node coordinate processing is carried out to the markd grid unit of institute, forms plane Corner-point Grids
S40：Row interpolation is entered to the formation thickness data, obtains thickness regular grid, and by thickness regular grid
Retained and rejected, at least so that the grid quantity of the thickness regular grid, distribution mode and the revised stratum
Top bottom surface regular grid is consistent；
Thickness degrees of data enters row interpolation over the ground, create-rule grid, and Method And Principle is consistent with step S30.Needed during gridding
Ensure that number of grid in the x, y direction is consistent with the revised stratum top bottom surface regular grid.In numerical simulation, net
Lattice minimum unit typically has limit value, but formation thickness wouldn't be entered by the requirement of minimum grid unit when regular grid is established
Row is divided, and the stratum total thickness value in current grid unit coordinates regional, follow-up three-dimensional angle point are deposited in formation thickness grid
Divided again when grid is established.Thickness grid is retained and rejected according to stratum top bottom surface border grid unit,
Make final number of grid consistent.
About " thickness regular grid " and " stratum top bottom surface regular grid ", using the present invention modeling method when, carry out
Stratum modeling needs to ensure that the quantity palpus of the two grid is consistent, coordinate is consistent.This is based on, current scheme interrupting layer is tied up
Surely can be in any one progress of two grid, the coordinate of another grid makes corresponding modification, after so bringing
Fruit is that the display effect of tomography is relatively simple.In the possible improvement of the present invention tomography may be carried out respectively for two grid
Binding, so as to improve the display effect in the three-dimensional model of tomography.And in the elaboration of foregoing embodiments, mainly with to described
The process that stratum top bottom surface regular grid carries out tomography binding is illustrated, in optional scheme of the invention, in the step
In S30 and S40, in addition to：
Tomography binding is carried out at least one the stratum top bottom surface regular grid and thickness regular grid, while causes institute
The quantity of stratum top bottom surface regular grid and thickness regular grid, coordinate, arrangement mode is stated to be consistent；
S5：Depth assignment is carried out to the stratum top bottom surface regular grid using the thickness regular grid, so as to carry out
Modeling, obtains three-dimensional Corner-point Grids；
In the most of embodiments of the present invention, the depth assignment in the step S5 is specially：Utilize thickness regular grid
In thickness determine the corresponding modeling depth of the three-dimensional Corner-point Grids, and then according to the depth modelling.
Furthermore, it is understood that in optional embodiment of the present invention, in the S5：
If the thickness of grid unit is more than the minimum thickness unit of modeling demand, first by the thickness regular grid edge
Thickness direction splits into some cell cubes in units of the minimum thickness unit, then in modeling, in net corresponding to top surface
Modeling is formed and the cell cube quantity, some three-dimensional grids of thickness identical under lattice unit or corresponding to bottom surface on grid cell
Unit；
If the thickness of grid unit is not more than the minimum thickness unit of modeling demand, in plane angle point corresponding to top surface
Modeling forms thickness and minimum thickness list under the grid cell of grid or corresponding to bottom surface on the grid cell of plane Corner-point Grids
First identical three-dimensional grid unit.
In optional embodiment of the present invention, in the step S5, in addition to the three dimensional angular after depth assignment
Dot grid carries out the process of attribute assignment.The attribute assignment is comprised at least to each grid list in the three-dimensional Corner-point Grids
The process of first porosity assignment and permeability assignment.
Specifically, comprise the following steps：
The first step：When formation thickness is more than minimum unit requirement thickness, formation thickness grid unit is marked and thick
Degree is split.For example the formation thickness value in certain grid unit is 15m, the minimum unit thickness of modeling demand is 3m, then marks this
Unit is thickness split cells, and it is split as into 5 cell cubes in Z-direction, i.e., final three-dimensional formation Corner-point Grids Z-direction can
It is of five storeys.
Second step：Based on stratum top bottom surface absolute altitude regular grid, Corner-point Grids assignment is carried out.
1. in the modeling of three-dimensional Corner-point Grids, if the top bottom surface absolute altitude regular grid knot using a certain stratum of COORD arrays
Point coordinates, recorded by the row of stratum top bottom surface absolute altitude regular grid；The depth of Corner-point Grids is recorded using ZCORN arrays,
The Z coordinate of i.e. all Corner-point Grids units, arranged in the way of the row, column, layer by stratum regular grid；It can be used in addition
PROPERTY arrays, the POROSITY arrays class record attribute information of Corner-point Grids unit, porosity value, using PERMX,
PERMY, PERMZ array record the permeability in the X, Y, Z direction of grid cell respectively.
2. carrying out assignment to stratum top bottom surface regular grid unit mainly includes 2 parts：The depth of Corner-point Grids unit is assigned
Value and attribute assignment.Stratum top bottom surface regular grid unit is traveled through, if formation thickness is not split, by its grid cell
ZCORN values are set to the value of the corresponding formation thickness regular grid unit of COORD values therewith；If formation thickness is split, if
Using stratum top surface grid, then ZCORN values are done into downwards multiple extension, then upwardly extended if bottom surface, formation is multiple to have phase
Unit grid with COORD values stacks.To its attribute assignment, such as porosity, permeability assignment while traveling through grid unit.
S6：Threedimensional model is obtained using the three-dimensional Corner-point Grids.
The present invention is specifically addressed with a concrete practice scheme below in conjunction with Fig. 6 to Figure 10, the block geological structure
Simply, seam inclination is small, and seat earth relief intensity is small, average coal seam thickness 3.15m.Use the seat earth absolute altitude in block
Data carry out gridding and obtain No. 3 regular grid data of seat earth 100 × 75, and its grid three-dimensional surface can be shown as Fig. 6.
Assuming that grid cell Z-direction minimum value is 5m during numerical simulation, can Xu Dui areas according to coal seam average thickness in research area
Coal seam thickness is split in domain, i.e., final three-dimensional Corner-point Grids are Z-direction single layer mesh.It is right according to the step of realizing of the present invention
Seat earth grid carries out tomography binding and thickness assignment, obtains Fig. 7.
Fig. 7 Corner-point Grids are carried out with attribute assignment to go forward side by side style of writing reason stick picture disposing, can obtain three-dimensional coal seam mould true to nature
Type, as shown in Figure 8.
Fig. 9 and Figure 10 also contrasts the schematic diagram for showing and being modified using block border to grid.
Existing three dimensional contour line method is artificial progress graphic plotting mostly, similar to be drawn using AutoCAD, than
Capture function as the line between point and point may need user to be provided by software is connected manually, for the ground of complexity
For layer modeling, the workload of this work is very big, and easily malfunctions.In more advanced software, brill can be first imported
Hole data, the three-dimensional curved surface (each face can be considered a stratigraphic model) of interpolation generation, by manually making between multiple curved surfaces
Fixed rule is overlapped, connected, Boolean calculation, and so as to generate the model of three-dimensional, workload is also very big, error rate
It is that comparison is high.The present invention uses what is write in advance primarily directed to the manual operation part in above-mentioned several three-dimensional modelings
Code regulation, by the attribute of coordinate, grid, row into threedimensional model, has similar use " stretching " to order automatically on backstage in program operation
The meaning of the threedimensional model formed is made, thus should can reduce the intensity of hand labor.Moreover, from manually to automatic
Change also necessarily brings whole work-based logic and flow to change, so, those skilled in the art can not possibly be according to artificial work
The method of work just can obviously obtain the workflow of the present invention.
In summary, the elevation data of stratum top bottom interface of the present invention first in borehole data, is calculated by interpolation
Method generates the regular grid of a certain survey region stratum top bottom surface, meanwhile, entered using the formation thickness data in borehole data
Row interpolation, the regional stratum thickness regular grid is generated, keep the consistent of grid quantity；Then according to the boundary coordinate in research area
With stratal pinch-outs line number evidence, border file is established, rejects the grid beyond border in grid file；Secondly, according to layer data
Tomography binding is carried out to stratum top bottom surface regular grid；It is finally based on two kinds of regular grids and realizes stratum using Corner-point Grids technology
Three-dimensional modeling.
Data format of the present invention has advantage.Stratum top bottom surface, formation thickness regular grid are turned to three
The data pre-processing of stratum modeling is tieed up, can be carried out by existing interpolation method and software, use the public affairs of ground popular software
Data format is opened, existing mature technology can be used, quickly realizes that data prepare；The regularization Grid square of stratum top bottom surface
In itself contained the absolute altitude (Z coordinate) and geodetic coordinates (X, Y-coordinate) on stratum, can be with to the three-dimensional geological model of foundation
Reach with true stratum carry out corresponding to effect.
The present invention is skillfully constructed using regular grid and the advantage of Corner-point Grids technology, modeling method.
Rule-based grid carries out Corner-point Grids in the present invention, can be carried out, subtracted using computer program rapid batch
Few manual intervention, improves modeling efficiency, and the workload of Man Graphics is greatly reduced.
Claims (10)
1. a kind of three dimensional contour line method of rule-based grid and Corner-point Grids technology, it is characterised in that：Including：
S1：For the research area of three dimension strata model to be established, extraction stratigraphic boundary, pinching line coordinates, tomography coordinate and drilling
Data；
S2：Stratum top bottom surface elevation data and formation thickness data are obtained according to the borehole data；Respectively enter step S30 and
Step S40；
S30：Row interpolation is entered to the stratum top bottom surface elevation data, obtains stratum top bottom surface regular grid；Then stratum is utilized
Stratum top bottom surface regular grid is modified described in border and pinch-out boundary coordinate pair；
S40：Row interpolation is entered to the formation thickness data, obtains thickness regular grid；
In the step S30 and S40, in addition to：
Tomography binding is carried out at least one the stratum top bottom surface regular grid and thickness regular grid, while is caused describedly
The quantity of layer top bottom surface regular grid and thickness regular grid, coordinate, arrangement mode are consistent；
S5：Depth assignment is carried out to the stratum top bottom surface regular grid using the thickness regular grid, so as to be modeled,
Obtain three-dimensional Corner-point Grids；
S6：Threedimensional model is obtained using the three-dimensional Corner-point Grids.
2. the three dimensional contour line method of rule-based grid as claimed in claim 1 and Corner-point Grids technology, its feature exist
In：When being modified in the step S30 to the stratum top bottom surface regular grid, the stratum top bottom surface rule lattice are removed
In net outside pinch-out boundary coordinate connecting line, and the grid unit disjoint with the pinching line coordinates connecting line.
3. the three dimensional contour line method of rule-based grid as claimed in claim 1 and Corner-point Grids technology, its feature exist
In：In the step S30, also tomography is carried out using the revised stratum top bottom surface regular grid of tomography coordinate pair
The process of binding；In the step S40, in addition to by the way that thickness regular grid is retained and rejected, at least so that institute
Grid quantity, coordinate and the arrangement mode for stating thickness regular grid are consistent with the revised stratum top bottom surface regular grid.
4. the three dimensional contour line method of rule-based grid as claimed in claim 3 and Corner-point Grids technology, its feature exist
In：The process of tomography binding includes：
S321：Fault line is obtained according to the tomography coordinate,
S322：Two nearest network computations of end points of detachment layer line are moved to endpoint location；
S323：Judge other network computations of near fault lines and the vertical range of the fault line, selection wherein vertical range
Nearest some network computations, by it along being moved to perpendicular to the direction of fault line on the fault line, after forming tomography binding
The stratum top bottom surface regular grid.
5. the three dimensional contour line method of rule-based grid as claimed in claim 4 and Corner-point Grids technology, its feature exist
In：The step S323 further comprises：According to the linear equation of fault line, the grid that tomography is passed through linearly is searched for
Unit, then each node of these grid units and fault line are carried out to the judgement of vertical range one by one, then, selection vertically away from
From some nodes of minimum, for wherein each node, the vertical line of the fault line of the node was utilized in the fault line
The coordinate of intersection point point replace the former coordinate of the node, the coordinate for completing all some nodes is replaced, after processing, is completed tomography and is tied up
It is fixed.
6. the three dimensional contour line method of rule-based grid as claimed in claim 1 and Corner-point Grids technology, its feature exist
In：Depth assignment in the step S5 is specially：The three-dimensional Corner-point Grids are determined using the thickness in thickness regular grid
Corresponding modeling depth, and then according to the depth modelling.
7. the three dimensional contour line method of the rule-based grid and Corner-point Grids technology as described in claim 1 or 6, its feature
It is：In the S5：
If the thickness of grid unit is more than the minimum thickness unit of modeling demand, first by the thickness regular grid along thickness
Direction splits into some cell cubes in units of the minimum thickness unit, then in modeling, in grid list corresponding to top surface
Modeling is formed and the cell cube quantity, some three-dimensional grid lists of thickness identical under member or corresponding to bottom surface on grid cell
Member.
8. the three dimensional contour line method of rule-based grid as claimed in claim 7 and Corner-point Grids technology, its feature exist
In：In the S5：
If the thickness of grid unit is not more than the minimum thickness unit of modeling demand, in plane Corner-point Grids corresponding to top surface
Grid cell under or bottom surface corresponding to plane Corner-point Grids grid cell on modeling form thickness and minimum thickness unit phase
Same three-dimensional grid unit.
9. the three dimensional contour line method of rule-based grid as claimed in claim 1 and Corner-point Grids technology, its feature exist
In：In the step S6, in addition to the three-dimensional Corner-point Grids after depth assignment are carried out with the process of attribute assignment.
10. the three dimensional contour line method of rule-based grid as claimed in claim 9 and Corner-point Grids technology, its feature exist
In：The attribute assignment comprises at least to be assigned to each the grid cell porosity assignment and permeability in the three-dimensional Corner-point Grids
The process of value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510631247.8A CN105184867B (en) | 2015-09-29 | 2015-09-29 | The three dimensional contour line method of rule-based grid and Corner-point Grids technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510631247.8A CN105184867B (en) | 2015-09-29 | 2015-09-29 | The three dimensional contour line method of rule-based grid and Corner-point Grids technology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105184867A CN105184867A (en) | 2015-12-23 |
| CN105184867B true CN105184867B (en) | 2017-11-10 |
Family
ID=54906917
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510631247.8A Expired - Fee Related CN105184867B (en) | 2015-09-29 | 2015-09-29 | The three dimensional contour line method of rule-based grid and Corner-point Grids technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105184867B (en) |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107991702B (en) * | 2016-10-26 | 2019-05-17 | 中国石油化工股份有限公司 | The simple accurate determination method of normal fault timing |
| CN108074282B (en) * | 2016-11-08 | 2021-02-23 | 中国石油化工股份有限公司 | Thickness correction method of physical property parameters based on three-dimensional geological modeling software |
| CN106652033B (en) * | 2016-12-05 | 2020-02-14 | 中国石油天然气股份有限公司 | Method for subdividing natural grid of geological profile |
| CN106897445A (en) * | 2017-03-01 | 2017-06-27 | 武汉科技大学 | It is a kind of that the novel method that operation delimited on border is improved based on ArcGIS |
| CN109215124A (en) * | 2018-09-25 | 2019-01-15 | 中国科学院武汉岩土力学研究所 | The construction method of underground engineering 3D grid model under a kind of complex geological condition |
| CN110838181B (en) * | 2019-11-08 | 2023-09-26 | 中煤航测遥感集团有限公司 | Horizontal cutting method, device and equipment for triangular net model and storage medium |
| CN111127644B (en) * | 2019-12-12 | 2021-03-02 | 北京龙软科技股份有限公司 | Method and device for drawing coal seam thickness contour line of open pit coal mine |
| CN111161403A (en) * | 2019-12-20 | 2020-05-15 | 核工业北京地质研究院 | Method for acquiring top and bottom surfaces of sedimentary stratum based on spatial interpolation of drilling data |
| CN111401451B (en) * | 2020-03-17 | 2023-06-30 | 中国石油化工股份有限公司 | Automatic identification method for reservoir configuration interface in reservoir geological model |
| CN111612903B (en) * | 2020-04-29 | 2023-03-21 | 中冶沈勘工程技术有限公司 | Geological data visualization method based on mixed data model |
| CN114114405A (en) * | 2020-08-26 | 2022-03-01 | 中国石油天然气股份有限公司 | Stratum true thickness determination method and device based on discrete data gridding |
| CN112461728B (en) * | 2020-11-12 | 2021-07-02 | 成都理工大学 | Method and system for determining whole core permeability and pore throat radius distribution |
| CN112634450A (en) * | 2021-01-06 | 2021-04-09 | 安徽理工大学 | Method for integrally constructing local complex fault through corner grid |
| CN113409463B (en) * | 2021-06-29 | 2022-06-07 | 中国地质大学（武汉） | Three-dimensional geological model construction method and device including pinch-out treatment |
| CN117671160B (en) * | 2024-01-30 | 2024-04-23 | 北京星天地信息科技有限公司 | Multi-source data collaborative coal seam modeling method, device, equipment and storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102194253A (en) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | Method for generating tetrahedron gridding for three-dimensional geological structure |
| CN102194252A (en) * | 2011-05-17 | 2011-09-21 | 北京航空航天大学 | Geological-stratum-structure-based method for generating triangular lattice grids |
| CN103514631A (en) * | 2013-10-21 | 2014-01-15 | 重庆市勘测院 | Gridding type parallel construction method for three-dimensional geological model based on massive borehole data |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102612682B (en) * | 2009-11-12 | 2016-04-27 | 埃克森美孚上游研究公司 | For the method and apparatus of reservoir modeling and simulation |
-
2015
- 2015-09-29 CN CN201510631247.8A patent/CN105184867B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102194252A (en) * | 2011-05-17 | 2011-09-21 | 北京航空航天大学 | Geological-stratum-structure-based method for generating triangular lattice grids |
| CN102194253A (en) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | Method for generating tetrahedron gridding for three-dimensional geological structure |
| CN103514631A (en) * | 2013-10-21 | 2014-01-15 | 重庆市勘测院 | Gridding type parallel construction method for three-dimensional geological model based on massive borehole data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105184867A (en) | 2015-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105184867B (en) | The three dimensional contour line method of rule-based grid and Corner-point Grids technology | |
| US10339709B1 (en) | Systems and methods for coordinated editing of seismic data in dual model | |
| Guo et al. | Three-dimensional geological modeling and spatial analysis from geotechnical borehole data using an implicit surface and marching tetrahedra algorithm | |
| CN111899341B (en) | Three-dimensional geological rapid modeling method for long and large banded linear engineering | |
| Natali et al. | Modeling Terrains and Subsurface Geology. | |
| US12061306B2 (en) | Constructing structural models of the subsurface | |
| Wu et al. | An approach to computer modeling and visualization of geological faults in 3D | |
| CN104574511B (en) | A kind of quick progressive three-dimensional geological modeling method | |
| CN105469443B (en) | Based on geotraverse（PRB）The method that process Data fuse generates three-dimensional geological figure | |
| CN102609982B (en) | Topology discovery method of space geological data based on unstructured mode | |
| CN103279986A (en) | Three-dimensional horizontal geological section map manufacturing method and applications thereof | |
| CA2713019A1 (en) | Gridless geological modeling | |
| CN107944086A (en) | A kind of fast modeling method based on borehole data | |
| GB2474740A (en) | Gridless geological modeling of a structural framework | |
| CN105593908A (en) | Identifying matching properties between a group of bodies representing a geological structure and a table of properties | |
| CN106709987B (en) | Dynamic construction method of three-dimensional geological profile model | |
| Wu et al. | On three-dimensional geological modeling and visualization | |
| CN110990919A (en) | Three-dimensional GIS route selection design analysis method for urban rail transit | |
| CN112231818B (en) | Mountain city three-dimensional geological modeling method based on BIM | |
| CN111445569B (en) | Sedimentary geological evolution dynamic simulation method | |
| CN103218850A (en) | True three-dimensional mining blast unit body modeling method | |
| CN117173357A (en) | Mine three-dimensional geological modeling and layered cutting method | |
| CN106097446A (en) | A kind of three dimensional geological curved surface method for reconstructing based on extensive geological data | |
| CN115035258A (en) | Efficient urban three-dimensional geological modeling method based on CAD (computer-aided design) drilling histogram | |
| CN106097447A (en) | A kind of curve reestablishing method of extensive geological data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171110 Termination date: 20180929 |