CN103236088B - Three-dimensional geological structure modeling method and system based on composite grid - Google Patents

Three-dimensional geological structure modeling method and system based on composite grid Download PDF

Info

Publication number
CN103236088B
CN103236088B CN201310130973.2A CN201310130973A CN103236088B CN 103236088 B CN103236088 B CN 103236088B CN 201310130973 A CN201310130973 A CN 201310130973A CN 103236088 B CN103236088 B CN 103236088B
Authority
CN
China
Prior art keywords
unit
grid
tetrahedral grid
tetrahedral
initial
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.)
Active
Application number
CN201310130973.2A
Other languages
Chinese (zh)
Other versions
CN103236088A (en
Inventor
赵明伟
袁钢辉
矫树春
齐建军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BEIJING GOLDENSUN PETROLEUM TECHNOLOGIES INC
Original Assignee
BEIJING GOLDENSUN PETROLEUM TECHNOLOGIES INC
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 BEIJING GOLDENSUN PETROLEUM TECHNOLOGIES INC filed Critical BEIJING GOLDENSUN PETROLEUM TECHNOLOGIES INC
Priority to CN201310130973.2A priority Critical patent/CN103236088B/en
Publication of CN103236088A publication Critical patent/CN103236088A/en
Application granted granted Critical
Publication of CN103236088B publication Critical patent/CN103236088B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention embodiment discloses a three-dimensional geological structure modeling method and system based on a composite grid. The three-dimensional geological structure modeling method includes selecting a grid system area of structure modeling, selecting an origin of the grid system and the grid step size and number on the coordinate direction, and generating an initial hexahedral grid system; dividing each initial hexahedral grid unit into multiple tetrahedral grid units, and generating a tetrahedral grid system; numbering the multiple tetrahedral grid units, recording the number of the multiple tetrahedral grid units contained in each initial hexahedral grid unit, and generating topology relationship of the tetrahedral grid system; adding an actual boundary of a set geological body model, redividing the tetrahedral grid system, and generating a composite grid system; and adding well location data and cross section data of the three-dimensional geological body model into the composite grid system, and performing layer tracking and three-dimensional geological grid construction.

Description

Based on three-dimensional geologic structure modeling method and the system of composite gridding
Technical field
The present invention relates to computational geometry and three-dimensional geological modeling field, particularly relate to a kind of three-dimensional geologic based on composite gridding structure modeling method and system, a kind of method and system generating complementary two nested grid systems concretely, landform complicated arbitrarily can be simulated, and grid fast query and the location of given coordinate position can be carried out.
Background technology
The generation technique of computing grid always is the important branch of oil exploration and development fields, and it has established the basis of a lot of work, such as: isoline algorithm, and three-dimensional geologic structure modeling, numerical reservoir simulation etc.Computing grid is divided into two large classes generally: structured grid and unstrctured grid.Conventional structured grid comprises rectangular node, hexahedral mesh etc., and conventional unstrctured grid then comprises triangular mesh, tetrahedral grid etc.For fairly simple geometric shape process, it is fairly perfect that structured grid has developed at present, and the effect of application is also all well and good.But for the geometric shape of more complicated, conventional structured grid often runs into very large difficulty when processing these profiles, and the effect brought is often also barely satisfactory.The proposition of unstructured grid and application, successfully solve the above-mentioned shortcoming of structured grid, greatly advances the application of computing method at petroleum exploration and development.Owing to not needing the structural and orthogonality considering grid in the generation of unstrctured grid, the size shape of grid cell and the position of net point is just made to be more prone to control, therefore unstrctured grid has just had very large dirigibility, and the description for complex appearance there has also been universality.Just because of there is this series of advantage, unstrctured grid technology receives the extensive attention of domestic and international oil exploration and development fields, and obtains and develop rapidly and apply.
But such as isoline algorithm in actual application, three-dimensional geological modeling, numerical reservoir simulation etc., the most basic requirement is that generated grid perfectly must express any given region, and this just means that it is not all right for being used alone structured grid; Need again in generated grid system, add new well location/layer data simultaneously, wherein related to grid fast query and the location of given coordinate position, and be used alone the obvious inefficiency of unstrctured grid.
In petroleum exploration and development, structured grid relatively conventional is at present mainly basic grid unit with rectangle in two-dimensional case, is mainly basic grid unit in three-dimensional case with hexahedron.The generation of structured grid is fairly simple, only needs the initial point determining grid system, along mesh spacing and the meshes number of all directions.The advantage of structured grid system is to generate very simple, and is easy to fast query and the location of coordinate position.But structured grid is very difficult to the simulation of complicated landform, this is because the unit (two-dimensional rectangle/three-dimensional hexahedron) of structured grid does not have very large dirigibility.When the border more complicated of given area, the effect that structured grid brings can be very poor.
In addition, in petroleum exploration and development, unstrctured grid relatively conventional is at present mainly basic grid unit with triangle in two-dimensional case, is mainly basic grid unit in three-dimensional case with tetrahedron.The generation method mainly Delaunay method of conventional unstrctured grid: be first covered with initial grid cell in space, then constantly introduce new point according to certain criterion, forms new grid cell by node restructuring.Here criterion refers to Delaunay sphere criterion, namely requires that the circumscribed circle of each triangular mesh unit does not comprise other nodes in two-dimensional case, namely requires that the circumsphere of each tetrahedral grid unit does not comprise other grid nodes in three-dimensional case.The advantage of Delaunay method is just that the speed of mess generation is fast, and mesh quality is high, can carry out two-dimensional/three-dimensional numerical evaluation well on institute's generating mesh basis.
But the shortcoming of Delaunay grid mainly contains 2 points: one is the unstrctured grid subdivision needing to be provided in advance given area.This is because require initial unstrctured grid subdivision in generation Delaunay grid process, then add some new points to realize the displacement to original point according to Delaunay sphere rule.In actual application, this will inevitably bring extra calculated amount, and the quality of initial mesh quality also can have impact to the speed of the final Delaunay grid generated, and this is uppity.Two is that generated Delaunay grid is unfavorable for Search and Orientation.In actual algorithm design, often need to search the grid meeting specified conditions and number and position, such as: the well newly added in which grid or tomography pass which grid.Due to Delaunay coarse gridding is the topological structure of point, line, surface, body, and must depend on circulation when carrying out the searching of given coordinate position, this greatly can have influence on the efficiency of involved algorithm.
Summary of the invention
The object of the invention is the deficiency in order to overcome structured grid and the unstrctured grid existed in prior art, there is provided one can simulate any complicated landform, and formation speed is fast, quality is good and be convenient to carry out the fast query of given coordinate points and the composite gridding of location, for the structure modeling of three-dimensional geologic.
In order to achieve the above object, the embodiment of the invention discloses a kind of three-dimensional geologic based on composite gridding structure modeling method, described method comprises: steps A, the grid system region of selected structure modeling, and mesh spacing in described grid system region on the initial point of selected grid system and coordinate direction and number, generate initial hexahedral mesh system, and record the numbering of each initial hexahedral mesh unit in described initial hexahedral mesh system; Described each initial hexahedral mesh dividing elements is multiple tetrahedral grid unit by step B, generates tetrahedral grid system; Step C, is numbered described multiple tetrahedral grid unit, records the numbering of the tetrahedral grid unit comprised in described each hexahedral mesh unit, generates the topological relation of described tetrahedral grid system; Step D, adds the actual boundary of the model of geological structure body of setting, carries out subdivision again to described tetrahedral grid system, generates composite gridding system; Well location data in described three-dimensional geological body Model and profile data are joined in described composite gridding system, carry out the structure of layer auto-tracking and three-dimensional geological volume mesh by step e.
In order to achieve the above object, the embodiment of the invention also discloses a kind of three-dimensional geologic based on composite gridding structure modeling, comprise: initial hexahedral mesh system generation unit, for the grid system region of selected structure modeling, and mesh spacing in described grid system region on the initial point of selected grid system and coordinate direction and number, generate initial hexahedral mesh system, and record the numbering of each initial hexahedral mesh unit in described initial hexahedral mesh system; Tetrahedral grid system generation unit, for being multiple tetrahedral grid unit by described each initial hexahedral mesh dividing elements, generates tetrahedral grid system; Topology generation unit, for being numbered described multiple tetrahedral grid unit, recording the numbering of the tetrahedral grid unit comprised in described each hexahedral mesh unit, generating the topological relation of described tetrahedral grid system; Composite gridding system generation unit, adds the actual boundary of the model of geological structure body of setting, carries out subdivision again to described tetrahedral grid system, generates composite gridding system; Three-dimensional geologic grid enable unit, for the well location data in described three-dimensional geological body Model and profile data being joined in described composite gridding system, carries out the structure of layer auto-tracking and three-dimensional geological volume mesh.
Three-dimensional geologic based on the composite gridding structure modeling method of the embodiment of the present invention and system, can generate a kind of composite gridding system with complementary two nested grid systems, can simulate landform complicated arbitrarily.Further, input needs the region of modeling, can generate the three-dimensional composite gridding system in this region fast.Given arbitrary coordinate position (being positioned at region), can the numbering of tetrahedral grid unit at this coordinate place of quick position.This composite gridding system can add the new node (coordinate is positioned at given area) of any amount very easily, only need carry out Local grid subdivision (at new coordinate points place) and topological relation renewal to original grid system.Therefore, the present invention can also realize grid fast query and the location of given coordinate position.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those skilled in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.
Fig. 1 is the method flow diagram of the structure of the three-dimensional geologic based on the composite gridding modeling method of the embodiment of the present invention;
Fig. 2 is the embodiment illustrated in fig. 1 method flow diagram adding well location data;
Fig. 3 is the structural representation of the structure of the three-dimensional geologic based on the composite gridding modeling of the embodiment of the present invention;
Fig. 4 be embodiment illustrated in fig. 3 in the structural representation of composite gridding system generation unit 104;
Fig. 5 be embodiment illustrated in fig. 3 in the structural representation of geologic body grid enable unit 105;
Fig. 6 illustrates the method flow diagram carrying out the specific embodiment of three-dimensional geologic structure based on composite gridding of the present invention;
Fig. 7 be embodiment illustrated in fig. 6 in the schematic diagram of initial hexahedral mesh system;
Fig. 8 be embodiment illustrated in fig. 6 in single hexahedral mesh unit be the schematic diagram of multiple tetrahedral grid unit by subdivision;
Fig. 9 be embodiment illustrated in fig. 8 in gone out by subdivision five tetrahedral grid unit;
Figure 10 be embodiment illustrated in fig. 6 in the schematic diagram of tetrahedral grid system of generation;
Figure 11 is the schematic diagram of the point being added on tetrahedron inside;
Figure 12 is the schematic diagram of the point be added in tetrahedron one side;
Figure 13 is the schematic diagram of the point be added on tetrahedron one side;
Figure 14 adds the schematic diagram of starting point on the limit on a certain summit of tetrahedron;
Figure 15 is the schematic diagram adding the limit that starting point does not overlap with tetrahedron top;
Figure 16 is the schematic diagram in the inner interpolation face of tetrahedron.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.
The composite gridding of the embodiment of the present invention does not refer to that the final grid generated comprises SST hexahedral meshes and tetrahedron unstrctured grid simultaneously, and refer to complementary two nested grid systems, a set of is SST hexahedral meshes system, for covering the initial SST hexahedral meshes system in this region at given area rapid build; Another set of is tetrahedron unstrctured grid system, inherits from initial SST hexahedral meshes system, and the boundary in given area is revised this tetrahedron unstrctured grid system.What real participation calculated in actual applications is tetrahedron unstrctured grid system, and what SST hexahedral meshes system played is effect tetrahedron unstrctured grid system being carried out to fast query location.
Fig. 1 is the method flow diagram of the structure of the three-dimensional geologic based on the composite gridding modeling method of the embodiment of the present invention.As shown in the figure, the method comprises:
Step S101, the grid system region of selected structure modeling, and mesh spacing in described grid system region on the initial point of selected grid system and coordinate direction and number, generate initial hexahedral mesh system, and record the numbering of each initial hexahedral mesh unit in described initial hexahedral mesh system; Described each initial hexahedral mesh dividing elements is multiple tetrahedral grid unit by step S102, generates tetrahedral grid system; Step S103, is numbered described multiple tetrahedral grid unit, records the numbering of the tetrahedral grid unit comprised in described each hexahedral mesh unit, generates the topological relation of described tetrahedral grid system; Step S104, adds the actual boundary of the model of geological structure body of setting, carries out subdivision again to described tetrahedral grid system, generates composite gridding system; Well location data in described three-dimensional geological body Model and profile data are joined in described composite gridding system, carry out the structure of layer auto-tracking and three-dimensional geological volume mesh by step S105.
In the step S101 of the present embodiment, the grid system region of structure modeling is hexahedron/rectangular parallelepiped region.The initial mesh system that this step generates is SST hexahedral meshes, if determine grid system initial point, along the mesh spacing of all directions and meshes number.The generation of SST hexahedral meshes belongs to the common practise of this area, therefore repeats no more in the embodiment of the present invention.
In the step S102 of the present embodiment, each initial hexahedral mesh unit in described initial mesh system is divided into multiple tetrahedral grid unit according to unified order, generates tetrahedral grid system.In the present embodiment, initial mesh unit is hexahedron, can be divided into five tetrahedrons.
In the step S103 of the present embodiment, the topological relation of the described tetrahedral grid system of generation comprises: the shared relationship of joint points, edges, faces, body.Initial mesh unit is hexahedron, can be divided into five tetrahedrons, then its topological relation is the shared relationship of joint points, edges, faces and body.
In the step S104 of the present embodiment, add the actual boundary of the model of geological structure body of setting, again subdivision is carried out to described tetrahedral grid system, comprising: in multiple tetrahedral grid unit of described tetrahedral grid system, add point, limit or face, generate new tetrahedral grid unit.
For tetrahedron, when adding point operation, if this point is inner at tetrahedron, then connecting three summits of this point with tetrahedral four faces successively, forming four new tetrahedrons; If this point drops on the inside in tetrahedron face, then three summits first connecting this summit and this face form three line segments, then connect this summit and this face to point, form three new tetrahedrons; If this point drops on the inside on a tetrahedron limit, then two summits first connecting this summit and this limit form two line segments, then connect this summit with will not this summit coplanar to point, form two new tetrahedrons; If this point drops on certain summit of tetrahedron, then disregard.
For tetrahedron, when adding limit operation, if this limit overlaps with tetrahedral certain limit, then disregard; If this limit does not all overlap with tetrahedral arbitrary limit, then distinguish following two kinds of situations: if the summit on this limit is certain summit tetrahedral, the operation then adding limit just can transform into a particular ray (starting point by the summit on interpolation limit, direction is along the terminal on added limit) with the operation of this tetrahedron find intersection, the operation adding point is carried out to tried to achieve intersection point, then completes the operation this tetrahedron being added to limit.For sake of convenience, if represent this tetrahedron with Tetra. the summit on this limit is not certain summit of Tetra, then first ask particular ray (starting point by the summit on interpolation limit, direction is along the terminal on added limit) first the intersection point A crossing with Tetra, then antinode A carries out the operation adding point in Tetra, and such Tetra is just split into as N number of sub-tetrahedron.Travel through this N number of sub-tetrahedron, find the sub-tetrahedron of this ray and which crossing and try to achieve intersection points B (representing this sub-tetrahedron with STetra), and then antinode B carries out the operation adding point in STetra, then Tetra subdivision completes.
For tetrahedron, when adding face operation, if this face overlaps with certain face tetrahedral, then disregard; If this face does not all overlap with tetrahedral any surface, then the operation adding face is converted into the operation that joining is asked on this face and this tetrahedral six limits, and concrete operations are as follows:
1) if this limit does not exist intersection point with the face of adding, next limit of traversal is continued;
2) if there is intersection point, then operation front carrying out adding point to this intersection point is being worked as;
3) tetrahedral six limits have all traveled through, and terminate.
In the step S104 of the present embodiment, after generating new tetrahedral grid unit, find the hexahedral mesh unit that these new subdivision grid cells are corresponding, upgrade numbering and the topological relation of the subdivision tetrahedral grid unit that hexahedral mesh unit inside comprises, generate composite gridding system.
In the step S105 of the present embodiment, for the method flow adding well location data shown in Fig. 2, illustrate and the well location data in described model of geological structure body are joined in described composite gridding system, carry out the step of the structure of layer auto-tracking and geologic body grid, it comprises: step S1051, searches the initial hexahedral mesh unit at the coordinate place of described well location data; Step S1052, all tetrahedral grid unit that the initial hexahedral mesh unit traveling through the coordinate place of described well location data comprises, search the tetrahedral grid unit comprising the coordinate of described well location data; Step S1053, adds described well location in the tetrahedral grid unit of coordinate comprising described well location data; Step S1054, upgrades the topological relation of the tetrahedral grid system after adding well location data.
Utilize the method for the structure of the three-dimensional geologic based on the composite gridding modeling of above-described embodiment, computation complexity when determining the tetrahedral grid unit at coordinate points place of well location data is O(1), if and the computation complexity adopting general unstrctured grid method (Delaunay method) is O(M), the total number of grid of M for comprising in unstrctured grid herein.Therefore, the composite gridding system in the present embodiment is faster in speed, will be effectively than general unstrctured grid method for the well location data of new interpolation and the process of profile data.
Fig. 3 is the structural representation of the structure of the three-dimensional geologic based on the composite gridding modeling of the embodiment of the present invention.As shown in the figure, three-dimensional geologic based on the composite gridding structure modeling of this embodiment comprises: initial hexahedral mesh system generation unit 101, for the grid system region of selected structure modeling, and mesh spacing in described grid system region on the initial point of selected grid system and coordinate direction and number, generate initial hexahedral mesh system, and record the numbering of each initial hexahedral mesh unit in described initial hexahedral mesh system; Tetrahedral grid system generation unit 102, for being multiple tetrahedral grid unit by described each initial hexahedral mesh dividing elements, generates tetrahedral grid system; Topology generation unit 103, for being numbered described multiple tetrahedral grid unit, recording the numbering of the tetrahedral grid unit comprised in described each hexahedral mesh unit, generating the topological relation of described tetrahedral grid system; Composite gridding system generation unit 104, adds the actual boundary of the model of geological structure body of setting, carries out subdivision again to described tetrahedral grid system, generates composite gridding system; Three-dimensional geologic grid enable unit 105, for the well location data in described three-dimensional geological body Model and profile data being joined in described composite gridding system, carries out the structure of layer auto-tracking and three-dimensional geological volume mesh.
In the present embodiment, the grid system region of the structure modeling that initial hexahedral mesh system generation unit 101 is selected is hexahedron/rectangular parallelepiped region.Initial mesh system in the embodiment of the present invention is SST hexahedral meshes, if determine grid system initial point, along the mesh spacing of all directions and meshes number.The generation of SST hexahedral meshes belongs to the common practise of this area, therefore repeats no more in the embodiment of the present invention.
In the present embodiment, each initial hexahedral mesh unit in described initial hexahedral mesh system is divided into multiple tetrahedral grid unit according to unified order by tetrahedral grid system generation unit 102, generates tetrahedral grid system.In the present embodiment, initial mesh unit is hexahedron, then can be divided into five tetrahedrons.
In the present embodiment, the topological relation of the described tetrahedral grid system of described Topology generation unit 103 generation comprises: the shared relationship of joint points, edges, faces, body.Initial mesh unit is hexahedron, can be divided into five tetrahedrons, then its topological relation is the shared relationship of joint points, edges, faces and body.
In the present embodiment, as shown in Figure 4, described composite gridding system generation unit 104 comprises: adding point module 1041, for adding a little in multiple tetrahedral grid unit of described tetrahedral grid system, generating new tetrahedral grid unit; Adding limit module 1042, for adding limit in multiple subdivision tetrahedral grid unit of described tetrahedral grid system, generating new tetrahedral grid unit; Adding face mould block 1043, for adding face in multiple tetrahedral grid unit of described tetrahedral grid system, generating new tetrahedral grid unit.In the present embodiment, as shown in Figure 4, described composite gridding system generation unit 104 also comprises: update module 1044, for after generating new tetrahedral grid unit, find the initial hexahedral mesh unit that new tetrahedral grid unit is corresponding, the numbering of the tetrahedral grid unit that the initial hexahedral mesh unit inside upgrading described correspondence comprises and topological relation, generate composite gridding system.
In the present embodiment, well location data in described model of geological structure body and profile data join in described composite gridding system by described geologic body grid enable unit 105, carry out the structure mould of layer auto-tracking and geologic body grid, as shown in Figure 5, it comprises: initial mesh unit searches unit 1051, for searching the initial hexahedral mesh unit at the coordinate place of described well location data; Tetrahedron element grid-search unit 1052, all tetrahedral grid unit that the initial hexahedral mesh unit for traveling through the coordinate place of described well location data comprises, search the tetrahedral grid unit comprising the coordinate of described well location data; Well location data adding device 1053, for adding described well location in described comprising in the tetrahedral grid unit of the coordinate of described well location data; Three-dimensional geologic mess generation unit 1054, for upgrading the topological relation of the tetrahedral grid system after adding well location data, generates three-dimensional geological volume mesh.
Utilize the system of the structure of the three-dimensional geologic based on the composite gridding modeling of above-described embodiment, computation complexity when determining the tetrahedral grid unit at coordinate points place of well location data is O(1), if and the computation complexity adopting general unstrctured grid method (Delaunay method) is O(M), the total number of grid of M for comprising in unstrctured grid herein.Therefore, the composite gridding system in the present embodiment is faster in speed, will be effectively than general unstrctured grid method for the well location data of new interpolation and the process of profile data.
Specific embodiment one:
As shown in Figure 6, be modeled as specific embodiment with three-dimensional geologic structure and the method flow carrying out three-dimensional geologic structure based on composite gridding of the present invention is described.
Step S601, the given hexahedron region needing structure modeling.
Step S602, respective mesh spacing and meshes number on the initial point of selected good hexahedral mesh system and x, y, z direction, generate the initial hexahedral mesh system in enough large this region of covering, as shown in Figure 7, and record respective mesh spacing, meshes number on initial point and x, y, z direction.
Step S603, is divided into five tetrahedrons by each hexahedral mesh unit of the initial hexahedral mesh system in step S602 by unified order, as shown in Figure 8, and records the numbering of the tetrahedron element that each hexahedral mesh unit inside comprises.Fig. 9 is subdivision five tetrahedrons out.
Step S604, through step S602, hexahedral mesh system just becomes tetrahedral grid system, as shown in Figure 10.To its topological relation of tetrahedral grid system record generated, comprise the shared relationship of joint points, edges, faces, body.
Step S605, on the basis of step S604, adds the border of given three-dimensional geological body region, often adds a segment boundary and just carries out subdivision to regular tetrahedral grid system.Concrete subdivision method relates separately to successively in single tetrahedron and adds a little, adds limit, adds the operation in face:
1) point is added: with this point for starting point
11) if this point is inner at tetrahedron, connect three summits of this point with tetrahedral four faces successively, form four new tetrahedrons, as shown in figure 11.
12) if to drop on tetrahedron face inner for this point, then three summits first connecting this summit and this face form three line segments, then connect this summit and this face to point, form three new tetrahedrons, as shown in figure 12.
13) whether if this point drops on the inside on a tetrahedron limit, then two summits first connecting this summit and this limit form two line segments, then connect this summit with or not coplanar to point this summit, form two new tetrahedrons, as shown in figure 13.
14) if this point drops on certain summit of tetrahedron, then disregard.
2) limit is added:
21) if this limit overlaps with tetrahedral certain limit, disregard.
22) if this limit does not all overlap with tetrahedral arbitrary limit:
221) if the summit on this limit is certain summit tetrahedral, the operation then adding limit just can transform into a particular ray (starting point by the summit on interpolation limit, direction is along the terminal on added limit) with the operation of this tetrahedron find intersection, the operation adding point is carried out to tried to achieve intersection point, then complete the operation this tetrahedron being added to limit, as shown in figure 14.
222) if the summit on this limit is not certain summit tetrahedral, for sake of convenience, this tetrahedron is represented with Tetra.If the summit on this limit is not certain summit of Tetra, then first ask particular ray (starting point by the summit on interpolation limit, direction is along the terminal on added limit) first the intersection point A crossing with Tetra, then antinode A carries out the operation adding point in Tetra, and such Tetra is just split into as N number of sub-tetrahedron.Travel through this N number of sub-tetrahedron, find the sub-tetrahedron of this ray and which crossing and try to achieve intersection points B (representing this sub-tetrahedron with STetra), and then antinode B carries out the operation adding point in STetra, then Tetra subdivision completes, as shown in figure 15.
3) face is added:
31) if this face overlaps with certain face tetrahedral, then disregard.
32) if this face does not all overlap with tetrahedral any surface; The operation in interpolation face is converted into the operation that joining is asked on this face and this tetrahedral six limits, specifically, from this tetrahedral Article 1 limit:
321) if this limit does not exist intersection point with the face of adding, next limit of traversal is continued;
322) if there is intersection point, then when on front to this intersection point carry out based on add point operation.
Until tetrahedral six limits have all traveled through, process has terminated, to complete the operation this tetrahedron being added to face, as shown in figure 16.
Step S606, to new subdivision tetrahedron element out, finds the hexahedral mesh unit comprising them, upgrades the numbering of the tetrahedron element that this hexahedral inner comprises, and be recorded in hexahedral mesh system.Further, upgrade the topological relation (shared relationship of joint points, edges, faces, body) of new subdivision tetrahedron element out, and this relation of inclusion is recorded in tetrahedral grid system.To tetrahedral grid system, only retain the topological structure of the mesh generation of given intra-zone.
Step S607, joins in composite gridding system with reference to point (well location) data and profile data.Value according to reference point carries out interpolation to all tetrahedral grid cell nodes in composite gridding system.
The enforcement of this step can embody the advantage of composite gridding system described in the invention well, first its embodiment is described in detail (here only to add reference point data instance, add data class on section seemingly, only to need section to be considered as a string discrete point);
Hypothetical reference point data is x1, x2 ..., xn, data total amount is n;
1) from first data x1, the hexahedral mesh unit at the coordinate place of x1 is found;
2) all tetrahedral grid unit that this hexahedral mesh unit comprises are determined;
3) travel through 2 successively) in all tetrahedral grid unit, until find the tetrahedral grid unit comprising x1 coordinate points;
4) to 3) in the tetrahedral grid unit determined carry out the operation adding point, the operation of this interpolation point is identical with the operation of the interpolation point in step S104;
5) by 4) in new subdivision tetrahedron element out add 1 to) in the hexahedral element determined;
6) topological structure of tetrahedral grid system is upgraded.
To x2 ..., xn repeats 1)--6) step, then reference point data have been added.
As can be seen from above specific embodiment, computation complexity when determining the subdivision tetrahedral grid unit at coordinate points place of well location data is O(1), if and the computation complexity adopting general unstrctured grid method (Delaunay method) is O(M), the total number of grid of M for comprising in unstrctured grid herein.Therefore, the composite gridding system in the present embodiment is faster in speed, will be effectively than general unstrctured grid method for the well location data of new interpolation and the process of profile data.
Three-dimensional geologic based on the composite gridding structure modeling method of the embodiment of the present invention and system, can generate a kind of composite gridding system with complementary two nested grid systems, can simulate landform complicated arbitrarily.Further, input needs the region of modeling, can generate the three-dimensional composite gridding system in this region fast.Given arbitrary coordinate position (being positioned at region), can the numbering of tetrahedral grid unit at this coordinate place of quick position.This composite gridding system can add the new node (coordinate is positioned at given area) of any amount very easily, only need carry out Local grid subdivision (at new coordinate points place) and topological relation renewal to original grid system.Therefore, the present invention can also realize grid fast query and the location of given coordinate position.
Those skilled in the art can also recognize the various illustrative components, blocks (illustrative logical block) that the embodiment of the present invention is listed, unit, and step can pass through electronic hardware, computer software, or both combinations realize.For the replaceability (interchangeability) of clear displaying hardware and software, above-mentioned various illustrative components (illustrative components), unit and step have universally described their function.Such function is the designing requirement realizing depending on specific application and whole system by hardware or software.Those skilled in the art for often kind of specifically application, can use the function described in the realization of various method, but this realization can should not be understood to the scope exceeding embodiment of the present invention protection.
The software module that method described in the embodiment of the present invention or the step of algorithm directly can embed hardware, processor performs or the combination of both.Software module can be stored in the storage medium of other arbitrary form in RAM storer, flash memory, ROM storer, eprom memory, eeprom memory, register, hard disk, moveable magnetic disc, CD-ROM or this area.Exemplarily, storage medium can be connected with processor, with make processor can from storage medium reading information, and write information can be deposited to storage medium.Alternatively, storage medium can also be integrated in processor.Processor and storage medium can be arranged in ASIC, and ASIC can be arranged in user terminal.Alternatively, processor and storage medium also can be arranged in the different parts in user terminal.
In one or more exemplary design, the above-mentioned functions described by the embodiment of the present invention can realize in the combination in any of hardware, software, firmware or this three.If realized in software, these functions can store on the medium with computer-readable, or are transmitted on the medium of computer-readable with one or more instruction or code form.Computer readable medium comprises computer storage medium and is convenient to make to allow computer program transfer to the telecommunication media in other place from a place.Storage medium can be that any general or special computer can the useable medium of access.Such as, such computer readable media can include but not limited to RAM, ROM, EEPROM, CD-ROM or other optical disc storage, disk storage or other magnetic storage device, or other anyly may be used for carrying or store the medium that can be read the program code of form with instruction or data structure and other by general or special computer or general or special processor.In addition, any connection can be properly termed computer readable medium, such as, if software is by a concentric cable, fiber optic cables, twisted-pair feeder, Digital Subscriber Line (DSL) or being also comprised in defined computer readable medium with wireless way for transmittings such as such as infrared, wireless and microwaves from a web-site, server or other remote resource.Described video disc (disk) and disk (disc) comprise Zip disk, radium-shine dish, CD, DVD, floppy disk and Blu-ray Disc, and disk is usually with magnetic duplication data, and video disc carries out optical reproduction data with laser usually.Above-mentioned combination also can be included in computer readable medium.
Those skilled in the art should understand, embodiments of the invention can be provided as method, system or computer program.Therefore, the present invention can adopt the form of complete hardware embodiment, completely software implementation or the embodiment in conjunction with software and hardware aspect.And the present invention can adopt in one or more form wherein including the upper computer program implemented of computer-usable storage medium (including but not limited to magnetic disk memory, CD-ROM, optical memory etc.) of computer usable program code.
The present invention describes with reference to according to the process flow diagram of the method for the embodiment of the present invention, equipment (system) and computer program and/or block scheme.Should understand can by the combination of the flow process in each flow process in computer program instructions realization flow figure and/or block scheme and/or square frame and process flow diagram and/or block scheme and/or square frame.These computer program instructions can being provided to the processor of multi-purpose computer, special purpose computer, Embedded Processor or other programmable data processing device to produce a machine, making the instruction performed by the processor of computing machine or other programmable data processing device produce device for realizing the function of specifying in process flow diagram flow process or multiple flow process and/or block scheme square frame or multiple square frame.
These computer program instructions also can be stored in can in the computer-readable memory that works in a specific way of vectoring computer or other programmable data processing device, the instruction making to be stored in this computer-readable memory produces the manufacture comprising command device, and this command device realizes the function of specifying in process flow diagram flow process or multiple flow process and/or block scheme square frame or multiple square frame.These computer program instructions also can be loaded in computing machine or other programmable data processing device, make on computing machine or other programmable devices, to perform sequence of operations step to produce computer implemented process, thus the instruction performed on computing machine or other programmable devices is provided for the step realizing the function of specifying in process flow diagram flow process or multiple flow process and/or block scheme square frame or multiple square frame.
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (8)

