CN100465998C - Method for finding mine with cubic forecast model based on three-dimensional modeling - Google Patents

Abstract

The invention discloses a mine finding method based on a cube forecasting model with three-dimensional modeling, characterised in that having following steps: (1) collecting data; (2) establishing a database; (3) establishing a three-dimensional mine model; (4) establishing a cube forecasting model; (5) validating a forecasting result. The invention can find a mine in a large scale study region where a twodimensional mine finding method can't realized by using the method of three-dimensional modeling during a mine finding process, the obtained forecasting result can reflect an actual situation more factually. With the cube forecasting model with three-dimensional modeling, various three-dimensional space data can be used to establish a three-dimensional geological model, a three-dimensional mine model and a three-dimensional mine finding model of the study region, the forecasting result is qualified in position, quantity and probability. Accordingly, the ''cube forecasting model'' mine finding method realizes a development from the two-dimensional mine finding to the three-dimensional mine finding, has important theory value and application value.

Description

Method for finding mine with cubic forecast model based on three-dimensional modeling

Technical field

The present invention relates to a kind of method of looking for the ore deposit, specifically a kind of method for finding mine with cubic forecast model based on three-dimensional modeling.

Background technology

At present, China many super large metallic ores Ji Qu is because neoteric and modern large scale mining, and the resource reserve of verifying seriously consumes, and the mining industry situation is very severe, and the economy and the society in these areas caused tremendous influence.Seek the new matter of utmost importance that resource has become these regional mining industry survival and development to solve of taking over.

The evaluation of small scale SOLID MINERAL RESOURCES is summed up from the simple regularity of ore formation, instruct and be summarized as ore deposit geologic background and the regularity of ore formation and be support with 3S technology (GIS, GPS, RS), geomathematics, computer technology etc. towards using modern metallogenic theory, the important meals mineral resource assessment direction of utilizing the Multielement Information System synthetic study to carry out different levels develops.Mineral resource prediction under the large scale condition of mining area then is the key that can really solve the mineral products urban sustainable development.

Traditional two-dimentional method of prospecting be on the basis of metallogenic prognosis theory in conjunction with geophysics, geochemistry and the mineral exploration engineering data on the face of land, reconnoitre the space geometry form and the scale of mineral resources.These mine locating methods have developed comparative maturity at present, from geophysical prospecting method, natural method, artificial field method, astable method are arranged, and gravitational method, magnetic method, electrical method, seismic method, radioactivity method are arranged, and ground and underground survey etc. are arranged; From sampling, surface occurence, exploratory trench, tunnel, drill hole sampling method etc. are arranged.The key issue that present two dimension looks for the ore deposit to solve is how to reject noise from so a large amount of two-dimensional signals, distinguish non-ore anomaly and cause ore anomaly, therefore in information processing, use mathematical method, formed a series of information processing technologies such as mineral deposit statistical forecast, collective geological prospecting quantity of information method, evidence weight method.These technical methods are looked in the process of ore deposit in middle small scale zone and have been played important role.

Yet the mineral deposit is not the plane, and it is to be present in the geologic body that three dimensions has Special Significance, and the ore deposit information of looking for that the face of land obtains is the stack of underground three-dimensional spatial information.The essence of the disposal route of two-dimensional signal above-mentioned is that the three dimensions with the two dimensional surface that is added to becomes ore deposit information to extract, and extracting to such an extent that the distortion or the loss of information can take place in the process unavoidably, causes incorrect conclusion.Therefore directly it is extremely necessary carrying out the mineral deposit prediction from three parallactic angle degree, also is an important directions of geological prospecting from now on.

Summary of the invention

The objective of the invention is to remedy and break through traditional two-dimentional method of prospecting, propose a kind of method for finding mine with cubic forecast model based on three-dimensional modeling.

The object of the present invention is achieved like this: the three-dimensional method of prospecting of the present invention is directly to use the three-dimensional space data of gathering under the face of land, use appropriate information to merge and extracting method then, extract the ore anomaly that causes in these data, in conjunction with the deposit genesis analysis, set up the three-dimensional geological of study area and look for the ore deposit model, the target area, ore deposit is looked in final proposition.Set up cubic forecast model on the basis of three-dimensional modeling, it can use various three-dimensional space datas, set up study area three-dimensional geological model, three-dimensional become ore deposit model and three-dimensional to look for the ore deposit model, it predicts the outcome and can accomplish location, quantitative and decide probability.

Carry out the detailed description of concrete technical scheme below, so that the clearer understanding of those skilled in the art content of the present invention:

1, Data acquisition;

Should collect mining geology, geochemistry, mine engineering data.Wherein, mining geology and mine engineering data obtain by collecting stage casing, mine geological map, mine centre section project arrangenent diagram and exploratory grid cross section (CAD drawing), and go deep into mine and understand mine engineering operation situation; Obtain in the common sample resources data file of element anomalies data from borehole data storehouse, mine.

2, Database Construction:

Set up borehole data storehouse, mining area.Wherein, can set up mining area element anomalies three-dimensional model by the common sample resources data file in the borehole data storehouse.The borehole data storehouse is the integrated foundation of Micromine KAN TAN 3D by the research and development of Australian MICROMINE company, and can carry out unified management.

The concrete establishment step in borehole data storehouse is as follows:

(1) collect mining area borehole data (drawing or Access database), the database file that need use mainly contains: aperture resources data file, deviational survey resources data file, common sample resources data file.Software platform is: Microsoft Access.

(2) the necessary attribute of using in the resources data file of aperture has: boring numbering, three-dimensional geographic coordinate (X, Y, Z), design inclination angle, design orientation, design hole depth, whole hole depth; Other attributes comprise: unit in charge of construction, ore deposit section, stage casing, tunnel title, ore body number, classification, perforate date, whole hole date, remarks, whether the emphasis drilling depth, edit and record people, ID.

(3) the necessary attribute of using in the deviational survey resources data file has: boring numbering, the deviational survey degree of depth, inclination angle, position angle, ID.

Can in subsequent step, set up the three-dimensional model of boring form according to above three data files.

(4) the necessary attribute of using in the common sample resources data file has: boring numbering, sample number into spectrum, sampling depth, sample length, long, the required element grade of core; Other attributes comprise: ore type, remarks, three-dimensional geographic coordinate, ID.

Can set up mining area element anomalies three-dimensional model according to the common sample resources data file in subsequent step, concrete steps are seen below.

3, the foundation of three-dimensional mine model;

The mine three-dimensional model that needs to set up comprises face of land DEM model, stratigraphic model, fracture model, granite model, all types of known ore body model, element anomalies three-dimensional model.The concrete grammar division is as follows:

(1) DEM model in the face of land is on the basis of collecting mining area surface elevation drawing (tiff form), use the MAPGIS Geographic Information System (GIS) software of China Geological Univ. Wuhan's independent research and carry out vector quantization, face of land elevation is changed into vector file, import to the foundation that realizes among the three-dimensional modeling visual software Micromine face of land three-dimensional entity model.

(2) method for building up of stratigraphic model, fracture model, granite and all types of known ore body models is roughly the same.

The concrete steps that three-dimensional mine model is set up are as follows:

1), collects section

At first collect stage casing, mining area geological map and exploratory grid cross section (CAD drawing, the dwg form), under AutoCAD 2006 software platforms, planimetric map collection work district section data is reconnoitred in the contrast workspace, and checks that collected section is whether corresponding with the section arranged in the planimetric map.

2), stack, adjustment section

1., in AutoCAD, open the platform master plan of certain sea level elevation.

2., in planimetric map, insert sectional view.

3., adjust profile position.

3), divide dissection surface

All sections are all decomposed:

In vertical view, choose profile line → modification → decomposition.How to decompose several times, till can not decomposing again.

Then file is preserved:

File → save as saves as AutoCAD R12/LT2 DXF form → preservation with file.

4), import Micromine

Under Micromine KAN TAN 3D software platform, section file that will pre-service is good in AutoCAD imports, and realizes the geologic body d solid modeling in this software.

1., a newly-built engineering;

2., a section is imported in the engineering;

3., according to step method 2. all sections are all imported among the Micromine.

5), editor's section

1., access section;

2., delimit the modeling scope;

3., editor's section.

6), three-dimensional modeling

1., newly-built type;

2., newly-built wire frame;

3., modeling;

4., set up the three-dimensional stereo model of modeling scope.

According to the method that the atural object three-dimensional model is set up, set up the three-dimensional stereo model of a modeling scope.

5., the three-dimensional Boolean calculation of atural object entity

Boolean (Boolean calculation) is by plural object is carried out the computing of union, difference set, common factor, thereby obtains a kind of operational method of new object form.

4 kinds of boolean operation methods are provided among the Micromine KAN TAN 3D.Be respectively that A hands over B (Ain B), A to subtract B (A out B), B to hand over A (B in A), B to subtract A (B out A).The logical meaning of their representatives is respectively:

A hands over B (A in B): output A hands over the part in B;

A subtracts B (A out B): the part of output A outside B;

B hands over A (B in A): output B hands over the part in A;

B subtracts A (B out A): the part of output B outside A.

More than four kinds of fundamental operation methods more can be incorporated into more complex calculations.For example:

A, A represent the modeling scope, and B represents the face of land.

Operational method: A hands over B, B to hand over A.

The model that obtains is the stack model on the modeling scope and the face of land:

B, A represent the stack model on the modeling scope and the face of land, the lower surface on first stratum that B representative and surface distance are nearest.

Operational method is: A subtracts B, B and hands over A.

The model that obtains is: the uppermost stratum of crop out.

The ore body of c, A representative and grouan meet, B represents the grouan solid model.

Operational method is: A subtracts B, B and hands over A.

The model that obtains is that ore body is exposed at the outer department pattern of granite.

(3) foundation of boring model

The boring model is made up of 3 files of raw data base, be respectively drilling orifice three-dimensional coordinate file, boring skew back file and analysis data file, wherein, preceding two files are used to set up the three-dimensional model of boring form, and the 3rd common sample resources data file is used to set up the element anomalies three-dimensional model.Study area horse drinking water area is positioned at No. 10 ore cluster south of reed dyke, and the research of this ore cluster is historical long, and the boring data with existing is more comprehensive, sets up so the foundation in borehole data storehouse can make full use of existing borehole data.Master file has: reed dyke drilling orifice data, deviational survey data, common sample resources data file.

Concrete steps are:

1) the Access database is imported among the Micromine

2) calibration database

3) create new database

4) load boring

4, the foundation of cubic forecast model;

1) extraction of geologic body three-dimensional cube

According to the announcement of geologic information to ore body, the particularly distribution of exploratory line, determine the scope and the basic parameter of modeling in conjunction with form, trend, tendency and the spatial distribution characteristic of ore body, the scope model is divided into the cell block of row * row * layer for 20m * 20m * 20m, adopt time piecemeal in edge of model, its specification is 10m * 10m * 10m.

Behind the three-dimensional grid model of setting up, can will look for digital model determined Prediction Parameters in ore deposit to compose to each cell block as attribute.Use the stratum solid model that 3D grid is limited, mark off the grid cell that Different Strata comprises, as the stratum variable in the prediction of mineral deposit; Use the fracture solid model that 3D grid is limited, mark off the grid cell that different fractures are comprised, as the fracture variable in the prediction of mineral deposit.Use known ore body solid model that 3D grid is limited, mark off the grid cell that dissimilar known ore body comprises, as the priori conditions in the prediction of mineral deposit.

Three-dimensional cube extracts concrete steps:

1., set up empty scope nugget model:

2., load piece

3., give the cell block ascription

After setting up this study area three-dimensional grid model, can utilize common sample data in the borehole data storehouse to analyze three-dimensional unusual distribution of element of these cell blocks.Required constituent content to sample in the boring is analyzed, and the weighting of service range inverse ratio is carried out required constituent content interpolation to the unknown elements piece.

Analytical procedure is as follows:

1. the statistical study of grade data: common sample data statistics is analyzed the mathematical distribution rule that obtains the required element in this area, and find extra-high-speed value, minimum value, mean value and the variance of these elements.

2. the combined treatment of sample: the back of averaging of each sample grade weighting in the combined segment is formed new boring.Getting sample length is 0.3, cutoff grade 0.2, and the extra-high-speed value gets 10%.

3. limit boring model effective coverage with the empty block models file of study area scope.

4. the length combination of sample: for the ease of interpolation, decompose the sample that sample is a regular length again, get that the combination sample is long to be 10m.

5. sample interpolation: with the common sample data after the length combination is source data, service range inverse ratio weighted method, it is 2 that the power inverse is set, minor increment 5m, the ellipsoid that definition search is oval, 50m successively, 100m, 150m, 200m are that radius search known sample point comes the required element of unknown nugget is carried out interpolation.

The three-dimensional unusual figure of analytical element, and the anomaly threshold of required element is carried out value according to the unusual figure of three-dimensional.

It is as follows that the element anomalies three-dimensional model is set up concrete steps:

A, set up study area boring solid model according to the method for building up of aforementioned boring model;

B, the different elements that can study are as required set up different element anomalies three-dimensional models.

2) calculating of geological prospecting variable and extraction

According to the data that study area was collected, confirmable types of variables comprises stratum, fracture, rock mass and geochemical exploration, respectively this four class variable is carried out the division of cubic units and extracts to such an extent that step is as follows.

Grasp the situation on following latent stratum, the study area face of land, with all latent stratum respectively as a variable;

Getting rupture surface both sides 150m (this value is decided on different deposit types) is buffer zone, sets up fracture buffer zone variable.

With the above 150m of rock mass end face (this value is decided on different deposit types) is buffer zone, sets up the rock mass buffer variable.

Required research constituent content to sample in the element anomalies three-dimensional model is analyzed, and obtains the constituent content of each cubic units, and each is studied element respectively as a predictive variable.

3) the collective geological prospecting information of calculating each geologic unit and being comprised

When calculating the collective geological prospecting information that each geologic unit comprised, the software platform of using is the Access database software in Microsoft Office 2003 software packages, to import among the Access by the autonomous integrated database of Micromine KAN TAN 3D, Access provides the several data statistical query that comprises SQL ad hoc inquiry function, can satisfy worker's different statistical demand.Concrete steps are as follows:

1., scrutinizing on main ore-controlling factors and the favourable indicator for deposit basis, choose geology, the geochemical exploration Information sign of statistical study variable, and each sign of agreement to have value in the unit be 1, do not exist and get 0, each sign of statistics is in the distribution of each unit.In the computation process, the high size of nugget length and width is unified for 10m * 10m * 10m (alleged unit sizes promptly), calculate then.

2., look for the ore deposit quantity of information to calculate

The ore deposit quantity of information of looking for of certain indicator for deposit is calculated with conditional probability, that is:

$I_{A\left(B\right)}=\lg \frac{P(A/B)}{P\left(A\right)}---\left(1\right)$

I in the formula _{A (B)}The quantity of information that the B ore deposit is arranged for the A sign; P (A/B) is for the probability of A occurring under the existence condition of the known B of having ore deposit; P (A) is for the probability of sign A occurring in study area.

Because the difficulty on the probability estimate is come the estimated probability value with frequency values.At this moment

$I_{A\left(B\right)}=\lg \frac{\left(\frac{N_j}{N}\right)}{\left(\frac{S_j}{S}\right)}---\left(2\right)$

N in the formula _jContain the ore deposit unit number for what have sign A in the study area; N contains the ore deposit unit number in the study area; S _jFor having the unit number of sign A in the study area; S is the unit sum of study area.

General Two-Dimensional looks for ore deposit quantity of information method to determine to look for the favourable sign in ore deposit with formula 3.

${\mathrm{AI}}^{+}=K\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{I}_j---\left(3\right)$

Under the situation of variable negligible amounts, just do not calculate the accumulative total critical value of useful information, but directly use all on the occasion of the information number.

Utilize each indicator for deposit quantity of information to calculate the quantity of information summation of each elementary cell respectively, determine to look for ore deposit quantity of information critical value according to the subjective probability method then, the distant view unit of determining, ore deposit that contains comprises known ore body unit, contains therefore that the prospecting evaluation probability is in the distant view unit, ore deposit: contain the known ore body unit number/known ore body unit number that comprises in the distant view unit, ore deposit

4) the ore deposit probability calculation is distinguished in prediction

With a regional classification is some unit, has four kinds of situations:

(1) blank;

(2) x unit, expression have certain geology sign to exist;

(3) xy unit, expression have this kind geology sign to have the ore deposit to exist again;

(4) y unit, expression have the ore deposit to exist, but do not have this geology sign.

Corresponding frequency is respectively: f ₀, f _x, f _Xy, f _y

Appoint in the unit that has certain geology sign to exist and get 1, it may contain the probability P in 1 mineral deposit _xFor:

$P_x=\frac{f_{\mathrm{xy}}}{(f_x+f_{\mathrm{xy}})}---\left(4\right)$

Appoint in the unit that does not have certain geology sign to exist and get 1, it may contain the probability P in 1 mineral deposit ₀For:

${\mathrm{<figure-callout\; id="0"\; label="P"\; filenames="C200710098940E00561.png,C200710098940E00802.png"\; state="\{\{state\}\}">P</figure-callout>}}_0=\frac{f_y}{(f_y+f_0)}---\left(5\right)$

Setting looks for ore deposit informational capacity lower limit as the comprehensive geology sign, utilizes probability estimate method Estimation Study district ore body unit output number to calculate f _x+ f _XyAnd f _Xy, and obtain P according to (4) formula _xCalculate f ₀+ f _yAnd f _y, and obtain P according to (5) formula ₀

Usually the calculating here not have to consider beyond known mineral deposit the also mineral deposit of possible output, needs the control zone of selecting working level high for this reason, obtains a correction factor, is used for revising P _xAnd P ₀, get P _x' and P ₀', be used for the prediction district with revised probability, obtain predicting the estimation of district's mineral deposit number.Concrete steps are as follows:

The study area of selecting certain elevation interval is that the control zone is revised, unit, calculation control district sum and gross information content are positioned at the unit number on the lower limit, correction coefficient is=the control zone ore body unit number of the known ore body unit number in control zone/obtained by the output probability that calculates above.

Proofread and correct the probability P after obtaining proofreading and correct with correction coefficient _x' and P ₀', f _x* P _x'+f ₀* P ₀' can be used for estimating the ore body unit output number in distant view prediction district.

5, the checking that predicts the outcome.

After obtaining the one-tenth ore deposit probability of workspace, engineering verification is carried out in the work that will carry out exactly on the spot to the field.

Advantage of the present invention and benefit:

The method that adds three-dimensional modeling in looking for the ore deposit process can be looked for the ore deposit at the large scale study area that two dimension looks for ore deposit institute not reach.Three-dimensional modeling can play a significant role in the metallogenic prognosis process of mineral deposit yardstick, can progressively deepen the mineralization understanding to study area in the three-dimensional modeling method application process, obtains predicting the outcome and can reflect actual conditions more really.

Cubic forecast model based on three-dimensional modeling, can make various three-dimensional space datas, set up the three-dimensional geological model of study area, three-dimensional one-tenth ore deposit model and three-dimensional and look for the ore deposit model, it predicts the outcome and can accomplish location, quantitative and decide probability: locate and can determine to predict the three dimensional space coordinate of target area exactly, quantitatively be the amount of metal that can determine that the target area is possible, deciding probability is exactly the one-tenth ore deposit possibility size of determining in the study area scope.Therefore, " cubic forecast model " method of prospecting has realized the breakthrough of looking for the ore deposit to look for the ore deposit to three-dimensional from two dimension having great value in theory and using value.

Description of drawings

Fig. 1 is " input MDB file " dialog box synoptic diagram;

Fig. 2 is " database successfully imports prompting " dialog box synoptic diagram;

Fig. 3 is " boring verification " dialog box synoptic diagram;

Fig. 4 is " well head territory " dialog box synoptic diagram;

Fig. 5 is " measurement territory " dialog box synoptic diagram;

Fig. 6 is " interval territory " dialog box synoptic diagram;

Fig. 7 is " creating the borehole data storehouse " dialog box synoptic diagram;

Fig. 8 is " borehole data storehouse " dialog box synoptic diagram;

Fig. 9 is " adding interval file " dialog box synoptic diagram;

Figure 10 is " loading boring " dialog box synoptic diagram;

Figure 11 is " color sign indicating number " dialog box synoptic diagram;

Figure 12 absorbs water mining area effect figure for the boring method for establishing model horse;

Figure 13 is that face of land DEM method for establishing model is at horse drinking water mining area effect figure (1);

Figure 14 is that face of land DEM method for establishing model is at horse drinking water mining area effect figure (2);

Figure 15 is the Insert dialog box synoptic diagram;

Figure 16 is " view " toolbar synoptic diagram;

Figure 17 is " new construction " dialog box synoptic diagram;

Figure 18 is " importing the DXF file " dialog box synoptic diagram;

Figure 19 is that the section introduction method is at horse drinking water mining area effect figure;

Figure 20 is " line " properties dialog synoptic diagram;

Figure 21 is " new file " dialog box synoptic diagram;

Figure 22 is " scope 1234.STR " attribute representation intention;

Figure 23 is regular closed wire frame schematic top plan view;

Figure 24 is " seeing eastwards " synoptic diagram;

Figure 25 is " seeing eastwards " 3 d effect graph;

Figure 26 is " line " synoptic diagram;

Figure 27 is " line translation " dialog box synoptic diagram;

Figure 28 is the 3-D effect synoptic diagram after closed wire frame moves forward backward;

Figure 29 is " definition triangular mesh graph type " dialog box synoptic diagram;

Figure 30 is " selection wire box type " dialog box synoptic diagram;

Figure 31 is " network of triangle trrellis diagram attribute " dialog box synoptic diagram;

Figure 32 is a modeling procedure A synoptic diagram (1)

Figure 33 is a modeling procedure A synoptic diagram (2)

Figure 34 is " verification network of triangle trrellis diagram " dialog box synoptic diagram;

Figure 35 is " triangular mesh surface verification situation " dialog box synoptic diagram (1);

Figure 36 is " triangular mesh surface verification situation " dialog box synoptic diagram (2);

Figure 37 is the three-dimensional entity model synoptic diagram;

Figure 38 is the three-dimensional entity model of modeling scope;

Figure 39 is three-dimensional Boolean calculation synoptic diagram (the 1) (A: modeling scope of atural object entity; The B face of land);

Figure 40 is the Boolean calculation results model on the modeling scope and the face of land;

Figure 41 is three-dimensional Boolean calculation synoptic diagram (the 2) (A: the scope face of land of atural object entity; B: the stratum);

Figure 42 is the uppermost stratum three-dimensional model synoptic diagram of crop out;

Figure 43 is three-dimensional Boolean calculation synoptic diagram (the 3) (A: ore body of atural object entity; B: the grouan entity);

Figure 44 is exposed at the outer department pattern synoptic diagram of granite for ore body;

Figure 45 is following latent stratum of mining area surface and rock mass three-dimensional model;

Figure 46 is the following buried faults model of mining area surface;

Figure 47 is the following latent known oxide ore ore body model of mining area surface;

Figure 48 is " distance inverse ratio weighting " dialog box synoptic diagram;

Figure 49 is " input domain " dialog box synoptic diagram;

Figure 50 is " module file " dialog box synoptic diagram;

Figure 51 is " particular value (nonnumeric) " dialog box synoptic diagram;

Figure 52 is " search data " dialog box synoptic diagram;

Figure 53 is " producing blank block models " dialog box synoptic diagram;

Figure 54 is " limiting with the network of triangle trrellis diagram " dialog box synoptic diagram;

Figure 55 is " ascription " dialog box synoptic diagram;

Figure 56 is " calling the network of triangle trrellis diagram " dialog box synoptic diagram;

Figure 57 is " network of triangle trrellis diagram attribute " dialog box synoptic diagram;

Figure 58 is " data information " dialog box synoptic diagram;

Figure 59 is " piece definition " dialog box synoptic diagram;

Figure 60 is " new file " dialog box synoptic diagram;

Figure 61 is a fracture properties form synoptic diagram;

Figure 62 is " interpolation territory " dialog box synoptic diagram;

Figure 63 is " loading piece " dialog box synoptic diagram;

Figure 64 is regional empty piece grid model for horse absorbs water;

Figure 65 is empty piece grid model configuration of surface figure;

Figure 66 is " wire frame New anyAttribute " dialog box synoptic diagram;

Figure 67 is " wire frame ascription " dialog box synoptic diagram;

Figure 68 is " wire frame assignment " dialog box synoptic diagram;

Figure 69 is " ascription " dialog box synoptic diagram;

Figure 70 is latent stratum, horse drinking water area three-dimensional grid model synoptic diagram;

Figure 71 is regional buried faults 3D grid synoptic diagram for horse absorbs water;

Figure 72 is the latent known oxide ore ore body 3D grid synoptic diagram in horse drinking water area;

Figure 73 is the three-dimensional unusual figure of horse drinking water mining area element S n;

Figure 74 is the three-dimensional unusual figure of horse drinking water mining area elements Pb;

Figure 75 is three-dimensional unusual figure (the A:Pb element anomalies of horse drinking water mining area element S n, Pb; The B:Sn element anomalies; C: shallow; D: dark);

Figure 76 is " distance inverse ratio weighting " dialog box synoptic diagram;

Figure 77 is " input domain " dialog box synoptic diagram;

Figure 78 is " module file " dialog box synoptic diagram;

Figure 79 is " particular values (nonnumeric) " dialog box synoptic diagram;

Figure 80 is " search data " dialog box synoptic diagram;

Figure 81 is fracture buffer zone variogram (1);

Figure 82 is fracture buffer zone variogram (2);

Figure 83 is rock mass buffer variable diagram;

Figure 84 is integrated traffic unit number statistics histogram;

Figure 85 looks for the ore deposit to contain much information in 1.88 cell cube distribution plan;

Figure 86 predicts the outcome for looking for the ore deposit quantity of information;

Figure 87 is 1720 a stage casings prediction vertical view;

Figure 88 is 1620 a stage casings prediction vertical view;

Figure 89 is 1520 a stage casings prediction vertical view;

Figure 90 is study area engineering verification district and known ore body graph of a relation;

Figure 91 is study area emphasis prediction district and buried faults and known ore body graph of a relation.

Embodiment

Utilize method of the present invention that Sector East, Geju City geology is analyzed and look for the ore deposit, concrete steps are as follows:

One, geology background:

Mining area, Geju City is positioned at the southwestern margin of south China fold system on macrotectonics, the northwest is the health Yunnan earth's axis, southwest is sad labor mountain metamorphic rocks belt, northern more ancient land is faced in the southeast, important component part for many metals metallogenic belt, Nancy, the eastern regions of the Yunnan Province, be global two huge TIN ORE BELT, the i.e. joint of huge TIN ORE BELT of tethys and the huge TIN ORE BELT in Pacific Rim.

The mining area has experienced repeated tectonic movements, forms the composite factor in current tectonic structure general layout and ultra-large type tin-polymetallic mineral deposit:

(1) synthem, Lvliang City's synthem (sad firm mountain group, precious jade mountain group and bright red mountain group), Yang Zi synthem (Kunyang group) and the tin of Caledon synthem, the plumbous abundance value of preceding li west phase formation exceed 2～3 times of Clark values, constitute into the basis, ore deposit.

(2) a seal rapid sedimentation of early metaphase trough, deposited the huge thick stratum of a cover based on carbonate, simultaneously submarine volcano taking place repeatedly sprays and spoils over, not only bring a large amount of copper metallogenic materials, the block sulphuretted copper ore bed that takes on a certain scale, and seabed hot water activity thereafter brings a large amount of tin metallogenic materials from substrate, causes the high background value of tin in the stratum, even may just form " interlayer tin ore deposit " in some layer position.

(3) in the Yanshan Mountain, strong acid magma intruder operations in late period, three big remelting type granitic batholiths that the spacing that has caused the regional northern west in south, the eastern regions of the Yunnan Province approximately to be 90km on the zone distributes, and finally form ultra-large type world-famous Geju City tin ore deposit.

(4) the happiness tectonic movement of mountain since the phase makes the further lifting of the earth's crust, most interlayer ore body and part grouan contact zones ore body suffer oxidation and become loose earthy oxide ore, the face of land has formed huge placer tin mineral deposit, has also formed the spelean deposit placer and the underground streamside placer of certain scale in the ancient solution cavity.

The stratum:

Exposure stratum, mining area and to hold ore bed be system Geju City group in the Triassic system is made up of carbonate and a little petroclastic rock, contains cream salt, rich algae carbonate rock stratum for tidal flat lagoon-sabkha environment deposit.Geju City group tin average content is 5.36 * 10 ^-6, be 5 times of world's carbonatite.Group stratum, Geju City can be divided into three lithological combination sections by its lithological combination.Card room section (T ₂g ₁): thick 500-575m.Limestone and dolomitic limestone or contain the argillaceous limestone alternating layers, bedding is grown, and rich algae also has Irish touchstone, diabase to distribute.Card room section is that ore horizon is mainly composed in the mining area, and this section tax is deposited ore body and accounted for 82.5% of the ore body verified in the mining area, and the tin reserves account for 90%, and the copper reserves account for 96%, and plumbous reserves account for 44%, are the main ore-controlling layer in mining area, Geju City.Card room section is divided into six lithology layers, i.e. T again ₂g ₁ ¹, T ₂g ₁ ², T ₂g ₁ ³, T ₂g ₁ ⁴, T ₂g ₁ ⁵, T ₂g ₁ ⁶, tin ore concentrates on T especially ₂g ₁ ⁵And T ₂g ₁ ⁶, T wherein ₂g ₁ ⁶Be the best ore bed that holds in the mining area, thick 15-100m, its lithology is grey medium bed crystallite-sand formation cuttings limy dolomite and microcrystalline limestone alternating layers, sees algae bonding grumeleuse, algae group, laminated striation stone in the limy dolomite.

Horse traction lattice section (T ₂g ₂): thick 21-342m, take in pierite as the combination of main folder limestone, the alternating layers band is than T ₂g ₁Difference, for mining area second holds ore bed, the tin reserves account for 5% of mining area, and copper accounts for 1%, and lead accounts for 47%.The interlayer ore deposit is far away from T ₂g ₁Grow.

White clay hole section (T ₂g ₃): thick 28-155m, take the limestone as the combination of main folder pierite, lithology is single, the agensis of interlayer ore deposit.

Structure

East northeast-north north east is mainly to hold the ore deposit ore-controlling structure to the fold fracture in the mining area.Geju City whole district is passed through in Geju City fracture in north-south, and the mining area is divided into Sector East and Westbound two parts, and the Sector East is main many metals of tin copper metallogenic province.The wuzishan compound anticline is a mining area one-level large-scale fold structure, prolong and reach 30km, be wider than 15km, the control of thing two ends across the nearly East and West direction rift structure on the compound anticline of the Wuzhi Mountain because of ruptured by nearly north-south, caused " ladder lattice " formula, these " ladder lattice " and wuzishan compound anticline are being controlled five big ore fields of Sector East, Geju City, i.e. horse traction lattice, pine tree pin, Gao Song, Lao Chang, card room ore field.

Magmatite

The magmatite outcrop area is little in the mining area, generally be latent state output, its type has geochemistry of Yanshanian granitoids and Indo-Chinese epoch alkaline basalt, wherein but the grouan invasion is to Middle Triassic series top, be shell source remelting biotite granite, the aluminium supersaturation, the rich alkali (potassium is greater than sodium) of Silicon-rich, stanniferous (10～29) * 10 ^-6, average stanniferous 17.5 * 10 ^-6The little stock of big batholith upper process has caused anticline (arched roof), down grouan strain projection is arranged, and is best ore-forming structure magma combination in the district.Alkaline basalt is mainly composed and is stored in section bottom, group card room, Geju City (T ₂g ₁ ¹) in the stratum, by 1～2 or multilayer constitute, thick several meters to over one hundred rice, stanniferous (9～51) * 10 ⁶, the enrichment multiple is 4.7～34, visible once in a while cassiterite.Ore fields such as the best Cheng Kuang location Lao Chang in Geju City, card room, pine tree pin by chance are positioned at Geju City remelting grouan and T ₂g ₁Volcano rock stratum grow the location overlapping.

In sum, mining area overwhelming majority tin-polymetallic ore body all originates in Geju City group, and tin ore mainly is distributed in Geju City group T ₂g ₁ ⁶Section has control ore deposit, tangible stratum feature.Above-mentioned rock mass, stratum stanniferous be higher to show that granite can bring mineral resources and in the invasion process, can assimilate, extract the mineral in the country rock, then extract tin in the volcanics at grouan and volcanics overlapping, group stratum, Geju City and volcanics also play the effect of source bed respectively.In the mining area, stratum, structure, the just right Geju City ultra-large type tin ore that is configured to of magmatite provide superior control ore deposit geologic condition.

Two, utilization is looked for the ore deposit based on the method for finding mine with cubic forecast model of three-dimensional modeling

1, Data acquisition;

Mining geology, geochemistry, mine engineering data are collected by required mine engineering owned enterprise or production technology department of geologic prospecting party.Wherein, mining geology and mine engineering data obtain by collecting stage casing, mine geological map, mine centre section project arrangenent diagram and exploratory grid cross section (CAD drawing or .dwg formatted file), and go deep into mine and understand mine engineering operation situation.Wherein, stage casing, mine geological map should include particular location, form, numbering, regional tectonics (main faults) occurrence, the title of milimeter grid, exploratory line; Should include following latent formation beds form of latent absolute altitude, milimeter grid, the face of land and position in each sectional view, main faults is extended form and position, rock mass form and position, all types of known ore body form and position downwards from the face of land.Geochemical data (being the geochemical exploration data) is provided from mine element chemistry grade database (being the borehole data storehouse) by bargh or geologic prospecting party too;

2, the construction in borehole data storehouse;

Set up borehole data storehouse, mining area.Wherein, can set up mining area element anomalies three-dimensional model by the common sample resources data file in the borehole data storehouse.The borehole data storehouse is the integrated foundation of Micromine KAN TAN 3D by the research and development of Australian MICROMINE company, and can carry out unified management.

The concrete establishment step in borehole data storehouse is as follows:

(1) collect mining area borehole data (drawing or Access database), the database file that need use mainly contains: aperture resources data file, deviational survey resources data file, common sample resources data file.Software platform is: Microsoft Access.

(2) the necessary attribute of using in the resources data file of aperture has: boring numbering, three-dimensional geographic coordinate (X, Y, Z), design inclination angle, design orientation, design hole depth, whole hole depth; Other attributes comprise: unit in charge of construction, ore deposit section, stage casing, tunnel title, ore body number, classification, perforate date, whole hole date, remarks, whether the emphasis drilling depth, edit and record people, ID.

(3) the necessary attribute of using in the deviational survey resources data file has: boring numbering, the deviational survey degree of depth, inclination angle, position angle, ID.

Can in subsequent step, set up the three-dimensional model of boring form according to above three data files.

(4) the necessary attribute of using in the common sample resources data file has: boring numbering, sample number into spectrum, sampling depth, sample length, long, the required element grade of core; Other attributes comprise: ore type, remarks, three-dimensional geographic coordinate, ID.

Can in subsequent step, set up mining area element anomalies three-dimensional model according to the common sample resources data file.

The boring model is made up of 3 files of raw data base, is respectively drilling orifice three-dimensional coordinate file, boring skew back file and common sample resources data file.Study area horse drinking water area is positioned at No. 10 ore cluster south of reed dyke, and the research of this ore cluster is historical long, and the boring data with existing is more comprehensive, sets up so the foundation in borehole data storehouse can make full use of existing borehole data.Master file has: reed dyke drilling orifice data, deviational survey data, common sample resources data file.

Concrete steps are:

1) the Access database is imported among the Micromine:

File → importing → MDB (Microsoft Access) database → enter " input MDB file " dialog box is seen Fig. 1.

Click " selection database " and " option table " → respectively and select " aperture data " table, deviational survey MSDS, common sample MSDS in the same database

In " target ", be respectively aperture data, deviational survey data and common sample data " file " name, other acquiescence → inputs with the prompting that data have successfully imported occurring, are seen Fig. 2.

The right-click file name can be edited this document.

2) calibration database

Boring → verification → enter " boring verification " dialog box is seen Fig. 3.

A, " boring well head "

Double-click selection well head file → " type " selection " data " → click in " well head territory " → enter " well head territory " dialog box is seen Fig. 4.

" boring title " selects field " boring numbering "; " eastern coordinate " selects field " X "; " northern coordinate " selects field " Y "; " elevation " selects field " Z "; " total depth " selection field " whole hole depth " → close.

B, " borehole survey "

Double-click selection deviational survey file → " type " selection " data " → click in " measurement territory " → enter " measurement territory " dialog box is seen Fig. 5.

" boring " selects field " boring numbering "; " fathom " and select field " the deviational survey degree of depth "; " position angle " selects field " position angle "; " inclination angle " selection field " inclination angle " → close.

C, " between bore area "

" interval file " selects common sample data → " type " selection " data ".

Field Count → OK that right-click " interval file " title → revise → enter " file modification " dialog box → record → insertion record → enter " insert record " dialog box → selections will be inserted → get back to " file modification " dialog box → in the field of heart interpolation, insert field name, type, character duration, decimal digits etc. → close

Click " interval territory " and enter " interval territory " dialog box, see Fig. 6.

" boring name " selects field " boring numbering "; " from " selection field " sampling depth "; " to " selection " is sampled to " → closes.

D, be the report file name.

E, operation

Right-click report file → edit → enter report file editing dialog frame, what occur in " warning " territory promptly is assay.The mistake that occurs roughly has following a few class:

Lost mouth coordinate

Solution: this boring of deletion in the file of aperture;

The deviational survey data are not since 0

Solution: earlier the peaceful preface of deviational survey data boring numbering is arranged, above this holes first line data, adds data line, allow its skew back degree of depth be 0, all the other with first go identical;

Discontinuous deviational survey data

Solution: replenish the deviational survey data, make it continuous;

Lost or incorrect position angle

Solution: replenish, delete or correct bearing data;

Do not have from 0 from field value

Solution: adding 0 value from field;

From the value of field more than or equal to value to field

Solution:, be recalculated to field with the long polishing of sample;

Interval overlapping

Solution: revise artificial property improper value in the interval;

Lost from or to field value

Solution: replenish " from " or " to " field value;

The interval has surpassed hole depth

Solution: revise the artificial property mistake that occurs in the database typing;

The record of this boring is discontinuous

Solution: replenish this borehole log.

3) create new database

Dialog box that the borehole data storehouse " is created " in boring → borehole data storehouse → create → enter is seen Fig. 7.

Choose in " borehole data storehouse " → for database name → create → enter " borehole data storehouse " dialog box, see Fig. 8

" definition well head file ": double-click and select the well head file

" definition is measured ": choose " deviational survey in the well " → double-click selection deviational survey file

Click " interval file " and enter " interval file " dialog box → add → enter dialog box that " adds interval file " → double-click and select " common sample data " file → close → determine → determine, see Fig. 9.

4) load boring

Forms collection → three-dimensional → boring → " loading boring " dialog box is seen Figure 10.

A, " boring is provided with "

Choose in " borehole data storehouse " → double-click the selection database

B, drilling track show

Track rugosity: can select the drilling track rugosity as required;

Color code track: can be drilling track and give color;

Click color coding → enter " color sign indicating number " dialog box, select interval file and color gamut → close, see Figure 11.

Show the boring name: can make amendment to color

Mark well head: can carry out mark to pithead position

Determine.

The effect that the above boring method for establishing model is used in horse drinking water mining area as shown in figure 12.

3, the foundation of three-dimensional mine model;

The mine three-dimensional model that needs to set up comprises face of land DEM model, stratigraphic model, fracture model, granite model, geochemical anomaly three-dimensional model.The concrete grammar division is as follows:

(1) DEM model in the face of land is on the basis of collecting mining area surface elevation drawing (tiff form), use the MAPGIS Geographic Information System (GIS) software of China Geological Univ. Wuhan's independent research and carry out vector quantization, face of land elevation is changed into vector file, import to the foundation that realizes among the three-dimensional modeling visual software Micromine face of land three-dimensional entity model.

The effect that the above face of land DEM method for establishing model is used in horse drinking water mining area is shown in Figure 13,14.

(2) method for building up of stratigraphic model, fracture model, granite model is roughly the same.

1), collects section

At first collect stage casing, mining area geological map and exploratory grid cross section (CAD drawing, the dwg form), under AutoCAD 2006 software platforms, planimetric map collection work district section data is reconnoitred in the contrast workspace, and checks that collected section is whether corresponding with the section arranged in the planimetric map.

2), stack, adjustment section

1., in AutoCAD, open the platform master plan of certain sea level elevation.

2., in planimetric map, insert sectional view.Concrete steps are as follows:

A, click " insertion " → piece in the menu bar, enter the Insert dialog box, see Figure 15.

B, click " the browsing " of inserting in the dialog box, need to select the section that inserts, and beat " √ " in " specifying the insertion point on screen " preceding square frame, other option is given tacit consent to and is got final product click " determining ".

C, in planimetric map, specify an insertion point of inserting section, click getting final product.

3., adjust profile position.Concrete steps are as follows:

A, click " view " → 3-D view in the menu bar, select " left side is looked ".Perhaps right-click status bar blank space is chosen " view ", and " view " toolbar is accessed, and sees Figure 16.

B, in left view, erect section (be about to section half-twist).

Click left button and choose section → click right button → rotation, specify a rotation center (is example with the profile line right endpoint), click the profile line right endpoint, the input anglec of rotation (90 °) can be upright with section.

C, in vertical view, section moved on to profile line position corresponding on the platform master plan.

Left button is clicked section and it is chosen → click right button → move, and clicks profile line left end point → the click corresponding end points of identical numbering profile line on the master plan of plane.

D, section is rotated to and the identical angle of section on the planimetric map.

Left button is clicked section and it is chosen → click right button → rotation → click section right endpoint → input alphabet " R " expression reference → click section right endpoint → click section left end point to determine rotating another end points of corresponding section on initial position → click planimetric map to determine rotation final position complete operation.

E, in left view be that the line of 1800m aligns with planimetric map with sea level elevation on the section.

F, section is accurately located.

