CN105321204B - A kind of three-dimensional geological grid body Model construction method - Google Patents
A kind of three-dimensional geological grid body Model construction method Download PDFInfo
- Publication number
- CN105321204B CN105321204B CN201410796339.7A CN201410796339A CN105321204B CN 105321204 B CN105321204 B CN 105321204B CN 201410796339 A CN201410796339 A CN 201410796339A CN 105321204 B CN105321204 B CN 105321204B
- Authority
- CN
- China
- Prior art keywords
- point
- cell cube
- calculating
- geological
- line
- 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
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
The invention discloses a kind of three-dimensional geological grid body Model construction method, belongs to structure three-dimensional geological body Model field.A kind of three-dimensional geological grid body Model construction method of the present invention, idiographic flow are：Determine that the standard engineering planar range → standard rectangular region G rasterizings → calculating standard rectangular region G earth's surface triangulation network and the ground end triangulation network → calculating standard rectangular region G unit grids calculate line and merge cell cube according to the calculating information and positional information of cell cube with grid lines → equidistant partition boost line with calculating line → calculating discrete point geological property → foundation calculating point is classified with corresponding relation between auxiliary magnet → structure three-dimensional geological grid cell body → cell cube →, totally ten steps, by spatial spreading into series of points, and by calculating the geological property of point, then relevance between points is analyzed, and then merge space geometric from where discrete point, so as to calculate the overall geology distribution situation in space, precision height can be realized, limited space is small, efficiency high, the strong purpose of application.
Description
Technical field
The present invention relates to structure three-dimensional geological body Model field, is specifically related to a kind of three-dimensional geological grid body Model structure
Method.
Background technology
Data acquisition is mainly using the collection of discrete scatterplot and single-point line string collection two in existing three-dimensional geologic structure
Kind pattern.Discrete scatterplot collection is exactly directly to carry out discrete single-point sampling, and the connection not gathered between sampled point and sampled point is closed
It is information, such as based on topographic map collection discrete elevation point, field acquisition discrete elevation point, Laser Radar Scanning collection.Single-point
The collection of line string is exactly to be gathered and organized data in a manner of continuous coordinate string, in the collection, digital photogrammetry such as contour
Collection of property line, etc..The common feature of above-mentioned sampling is the spatial positional information for only gathering sampled point, does not gather sampled point work
With the relief form attribute information in region.In fact, each sampled point of topographical surface, it is not only a space orientation number
According to, it is often more important that its form properties information, whether such as Feature line, it is smooth company just to have in tangential direction
Whether continuous curve or broken line, smooth treatment, a series of attributes such as smooth degree is how many can be carried out in vertical section direction.
In theory, the relief form of each sampled point subrange can be reconstructed by combining its periphery sampled point, but in reality
In the application of border, because the complexity of landform is far beyond the analog capability of existing mathematics and computer science, therefore work as sampled point
Local form attribute information lack when, especially the form properties loss of learning of feature terrain when, in practice be difficult high-fidelity weight
Its terrain surface of structure.In fact this is also to cause existing three-dimensional geological body Model mesorelief expression the distortion of critical regions property to be present to show
As the main reason for one of.
The main method of current three-dimensional geologic structure is to be entered with Arcgis, Surpac, QvCoalMine etc. for platform
Row fitting builds three-dimensional geologic mould with manually repairing the method structure three-dimensional geological body Model that woollen yarn knitting closes in these traditional methods
Type is, it is necessary to substantial amounts of artificial, and its accuracy is subject to applied platform technology, in all big discounting of precision, efficiency etc.
Button, so, need one kind badly and be based on a wide range of geologic data, the three-dimensional geological object model skill using a wide range of geological form as background
Art.
2009, Paradigm Geophysical Corp. of the U.S. proposed it is the spatial spread by column, and then instead
The preparation method for drilling the structure three-dimensional geological model of the distribution of space geologic body, but it calculates in space optional position and limited to, mould
Type limited precision.
Chinese Patent Application No.：201310238811, publication date：In September, 2013 patent document of 4 days discloses one kind three
Tie up horizontal geologic profile figure preparation method and application thereof.This method is with 1. according to energy and mines geological and engineering geology correlation technique
Code requirement, obtain field geology firsthand information, include without being confined to drilling geologic information, drill site three dimensional local information,
Geology geophysical prospecting informations, and the 10 steps step such as associated chamber is interior, outdoor test data are completed.The invention is used for：Mapping of Engineering Geology
Energy and mines geological drawing of the present invention；In engineering geology and energy and mines geological in geological mapping geological mapping achievement expression；
Include without the expression for being confined to related investigation result in engineering investigation, design and analysis, construction, monitoring, operation link；Including
And it is not limited to energy minerals planning, exploration, design, monitoring, exploitation, operation, the expression for repairing related investigation result in link.
The patent document determines that block diagram carries out the distribution structure of inverting space geologic body dimensionally by existing original documents and materials
Matter model instructs geological mapping to work, due to firsthand information by the age, exploring equipment was horizontal at that time, technical staff specialty water
The limitation of flat, geology change etc., the reliability of the precision of constructed three-dimensional geological model allow people to suspect unavoidably.
The content of the invention
1. technical problems to be solved
For being currently based on a wide range of geologic data, the three-dimensional geological object model using a wide range of geological form as background, essence
The problem of degree, efficiency limitation, of the invention provides a kind of three-dimensional geological grid body Model construction method.The construction method is empty
Between be separated into a series of point, by calculate point geological property, then analyze relevance between points, so merge from
Space geometric from where scatterplot, so as to calculate the overall geology distribution situation in space, it is possible to achieve precision is high, limited space
It is small, efficiency high, the strong purpose of application.
2. technical scheme
To reach above-mentioned purpose, technical scheme provided by the invention is：
A kind of three-dimensional geological grid body Model construction method, idiographic flow are as follows：
Step 1, determine standard engineering planar range：Known standard engineering planar range is rectangular area [(x1, x2), (y1,
y2)], the unit of grid body be：u[ux, uy, uz], obtain standard rectangular region G [(xmin, xmax), (ymin, ymax)]；
Step 2, standard rectangular region G rasterizings：With point O (0,0) for origin, spacing is [ux, uy], to standard rectangular area
Domain G carries out rasterizing, and standard rectangular region G is divided into nx×nyIndividual continuous grids unit, nxFor
ForGrid cell centers point is taken as grid computing point, then corresponding four grid nodes of each grid cell and
One grid computing point；Then standard rectangular region G is expressed as (nx+1)×(ny+ 1) individual grid node and nx×nyIndividual grid computing
The point set P of point composition；Will point pi[(spatial value), point mark] enters point set P by ranks ordered storage；pi∈ P, point mark are
Distinguish piIt is grid computing point or grid node；
Step 3, the earth's surface triangulation network and the ground end triangulation network for calculating standard rectangular region G：
According to drilling top coordinate structure earth's surface triangulation network TRIT；According to drilling bottom coordinate structure ground end triangle T RIB, and
The triangulation network subdivision of each big layer；
Step 4, the calculating line and grid lines for calculating standard rectangular region G unit grids：
Calculate point piIn triangulation network TRIT, TRIBOn subpoint：By dividing equally cloth standard rectangular area in point set P a little
In the G of domain, so in triangulation network TRIT, TRIBOn can find the point corresponding to subpoint, two subpoints for connecting the point obtain line
li；Such as fruit dot piFor grid computing point, then liTo calculate line；Such as fruit dot piFor grid node, then liFor boost line；According to triangle
Identified plane equation, calculate point piIn triangulation network TRIT、TRIBOn subpoint Z values be respectively Zt、Zb；
Then Current standards works scope is expressed as nx×nyBar calculates line and (nx+1)×(ny+ 1) bar boost line, by grid
Node and grid computing point correspondence understand that every calculates corresponding four boost lines of line；
Step 5, equidistant partition boost line are with calculating line：
WithOn the basis of point, with uzGrown for unit, averagely divide boost line, division points are referred to as auxiliary magnet；
The point on the basis of Z=0, with uzGrown for unit, average division calculates line, and division points are referred to as to calculate point；
Understand that calculates point eight auxiliary magnets of corresponding surrounding with calculating line mechanism by equidistant partition boost line, it is corresponding auxiliary
Help the cell cube being a little built into where the calculating point；
Step 6, discrete point geological property is calculated, it is determined that the cell cube where each discrete point, and the ground of confirmation unit body
Matter attribute；
The borehole space positional information collected, the distribution of strata situation of drilling, geologic setting parameter, geology physical prospecting letter
Breath, the comprehensive layering data of drilling, space of discrete points position are as Given information, by fuzzy, probability, random process uncertain
Property method, to discrete point carry out metric calculation, obtain probability distribution situation of every kind of geological property in the discrete point；Take probability set
The geological property for closing the maximum in L is the geological property of the discrete point；Take the probability in set L impartial, be considered as the discrete point
Boundary in geologic body；Take the geological property that any one geological property is the discrete point；
Step 7, establish the corresponding relation calculated between point and auxiliary magnet；
A：It is numbered to calculating line with boost line, coding rule：[row number, line number], is designated as [I, J]
B：The corresponding relation for calculating line and boost line is established, numbering is the numbering for calculating boost line corresponding to line of [I, J]
For：[I, J], [I+1, J], [I+1, J+1], [I, J+1]；
C：Successively to auxiliary magnet in the boost line of every, it is numbered from top to bottom with natural number order；Successively to every
Calculating line on calculate a little, be numbered from top to bottom with natural number order；
D：From top to bottom, auxiliary magnet corresponding to each calculating point is found out successively；Its corresponding principle is：Auxiliary magnet Z values point
Cloth is positioned at [calculating point Z values-uz, calculate point Z values+uz] in section；
Step 8, structure three-dimensional geological grid cell body：Current standards works scope is expressed as to include geology category by multiple
Property, joint account information, the cell cube of positional information assemble T three-dimensional geological model, wherein cell cube Ti,j,k∈T；
One calculates point eight auxiliary magnets of corresponding surrounding, and corresponding auxiliary magnet is built into the cell cube where the calculating point,
Three-dimensional geological grid cell body is built, calculates the geological property of geological property cell cube where the calculating point a little；
Step 9, Column Row, according to computation sequence and cell cube positional information, cell cube is classified；Pass through construction unit
Body, each calculates cell cube corresponding to calculating point on line, forms a cell cube row, line number and columns one with calculating line
Cause；
Step 10, according to calculating information, the positional information of cell cube cell cube is merged；Merge geological property identical list
First body, and obtain the distribution entity of every kind of geology；
Preferably, in described step 5, auxiliary magnet Z Distribution values are Z=[n2,{n},n1]×uz, wherein N is section [n2, n1] between natural number set；Calculating point Z Distribution values is：Z
=[m2,{m},m1]×uz, whereinM is section [m2, m1] between natural number set.
Preferably, in described step 6, the geological property calculating process of discrete point is：
A：If sharing m kind geological properties in borehole data, then its probability distribution is set L：{ l1, l2, l3 ..., lm }, its
In
B：The maximum in set L is taken, writes down and is designated as i, then Li>Lj, wherein j=1,2 ..., i-1, i+1 ... m；Then i-th
Kind geological property is the geological property of the discrete point；
C：The maximum in set L is taken, writes down and is designated as i1, i2 ..., ik common k, and Li1=Li2=...=Lik>Lj,
Wherein j ≠ i1, i2 ..., ik；Then the discrete point is in the boundary of geologic body, and it is i-th 1, i2 to remember the geological property ..., ik
In any one geological property.
Preferably, three-dimensional geological grid body Model construction step is in described step 8：
A：Understand that calculates point eight auxiliary magnets of corresponding surrounding with calculating line mechanism by equidistant partition boost line, it is corresponding
Auxiliary magnet is built into the cell cube where the calculating point；
B：The face of the cell cube is numbered, counts the type in all faces of the cell cube, by face be divided into vertical plane with
Horizontal plane, the point numbering that each bread contains is counted, subscript of the index face in mesh generation, is the unit by gained statistical information
Body joint account information；
C：Remember the geological property of geological property cell cube where the calculating point of the calculating point；
Then current works scope is expressed as by multiple units for including geological property, joint account information, positional information
Body set T three-dimensional geological grid body, wherein cell cube Ti,j,k∈T。
Preferably, cell cube sorting technique is in described step 9：
With cell cube Ti,j,kExemplified by, its positional information [i, j, k], i be column ordinal number i ∈ 1,2 ... .., nx, j
By the ordinal number j ∈ that are expert at 1,2 ... .., ny, K is vertical direction ordinal number；
Class1：Do not intersect with other any cell cubes；
Type 2：Intersect with a cell cube；
Type 3：Intersect with two cell cubes；
Type 4：Intersect simultaneously with three cell cubes.
Preferably, the method for cell cube merging is in described step 10：
A：The cell cube not intersected with other any cell cubes merges：Determine cell cube T0,0,0To originate geologic body, and the ground
The geological property of plastid and cell cube T0,0,0Unanimously；
B：The cell cube intersected with a cell cube merges：Such as the cell cube is different from the geological property of intersecting cell cube,
Nonjoinder；Conversely, find out the place geologic body of intersecting cell cube, deleted from geologic body the cell cube intersection surface and its
Place geologic data, and remaining 5 face is added into the place geologic body of intersecting cell cube；The cell cube is with intersecting cell cube
It is incorporated into same geologic body；
C：The cell cube intersected with two cell cubes merges：Such as cell cube hands over the geological property of cell cube not with two-phase
Together, equal nonjoinder；It is as identical in only intersected cell cube geological property with wherein one, then merged according to B modes；Such as three lists
First body geological property is identical, and cell cube is added into geologic body where adjacent cells body, and need to merge intersecting cell cube place
Geologic body；
D：The cell cube intersected with three cell cubes merges：Such as cell cube hands over the geological property of cell cube not with two-phase
Together, equal nonjoinder；It is as identical in only intersected cell cube geological property with wherein one, then merged according to C modes；Such as four lists
First body geological property is identical, merges four cell cube Cheng Xin cell cube.
3. beneficial effect
Compared to immediate prior art, the present invention is advantageous in that：
(1) a kind of three-dimensional geological grid body Model construction method of the invention, three-dimensional geological grid body Model architecture
Three-dimensional geological model, straightforward discloses abundant stratum space distribution information, and the original geology letter of direct correlation
Breath, makes technical staff quickly form complete, accurate geology three-dimensional spatial distribution form concept, is provided more for applied geology information
For effective, comprehensive, clearly technical support；
(2) a kind of three-dimensional geological grid body Model construction method of the invention, the whole engineering geology of concentrated expression are three-dimensional
Stratum space distribution rule and variation tendency；
(3) a kind of three-dimensional geological grid body Model construction method of the invention, engineering geology and energy ore deposit are accurately reflected
The geological condition of place of production matter, can be the related design and analysis of Geotechnical Engineering, the construction design related to energy minerals, exploitation,
Repair and directly perceived, accurate, comprehensive technical support is provided, greatly improve engineering efficiency and reduce engineering cost.
Brief description of the drawings
Fig. 1 is standard area G grid nodes and grid computing point distribution map；
Fig. 2 is works scope internal drilling plane distribution schematic diagram；
Fig. 3 is calculating point and auxiliary magnet plane distribution schematic diagram；
Fig. 4 is three-dimensional geological model schematic diagram；
Fig. 5 is that three-dimensional geological cell cube lists intention.
Embodiment
Explanation is further explained to technical scheme with reference to embodiment and accompanying drawing.
Embodiment 1
First related notion is illustrated below：
Geologic body：Certain geology spatial distribution geometric shape, in particular to similar stair-stepping geometric shape.
Geologic body specifying information has geologic body where cell cube numbering, geologic body numbering and cell cube containing cell cube to compile
Number.
A kind of three-dimensional geological grid body Model construction method of the present embodiment, idiographic flow are as follows：
Step 1, determine standard engineering planar range：Known standard engineering planar range is rectangular area [(x1, x2), (y1,
y2)], the unit of grid body be：u[ux, uy, uz], obtain standard rectangular region G [(xmin, xmax), (ymin, ymax)]；Then G is calculated
Method is：
Step 2, standard rectangular region G rasterizings, ensure the geology distribution uniformity of any two adjacent area：
With point O (0,0) for origin, spacing is [ux, uy], rasterizing is carried out to standard rectangular region G, by standard rectangular area
Domain G is divided into nx×nyIndividual continuous grids unit, nxFornyForTake grid cell centers
Point is used as grid computing point, as shown in figure 1, each corresponding four grid nodes of grid cell and a grid computing point；Then mark
Quasi- rectangular area G is expressed as (nx+1)×(ny+ 1) individual grid node and nx×nyThe point set P of individual grid computing point composition；Will point pi
[(spatial value), point mark] enters point set P by ranks ordered storage；pi∈ P, point mark are to discriminate between piBe grid computing point or
It is grid node；
Step 3, the earth's surface triangulation network and the ground end triangulation network for calculating standard rectangular region G：
According to drilling top coordinate structure earth's surface triangulation network TRIT：
As shown in Fig. 2 the node in figure in the range of black wire frame is the drilling plane distribution schematic diagram of certain engineering, by drilling
Top and bottom coordinate, respectively interpolation go out tetra- point [x of standard rectangular area Gmin,ymin]、[xmin,ymax]、[xmax,ymin]、
[xmax,ymax] top and bottom elevation, be designated as [zt(xmin,ymin),zb(xmin,ymin)]、[zt(xmin,ymax),zb(xmin,ymax)]、
[zt(xmax,ymin),zb(xmax,ymin)]、[ztxmax,ymax,zbxmax,ymax]；Note drilling apex coordinate and point [xmin,ymin,
zt(xmin,ymin)]、[xmin,ymax,zt(xmin,ymax)]、[xmax,ymin,zt(xmax,ymin)]、[xmax,ymax,ztxmax,ymax] it is set
ST；Calculated using Delaunay Triangulation and carry out the structure triangulation network：
A：A wide range of triangle is built, comprising all scatterplots in all set ST, big Atria summit is not wrapped
It is contained in set ST；
B：It is any from set ST to choose a bit, first Delaunay Triangulation is built, altogether comprising three triangles；
C：Still uncalculated point in set of computations ST successively, finds out the public portion for being located at adjacent two triangle circumscribed circle
Two adjacent triangles in being divided to, and the common edge of adjacent triangle is deleted, and 2 points of the point and residue are connected, formed new
Delaunay Triangulation；
D：Delete the triangle containing a wide range of triangular apex, you can obtain Delaunay Triangulation to the end；
Similarly according to above method, structure ground end triangle T RIB, and the triangulation network subdivision of each big layer；
Step 4, the calculating line and grid lines for calculating standard rectangular region G unit grids：
Calculate point piIn triangulation network TRIT, TRIBOn subpoint：By dividing equally cloth standard rectangular area in point set P a little
In the G of domain, so in triangulation network TRIT, TRIBOn can find the point corresponding to subpoint, two subpoints for connecting the point obtain line
li；Such as fruit dot piFor grid computing point, then liTo calculate line；Such as fruit dot piFor grid node, then liFor boost line；According to triangle
Identified plane equation, calculate point piIn triangulation network TRIT、TRIBOn subpoint Z values be respectively Zt、Zb；
Then Current standards works scope is expressed as nx×nyBar calculates line and (nx+1)×(ny+ 1) bar boost line, by grid
Node and grid computing point correspondence understand that every calculates corresponding four boost lines of line；
Step 5, equidistant partition boost line are with calculating line：
With Z=On the basis of point, with uzGrown for unit, averagely divide boost line, division points are referred to as auxiliary magnet, then aid in
Point Z Distribution values are Z=[n2,{n},n1]×uz, whereinN is section [n2, n1]
Between natural number set；
The point on the basis of Z=0, with uzGrown for unit, average division calculates line, and division points are referred to as to calculate point, then calculate point Z
Distribution value is：Z=[m2,{m},m1]×uz, whereinM is section [m2, m1] between natural manifold
Close；
As shown in figure 3, understand that calculates corresponding eight auxiliary of surrounding with calculating line mechanism by equidistant partition boost line
Point, corresponding auxiliary magnet are built into the cell cube where the calculating point；
Step 6, discrete point geological property is calculated, it is determined that the cell cube where each discrete point, and the ground of confirmation unit body
Matter attribute；
The borehole space positional information collected, the distribution of strata situation of drilling, geologic setting parameter, geology physical prospecting letter
Breath, the comprehensive layering data of drilling, space of discrete points position are as Given information, by fuzzy, probability, random process uncertain
Property method, to discrete point carry out metric calculation, obtain probability distribution situation of every kind of geological property in the discrete point；
The geological property calculating process of discrete point is：
A：If sharing m kind geological properties in borehole data, then its probability distribution is set L：{ l1, l2, l3 ..., lm }, its
In
B：The maximum in set L is taken, writes down and is designated as i, then Li>Lj, wherein j=1,2 ..., i-1, i+1 ... m；Then i-th
Kind geological property is the geological property of the discrete point；
C：The maximum in set L is taken, writes down and is designated as i1, i2 ..., ik common k, and Li1=Li2=...=Lik>Lj,
Wherein j ≠ i1, i2 ..., ik；Then the discrete point is in the boundary of geologic body, and it is i-th 1, i2 to remember the geological property ..., ik
In any one geological property.
Step 7, establish the corresponding relation calculated between point and auxiliary magnet；
A：It is numbered to calculating line with boost line, coding rule：[row number, line number], is designated as [I, J]
B：The corresponding relation for calculating line and boost line is established, numbering is the numbering for calculating boost line corresponding to line of [I, J]
For：[I, J], [I+1, J], [I+1, J+1], [I, J+1]；
C：Successively to auxiliary magnet in the boost line of every, it is numbered from top to bottom with natural number order；Successively to every
Calculating line on calculate a little, be numbered from top to bottom with natural number order；
D：From top to bottom, auxiliary magnet corresponding to each calculating point is found out successively；Its corresponding principle is：Auxiliary magnet Z values point
Cloth is positioned at [calculating point Z values-uz, calculate point Z values+uz] in section；
Step 8, structure three-dimensional geological grid cell body：
A：One calculates point eight auxiliary magnets of corresponding surrounding, and corresponding auxiliary magnet is built into the unit where the calculating point
Body；
B：The face of the cell cube is numbered, counts the type in all faces of the cell cube, by face be divided into vertical plane with
Horizontal plane, the point numbering that each bread contains is counted, subscript of the index face in mesh generation, is the unit by gained statistical information
Body joint account information；
C：Remember the geological property of geological property cell cube where the calculating point of the calculating point；
As shown in figure 4, then current works scope is expressed as by multiple comprising geological property, joint account information, position
Three-dimensional geological the grid body T, wherein cell cube T of the cell cube set of informationi,j,k∈T；
Step 9, Column Row, according to computation sequence and cell cube positional information, cell cube is classified；Pass through construction unit
Body, each calculates cell cube corresponding to calculating point on line and forms a cell cube row, consistent with columns with the line number for calculating line,
As shown in Figure 5；
With cell cube Ti,j,kExemplified by, its positional information [i, j, k], i be column ordinal number i ∈ 1,2 ... .., nx, j
By the ordinal number j ∈ that are expert at 1,2 ... .., ny, K is vertical direction ordinal number；
Class1：Do not intersect with other any cell cubes；
Type 2：Intersect with a cell cube；
Type 3：Intersect with two cell cubes；
Type 4：Intersect simultaneously with three cell cubes；
Step 10, according to calculating information, the positional information of cell cube cell cube is merged；Merge geological property identical list
First body, and obtain the distribution entity of every kind of geology：
The cell cube not intersected with other any cell cubes is merged：Determine cell cube T0,0,0To originate geologic body, and the ground
The geological property of plastid and cell cube T0,0,0Unanimously；
Embodiment 2
A kind of three-dimensional geological grid body Model construction method of the present embodiment, substantially the same manner as Example 1, difference exists
It is different in the cell cube type of merging：
The cell cube intersected with a cell cube is merged：Such as the cell cube is different from the geological property of intersecting cell cube,
Nonjoinder；Conversely, find out the place geologic body of intersecting cell cube, deleted from geologic body the cell cube intersection surface and its
Place geologic data, and remaining 5 face is added into the place geologic body of intersecting cell cube；The cell cube is with intersecting cell cube
It is incorporated into same geologic body.
Embodiment 3
A kind of three-dimensional geological grid body Model construction method of the present embodiment, substantially the same manner as Example 1, difference exists
It is different in the cell cube type of merging：
The cell cube intersected with two cell cubes is merged：Such as cell cube hands over the geological property of cell cube not with two-phase
Together, equal nonjoinder；It is as identical in only intersected cell cube geological property with wherein one, then merged according to the mode of embodiment 2；
Such as three cell cube geological properties are identical, cell cube added into geologic body where adjacent cells body, and need to merge intersecting unit
Geologic body where body；
Embodiment 4
A kind of three-dimensional geological grid body Model construction method of the present embodiment, substantially the same manner as Example 1, difference exists
It is different in the cell cube type of merging：
The cell cube intersected with three cell cubes is merged：Such as cell cube hands over the geological property of cell cube not with two-phase
Together, equal nonjoinder；It is as identical in only intersected cell cube geological property with wherein one, then merged according to the mode of embodiment 3；
Such as four cell cube geological properties are identical, merge four cell cube Cheng Xin cell cube.
Schematically the invention and embodiments thereof are described above, actual technical scheme is not limited to
In this.So if one of ordinary skill in the art is enlightened by it, in the case where not departing from this creation objective, without wound
The protection domain designed the embodiment similar to the technical scheme, this patent all should be belonged to for the property made.
Claims (5)
1. a kind of three-dimensional geological grid body Model construction method, method flow are as follows：
Step 1, determine standard engineering planar range：Known standard engineering planar range is rectangular area [(x1, x2), (y1, y2)]、
The unit of grid body is：u[ux, uy, uz], obtain standard rectangular region G [(xmin, xmax), (ymin, ymax)]；
Step 2, standard rectangular region G rasterizings：With point O (0,0) for origin, spacing is [ux, uy], standard rectangular region G is entered
Row rasterizing, standard rectangular region G is divided into nx×nyIndividual continuous grids unit, nxFornyForGrid cell centers point is taken as grid computing point, then corresponding four grid nodes and one of each grid cell
Individual grid computing point；Then standard rectangular region G is expressed as (nx+1)×(ny+ 1) individual grid node and nx×nyIndividual grid computing point
The point set P of composition；Will point pi[(spatial value), point mark] enters point set P by ranks ordered storage；pi∈ P, point mark is area
Divide piIt is grid computing point or grid node；
Step 3, the earth's surface triangulation network and the ground end triangulation network for calculating standard rectangular region G：
According to drilling top coordinate structure earth's surface triangulation network TRIT；According to drilling bottom coordinate structure ground end triangle T RIB, and each
The triangulation network subdivision of big layer；
Step 4, the calculating line and grid lines for calculating standard rectangular region G unit grids：
Calculate point piIn triangulation network TRIT, TRIBOn subpoint：By dividing equally cloth standard rectangular region G in point set P a little
It is interior, so in triangulation network TRIT, TRIBOn can find the point corresponding to subpoint, two subpoints for connecting the point obtain line li；
Such as fruit dot piFor grid computing point, then liTo calculate line；Such as fruit dot piFor grid node, then liFor boost line；According to triangle institute
The plane equation of determination, calculate point piIn triangulation network TRIT、TRIBOn subpoint Z values be respectively Zt、Zb；
Then Current standards works scope is expressed as nx × ny bars and calculates line and (nx+1)×(ny+ 1) bar boost line, by grid node
Understand that every calculates corresponding four boost lines of line with grid computing point correspondence；
Step 5, equidistant partition boost line are with calculating line：
WithOn the basis of point, with uzGrown for unit, averagely divide boost line, division points are referred to as auxiliary magnet；
The point on the basis of Z=0, with uzGrown for unit, average division calculates line, and division points are referred to as to calculate point；
Understand that one calculates point eight auxiliary magnets of corresponding surrounding, corresponding auxiliary magnet with calculating line mechanism by equidistant partition boost line
The cell cube being built into where the calculating point；
Step 6, discrete point geological property is calculated, it is determined that the cell cube where each discrete point, and the geology category of confirmation unit body
Property：The borehole space positional information collected, the distribution of strata situation of drilling, geologic setting parameter, geology geophysical prospecting informations, brill
Hole synthesis layering data, space of discrete points position pass through fuzzy, probability, random process uncertainty side as Given information
Method, metric calculation is carried out to discrete point, obtains probability distribution situation of every kind of geological property in the discrete point；Take Making by Probability Sets L
In maximum geological property be the discrete point geological property, when taking the probability equalization of the maximum in set L, be considered as
The discrete point is in the boundary of geologic body, takes the geological property that any one geological property is the discrete point；
Step 7, establish the corresponding relation calculated between point and auxiliary magnet；
A：It is numbered to calculating line with boost line, coding rule：[row number, line number], is designated as [I, J]
B：The corresponding relation for calculating line and boost line is established, numbering is that calculating numbering of boost line corresponding to line of [I, J] is：
[I, J], [I+1, J], [I+1, J+1], [I, J+1]；
C：Successively to auxiliary magnet in the boost line of every, it is numbered from top to bottom with natural number order；Meter to every successively
Calculate and calculated a little on line, is numbered from top to bottom with natural number order；
D：From top to bottom, auxiliary magnet corresponding to each calculating point is found out successively；Its corresponding principle is：Auxiliary magnet Z Distribution values position
In [calculating point Z values-uz, calculate point Z values+uz] in section；
Step 8, structure three-dimensional geological grid cell body, step are：
A：One calculates point eight auxiliary magnets of corresponding surrounding, and corresponding auxiliary magnet is built into the cell cube where the calculating point；
B：The face of the cell cube is numbered, counts the type in all faces of the cell cube, face is divided into vertical plane and level
Face, the point numbering that each bread contains is counted, subscript of the index face in mesh generation, is that the cell cube is closed by gained statistical information
And calculate information；
C：Remember the geological property of geological property cell cube where the calculating point of the calculating point；
Then Current standards works scope is expressed as by multiple units for including geological property, joint account information, positional information
Body assembles T three-dimensional geological model, wherein cell cube Ti,j,k∈T；
Step 9, Column Row, according to computation sequence and cell cube positional information, cell cube is classified；By construction unit body,
Each calculates cell cube corresponding to calculating point on line and forms a cell cube row, consistent with columns with the line number for calculating line；
Step 10, according to calculating information, the positional information of cell cube cell cube is merged；Merge geological property identical unit
Body, and obtain the distribution entity of every kind of geology.
A kind of 2. three-dimensional geological grid body Model construction method according to claim 1, it is characterised in that：Described step
In 5, auxiliary magnet Z Distribution values are Z=[n2,{n},n1]×uz, wherein N is section [n2,
n1] between natural number set；Calculating point Z Distribution values is：Z=[m2,{m},m1]×uz, whereinM is
Section [m2, m1] between natural number set.
A kind of 3. three-dimensional geological grid body Model construction method according to claim 1 or 2, it is characterised in that：Described
In step 6, the geological property calculating process of discrete point is：
A：If sharing m kind geological properties in borehole data, then its probability distribution is set L：{ l1, l2, l3 ..., lm }, wherein
B：The maximum in set L is taken, writes down and is designated as i, then Li>Lj, wherein j=1,2 ..., i-1, i+1 ... m；Then i-th kind of ground
Matter attribute is the geological property of the discrete point；
C：The maximum in set L is taken, writes down and is designated as i1, i2 ..., ik common k, and Li1=Li2=...=Lik>Lj, wherein j
≠i1,i2,…,ik；Then the discrete point is in the boundary of geologic body, and it is i-th 1, i2 to remember the geological property ..., in ik
Any one geological property.
A kind of 4. three-dimensional geological grid body Model construction method according to claim 3, it is characterised in that：Described step
Cell cube sorting technique is in 9：
With cell cube Ti,j,kExemplified by, its positional information [i, j, k], i be column ordinal number i ∈ 1,2 ... .., nx, j is institute
Be expert at ordinal number j ∈ 1,2 ... .., ny, k is vertical direction ordinal number；
Class1：Do not intersect with other any cell cubes；
Type 2：Intersect with a cell cube；
Type 3：Intersect with two cell cubes；
Type 4：Intersect simultaneously with three cell cubes.
A kind of 5. three-dimensional geological grid body Model construction method according to claim 4, it is characterised in that：Described step
The method of cell cube merging is in 10：
A：The cell cube not intersected with other any cell cubes merges：Determine cell cube T0,0,0To originate geologic body, and the geologic body
Geological property and cell cube T0,0,0Unanimously；
B：The cell cube intersected with a cell cube merges：Such as the cell cube is different from the geological property of intersecting cell cube, does not conform to
And；Conversely, finding out the place geologic body of intersecting cell cube, intersection surface and its place of the cell cube are deleted from geologic body
Geologic data, and remaining 5 face is added into the place geologic body of intersecting cell cube；The cell cube merges with intersecting cell cube
Into same geologic body；
C：The cell cube intersected with two cell cubes merges：Such as cell cube hands over the geological property of cell cube different from two-phase,
Equal nonjoinder；It is as identical in only intersected cell cube geological property with wherein one, then merged according to B modes；Such as Unit three
Body geological property is identical, and cell cube is added into geologic body where adjacent cells body, and where need to merging intersecting cell cube
Geologic body；
D：The cell cube intersected with three cell cubes merges：Such as cell cube hands over the geological property of cell cube different from two-phase,
Equal nonjoinder；It is as identical in only intersected cell cube geological property with wherein one, then merged according to C modes；Such as Unit four
Body geological property is identical, merges four cell cube Cheng Xin cell cube.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410796339.7A CN105321204B (en) | 2014-12-18 | 2014-12-18 | A kind of three-dimensional geological grid body Model construction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410796339.7A CN105321204B (en) | 2014-12-18 | 2014-12-18 | A kind of three-dimensional geological grid body Model construction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105321204A CN105321204A (en) | 2016-02-10 |
| CN105321204B true CN105321204B (en) | 2017-11-14 |
Family
ID=55248515
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410796339.7A Active CN105321204B (en) | 2014-12-18 | 2014-12-18 | A kind of three-dimensional geological grid body Model construction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105321204B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106920275B (en) * | 2017-01-24 | 2021-05-28 | 天衍智（北京）科技有限公司 | Complex attribute boundary three-dimensional vector iteration method and application system |
| CN108510134B (en) * | 2017-02-27 | 2022-02-01 | 中国石油化工股份有限公司 | Method and system for determining geological variable quantile value |
| CN107424210A (en) * | 2017-06-16 | 2017-12-01 | 中国电力科学研究院 | The processing method and processing device that mountain area mima type microrelief boundary profile smoothly transits |
| CN111414654B (en) * | 2020-03-16 | 2024-02-23 | 南京交通职业技术学院 | Method for retrieving continuous numerical model of division |
| CN112102484B (en) * | 2020-08-12 | 2022-11-25 | 特雷西能源科技(杭州)有限公司 | Geological model parameter field adjusting method and device |
| CN111986325B (en) * | 2020-08-25 | 2023-07-18 | 东北大学 | Implicit sectioning method for three-dimensional geologic body model |
| CN113393578B (en) * | 2021-06-02 | 2024-08-23 | 南京师范大学 | Method and device for constructing geological section three-dimensional model |
| CN117392335B (en) * | 2023-09-26 | 2024-08-20 | 深圳市地质环境研究院有限公司 | Conversion method from three-dimensional surface element geologic model to three-dimensional surface element geologic model |
| CN117541741B (en) * | 2024-01-10 | 2024-03-19 | 自然资源部第三地理信息制图院 | Construction method of three-dimensional model of geologic body and electronic equipment |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102194253A (en) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | Method for generating tetrahedron gridding for three-dimensional geological structure |
| CN103500468A (en) * | 2013-09-09 | 2014-01-08 | 电子科技大学 | Space curved surface reconstruction method based on boundary constraint |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101582173B (en) * | 2009-06-24 | 2012-07-11 | 中国石油天然气集团公司 | Block model building method for complex geological structure |
-
2014
- 2014-12-18 CN CN201410796339.7A patent/CN105321204B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102194253A (en) * | 2011-06-15 | 2011-09-21 | 北京航空航天大学 | Method for generating tetrahedron gridding for three-dimensional geological structure |
| CN103500468A (en) * | 2013-09-09 | 2014-01-08 | 电子科技大学 | Space curved surface reconstruction method based on boundary constraint |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105321204A (en) | 2016-02-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105321204B (en) | A kind of three-dimensional geological grid body Model construction method | |
| CN103279986B (en) | Three-dimensional horizontal geologic profile figure preparation method and application thereof | |
| CN103514630B (en) | Fault structure three-dimensional modeling method | |
| CN106023313B (en) | A kind of three-dimensional geological number exploration report generation method | |
| CN116152461B (en) | Geological modeling method, device, computer equipment and computer readable storage medium | |
| CN102609982B (en) | Topology discovery method of space geological data based on unstructured mode | |
| Dou et al. | 3D geological suitability evaluation for urban underground space development–A case study of Qianjiang Newtown in Hangzhou, Eastern China | |
| CN107944086B (en) | Rapid modeling method based on drilling data | |
| CN111339691A (en) | Intelligent geotechnical engineering parameter three-dimensional analysis and evaluation system and method based on voxler software | |
| CN102622526A (en) | Digital mine tunneling search prediction method | |
| CN112446956B (en) | Geological three-dimensional model construction method based on drilling data | |
| CN104809266A (en) | Spline based accurate predicating method for face coal seam occurrence condition | |
| CN109507749A (en) | A kind of heavy magnetic is from constraint 3-d inversion and joint interpretation method | |
| CN107886575A (en) | A kind of method that open-pit mine stope triangular mesh cuts coal seam quadrilateral mesh | |
| He et al. | Three-dimensional hydrogeological modeling method and application based on TIN-GTP-TEN | |
| Leoni et al. | GIS methodology to assess landslide susceptibility: application to a river catchment of Central Italy | |
| Zhang et al. | Spatial patterns and controlling factors of the evolution process of karst depressions in Guizhou province, China | |
| CN109949415B (en) | System and method for modeling three-dimensional earth surface and geological body model in topological consistency | |
| WO2018196214A1 (en) | Statistics system and statistics method for geographical influence on vernacular architectural form | |
| CN112489162A (en) | Large-range micro-unit coal seam geological prediction and profile drawing method | |
| CN104462649B (en) | A kind of automatic update method of ore body block segment model reserves | |
| Shen et al. | Three-dimensional modeling of loose layers based on stratum development law | |
| CN117830561A (en) | Gold ore prospecting method based on three-dimensional geological modeling | |
| CN106504319B (en) | Reservoir Three Dimensional Contrast map generalization method and device between well | |
| CN115831276A (en) | Space classification and partition method and system for soil materials |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |