CN105321204A - Three-dimensional geological raster body model construction method - Google Patents

Abstract

The present invention discloses a three-dimensional geological raster body model construction method, and belongs to the field of constructing three-dimensional geological body models. The three-dimensional geological raster body model construction method disclosed by the present invention has a particular procedure that comprises the following ten steps: determining a standard project plane range; rasterizing a standard rectangle region G; calculating an earth-surface triangulation net and an underground triangulation net of the standard rectangle region G; calculating a unit grid calculation line and a grid line of the standard rectangle region G; equidistantly dividing an auxiliary line and the calculation line; calculating geological attributes of discrete points; establishing a correspondence relationship between a calculation point and an auxiliary point; constructing three-dimensional geological raster unit bodies; classifying the unit bodies; and combining the unit bodies according to calculation information and position information of the unit bodies. A space is dispersed into a series of points, and the geological attributes of the points are calculated, and then association between the points is analyzed, and thus spatial geometrical bodies where the discrete points are located are combined, so that an overall geological distribution condition of the space is calculated, and the goals of high precision, little space limitation, high efficiency and strong applicability can be achieved.

Description

A kind of three-dimensional geological grid body Model construction method

Technical field

The present invention relates to and build three-dimensional geological body Model field, be specifically related to a kind of three-dimensional geological grid body Model construction method.

Background technology

During existing three-dimensional geologic builds, data acquisition mainly adopts discrete loose collection and single-point line string to gather two kinds of patterns.The collection of discrete loose point is exactly directly carry out discrete single-point sampling, does not gather the connection relation information between sampled point and sampled point, as gathered discrete elevation point, field acquisition discrete elevation point, Laser Radar Scanning collection etc. based on topomap.The collection of single-point line string is exactly gather and organising data in the mode of continuous coordinate string, as the collection of geography line in isocontour collection, digital photogrammetry, etc.The common feature of above-mentioned sampling is the spatial positional information only gathering sampled point, does not gather the relief form attribute information of the sampled point zone of action.In fact, each sampled point of topographical surface, it is not only space orientation data, the more important thing is its form properties information, as for Feature line, whether be smooth continuous curve or broken line, whether can carry out smooth treatment, a series of attributes such as smooth degree is how many in direction, longitudinal profile if just having in tangential direction.In theory, the relief form of each sampled point subrange can be reconstructed by its periphery sampled point of associating, but in actual applications, because the complexity of landform is far beyond the analog capability of existing mathematics and computer science, therefore when sampled point local form attribute information lacks, during the form properties loss of learning of especially feature terrain, be difficult to high-fidelity in reality and reconstruct its terrain surface.In fact this is also cause existing three-dimensional geological body Model mesorelief to express the one of the main reasons that there is critical regions property distortion phenomenon.

The main method that current three-dimensional geologic builds is for platform with Arcgis, Surpac, QvCoalMine etc., carry out matching with manually repairing and knit the method for closing and build three-dimensional geological body Model, three-dimensional geological body Model is built with these traditional methods, need a large amount of artificial, its accuracy is also limited by applied platform technology, all have a greatly reduced quality in precision, efficiency etc., so, need badly a kind of based on geologic data on a large scale, the three-dimensional geological object model technology being background with geological form on a large scale.

2009, Paradigm Geophysical Corp. of the U.S. proposed to be spatial spread by column, and then the method for making of the inverting structure three-dimensional geological model of the distribution of plastid spatially, but it is in optional position, space calculation limitation, and model accuracy is limited.

Chinese Patent Application No.: 201310238811, publication date: on September 4th, 2013 patent document discloses a kind of three-dimensional horizontal geologic profile figure method for making and uses thereof.The method 1. foundation energy and mines geological and requirements of engineering geology associated specifications, obtain field geology source book, comprise and be not limited to probing geologic information, probing point three dimensional local information, a geology geophysical prospecting informations, and in associated chamber, 10 steps such as outdoor test data have walked.This invention is used for: Mapping of Engineering Geology energy and mines geological drawing of the present invention; The expression of geological mapping achievement in geological mapping in engineering geology and energy and mines geological; Comprise and to be not limited in engineering investigation, design and analysis, construction, monitoring, operation link the expression of relevant investigation result; Comprise and be not limited to energy minerals planning, exploration, design, monitoring, exploitation, operation, repair the expression of relevant investigation result in link.This patent documentation by existing original documents and materials determination histogram carry out inverting spatially plastid distribution build three-dimensional geological model instruct geological mapping work, because source book is by the limitation of age, at that time exploring equipment level, technician's professional standards, geology change etc., the reliability of the precision of constructed three-dimensional geological model allows people suspect unavoidably.