View → dynamic observe → freely dynamic observe, find in sectional view and the planimetric map can be corresponding mutually significant structure (as fracture), the intersection point of 1800m altitude line and fracture on the section is moved to the intersection point place of rupturing on sectional view 1800m altitude line and the master plan, section accurately can be located.

G, inspection section.

Get back to left view, check whether section 1800m line still overlaps with master plan, if do not overlap repeating step e～f.

H, preservation.If the 1800m line overlaps with master plan in left view, then can get back in the vertical view this section is preserved.

Clicking after choosing this section and clicking " file " → output → name → determine → get back to principal plane figure, principal plane figure can finish preservation work.

After according to above step all sections all being carried out, section is saved as AutoCADR12/LT2 DXF form.

3), divide dissection surface

All sections are all decomposed:

In vertical view, choose profile line → modification → decomposition.How to decompose several times, till can not decomposing again.

Then file is preserved:

File → save as saves as AutoCAD R12/LT2 DXF form → preservation with file.

4), import Micromine

Under Micromine KAN TAN 3D software platform, section file that will pre-service is good in AutoCAD imports, and realizes the geologic body d solid modeling in this software.

1., a newly-built engineering

Click " file " → engineering → newly-built → enter " new construction dialog box "

The input engineering name is also set engineered paths, fills in title, and the graphic file that generates in Micromine later will all be saved in this document folder, see Figure 17.

2., section is imported in the engineering

" √ " → move beaten in " incoming line " front of clicking file → importing → DXF file → enter dialog box → input file path that " imports the DXF file " → in " output ", sees Figure 18.

3., according to step method 2. all sections are all imported among the Micromine.

The effect that the above section introduction method is used in horse drinking water mining area as shown in figure 19.

5), editor's section

1., access section

Forms collection → View Manager → double-click " line ", enter " line " dialog box → in " file " of display file, fill in section file name or double-click browser document select → type selecting " line " → eastern territory fills in " EAST ", northern territory and fills in " NORTH ", Z territory and fill in " RL ", link field and fill in " JOINT " → determine, can access this section file, see Figure 20.

2., delimit the modeling scope

A, newly-built scope file.

After all sections are all accessed, the white bars frame of forms collection → View Manager → line → enter " line " dialog box → right-click " file " back → newly-built → enter " new file " dialog box, as shown in figure 21.

" file " located to fill in scope name, " type " and located to select " line " → determine → enter ew line attributes edit dialog box.Fill in the title of eastern territory, northern territory, Z territory, line territory, link field in " field name " hurdle respectively, type selecting N---numeric type is set live width and decimal digits, sees Figure 22.

Close → system prompt " whether newly-built line file " → determine → enter file attribute editing dialog frame → close.Enter " line " dialog box after closing, click OK adjustablely goes out newly-built scope file.

B, draw the line

Display panel → View Manager → right-click scope file → editor, the scope file of making enters editing mode.

Access all section files, click " newly-built line " icon in the toolbar → regular closed wire frame smaller slightly of picture, see Figure 23 than section scope.

Click and " see eastwards " in the view's toolbar that icon → scope is duplicated translation downwards, and (see that eastwards design sketch sees Figure 24,3 d effect graph is seen Figure 25.Display panel → View Manager → right-click scope file → renewal is preserved the scope file.

3., editor's section

A, enter editing mode:

Section sequence number → editor that display panel → right-click will be edited.Click " section " icon in the toolbar

→ along section with standardized line of cross star → enter editing mode that conforms to profile position.

B, line

Need to determine the line drawing shown in the atural object (as stratum, fracture, ore body etc.) of modeling → click " newly-built line " icon in the toolbar → on the section to go out the atural object profile.

To defer to following principles during setting-out:

On the line too many point is not arranged in the time of line, in order to avoid mistake when next step three-dimensional modeling, occurs.

The place that atural object intersects, line will be drawn longerly.Fracture and the granitic mass place of intersecting for example will be with the line drawing of rupturing in granitic mass, as shown in figure 26.

For single ore body, ore body is drawn as the wire frame of a sealing, and with line segment translation one segment distance forward and backward respectively.Concrete operations are as follows: after the profile of this ore body is described, click left button and choose this ore body → click right button → move → enter " line translation " dialog box → selection " vertical range " pattern, and in the square frame of " duplicating " front, beat " √ " → fill in the distance (moving forward) that will move → determine, see Figure 27 on the occasion of moving backward, duplicating.

3-D effect after moving as shown in figure 28.

Beyond inc atural object outline line (stratum, rock mass, fracture, known ore body etc.) the modeling scope of as far as possible all drawing.Carry out rational Boolean calculation when being convenient to three-dimensional modeling.

According to the above four principle, delineated the atural object profile of the need modeling on all sections.

6), three-dimensional modeling

Forms collection → three-dimensional

1., newly-built type.

Wire frame → type → enter " definition triangular mesh graph type " dialog box → newly-built → enter " newtype " dialog box → fill in newtype title (as stratum, fracture, grouan, level line etc.) → determine can be set up a newtype, sees Figure 29.

2., newly-built wire frame

Forms collection → three-dimensional → line accesses a section.

Wire frame → newly-built promptly can be at the new wire frame of display panel → sky of " three-dimensional " lining appearance.Double-click this wire frame → enter " selection wire box type " dialog box, see Figure 30.

Select a wire frame type → select → enter " network of triangle trrellis diagram attribute " dialog box, fill in title, code, title, select color → determine, see Figure 31.

At this moment get back to principal plane figure, can edit wire frame.

3., modeling

A, forms collection → three-dimensional → line can access profile line.In each section, seek corresponding line, click and connect, see Figure 32,33.

B, verification

After wire frame connected, click right button → verification → enter " verification network of triangle trrellis diagram " dialog box, and beat " √ " in the square frame of " inspection " in " compression point " in " optimization option " and " compression triangle " and " the verification option " and " inspection closure " front from intersecting, activate color can be set to more eye-catching voluntarily and color that contrast is bigger → determine, see Figure 34.

Can occur triangular mesh surface verification situation after determining,, generally can point out 1 ~ 2 opening for inc surfaces such as stratum, large-scale fractures; For closure surfaces such as single ore bodies, two openings generally can have been pointed out.It should be noted that necessarily can not have and intersect triangle and occur,, will make amendment by the diabolo surface mesh, and then verification,, see Figure 35 up to not intersecting till the triangle if intersect triangle.

For single ore body, need seal its surface, make it become the body of a sealing.Concrete steps are: click the triangular mesh surface opening that verification comes out → click right button → sealing.By the time after opening all seals, also to carry out verification,, see Figure 36 up to both not having invalid connection and opening also not to intersect till the triangle.Check errorless after, click right button → withdraw from editing machine, the selection wire box type is inserted title, code and title in the triangular mesh properties dialog, click the three-dimensional entity model of determining to obtain atural object.Setting up good entity three-dimensional model can access by forms collection → three-dimensional → network of triangle trrellis diagram, sees Figure 37.

4., set up the three-dimensional stereo model of modeling scope

According to the method that the atural object three-dimensional model is set up, set up the three-dimensional stereo model of a modeling scope, see Figure 38.

5., the three-dimensional Boolean calculation of atural object entity

Boolean (Boolean calculation) is by plural object is carried out the computing of union, difference set, common factor, thereby obtains a kind of operational method of new object form.

4 kinds of boolean operation methods are provided among the Micromine KAN TAN 3D.Be respectively that A hands over B (Ain B), A to subtract B (A out B), B to hand over A (B in A), B to subtract A (B out A).The logical meaning of their representatives is respectively:

A hands over B (A in B): output A hands over the part in B;

A subtracts B (A out B): the part of output A outside B;

B hands over A (B in A): output B hands over the part in A;

B subtracts A (B out A): the part of output B outside A.

More than four kinds of fundamental operation methods more can be incorporated into more complex calculations.For example:

A, A represent the modeling scope, and B represents the face of land, see Figure 39 (operational method: A hands over B, B to hand over A).

The model that obtains is the stack model on the modeling scope and the face of land, sees Figure 40.

B, A represent the stack model (promptly last figure) on the modeling scope and the face of land, and the lower surface on first stratum that B representative and surface distance are nearest is seen Figure 41 (operational method is: A subtracts B, B and hands over A).

The model that obtains is: the uppermost stratum of crop out, see Figure 42.

The ore body of c, A representative and grouan meet, B represents the grouan solid model, sees Figure 43 (operational method is: A subtracts B, B and hands over A).The model that obtains is that ore body is exposed at the outer department pattern of granite, sees Figure 44.

The effect that the above three-dimensional mine model method for building up is used in horse drinking water mining area such as Figure 45 (1.-T ₂g ₂ ³2.-T ₂g ₂ ²3.-T ₂g ₂ ¹4.-T ₂g ₁ ⁶5.-T ₂g ₁ ⁵6.-grouan; ), Figure 46 (1.-big Qing Dong fracture; 2.-fracture of kylin mountain; 3.-fracture of reed dyke; 4.-No. 1 fracture; 5.-Da Qing South Mountain fracture), Figure 47 (1.-steep dip rake shape ore body; 2.-moderate dip is like bank of ore) shown in, be followed successively by: following latent stratum of mining area surface and rock mass three-dimensional model, fracture model, known oxide ore ore body model.

Geochemical exploration modelling concrete steps are as follows:

A, set up study area boring solid model according to the method for building up of aforementioned boring model.

B, set up the geochemical exploration model according to the GEOCHEMICAL DATA file that provides in the borehole data storehouse

Figure 48 is seen in modeling → three-dimensional nugget estimation → distance inverse ratio weighting → enter " distance inverse ratio weighting " dialog box.

A, input:

Double-click selection sample analysis file name → " type " selection " data " → click " input domain " → enter " input domain " dialog box is seen Figure 49.

In must filling out territory " input domain 1 ", double-click the element that selection will show

Double-click and select " eastern territory ", " northern territory ", " elevation "

Close

B, from the document definition piece

Choose " from the document definition piece "

Click " more " → enter " module file " and see Figure 50

Documentation requirements: double-click and select empty block models file → " type " selection " data "

The module territory is set: double-click and insert " east ", " north ", " elevation "

Close

C, particular values

Click " particular values (nonnumeric) " → enter " particular values (nonnumeric) " dialog box, see Figure 51.

Choose " ignoring character " and " skipping blanks " → close.

D, model

" method " selection " distance inverse ratio " → " data search " selection engineering nugget model file → right-click this document name → enter " search data " dialog box is seen Figure 52.

" radius " inserts interpolation search ellipsoid radius value → selection sector → determine.

E, output

" video data " selection " nothing "

Select output file

" type " selection " data "

F, click " OK " and carry out.

Use the interpretation of result after modeling is carried out in Sector East, Geju City, Yunnan horse drinking water mining area of above-mentioned modeling method:

The mining area of this task is positioned at the big cyanines in high loose ore field middle part-A Xi stockaded village to the oblique wing one side by north, link to each other with reed dyke ore cluster with big cyanines east, this area ore cluster, the face of land does not have the exposure of magma rock mass in the district, institute exposure stratum is system Geju City group carbonate formation in the Triassic system, rift structure is extremely grown, the location of ore body and output form are subjected to formation lithology and architectonic control fairly obvious, be mining area, Geju City to tiltedly and east northeast jointly control the interlevel oxidation mineral deposit and be one of special location of multi-level intensive output to rift structure.

The stratum that distributes in the mining area is system Geju City group card room section and horse traction lattice section (T in the Triassic system ₂g ₁ ⁵-T ₂g ₂ ³) carbonate formation, it is different that different formation lithologies plays the mineralization intensity of ore-forming element.Its rock signature is described below:

T ₂g ₂ ³: press from both sides a small amount of dolomitic limestone and limestone phacoid to thick-layer shape pierite in grey, the light gray, microbedding haircut is educated, and the part contains the crinoid stem fossil, thick 97-440 rice.

T ₂g ₂ ²: grey medium bed shape contains limy dolomite, and dolomitic limestone and limestone phacoid that interrupted folder differs in size contain several layers of crinoid stem fossil, thick 301.8 meters.

T ₂g ₂ ¹: Dark grey, grey thick-layer shape pierite, face of land mineralising is more weak, and the following mineralising of 300-400 rice strengthens under the face of land, and compose in local location has the rich lead ore body, thick 189.6 meters.

T ₂g ₁ ⁶: grey, light grey medium bed shape limestone and limy dolomite alternating layers, variation in thickness is big, for the surplus rice of 41.5-200, is the main ore--bearing horizon of high loose ore field.Because the two chemistry, physical property of limestone and limy dolomite alternating layers is very greatly different, and be alternately distributed, under the tectonic stress effect, easily produce splitting and interlayer fragmentation, form good space, storage ore deposit, so the mineralization intensity of tin, copper, lead, zinc, silver is higher than other layer position significantly in this layer, the overwhelming majority of verifying in the mining area promptly results from this layer like layered minerals ore body and veiny ore deposit shape.

T ₂g ₁ ⁵: grey, light grey medium bed shape ls contain shale, tool corrugated and tiger fur columnar structure.Top folder layer 2-3 ribbon limy dolomite extends stablely, and the bottom is for closing chert nodule thin layer limestone, and accidental have gastropod, lamellibranch, crinoid stem fossil, thick 336-662.7 rice.This floor is for another composes ore horizon in the district, and output vein orebody has been main, has a small amount of like the bank of ore output.

The fracture model analysis result:

Fold and fracture army grow in the mining area, and bending phenomenon is obvious, but the joint that is mostly rift structure and derivation thereof, crack that the output and the control of ore body is played a major role.

Rift structure is extremely grown in the loose ore field of whole height, and is crisscross, on the basis of oblique structure, forms the tectonic framework that seems to be " chessboard " form on big cyanines-A Xi stockaded village, with the mineralization relation significantly mainly be east northeast side and nearly west to two groups of fractures.Main faults is configured to the fracture of kylin mountain, the fracture of stove dyke, a fracture and big Jing Dong fracture in this study area, and each fracture characteristic is summarized as follows:

1. big Jing Dong fracture: northwest (NW) moves towards fracture, be distributed in study area the north, in whole mining area, move towards long 2-3 km, broken bandwidth 5-10 rice, the part reaches 30 meters, be northern 40-50 degree west trend, east northeast tilts, inclination angle 66-88 degree, and the dust rock zone is obvious, visible broken porphyry of while and particle rock, the weak haematization of face of land tool.Past its trend of east southeast of this fracture fades to East and West direction and lies concealed in underground, and face of land sign is not obvious, and ore-bearing potential is strengthened gradually downwards, and more obvious below 1720 stage casings, the lode tin ore body of loose ore deposit maximum promptly results from this No. 131.

2., kylin mountain fracture: east-west single-row, come across high loose ore field middle part, trend is grown up in 6 kms, and oblique extension is greater than 800 meters, broken bandwidth 6-30 rice, the southern side has several parallel shatter belts to produce fracture tendency east northeast, inclination angle 70-83 degree, this obvious bad break east northeast that ruptures is to fracture, in its deep mineralising crescendo, the part has the tin ore bodies to deposit, and it is positioned at the north of study area.

3. reed dyke fracture: east northeast move towards fracture, for the trunk structure and derive from a series of secondary fractures and crack, as a fracture etc., is the main structure that local area is controlled the ore body output.High loose ore field middle part is tiltedly worn in the fracture of reed dyke, be northern 36-45 degree east trend, tilt to northwest (NW), inclination angle 46-88 degree, trend is grown up in 8 kms, broken bandwidth 5-30 rice, and the part reaches 50-60 rice, form by pressure break rock, crushed rock, broken porphyry, particle rock, mylonite and rubblerock, show as the cleavageization toward the deep extension.This fracture is derived to have and is lain concealed in underground parallel with it or be the secondary fracture of east northeast direction, general and reed dyke fracture formation " going into " font on section of parallel output.

4. a fracture: be positioned on the reed dyke fracture northwest (NW) side and coil, derive from the buried faults of branch for the fracture of reed dyke, intersect gradually with the fracture of reed dyke from north orientation south on the plane, separate gradually with the fracture of reed dyke from top to bottom on the section, its feature is similar to the fracture of reed dyke, become northern 25-35 degree east trend, the tendency northwest (NW), inclination angle 45-80 degree, trend is grown up in 2 kms, tilt to extend greater than 500 meters broken bandwidth 0.2-10 rice.This fracture constitutes this Qu Chengkuang with the fracture of reed dyke, and a structure typical case in control ore deposit promptly at the retaining part of fracture, is controlling well stratiform, veiny, the mutual shape ore body of layer arteries and veins output.In the secondary little fracture of above-mentioned two fracture sides, Fu Xi vein orebody output is also arranged.

The granite model result is analyzed: the face of land does not have the magmatite exposure in the scope of mining area, but there is latent granitic mass to distribute in the deep, according to the latent grouan solid model in the mine of building up, the latent table of this district's granite is high about 900-1200 rice, apart from face of land 1200-1400 rice, this study area is in the depression position of granitic mass, its configuration of surface be east northeast to groove, three in east southeast west is higher, and the north is lower; This rock mass lithology is similar to old card rock mass, is middle particulate biotite granite, and this granite is the parent rock body that becomes the ore deposit in the district.

4, the foundation of cubic forecast model;

1) extraction of geologic body three-dimensional cube

According to the announcement of geologic information to ore body, the particularly distribution of exploratory line, determine the scope and the basic parameter of modeling in conjunction with form, trend, tendency and the spatial distribution characteristic of ore body, the scope model is divided into the cell block of row * row * layer for 20m * 20m * 20m, adopt time piecemeal in edge of model, its specification is 10m * 10m * 10m.

Behind the three-dimensional grid model of setting up, can will look for digital model determined Prediction Parameters in ore deposit to compose to each cell block as attribute.Use the stratum solid model that 3D grid is limited, mark off the grid cell that Different Strata comprises, as the lithology variable in the prediction of mineral deposit; Use the fracture solid model that 3D grid is limited, mark off the grid cell that different fractures are comprised, as the fracture variable in the prediction of mineral deposit.Use known ore body solid model that 3D grid is limited, mark off the grid cell that different ore body comprises, as the priori conditions in the prediction of mineral deposit.

Three-dimensional cube extracts concrete steps:

1. set up empty scope nugget model:

Dialog box that the new nugget model of modeling → measurements of the chest, waist and hips nugget estimation → establishment → enter " produces blank block models " is seen Figure 53.

A, choose " limiting empty block models " and " network of triangle trrellis diagram " → click " more " → enter " limiting " dialog box, see Figure 54 with the network of triangle trrellis diagram.

Select " single ", the attribute of single predictive variable can be composed to cell block; Select " storehouse " the storehouse attribute can be composed and enter " ascription " dialog box → close, see Figure 55 to cell block → click " ascription ".

" limit " dialog box → in " block models ", choose " inferior piecemeal " with the network of triangle trrellis diagram

" limit " regulation of dialog box → in " inferior blocking " " east time block count ", " north time block count " and " the inferior block count of elevation " → close with the network of triangle trrellis diagram.

B, piece definition

Forms collection → three-dimensional → network of triangle trrellis diagram → empty the block models of the scope of opening

Dialog box that this document → attribute in the right-click display panel → enter " calls the network of triangle trrellis diagram ", see Figure 56 → right-click title → enter " network of triangle trrellis diagram attribute " dialog box, see and Figure 57 → click " data information " → enter " data information " dialog box see Figure 58.

" X minimum ", " Y minimum ", " Z minimum ", " X maximum ", " Y maximum ", " Z maximum " in the record " limit ".

In " producing blank block models " dialog box, click " piece definition " → enter " piece definition " dialog box, see Figure 59 → insert respectively " " east " of starting block " center " (X minimum), " north " (Y minimum), " elevation " (Z minimum) coordinate and " stop " east " (X maximum), " north " (Y maximum), " elevation " (Z maximum) coordinate of piece " center ", insert cell block center " spacing " → click " piece number " system with automatic computing block number → close.

C, output

Right-click " file " → newly-built → enter " new file " dialog box, see Figure 60 → insert file name → selection type (data) → can select to use masterplate → determine → enter " new file editor " dialog box, can edit field name, type, character duration, decimal digits, see Figure 61.

" produce blank block models " dialog box → insert " eastern territory " (EAST), " northern territory " (NORTH) and " elevation " (RL).

Figure 62 → can make amendment to field name, type, character duration, decimal digits in the neofield is seen in " produce blank block models " dialog box → " interpolation territory " → enter " interpolation territory " dialog box.

Operation.

2., load piece

Forms collection → three-dimensional → nugget model → enter " loading piece " dialog box is seen Figure 63.

Right-click " nugget model file " → edit → enter nugget model file database can be checked and edits → close this database

" loading piece " dialog box → selection type (data) → insert " eastern territory " (EAST), " northern territory " (NORTH) and " elevation " (RL) → determine.

The effect that the method for the empty block models of above-mentioned foundation is used in horse drinking water mining area is shown in Figure 64,65;

3., give the cell block ascription

A, see Figure 66 for wire frame New anyAttribute (is example with the stratum solid model of having built up).

Type under the layer entity model of wire frame → type → selectively → editor → newly-built → enter " NewAttribute " dialog box → input attributes title (can set a plurality of attributes of different purposes simultaneously---set storehouse attribute and layer attribute in the work) → determine → close → close.

B, be the wire frame ascription, see Figure 67.

A stratum solid model file of forms collection → three-dimensional → network of triangle trrellis diagram → open.

This stratum solid model file → attribute of right-click in the control panel → enter " calling the network of triangle trrellis diagram " dialog box → right-click entity stratum title → enter " network of triangle trrellis diagram attribute " dialog box → click " user-defined characteristic " → enter " user-defined network of triangle trrellis diagram attribute " dialog box → behind Property Name, insert respectively property value → close → determine → determine.

C, be the cell block ascription

Modeling → assignment → network of triangle trrellis diagram → enter " wire frame assignment " dialog box is seen Figure 68.

" input " → type → selection " block models "

The empty block models of " file " range of choice → insert " type " (data), " eastern territory " (EAST), " northern territory " (NORTH), " elevation " (RL).

Can select in " network of triangle trrellis diagram " to give cell block ascription → click ascription " more " → enter " ascription " dialog box, see that Figure 69 → selections is with the empty block models of " attribute " of wire frame " compose to " with single wire frame or with the form in storehouse.

Select " rewriting aiming field "

" block models " selection " blocking factor " → input " east time block count " value, " north time block count " value and " elevation time block count " value → input " report file " title → operation.

The above be the effect used in horse drinking water area of the method for cell block ascription shown in Figure 70,71,72,73,74,75, i.e. the three-dimensional unusual figure of the three-dimensional unusual figure of stratum three-dimensional grid model, fracture three-dimensional grid model, ore body three-dimensional grid model, horse drinking water mining area element S n, horse drinking water mining area elements Pb, the three-dimensional unusual stacking diagram of horse drinking water mining area element S n, Pb.

After setting up this study area three-dimensional grid model, can utilize common sample data spare in the borehole data storehouse to analyze three-dimensional unusual distribution of element of these cell blocks.Required constituent content to sample in the boring is analyzed, and the weighting of service range inverse ratio is carried out required constituent content interpolation to the unknown elements piece.

Analytical procedure is as follows:

1. the statistical study of grade data: common sample data statistics is analyzed the mathematical distribution rule that obtains the required element in this area, and find extra-high-speed value, minimum value, mean value and the variance of these elements.

2. the combined treatment of sample: the back of averaging of each sample grade weighting in the combined segment is formed new boring.Getting sample length is 0.3, cutoff grade 0.2, and the extra-high-speed value gets 10%.

3. limit boring with the empty block models file of study area scope.

4. the length combination of sample: for the ease of interpolation, decompose the sample that sample is a regular length again, get that the combination sample is long to be 10m.

5. sample interpolation: with the sample analysis data after the length combination is source data, service range inverse ratio weighted method, it is 2 that the power inverse is set, minor increment 5m, the ellipsoid that definition search is oval, 50m successively, 100,150m, 200m are that radius search known sample point comes the required element of unknown nugget is carried out interpolation.

The three-dimensional unusual figure of analytical element, and the anomaly threshold of required element is carried out value according to the unusual figure of three-dimensional.

Element anomalies modelling concrete steps are as follows:

A, set up study area boring solid model according to the method for building up of aforementioned boring model.

B, set up the element anomalies model

Figure 76 is seen in modeling → three-dimensional nugget estimation → distance inverse ratio weighting → enter " distance inverse ratio weighting " dialog box.

A, input:

Double-click selection sample analysis file name → " type " selection " data " → click " input domain " → enter " input domain " dialog box is seen Figure 77.

In must filling out territory " input domain 1 ", double-click the element that selection will show

Double-click and select " eastern territory ", " northern territory ", " elevation "

Close

B, from the document definition piece

Choose " from the document definition piece "

Click " more " → enter " module file " and see Figure 78.

Documentation requirements: double-click and select empty block models file → " type " selection " data "

The module territory is set: double-click and insert " east ", " north ", " elevation "

Close

C, particular values

Click " particular values (nonnumeric) " → enter " particular values (nonnumeric) " dialog box, see Figure 79.

Choose " ignoring character " and " skipping blanks " → close.

D, model

" method " selection " distance inverse ratio " → " data search " selection engineering nugget model file → right-click this document name → enter " search data " dialog box is seen Figure 80.

" radius " inserts interpolation search ellipsoid radius value → selection sector → determine.

E, output

" video data " selection " nothing "

Select output file

" type " selection " data "

F, click " OK " and carry out.

2) calculating of geological prospecting variable and extraction

According to the data that study area was collected, confirmable types of variables comprises stratum, fracture, rock mass and element anomalies, respectively this four class variable is carried out the division of cubic units and extracts to such an extent that step is as follows.

Grasp the situation on following latent stratum, the study area face of land, with all latent stratum respectively as a variable;

Getting rupture surface both sides 150m (this value is decided on different deposit types) is buffer zone, sets up fracture buffer zone variable.

With the above 150m of rock mass end face (this value is decided on different deposit types) is buffer zone, sets up the rock mass buffer variable.

In the boring model the required research constituent content of sample is analyzed, obtained the constituent content of each cubic units, all are studied elements respectively as a predictive variable.

Horse absorbs water regional stratum T ₂g ₂ ³, T ₂g ₂ ², T ₂g ₂ ¹, T ₂g ₁ ⁶, T ₂g ₁ ⁵, the stratum of five groups is respectively as 5 variablees;

Getting rupture surface both sides 150m is buffer zone, sets up fracture buffer zone variable, shown in Figure 81,82.

Above 150m is a buffer zone with the rock mass end face, sets up the rock mass buffer variable, shown in Figure 83.

Boring is analyzed Sn, the Pb content of sample in the model, obtains the constituent content of each cubic units, with Sn, Pb as two predictive variables.

Behind the three-dimensional grid model of setting up, can will look for digital model determined Prediction Parameters in ore deposit to compose to each cell block as attribute.

Use the stratum solid model that 3D grid is limited, the grid cell that division place Different Strata is comprised is as the lithology variable in the prediction of mineral deposit.

Use the fracture solid model that 3D grid is limited, mark off the grid cell that different fractures are comprised, as the fracture variable in the prediction of mineral deposit.

Use known ore body solid model that 3D grid is limited, mark off the grid cell that different ore body comprises, as the priori conditions in the prediction of mineral deposit.

3) the collective geological prospecting information of calculating each geologic unit and being comprised

When calculating the collective geological prospecting information that each geologic unit comprised, the software platform of using is the Access database software in Microsoft Office 2003 software packages, to import among the Access by the autonomous integrated database of Micromine KAN TAN 3D, Access provides the several data statistical query that comprises SQL ad hoc inquiry function, can satisfy worker's different statistical demand.Concrete steps are as follows:

1., scrutinizing on main ore-controlling factors and the favourable indicator for deposit basis, choose geology, the element anomalies sign of statistical study variable, and each sign of agreement to have value in the unit be 1, do not exist and get 0, each sign of statistics is in the distribution of each unit.In the computation process, the high size of nugget length and width is unified for 10m * 10m * 10m (alleged unit sizes promptly), calculate then.

2., look for the ore deposit quantity of information to calculate

The ore deposit quantity of information of looking for of certain indicator for deposit is calculated with conditional probability, that is:

$I_{A\left(B\right)}=\lg \frac{P(A/B)}{P\left(A\right)}---\left(1\right)$

I in the formula _{A (B)}The quantity of information that the B ore deposit is arranged for the A sign; P (A/B) is for the probability of A occurring under the existence condition of the known B of having ore deposit; P (A) is for the probability of sign A occurring in study area.

Because the difficulty on the probability estimate is come the estimated probability value with frequency values.At this moment

$I_{A\left(B\right)}=\lg \frac{\left(\frac{N_j}{N}\right)}{\left(\frac{S_j}{S}\right)}---\left(2\right)$

N in the formula _jContain the ore deposit unit number for what have sign A in the study area; N contains the ore deposit unit number in the study area; S _jFor having the unit number of sign A in the study area; S is the unit sum of study area.

General Two-Dimensional looks for ore deposit quantity of information method to determine to look for the favourable sign in ore deposit with formula 3.

${\mathrm{AI}}^{+}=K\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{I}_j---\left(3\right)$

Under the situation of variable negligible amounts, just do not calculate the accumulative total critical value of useful information, but directly use all on the occasion of the information number.

Utilize each indicator for deposit quantity of information to calculate the quantity of information summation of each elementary cell respectively, determine to look for ore deposit quantity of information critical value according to the subjective probability method then, the distant view unit of determining, ore deposit that contains comprises known ore body unit, contains therefore that the prospecting evaluation probability is in the distant view unit, ore deposit: contain the known ore body unit number/known ore body unit number that comprises in the distant view unit, ore deposit

On the basis of scrutinizing main empty ore deposit factor and favourable indicator for deposit, choose 5 of geology, the element anomalies signs of statistical study variable, and the agreement case marker to aim at existing in the unit value be 1, not having value is 0, each sign of statistics is in the distribution of each unit, that calculates each sign looks for the ore deposit quantity of information, the results are shown in following table.

Look for ore deposit quantity of information result of calculation table

	Nj	Sj	(Nj/n)/(Sj/S)	Ij
					T 2g 2 1	66	231609	0.31801029	-0.497558827
T 2g 1 6	753903	322611	2.604764557	0.415768474
					T 2g 1 5	903	623846	1.615336567	0.208263025
Fracture	1788	1098987	1.815631163	0.259027628
					Rock mass	30	238356	0.140458438	-0.852452165

Sn is unusual	1068	141085	8.447789413	0.926743079
					Pb is unusual	348	203524	1.908166572	0.280616283

4) the ore deposit probability calculation is distinguished in prediction

With a regional classification is some unit, has four kinds of situations:

(1) blank; (2) x unit, expression have certain geology sign to exist; (3) xy unit, expression have this kind geology sign to have the ore deposit to exist again; (4) y unit, expression have the ore deposit to exist, but do not have this geology sign.Corresponding frequency is: f ₀, f _x, f _Xv, f _y

Appoint in the unit that has certain geology sign to exist and get 1, it may contain the probability P in 1 mineral deposit _xFor:

$P_x=\frac{f_{\mathrm{xy}}}{(f_x+f_{\mathrm{xy}})}---\left(4\right)$

Appoint in the unit that does not have certain geology sign to exist and get 1, it may contain the probability P in 1 mineral deposit ₀For:

${\mathrm{<figure-callout\; id="0"\; label="P"\; filenames="C200710098940E00561.png,C200710098940E00802.png"\; state="\{\{state\}\}">P</figure-callout>}}_0=\frac{f_y}{(f_y+f_0)}---\left(5\right)$

Setting looks for ore deposit informational capacity lower limit as the comprehensive geology sign, utilizes probability estimate method Estimation Study district ore body unit output number.Calculate f _x+ f _XyAnd f _Xy, and obtain P according to (4) formula _xCalculate f ₀+ f _yAnd f _y, and obtain P according to (5) formula ₀

Usually the calculating here not have to consider beyond known mineral deposit the also mineral deposit of possible output, needs the control zone of selecting working level high for this reason, obtains a correction factor, is used for revising P _xAnd P ₀, get P _x' and P ₀', be used for the prediction district with revised probability, obtain predicting the estimation of district's mineral deposit number.

The study area of selecting certain elevation interval is that the control zone is revised, unit, calculation control district sum and gross information content are positioned at the unit number on the lower limit, must correction coefficient be: the control zone ore body unit number of the known ore body unit number in control zone/obtained by the output probability that calculates above.

Proofread and correct the probability P after obtaining proofreading and correct with correction coefficient _x' and P ₀', f _x* P _x'+f ₀* P ₀' can be used for estimating the ore body unit output number in distant view prediction district.

The ore deposit quantity of information of looking for of certain sign is calculated with conditional probability in the table, that is:

I _A(B)＝lgP(A/B)/P(A)

Several I in the formula _{A (B)}The quantity of information that the B ore deposit is arranged for the A sign; P (A/B) is for the probability of A occurring under the existence condition of the known B of having ore deposit; P (A) is for the probability of sign A occurring in study area.

Because the difficulty on the probability estimate is come the estimated probability value with frequency values.At this moment,

I _A(B)＝lg(Nj/N)/(Sj/S)

In the formula Nj be have sign A in the study area contain the ore deposit unit number; N contains the ore deposit unit number in the study area; Sj is the unit number that has sign A in the study area; S is the unit sum of study area.

In the computation process, the high size unification of nugget length and width is 10m * 10m * 10m (promptly calling unit sizes in the following text), calculates then.The study area cube adds up to 1921701 units, and comprise the nugget cube and add up to 1722 unit,

Forecasting institute is as follows with variable and statistics thereof:

1. T ₂g ₂ ³: contain 0 unit in ore deposit, shared total cube number is 363022 units.

2. T ₂g ₂ ²: contain 0 unit of ore deposit cube, shared total cube number is 380613 units.

3. T ₂g ₂ ¹: contain 66 units of ore deposit cube, shared total cube number is 231609 units.

4. T ₂g ₁ ⁶: contain 753 units of ore deposit cube, shared total cube number is 322611 units.

5. T ₂g ₁ ⁵: contain 903 units of ore deposit cube, shared total cube number is 623846 units.

6. fracture: contain 1722 units of ore deposit cube, shared total cube number is 1098987 units.

7. rock mass: contain 30 units of ore deposit cube, shared total cube number is 238356 units.

8. Sn unusual (non-null value): the result of laboratory test of normal state statistics Sn is an average 0.22999, variance 0.0746407, and standard deviation 0.2732044 0.18 to be anomaly threshold, then contains 1068 units of ore deposit cube, 141085 units of shared total cube.

Pb unusual (non-null value): the result of laboratory test of normal state statistics Pb is an average 0.15698, variance 0.0243832, and standard deviation 0.1561512 0.1 to be anomaly threshold, then contains 348 units of ore deposit cube, 203524 units of shared total cube.

5, the checking that predicts the outcome.

The effect that the above Forecasting Methodology is used in horse drinking water mining area is as follows:

Utilize each indicator for deposit quantity of information to calculate the quantity of information summation of each elementary cell respectively, determine to look for ore deposit quantity of information critical value to be according to the subjective probability method then 〉=0.941 o'clock, the distant view unit of determining, ore deposit that contains comprises 1068 unit of known ore body, therefore account for the 1068/1722=62.02% of known ore body, contain that the prospecting evaluation probability is 62.02% in the distant view unit, ore deposit.

Figure 84 is the statistic of classification of integrated traffic unit number, calculates integrated traffic 〉=0.67 cell cube accounts for 19.7% of total unit number, and 0.95 unit number accounts for 10.06% of total unit number.To be integrated traffic distribute greater than 1.88 cell cube Figure 85, can see mainly being distributed in the northern position near reed dyke ore cluster of study area, more identical with actual conditions.Figure 86 predicts the outcome for looking for the ore deposit quantity of information to send out.Figure 87, Figure 88, Figure 89 are respectively the prediction sectional views in 1720,1620,1,520 three stage casings.

The engineering verification district

According to the element anomalies result of study area, determine three engineering verification district I, II, the III of study area, see Figure 90, coordinate is respectively X=124820, Y=5581700, Z=1650; X=124890, Y=5581820, Z=1650; X=125200, Y=5582000, Z=1680.

The cube piece number in checking I district is 373, and the prediction amount of metal is 53864.99 tons.

The cube piece number in checking II district is 414, and the prediction amount of metal is 48374.16 tons.

The cube piece number in checking III district is 127, and the prediction amount of metal is 7157.49 tons.

The total syncube number of engineering verification Qu is 915, cumulative volume 3536000m3, and with oxidized ore average weight 2.451, the prediction total metal content that calculates engineering verification district Sn is 109396.63 tons.

Emphasis prediction district

According to predicting the outcome among quantity of information result of calculation and Fig. 8, determine to look for the ore deposit quantity of information=0.941 legislation body is as study area emphasis prediction district, its coordinate range such as Figure 91.As can be seen, the position in prediction district is subjected to the strictness control of stratum and fracture

According to above result, can determine that stratum and fracture are the main control ore deposit factors of study area, therefore (do not contain T with fracture 150m buffering ₂g ₁ ⁶, T ₂g ₁ ⁵), T ₂g ₁ ⁶, T ₅g ₁ ⁵The stratum scope obtains predicting that the total unit number of prospective area is 1634576, wherein with quantity of information as distant view prediction district 〉=0.941 as the comprehensive geology sign, be masked as 1 unit number f _x=164697, be masked as 0 unit number f ₀=1469879, in conjunction with the P that determines _x' and P ₀', the prediction ore body element number that can obtain study area prediction prospective area is 2324, known study area is found 1788 of ore body unit number, therefore also has 536 ore body unit not to be found.

Table 1 distant view prediction district ore body unit number

After obtaining the one-tenth ore deposit probability of workspace, engineering verification is carried out in the work that will carry out exactly on the spot to the field.Facts have proved that use the prediction that the method for the invention is carried out, real result is effective.

Claims (5)

1, a kind of method for finding mine with cubic forecast model based on three-dimensional modeling is characterized in that concrete steps are as follows:

(1) Data acquisition;

(2) Database Construction;

(3) foundation of three-dimensional mine model;

(4) foundation of cubic forecast model;

(5) checking that predicts the outcome;

The concrete steps of setting up of described cubic forecast model are:

(1) extraction of geologic body three-dimensional cube: 1., set up empty scope nugget model; 2., load piece; 3., give the cell block ascription; Required constituent content to sample in the boring is analyzed, and the weighting of service range inverse ratio is carried out required constituent content interpolation to the unknown elements piece; The three-dimensional unusual figure of analytical element, and the anomaly threshold of required element is carried out value according to the unusual figure of three-dimensional;

(2) calculating of geological prospecting variable and extraction: stratum, fracture, rock mass and four typess of variables of geochemical exploration, respectively these four typess of variables are carried out the division and the extraction of cubic units;

(3) the collective geological prospecting information of calculating each geologic unit and being comprised: the software platform of application is the Access database software in Microsoft Office 2003 software packages, will be imported among the Access by the autonomous integrated database of MicromineKAN TAN 3D; Step is: 1., on research main ore-controlling factors and favourable indicator for deposit basis, choose geology, the geochemical exploration Information sign of statistical study variable, and each sign of agreement to have value in the unit be 1, do not exist and get 0, each sign of statistics is in the distribution of each unit; In the computation process, the high size unification of nugget length and width is 10m * 10m * 10m, calculates then; 2., look for the ore deposit quantity of information to calculate: to look for the ore deposit quantity of information to calculate with conditional probability; Utilize each indicator for deposit quantity of information to calculate the quantity of information summation of each elementary cell respectively, determine to look for ore deposit quantity of information critical value according to the subjective probability method then;

(4) the ore deposit probability calculation is distinguished in prediction:

With a regional classification is some unit, has four kinds of situations:

1. blank;

2. the x unit represents to have certain geology sign to exist;

3. the xy unit represents to have this kind geology sign to have the ore deposit to exist again;

4. the y unit is represented to have the ore deposit to exist, but is not had this geology sign;

Corresponding frequency is respectively: f ₀, f _x, f _Xy, f _y

Appoint in the unit that has the geology sign to exist and get 1, it may contain the probability P in 1 mineral deposit _xFor:

$P_x=\frac{f_{\mathrm{xy}}}{(f_x+f_{\mathrm{xy}})}---\left(4\right)$

Appoint in the unit that does not have the geology sign to exist and get 1, it may contain the probability P in 1 mineral deposit ₀For:

$P_0=\frac{f_y}{(f_y+f_0)}---\left(5\right)$

Setting looks for ore deposit informational capacity lower limit as the comprehensive geology sign, utilizes probability estimate method Estimation Study district ore body unit output number to calculate f _x+ f _XyAnd f _Xy, and obtain P according to (4) formula _xCalculate f ₀+ f _vAnd f _v, and obtain P according to (5) formula ₀

Select the high control zone of working level, obtain correction factor, be used for revising P _xAnd P ₀, get P _x' and P ₀', be used for the prediction district with revised probability, obtain predicting the estimation of district's mineral deposit number; Concrete steps are as follows:

Selecting the study area in elevation interval is that the control zone is revised, unit, calculation control district sum and gross information content are positioned at the unit number on the lower limit, correction coefficient is=the control zone ore body unit number of the known ore body unit number in control zone/obtained by the output probability of aforementioned calculation;

Proofread and correct the probability P after obtaining proofreading and correct with correction coefficient _x' and P ₀'.

2, a kind of method for finding mine with cubic forecast model based on three-dimensional modeling according to claim 1 is characterized in that described Data acquisition, comprises: mining geology, geochemistry, mine engineering data.

3, a kind of method for finding mine with cubic forecast model based on three-dimensional modeling according to claim 1 is characterized in that described database establishment comprises: the construction in borehole data storehouse.

4, a kind of method for finding mine with cubic forecast model according to claim 1 based on three-dimensional modeling, it is characterized in that the foundation of described three-dimensional mine model comprises: the foundation of face of land DEM model, stratigraphic model, fracture model, granite model, boring model, element anomalies three-dimensional model.

5, a kind of method for finding mine with cubic forecast model based on three-dimensional modeling according to claim 1 is characterized in that, described predicting the outcome is verified as the field and verifies on the spot or engineering verification.

Images