1., based on a three-dimensional geologic structure modeling method for composite gridding, it is characterized in that, described method comprises:
Steps A, the grid system region of selected structure modeling, and mesh spacing in described grid system region on the initial point of selected grid system and coordinate direction and number, generate initial hexahedral mesh system, and record the numbering of each initial hexahedral mesh unit in described initial hexahedral mesh system;
Described each initial hexahedral mesh dividing elements is multiple tetrahedral grid unit by step B, generates tetrahedral grid system;
Step C, is numbered described multiple tetrahedral grid unit, records the numbering of the tetrahedral grid unit comprised in described each hexahedral mesh unit, generates the topological relation of described tetrahedral grid system;
Step D, adds the actual boundary of the model of geological structure body of setting, carries out subdivision again to described tetrahedral grid system, generates composite gridding system;
Well location data in described three-dimensional geological body Model and profile data are joined in described composite gridding system, carry out the structure of layer auto-tracking and three-dimensional geological volume mesh by step e;
In described step e, the well location data in described model of geological structure body and profile data are joined in described composite gridding system, carry out the structure of layer auto-tracking and geologic body grid, comprising:
Search the initial hexahedral mesh unit at the coordinate place of described well location data;
All tetrahedral grid unit that the initial hexahedral mesh unit traveling through the coordinate place of described well location data comprises, search the tetrahedral grid unit comprising the coordinate of described well location data;
Described well location data are added in described comprising in the tetrahedral grid unit of described well location data coordinates;
Upgrade the topological relation of the tetrahedral grid system after adding well location data.
2. the structure of the three-dimensional geologic based on composite gridding modeling method according to claim 1, it is characterized in that, in described step C, the topological relation of the described tetrahedral grid system of generation comprises: the shared relationship of joint points, edges, faces, body.
3. the structure of the three-dimensional geologic based on composite gridding modeling method according to claim 1, is characterized in that, in described step D, add the actual boundary of the model of geological structure body of setting, carry out subdivision again, comprising described tetrahedral grid system:
In multiple tetrahedral grid unit of described tetrahedral grid system, add point, limit or face, generate new subdivision tetrahedral grid unit.
4. the structure of the three-dimensional geologic based on composite gridding modeling method according to claim 3, it is characterized in that, in described step D, after generating new subdivision tetrahedral grid unit, find the initial hexahedral mesh unit that new subdivision tetrahedral grid unit is corresponding, the numbering of the tetrahedral grid unit that the initial hexahedral mesh unit inside upgrading described correspondence comprises and topological relation, generate composite gridding system.
5., based on a three-dimensional geologic structure modeling for composite gridding, it is characterized in that, described system comprises:
Initial hexahedral mesh system generation unit, for the grid system region of selected structure modeling, and mesh spacing in described grid system region on the initial point of selected grid system and coordinate direction and number, generate initial hexahedral mesh system, and record the numbering of each initial hexahedral mesh unit in described initial hexahedral mesh system;
Tetrahedral grid system generation unit, for being multiple tetrahedral grid unit by described each initial hexahedral mesh dividing elements, generates tetrahedral grid system;
Topology generation unit, for being numbered described multiple tetrahedral grid unit, recording the numbering of the tetrahedral grid unit comprised in described each hexahedral mesh unit, generating the topological relation of described tetrahedral grid system;
Composite gridding system generation unit, adds the actual boundary of the model of geological structure body of setting, carries out subdivision again to described tetrahedral grid system, generates composite gridding system;
Three-dimensional geologic grid enable unit, for the well location data in described three-dimensional geological body Model and profile data being joined in described composite gridding system, carries out the structure of layer auto-tracking and three-dimensional geological volume mesh;
Described three-dimensional geologic grid enable unit comprises:
Initial hexahedral mesh unit searches unit, for searching the initial hexahedral mesh unit at the coordinate place of described well location data;
Tetrahedron element grid-search unit, all tetrahedral grid unit that the initial hexahedral mesh unit for traveling through the coordinate place of described well location data comprises, search the tetrahedral grid unit comprising the coordinate of described well location data;
Well location data adding device, for adding described well location data in described comprising in the tetrahedral grid unit of described well location data coordinates;
Three-dimensional geologic mess generation unit, for upgrading the topological relation of the tetrahedral grid system after adding well location data, generates three-dimensional geological volume mesh.
6. the structure of the three-dimensional geologic based on composite gridding modeling according to claim 5, it is characterized in that, the topological relation of the tetrahedral grid system that described Topology generation unit generates comprises: the shared relationship of joint points, edges, faces, body.
7. the structure of the three-dimensional geologic based on composite gridding modeling according to claim 5, it is characterized in that, described composite gridding system generation unit comprises:
Adding point module, for adding a little in multiple tetrahedral grid unit of described tetrahedral grid system, generating new tetrahedral grid unit;
Adding limit module, for adding limit in multiple tetrahedral grid unit of described tetrahedral grid system, generating new tetrahedral grid unit;
Adding face mould block, for adding face in multiple tetrahedral grid unit of described tetrahedral grid system, generating new tetrahedral grid unit.
8. the structure of the three-dimensional geologic based on composite gridding modeling according to claim 7, it is characterized in that, described composite gridding system generation unit also comprises:
Update module, for after generating new tetrahedral grid unit, find the initial hexahedral mesh unit that new tetrahedral grid unit is corresponding, the numbering of the tetrahedral grid unit that the initial hexahedral mesh unit inside upgrading described correspondence comprises and topological relation, generate composite gridding system.
CN201310130973.2A 2013-04-16 2013-04-16 Three-dimensional geological structure modeling method and system based on composite grid Active CN103236088B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310130973.2A CN103236088B (en) 2013-04-16 2013-04-16 Three-dimensional geological structure modeling method and system based on composite grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310130973.2A CN103236088B (en) 2013-04-16 2013-04-16 Three-dimensional geological structure modeling method and system based on composite grid