Summary of the invention

1. the technical matters that will solve

For at present based on geologic data on a large scale, the three-dimensional geological object model being background with geological form on a large scale, the problem of precision, efficiency limitation, of the present inventionly provides a kind of three-dimensional geological grid body Model construction method.This construction method becomes a series of point spatial spreading, by the geological property of calculation level, then relevance is between points analyzed, and then merge discrete point place space geometric from, thus calculate the geology distribution situation of space entirety, can realize precision high, limited space is little, efficiency is high, the object that application is strong.

2. technical scheme

For achieving the above object, technical scheme provided by the invention is:

A kind of three-dimensional geological grid body Model construction method, idiographic flow is as follows:

Step 1, confirmed standard engineering planar range: known standard engineering planar range is rectangular area [(x ₁, x ₂), (y ₁, y ₂)], the unit of grid body is: u [u _x, u _y, u _z], obtain standard rectangular region G [(x _min, x _max), (y _min, y _max)];

Step 2, standard rectangular region G rasterizing: with an O (0,0) for initial point, spacing is [u _x, u _y], rasterizing is carried out to standard rectangular region G, standard rectangular region G is divided into n _x× n _yindividual continuous grids unit, n _xfor for get grid cell centers point as grid computing point, then corresponding four grid nodes of each grid cell and a grid computing point; Then standard rectangular region G is expressed as (n _x+ 1) × (n _y+ 1) individual grid node and n _x× n _ythe point set P of individual grid computing point composition; P will be put _i[(spatial value), some mark] enters point set P by ranks ordered storage; p _i∈ P, some mark distinguishes p _igrid computing point or grid node;

Step 3, the earth's surface triangulation network calculating standard rectangular region G and the ground end triangulation network:

Earth's surface triangulation network TRI is built according to boring top coordinate _t; Ground end triangle T RI is built according to boring end coordinate _b, and the triangulation network subdivision of each large layer;

Step 4, the calculating line calculating standard rectangular region G unit grid and mesh lines:

Calculation level p _iat triangulation network TRI _t, TRI _bon subpoint: owing to dividing equally a little in point set P in cloth standard rectangular region G, so at triangulation network TRI _t, TRI _bon all can find the subpoint of this some correspondence, two subpoints connecting this point obtain line l _i; As fruit dot p _ifor grid computing point, then l _ifor calculating line; As fruit dot p _ifor grid node, then l _ifor boost line; According to the determined plane equation of triangle, calculate a p _iat triangulation network TRI _t, TRI _bon the Z value of subpoint be respectively Z _t, Z _b;

Then Current standards works scope is expressed as n _x× n _ybar calculates line and (n _x+ 1) × (n _y+ 1) bar boost line, calculates corresponding four boost lines of line by grid node and the known every bar of grid computing point correspondence;

Step 5, equidistant partition boost line and calculating line:

With for reference point, with u _zfor unit head, on average divide boost line, title division points is auxiliary point;

Be reference point with Z=0, with u _zfor unit head, average division calculates line, and title division points is calculation level;

Eight auxiliary points around corresponding with the known calculation level of calculating line mechanism by equidistant partition boost line, corresponding auxiliary point is built into the cell cube at this calculation level place;

Step 6, calculating discrete point geological property, determine the cell cube at each discrete point place, and the geological property of confirmation unit body;

Using the stratigraphic distribution situation of the borehole space positional information collected, boring, geologic background parameter, geology geophysical prospecting informations, hole comprehensive layering data, space of discrete points position as Given information, by the uncertain method of fuzzy, probability, stochastic process, metric calculation is carried out to discrete point, obtains the probability distribution situation of often kind of geological property at this discrete point; The geological property of getting the maximal value in Making by Probability Sets L is the geological property of this discrete point; The probability got in set L is impartial, is considered as the boundary that this discrete point is in geologic body; Getting any geological property is the geological property of this discrete point;

Step 7, set up corresponding relation between calculation level and auxiliary point;

A: calculating line and boost line are numbered, coding rule: [row number, line number], is designated as [I, J]

B: set up the corresponding relation calculating line and boost line, is numbered being numbered of boost line corresponding to the calculating line of [I, J]: [I, J], [I+1, J], [I+1, J+1], [I, J+1];

C: successively to auxiliary point in the boost line of every bar, be numbered with natural number order from top to bottom; Successively to calculation level on the calculating line of every bar, be numbered with natural number order from top to bottom;

D: from top to bottom, finds out the auxiliary point that each calculation level is corresponding successively; The principle of its correspondence is: auxiliary point Z Distribution value is positioned at [calculation level Z value-u _z, calculation level Z value+u _z] in interval;

Step 8, build three-dimensional geological grid cell body: Current standards works scope is expressed as by multiple comprise geological property, joint account information, positional information cell cube assemble the three-dimensional geological model of T, wherein cell cube T _{i, j, k}∈ T;

Eight auxiliary points around a calculation level correspondence, corresponding auxiliary point is built into the cell cube at this calculation level place, and build three-dimensional geological grid cell body, the geological property of calculation level is the geological property of this calculation level place cell cube;

Step 9, Column Row, according to computation sequence and cell cube positional information, classify cell cube; By construction unit body, each calculates the cell cube that on line, calculation level is corresponding, forms cell cube row, consistent with the line number and columns calculating line;

Step 10, according to the computing information of cell cube, positional information, cell cube to be merged; Merge the cell cube that geological property is identical, and obtain the distribution entity of often kind of geology;

Preferably, in described step 5, auxiliary point Z Distribution value is Z=[n ₂, { n}, n ₁] × u _z, wherein n is interval [n ₂, n ₁] between natural number set; Calculation level Z Distribution value is: Z=[m ₂, { m}, m ₁] × u _z, wherein m is interval [m ₂, m ₁] between natural number set.

Preferably, in described step 6, the geological property computation process of discrete point is:

A: establish total m kind geological property in borehole data, then its probability distribution is set L:{l1, l2, l3 ..., lm}, wherein $\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{l}_i=1;$

B: get the maximal value in set L, write down and be designated as i, then Li>Lj, wherein j=1,2 ..., i-1, i+1 ... m; Then i-th kind of geological property is the geological property of this discrete point;

C: get the maximal value in set L, write down and be designated as i1, i2 ..., ik k altogether, and Li1=Li2=...=Lik>Lj, wherein j ≠ i1, i2 ..., ik; Then this discrete point is in the boundary of geologic body, remembers that this geological property is the i-th 1, i2 ..., any one geological property in ik.

Preferably, in described step 8, three-dimensional geological grid body Model construction step is:

A: eight auxiliary points around corresponding with the known calculation level of calculating line mechanism by equidistant partition boost line, corresponding auxiliary point is built into the cell cube at this calculation level place;

B: the face of this cell cube is numbered, add up the type of all of this cell cube, face is divided into vertical plane and surface level, add up the some numbering that each bread contains, the subscript of index face in stress and strain model, by gained statistical information i.e. this cell cube joint account information;

C: remember that the geological property of this calculation level is the geological property of this calculation level place cell cube;

Then current works scope is expressed as by multiple three-dimensional geological grid body of cell cube set T comprising geological property, joint account information, positional information, wherein cell cube T _{i, j, k}∈ T.

Preferably, in described step 9, cell cube sorting technique is:

With cell cube T _{i, j, k}for example, its positional information [i, j, k], i be column Ser.No. i ∈ 1,2 ... .., n _x, j by the Ser.No. j ∈ that is expert at 1,2 ... .., n _y, K is vertical direction Ser.No.;

Class1: not any with other cell cube is crossing;

Type 2: crossing with a cell cube;

Type 3: crossing with two cell cubes;

Type 4: simultaneously crossing with three cell cubes.

Preferably, the method that in described step 10, cell cube merges is:

A: the crossing cell cube of not any with other cell cube merges: cell body T _0,0,0for starting point plastid, and the geological property of this geologic body and cell cube T _0,0,0unanimously;

B: the cell cube intersected from a cell cube merges: the geological property as this cell cube and crossing cell cube is different, nonjoinder; Otherwise, find out the location plastid of crossing cell cube, from geologic body, delete intersection surface and its location prime number certificate of this cell cube, and all the other 5 faces are added in the location plastid of crossing cell cube; This cell cube and crossing cell cube are incorporated in same geologic body;

C: the cell cube intersected from two cell cubes merges: as this cell cube and two-phase hand over the geological property of cell cube all different, equal nonjoinder; As only with wherein a crossing cell cube geological property is identical, then carry out merging according to B mode; As identical in three cell cube geological property, cell cube is added in the plastid of adjacent cells body location, and the geologic body at crossing cell cube place need be merged;

D: the cell cube intersected from three cell cubes merges: as this cell cube and two-phase hand over the geological property of cell cube all different, equal nonjoinder; As only with wherein a crossing cell cube geological property is identical, then carry out merging according to C mode; As identical in four cell cube geological property, merge the cell cube of four cell cube Cheng Xin.

3. beneficial effect

Compared to immediate prior art, usefulness of the present invention is:

(1) a kind of three-dimensional geological grid body Model construction method of the present invention, the three-dimensional geological model of three-dimensional geological grid body Model architecture, straightforward discloses abundant stratum space distribution information, and the original geological information of direct correlation, make technician form complete, geology three-dimensional spatial distribution form concept accurately fast, for applied geology information provide more effectively, comprehensively, technical support clearly;

(2) a kind of three-dimensional geological grid body Model construction method of the present invention, concentrated expression whole engineering geology three-dimensional formation space distribution rule and variation tendency;

(3) a kind of three-dimensional geological grid body Model construction method of the present invention, reflect the geological condition of engineering geology and energy and mines geological accurately, the design and analysis can be correlated with for Geotechnical Engineering, design that construction is relevant with energy minerals, exploitation, reparation provide intuitively, accurately, comprehensive technical support, greatly improve engineering efficiency and reduce engineering cost.

Accompanying drawing explanation

Fig. 1 is standard area G grid node and grid computing point distribution plan;

Fig. 2 is works scope internal drilling plane distribution schematic diagram;

Fig. 3 is calculation level and auxiliary point 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

Below in conjunction with embodiment and accompanying drawing, explanation is further explained to technical scheme of the present invention.

Embodiment 1

First related notion is illustrated below:

Geologic body: certain geology, in the geometric shape of space distribution, refers in particular to similar stair-stepping geometric shape here.

Geologic body specifying information has the numbering and cell cube location plastid numbering that cell cube is numbered, geologic body contains cell cube.

A kind of three-dimensional geological grid body Model construction method of the present embodiment, idiographic flow is as follows:

Step 1, confirmed standard engineering planar range: known standard engineering planar range is rectangular area [(x ₁, x ₂), (y ₁, y ₂)], the unit of grid body is: u [u _x, u _y, u _z], obtain standard rectangular region G [(x _min, x _max), (y _min, y _max)]; Then G computing method are:

Step 2, standard rectangular region G rasterizing, ensure the geology distribution consistance of any two adjacent areas:

With an O (0,0) for initial point, spacing is [u _x, u _y], rasterizing is carried out to standard rectangular region G, standard rectangular region G is divided into n _x× n _yindividual continuous grids unit, n _xfor n _yfor get grid cell centers point as grid computing point, as shown in Figure 1, corresponding four grid nodes of each grid cell and a grid computing point; Then standard rectangular region G is expressed as (n _x+ 1) × (n _y+ 1) individual grid node and n _x× n _ythe point set P of individual grid computing point composition; P will be put _i[(spatial value), some mark] enters point set P by ranks ordered storage; p _i∈ P, some mark distinguishes p _igrid computing point or grid node;

Step 3, the earth's surface triangulation network calculating standard rectangular region G and the ground end triangulation network:

Earth's surface triangulation network TRI is built according to boring top coordinate _t:

As shown in Figure 2, the node in figure within the scope of black wire frame is the drilling plane distribution schematic diagram of certain engineering, and by the top and bottom coordinate of holing, interpolation goes out standard rectangular area G tetra-point [x respectively _min, y _min], [x _min, y _max], [x _max, y _min], [x _max, y _max] 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)}], [zt _{xmax, ymax}, zb _{xmax, ymax}]; Note boring apex coordinate and point [x _min, y _min, zt _{(xmin, ymin)}], [x _min, y _max, zt _{(xmin, ymax)}], [x _max, y _min, zt _{(xmax, ymin)}], [x _max, y _max, zt _{xmax, ymax}] be S set T; Adopt Delaunay triangulation to calculate and carry out the structure triangulation network:

A: build triangle on a large scale, comprise all loose point in all S set T, a large Atria summit is not included in S set T;

B: choose arbitrarily a bit from S set T, builds first Delaunay triangulation, comprises three triangles altogether;

C: the point not yet calculated in set of computations ST successively, find out two adjacent triangle in the public part being located at adjacent two triangle circumscribed circle, and delete the common edge of adjacent triangle, and connect this point and residue 2 points, form new Delaunay triangulation;

D: delete the triangle containing triangular apex on a large scale, last Delaunay triangulation can be obtained;

In like manner according to above method, build ground end triangle T RI _b, and the triangulation network subdivision of each large layer;

Step 4, the calculating line calculating standard rectangular region G unit grid and mesh lines:

Calculation level p _iat triangulation network TRI _t, TRI _bon subpoint: owing to dividing equally a little in point set P in cloth standard rectangular region G, so at triangulation network TRI _t, TRI _bon all can find the subpoint of this some correspondence, two subpoints connecting this point obtain line l _i; As fruit dot p _ifor grid computing point, then l _ifor calculating line; As fruit dot p _ifor grid node, then l _ifor boost line; According to the determined plane equation of triangle, calculate a p _iat triangulation network TRI _t, TRI _bon the Z value of subpoint be respectively Z _t, Z _b;

Then Current standards works scope is expressed as n _x× n _ybar calculates line and (n _x+ 1) × (n _y+ 1) bar boost line, calculates corresponding four boost lines of line by grid node and the known every bar of grid computing point correspondence;

Step 5, equidistant partition boost line and calculating line:

With Z= for reference point, with u _zfor unit head, on average divide boost line, title division points is auxiliary point, then auxiliary point Z Distribution value is Z=[n ₂, { n}, n ₁] × u _z, wherein n is interval [n ₂, n ₁] between natural number set;

Be reference point with Z=0, with u _zfor unit head, average division calculates line, and title division points is calculation level, then calculation level Z Distribution value is: Z=[m ₂, { m}, m ₁] × u _z, wherein m is interval [m ₂, m ₁] between natural number set;

As shown in Figure 3, eight auxiliary points around corresponding with the known calculation level of calculating line mechanism by equidistant partition boost line, corresponding auxiliary point is built into the cell cube at this calculation level place;

Step 6, calculating discrete point geological property, determine the cell cube at each discrete point place, and the geological property of confirmation unit body;

Using the stratigraphic distribution situation of the borehole space positional information collected, boring, geologic background parameter, geology geophysical prospecting informations, hole comprehensive layering data, space of discrete points position as Given information, by the uncertain method of fuzzy, probability, stochastic process, metric calculation is carried out to discrete point, obtains the probability distribution situation of often kind of geological property at this discrete point;

The geological property computation process of discrete point is:

A: establish total m kind geological property in borehole data, then its probability distribution is set L:{l1, l2, l3 ..., lm}, wherein $\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{l}_i=1;$

B: get the maximal value in set L, write down and be designated as i, then Li>Lj, wherein j=1,2 ..., i-1, i+1 ... m; Then i-th kind of geological property is the geological property of this discrete point;

C: get the maximal value in set L, write down and be designated as i1, i2 ..., ik k altogether, and Li1=Li2=...=Lik>Lj, wherein j ≠ i1, i2 ..., ik; Then this discrete point is in the boundary of geologic body, remembers that this geological property is the i-th 1, i2 ..., any one geological property in ik.

Step 7, set up corresponding relation between calculation level and auxiliary point;

A: calculating line and boost line are numbered, coding rule: [row number, line number], is designated as [I, J]

B: set up the corresponding relation calculating line and boost line, is numbered being numbered of boost line corresponding to the calculating line of [I, J]: [I, J], [I+1, J], [I+1, J+1], [I, J+1];

C: successively to auxiliary point in the boost line of every bar, be numbered with natural number order from top to bottom; Successively to calculation level on the calculating line of every bar, be numbered with natural number order from top to bottom;

D: from top to bottom, finds out the auxiliary point that each calculation level is corresponding successively; The principle of its correspondence is: auxiliary point Z Distribution value is positioned at [calculation level Z value-u _z, calculation level Z value+u _z] in interval;

Step 8, structure three-dimensional geological grid cell body:

Eight auxiliary points around A: one calculation level correspondence, corresponding auxiliary point is built into the cell cube at this calculation level place;

B: the face of this cell cube is numbered, add up the type of all of this cell cube, face is divided into vertical plane and surface level, add up the some numbering that each bread contains, the subscript of index face in stress and strain model, by gained statistical information i.e. this cell cube joint account information;

C: remember that the geological property of this calculation level is the geological property of this calculation level place cell cube;

As shown in Figure 4, then current works scope is expressed as by multiple the three-dimensional geological grid body T of cell cube set comprising geological property, joint account information, positional information, wherein cell cube T _{i, j, k}∈ T;

Step 9, Column Row, according to computation sequence and cell cube positional information, classify cell cube; By construction unit body, each calculates the cell cube that on line, calculation level is corresponding and forms cell cube row, consistent with the line number and columns calculating line, as shown in Figure 5;

With cell cube T _{i, j, k}for example, its positional information [i, j, k], i be column Ser.No. i ∈ 1,2 ... .., n _x, j by the Ser.No. j ∈ that is expert at 1,2 ... .., n _y, K is vertical direction Ser.No.;

Class1: not any with other cell cube is crossing;

Type 2: crossing with a cell cube;

Type 3: crossing with two cell cubes;

Type 4: simultaneously crossing with three cell cubes;

Step 10, according to the computing information of cell cube, positional information, cell cube to be merged; Merge the cell cube that geological property is identical, and obtain the distribution entity of often kind of geology:

Cell cube crossing for not any with other cell cube is merged: cell body T _0,0,0for starting point plastid, and the geological property of this geologic body and cell cube T _0,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 is that the cell cube type merged is different:

Merge from the cell cube that a cell cube intersects: the geological property as this cell cube and crossing cell cube is different, nonjoinder; Otherwise, find out the location plastid of crossing cell cube, from geologic body, delete intersection surface and its location prime number certificate of this cell cube, and all the other 5 faces are added in the location plastid of crossing cell cube; This cell cube is incorporated in same geologic body with crossing cell cube.

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 is that the cell cube type merged is different:

Merge from the cell cube that two cell cubes intersect: as this cell cube and two-phase hand over the geological property of cell cube all different, equal nonjoinder; As only with wherein a crossing cell cube geological property is identical, then carry out merging according to embodiment 2 mode; As identical in three cell cube geological property, cell cube is added in the plastid of adjacent cells body location, and the geologic body at crossing cell cube place need be merged;

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 is that the cell cube type merged is different:

Merge from the cell cube that three cell cubes intersect: as this cell cube and two-phase hand over the geological property of cell cube all different, equal nonjoinder; As only with wherein a crossing cell cube geological property is identical, then carry out merging according to embodiment 3 mode; As identical in four cell cube geological property, merge the cell cube of four cell cube Cheng Xin.

Below be schematically described the invention and embodiment thereof, actual technical scheme is not limited thereto.So, if those of ordinary skill in the art enlightens by it, when not departing from this creation aim, designing the embodiment similar to this technical scheme without creationary, the protection domain of this patent all should be belonged to.

Claims (6)

1. a three-dimensional geological grid body Model construction method, method flow is as follows:

Step 1, confirmed standard engineering planar range: known standard engineering planar range is rectangular area [(x ₁, x ₂), (y ₁, y ₂)], the unit of grid body is: u [u _x, u _y, u _z], obtain standard rectangular region G [(x _min, x _max), (y _min, y _max)];

Step 2, standard rectangular region G rasterizing: with an O (0,0) for initial point, spacing is [u _x, u _y], rasterizing is carried out to standard rectangular region G, standard rectangular region G is divided into n _x× n _yindividual continuous grids unit, n _xfor n _yfor get grid cell centers point as grid computing point, then corresponding four grid nodes of each grid cell and a grid computing point; Then standard rectangular region G is expressed as (n _x+ 1) × (n _y+ 1) individual grid node and n _x× n _ythe point set P of individual grid computing point composition; P will be put _i[(spatial value), some mark] enters point set P by ranks ordered storage; p _i∈ P, some mark distinguishes p _igrid computing point or grid node;

Step 3, the earth's surface triangulation network calculating standard rectangular region G and the ground end triangulation network:

Earth's surface triangulation network TRI is built according to boring top coordinate _t; Ground end triangle T RI is built according to boring end coordinate _b, and the triangulation network subdivision of each large layer;

Step 4, the calculating line calculating standard rectangular region G unit grid and mesh lines:

Calculation level p _iat triangulation network TRI _t, TRI _bon subpoint: owing to dividing equally a little in point set P in cloth standard rectangular region G, so at triangulation network TRI _t, TRI _bon all can find the subpoint of this some correspondence, two subpoints connecting this point obtain line l _i; As fruit dot p _ifor grid computing point, then l _ifor calculating line; As fruit dot p _ifor grid node, then l _ifor boost line; According to the determined plane equation of triangle, calculate a p _iat triangulation network TRI _t, TRI _bon the Z value of subpoint be respectively Z _t, Z _b;

Then Current standards works scope is expressed as n _x× n _ybar calculates line and (n _x+ 1) × (n _y+ 1) bar boost line, calculates corresponding four boost lines of line by grid node and the known every bar of grid computing point correspondence;

Step 5, equidistant partition boost line and calculating line:

With for reference point, with u _zfor unit head, on average divide boost line, title division points is auxiliary point;

Be reference point with Z=0, with u _zfor unit head, average division calculates line, and title division points is calculation level;

Eight auxiliary points around corresponding with the known calculation level of calculating line mechanism by equidistant partition boost line, corresponding auxiliary point is built into the cell cube at this calculation level place;

Step 6, calculating discrete point geological property, determine the cell cube at each discrete point place, and the geological property of confirmation unit body: using the stratigraphic distribution situation of the borehole space positional information collected, boring, geologic background parameter, geology geophysical prospecting informations, hole comprehensive layering data, space of discrete points position as Given information, by the uncertain method of fuzzy, probability, stochastic process, metric calculation is carried out to discrete point, obtains the probability distribution situation of often kind of geological property at this discrete point; The geological property of getting the maximal value in Making by Probability Sets L is the geological property of this discrete point; The probability got in set L is impartial, is considered as the boundary that this discrete point is in geologic body; Getting any geological property is the geological property of this discrete point;

Step 7, set up corresponding relation between calculation level and auxiliary point;

A: calculating line and boost line are numbered, coding rule: [row number, line number], is designated as [I, J]

B: set up the corresponding relation calculating line and boost line, is numbered being numbered of boost line corresponding to the calculating line of [I, J]: [I, J], [I+1, J], [I+1, J+1], [I, J+1];

C: successively to auxiliary point in the boost line of every bar, be numbered with natural number order from top to bottom; Successively to calculation level on the calculating line of every bar, be numbered with natural number order from top to bottom;

D: from top to bottom, finds out the auxiliary point that each calculation level is corresponding successively; The principle of its correspondence is: auxiliary point Z Distribution value is positioned at [calculation level Z value-u _z, calculation level Z value+u _z] in interval;

Step 8, build three-dimensional geological grid cell body: Current standards works scope is expressed as by multiple comprise geological property, joint account information, positional information cell cube assemble the three-dimensional geological model of T, wherein cell cube T _{i, j, k}∈ T;

Eight auxiliary points around a calculation level correspondence, corresponding auxiliary point is built into the cell cube at this calculation level place, and build three-dimensional geological grid cell body, the geological property of calculation level is the geological property of this calculation level place cell cube;

Step 9, Column Row, according to computation sequence and cell cube positional information, classify cell cube; By construction unit body, each calculates the cell cube that on line, calculation level is corresponding and forms cell cube row, consistent with the line number and columns calculating line;

Step 10, according to the computing information of cell cube, positional information, cell cube to be merged; Merge the cell cube that geological property is identical, and obtain the distribution entity of often kind of geology.

