CN107958483A - A kind of Stratum Modeling - Google Patents
A kind of Stratum Modeling Download PDFInfo
- Publication number
- CN107958483A CN107958483A CN201610898488.3A CN201610898488A CN107958483A CN 107958483 A CN107958483 A CN 107958483A CN 201610898488 A CN201610898488 A CN 201610898488A CN 107958483 A CN107958483 A CN 107958483A
- Authority
- CN
- China
- Prior art keywords
- triangle
- intersecting
- vertex
- crosspoint
- triangles
- 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
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 present invention provides a kind of Stratum Modeling, including：It will explain that the geological surface of reconstruct is decomposed into continuous space triangle based on seismotectonics；Friendship is asked to generate crosspoint two-by-two on the triangle of there may exist overlapping relation；Using triangle as foundation, the intersecting lens in geological surface is generated based on crosspoint；Subdivision is carried out again to the triangle comprising crosspoint based on intersecting lens to update geological surface；Dough sheet division is carried out to the geological surface of renewal using the side of intersecting lens and geological surface as border；Using intersecting lens as suture, the geological surface dough sheet generation geological surface closed shape based on division.The present invention can reduce the man-machine interactively operation during Geologic modeling, can effectively improve the precision of geological model.
Description
Technical field
The invention belongs to the three-dimensional geological modeling technical field in geophysical prospecting for oil, specifically, more particularly to
A kind of Stratum Modeling based on space curved surface topological analysis.
Background technology
Three-dimensional geological modeling is a kind of under three-dimensional environment, uses the information management of computer technology general space, geology solution
Release, spatial analysis and prediction, soil statistics, the multiple types of tools such as physical contents analysis and graph visualization and multinomial technology carry out
The leading edge integrated technology of geological analysis.
Stratum object will be expressed and shown to three-dimensional geological modeling not only in itself, but also express between the internal object of stratum
Correlation (such as topological relation and semantic relation), in addition also wants the attribute information of storage object.At present, national system is not only lacked
One standard stratigraphic model, but also can be by theory and model that most people receives without being formed.Therefore, three-dimensional geological is built
Mould still need further research with it is perfect.It is also, also more, smart there are man-machine interactively in existing Geologic modeling construction modeling process
Spend the problem of error is big.
The content of the invention
To solve problem above, the present invention provides a kind of Stratum Modeling, during can reducing Geologic modeling
Man-machine interactively operates, and can effectively improve the precision of geological model.
According to one embodiment of present invention, there is provided a kind of Stratum Modeling, including：
It will explain that the geological surface of reconstruct is decomposed into continuous space triangle based on seismotectonics；
Friendship is asked to generate crosspoint two-by-two on the triangle of there may exist overlapping relation；
Using the triangle as foundation, the intersecting lens in geological surface is generated based on the crosspoint；
Subdivision is carried out again to the triangle comprising the crosspoint based on the intersecting lens to update geological surface；
Dough sheet division is carried out to the geological surface of renewal using the side of the intersecting lens and geological surface as border；
Using the intersecting lens as suture, the geological surface dough sheet generation geological surface closed shape based on division.
According to one embodiment of present invention, friendship is asked two-by-two to the triangle of there may exist overlapping relation
Further comprised with generating crosspoint：
Judgement ask whether two triangles of friendship overlap, such as it is misaligned, calculate one of triangle three vertex and
The directed distance of plane where another triangle；
The contact type and contact relation of crosspoint are judged according to the directed distance.
According to one embodiment of present invention, using the triangle as foundation, geological surface is generated based on the crosspoint
In intersecting lens further comprise：
In individual geological surface, friendship triangle is asked comprising identical to asking friendship triangle to be arranged and connected according to crosspoint
The crosspoint of shape is to generate intersecting line segment；
Intersecting line segment is arranged by endpoint order, merges the intersecting line segment comprising same endpoints to generate intersecting lens；
Intersecting lens is arranged by endpoint order, merges the intersecting lens comprising same endpoints to generate new intersecting lens；
All intersecting line segments and intersecting lens in all geological surfaces are traveled through until new intersecting lens cannot be produced.
According to one embodiment of present invention, the triangle comprising the crosspoint is carried out again based on the intersecting lens
Subdivision is further comprised with updating geological surface：
In individual geological surface, the triangle of friendship is sought according to position division of the crosspoint in friendship triangle is sought；
It is that vertex of a triangle and intersecting lens are the side of triangle for constraints using the crosspoint, to the ground of reconstruct
Matter curved surface carries out triangle subdivision.
According to one embodiment of present invention, the triangle of friendship is sought according to position division of the crosspoint in the triangle of friendship is sought
Shape further comprises：
Such as crosspoint is in a line of triangle for asking friendship, then using the crosspoint as vertex, by the triangle point
Split for two common side triangles.
According to one embodiment of present invention, the triangle of friendship is sought according to position division of the crosspoint in the triangle of friendship is sought
Shape further comprises：
If the crosspoint is in the inside for the triangle for asking friendship, then using the crosspoint as vertex, which is split into
Three triangles.
According to one embodiment of present invention, it is bent to the geology of renewal as border using the side of the intersecting lens and geological surface
Face carries out dough sheet division and further comprises：
The intersection edges for all triangles for asking friendship for extracting an intersecting lens and obtaining on the intersecting lens；
An intersection edges on the intersecting lens are marked, the associated intersecting triangles of the intersection edges is marked and travels through associated phase
Hand over all sides of triangle；
According to associated intersecting triangles while attribute mark this while associated triangle, it is more until that can not mark
Triangle is so as to obtain a geological surface dough sheet；
Travel through all intersecting lenses and mark the associated intersecting triangles of intersection edges of all triangles for asking friendship, to obtain
All dough sheets of geological surface.
According to one embodiment of present invention, according to associated intersecting triangles while attribute mark this while associated three
It is angular to further comprise：
Side such as associated intersecting triangles is an intersecting line segment, stops traversal；
Frontier juncture such as associated intersecting triangles joins a triangle, then the triangle is located at geological surface border, marks
The triangle simultaneously stops traveling through；
Side such as associated intersecting triangles is not intersecting line segment and associates two triangles, then marks another not visit
The triangle asked, and travel through all sides for not accessing triangle.
According to one embodiment of present invention, using the intersecting lens as suture, the geological surface dough sheet life based on division
Further comprise into geological surface closed shape：
Select one not access geological surface dough sheet as front piece is worked as, mark dough sheet access order and dough sheet is put into access number
Group；
Selection finds out next dough sheet when a unmarked intersecting lens of front piece according to the dough sheet access order；
The access direction of next dough sheet is determined according to the direction of current intersecting lens；
If next dough sheet is not put into access array, which is put into access array；
It is when the inside line segment of front piece and next dough sheet to mark the intersecting lens；
Judge whether all access orders of next dough sheet are all present in access in array, if all existed, lower one side
Piece is outstanding dough sheet；
Unmarked intersecting lens is traveled through to find next dough sheet successively, until when all intersecting lenses of front piece are labeled as inside
Line segment, removes all outstanding dough sheets and merges all dough sheets, completes the construction of a closed shape；
All geological surface dough sheets that do not access are traveled through until all geological surface dough sheets are collectively labeled as having accessed, then institute
There is the generation of geological surface closed shape to finish.
According to one embodiment of present invention, the contact type and contact relation of crosspoint are judged according to the directed distance
Further comprise：
As the directed distance on three vertex is same, just or with bearing, then two triangles are non-intersect, and two triangles are not intersecting
Point；
Directed distance such as three vertex is all zero, then two triangles are coplanar；
If the directed distance on two vertex is zero, the directed distance on another vertex is just or is negative, a then triangle
There are a line another triangle institute in the plane, planar judge two triangles in another triangle
Intersecting situation；
Directed distance such as a vertex is zero, and the directed distance on another two vertex is same just or with bearing, then a triangle
There are a vertex another triangle institute in the plane, another triangle planar judge the vertex with
The intersecting situation of another triangle；
Directed distance such as a vertex is zero, and the directed distance on another two vertex is negative for one positive one, then a triangle
There are a vertex another triangle institute in the plane, while the directed distance be zero vertex opposite side and another
The place plane of triangle judges the intersection point and another triangle in the plane there are intersection point in another triangle
Intersecting situation；
If the directed distance on vertex is just, the directed distance on another two vertex is with negative, or vertex is oriented
Distance is negative, and the directed distance on another two vertex is with just, then there are two sides and the place of another triangle for a triangle
Plane intersects, side where taking two vertex of directed distance contrary sign, the intersection point of plane where calculating the side and another triangle,
And in the intersecting situation for judging the intersection point and another triangle in the plane of another triangle.
Beneficial effects of the present invention：
Stratum Modeling provided by the invention based on space curved surface topological analysis, during can reducing Geologic modeling
Man-machine interactively operation, the precision of geological model can be effectively improved.Calculated by geological surface space contact relation and stratum is closed
Fit generation, and space representation is carried out to stratum with visualization technique, it can be used in oil and gas reservoir evaluation of trap and reservoir commented
Valency, instructs identification of hydrocarbon to hide Distribution Characteristics and carries out Beneficial Zones of Exploring selection.
Other features and advantages of the present invention will be illustrated in the following description, also, partly becomes from specification
Obtain it is clear that or being understood by implementing the present invention.The purpose of the present invention and other advantages can be by specification, rights
Specifically noted structure is realized and obtained in claim and attached drawing.
Brief description of the drawings
To describe the technical solutions in the embodiments of the present invention more clearly, it is required in being described below to embodiment
Attached drawing does simple introduction：
Fig. 1 is method flow diagram according to an embodiment of the invention；
Fig. 2 is triangle topology relation schematic diagram according to an embodiment of the invention；
Fig. 3 is intersecting line segment and triangle topology relationship analysis schematic diagram according to an embodiment of the invention；
Fig. 4 is intersecting lens generation schematic diagram according to an embodiment of the invention；
Fig. 5 is geological surface dough sheet division schematic diagram according to an embodiment of the invention；
Fig. 6 is four face schematic diagrames used in closed shape organigram according to an embodiment of the invention；
Fig. 7 is that four curved surfaces intersect schematic diagram two-by-two in Fig. 6 according to an embodiment of the invention；
Fig. 8 is reconstruct and the dough sheet division schematic diagram of four curved surfaces according to an embodiment of the invention；
Fig. 9 is geological model closed shape generation schematic diagram according to an embodiment of the invention；
Figure 10 a are geological surface topological analysis schematic diagrames according to an embodiment of the invention；
Figure 10 b are single stratum closed shape schematic diagrames according to an embodiment of the invention；
Figure 10 c are stratigraphic model schematic diagrames according to an embodiment of the invention.
Embodiment
Carry out the embodiment that the present invention will be described in detail below with reference to accompanying drawings and embodiments, how the present invention is applied whereby
Technological means solves technical problem, and that reaches technique effect realizes that process can fully understand and implement according to this.Need to illustrate
As long as not forming conflict, each embodiment in the present invention and each feature in each embodiment can be combined with each other,
The technical solution formed is within protection scope of the present invention.
Current geological model becomes the digitized important means in oil field, and learns information technology field as current position
One of most dynamic research direction.In the geologic structure interpretation stage, it is main that three-dimensional geological modeling can make up structure interpretation
Dependent on plan and profile and need to project to the stratum in three dimensions, construction, landforms and other geological phenomenons
Deficiency in a certain plane, by means of visualization technique, directly goes to understand and expresses geologic body and ground from the angle of three dimensions
Matter environment, effectively compensate for the loss of spatial information in conventional method and distortion, drawing course is numerous and diverse and information updating is difficult
The defects of, complicated geological structure can be clearly descended in a more intuitive way, it is effective to improve geological personnel to block
The understanding that texture is made.
Meanwhile in the Reservoir Development stage, integrated geological interpretation task is extremely complex, and it is all to be related to geology, earthquake, well logging etc.
It is multidisciplinary, and three-dimensional geological model can be explained with synthesis structure, stratigraphic interpretation, comprehensive log interpretation, seismic attributes analysis, reservoir
The Data Analysis Results such as comprehensive characterization (pool description), carry out reservoir frame model, Reservoir Parameter Models and oil-gas reservoir change
Dynamic analog etc. research, provide decision-making foundation for the exploitation of oil-gas reservoir.It is also, continuous with oil reservoir development degree
Improve, establish fine hydrocarbon reservoir geologic model, be also of great significance for improving oil reservoir understanding.
Therefore, the present invention provides a kind of Stratum Modeling based on space curved surface topological analysis.The present invention is by right
Form geological surface spatial triangle and asked friendship two-by-two, obtain all crosspoints, and close according to crosspoint and triangle position
System, generates the intersecting lens between any two geological surfaces, realizes the space topology analyzing of geological surface, and utilize topological analysis
As a result, carrying out burst and stitching processing to geological surface, realize stratum block closed shape generation, realize stratigraphic model.The present invention
The achievement of seismotectonics explanation is taken full advantage of, error caused by avoiding man-machine interactively means, farthest remains ground
The description precision that texture is made, epistemically can descend complicated geological structure to provide more accurate model and foundation for geological personnel.
It is the method flow diagram according to one embodiment of the present of invention as shown in Figure 1, to carry out specifically the present invention below with reference to Fig. 1
It is bright.
It is step S100 first, will explains that the geological surface of reconstruct is decomposed into continuous space triangle based on seismotectonics.
Point is explained specifically, being explained to obtain by seismotectonics, by explaining point reconstruct geological surface.Then by the geological surface after reconstruct
It is decomposed into continuous spatial triangle structure.Contained in geological surface information after decomposition spatial triangle vertical array and
Vertex of a triangle indexes.
Followed by step S200, asked friendship intersecting to generate two-by-two on the triangle of there may exist overlapping relation
Point.
Triangle there may be overlapping relation can obtain in the following manner.First successively to the geological surface of reconstruct
It is detected according to spatial dimension, judges the curved surface there are overlapping relation.Then to all on the curved surface there are overlapping relation
Triangle is detected according to spatial dimension, judges the triangle there are overlapping relation.
Crosspoint is obtained by calculating the intersection point of two triangles.Contained in the information of crosspoint crosspoint position,
The topological relation of triangle information representing and crosspoint and triangle where crosspoint.Crosspoint position refers to the sky of crosspoint
Between coordinate.Triangle information representing refers to triangular apex and triangle edges.Topological relation refers to that the crosspoint is located at two three respectively
Position on angular, for single triangle, including crosspoint is located at vertex of a triangle, side and positioned at triangle interior three
Kind position relationship.
The crosspoint for calculating two triangles specifically includes procedure below.First determine whether to ask whether two triangles of friendship weigh
Close, if overlapped, present invention determine that method is invalid；As misaligned, three of one of triangle (being set to A ') are calculated
Vertex and the directed distance of plane (being set to P ') where another triangle (being set to B ').Then, two are judged according to directed distance
The contact type and contact relation of the crosspoint of a triangle.The contact type and contact relation of crosspoint are specifically divided into following
Six kinds of situations.
Situation 1：The directed distance on three vertex is with just or with bearing, then two triangles A ' and B ' are non-intersect, two triangles
Shape A ' and B ' do not have crosspoint, as a parts show the contact relation of two triangles in Fig. 2.It is empty in the attached drawing of the present invention
Line represent projection relation ,+number ,-number and 0 represent directed distance, this method directly using geometric algorithm judge topological relation, than
The position relationship on vertex and vertex, vertex and side more can be accurately judged using algebraic equation.Other accompanying drawings are using similar
Geometric algorithm.
Situation 2：The directed distance on three vertex is all zero, then two triangles A ' and B ' are coplanar, such as b parts institute in Fig. 2
Show.
Situation 3：The directed distance on two vertex is zero, and the directed distance on another vertex is just or is negative, then triangle
B ' on plane P ', judges the intersecting situation of two triangles, as shown in c parts in Fig. 2 there are a line in plane P '.
Situation 4：The directed distance on one vertex is zero, and the directed distance on another two vertex is same just or with bearing, then triangle
B ' on plane P ', judges the vertex and the intersecting situation of triangle A ' there are a vertex on plane P ', such as d portions in Fig. 2
Shown in point.
Situation 5：The directed distance on one vertex is zero, and the directed distance on another two vertex is negative for one positive one, then triangle
B ' there are a vertex on plane P ', while the directed distance be zero vertex opposite side and plane P ' there are intersection point, flat
The intersection point and the intersecting situation of triangle A are judged on the P ' of face, as shown in e parts in Fig. 2.
Situation 6：The directed distance on one vertex is just, and the directed distance on another two vertex is with negative, alternatively, vertex
Directed distance is negative, and the directed distance on another two vertex is with just, then B ' takes directed distance different there are two sides and plane P ' are intersecting
Number two vertex where side, calculate the intersection point of the side and plane P ', and judge on plane P ' intersection point and triangle A '
Intersecting situation, as shown in f parts in Fig. 2.
Followed by step S300, using triangle as foundation, the intersecting lens in geological surface is generated based on crosspoint.
Intersecting lens in geological surface is formed by continuously intersecting line segment, and intersecting line segment is intersected by two triangles
What obtained intersection point was formed by connecting.Generation intersecting lens comprises the following steps that.
First, in individual geological surface, asked according to crosspoint comprising identical asking friendship triangle to be arranged and connected
The crosspoint of triangle is handed over to generate intersecting line segment, as shown in b parts in Fig. 4, a parts show two triangles in Fig. 4
Crosspoint.Specifically, corresponding to the position relationship of two triangles in Fig. 2, the intersecting line segment formed by two triangles
Position relationship is as shown in Figure 3.Wherein, the triangle table in Fig. 3 is shown as corresponding ask of the plane and hands over triangle A ', solid black line
Segment table show be another triangle B ' a line projection in the plane, horizontal dotted line represented where the Projection Line Segment on the side
Straight line, vertical dotted line represent projection straight line and the scope R there may be crosspoint of triangle A '.A plurality of dotted line in Fig. 3
What is represented is projection side and region R different position relationships that may be present and specific intersection topology situation.Wherein region R
It is the scope blocked by delta-shaped region of straight line projected where side, each several part in corresponding diagram 3 is analyzed, and there are following feelings
Condition：(1) in b parts region R corresponding flat triangles A ' a line, region R corresponding flat triangles A ' in (2) c parts
One vertex, region R are reduced to a point, and in (3) d parts a vertex and the vertex of region R corresponding flat triangles A '
On opposite side a little between line segment, between two points in (4) e parts on two sides of region R corresponding flat triangles A '
Line segment.Following analysis method can be taken for different region R, the side entered laterally from left to right using Projection Line Segment
Formula is understood, illustrates that specific intersection point only need to be oriented by what is obtained in Fig. 2 mainly for obtaining complete topology status below
Distance, can be obtained by, projecting side and the region R situations of intersecting has following state by calculating the point that directed distance is 0：(1) throw
Shadow side does not have intersection point in left side and with region R；(2) the left margin point for the right endpoint and region R for projecting side intersects, (3) projection side
Internal point and the left margin point of region R intersect, the left end point on (4) projection side and the left margin point of region R intersect, (5) projection
While falling inside the R of region, that is, project side and triangle interior is met at a bit, (6) project the right endpoint and region R right margins point on side
Intersecting, the internal point on (7) projection side and the right margin point of region R intersect, the left end point on (8) projection side and the right margin of region R
Point is intersecting, and (9) projection side does not have intersection point on right side and region R.Because project side it is corresponding be triangle B ' a line, because
This, the overlapping relation for the triangle A ' and B ' that can know for sure.It can do and as above divide for the projection side of each figure in Fig. 2
Analysis, so as to obtain a-e parts in Fig. 3.It can be obtained by the following topological relation of crosspoint by Fig. 3：By in side and triangle
Portion intersects to obtain, by two while intersecting to obtain (positioned at while inside, not comprising vertex), intersect to obtain by side and vertex or by
Vertex and vertex intersect to obtain.And complete triangle, side and vertex information can be recorded, so as to obtain complete topology letter
Breath.Herein on basis, crosspoint can be connected into intersecting line segment according to these information.
Then, intersecting line segment is arranged by endpoint order, merges the intersecting line segment comprising same endpoints to generate phase
Intersection.Then, intersecting lens is arranged by endpoint order, merges the intersecting lens comprising same endpoints to generate new intersect
Line.Finally, all intersecting line segments and intersecting lens in all geological surfaces are traveled through until new intersecting lens cannot be produced.
Followed by step S400, subdivision is carried out again to the triangle comprising crosspoint based on intersecting lens to update geology
Curved surface.Specifically, in individual geological surface, the topology information that is included according to crosspoint, the triangle to producing crosspoint
Reconstructed again.
First, in individual geological surface, friendship triangle is asked according to position division of the crosspoint in friendship triangle is sought.Root
According to the position of crosspoint in the triangles, friendship triangle is sought in division in two kinds of situation.If crosspoint is in a line of triangle,
Then using the crosspoint as vertex, which is split into two common side triangles.For example, triangle ABC includes crosspoint, A,
B, C is three vertex of triangle, and P is crosspoint, and triangle ABC is then split into triangle by P on the BC sides of triangle ABC
Shape ABP and the triangle on the common sides of triangle ACP two.If crosspoint is in the inside of triangle, then using the crosspoint as vertex, general
The triangle is split into three triangles.For example, P is split into triangle at the inside of triangle ABC, by triangle ABC
ABP, triangle ACP and triangle BCP.
Then, it is that vertex of a triangle and intersecting lens are the side of triangle for constraints using crosspoint, to reconstruct ground
Matter curved surface carries out triangle subdivision, and the geological surface after subdivision can so optimize the geometry of geological surface as shown in c in Fig. 4
Form.
Followed by step S500, using the side of intersecting lens and original reconstruct geological surface as border, to the geology of renewal
Curved surface carries out geological surface dough sheet division.Geological surface dough sheet refers to the connected region by intersecting lens and surface boundary encirclement, is
One independent space polygon region.The geological surface obtained in step S400 contains the vertex of a triangle for forming curved surface
Indexed with triangular apex.According to geometric theory, a triangle should associate three sides, and triangle can be established according to the theory
Shape and the topology information on side, complete the mapping on the side from triangle to triangle.For the associated triangle in side of triangle, deposit
In two kinds of situations：Such as a line one triangle of association, then the side is the border of geological surface；As a line associates two triangles
Shape, then the side is the common edge of two triangles.According to both of side and triangle incidence relations, side and triangle are established
Topology information, completes the mapping from the side of triangle to triangle.
In addition, understood according to the generation of intersecting lens description：Intersecting lens is made of intersecting line segment, intersects line segment by crosspoint group
Into crosspoint contains the information of intersecting triangles；The description being updated according to the crosspoint using generation to geological surface
Understand：Crosspoint must be the triangular apex of the geological surface of renewal, and intersecting line segment must be the side of triangle, and be phase
Hand over the common edge of triangle.According to the above description, the topology information of intersecting line segment and intersecting triangles can be established, is completed intersecting
Mapping of the side to triangle.
Information from the description above, the partition process of geological surface dough sheet include following steps.First, one is extracted
Intersecting lens simultaneously obtains all intersection edges on the intersecting lens.Then an intersection edges therein are marked, according to intersection edges mark pair
Answer associated intersecting triangles (as shown in a in Fig. 5, the triangle of an adjoining is found from an intersection edges) and travel through corresponding association
All sides of triangle.For each side of triangle, the situation that triangle is joined according to each frontier juncture determines how progress
Traversal.Specifically, such as side is an intersection edges, then stop traversal；Such as the frontier juncture joins a triangle, then the triangle is located at
Geological surface border, stops traversal；Such as the frontier juncture joins two triangles, the then triangle for marking another not access, and travels through
All sides of the triangle not accessed.Then the adjacent edge of triangle, base are found from all sides of corresponding association triangle
In adjacent edge and intersecting wire tag triangle until can not mark more triangles (as shown in b-e in Fig. 5, wherein, b illustrates root
According to the triangle found, find one and do not access side, and another adjoining triangle is found according to the side, so as to complete once
Diffusion, c-e illustrate multiple diffusion process).At this time markd triangle be exactly an internal connection region, be one
Dough sheet, as shown in f in Fig. 5.Above step is circulated, all intersecting lenses is traveled through and marks the intersection edges in all intersection edges to associate
Intersecting triangles, then can obtain all dough sheets on the geological surface.
It is finally step S600, using intersecting lens as suture, the geological surface dough sheet generation geological surface based on division closes
It is fit.Geological surface closed shape refers to the space geometry form surrounded by dough sheet, and space is divided into disconnected inside and outer
Two, portion region.
Repeated accesses and sequence consensus are, it is necessary to define direction in order to prevent.Present invention agreement uses right hand rule, this definition
It is as follows：The order for defining side (line segment) is the order on vertex；According to right hand rule, triangle is defined using the order on the side of triangle
The outer normal direction (positive direction) of shape；The direction of dough sheet is consistent with triangle direction；The direction of closed surface is consistent with dough sheet direction；
, it is necessary to respectively be accessed once along (or positive and negative) direction inside and outside dough sheet in closed shape generating process.
The geological surface dough sheet obtained by step S500, can obtain following information：Complete from dough sheet to triangle
The mapping of mapping, the mapping of dough sheet to intersecting lens and intersecting lens to dough sheet；One intersecting lens can associate multiple dough sheets, therefore need
The access order of multiple dough sheets is defined, by following agreement：With the vertex index order of intersecting lens, extract one therein and intersect
Line segment does perspective plane for normal direction, this is intersected the associated triangle of line segment projects to the perspective plane, with the side counterclockwise of projection line
To the access order to put in order as associated multiple dough sheets.
Closed shape generating process is described in detail below, and represents that positive direction accesses with (+), (-) represents that negative direction accesses.
One is selected not access dough sheet A (+) as front piece is worked as first, mark access order is positive direction, by the dough sheet
It is put into access array.If two access orders of all dough sheets are all labeled, closed shape generation terminates.
Then selection is found out next when a unlabelled intersecting lens L (+) of front piece A (+) according to dough sheet access order
A dough sheet B.If can not find next dough sheet, B=A (-) is marked, i.e., is accessed along negative direction and works as front piece, while set intersecting
The access order of line is negative direction, i.e. L (-).
The access direction of dough sheet B is determined then according to the direction of current intersecting lens.
Then, if dough sheet B is not put into access array, which is put into access array.
Then, the inside line segment that intersecting lens L is dough sheet A and dough sheet B is marked.
Next, it is determined that whether all access orders (present invention is by taking two access orders as an example) of dough sheet B are all present in visit
Ask in array, if all existed, dough sheet B is outstanding dough sheet.
Then travel through unmarked intersecting lens and find next dough sheet successively, until all intersecting lenses as front piece A are labeled as
Internal line segment, removes all outstanding dough sheets, and merges all dough sheets, then completes the construction of a closed shape.
All dough sheets that do not access are finally traveled through until all dough sheets are collectively labeled as having accessed, then all geological surfaces close
Zoarium generation finishes.
Illustrate the generating process of closed shape below by way of a specific example.Fig. 6 is participate in closed shape generation 4
Original face, is respectively top surface s1, above bottom surface s2, s3 and back side s4.Fig. 7 is that signal intersecting two-by-two is carried out to space curved surface
Figure, wherein intersection point (line) use the set expression of two curved surfaces, as s13 represents the intersection point (line) of s1 and s3.Fig. 8 is each face root
Intersection, and subdivision again are generated according to intersection point, ultimately generates the process of dough sheet, wherein dough sheet is represented using curved surface and index, such as
S1p1 represents first dough sheet of curved surface s1.Fig. 9 is the generating process of closed shape stratigraphic model, and it is each to have listed closed shape in figure
The anabolic process in face, wherein, 1.s1p3 (+) expression is the positive direction access order of the 3rd dough sheet of curved surface s1, wherein+represent
Positive direction access order, the access order of multiple dough sheets, is indicated in figure by counterclockwise arrow；(2.s1p3 -) table
Show the negative direction access order for the 3rd dough sheet for being curved surface s1, wherein-represent positive direction access order, because have accessed at the same time
S1p3 both forward and reverse directions, so s1p3 is to hang dough sheet, it is necessary to remove.
As shown in table 1, institute's orderliness in Fig. 9 is followed successively by：
Table 1
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
| s1p3(+) | s1p3(-) | s3p1(+) | s2p3(+) | s2p3(-) | s3p2(+) | s3p2(-) | s2p2(+) | s4p3(+) |
| 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 |
| s4p3(-) | s2p1(+) | s2p1(-) | s4p2(+) | s1p1(+) | s1p1(-) | s4p1(+) | s4p1(-) | s1p2(+) |
By deleting invalid outstanding dough sheet, the closed shape finally obtained is by tetra- dough sheet structures of s3p1, s2p2, s4p2, s1p2
Into.
It is the practical application of the closed shape generation based on space curved surface topological analysis as shown in Figure 10 a-10c, utilizes somewhere
The practice interpretation data in area have carried out closed shape generation, establish stratigraphic model.Wherein, Figure 10 a are carried out according to geological surface
For topological analysis as a result, one of closed shape in the final result that Figure 10 b are, Figure 10 c are the closed shapes on whole stratum
The common display of model and original geological surface, illustrates that topological coherence is effectively maintained.
Stratum modeling technology provided by the invention based on space curved surface topological analysis, during can reducing Geologic modeling
Man-machine interactively operation, the precision of geological model can be effectively improved.By to the calculating of geology curved space contact relation and stratum
The generation of closed shape, and space representation is carried out to stratum with visualization technique, it can be used in oil and gas reservoir evaluation of trap and reservoir
Evaluation, and guides identification of hydrocarbon to hide Distribution Characteristics and carries out Beneficial Zones of Exploring selection.
While it is disclosed that embodiment as above, but the content only to facilitate understand the present invention and adopt
Embodiment, is not limited to the present invention.Any those skilled in the art to which this invention pertains, are not departing from this
On the premise of the disclosed spirit and scope of invention, any modification and change can be made in the implementing form and in details,
But the scope of patent protection of the present invention, still should be subject to the scope of the claims as defined in the appended claims.
Claims (10)
1. a kind of Stratum Modeling, including：
It will explain that the geological surface of reconstruct is decomposed into continuous space triangle based on seismotectonics；
Friendship is asked to generate crosspoint two-by-two on the triangle of there may exist overlapping relation；
Using the triangle as foundation, the intersecting lens in geological surface is generated based on the crosspoint；
Subdivision is carried out again to the triangle comprising the crosspoint based on the intersecting lens to update geological surface；
Dough sheet division is carried out to the geological surface of renewal using the side of the intersecting lens and geological surface as border；
Using the intersecting lens as suture, the geological surface dough sheet generation geological surface closed shape based on division.
2. according to the method described in claim 1, it is characterized in that, to the triangle of there may exist overlapping relation into
Row asks friendship to further comprise to generate crosspoint two-by-two：
Judgement asks whether two triangles of friendship overlap, as misaligned, calculate three vertex of one of triangle with it is another
The directed distance of plane where a triangle；
The contact type and contact relation of crosspoint are judged according to the directed distance.
3. method according to claim 1 or 2, it is characterised in that using the triangle as foundation, based on the crosspoint
Intersecting lens in generation geological surface further comprises：
In individual geological surface, triangle is handed over comprising identical ask to asking to hand over triangle to be arranged and connect according to crosspoint
Crosspoint is to generate intersecting line segment；
Intersecting line segment is arranged by endpoint order, merges the intersecting line segment comprising same endpoints to generate intersecting lens；
Intersecting lens is arranged by endpoint order, merges the intersecting lens comprising same endpoints to generate new intersecting lens；
All intersecting line segments and intersecting lens in all geological surfaces are traveled through until new intersecting lens cannot be produced.
4. method according to any one of claim 1-3, it is characterised in that based on the intersecting lens to including the phase
The triangle of intersection point carries out again subdivision and is further comprised with updating geological surface：
In individual geological surface, the triangle of friendship is sought according to position division of the crosspoint in friendship triangle is sought；
It is that vertex of a triangle and intersecting lens are the side of triangle for constraints using the crosspoint, it is bent to the geology of reconstruct
Face carries out triangle subdivision.
5. according to the method described in claim 4, it is characterized in that, divided according to position of the crosspoint in the triangle of friendship is sought
The triangle of friendship is asked to further comprise：
Such as crosspoint, then using the crosspoint as vertex, which is split into a line of triangle for asking friendship
Two common side triangles.
6. according to the method described in claim 4, it is characterized in that, divided according to position of the crosspoint in the triangle of friendship is sought
The triangle of friendship is asked to further comprise：
If the crosspoint is in the inside for the triangle for asking friendship, then using the crosspoint as vertex, which is split into three
Triangle.
7. according to the method any one of claim 1-6, it is characterised in that with the side of the intersecting lens and geological surface
Dough sheet division is carried out for border to the geological surface of renewal to further comprise：
The intersection edges for all triangles for asking friendship for extracting an intersecting lens and obtaining on the intersecting lens；
An intersection edges on the intersecting lens are marked, the associated intersecting triangles of the intersection edges is marked and travels through associated intersecting three
Angular all sides；
According to associated intersecting triangles while attribute mark this while associated triangle, until more triangles can not be marked
Shape is so as to obtain a geological surface dough sheet；
Travel through all intersecting lenses and mark the associated intersecting triangles of intersection edges of all triangles for asking friendship, to obtain geology
All dough sheets of curved surface.
8. the method according to the description of claim 7 is characterized in that should according to the attribute on the side of associated intersecting triangles mark
The associated triangle in side further comprises：
Side such as associated intersecting triangles is an intersecting line segment, stops traversal；
As associated intersecting triangles frontier juncture join a triangle, then the triangle be located at geological surface border, mark this three
It is angular and stop traveling through；
Side such as associated intersecting triangles is not intersecting line segment and associates two triangles, then marks what another was not accessed
Triangle, and travel through all sides for not accessing triangle.
9. according to the method any one of claim 1-8, it is characterised in that using the intersecting lens as suture, be based on
The geological surface dough sheet generation geological surface closed shape of division further comprises：
Select one not access geological surface dough sheet as front piece is worked as, mark dough sheet access order and dough sheet is put into access array；
Selection finds out next dough sheet when a unmarked intersecting lens of front piece according to the dough sheet access order；
The access direction of next dough sheet is determined according to the direction of current intersecting lens；
If next dough sheet is not put into access array, which is put into access array；
It is when the inside line segment of front piece and next dough sheet to mark the intersecting lens；
Judge whether all access orders of next dough sheet are all present in access in array, if all existed, next dough sheet is
Outstanding dough sheet；
Unmarked intersecting lens is traveled through to find next dough sheet successively, until all intersecting lenses when front piece are labeled as inner wire
Section, removes all outstanding dough sheets and merges all dough sheets, completes the construction of a closed shape；
All geological surface dough sheets that do not access are traveled through until all geological surface dough sheets are collectively labeled as having accessed, then allly
The generation of matter curved surface closed shape finishes.
10. according to the method described in claim 2, it is characterized in that, the contact class of crosspoint is judged according to the directed distance
Type and contact relation further comprise：
The directed distance on such as three vertex is with just or with bearing, then two triangles are non-intersect, and two triangles do not have crosspoint；
Directed distance such as three vertex is all zero, then two triangles are coplanar；
Directed distance such as two vertex is zero, and the directed distance on another vertex is just or is negative, then a triangle exists
A line another triangle institute in the plane, in the phase for planar judging two triangles of another triangle
Friendship condition；
Directed distance such as a vertex is zero, and the directed distance on another two vertex is the same as just or the same as negative, a then triangle presence
One vertex another triangle institute in the plane, planar judge the vertex and another in another triangle
The intersecting situation of a triangle；
Directed distance such as a vertex is zero, and the directed distance on another two vertex is negative for one positive one, then a triangle exists
One vertex another triangle institute in the plane, while the directed distance be zero vertex opposite side and another triangle
The place plane of shape judges the intersecting of the intersection point and another triangle in the plane there are intersection point, in another triangle
Situation；
If the directed distance on a vertex is just, the directed distance on another two vertex is same negative, or the directed distance on a vertex
It is negative, the directed distance on another two vertex is with just, then there are two sides and the place plane of another triangle for a triangle
It is intersecting, side where taking two vertex of directed distance contrary sign, the intersection point of plane where calculating the side and another triangle, and
The intersecting situation for judging the intersection point and another triangle in the plane of another triangle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610898488.3A CN107958483B (en) | 2016-10-14 | 2016-10-14 | Stratum modeling method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610898488.3A CN107958483B (en) | 2016-10-14 | 2016-10-14 | Stratum modeling method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107958483A true CN107958483A (en) | 2018-04-24 |
| CN107958483B CN107958483B (en) | 2021-05-25 |
Family
ID=61953395
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610898488.3A Active CN107958483B (en) | 2016-10-14 | 2016-10-14 | Stratum modeling method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107958483B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112836262A (en) * | 2019-11-22 | 2021-05-25 | 久瓴(上海)智能科技有限公司 | Method for acquiring intersection state between planes, model generation method and related product |
| CN113689557A (en) * | 2021-04-27 | 2021-11-23 | 电子科技大学 | Theoretical three-dimensional construction model construction method based on wire frame model |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103413297A (en) * | 2013-07-16 | 2013-11-27 | 南京师范大学 | Cutting method based on integrated three-dimensional GIS model |
| CN105844710A (en) * | 2016-04-11 | 2016-08-10 | 北京石油化工学院 | Data detection method for use in the process of geologic body meshing |
| CN105913491A (en) * | 2016-03-29 | 2016-08-31 | 浙江大学城市学院 | Gridding data processing method for three-dimensional geological surface model |
| CN105931297A (en) * | 2016-03-29 | 2016-09-07 | 浙江大学 | Data processing method applied to three-dimensional geological surface model |
-
2016
- 2016-10-14 CN CN201610898488.3A patent/CN107958483B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103413297A (en) * | 2013-07-16 | 2013-11-27 | 南京师范大学 | Cutting method based on integrated three-dimensional GIS model |
| CN105913491A (en) * | 2016-03-29 | 2016-08-31 | 浙江大学城市学院 | Gridding data processing method for three-dimensional geological surface model |
| CN105931297A (en) * | 2016-03-29 | 2016-09-07 | 浙江大学 | Data processing method applied to three-dimensional geological surface model |
| CN105844710A (en) * | 2016-04-11 | 2016-08-10 | 北京石油化工学院 | Data detection method for use in the process of geologic body meshing |
Non-Patent Citations (2)
| Title |
|---|
| 安聪荣 等: "基于层面结构的地质块体拓扑关系的自动建模", 《测绘学报》 * |
| 魏嘉 等: "三维地质构造建模技术研究", 《石油勘探》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112836262A (en) * | 2019-11-22 | 2021-05-25 | 久瓴(上海)智能科技有限公司 | Method for acquiring intersection state between planes, model generation method and related product |
| CN113689557A (en) * | 2021-04-27 | 2021-11-23 | 电子科技大学 | Theoretical three-dimensional construction model construction method based on wire frame model |
| CN113689557B (en) * | 2021-04-27 | 2023-03-14 | 电子科技大学 | Theoretical three-dimensional construction model construction method based on wire frame model |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107958483B (en) | 2021-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105825550B (en) | Take the complex three-dimensional building model cutting modeling method of consistency into account | |
| CN109815604B (en) | BIM indoor space connection diagram construction method based on topological relation of building elements | |
| CN101582173B (en) | Block model building method for complex geological structure | |
| US6907392B2 (en) | Method of generating a hybrid grid allowing modelling of a heterogeneous formation crossed by one or more wells | |
| CN106127857B (en) | The on-board LiDAR data modeling method of integrated data driving and model-driven | |
| CN105022865B (en) | A kind of fuel tanker inner surface model extraction method based on STL model boolean operations | |
| Chiang et al. | Simple and optimal output-sensitive construction of contour trees using monotone paths | |
| CN103413297A (en) | Cutting method based on integrated three-dimensional GIS model | |
| CN106981093A (en) | A kind of three-dimensional formation parallel modeling method of subregion constraint coupling | |
| CN107423518B (en) | Comprehensive pipe gallery full-line BIM collaborative design method combining Civil3D and Revit | |
| CN102495427B (en) | Interface perception ray tracing method based on implicit model expression | |
| CN104635262B (en) | A kind of positive reversed fault isopleth automatic generation method based on extended rectangular grid | |
| Tang et al. | An application-driven LOD modeling paradigm for 3D building models | |
| CN102609982B (en) | Topology discovery method of space geological data based on unstructured mode | |
| Qian et al. | Automatic unstructured all-hexahedral mesh generation from B-Reps for non-manifold CAD assemblies | |
| CN105205864B (en) | Geologic structure face threedimensional model method for automatic modeling and system based on multi-source data | |
| CN102279980A (en) | Geological exploration ore body three-dimensional modeling method and device thereof | |
| CN102314711A (en) | Three-dimensional visualization method and device for mineral resource evaluation information | |
| CN104318618B (en) | Three-dimensional sectioning method of generalized tri-prism spatial data model | |
| Xie et al. | A semantics-constrained profiling approach to complex 3D city models | |
| Ying et al. | Building 3D cadastral system based on 2D survey plans with SketchUp | |
| CN106447777A (en) | Three-dimensional topological relation expressing and mapping achieved under support of Boolean operation | |
| CN112150582B (en) | Multi-modal data-oriented geological profile approximate expression method | |
| CN107958483A (en) | A kind of Stratum Modeling | |
| Li et al. | 3D geological implicit modeling method of regular voxel splitting based on layered interpolation data |
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 |