Publications (2)

Publication Number Publication Date
CN103236088A CN103236088A (en) 2013-08-07
CN103236088B true CN103236088B (en) 2015-07-15

Family

ID=48884127

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310130973.2A Active CN103236088B (en) 2013-04-16 2013-04-16 Three-dimensional geological structure modeling method and system based on composite grid

Country Status (1)

Country Link
CN (1) CN103236088B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103824331A (en) * 2014-01-28 2014-05-28 浙江大学 Iteration generation method of hexahedral mesh model
CN104504756A (en) * 2015-01-20 2015-04-08 西安煤航卫星数据应用有限公司 Dynamic updating method of ore deposit three dimension geology model
CN105551085B (en) * 2015-12-30 2019-02-19 中国科学院东北地理与农业生态研究所 A kind of construction method of doubly ruled surface geologic body
CN109584369B (en) * 2018-10-24 2021-09-28 中国石油天然气股份有限公司 Actual stratum full hexahedron grid generation method and device
CN109658505B (en) * 2018-11-05 2023-05-09 武汉中地数码科技有限公司 Three-dimensional space conformation method and system with fault plane
CN111815769A (en) * 2019-04-12 2020-10-23 中国石油天然气股份有限公司 Modeling method, computing device and storage medium for thrust-driven tectonic belt structure
CN111414657B (en) * 2020-03-09 2023-02-10 中煤科工集团沈阳设计研究院有限公司 Method for determining coal seam bifurcation position and establishing refined three-dimensional solid model
CN111414654B (en) * 2020-03-16 2024-02-23 南京交通职业技术学院 Method for retrieving continuous numerical model of division
CN111985124B (en) * 2020-07-15 2023-07-07 北京金阳普泰石油技术股份有限公司 Single sand body model construction method, system, equipment and readable storage medium
CN114491782B (en) * 2022-04-18 2022-07-05 煤炭科学研究总院有限公司 Underground engineering surrounding rock modeling method based on combined mesh subdivision

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8150663B2 (en) * 2007-03-30 2012-04-03 Paradigm Geophysical (Luxembourg) S.A.R.L. Partitioning algorithm for building a stratigraphic grid
CN102194253B (en) * 2011-06-15 2013-09-18 北京航空航天大学 Method for generating tetrahedron gridding for three-dimensional geological structure

Also Published As

Publication number Publication date
CN103236088A (en) 2013-08-07

Similar Documents

Publication Publication Date Title
CN103236088B (en) Three-dimensional geological structure modeling method and system based on composite grid
CN102306396B (en) Three-dimensional entity model surface finite element mesh automatic generation method
Ma et al. Associative feature modeling for concurrent engineering integration
Crocker et al. An editable nonmanifold boundary representation
KR102040988B1 (en) Creating a surface from a plurality of 3d curves
WO2011163166A2 (en) Machine, program product, and computer-implemented method to simulate reservoirs as 2.5d unstructured grids
CN105761303A (en) Creation Of Bounding Boxes On 3d Modeled Assembly
CN102682476B (en) The boolean operation method of triangle grid data and system thereof
CN105118091B (en) A kind of method and system constructing the non-homogeneous geology mesh surface model of more precision
CN102819865A (en) Modeling method for magnetotelluric three-dimensional geologic structure model
CN105022865A (en) STL model Boolean calculation based method for extracting inner surface model of fuel tank in airplane
CN102194252A (en) Geological-stratum-structure-based method for generating triangular lattice grids
CN102737407A (en) Fitting optimization method of triangular mesh data and system for achieving fitting optimization method
CN105243237B (en) A kind of construction method of assembly technology structure tree
CN106649817A (en) Three-dimensional pipe model establishment lofting method and device for geography information system
CN105225272A (en) A kind of tri-dimensional entity modelling method based on the reconstruct of many outline lines triangulation network
Penninga et al. A tetrahedronized irregular network based DBMS approach for 3D topographic data modeling
CN104965938A (en) Designing industrial products by using geometries connected by geometrical constraints
CN103279602B (en) The delicate tissues method of the geographical entity space geometric data of a kind of rule
CN104700453A (en) Method for generating three-dimensional model of cable pipe network based on GIS
CA2803144A1 (en) Machine, program product, and computer-implemented method to simulate reservoirs as 2.5d unstructured grids
Li et al. Immersion of self-intersecting solids and surfaces
Aubry et al. A surface remeshing approach
CN104090945B (en) Geographic space entity constructing method and system
Ohori et al. Modeling and manipulating spacetime objects in a true 4D model

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant