CN110992488A - Inclined crack grid dividing method based on embedded discrete crack model - Google Patents
Inclined crack grid dividing method based on embedded discrete crack model Download PDFInfo
- Publication number
- CN110992488A CN110992488A CN201911328526.1A CN201911328526A CN110992488A CN 110992488 A CN110992488 A CN 110992488A CN 201911328526 A CN201911328526 A CN 201911328526A CN 110992488 A CN110992488 A CN 110992488A
- Authority
- CN
- China
- Prior art keywords
- crack
- fracture
- grid
- projection
- matrix
- 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
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/05—Geographic models
Abstract
The invention provides an inclined fracture grid dividing method based on an embedded discrete fracture model, which belongs to the technical field of oil reservoir simulation and comprises the following steps: establishing a matrix grid system; acquiring the boundary of the crack surface and the projection surface of the matrix grid system on each coordinate axis plane; inverting the projection of the matrix grid penetrated by the crack surface on each projection surface to obtain the matrix grid penetrated by the crack surface; dividing fracture surfaces by using the boundaries of the matrix grids to obtain initial fracture grids; the method and the device have the advantages that the fracture surface boundary is used as a constraint condition, the part, located in the fracture boundary, in the preliminary fracture grid is reserved, and the final fracture grid is obtained.
Description
Technical Field
The invention belongs to the technical field of oil reservoir simulation, and particularly relates to an inclined fracture grid subdivision method based on an embedded discrete fracture model.
Background
With the development of conventional oil and gas resources, conventional reservoirs with high porosity and high permeability have been produced almost completely. To meet the increasing demand for oil and gas resources, the development of unconventional low permeability reservoirs has become imperative. The mine practice shows that: the low-permeability compact reservoir has the characteristic of low permeability, so that a fracturing horizontal well technology is generally adopted. Because the conventional oil reservoir engineering method is difficult to deal with the complex fractured reservoir seepage, a numerical simulation method is mostly adopted to research the complex fractured reservoir seepage. In the numerical reservoir simulation method, the embedded discrete fracture model does not need to encrypt the grid, so that the embedded discrete fracture model has higher operation speed and is favored by more and more scholars. However, in the embedded discrete fracture model, the fracture mesh is obtained by dividing the boundary of the structured matrix mesh, and the intersection mode of the fracture plane and the matrix mesh in the three-dimensional model is very complex, so that the difficulty in dividing the inclined fracture mesh is very high, most students still adopt a simplified model of a vertical fracture model, but the simplified model has a large difference from the actual situation, and the simulation result cannot reflect the real stratum situation. The invention aims to provide an inclined crack grid dividing method based on an embedded discrete crack model, and solves the problem that the inclined crack is difficult to deal with by applying the embedded discrete crack model.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention provides an inclined crack grid subdivision method based on an embedded discrete crack model.
The technical scheme of the invention is as follows:
the invention provides an inclined crack mesh generation method based on an embedded discrete crack model, which comprises the following steps:
s1, subdividing the matrix mesh according to the basic geometric parameters of the mesh to establish a matrix mesh system;
s2, projecting the crack surface boundary and the matrix grid system to 3 coordinate axis planes respectively to obtain projection surfaces of the matrix grid penetrated by the crack surface in the 3 coordinate axis planes;
s3, inverting the projection of the matrix grid penetrated by the crack surface in each projection plane to obtain the matrix grid penetrated by the crack surface;
s4, dividing a fracture plane by using the boundary of the matrix grid penetrated by the fracture surface to obtain a primary fracture grid;
and S5, taking the central point in the crack surface, taking the crack boundary as a constraint condition, obtaining the point in the crack boundary in the preliminary crack grid by a vector point-product judgment method, and connecting all the points to compile a final crack grid.
Preferably, the specific steps of step S2 are as follows:
s21, projecting each edge of the crack surface and the matrix grid to the same coordinate plane to obtain a projection surface;
s22, obtaining each intersection point of the projection of each side of the crack surface and the projection of the matrix grid boundary on the projection surface, and sequentially arranging the intersection points and the first and last points of the projection of each side of the crack surface around the boundary of the crack surface to form a point set;
s23, finding the middle point of each pair of adjacent points in the point set on the projection surface and finding the matrix grid projection where the middle point is located, connecting the matrix grid projections where the middle point is located in sequence to form a closed graph, wherein the closed graph and the matrix grid projection inside the closed graph are the matrix grid projection through which the crack surface on the projection surface passes;
s24, repeating the steps S21-S23, and projecting to x-y, x-y and y-z3 coordinate planes respectively to obtain projection planes of the matrix grid penetrated by the crack plane in 3 coordinate planes.
Preferably, the specific steps of step S3 are as follows:
s31, recording the positions of matrix grid projections penetrated by the crack surfaces in the projection plane of the coordinate plane in the form of column vectors respectively, wherein the set of the projection positions on the X-y coordinate plane is { Xxy，Yxy}，The set of its projection positions on the X-z coordinate plane is { Xxz，ZxzA set of projection positions thereof on the Y-z coordinate plane is { Y }yz，Zyz}；
S32, mixing XxyAnd XxzWherein the same element constitutes a vector Xxy+xzAnd recording the position i of the same element within the respective vectorxyAnd ixzFrom YxyThe selected position is ixyOf (d) constitutes a new vector Yxy+xzFrom Z to ZxzIn the selected position is ixzOf (c) constitutes a new vector Zxy+xzThe above three new vectors are formed into a coordinate set { X }xy+xz，Yxy+xz，Zxy+xz}；
S33, in coordinate set { Xxy+xz，Yxy+xz，Zxy+xzSelect its two columns and Yyz，ZyzThe same coordinate lines constitute a coordinate set { X }x，Yy，Zz}；
S34, in coordinate set { Xx，Yy，ZzAnd (4) screening the matrix grids penetrated by the fracture surfaces, wherein the screening method comprises the following specific steps:
determining the fracture surface equation as:
f(x,y,z)＝ax+by+cz+d＝0，
wherein a, b, c and d are constants related to the geometrical shape of the crack surface,
for coordinate set { Xx，Yy，ZzDetermining the vertex coordinates of a single substrate mesh if there is a pair of adjacent vertices (x) in the substrate mesh1,y1,z1)、(x2,y2,z2) Such that:
f(x1,y1,z1)·f(x2,y2,z2)＜0
the fracture faces cross the matrix lattice, whereas they do not.
Preferably, the specific steps of step S4 are as follows:
s41, determining the intersection point position of the matrix grid penetrated by the crack surface and the crack surface, wherein the specific method comprises the following steps:
s411, determining the equation of the crack surface as follows:
f(x,y,z)＝ax+by+cz+d＝0
wherein a, b, c and d are constants related to the geometrical shape of the crack surface,
s412, each pair of adjacent vertexes (x) of the single matrix mesh through which the fracture surface passes1,y1,z1)、(x2,y2,z2) Substituting into the equation of the fracture surface if it satisfies
f(x1,y1,z1)·f(x2,y2,z2)≤0，
Then, a straight line passing through the pair of adjacent vertexes is obtained, and the straight line equation and the crack surface equation are combined to obtain an intersection point of the straight line equation and the crack surface equation;
s413, repeating the step S412, and obtaining intersection points of the matrix grids and the crack surfaces, through which all the crack surfaces penetrate;
and S42, respectively connecting the intersection points in each matrix grid in the S413 into convex polygons by using a convex hull algorithm to form a primary fracture grid.
Preferably, the specific steps in step S5 are as follows:
s51, determining the end points of all edges on the crack boundary, forming an end point set { X } and determining a central point C in the crack plane;
s52, selecting a point D in the preliminary crack grid, if any two adjacent points A, B in the endpoint set { X } satisfy
If the D point is positioned in the crack boundary, otherwise, the D point is positioned in the primary crack grid and is positioned outside the crack boundary;
s53, if a point in a certain crack grid is positioned outside the crack boundary in the step S52, removing the point, and obtaining a new point on the boundary by taking the crack boundary as constraint to obtain a final crack grid positioned in the boundary;
and S54, repeating the steps S52-S53, removing points outside the crack boundaries in all the initial crack grids, and obtaining the final crack grids.
Has the advantages that: the invention provides an embedded discrete fracture model-based inclined fracture mesh generation method, which is used for generating a generation of an inclined arbitrary polygonal fracture surface in a stratum by using the embedded discrete fracture model and is beneficial to researching the seepage rule in the stratum with complex fracture distribution.
Drawings
FIG. 1 is a diagram of a matrix grid system;
FIG. 2 is a diagram showing the relative positions of the fracture faces and the substrate;
FIG. 3 is a plane of projection of a matrix lattice through which fracture planes pass in the plane of the x-z coordinate axes;
FIG. 4 is a plane of projection of the matrix lattice through which the fracture planes pass in the plane of the x-y coordinate axes;
FIG. 5 is a plane of projection of the matrix lattice through which the fracture planes pass in the plane of the y-z coordinate axes;
FIG. 6 is a matrix grid map with fracture faces passing through;
FIG. 7 is a preliminary fracture grid diagram;
fig. 8 is a final fracture grid.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be described clearly and completely with reference to the accompanying drawings of the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
The invention provides an inclined crack mesh generation method based on an embedded discrete crack model, which comprises the following steps:
the invention provides an inclined crack mesh generation method based on an embedded discrete crack model, which comprises the following steps:
s1, subdividing the matrix mesh according to the basic geometric parameters of the mesh to establish a matrix mesh system;
the reservoir with the size of 100 x 100 establishes a matrix grid system according to the grid size of 10 x 10, as shown in figure 1, wherein the reservoir has a quadrilateral fracture surface with the coordinates of end points of (40, 100, 0); (15, 9, 30); (15, 35, 90); (25, 80, 97.84) and the relative positions of the matrix lattice and the fracture faces are shown in FIG. 2.
S2, projecting the crack surface boundary and the matrix grid system to 3 coordinate axis planes (x-z, y-z, x-y) respectively to obtain projection surfaces of the matrix grid passed by the crack surface in the 3 coordinate planes, as shown in figures 3-5, wherein the shadow part is the matrix grid projection passed by the crack surface on the projection surface, and the specific steps are as follows:
s21, projecting each side of the crack surface and the matrix grid to an x-z coordinate plane to obtain a projection surface, wherein as shown in FIG. 3, the projection of the matrix grid is a square part in the drawing, and the projection of each side of the crack surface is a non-irregular quadrangle in the square graph;
s22, calculating the intersection point of the projection of each side of the crack surface and the projection of the matrix grid boundary on the projection surface, and sequentially arranging the intersection point and the first and last points of the projection of each side of the crack surface around the boundary of the crack surface to form a point set;
s23, finding the middle point of each pair of adjacent points in the point set on the projection surface and finding the matrix grid projection where the middle point is located, as shown in FIG. 3, the slash shadow part in the graph is the matrix grid projection where the middle point is located, connecting the matrix grid projections where the middle point is located in sequence to form a closed graph, as shown in the slash shadow part in FIG. 3, wherein the closed graph (the slash shadow part in FIG. 3) and the matrix grid projection inside the closed graph (the cross shadow part in FIG. 3) are the matrix grid projection where the crack surface passes through on the x-z coordinate plane;
s24, repeating the steps S21-S23, and projecting to the x-y coordinate plane and the y-z coordinate plane respectively to obtain projection planes of the matrix grid penetrated by the crack plane in the x-y coordinate plane and the y-z coordinate plane, which are respectively shown in the figure 4 and the figure 5.
S3, in each projection plane, inverting the projection of the matrix mesh passed through by the fracture plane to obtain the matrix mesh passed through by the fracture plane, as shown in fig. 6, the specific steps are as follows:
s31, recording the positions of matrix grid projections penetrated by the crack surfaces in the projection plane of the coordinate plane in the form of column vectors respectively, wherein the set of the projection positions on the X-y coordinate plane is { Xxy，YxyA set of projection positions thereof on an X-z coordinate plane is { X }xz，ZxzA set of projection positions thereof on the Y-z coordinate plane is { Y }yz，Zyz}；
S32, mixing XxyAnd XxzWherein the same element constitutes a vector Xxy+xzAnd recording the position i of the same element within the respective vectorxyAnd ixzFrom YxyThe selected position is ixyOf (d) constitutes a new vector Yxy+xzFrom Z to ZxzIn the selected position is ixzOf (c) constitutes a new vector Zxy+xzThe above three new vectors are formed into a coordinate set { X }xy+xz，Yxy+xz，Zxy+xz}；
S33, in coordinate set { Xxy+xz，Yxy+xz，Zxy+xzSelect its two columns and Yyz，ZyzThe same coordinate lines constitute a coordinate set { X }x，Yy，Zz}；
S34, in coordinate set { Xx，Yy，ZzAnd (4) screening the matrix grids penetrated by the fracture surfaces, wherein the screening method comprises the following specific steps:
determining the fracture surface equation as:
f(x,y,z)＝ax+by+cz+d＝0，
wherein a, b, c and d are constants related to the geometrical shape of the crack surface,
for coordinate set { Xx，Yy，ZzDetermining the vertex coordinates of a single substrate mesh if there is a pair of adjacent vertices (x) in the substrate mesh1,y1,z1)、(x2,y2,z2) Such that:
f(x1,y1,z1)·f(x2,y2,z2)＜0
the fracture faces cross the matrix lattice, whereas they do not.
S4, splitting a fracture plane by using the boundary of the matrix mesh penetrated by the fracture surface to obtain a primary fracture mesh, as shown in FIG. 7, the method comprises the following specific steps:
s41, determining the intersection point position of the matrix grid penetrated by the crack surface and the crack surface, wherein the specific method comprises the following steps:
s411, determining the equation of the crack surface as follows:
f(x,y,z)＝ax+by+cz+d＝0
wherein a, b, c and d are constants related to the geometrical shape of the crack surface,
s412, each pair of adjacent vertexes (x) of the single matrix mesh through which the fracture surface passes1,y1,z1)、(x2,y2,z2) Substituting into the equation of the fracture surface if it satisfies
f(x1,y1,z1)·f(x2,y2,z2)≤0，
Then, a straight line passing through the pair of adjacent vertexes is obtained, and the straight line equation and the crack surface equation are combined to obtain an intersection point of the straight line equation and the crack surface equation;
s413, repeating the step S412, and obtaining intersection points of the matrix grids and the crack surfaces, through which all the crack surfaces penetrate;
and S42, respectively connecting the intersection points in each matrix grid in the S413 into convex polygons by using a convex hull algorithm to form a primary fracture grid.
S5, removing points outside the fracture boundary in the preliminary fracture grid with the fracture boundary as a constraint condition to obtain a final fracture grid, as shown in fig. 8, the specific steps are as follows:
s51, determining the end points of all edges on the crack boundary, forming an end point set { X } and determining a central point C in the crack plane;
s52, selecting a point D in the preliminary crack grid, if any two adjacent points A, B in the endpoint set { X } satisfy
If the D point is positioned in the crack boundary, otherwise, the D point is positioned in the primary crack grid and is positioned outside the crack boundary;
s53, if a point in a certain crack grid is positioned outside the crack boundary in the step S52, removing the point, and obtaining a new point on the boundary by taking the crack boundary as constraint to obtain a final crack grid positioned in the boundary;
and S54, repeating the steps S52-S53, removing points outside the crack boundaries in all the initial crack grids, and obtaining the final crack grids.
Although the present invention has been described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from the scope of the present invention.
Claims (5)
1. An inclined fracture mesh generation method based on an embedded discrete fracture model is characterized by comprising the following steps:
s1, subdividing the matrix mesh according to the basic geometric parameters of the mesh to establish a matrix mesh system;
s2, projecting the crack surface boundary and the matrix grid system to 3 coordinate axis planes respectively to obtain projection surfaces of the matrix grid penetrated by the crack surface in the 3 coordinate axis planes;
s3, inverting the projection of the matrix grid penetrated by the crack surface in each projection plane to obtain the matrix grid penetrated by the crack surface;
s4, dividing a fracture plane by using the boundary of the matrix grid penetrated by the fracture surface to obtain a primary fracture grid;
and S5, taking the central point in the fracture surface, taking the fracture boundary as a constraint condition, and removing the points outside the fracture boundary in the preliminary fracture grid by a vector point-product judgment method to obtain the final fracture grid in the boundary.
2. The oblique fracture mesh splitting method based on the embedded discrete fracture model as claimed in claim 1, wherein the specific steps of step S2 are as follows:
s21, projecting each edge of the crack surface and the matrix grid to the same coordinate plane to obtain a projection surface;
s22, obtaining intersection points of the projection of each side of the crack surface and the projection of the matrix grid boundary on the projection surface, and sequentially arranging the intersection points and the first and last points of the projection of each side of the crack surface according to the adjacent sequence to form a point set;
s23, finding the middle point of each pair of adjacent points in the point set on the projection surface and finding the matrix grid projection where the middle point is located, connecting the matrix grid projections where the middle point is located in sequence to form a closed graph, wherein the closed graph and the matrix grid projection inside the closed graph are the matrix grid projection through which the crack surface on the projection surface passes;
s24, repeating the steps S21-S23, and projecting to 3 coordinate planes x-y, x-y and y-z respectively to obtain projection planes of the matrix grid penetrated by the crack plane in the 3 coordinate planes.
3. The oblique fracture mesh splitting method based on the embedded discrete fracture model as claimed in claim 1, wherein the specific steps of step S3 are as follows:
s31, recording the positions of matrix grid projections penetrated by the crack surfaces in the projection plane of the coordinate plane in the form of column vectors respectively, wherein the set of the projection positions on the X-y coordinate plane is { Xxy，YxyA set of projection positions thereof on an X-z coordinate plane is { X }xz，ZxzA set of projection positions thereof on the Y-z coordinate plane is { Y }yz，Zyz}；
S32, mixing XxyAnd XxzWherein the same element constitutes a vector Xxy+xzAnd recording the position i of the same element within the respective vectorxyAnd ixzFrom YxyThe selected position is ixyOf (d) constitutes a new vector Yxy+xzFrom Z to ZxzIn the selected position is ixzOf (c) constitutes a new vector Zxy+xzThe above three new vectors are formed into a coordinate set { X }xy+xz，Yxy+xz，Zxy+xz}；
S33, in coordinate set { Xxy+xz，Yxy+xz，Zxy+xzSelect its two columns and Yyz，ZyzThe same coordinate lines constitute a coordinate set { X }x，Yy，Zz}；
S34, in coordinate set { Xx，Yy，ZzAnd (4) screening the matrix grids penetrated by the fracture surfaces, wherein the screening method comprises the following specific steps:
determining the fracture surface equation as:
f(x,y,z)＝ax+by+cz+d＝0，
wherein a, b, c and d are constants related to the geometrical shape of the crack surface,
for coordinate set { Xx，Yy，ZzDetermining the vertex coordinates of a single substrate mesh if there is a pair of adjacent vertices (x) in the substrate mesh1,y1,z1)、(x2,y2,z2) Such that:
f(x1,y1,z1)·f(x2,y2,z2)＜0
the fracture faces cross the matrix lattice, whereas they do not.
4. The oblique fracture mesh splitting method based on the embedded discrete fracture model as claimed in claim 1, wherein the specific steps of step S4 are as follows:
s41, determining the intersection point position of the matrix grid penetrated by the crack surface and the crack surface, wherein the specific method comprises the following steps:
s411, determining the equation of the crack surface as follows:
f(x,y,z)＝ax+by+cz+d＝0
wherein a, b, c and d are constants about the geometrical form of the crack surface;
s412, each pair of adjacent vertexes (x) of the single matrix mesh through which the fracture surface passes1,y1,z1)、(x2,y2,z2) Substituting into the equation of the fracture surface if it satisfies
f(x1,y1,z1)·f(x2,y2,z2)≤0，
Then, a straight line passing through the pair of adjacent vertexes is obtained, and the straight line equation and the crack surface equation are combined to obtain an intersection point of the straight line equation and the crack surface equation;
s413, repeating the step S412, and obtaining intersection points of the matrix grids and the crack surfaces, through which all the crack surfaces penetrate;
and S42, respectively connecting the intersection points in each matrix grid in the S413 into convex polygons by using a convex hull algorithm to form a primary fracture grid.
5. The oblique fracture mesh splitting method based on the embedded discrete fracture model as claimed in claim 1, wherein the specific steps of step S5 are as follows: :
s51, determining the end points of all edges on the crack boundary, forming an end point set { X } and determining a central point C in the crack plane;
s52, selecting a point D in the preliminary crack grid, if any two adjacent points A, B in the endpoint set { X } satisfy
If the D point is positioned in the crack boundary, otherwise, the D point is positioned in the primary crack grid and is positioned outside the crack boundary;
s53, if a point in a certain crack grid is positioned outside the crack boundary in the step S52, removing the point, and obtaining a new point on the boundary by taking the crack boundary as constraint to obtain a final crack grid positioned in the boundary;
and S54, repeating the steps S52-S53, removing points outside the crack boundaries in all the initial crack grids, and obtaining the final crack grids.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911328526.1A CN110992488B (en) | 2019-12-20 | 2019-12-20 | Inclined crack grid dividing method based on embedded discrete crack model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911328526.1A CN110992488B (en) | 2019-12-20 | 2019-12-20 | Inclined crack grid dividing method based on embedded discrete crack model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110992488A true CN110992488A (en) | 2020-04-10 |
| CN110992488B CN110992488B (en) | 2020-09-18 |
Family
ID=70073679
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911328526.1A Active CN110992488B (en) | 2019-12-20 | 2019-12-20 | Inclined crack grid dividing method based on embedded discrete crack model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110992488B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8695476B2 (en) * | 2011-03-14 | 2014-04-15 | The United States Of America, As Represented By The Secretary Of The Navy | Armor plate with shock wave absorbing properties |
| CN105653746A (en) * | 2014-11-10 | 2016-06-08 | 中国石油化工股份有限公司 | Embedded discrete fracture model-based fractured well modeling and simulation method |
| CN108595782A (en) * | 2018-04-03 | 2018-09-28 | 中国石油大学(北京) | A kind of discrete fractures mesostroma and interstitial mass transfercomputation method |
-
2019
- 2019-12-20 CN CN201911328526.1A patent/CN110992488B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8695476B2 (en) * | 2011-03-14 | 2014-04-15 | The United States Of America, As Represented By The Secretary Of The Navy | Armor plate with shock wave absorbing properties |
| CN105653746A (en) * | 2014-11-10 | 2016-06-08 | 中国石油化工股份有限公司 | Embedded discrete fracture model-based fractured well modeling and simulation method |
| CN108595782A (en) * | 2018-04-03 | 2018-09-28 | 中国石油大学(北京) | A kind of discrete fractures mesostroma and interstitial mass transfercomputation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110992488B (en) | 2020-09-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ito et al. | Octree‐based reasonable‐quality hexahedral mesh generation using a new set of refinement templates | |
| US7800627B2 (en) | Mesh quilting for geometric texture synthesis | |
| Lo | Finite element mesh generation and adaptive meshing | |
| CN110543716A (en) | three-dimensional overhead cable level power grid optimization method and device and computer equipment | |
| CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
| CN113593030A (en) | Geological profile generation method, system, terminal and medium based on three-dimensional model | |
| CN104851126A (en) | Three-dimensional model decomposition method and three-dimensional model decomposition device based on generalized cylinder | |
| CN110334365B (en) | Numerical simulation method and system for flow of reservoir after heterogeneous fracturing | |
| CN104504754A (en) | Multipoint statistic modeling method and device | |
| CN109461209B (en) | Novel structural grid generation method | |
| CN106384386A (en) | Grid processing method for LOD model generation and grid processing system thereof and 3D reconstruction method and system | |
| CN105719342B (en) | A kind of the three-dimensional modeling method for visualizing and device of ground fissure geologic body | |
| CN110826128A (en) | Design method for rapid forming of dredging and trenching of bottom surface with any shape | |
| CN108595782B (en) | Calculation method for mass transfer between matrix and cracks in discrete cracks | |
| CN110992488B (en) | Inclined crack grid dividing method based on embedded discrete crack model | |
| CN115937467B (en) | Method and system for dividing random three-dimensional fracture in upscale model grid | |
| CN115908733A (en) | Real-time segmentation and three-dimensional visualization method for angular point grid geological model | |
| CN116664758A (en) | Tetrahedral grid self-adaptive hierarchical refinement method based on fault plane | |
| Remacle et al. | Fast and robust mesh generation on the sphere–application to coastal domains | |
| CN107767451A (en) | A kind of 3DEC complicated landform modeling methods based on GOCAD Raster Data Models | |
| Tian | 3D modeling and digital preservation of ancient architectures based on autoCAD and 3Dmax | |
| Zhang et al. | Conformal adaptive hexahedral-dominant mesh generation for CFD simulation in architectural design applications | |
| CN112507551A (en) | Unstructured dynamic mesh generation method and device | |
| Li et al. | Generalized prism grid: a pillar-based unstructured grid for simulation of reservoirs with complicated geological geometries | |
| CN111951394A (en) | Fault structure unit three-dimensional model construction method and device based on geological map |
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 |