CN110968930B - Geological variable attribute interpolation method and system - Google Patents
Geological variable attribute interpolation method and system Download PDFInfo
- Publication number
- CN110968930B CN110968930B CN201811159460.3A CN201811159460A CN110968930B CN 110968930 B CN110968930 B CN 110968930B CN 201811159460 A CN201811159460 A CN 201811159460A CN 110968930 B CN110968930 B CN 110968930B
- Authority
- CN
- China
- Prior art keywords
- line segments
- mesh
- points
- new
- tetrahedral
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
Landscapes
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Generation (AREA)
Abstract
The invention discloses a geological variable attribute interpolation method and a geological variable attribute interpolation system, wherein the method comprises the following steps: loading the explained geological body result data in the set work area range; constructing a surface grid model of the geologic body according to the control point data of the geologic body; taking the triangular grid data of the block geometric surface as a limiting condition, and carrying out regularization processing on the limiting condition to obtain regular limiting points, limiting line segments and limiting surfaces and obtain a regularized surface grid model; adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency; and performing attribute interpolation on the tetrahedral mesh by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result. The variable attribute interpolation method and system for the geologic body realize variable attribute interpolation calculation of the complex geologic body, have high modeling efficiency and strong practicability, achieve the aim of rapid attribute modeling, and are suitable for parallel calculation of large data body attribute models.
Description
Technical Field
The invention belongs to the field of seismic forward attribute modeling, and particularly relates to a geological variable attribute interpolation method and system.
Background
The seismic forward modeling technology is widely applied to acquisition, processing and interpretation of seismic exploration, and plays an important role in design optimization of an observation system, processing parameter extraction and verification of an interpretation scheme. However, the quality of the result of the seismic forward modeling depends on multiple aspects such as wavelet selection, observation system definition, attribute modeling, a forward modeling numerical solution and the like. The establishment of reliable attribute models (including compressional wave velocity, shear wave velocity, density, saturation, porosity and the like) is a crucial and inevitable link in seismic forward simulation and numerical analysis. At present, most of adopted attribute modeling software is realized on the basis of layers, but only one constant can be used for filling the layers, and the situation that the field real complex geological structure is completely not met. Therefore, a modeling method capable of adapting to the field real complex geological structure condition is particularly needed.
Disclosure of Invention
The invention aims to provide a geologic body variable attribute interpolation method and system which have higher modeling efficiency and are suitable for complex geologic structure conditions.
In order to achieve the above object, the present invention provides a geological variable property interpolation method, which includes:
step 1: loading the explained geological body result data in the set work area range;
step 2: constructing a surface grid model of the geologic body according to the control point data of the geologic body;
and step 3: taking the triangular grid data of the block geometric surface as a limiting condition, and carrying out regularization processing on the limiting condition to obtain regular limiting points, limiting line segments and limiting surfaces and obtain a regularized surface grid model;
and 4, step 4: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency;
and 5: and (5) performing attribute interpolation on the tetrahedral mesh obtained in the step (4) by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result.
Preferably, the surface mesh model of the complex geologic body is constructed by using a geometric modeling method, wherein the geometric modeling method comprises the following steps: two-dimensional limited Delaunay subdivision, surface intersection, surface stitching and complex geologic body rapid tracking method.
Preferably, step 3 specifically comprises:
step 301: setting a constraint line set and a discrete point set according to the triangular grid data;
step 302: recording the end points of all the line segments contained in the constraint line segment set into the discrete point set;
step 303: selecting any two line segments in the constraint line segment set, judging whether the two line segments are intersected or not, if the two line segments are intersected and collinear, recording collinear parts into the constraint line segment set, deleting the two line segments in the constraint line segment set, and if the two line segments are only intersected and not collinear, recording intersection points into the discrete point set;
step 304: if the triangular mesh does not meet the regularization requirement, adding a new point, judging whether the new point belongs to the discrete point set, and if not, recording the new point into the discrete point set; judging whether the line segments in the constraint line segment set contain the new points or not, if so, dividing the line segments containing the new points into two new line segments by taking the new points as end points, recording the new line segments into the constraint line segment set, and deleting the line segments containing the new points;
step 305: and (6) repeating the steps 303-304 until all the triangular meshes meet the regularization requirement, obtaining regular limiting points, limiting line segments and limiting surfaces, and obtaining the regularization surface mesh model consisting of the regular limiting points, the regular line segments and the limiting surfaces.
Preferably, whether the regularization requirement is satisfied is determined based on the area of the triangle and the minimum angle of the triangle.
Preferably, step 4 specifically includes:
step 401: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain a tetrahedral mesh;
step 402: checking whether each point, each edge and each triangle in the triangle mesh data of the limited condition are in the subdivision result or not, and if not, respectively adding the intersection points into the current triangle mesh and the tetrahedral mesh;
step 403: and optimizing the tetrahedral mesh by adopting a variable mesh method to obtain the tetrahedral mesh with topological consistency.
Preferably, the mesh changing method is to add auxiliary points according to the geologic body construction form to construct a tetrahedral mesh.
Preferably, attribute interpolation is performed on the complex geologic body by using a gradient-based variable attribute interpolation algorithm, wherein a formula used by the variable attribute interpolation algorithm is as follows:
SampleValue[i]＝mBaseValue+mGX*gx+mGY*gy+mGZ*gz (1)
wherein, sampleValue [ i ] represents the attribute interpolation result of the ith tetrahedron, mBaseValue represents the basic attribute value of the block where the tetrahedron to be inserted is located, mGX, mGY and mGZ respectively represent the attribute change gradient in the X, Y, Z coordinate direction, and gx, gy and gz respectively represent the coordinate increment in the X, Y, Z coordinate direction.
The invention also provides a geological variable attribute interpolation system, which comprises: a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
step 1: loading the explained geologic body result data in a set work area range;
step 2: constructing a surface grid model of the geologic body according to the control point data of the geologic body;
and step 3: taking the triangular grid data of the block geometric surface as a limiting condition, and carrying out regularization processing on the limiting condition to obtain regular limiting points, limiting line segments and limiting surfaces and obtain a regularized surface grid model;
and 4, step 4: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency;
and 5: and (4) performing attribute interpolation on the tetrahedral mesh obtained in the step (4) by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result.
Preferably, step 3 specifically comprises:
step 301: setting a constraint line set and a discrete point set according to the triangular grid data;
step 302: recording the end points of all line segments contained in the constraint line segment set into the discrete point set;
step 303: selecting any two line segments in the constraint line segment set, judging whether the two line segments are intersected or not, if the two line segments are intersected and collinear, recording collinear parts into the constraint line segment set, deleting the two line segments in the constraint line segment set, and if the two line segments are only intersected and not collinear, recording intersection points into the discrete point set;
step 304: if the triangular mesh does not meet the regularization requirement, adding a new point, judging whether the new point belongs to the discrete point set, and if not, recording the new point into the discrete point set; judging whether the line segments in the constraint line segment set contain the new points or not, if so, dividing the line segments containing the new points into two new line segments by taking the new points as end points, recording the new line segments into the constraint line segment set, and deleting the line segments containing the new points;
step 305: and repeating the steps 303-304 until all the triangular meshes meet the regularization requirement, obtaining regular limit points, limit line segments and limit surfaces, and obtaining a regularization surface mesh model consisting of the regular limit points, the limit line segments and the limit surfaces.
Preferably, step 4 specifically includes:
step 401: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain a tetrahedral mesh;
step 402: checking whether each point, each edge and each triangle in the triangle grid data of the limited condition are in the subdivision result, if not, respectively adding the intersection points into the current triangle grid and the tetrahedral grid;
step 403: and optimizing the tetrahedral mesh by adopting a mesh changing method to obtain the tetrahedral mesh with topological consistency.
The invention has the beneficial effects that: according to the variable attribute interpolation method for the geological body, the triangular grid data of the geometric surface of the block body is used as a limiting condition, the limiting condition is subjected to regularization processing to obtain a regularized surface grid model, a three-dimensional Delaunay subdivision algorithm is adopted to subdivide the regularized surface grid model to obtain a tetrahedral grid with topological consistency, and attribute interpolation is performed on the tetrahedral grid, so that variable attribute interpolation calculation of the complex geological body is realized, the modeling efficiency is high, the practicability is high, the purpose of rapid attribute modeling is achieved, and the method is suitable for parallel calculation of the attribute model of the large data body.
The present invention has other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts throughout.
FIG. 1 shows a flowchart of a method for interpolating a geological variable attribute according to an embodiment of the invention.
FIG. 2 illustrates loading a complex geological control point map based on a work area according to one embodiment of the invention.
FIG. 3 illustrates a regularized surface mesh model of a method for geologic variable property interpolation according to one embodiment of the present invention.
FIG. 4 shows a tetrahedral mesh with topological consistency for a method of geologic variable property interpolation according to one embodiment of the present invention.
FIG. 5 is a diagram illustrating variable attribute interpolation results of a geological variable attribute interpolation method according to an embodiment of the present invention.
Detailed Description
Preferred embodiments of the present invention will be described in more detail below. While the following describes preferred embodiments of the present invention, it should be understood that the present invention may be embodied in various forms and should not be limited by the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
The geologic body variable attribute interpolation method comprises the following steps:
step 1: loading the explained geologic body result data in a set work area range;
and 2, step: constructing a surface grid model of the geologic body according to the control point data of the geologic body;
and 3, step 3: taking the triangular grid data of the block geometric surface as a limiting condition, and carrying out regularization processing on the limiting condition to obtain regular limiting points, limiting line segments and limiting surfaces and obtain a regularized surface grid model;
and 4, step 4: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency;
and 5: and (4) performing attribute interpolation on the surface volume grid obtained in the step (4) by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result.
Specifically, work area range control point coordinates (0,0,0), (16000, 32000,0), (16000, 32000 and 16000) are input to form a work area contour line, complex geologic body result data which are already explained by geologists in the work area range are loaded, the complex geologic body result data comprise horizon data, fault data, sand body data and the like, and basic data are provided for the geometric modeling of the complex geologic body. According to control point data of the complex geologic body, a surface mesh model of the complex geologic body is constructed by using a geometric modeling method, block geometric surface triangular mesh data obtained by geometric modeling is used as a limiting condition, the limiting condition is subjected to regularization processing to obtain regular limiting points, limiting line segments and limiting surfaces, a regularized surface mesh model is obtained, a three-dimensional limiting DELAUNAY subdivision algorithm is adopted to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency, and a gradient-based variable attribute interpolation algorithm is used to perform attribute interpolation on the tetrahedral meshes with topological consistency.
According to an exemplary implementation mode, the variable attribute interpolation method of the geological body is characterized in that triangular grid data of a block geometric surface are used as limiting conditions, the limiting conditions are subjected to regularization processing to obtain a regularized surface grid model, a three-dimensional Delaunay subdivision algorithm is adopted to subdivide the regularized surface grid model to obtain tetrahedral grids with topological consistency, and attribute interpolation is carried out on the tetrahedral grids, so that variable attribute interpolation calculation of the complex geological body is realized, the modeling efficiency is high, the practicability is high, the purpose of rapid attribute modeling is achieved, and the method is suitable for parallel calculation of large data body attribute models.
Preferably, the surface mesh model of the complex geologic body is constructed by using a geometric modeling method, wherein the geometric modeling method comprises the following steps: two-dimensional limited Delaunay subdivision, surface intersection, surface stitching and complex geologic body rapid tracking method.
As a preferred scheme, step 3 specifically comprises:
step 301: setting a constraint line set and a discrete point set according to the triangular grid data;
step 302: recording the end points of all line segments contained in the constraint line segment set into the discrete point set;
step 303: selecting any two line segments in the constraint line segment set, judging whether the two line segments are intersected or not, if the two line segments are intersected and collinear, recording collinear parts into the constraint line segment set, deleting the two line segments in the constraint line segment set, and if the two line segments are only intersected and not collinear, recording intersection points into a discrete point set;
step 304: if the triangular mesh does not meet the regularization requirement, adding a new point, judging whether the new point belongs to the discrete point set, and if not, recording the new point into the discrete point set; judging whether the line segments in the constraint line segment set contain new points or not, if so, dividing the line segments containing the new points into two new line segments by taking the new points as end points, recording the new line segments into the constraint line segment set, and deleting the line segments containing the new points;
step 305: and repeating the steps 303-304 until all the triangular meshes meet the regularization requirement, obtaining regular limit points, limit line segments and limit surfaces, and obtaining a regularization surface mesh model consisting of the regular limit points, the limit line segments and the limit surfaces.
Specifically, all line segments of triangular grid data of the surface grid model are extracted, all line segments exist together and are set as a constraint line segment set, all points of the triangular grid data of the surface grid model are extracted, all point segments exist together and are set as a discrete point set, the discrete point set comprises end points of all line segments in the constraint line segment set, the constraint line segment set and the discrete point set are processed according to the regularization requirement, any two line segments in the constraint line segment set are selected, if the two line segments are collinear, a collinear part is recorded into the constraint line segment set, the two line segments are deleted from the constraint line segment set, and if the two line segments are only intersected and not collinear, an intersection point is recorded into the discrete point set; if the triangular meshes do not meet the regularization requirement, adding a new point, determining the position of the new point according to the actual situation, judging whether the new point belongs to a discrete point set, if not, recording the new point into the discrete point set, then judging whether the line segments in the constraint line segment set contain the new point, if so, taking the new point as an end point, dividing the line segment containing the new point into two new line segments, recording the new line segments into the constraint line segment set, deleting the line segments containing the new point in the constraint line segment set, regularizing all the line segments in the constraint line segment set and the discrete point set, including the added new line segments and the new point, until all the triangular meshes meet the regularization requirement, obtaining regular limit points, limit line segments and limit surfaces, and obtaining a regularization surface mesh model consisting of the regular limit points, limit line segments and limit surfaces.
Preferably, whether the regularization requirement is satisfied is determined based on the area of the triangle and the minimum angle of the triangle.
As a preferred scheme, step 4 specifically comprises:
step 401: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes;
step 402: checking whether each point, each edge and each triangle in the triangle mesh data of the limited condition are in the subdivision result, if not, respectively adding the intersection points into the current triangle mesh and the tetrahedral mesh;
step 403: and optimizing the tetrahedral mesh by adopting a variable mesh method to obtain the tetrahedral mesh with topological consistency.
Specifically, a general tetrahedron subdivision is based on the geometric configuration of a point set, which can only ensure the existence of points in a Delaunay tetrahedron, but cannot ensure the existence of limit points, limit edges and limit surfaces in the Delaunay tetrahedron, so a regularized surface mesh model is subdivided by using a three-dimensional limit Delaunay subdivision algorithm to obtain tetrahedral meshes, then the existence of each edge in a triangle mesh under a limited condition is checked, the intersection points of the edges and the tetrahedrons which do not exist in the subdivision result are respectively added into the current triangle mesh and the tetrahedral mesh by using a Delaunay void method, then the existence of each triangle in the triangle mesh under the limited condition is checked, the intersection points of the triangles and the tetrahedrons which do not exist in the subdivision result are respectively added into the current triangle mesh and the tetrahedral mesh by using a Delaunay void method, the tetrahedral mesh is optimized by using a variable mesh method to obtain tetrahedral meshes with topological consistency, and finally the consistency of the limit conditions and the topological result on the geometry and topology is checked, and the tetrahedral meshes with the consistency on the limit conditions and the topological consistency of the subdivision result, and the tetrahedral meshes are obtained after the constraint conditions and the subdivision result reach the geometric consistency on the topological consistency of the topology.
Preferably, the variable mesh method is to add auxiliary points according to the geologic body structure form to construct a tetrahedral mesh.
Specifically, the geologic body complex structure is finely described by adding more auxiliary points at a place with a complex structure, and the geologic structure is described by adding a small number of auxiliary points at a place with a simpler structure, so that the grid quality is ensured, and the existence of the defined points, the defined edges and the defined surfaces in a Delaunay tetrahedron is ensured, thereby providing an accurate tetrahedral grid with topological consistency for the next step of complex geologic body variable attribute interpolation.
As a preferred scheme, a gradient-based variable attribute interpolation algorithm is adopted to perform attribute interpolation on the complex geologic body, and the formula used by the variable attribute interpolation algorithm is as follows:
SampleValue[i]＝mBaseValue+mGX*gx+mGY*gy+mGZ*gz (1)
wherein, sampleValue [ i ] represents the attribute interpolation result of the ith tetrahedron, mBaseValue represents the basic attribute value of the block body where the tetrahedron to be inserted is located, mGX, mGY and mGZ respectively represent the attribute change gradient in the X, Y, Z coordinate direction, and gx, gy and gz respectively represent the coordinate increment in the X, Y, Z coordinate direction.
The geologic variable attribute interpolation system according to the present invention comprises: a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
step 1: loading the explained geologic body result data in a set work area range;
step 2: constructing a surface grid model of the geologic body according to the control point data of the geologic body;
and step 3: taking the triangular grid data of the block geometric surface as a limiting condition, and carrying out regularization processing on the limiting condition to obtain regular limiting points, limiting line segments and limiting surfaces and obtain a regularized surface grid model;
and 4, step 4: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency;
and 5: and (4) performing attribute interpolation on the surface volume grid obtained in the step (4) by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result.
Specifically, work area range control point coordinates (0,0,0), (16000, 32000,0), (16000, 32000 and 16000) are input to form a work area contour line, complex geologic body result data which are already explained by geologists in the work area range are loaded, the complex geologic body result data comprise horizon data, fault data, sand body data and the like, and basic data are provided for the geometric modeling of the complex geologic body. According to control point data of the complex geologic body, a geometric modeling method is utilized to construct a surface mesh model of the complex geologic body, block geometric surface triangular mesh data obtained through geometric modeling is used as a limiting condition, the limiting condition is subjected to regularization processing, regular limiting points, limiting line segments and limiting surfaces are obtained, a regularized surface mesh model is obtained, a three-dimensional limiting DELAUNAY subdivision algorithm is adopted to subdivide the regularized surface mesh model, tetrahedral meshes with topological consistency are obtained, and a gradient-based variable attribute interpolation algorithm is utilized to carry out attribute interpolation on the tetrahedral meshes with topological consistency.
According to an exemplary implementation mode, the variable attribute interpolation system of the geological body takes triangular grid data of a block geometric surface as a limiting condition, the limiting condition is subjected to regularization processing to obtain a regularized surface grid model, a three-dimensional Delaunay subdivision algorithm is adopted to subdivide the regularized surface grid model to obtain a tetrahedral grid with topological consistency, and attribute interpolation is carried out on the tetrahedral grid, so that variable attribute interpolation calculation of a complex geological body is realized, the modeling efficiency is high, the practicability is high, the purpose of rapid attribute modeling is achieved, and the variable attribute interpolation system is suitable for parallel calculation of a large data body attribute model.
Preferably, the surface mesh model of the complex geologic body is constructed by using a geometric modeling method, wherein the geometric modeling method comprises the following steps: two-dimensional limited Delaunay subdivision, surface intersection, surface stitching and complex geologic body rapid tracking method.
As a preferred scheme, step 3 specifically comprises:
step 301: setting a constraint line set and a discrete point set according to the triangular grid data;
step 302: recording the end points of all the line segments contained in the constraint line segment set into the discrete point set;
step 303: selecting any two line segments in the constraint line segment set, judging whether the two line segments are intersected or not, if the two line segments are intersected and collinear, recording the collinear part into the constraint line segment set, deleting the two line segments in the constraint line segment set, and if the two line segments are only intersected and not collinear, recording the intersection point into the discrete point set;
step 304: if the triangular mesh does not meet the regularization requirement, adding a new point, judging whether the new point belongs to the discrete point set, and if not, recording the new point into the discrete point set; judging whether the line segments in the constraint line segment set contain new points or not, if so, dividing the line segments containing the new points into two new line segments by taking the new points as end points, recording the new line segments into the constraint line segment set, and deleting the line segments containing the new points;
step 305: and repeating the steps 303-304 until all the triangular meshes meet the regularization requirement, obtaining regular limit points, limit line segments and limit surfaces, and obtaining a regularization surface mesh model consisting of the regular limit points, the limit line segments and the limit surfaces.
Specifically, all line segments of triangular grid data of the surface grid model are extracted, all line segments exist together and are set as a constraint line segment set, all points of the triangular grid data of the surface grid model are extracted, all point segments exist together and are set as a discrete point set, the discrete point set comprises end points of all line segments in the constraint line segment set, the constraint line segment set and the discrete point set are processed according to the regularization requirement, any two line segments in the constraint line segment set are selected, if the two line segments are collinear, the collinear part is recorded into the constraint line segment set, the two line segments are deleted from the constraint line segment set, and if the two line segments are only intersected and not collinear, the intersection point is recorded into the discrete point set; if the triangular mesh does not meet the regularization requirement, adding a new point, determining the position of the new point according to the actual situation, judging whether the new point belongs to a discrete point set, if not, recording the new point into the discrete point set, then judging whether the line segments in the constraint line segment set contain the new point, if so, taking the new point as an end point, dividing the line segment containing the new point into two new line segments, recording the new line segments into the constraint line segment set, deleting the line segments containing the new point in the constraint line segment set, processing all the line segments and points in the constraint line segment set and the discrete point set according to the regularization requirement, including the added new line segments and new points, until all the triangular meshes meet the regularization requirement, obtaining regular limit points, limit line segments and limit surfaces, and obtaining a regularization surface mesh model consisting of the regular limit points, limit line segments and limit surfaces.
Preferably, whether the regularization requirement is satisfied is determined based on the area of the triangle and the minimum angle of the triangle.
As a preferred scheme, step 4 specifically comprises:
step 401: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes;
step 402: checking whether each point, each edge and each triangle in the triangle mesh data of the limited condition are in the subdivision result or not, and if not, respectively adding the intersection points into the current triangle mesh and the tetrahedral mesh;
step 403: and optimizing the tetrahedral mesh by adopting a variable mesh method to obtain the tetrahedral mesh with topological consistency.
Specifically, a general tetrahedron subdivision is based on the geometric configuration of a point set, which can only ensure the existence of points in a Delaunay tetrahedron, but cannot ensure the existence of limit points, limit edges and limit surfaces in the Delaunay tetrahedron, so that a regular surface mesh model is subdivided by adopting a three-dimensional limit Delaunay subdivision algorithm to obtain tetrahedral meshes, the interior of the surface mesh model is filled, then the existence of each edge in a condition-limited triangular mesh is checked, the intersection points of the edges and the tetrahedrons which do not exist in the subdivision result are respectively added into the current triangular mesh and the tetrahedral mesh by applying a Delaunay void method, then the existence of each triangle in the condition-limited triangular mesh is checked, the intersection points of the triangles and the tetrahedrons which do not exist in the subdivision result are respectively added into the current triangular mesh and the tetrahedral mesh by applying a Delaunay void method, the tetrahedral mesh is optimized by adopting a variable mesh method, the mesh model is compacted by subdivision and optimization, the interior is filled, the tetrahedral mesh with consistency is obtained, and finally the topology consistency of the constraint conditions and the tetrahedral mesh is obtained after the subdivision result is consistent with the topology and the topology on the limitation of the tetrahedral mesh.
Preferably, the variable mesh method is to add auxiliary points according to the geologic body structure form to construct a tetrahedral mesh.
Specifically, the geologic body complex structure is finely described by adding more auxiliary points at a place with a complex structure, and the geologic structure is described by adding a small number of auxiliary points at a place with a simpler structure, so that the grid quality is ensured, and the existence of the defined points, the defined edges and the defined surfaces in a Delaunay tetrahedron is ensured, thereby providing a precise tetrahedron grid with topological consistency for the next step of complex geologic body variable attribute interpolation.
As a preferred scheme, a gradient-based variable attribute interpolation algorithm is adopted to perform attribute interpolation on the complex geologic body, and the formula used by the variable attribute interpolation algorithm is as follows:
SampleValue[i]＝mBaseValue+mGX*gx+mGY*gy+mGZ*gz (1)
wherein, sampleValue [ i ] represents the attribute interpolation result of the ith tetrahedron, mBaseValue represents the basic attribute value of the block body where the tetrahedron to be inserted is located, mGX, mGY and mGZ respectively represent the attribute change gradient in the X, Y, Z coordinate direction, and gx, gy and gz respectively represent the coordinate increment in the X, Y, Z coordinate direction.
Examples
FIG. 1 shows a flowchart of a method for interpolating geological variable attributes according to an embodiment of the invention.
As shown in fig. 1, the method for interpolating geologic variable attributes includes:
step 1: loading the explained geologic body result data in a set work area range;
step 2: constructing a surface grid model of the geologic body according to the control point data of the geologic body;
the method comprises the following steps of constructing a surface mesh model of the complex geologic body by using a geometric modeling method, wherein the geometric modeling method comprises the following steps: two-dimensional limitation Delaunay subdivision, curved surface intersection, curved surface stitching and complex geologic body rapid tracking method;
and step 3: taking the triangular grid data of the block geometric surface as a limiting condition, and carrying out regularization processing on the limiting condition to obtain regular limiting points, limiting line segments and limiting surfaces and obtain a regularized surface grid model;
step 3 comprises steps 301-305:
step 301: setting a constraint line set and a discrete point set according to the triangular grid data;
step 302: recording the end points of all the line segments contained in the constraint line segment set into the discrete point set;
step 303: selecting any two line segments in the constraint line segment set, judging whether the two line segments are intersected or not, if the two line segments are intersected and collinear, recording collinear parts into the constraint line segment set, deleting the two line segments in the constraint line segment set, and if the two line segments are only intersected and not collinear, recording intersection points into a discrete point set;
step 304: if the triangular mesh does not meet the regularization requirement, adding a new point, judging whether the new point belongs to the discrete point set, and if not, recording the new point into the discrete point set; judging whether the line segments in the constraint line segment set contain new points or not, if so, taking the new points as end points, dividing the line segments containing the new points into two new line segments, recording the new line segments into the constraint line segment set, and deleting the line segments containing the new points;
step 305: repeating the steps 303-304 until all the triangular meshes meet the regularization requirement, obtaining regular limit points, limit line segments and limit surfaces, and obtaining a regularization surface mesh model consisting of the regular limit points, the limit line segments and the limit surfaces;
and 4, step 4: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency;
step 4 comprises steps 401-403:
step 401: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes;
step 402: checking whether each point, each edge and each triangle in the triangle mesh data of the limited condition are in the subdivision result or not, and if not, respectively adding the intersection points into the current triangle mesh and the tetrahedral mesh;
step 403: optimizing the tetrahedral mesh by adopting a mesh changing method to obtain the tetrahedral mesh with topological consistency;
the variable mesh method is to increase auxiliary points according to the geologic body structure form to construct a tetrahedral mesh;
and 5: and (5) performing attribute interpolation on the surface volume grid obtained in the step (4) by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result.
The attribute interpolation is carried out on the complex geologic body by adopting a gradient-based variable attribute interpolation algorithm, wherein the formula used by the variable attribute interpolation algorithm is as follows:
SampleValue[i]＝mBaseValue+mGX*gx+mGY*gy+mGZ*gz (1)
wherein, sampleValue [ i ] represents the attribute interpolation result of the ith tetrahedron, mBaseValue represents the basic attribute value of the block body where the tetrahedron to be inserted is located, mGX, mGY and mGZ respectively represent the attribute change gradient in the X, Y, Z coordinate direction, and gx, gy and gz respectively represent the coordinate increment in the X, Y, Z coordinate direction.
FIG. 2 illustrates loading a complex geological control point map based on a work area according to one embodiment of the invention. FIG. 3 illustrates a regularized surface mesh model of a method for geologic variable property interpolation according to one embodiment of the present invention. FIG. 4 shows a tetrahedral mesh with topological consistency for a method of geologic variable property interpolation according to one embodiment of the present invention. FIG. 5 is a diagram illustrating variable attribute interpolation results of a geological variable attribute interpolation method according to an embodiment of the present invention.
As shown in fig. 2, the coordinates of the work area range control points (0,0,0), (16000, 32000,0), (16000, 32000, 16000) are input to form a work area contour line, and a control point diagram formed by complex geological body result data which has been explained by geological personnel in the work area range is loaded.
As shown in fig. 3, after the regularization is performed by using the block geometric surface triangular mesh data as the limiting condition, the obtained regularized surface mesh model is relatively sparse, and only the surface of the model has lines.
As shown in fig. 4, a three-dimensional limited DELAUNAY subdivision algorithm and optimization processing are performed to obtain a tetrahedral mesh with topology consistency, the mesh is denser than a regular surface mesh model, and points and lines are subdivided and optimized inside the model to be full.
As shown in fig. 5, attribute interpolation is performed by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result, and the tetrahedral mesh after attribute difference is performed becomes more dense.
While embodiments of the present invention have been described above, the above description is illustrative, not exhaustive, and not limited to the disclosed embodiments. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.
Claims (7)
1. A geological variable attribute interpolation method is characterized by comprising the following steps:
step 1: loading the explained geologic body result data in a set work area range;
step 2: constructing a surface grid model of the geologic body according to the control point data of the geologic body;
and step 3: taking the triangular grid data of the block geometric surface as a limiting condition, and carrying out regularization processing on the limiting condition to obtain regular limiting points, limiting line segments and limiting surfaces and obtain a regularized surface grid model;
and 4, step 4: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency;
and 5: performing attribute interpolation on the tetrahedral mesh obtained in the step (4) by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result;
constructing a surface mesh model of the complex geologic body by using a geometric modeling method, wherein the geometric modeling method comprises the following steps: two-dimensional limitation Delaunay subdivision, curved surface intersection, curved surface stitching and complex geologic body rapid tracking method;
the step 3 specifically comprises:
step 301: setting a constraint line set and a discrete point set according to the triangular grid data;
step 302: recording the end points of all line segments contained in the constraint line segment set into the discrete point set;
step 303: selecting any two line segments in the constraint line segment set, judging whether the two line segments are intersected or not, if the two line segments are intersected and collinear, recording collinear parts into the constraint line segment set, deleting the two line segments in the constraint line segment set, and if the two line segments are only intersected and not collinear, recording intersection points into the discrete point set;
step 304: if the triangular mesh does not meet the regularization requirement, adding a new point, judging whether the new point belongs to the discrete point set, and if not, recording the new point into the discrete point set; judging whether the line segments in the constraint line segment set contain the new points or not, if so, dividing the line segments containing the new points into two new line segments by taking the new points as end points, recording the two new line segments into the constraint line segment set, and deleting the line segments containing the new points;
step 305: and repeating the steps 303-304 until all the triangular meshes meet the regularization requirement, obtaining regular limit points, limit line segments and limit surfaces, and obtaining a regularization surface mesh model consisting of the regular limit points, the limit line segments and the limit surfaces.
2. The method of claim 1, wherein the determination of whether the regularization requirement is satisfied is based on the area of the triangle and the minimum angle of the triangle.
3. The method of attribute-variant interpolation according to claim 1, wherein step 4 specifically comprises:
step 401: a three-dimensional limited DELAUNAY subdivision algorithm is adopted to subdivide the regularized surface mesh model to obtain a tetrahedral mesh;
step 402: checking whether each point, each edge and each triangle in the triangle mesh data of the limited condition are in the subdivision result or not, and if not, respectively adding the intersection points into the current triangle mesh and the tetrahedral mesh;
step 403: and optimizing the tetrahedral mesh by adopting a variable mesh method to obtain the tetrahedral mesh with topological consistency.
4. The method of claim 3, wherein the method of varying the mesh is to add auxiliary points according to the geologic structure to construct a tetrahedral mesh.
5. The method of claim 1, wherein the attribute interpolation is performed on the complex geologic body by using a gradient-based variable attribute interpolation algorithm, and the formula used by the variable attribute interpolation algorithm is as follows:
SampleValue[i] ＝ mBaseValue +mGX*gx +mGY*gy + mGZ*gz (1)
wherein, sampleValue [ i ] represents the attribute interpolation result of the ith tetrahedron, mBaseValue represents the basic attribute value of the block where the tetrahedron to be inserted is located, mGX, mGY and mGZ respectively represent the attribute change gradient in the X, Y, Z coordinate direction, and gx, gy and gz respectively represent the coordinate increment in the X, Y, Z coordinate direction.
6. A geological variable attribute interpolation system, comprising: a memory storing computer-executable instructions;
a processor executing computer executable instructions in the memory to perform the steps of:
step 1: loading the explained geologic body result data in a set work area range;
step 2: constructing a surface grid model of the geologic body according to the control point data of the geologic body;
and step 3: taking the triangular grid data of the block geometric surface as a limiting condition, and carrying out regularization processing on the limiting condition to obtain regular limiting points, limiting line segments and limiting surfaces and obtain a regularized surface grid model;
and 4, step 4: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain tetrahedral meshes with topological consistency;
and 5: performing attribute interpolation on the tetrahedral mesh obtained in the step 4 by using a gradient-based variable attribute interpolation algorithm to obtain an interpolation result;
the step 3 specifically comprises the following steps:
step 301: setting a constraint line set and a discrete point set according to the triangular grid data;
step 302: recording the end points of all the line segments contained in the constraint line segment set into the discrete point set;
step 303: selecting any two line segments in the constraint line segment set, judging whether the two line segments are intersected or not, if the two line segments are intersected and collinear, recording collinear parts into the constraint line segment set, deleting the two line segments in the constraint line segment set, and if the two line segments are only intersected and not collinear, recording intersection points into the discrete point set;
step 304: if the triangular mesh does not meet the regularization requirement, adding a new point, judging whether the new point belongs to the discrete point set, and if not, recording the new point into the discrete point set; judging whether the line segments in the constraint line segment set contain the new points or not, if so, dividing the line segments containing the new points into two new line segments by taking the new points as end points, recording the two new line segments into the constraint line segment set, and deleting the line segments containing the new points;
step 305: and repeating the steps 303-304 until all the triangular meshes meet the regularization requirement, obtaining regular limit points, limit line segments and limit surfaces, and obtaining a regularization surface mesh model consisting of the regular limit points, the limit line segments and the limit surfaces.
7. The system of claim 6, wherein step 4 comprises:
step 401: adopting a three-dimensional limited DELAUNAY subdivision algorithm to subdivide the regularized surface mesh model to obtain a tetrahedral mesh;
step 402: checking whether each point, each edge and each triangle in the triangle grid data of the limited condition are in the subdivision result, if not, respectively adding the intersection points into the current triangle grid and the tetrahedral grid;
step 403: and optimizing the tetrahedral mesh by adopting a variable mesh method to obtain the tetrahedral mesh with topological consistency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811159460.3A CN110968930B (en) | 2018-09-30 | 2018-09-30 | Geological variable attribute interpolation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811159460.3A CN110968930B (en) | 2018-09-30 | 2018-09-30 | Geological variable attribute interpolation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110968930A CN110968930A (en) | 2020-04-07 |
| CN110968930B true CN110968930B (en) | 2023-03-21 |
Family
ID=70028920
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811159460.3A Active CN110968930B (en) | 2018-09-30 | 2018-09-30 | Geological variable attribute interpolation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110968930B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113658333B (en) * | 2020-05-12 | 2024-03-01 | 中国石油化工股份有限公司 | Geologic body modeling method based on isosurface extraction |
| CN112037326B (en) * | 2020-08-11 | 2023-09-05 | 北京梦之岩科技有限公司 | Method for accurately describing adjacent spatial attribute condition of boundary of geologic body |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101105864A (en) * | 2007-08-07 | 2008-01-16 | 中国地质大学（武汉） | Tri-dimensional geologic body dynamic reconfiguration method based on sequence profile topology inference |
| WO2008029420A1 (en) * | 2006-09-04 | 2008-03-13 | Geosystem S.R.L. | Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data |
| CN102194253A (en) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | Method for generating tetrahedron gridding for three-dimensional geological structure |
| CN104240296A (en) * | 2014-07-24 | 2014-12-24 | 武汉地大信息工程股份有限公司 | Multi-tool collaborative three-dimensional geological modeling method and system for multi-source data |
| CN105388524A (en) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | Complex geologic body precise description taking triangular surface grid as limiting condition |
-
2018
- 2018-09-30 CN CN201811159460.3A patent/CN110968930B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008029420A1 (en) * | 2006-09-04 | 2008-03-13 | Geosystem S.R.L. | Method for building velocity models for pre-stack depth migration via the simultaneous joint inversion of seismic, gravity and magnetotelluric data |
| CN101105864A (en) * | 2007-08-07 | 2008-01-16 | 中国地质大学（武汉） | Tri-dimensional geologic body dynamic reconfiguration method based on sequence profile topology inference |
| CN102194253A (en) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | Method for generating tetrahedron gridding for three-dimensional geological structure |
| CN104240296A (en) * | 2014-07-24 | 2014-12-24 | 武汉地大信息工程股份有限公司 | Multi-tool collaborative three-dimensional geological modeling method and system for multi-source data |
| CN105388524A (en) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | Complex geologic body precise description taking triangular surface grid as limiting condition |
Non-Patent Citations (2)
| Title |
|---|
| 《地震成像块体建模技术研究》;孟祥宾;《油气地球物理》;全文 * |
| 《复杂地质体几何表面建模方法研究》;王鑫;《油气地球物理》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110968930A (en) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2869096B1 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| EP2653893B1 (en) | Faulted geological structures containing unconformities | |
| 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 | |
| US6052650A (en) | Enforcing consistency in geoscience models | |
| EP3371629B1 (en) | Representing structural uncertainty in a mesh representing a geological environment | |
| CA2918493C (en) | Generating a custom reservoir from multiple compartments representing one or more geological structures | |
| US20130218539A1 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
| NO20171868A1 (en) | Geologic stratigraphy via implicit and jump functions | |
| CN106687827B (en) | Stratum modeling method for fault | |
| EP2631685A2 (en) | Building faulted grids for a sedimentary basin including structural and stratigraphic interfaces | |
| KR20150093391A (en) | Method for rendering terrain | |
| CN110968930B (en) | Geological variable attribute interpolation method and system | |
| CN111445569B (en) | Sedimentary geological evolution dynamic simulation method | |
| CN115795974A (en) | Air-ground integrated monitoring method for evolution of underground water seepage field of deep-buried tunnel | |
| CN111951394A (en) | Fault structure unit three-dimensional model construction method and device based on geological map | |
| CN106875484A (en) | A kind of geology accumulation body Fast Fitting modeling method based on dimensional topography | |
| CN107886573B (en) | Slope three-dimensional finite element grid generation method under complex geological conditions | |
| CN113658333B (en) | Geologic body modeling method based on isosurface extraction | |
| AU2016223281B2 (en) | Method for parameterizing a 3D domain with discontinuities | |
| CN107767451B (en) | 3DEC complex terrain modeling method based on GOCAD grid model | |
| Karimi-Fard | Growing region technique applied to grid generation of complex fractured porous media | |
| US10310116B2 (en) | Method for determining a stacked model describing architectural elements | |
| CN114491897A (en) | Seismic wave numerical simulation method, device, medium and electronic equipment | |
| EP1070267B1 (en) | Enforcing consistency in geoscience models | |
| CN117518245A (en) | Method, device, computer equipment and storage medium for establishing speed model |
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 |