CN103236087A - Triangular prism-shaped geological model construction method - Google Patents
Triangular prism-shaped geological model construction method Download PDFInfo
- Publication number
- CN103236087A CN103236087A CN2013101475408A CN201310147540A CN103236087A CN 103236087 A CN103236087 A CN 103236087A CN 2013101475408 A CN2013101475408 A CN 2013101475408A CN 201310147540 A CN201310147540 A CN 201310147540A CN 103236087 A CN103236087 A CN 103236087A
- Authority
- CN
- China
- Prior art keywords
- node
- layer
- grid
- nodes
- unit
- 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.)
- Granted
Links
Images
Landscapes
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a triangular prism-shaped geological model construction method. The method includes the following steps: (1) a model projection region is set and divided into planar grid units, and nodes in a model bottom control region are divided into nodes inside the control region, nodes outside the control region and nodes on the control boundary; (2) each layer of geological materials are divided into a grid layer, whether the nodes are located on the boundary surface of the materials of the layer is judged, the coordinates of the nodes of the materials of the layer higher than the boundary surface of the layer are vertically moved down to the boundary surface of the materials, the nodes of each layer are numbered, and the control boundary and the region interior are specially processed; (3) triangular prism units k, l, m, k plus np0, 1 plus np0 and m plus np0 are sequentially constructed on virtual nodes corresponding to the ith layer and the i plus 1th layer until the ground surface; (4) real node numbering is carried out on the nodes from the first layer, whether the nodes are superposed with the corresponding nodes of the lower layer is judged, if so, then the nodes are numbered by adopting the corresponding real nodes of the lower layer, or else the nodes are new real nodes; and real node numbering is carried out; (5) the virtual nodes forming the units are replaced by the real nodes, and thereby a geological model is generated.
Description
Technical field
The present invention relates to a kind of construction method of triangular prism geologic model.
Background technology
Numerical simulation technology is the important means of civil engineering work mechanical analysis, and it then is the chief component of numerical simulation pre-treatment that geologic model makes up.Yet, under the influence of the complicated geological origin cause of formation, the grid irregular appears to the geologic model that stratification distributes in normal existence because of the crucial weak formation of part when carrying out the finite element grid subdivision, the thin layer stratum is difficult to the accurately phenomenon of structure, causes model to realize difficulty, result of calculation distortion.Simultaneously, because complex layered stratum can only adopt tetrahedral grid to divide usually, its computing velocity, precision far are inferior to hexahedron, triangular prism grid etc., therefore, optimize the structure model method and have the engineering application value.
Summary of the invention
Goal of the invention: the objective of the invention is to the layered geology model object at fields such as civil engineering works, provide a kind of can be from bottom to top, based on the corresponding construction method that generates the triangular prism geologic model fast of node actual situation.
Technical scheme: the construction method of triangular prism geologic model of the present invention comprises the steps:
Step 1 arranges control area, model bottom surface, and be the triangular mesh unit (Fig. 1) that is uneven in length with model bottom surface subdivision, then with crucial weak formation range of control to the projection of model bottom surface, thereby control area, model bottom surface interior nodes is divided into 3 classes: crucial control zone interior nodes (Fig. 2), crucial control zone exterior node and pass key control boundary node (Fig. 3), control area, described model bottom surface node total number is np0, and the unit total number is ne0;
Step 2 generates the boundary surface reference mark that the level geologic model needs known models bottom surface elevation and layers of material, yet in actual geology modeling process, complicated geo-logical terrain boundary surface is difficult to determine, adopts the geotechnical boring exploration usually, and interpolation obtains the earth formation material interphase according to exploring as a result.The present invention sets up on known layer boundary surface basis.Generate the triangular prism grid node according to formation information by layer, and give dummy node to it and number, definition model bottom surface elevation z0, determine each earth formation material interphase and ground surface from bottom to top successively respectively, then the 1st layer of stratum be between model bottom surface and the 1st layer of strata interface, the 2nd layer material between the 1st layer of strata interface and the 2nd layer of strata interface ... and the like, N layer stratum is between N-1 layer strata interface and N layer strata interface; Adopt length d z1 that the triangular prism between the 1st layer is carried out dividing elements, if the triangular prism node maximum elevation difference of stratigraphic boundary correspondence is dzz1 up and down, the grid number of plies N1=dzz1/dz1 of this layer stratigraphic division then, N1 adopts into, and a method rounds, so the 1st layer of stratum can be controlled by N1 layer grid cell, each layer grid node number is np0, and number of unit is ne0, so the numbering of j node of i layer can be designated as (i-1) * np0+j in this stratum grid; Travel through all grid nodes of the 1st layer material, each node and current stratigraphic boundary are compared, if node is positioned at the top, stratigraphic boundary, node is moved on the current strata interface along its place triangular prism rib direction; In like manner, when carrying out the grid division for k layer stratum, be the bottom surface with coboundary, k-1 layer stratum, adopt length d zk that the triangular prism between the k layer is carried out dividing elements, if the triangular prism node maximum elevation difference of stratigraphic boundary correspondence is dzzk up and down, the grid number of plies Nk=dzzk/dzk of k layer stratigraphic division then, Nk adopts into, and a method rounds, so k layer stratum can be controlled by Nk layer grid cell, because each layer grid node number is similarly np0, number of unit also is similarly ne0, so the numbering of m node of l layer can be designated as (l-1) * np0+m in the grid of k layer stratum; Travel through all grid node of k layer, ratio carried out in each node and coboundary, k layer stratum, if node is positioned at coboundary, k layer stratum top, with node along its place triangular prism rib direction translation to current strata interface; From traveling through all stratum down, can generate geologic model successively;
When crucial weak formation being carried out the grid cell structure, need grid node outside grid node, the crucial control zone in the crucial control zone and close the key control boundary node to handle respectively, it is characterized in that, in the crucial control zone during grid node division, crucial weak formation lower boundary is the coboundary of one deck under the bottom, the coboundary then is lower boundary or the ground surface on the top, outside the crucial control zone during grid node division, crucial weak formation lower boundary is identical with the coboundary, close key control boundary node position, upper and lower border, stratum is also identical; Be as the criterion with crucial weak formation control border and only carry out grid node at all weak formation triangular prisms in the control area and divide, can generate the thin layer element that thickness is d, namely control soft stratum;
Step 3 self model bottom surface z0 begins, and is controlled condition with model bottom surface triangular mesh, carries out cell formation according to i layer and i+1 layer dummy node, and advance to the face of land successively, if i layer plane unit grid below node is k, l, m, then the node of this grid cell correspondence is k, l, m, k+np0, l+np0, m+np0; Because the cell node coordinate has coincidence, forms so the unit can be regarded as by dummy node, the dummy node total number is (NT+1) * np0, and NT is the total number of plies of geologic model grid;
Step 4 self model bottom surface z0 begins, again the grid dummy node is carried out the physical node numbering, to j node of i layer, judge whether it vertically overlaps with the corresponding node of i-1 layer on the edge, this node is designated as the last layer node serial number if overlap then, otherwise this node is new physical node, is numbered, travel through all grid dummy nodes, finally can obtain physical node sum NP; Travel through all grid cells, each unit is carried out volume calculate, if unit volume is 0, then reject this unit, if it is volume greater than 0, then carries out element number to this unit, final that real unit total number is NE.
All grid cells of step 5 traversal are replaced the dummy node in each real unit with the corresponding physical node numbering of this point, then node unit information is pressed specific format output, then can generate corresponding triangular prism geologic model.
By adopting above-mentioned technology, the time that the present invention generates the complicated geological model can shorten to dozens of minutes by several days even several months, has saved a large amount of quality time.In addition, adopt the present invention, the very regular weak thin layer in generating mesh unit flexibly also, grid cell is triangular prism, has not only reduced the grid cell number and has also reduced lopsided grid cell, thereby guaranteed the reliability of numerical result.
The present invention compared with prior art, its beneficial effect is: 1, make up model with node actual situation corresponded manner, solved the problem that crucial weak formation model is realized difficulty, greatly improved modeling efficiency, and modeling is easy, quick; 2, grid cell adopts triangular prism, avoided needing situation with a large amount of tetrahedron modelings because of the cast material complex boundary, than the tetrahedral grid model of element, the about tetrahedral grid model of element of sum unit, its unit total half, thereby shorten computing time greatly; 3, this modeling method is not only applicable to the structure of geologic model, and suitable equally for the model construction on the material with complex border of other field, so its scope of application is extensive, controllability is stronger, and the model accuracy of foundation is higher.
In a word, the present invention from bottom to top, based on the corresponding geologic model construction method that generates the triangular prism grid cell fast of node actual situation, the geomechanics model grid that makes up is more regular than tetrahedron, it is convenient, accurate that the weak floor grid is realized, the node unit number falls sharply, computing velocity is fast, and the result is more accurate.
Description of drawings
Fig. 1 is model bottom surface control mesh figure.
The crucial weakness zone perspective view of Fig. 2.
Fig. 3 is model control area boundary graph.
Fig. 4 is simplified model terrestrial net lattice and control boundary graph.
Fig. 5 is that simplified model ground floor material generates figure.
Fig. 6 is simplified model ground floor material correction figure.
Fig. 7 is that figure is revised on simplified model materials at two layers border.
Fig. 8 is simplified model layered geology grid chart
Fig. 9 is complex model terrestrial net trrellis diagram.
Figure 10 is complex model control area boundary graph.
Figure 11 is complex model layered geology grid chart.
Figure 12 is crucial weak floor grid chart.
Figure 13 is layered geology model product process figure.
Embodiment
In order to understand the modeling technique of layered geology of the present invention easily, below respectively from adopting simple and complicated two model instances to describe.
Embodiment 1:(naive model)
Certain simplifies geologic model bottom surface elevation z0=0m step 1, its base surface area is a rectangle, long 100m, wide 60m, bold box is that crucial weak formation range of control is in the projection of model bottom surface in the rectangular area, at first the every 10m in model bottom surface is divided a square net, connect its diagonal line then successively, but then each little square net again subdivision be two triangular mesh, finally obtain np0=77 of model base surface area node sum, npe=120 of triangular mesh unit sum, as shown in Figure 4.
Step 2 is looked unfamiliar into the ground floor material according to the model bottom surface with the coboundary of ground floor material, clathrum spacing dz1=10m is set, but be dzz1=75m according to the maximum interlamellar spacing of known materials boundary surface computation model bottom surface and ground floor material interface, the number of plies N1=dzz1/dz1=75/10=7.5 that divides of ground floor material then, get N1=8 by advancing a method, the self model bottom surface upwards generates 8 layers of control mesh layer successively, every layer of grid node is 77, flat unit is 120, then the numbering of j node of i layer can be expressed as (i-1) * np0+j, and cell formation is carried out to i layer and i+1 layer dummy node successively in the self model bottom surface.Grid node as i layer plane trigonometry unit is k, l, m; The grid corresponding node of i+1 layer respective planes triangular unit is k+np0, l+np0, and m+np0, then this triangular prism grid cell is by node k, l, m, k+np0, l+np0, m+np0, ground floor web material lattice is built as shown in Figure 5.Compare by node and interface, ground floor material top for 8 layers of control mesh layer that generate, will be higher than the node at interface, material top, move to interface, material top straight down, as shown in Figure 6.
Step 3 is that the bottom surface generates second layer material according to the interface, top of ground floor material, clathrum spacing dz2 still is set to 10m, the maximal value of the corresponding coboundary with it of node, the maximum interlamellar spacing dzz2=48m(second layer material bottom surface of second layer material nodal distance as calculated), divide number of plies N2=dzz2/dz2=28/10=2.8, get N2=3, upwards generate 3 layers of control mesh layer successively from second layer material bottom surface, all nodes of traversal second layer material, the node that will be higher than its interface, top vertically is displaced downwardly on its interface, top, and net result as shown in Figure 7.
The generation of the crucial weak key-course of step 4, the crucial weak floor of this model is positioned at ground surface for simplicity, namely the 3rd layer.To in the crucial control zone, outside the crucial control zone and close the borderline grid node of key control and handle respectively, coboundary in this routine crucial control zone is ground surface, lower boundary is the coboundary of second layer material, overlap up and down with the boundary surface face that closes the key control border outside the crucial control zone, be second layer material coboundary.According to step 3, crucial weak formation is generating mesh key-course in crucial control zone only.
Step 5 self model bottom surface ground floor carries out the physical node numbering to grid node, to j node of i layer, judge whether it overlaps with the corresponding node of i-1 layer, if overlapping then, this node is designated as the real numbering of last layer node, otherwise this node is new physical node, be numbered, finally can obtain NP=952 of physical node sum; Travel through all earth formation material and ask the volume of each unit, reject this unit if unit volume is 0, if volume is this element greater than 0 and carries out element number, final NE=1164 of the real unit sum that get.
All earth formation material of step 6 traversal are replaced the dummy node in each unit with the corresponding physical node numbering of this point, at last node unit information is pressed specific format output, then generate corresponding triangular prism geologic model as shown in Figure 8.
Embodiment 2:(complex model)
Certain stratiform complicated geological model divides three layers, and the deposit on the displacement face of land is key-course, and its exposure limit, face of land such as Fig. 9 have the thick soft layer of one deck 0.2m in this deposit bottom.It is as follows to use modeling procedure of the present invention:
(1) control area, model bottom surface at first is set as shown in Figure 9, be the triangular mesh unit that is uneven in length with model bottom surface subdivision, can obtain np0=8000 of bottom layer node sum, npe=15642 of triangular element sum, then with crucial weak formation range of control to the projection of model bottom surface, be three classes according to view field with model bottom surface node division: crucial control zone interior nodes, crucial control zone boundary node and crucial control zone exterior node generate as shown in figure 10 with the control that makes things convenient for weak floor.
(2) look unfamiliar into earth formation material according to the known formation material boundary, every layer of clathrum spacing dzi=100m is set, the self model bottom surface makes progress the maximum interlamellar spacing of every layer material apart from being followed successively by 1495,192,2m, then the number of plies of Hua Fening is followed successively by 15,2,1, generate the control mesh layer successively from bottom to top, according to the boundary surface of known layers of material, in layer of material, travel through all unit with the height node at interface, layer material top therewith, move to straight down on the material boundary face.Carry out special processing for weak floor, namely generate the material of one deck thickness d d=20cm.
(3) self model bottom surface ground floor begins node is carried out the physical node numbering, for j node of i layer, judge whether it overlaps with the i-1 node layer, then this puts corresponding last layer node serial number if overlap, otherwise this node is new physical node, carries out the physical node numbering.Finally, can obtain 89739 of physical node number N P=; Travel through all unit to each unit cube, reject this unit if unit volume is 0, be real unit greater than 0 and carry out element number as if volume, final that real unit total number is 161934 of NE=.
(4) travel through all earth formation material, dummy node in each unit is replaced with the corresponding physical node numbering of this point, at last node unit information is pressed specific format output, then can generate corresponding triangular prism geologic model as shown in figure 11, the weak layer of crucial soft control as shown in figure 12.
Base area layer material boundary surface and model bottom surface elevation, be the triangular mesh unit that is uneven in length with model bottom surface subdivision, then with crucial weak formation range of control to the projection of model bottom surface, thereby control area, model bottom surface interior nodes is divided into 3 classes, obtain control area, model bottom surface node total number np0, unit total number ne0.The grid key-course is upwards set up in the self model bottom surface, until ground surface, base area layer material boundary surface will vertically be displaced downwardly to boundary surface with the node on the layer material face, then to the real numbering of layers of material node, the rejecting unit volume is 0 unit, and each unit dummy node is replaced with physical node in the model the most at last, and actual situation is mapped, can set up the geologic model that meets the following requirements: 1. grid cell controllable size is even and regular; 2. the generation of weak floor control is very flexible; 3. model unit is triangular prism, and lopsided unit is few, and computational accuracy is the four sides height; 4. modeling applicability is extensive, saves time, and model is more complicated, and it is more many to save time.The node unit information that obtains by this technology modeling can directly import software for calculation such as large-scale finite element ansys, finite difference flac3d by specific format output, for user modeling provides great convenience.
As mentioned above, although represented and explained the present invention that with reference to specific preferred embodiment it shall not be construed as the restriction to the present invention self.Under the spirit and scope of the present invention prerequisite that does not break away from the claims definition, can make various variations in the form and details to it.
Claims (2)
1. the construction method of a triangular prism geologic model is characterized in that comprising the steps:
(1) control area, model bottom surface is set, and be the triangular mesh unit that is uneven in length with model bottom surface subdivision, then with crucial weak formation range of control to the projection of model bottom surface, thereby control area, model bottom surface interior nodes is divided into 3 classes: crucial control zone interior nodes, crucial control zone exterior node and pass key control boundary node, control area, described model bottom surface node total number is np0, and the unit total number is ne0;
(2) generate the triangular prism grid node according to formation information by layer, and give the dummy node numbering to it, definition model bottom surface elevation z0 determines each earth formation material interphase and ground surface from bottom to top successively respectively; Adopt length d z1 that the triangular prism between the 1st layer is carried out dividing elements, if the triangular prism node maximum elevation difference of stratigraphic boundary correspondence is dzz1 up and down, the grid number of plies N1=dzz1/dz1 of this layer stratigraphic division then, N1 adopts into, and a method rounds, so the 1st layer of stratum can be controlled by N1 layer grid cell, each layer grid node number is np0, and number of unit is ne0, so the numbering of j node of i layer can be designated as (i-1) * np0+j in this stratum grid; Travel through all grid nodes of the 1st layer material, each node and current stratigraphic boundary are compared, if node is positioned at the top, stratigraphic boundary, node is moved on the current strata interface along its place triangular prism rib direction; In like manner, when carrying out the grid division for k layer stratum, be the bottom surface with coboundary, k-1 layer stratum, adopt length d zk that the triangular prism between the k layer is carried out dividing elements, if the triangular prism node maximum elevation difference of stratigraphic boundary correspondence is dzzk up and down, the grid number of plies Nk=dzzk/dzk of k layer stratigraphic division then, Nk adopts into, and a method rounds, so k layer stratum can be controlled by Nk layer grid cell, because each layer grid node number is similarly np0, number of unit also is similarly ne0, so the numbering of m node of l layer can be designated as (l-1) * np0+m in the grid of k layer stratum; Travel through all grid node of k layer, ratio carried out in each node and coboundary, k layer stratum, if node is positioned at coboundary, k layer stratum top, with node along its place triangular prism rib direction translation to current strata interface; From traveling through all stratum down, namely generate geologic model successively;
(3) self model bottom surface z0 begins, and is controlled condition with model bottom surface triangular mesh, carries out cell formation according to i layer and i+1 layer dummy node, and advance to the face of land successively, if i layer plane unit grid below node is k, l, m, then the node of this grid cell correspondence is k, l, m, k+np0, l+np0, m+np0; Because the cell node coordinate has coincidence, forms so the unit can be regarded as by dummy node, the dummy node total number is (NT+1) * np0, and NT is the total number of plies of geologic model grid;
(4) self model bottom surface z0 begins, again the grid dummy node is carried out the physical node numbering, to j node of i layer, judge whether it vertically overlaps with the corresponding node of i-1 layer on the edge, this node is designated as the last layer node serial number if overlap then, otherwise this node is new physical node, is numbered, travel through all grid dummy nodes, finally can obtain physical node sum NP; Travel through all grid cells, each unit is carried out volume calculate, if unit volume is 0, then reject this unit, if it is volume greater than 0, then carries out element number to this unit, final that real unit total number is NE;
(5) travel through all grid cells, the dummy node in each real unit is replaced with the corresponding physical node numbering of this point, then node unit information is pressed specific format output, then generate corresponding triangular prism geologic model.
2. the construction method of triangular prism geologic model according to claim 1, it is characterized in that: in the step (2), when crucial weak formation is carried out the grid cell structure, with grid node in the crucial control zone, grid node and pass key control boundary node are handled respectively outside the crucial control zone, in the crucial control zone during grid node division, crucial weak formation lower boundary is the coboundary of one deck under the bottom, the coboundary then is lower boundary or the ground surface on the top, outside the crucial control zone during grid node division, crucial weak formation lower boundary is identical with the coboundary, close key control boundary node position, on the stratum, lower boundary is also identical; Be as the criterion with crucial weak formation control border and only carry out grid node at all weak formation triangular prisms in the control area and divide, can generate the thin layer element that thickness is d, namely control soft stratum.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310147540.8A CN103236087B (en) | 2013-04-25 | 2013-04-25 | A kind of construction method of triangular prism geologic model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310147540.8A CN103236087B (en) | 2013-04-25 | 2013-04-25 | A kind of construction method of triangular prism geologic model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103236087A true CN103236087A (en) | 2013-08-07 |
| CN103236087B CN103236087B (en) | 2015-09-16 |
Family
ID=48884126
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310147540.8A Active CN103236087B (en) | 2013-04-25 | 2013-04-25 | A kind of construction method of triangular prism geologic model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103236087B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106373193A (en) * | 2016-09-22 | 2017-02-01 | 河海大学 | Box covering method-based slope assembly type discrete element model generating method |
| CN108197360A (en) * | 2017-12-22 | 2018-06-22 | 哈尔滨汽轮机厂有限责任公司 | The automatic dividing system of turbine rotor grid and method |
| CN111414654A (en) * | 2020-03-16 | 2020-07-14 | 南京交通职业技术学院 | Grid-division continuous numerical model retrieval method |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102063555A (en) * | 2011-01-26 | 2011-05-18 | 河海大学 | Finite element numerical model debugging method based on grid structure |
-
2013
- 2013-04-25 CN CN201310147540.8A patent/CN103236087B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102063555A (en) * | 2011-01-26 | 2011-05-18 | 河海大学 | Finite element numerical model debugging method based on grid structure |
Non-Patent Citations (4)
| Title |
|---|
| JING MING等: "GSIS: A 3D geological multi-body modeling system from netty cross-sections with topology", 《COMPUTERS & GEOSCIENCES》, vol. 36, 31 December 2010 (2010-12-31), pages 756 - 767, XP027059005 * |
| 周坤等: "三维地质构造的三角形格架式网格生成方法", 《计算机工程与设计》, vol. 32, no. 4, 31 December 2011 (2011-12-31) * |
| 孟宪海等: "三维地质体类三棱柱网格建模方法研究", 《计算机工程与应用》, vol. 46, no. 4, 31 December 2010 (2010-12-31) * |
| 芮小平等: "一种基于三棱柱的三维地质体可视化方法研究", 《中国矿业大学学报》, vol. 33, no. 5, 30 September 2004 (2004-09-30), pages 584 - 588 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106373193A (en) * | 2016-09-22 | 2017-02-01 | 河海大学 | Box covering method-based slope assembly type discrete element model generating method |
| CN106373193B (en) * | 2016-09-22 | 2019-03-01 | 河海大学 | A kind of side slope assembled discrete element analysis generation method based on box cladding process |
| CN108197360A (en) * | 2017-12-22 | 2018-06-22 | 哈尔滨汽轮机厂有限责任公司 | The automatic dividing system of turbine rotor grid and method |
| CN108197360B (en) * | 2017-12-22 | 2021-07-13 | 哈尔滨汽轮机厂有限责任公司 | Automatic division system and method for steam turbine rotor grids |
| CN111414654A (en) * | 2020-03-16 | 2020-07-14 | 南京交通职业技术学院 | Grid-division continuous numerical model retrieval method |
| CN111414654B (en) * | 2020-03-16 | 2024-02-23 | 南京交通职业技术学院 | Method for retrieving continuous numerical model of division |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103236087B (en) | 2015-09-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103514631B (en) | Magnanimity three-dimensional geological model grid type parallel constructing method based on borehole data | |
| US20140136171A1 (en) | Unstructured Grids For Modeling Reservoirs | |
| US8463586B2 (en) | Machine, program product, and computer-implemented method to simulate reservoirs as 2.5D unstructured grids | |
| 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 | |
| CN110529077B (en) | Modeling cross-faults and complex wellbores in reservoir simulation | |
| NO309548B1 (en) | Method for generating a three-dimensional lattice pattern corresponding to the geometry of a body, to obtain a model representing the body | |
| CN107944086B (en) | Rapid modeling method based on drilling data | |
| EA022006B1 (en) | Method of geophysical survey of prospective oil-bearing area | |
| CN105513131B (en) | A kind of three-dimensional geological structural model Corner-point Grids subdivision method taking fold into account | |
| US10174593B2 (en) | System and method for coarsening in reservoir simulation system | |
| EP3191874B1 (en) | A method of stratigraphic modeling of faults | |
| CN103226847B (en) | A kind of meandering river point bar 3-D data volume generation method and device | |
| US20120022837A1 (en) | Smoothing Of Stair-Stepped Geometry In Grids | |
| CN104504754B (en) | A kind of method and device of oil and gas reservoir multi-point statistic modeling | |
| CN105184867A (en) | Three-dimensional stratum modeling method based on regular grids and corner-point grid technology | |
| CN113313825A (en) | Three-dimensional geological modeling method and device for oil reservoir | |
| CN102194252A (en) | Geological-stratum-structure-based method for generating triangular lattice grids | |
| NO20121440A1 (en) | Rotations from gradient directions | |
| CN107481320A (en) | A kind of gridless routing of three-dimensional geological modeling | |
| CN103678899B (en) | A kind of oil-gas reservoir multi-point statistic modeling method based on variable effect ratio and equipment | |
| CN103236087A (en) | Triangular prism-shaped geological model construction method | |
| CN104535391B (en) | A kind of Geophysical Data Processing method based on layered geology model | |
| WO2014154538A1 (en) | Method of modelling a subsurface volume | |
| CN109215464B (en) | Trap and oil reservoir three-dimensional conceptual model construction method | |
| Fremming | 3D geological model construction using a 3D grid |
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 |