2. a kind of three-dimensional geological grid body Model construction method according to claim 1, it is characterized in that: in described step 5, auxiliary point Z Distribution value is Z=[n ₂, { n}, n ₁] × u _z, wherein $n_1=\frac{z_t-0.5\×u_z}{u_z},n_2=\frac{z_b-0.5\×u_z}{u_z},$ N is interval [n ₂, n ₁] between natural number set; Calculation level Z Distribution value is: Z=[m ₂, { m}, m ₁] × u _z, wherein m is interval [m ₂, m ₁] between natural number set.

3. a kind of three-dimensional geological grid body Model construction method according to claim 1 and 2, it is characterized in that: in described step 6, the geological property computation process of discrete point is:

A: establish total m kind geological property in borehole data, then its probability distribution is set L:{l1, l2, l3 ..., lm}, wherein $\underset{i=1}{\overset{m}{\mathrm{\Σ}}}{l}_i=1;$

B: get the maximal value in set L, write down and be designated as i, then Li>Lj, wherein j=1,2 ..., i-1, i+1 ... m; Then i-th kind of geological property is the geological property of this discrete point;

C: get the maximal value in set L, write down and be designated as i1, i2 ..., ik k altogether, and Li1=Li2=...=Lik>Lj, wherein j ≠ i1, i2 ..., ik; Then this discrete point is in the boundary of geologic body, remembers that this geological property is the i-th 1, i2 ..., any one geological property in ik.

4. a kind of three-dimensional geological grid body Model construction method according to claim 1 and 2, is characterized in that: in described step 8, three-dimensional geological grid body Model construction step is:

Eight auxiliary points around A: one calculation level correspondence, corresponding auxiliary point is built into the cell cube at this calculation level place;

B: the face of this cell cube is numbered, add up the type of all of this cell cube, face is divided into vertical plane and surface level, add up the some numbering that each bread contains, the subscript of index face in stress and strain model, by gained statistical information i.e. this cell cube joint account information;

C: remember that the geological property of this calculation level is the geological property of this calculation level place cell cube;

Then current works scope is expressed as by multiple three-dimensional geological grid body of cell cube set T comprising geological property, joint account information, positional information, wherein cell cube T _{i, j, k}∈ T.

5. a kind of three-dimensional geological grid body Model construction method according to claim 4, is characterized in that: in described step 9, cell cube sorting technique is:

With cell cube T _{i, j, k}for example, its positional information [i, j, k], i be column Ser.No. i ∈ 1,2 ... .., n _x, j by the Ser.No. j ∈ that is expert at 1,2 ... .., n _y, k is vertical direction Ser.No.;

Class1: not any with other cell cube is crossing;

Type 2: crossing with a cell cube;

Type 3: crossing with two cell cubes;

Type 4: simultaneously crossing with three cell cubes.

6. a kind of three-dimensional geological grid body Model construction method according to claim 5, is characterized in that: the method that in described step 10, cell cube merges is:

A: the crossing cell cube of not any with other cell cube merges: cell body T _0,0,0for starting point plastid, and the geological property of this geologic body and cell cube T _0,0,0unanimously;

B: the cell cube intersected from a cell cube merges: the geological property as this cell cube and crossing cell cube is different, nonjoinder; Otherwise, find out the location plastid of crossing cell cube, from geologic body, delete intersection surface and its location prime number certificate of this cell cube, and all the other 5 faces are added in the location plastid of crossing cell cube; This cell cube and crossing cell cube are incorporated in same geologic body;

C: the cell cube intersected from two cell cubes merges: as this cell cube and two-phase hand over the geological property of cell cube all different, equal nonjoinder; As only with wherein a crossing cell cube geological property is identical, then carry out merging according to B mode; As identical in three cell cube geological property, cell cube is added in the plastid of adjacent cells body location, and the geologic body at crossing cell cube place need be merged;

D: the cell cube intersected from three cell cubes merges: as this cell cube and two-phase hand over the geological property of cell cube all different, equal nonjoinder; As only with wherein a crossing cell cube geological property is identical, then carry out merging according to C mode; As identical in four cell cube geological property, merge the cell cube of four cell cube Cheng Xin.