CN109584369B - Actual stratum full hexahedron grid generation method and device - Google Patents
Actual stratum full hexahedron grid generation method and device Download PDFInfo
- Publication number
- CN109584369B CN109584369B CN201811241872.1A CN201811241872A CN109584369B CN 109584369 B CN109584369 B CN 109584369B CN 201811241872 A CN201811241872 A CN 201811241872A CN 109584369 B CN109584369 B CN 109584369B
- Authority
- CN
- China
- Prior art keywords
- layer
- layers
- straight line
- grid node
- actual stratum
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 92
- 238000012952 Resampling Methods 0.000 claims abstract description 30
- 238000004590 computer program Methods 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 5
- 238000002360 preparation method Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 description 23
- 238000010586 diagram Methods 0.000 description 18
- 239000011159 matrix material Substances 0.000 description 10
- 230000008901 benefit Effects 0.000 description 8
- 238000012545 processing Methods 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 6
- 238000013507 mapping Methods 0.000 description 6
- 239000007787 solid Substances 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000009941 weaving Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005686 electrostatic field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses a method and a device for generating a full hexahedral mesh of an actual stratum, wherein the method comprises the following steps: regular resampling is carried out on the actual stratum layer position, and grid node coordinates of the regular actual stratum layer position are obtained; determining the straight line of each grid node coordinate in each layer in the regularized actual stratum horizon and the corresponding grid node coordinate in the adjacent layer; determining the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction, and subdividing the space between each layer and adjacent layers in the regularized actual stratum layer according to the straight line and the number of subdivision layers to obtain the grid node coordinates of each subdivision layer between each layer and adjacent layers on the straight line; and generating the hexahedral mesh according to the mesh node coordinates of each layer, the adjacent layers and each split layer between the two layers. The method and the device can quickly generate the full hexahedral mesh in the actual stratum subdivision, and have high automation degree and good mesh subdivision quality.
Description
Technical Field
The invention relates to the field of reservoir modeling research in oil exploration and development, in particular to a method and a device for generating a full hexahedral mesh of an underground complex reservoir.
Background
Under the continuous efforts of numerous geophysicists, the finite element method is not only widely applied to the calculation of static fields in the geophysics, such as static magnetic fields, electrostatic fields, ground stress, geothermal fields and the like, but also applied to solving the geophysical problems changing along with time, such as seismic forward modeling, seismic inversion, numerical reservoir modeling and the like. The first step in carrying out a finite element calculation is to discretize a continuous medium in the earth physics, such as an actual earth formation, into a number of finite subdomains (elements), i.e., finite element mesh generation, and approximate the continuous medium by an assembly of these discrete elements. From the geometric shape of the generating unit, aiming at a two-dimensional continuous medium, the currently generated mesh mainly adopts triangular and quadrilateral units; for three-dimensional continuous media, the currently generated meshes mainly adopt tetrahedral and hexahedral units. Furthermore, for three-dimensional continuous media there are also some other geometrical elements like pyramidal elements, prismatic elements etc., but these are not commonly used. In the process of research on the finite element mesh generation technology, research on the tetrahedral (or triangular) mesh generation technology is mature, and the tetrahedral (or triangular) mesh generation technology is widely applied to commercial software. Compared with the tetrahedral unit, the hexahedral unit has the advantage of higher precision under the condition of the same degree of freedom, and is a hotspot and difficulty of the current finite element mesh generation technology research.
The method for generating the full hexahedron grid is commonly used by a mapping method, a grid-based method, a tetrahedron conversion method, a medial axial plane decomposition method, a paving and paving method, a whisker segment weaving method and a scanning method. The mapping method firstly adopts a manual or automatic method to decompose the region to be subdivided into hexahedral subregions (for three-dimensional regions), and then an isoparametric mapping or infinite mapping method is applied to each subregion to obtain a more detailed hexahedral mesh. The mapping method has the advantages of high quality and high speed of generating units, but the problems of automatic decomposition of a complex region, grid density transition, grid compatibility of different sub-regions and the like need to be solved. The grid-based method has the advantages of high subdivision automation degree and easiness in control of grid density, but the quality of grid units is poor when complex boundary conditions are processed, and the program implementation difficulty is high. The tetrahedral mesh generation algorithm is mature, and the hexahedral mesh generation result can be obtained by decomposing one tetrahedron into four hexahedrons, so the tetrahedral conversion algorithm is simple to implement and high in universality, but the number of meshes is large, a large number of irregular nodes are obtained in the conversion process, and the quality of the meshes is poor. The mid-axial plane decomposition method is suitable for subdivision of certain specific regular solid models and is not suitable for automatic subdivision of geometric regions in any shapes. The paving layering method gradually advances to divide grids from the boundary of the area to the inner layer, the quality of the generated grid boundary unit is good, but the quality of the grid in the area is poor. The hexahedral mesh generated by the whisker segment weaving method is high in quality, but the stability and the reliability are not verified. The scanning method is suitable for subdivision of a three-dimensional solid model with a simple and regular geometric shape, requires more man-machine interaction and has low automation degree.
The actual stratum is obtained by regular interpolation in the transverse direction according to the horizon explained by the seismic data and stacking in sequence from small to large according to time or depth values in the longitudinal direction, and has the characteristics of large transverse fluctuation change, strong irregularity and easy occurrence of voids at the boundary. For the subdivision object, the hexahedron subarea is difficult to divide by the mapping subdivision method; the top and bottom boundaries are processed inaccurately based on a grid method; the tetrahedron conversion method needs to firstly divide a tetrahedron and then ensures unit quality in the process of converting the tetrahedron into a hexahedron; the medial axial plane decomposition method, the paving and paving method, the whisker segment weaving method and the scanning method are not suitable for processing complex fluctuated outer boundaries and any inner boundaries, and have low automation degree.
Disclosure of Invention
The embodiment of the invention provides a method for generating full hexahedron grids of an actual stratum, which is used for quickly generating full hexahedron grids in actual stratum subdivision and has the advantages of high automation degree and good grid subdivision quality. The method comprises the following steps:
regular resampling is carried out on the actual stratum layer position, and grid node coordinates of the regular actual stratum layer position are obtained;
determining the straight line of each grid node coordinate in each layer in the regularized actual stratum horizon and the corresponding grid node coordinate in the adjacent layer;
determining the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction, and subdividing the space between each layer and adjacent layers in the regularized actual stratum layer according to the straight line and the number of subdivision layers to obtain the grid node coordinates of each subdivision layer between each layer and adjacent layers on the straight line;
generating hexahedral meshes according to the mesh node coordinates of each layer, the adjacent layers and each split layer between the two layers;
regular resampling is carried out on the actual stratum layer position to obtain the regular actual stratum layer position grid node coordinate, and the method comprises the following steps: obtaining the maximum value and the minimum value of the main survey line direction and the contact survey line direction of each layer in the actual stratum horizon; obtaining a rectangular area of each layer according to the maximum value and the minimum value of the main survey line direction and the contact survey line direction, dividing the rectangular area of each layer into a plurality of rectangular sub-areas at equal intervals, and determining the rectangular area of each layer as each layer after regularization; and resampling each layer after the regularization to obtain grid node coordinates of the regularized actual stratum layer.
The embodiment of the invention provides a device for generating a full hexahedral mesh of an actual stratum, which comprises:
the regular resampling module is used for carrying out regular resampling on the actual stratum layer position to obtain a grid node coordinate of the regular actual stratum layer position;
the grid node generation preparation module is used for determining the coordinates of each grid node in each layer in the regularized actual stratum horizon and the straight line where the coordinates of the corresponding grid node in the adjacent layer are located;
the grid node generation module is used for determining the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction, subdividing the space between each layer and adjacent layers in the regularized actual stratum layer according to the straight line and the number of subdivision layers, and obtaining the grid node coordinates of each subdivision layer between each layer and adjacent layers on the straight line;
the hexahedral mesh generation module is used for generating hexahedral meshes according to the mesh node coordinates of each layer, the adjacent layers and each split layer between the two layers;
the rule resampling module is specifically configured to: obtaining the maximum value and the minimum value of the main survey line direction and the contact survey line direction of each layer in the actual stratum horizon; obtaining a rectangular area of each layer according to the maximum value and the minimum value of the main survey line direction and the contact survey line direction, dividing the rectangular area of each layer into a plurality of rectangular sub-areas at equal intervals, and determining the rectangular area of each layer as each layer after regularization; and resampling each layer after the regularization to obtain grid node coordinates of the regularized actual stratum layer.
In the embodiment of the invention, the actual stratum layer is regularly resampled to obtain the grid node coordinates of the regular actual stratum layer, so that the grid node coordinates of the marginal blank area are supplemented; determining a straight line where each grid node in each layer and the corresponding grid node in the adjacent layer in the regularized actual stratum horizon are located, and the number of layers which are vertically split, and obtaining grid node coordinates between each layer and the adjacent layer; the hexahedral mesh is generated according to the mesh node coordinates between each layer and the adjacent layers, the whole subdivision process does not need to divide hexahedral subregions, and the subdivision process adopting the straight line in the solid analytic geometry is simple, high in speed and high in automation degree; in addition, the generated hexahedron grid has high quality because the subdivision process does not need intermediate conversion.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
FIG. 1 is a flow chart of a method for generating a full hexahedron mesh of an actual formation according to an embodiment of the present invention;
FIG. 2 is a diagram of a real-world layer bitmap before regularization according to an embodiment of the present invention;
FIG. 3 is a normalized actual floor level map according to an embodiment of the present invention;
fig. 4 is a regularized three actual stratum layer bitmaps to be divided according to the embodiment of the present invention;
FIG. 5 is a schematic diagram illustrating a subdivision principle, mesh node numbering and hexahedron unit numbering process according to an embodiment of the present invention;
FIG. 6 is a diagram of the result of generating a full hexahedral mesh of three actual stratigraphic horizons according to an embodiment of the present invention;
FIG. 7 is a diagram illustrating an evaluation of formation horizon full hexahedron mesh generation quality in accordance with an embodiment of the present invention;
fig. 8 is a schematic structural diagram of an actual stratum full hexahedron mesh generation device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
Fig. 1 is a flowchart of an actual stratum full hexahedron mesh generation method according to an embodiment of the present invention. As shown in fig. 1, the method for generating an actual stratum full hexahedron mesh of this embodiment includes:
103, determining the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction, and subdividing the space between each layer and adjacent layers in the regularized actual stratum layer according to the straight line and the number of subdivision layers to obtain the grid node coordinates of each subdivision layer between each layer and adjacent layers on the straight line;
and 104, generating a hexahedral mesh according to the mesh node coordinates of each layer, the adjacent layers and each split layer between the two layers.
As can be seen from the process shown in fig. 1, the embodiment of the present invention performs regular resampling on the actual stratum layer to obtain the grid node coordinates of the regular actual stratum layer, so as to supplement the grid node coordinates of the marginal blank area; determining a straight line where each grid node in each layer and the corresponding grid node in the adjacent layer in the regularized actual stratum horizon are located, and the number of layers which are vertically split, and obtaining grid node coordinates between each layer and the adjacent layer; the hexahedral mesh is generated according to the mesh node coordinates between each layer and the adjacent layers, the whole subdivision process does not need to divide hexahedral subregions, and the subdivision process adopting the straight line in the solid analytic geometry is simple, high in speed and high in automation degree; in addition, the generated hexahedron grid has high quality because the subdivision process does not need intermediate conversion.
Because the actual stratum horizon is manually explained on the basis of the seismic data processing result, the phenomenon of discontinuous explanation points can occur in the explanation process, and a blank area can occur at the boundary, the actual stratum horizon is firstly regularly re-sampled before the grid is generated.
In specific implementation, the specific process of performing regular resampling on the actual stratum layer to obtain the grid node coordinates of the regularized actual stratum layer may include:
obtaining the maximum value and the minimum value of the main survey line direction (inline) and the cross survey line direction (crossline) of each layer in the actual stratum horizon;
obtaining a rectangular area of each layer according to the maximum value and the minimum value of the main survey line direction and the contact survey line direction, dividing the rectangular area of each layer into a plurality of rectangular sub-areas at equal intervals, and determining the rectangular area of each layer as each layer after regularization;
and resampling each layer after the regularization to obtain grid node coordinates of the regularized actual stratum layer.
In specific implementation, there may be many ways to resample each regularized layer, for example, a nearest neighbor interpolation method may be adopted, but it should be understood that the above-mentioned method of using nearest neighbor interpolation is merely an example, and those skilled in the art may adopt other methods to resample each regularized layer according to actual needs, and related variations should fall within the scope of the present invention.
In an embodiment, determining the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction, and subdividing the space between each layer and adjacent layers in the regularized actual stratum layer according to the straight line and the number of subdivision layers to obtain the grid node coordinates of each subdivision layer between each layer and adjacent layers on the straight line may include:
determining the number of subdivision layers of each layer and adjacent layers in the vertical direction in the regularized actual stratum horizon;
obtaining the thickness of each layer and adjacent layers on the straight line according to the number of the subdivision layers and the straight line;
calculating the coordinate of each section layer between each layer and the adjacent layer in the straight line vertical direction according to the thickness;
and calculating to obtain the grid nodes of each split layer between each layer and the adjacent layer on the straight line according to the coordinates of each split layer on the vertical direction of the straight line and the straight line.
In an embodiment, generating the hexahedral mesh according to the mesh node coordinates of each layer and adjacent layers and each hierarchical layer therebetween may include:
numbering grid node coordinates of each layer, adjacent layers and each split layer between the adjacent layers;
forming a hexahedron unit by the eight adjacent grid nodes to obtain a hexahedron unit set;
and numbering the hexahedron units in the hexahedron unit set to generate hexahedron grids, and outputting grid node coordinates, grid node numbers, hexahedron unit numbers and grid node numbers forming the hexahedron units.
In specific implementation, there are many ways to number the mesh node coordinates of the actual stratigraphic layer position that is subdivided, for example, a numbering method of "from top to bottom, from left to right, from front to back" may be adopted to obtain all mesh node numbers, which may specifically include:
setting top grid nodes at the forefront and the leftmost of all grid nodes as 1; continuously numbering all grid nodes on a straight line of the top grid node along the vertical downward direction until all grid nodes on the straight line are numbered; continuously numbering the straight lines according to grid nodes on adjacent straight lines in the left-to-right direction until the numbering of the network nodes on all the straight lines in the left-to-right direction is completed; and continuously numbering the grid nodes on the straight line on the plane formed by the straight line and the adjacent straight line according to the sequence from front to back until all the grid nodes are numbered.
Of course, it should be understood that the above method of numbering "from top to bottom, left to right, and front to back" is only an example, and those skilled in the art may adopt other methods of numbering grid node coordinates according to actual needs, and all the related variations should fall within the scope of the present invention.
In one embodiment, after the hexahedral mesh is generated, the hexahedral mesh generation quality may be evaluated, and in particular, there are many methods for evaluating the hexahedral mesh generation quality.
In one embodiment, the specific process of evaluating the quality of the hexahedral mesh generation may include:
calculating the value of the unitized Jacobian matrix determinant of 8 grid nodes in each hexahedral unit in the hexahedral grid, and determining the value of the minimum Jacobian matrix determinant of each hexahedral unit, wherein the value range of the value of the minimum Jacobian matrix determinant is-1 to 1; and judging the degree that the value of the determinant of the minimum Jacobian matrix is close to 1, wherein the closer the value is to 1, the better the quality of the corresponding hexahedral unit is, and the value is less than zero, so that the corresponding hexahedral unit does not meet the quality requirement.
A specific embodiment is given below to illustrate a specific application of the method for generating a full hexahedral mesh of a practical stratum according to the present invention.
Please refer to fig. 2, fig. 3 and fig. 4. Fig. 2 is a layer bitmap of an actual horizon before regularization according to an embodiment of the present invention, where a selected horizon is a manually explained horizon of a certain western basin based on seismic processing data, and it can be seen that the horizon has a phenomenon of missing points at a boundary, and the explained points are more and are distributed irregularly.
Fig. 3 is a normalized actual layer bitmap according to this embodiment. In order to display the grid more clearly, the embodiment adopts a 100 × 100 grid to regularize the entire horizon, and the specific process includes: obtaining the maximum value and the minimum value of the main survey line direction and the contact survey line direction of each layer in the actual stratum horizon; obtaining a rectangular area of each layer according to the maximum value and the minimum value of the main survey line direction and the contact survey line direction, dividing the rectangular area of each layer into a plurality of rectangular sub-areas at equal intervals, and determining the rectangular area of each layer as each layer after regularization; and resampling each layer after the regularization to obtain grid node coordinates of the regularized actual stratum layer.
Fig. 4 is a layer bitmap of three actual layers to be divided after regularization in this embodiment, and subsequent mesh generation is implemented on the basis of the three layers.
Referring to fig. 5, fig. 5 is a schematic diagram illustrating a subdivision principle, mesh node numbering, and hexahedron cell numbering process according to an embodiment of the present invention.
In this embodiment, the level a and the level B represent eight regular matrix grids intercepted from the regularized adjacent levels, and the number of subdivision layers between the level a and the level B is determined to be 2, i.e. a new level C is inserted, where the grid nodes (x1, y1, z1) and the grid nodes (x2, y2, z2) are two corresponding points on the level a and the level B, respectively, and are two points of a straight line equation in a stereo-analytic geometry according to the following:
(x-x1)/(x2-x1)＝(y-y1)/(y2-y1)＝(z-z1)/(z2-z1)
grid node (x, y, z) is a coordinate point in the cartesian coordinate system of horizon C and is located on a straight line determined by grid node (x1, y1, z1) and grid node (x2, y2, z 2);
in this embodiment, the number of subdivision layers is 2, so the thickness zn of the horizon a and the horizon B on the straight line can be calculated by using the following expression:
zn＝|z2-z1|/2；
then, the coordinate z' of the horizon C in the vertical direction of the straight line is calculated by adopting the following expression:
z’＝z1+zn；
and substituting the coordinate z 'into two points of the linear equation, and calculating to obtain a coordinate x' and a coordinate y ', namely determining the grid node coordinate (x', y ', z') inserted by the horizon C.
The coordinate points of the horizon C can be sequentially obtained by adopting the method, and finally the coordinates of all grid nodes inserted into the horizon C are obtained.
Continuing with fig. 5, the mesh nodes of the actual stratum layer after being subdivided are numbered, in this embodiment, all mesh node numbers are obtained by using a method of "numbering from top to bottom, from left to right, from front to back", and the specific process is as follows:
the mesh node (x1, y1, z1) number is 1, the mesh node (x ', y ', z ') number is 2, and the mesh node (x2, y2, z2) number is 3. After the grid nodes (x1, y1, z1) and all the grid nodes on the straight line where the grid nodes (x ', y ', z ') are located are numbered, the straight line is numbered continuously from the adjacent straight line in the left-to-right direction, and after the grid nodes on all the straight lines of the front-most plane are numbered, the grid nodes on the straight line of the front-most adjacent plane are numbered continuously. For example, in this embodiment, eight adjacent grid nodes, namely grid nodes 1, 2, 4, 5, 10, 11, 13, and 14, constitute a unit E1, and eight adjacent grid nodes, namely grid nodes 2, 3, 5, 6, 11, 12, 14, and 15, constitute a hexahedron unit E2. Fig. 5 is a schematic diagram of a numbering principle, and the actual strata can be analogized according to the principle. And finally, outputting the grid node coordinates, the grid node numbers, the hexahedral unit numbers and the grid node numbers forming the hexahedral units.
FIG. 6 is a diagram of the result of generating a full hexahedral mesh of three actual stratigraphic horizons according to an embodiment of the present invention. The whole hexahedron grid generation process is realized by programming of a common desktop, the number of generated grid nodes is 90000, the number of hexahedron units is 101178, and the generation speed of the hexahedron grid meets the requirements of industrial application when the actual stratum horizon data and the hexahedron grid generation are read and the hexahedron grid generation result is output for 36 s.
FIG. 7 is a diagram of formation horizon full hexahedron mesh generation quality evaluation according to an embodiment of the present invention, in this embodiment, values of unitized Jacobian matrix determinant of 8 mesh nodes in each hexahedron unit in the hexahedron mesh are calculated, and a value of a minimum Jacobian matrix determinant of each hexahedron unit is determined, where a value range of the value of the minimum Jacobian matrix determinant is-1 to 1; and judging the degree that the value of the determinant of the minimum Jacobian matrix is close to 1, wherein the closer the value is to 1, the better the quality of the corresponding hexahedral unit is, and the value is less than zero, so that the corresponding hexahedral unit does not meet the quality requirement. Fig. 7 shows that the minimum value of the jacobian matrix determinant of all the hexahedral elements in the embodiment is greater than 0, and most values are greater than 0.6, which illustrates that the hexahedral elements obtained by the hexahedral mesh generation method provided by the embodiment of the present invention not only meet the requirement of finite element analysis, but also the generated hexahedral elements have very high quality.
In the embodiment of the invention, regular resampling is carried out on the actual stratum horizon to obtain the regular actual stratum horizon grid node coordinates; determining a straight line where each grid node in each layer and the corresponding grid node in the adjacent layer in the regularized actual stratum horizon are located, and the number of layers which are vertically split, and obtaining grid node coordinates between each layer and the adjacent layer; and generating a hexahedral mesh according to the mesh node coordinates between each layer and the adjacent layers. The method for generating the actual stratum full hexahedron mesh provided by the embodiment of the invention has the advantages of high automation degree, high speed and good mesh generation quality in the aspect of actual stratum generation, does not need human intervention in the generation process, and provides high-quality hexahedron mesh for industrial application of finite element forward simulation, reservoir inversion, ground stress calculation, oil reservoir numerical simulation and the like of an actual stratum. In the embodiment of the invention, the actual stratum layer is regularly resampled, so that the area of the actual stratum layer with blank area at the boundary is effectively processed; the method of the linear equation in the solid analytic geometry is adopted to obtain the grid node coordinates between each layer and the adjacent layer, the automation degree is high, the speed is high, the grid subdivision quality is good, and the subdivision process does not need human intervention. And hexahedron grids are generated according to the grid node coordinates between each layer and the adjacent layer, so that the automation degree is high and the speed is high.
Based on the same inventive concept, the embodiment of the invention also provides a device for generating the actual stratum full hexahedron grid, as described in the following implementation. Because the principles of solving the problems are similar to the actual stratum full hexahedron mesh generation method, the implementation of the device can refer to the implementation of the method, and repeated parts are not described in detail.
Fig. 8 is a schematic structural diagram of an actual stratum full hexahedron mesh generation apparatus in an embodiment of the present invention, and as shown in fig. 7, the apparatus includes:
the regular resampling module 801 is configured to perform regular resampling on an actual formation layer position to obtain a grid node coordinate of a regularized actual formation layer position;
a grid node generation preparation module 802, which determines the straight line where each grid node coordinate in each layer and the corresponding grid node coordinate in the adjacent layer are located in the regularized actual stratum layer;
the grid node generation module 803 determines the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction, and subdivides the space between each layer and adjacent layers in the regularized actual stratum layer according to the straight line and the number of subdivision layers to obtain the grid node coordinates of each subdivision layer between each layer and adjacent layers on the straight line;
and a hexahedral mesh generation module 804, configured to generate a hexahedral mesh according to the mesh node coordinates of each layer, the adjacent layers, and each split layer therebetween.
In an embodiment, the rule resampling module 801 may be specifically configured to:
obtaining the maximum value and the minimum value of the main survey line direction and the contact survey line direction of each layer in the actual stratum horizon;
obtaining a rectangular area of each layer according to the maximum value and the minimum value of the main survey line direction and the contact survey line direction, dividing the rectangular area of each layer into a plurality of rectangular sub-areas at equal intervals, and determining the rectangular area of each layer as each layer after regularization;
and resampling each layer after the regularization to obtain grid node coordinates of the regularized actual stratum layer.
In an embodiment, the grid node generating module 803 may be specifically configured to:
determining the number of subdivision layers of each layer and adjacent layers in the vertical direction in the regularized actual stratum horizon;
obtaining the thickness of each layer and adjacent layers on the straight line according to the number of the subdivision layers and the straight line;
calculating the coordinate of each section layer between each layer and the adjacent layer in the straight line vertical direction according to the thickness;
and calculating to obtain the grid nodes of each split layer between each layer and the adjacent layer on the straight line according to the coordinates of each split layer on the vertical direction of the straight line and the straight line.
In an embodiment, the hexahedral mesh generation module 804 may be specifically configured to:
numbering grid node coordinates of each layer, adjacent layers and each split layer between the adjacent layers;
forming a hexahedron unit by the eight adjacent grid nodes to obtain a hexahedron unit set;
and numbering the hexahedron units in the hexahedron unit set to generate hexahedron grids, and outputting grid node coordinates, grid node numbers, hexahedron unit numbers and grid node numbers forming the hexahedron units.
In an embodiment, the rule resampling module 801 may be further configured to: and resampling each layer after the regularization by adopting a nearest neighbor interpolation method to obtain grid node coordinates of the regularized actual stratum layer. Of course, the above nearest interpolation method is only an example, the grid node coordinate generation module may also adopt other methods to resample each layer after the regularization, and the military camp of the related variation example falls into the protection scope of the present invention.
In summary, the embodiment of the present invention performs regular resampling on the actual stratum horizon to obtain the regularized actual stratum horizon grid node coordinates; determining a straight line where each grid node in each layer and the corresponding grid node in the adjacent layer in the regularized actual stratum horizon are located, and the number of layers which are vertically split, and obtaining grid node coordinates between each layer and the adjacent layer; and generating a hexahedral mesh according to the mesh node coordinates between each layer and the adjacent layers. The method for generating the actual stratum full hexahedron mesh provided by the embodiment of the invention has the advantages of high automation degree, high speed and good mesh generation quality in the aspect of actual stratum generation, does not need human intervention in the generation process, and provides high-quality hexahedron mesh for industrial application of finite element forward simulation, reservoir inversion, ground stress calculation, oil reservoir numerical simulation and the like of an actual stratum. In the embodiment of the invention, the actual stratum layer is regularly resampled, so that the area of the actual stratum layer with blank area at the boundary is effectively processed; the method of the linear equation in the solid analytic geometry is adopted to obtain the grid node coordinates between each layer and the adjacent layer, the automation degree is high, the speed is high, the grid subdivision quality is good, and the subdivision process does not need human intervention. And hexahedron grids are generated according to the grid node coordinates between each layer and the adjacent layer, so that the automation degree is high and the speed is high.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A method for generating a full hexahedron grid of an actual stratum is characterized by comprising the following steps:
regular resampling is carried out on the actual stratum layer position, and grid node coordinates of the regular actual stratum layer position are obtained;
determining the straight line of each grid node coordinate in each layer in the regularized actual stratum horizon and the corresponding grid node coordinate in the adjacent layer;
determining the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction, and subdividing the space between each layer and adjacent layers in the regularized actual stratum layer according to the straight line and the number of subdivision layers to obtain the grid node coordinates of each subdivision layer between each layer and adjacent layers on the straight line;
generating hexahedral meshes according to the mesh node coordinates of each layer, the adjacent layers and each split layer between the two layers;
regular resampling is carried out on the actual stratum layer position to obtain the regular actual stratum layer position grid node coordinate, and the method comprises the following steps: obtaining the maximum value and the minimum value of the main survey line direction and the contact survey line direction of each layer in the actual stratum horizon; obtaining a rectangular area of each layer according to the maximum value and the minimum value of the main survey line direction and the contact survey line direction, dividing the rectangular area of each layer into a plurality of rectangular sub-areas at equal intervals, and determining the rectangular area of each layer as each layer after regularization; and resampling each layer after the regularization to obtain grid node coordinates of the regularized actual stratum layer.
2. The method for generating a full hexahedral mesh of an actual stratum according to claim 1, wherein the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction is determined, and the space between each layer and adjacent layers in the regularized actual stratum layer is subdivided according to the straight line and the number of subdivision layers to obtain the mesh node coordinates of each subdivision layer between each layer and adjacent layers on the straight line, and the method comprises the following steps:
determining the number of subdivision layers of each layer and adjacent layers in the vertical direction in the regularized actual stratum horizon;
obtaining the thickness of each layer and adjacent layers on the straight line according to the number of the subdivision layers and the straight line;
calculating the coordinate of each section layer between each layer and the adjacent layer on the straight line vertical direction according to the thickness;
and calculating to obtain the grid nodes of each section layer between each layer and the adjacent layer on the straight line according to the coordinates of each section layer on the straight line vertical direction and the straight line.
3. The method for generating a full hexahedral mesh of a practical stratigraphic layer according to claim 1, wherein generating the hexahedral mesh according to the mesh node coordinates of each layer and adjacent layers and each hierarchical layer therebetween, comprises:
numbering grid node coordinates of each layer, adjacent layers and each split layer between the adjacent layers;
forming a hexahedron unit by the eight adjacent grid nodes to obtain a hexahedron unit set;
and numbering the hexahedron units in the hexahedron unit set to generate hexahedron grids, and outputting grid node coordinates, grid node numbers, hexahedron unit numbers and grid node numbers forming the hexahedron units.
4. The method of claim 1, wherein the regularized layers are resampled by a nearest neighbor interpolation method to obtain the grid node coordinates of the regularized actual formation layer.
5. An actual formation full hexahedron mesh generation device, comprising:
the regular resampling module is used for carrying out regular resampling on the actual stratum layer position to obtain a grid node coordinate of the regular actual stratum layer position;
the grid node generation preparation module is used for determining the coordinates of each grid node in each layer in the regularized actual stratum horizon and the straight line where the coordinates of the corresponding grid node in the adjacent layer are located;
the grid node generation module is used for determining the number of subdivision layers of each layer and adjacent layers in the regularized actual stratum layer in the vertical direction, subdividing the space between each layer and adjacent layers in the regularized actual stratum layer according to the straight line and the number of subdivision layers, and obtaining the grid node coordinates of each subdivision layer between each layer and adjacent layers on the straight line;
the hexahedral mesh generation module is used for generating hexahedral meshes according to the mesh node coordinates of each layer, the adjacent layers and each split layer between the two layers;
the rule resampling module is specifically configured to: obtaining the maximum value and the minimum value of the main survey line direction and the contact survey line direction of each layer in the actual stratum horizon; obtaining a rectangular area of each layer according to the maximum value and the minimum value of the main survey line direction and the contact survey line direction, dividing the rectangular area of each layer into a plurality of rectangular sub-areas at equal intervals, and determining the rectangular area of each layer as each layer after regularization; and resampling each layer after the regularization to obtain grid node coordinates of the regularized actual stratum layer.
6. The actual stratigraphic full hexahedron mesh generation apparatus of claim 5, wherein the mesh node generation module is specifically configured to:
determining the number of subdivision layers of each layer and adjacent layers in the vertical direction in the regularized actual stratum horizon;
obtaining the thickness of each layer and adjacent layers on the straight line according to the number of the subdivision layers and the straight line;
calculating the coordinate of each section layer between each layer and the adjacent layer on the straight line vertical direction according to the thickness;
and calculating to obtain the grid nodes of each section layer between each layer and the adjacent layer on the straight line according to the coordinates of each section layer on the straight line vertical direction and the straight line.
7. The actual stratigraphic full hexahedron mesh generation apparatus of claim 5, wherein the hexahedron mesh generation module is specifically configured to:
numbering grid node coordinates of each layer, adjacent layers and each split layer between the adjacent layers;
forming a hexahedron unit by the eight adjacent grid nodes to obtain a hexahedron unit set;
and numbering the hexahedron units in the hexahedron unit set to generate hexahedron grids, and outputting grid node coordinates, grid node numbers, hexahedron unit numbers and grid node numbers forming the hexahedron units.
8. The actual stratigraphic full hexahedral mesh generation apparatus of claim 5, wherein the rule resampling module is further for: and resampling each layer after the regularization by adopting a nearest neighbor interpolation method to obtain grid node coordinates of the regularized actual stratum layer.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 4 when executing the computer program.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811241872.1A CN109584369B (en) | 2018-10-24 | 2018-10-24 | Actual stratum full hexahedron grid generation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811241872.1A CN109584369B (en) | 2018-10-24 | 2018-10-24 | Actual stratum full hexahedron grid generation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109584369A CN109584369A (en) | 2019-04-05 |
| CN109584369B true CN109584369B (en) | 2021-09-28 |
Family
ID=65920434
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811241872.1A Active CN109584369B (en) | 2018-10-24 | 2018-10-24 | Actual stratum full hexahedron grid generation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109584369B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110188395B (en) * | 2019-05-05 | 2021-06-15 | 广东海洋大学 | Dimension increasing type calculation fluid grid generation method based on linear surface body |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5729670A (en) * | 1996-01-16 | 1998-03-17 | Ford Global Technologies, Inc. | Method for producing a mesh of quadrilateral/hexahedral elements for a body to be analyzed using finite element analysis |
| US7671858B1 (en) * | 2005-09-06 | 2010-03-02 | Sandia Corporation | Unconstrained paving and plastering method for generating finite element meshes |
| CN103236088A (en) * | 2013-04-16 | 2013-08-07 | 北京金阳普泰石油技术股份有限公司 | Three-dimensional geological structure modeling method and system based on composite grid |
| CN105302951A (en) * | 2015-10-20 | 2016-02-03 | 北京工业大学 | Finite element mesh surface triangulation method |
| CN106251401A (en) * | 2016-07-06 | 2016-12-21 | 西安交通大学 | A kind of generation system and method for hexahedron 20 node unit grid |
| CN106652033A (en) * | 2016-12-05 | 2017-05-10 | 中国石油天然气股份有限公司 | Method for subdividing natural grid of geological profile |
-
2018
- 2018-10-24 CN CN201811241872.1A patent/CN109584369B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5729670A (en) * | 1996-01-16 | 1998-03-17 | Ford Global Technologies, Inc. | Method for producing a mesh of quadrilateral/hexahedral elements for a body to be analyzed using finite element analysis |
| US7671858B1 (en) * | 2005-09-06 | 2010-03-02 | Sandia Corporation | Unconstrained paving and plastering method for generating finite element meshes |
| CN103236088A (en) * | 2013-04-16 | 2013-08-07 | 北京金阳普泰石油技术股份有限公司 | Three-dimensional geological structure modeling method and system based on composite grid |
| CN105302951A (en) * | 2015-10-20 | 2016-02-03 | 北京工业大学 | Finite element mesh surface triangulation method |
| CN106251401A (en) * | 2016-07-06 | 2016-12-21 | 西安交通大学 | A kind of generation system and method for hexahedron 20 node unit grid |
| CN106652033A (en) * | 2016-12-05 | 2017-05-10 | 中国石油天然气股份有限公司 | Method for subdividing natural grid of geological profile |
Non-Patent Citations (2)
| Title |
|---|
| 基于典型焊接构件模型的有限元六面体网格自动划分算法研究;申刚;《中国优秀硕士学位论文全文数据库》;20140615;第5.3.2-5.3.3节 * |
| 基于几何分解的六面体网格生成算法的研究;王东风;《中国优秀硕士学位论文全文数据库》;20111215;第4.2-4.6节，图4.6 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109584369A (en) | 2019-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106327577B (en) | Dimensional topography optimal method based on local curvature's entropy and quad-tree structure | |
| RU2282242C2 (en) | Method for evening out a grid of nodes in accordance to characteristic features in digital image | |
| Natali et al. | Modeling Terrains and Subsurface Geology. | |
| EP0951652B1 (en) | Performing geoscience interpretation with simulated data | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| US10795053B2 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| US9390556B2 (en) | Systems and methods for generating a large scale polygonal mesh | |
| CN106687827B (en) | Stratum modeling method for fault | |
| WO2017007924A1 (en) | Improved geobody continuity in geological models based on multiple point statistics | |
| CN101582173A (en) | Block model construction method for complex geological structure | |
| CN104504754A (en) | Multipoint statistical modeling method and device | |
| CN103606191A (en) | Method for establishing complex geologic model quickly | |
| CN102867332A (en) | Multi-level subdivided mesh surface fitting method based on complicated boundary constraint | |
| CN116030218A (en) | Tetrahedral mesh dividing method, device, system and storage medium | |
| CN111445569B (en) | Sedimentary geological evolution dynamic simulation method | |
| CN109584369B (en) | Actual stratum full hexahedron grid generation method and device | |
| CN111968227A (en) | Three-dimensional geological fault network uncertainty analysis method, system and storage medium | |
| CN103207410A (en) | Rugged seabed aimed hybrid grid model building method | |
| WO2018158378A1 (en) | Computer implemented method for improving a velocity model for seismic imaging | |
| CN117635854A (en) | Knowledge-element-based geological model local updating method and system | |
| CN110968930B (en) | Geological variable attribute interpolation method and system | |
| US10254441B2 (en) | Method of modelling a subsurface volume | |
| CN113484909B (en) | Method for establishing fracture-cavity reservoir based on geometric gridding and parameter distribution | |
| Kononov et al. | A 3D tetrahedral mesh generator for seismic problems | |
| CN108876916B (en) | Braided river training image generation method and equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |