CN104573332B - A kind of mineral deposit volume determines method - Google Patents
A kind of mineral deposit volume determines method Download PDFInfo
- Publication number
- CN104573332B CN104573332B CN201410805000.9A CN201410805000A CN104573332B CN 104573332 B CN104573332 B CN 104573332B CN 201410805000 A CN201410805000 A CN 201410805000A CN 104573332 B CN104573332 B CN 104573332B
- Authority
- CN
- China
- Prior art keywords
- mathematical mesh
- volume
- mathematical
- physical grid
- mineral deposit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
Method is determined the invention discloses a kind of mineral deposit volume, wherein the determination method includes：Obtain all blocks that the mineral deposit includes；Physical grid is set up according to resulting block, wherein each of each of described physical grid and described block are corresponded and both overlapping margins；Mathematical mesh is set up, to cause the Mathematical mesh that the physical grid is completely covered；And calculate the volume of the Mathematical mesh to determine the volume of all blocks.It so, it is possible to realize the volume for quickly and accurately determining mineral deposit.
Description
Technical field
The present invention relates to fields such as geological prospecting, reserves calculating, mining designs, in particular it relates to which a kind of mineral deposit volume is true
Determine method.
Background technology
In geological prospecting and exploitation mineral, it is necessary to estimate border and its regularity of distribution in mineral deposit according to inspecting hole data,
And then mineral deposit volume is determined, but for the more complicated mineral deposit body of geological structure, be generally difficult to be accurately determined volume, thus nothing
Method is preferably designed and constructed.
At present, mineral deposit volume is generally determined by the following method：
Rhizoma Sparganii column method, this method is larger due to the fineness ratio of use, therefore there is larger error, especially for one
For the ore body of a little mineral ores, Rhizoma Sparganii column method is excessively coarse, and calculation error is larger；
Curved surface fitting method, this method belongs to the calculating in junior unit face, and its process is extremely complex.
The content of the invention
Method is determined it is an object of the invention to provide a kind of mineral deposit volume, the volume determines that method can be rapidly and precisely
Determine the volume in mineral deposit.
To achieve these goals, the present invention provides a kind of mineral deposit volume and determines method, and the determination method includes：Obtain
All blocks that the mineral deposit includes；Physical grid is set up according to resulting block, wherein every in the physical grid
Each of one and the block are corresponded and both overlapping margins；Mathematical mesh is set up, to cause the mathematics net
The physical grid is completely covered in lattice；And calculate the volume of the Mathematical mesh to determine the volume of all blocks.
Preferably, the Mathematical mesh includes Mathematical mesh unit, and the determination method also includes：It is determined that each block exists
The maximum and minimum value extended on three-dimensional；The maximum and minimum value extended based on all blocks on three-dimensional come
Determine maximum and minimum value that the mineral deposit extends on three-dimensional；Select the geomery of the Mathematical mesh unit；
The maximum and minimum value and the geomery of the Mathematical mesh unit extended based on the mineral deposit on three-dimensional is true
The quantity of the fixed Mathematical mesh unit, to determine the geometrical correspondence of the physical grid and the Mathematical mesh unit；
And determine that the Mathematical mesh unit is with the geometrical correspondence of the Mathematical mesh unit according to the physical grid
It is no to be completely covered by the physical grid, partly cover or do not cover completely.
Preferably, the determination method also includes：If the Mathematical mesh unit is completely covered by the physical grid
Or part is covered, then retains the Mathematical mesh unit；If the Mathematical mesh unit is not covered completely by the physical grid,
Then delete the Mathematical mesh unit；And set up the Mathematical mesh unit and the covering relation of physical grid.
Preferably, the Mathematical mesh includes Mathematical mesh unit, calculates the volume of the Mathematical mesh and includes：Judge every
Individual mathematics grid cell is completely covered by the physical grid overlapped with Block Boundary or part is covered；If the mathematics
Grid cell is completely covered by the physical grid overlapped with Block Boundary, then calculates the volume Ve of the Mathematical mesh；If
The physical grid part covering that the Mathematical mesh unit is overlapped with Block Boundary, then calculate the Mathematical mesh unit
With the volume Vp of physical grid intersection；Ve or Vp, and all Ve to calculating are calculated for each Mathematical mesh unit
With Vp summations, the volume of the Mathematical mesh of each block is directed to acquisition；And according to the mathematics net for each block
The volume of lattice obtains the volume of the Mathematical mesh for the mineral deposit.
Preferably for the Mathematical mesh unit covered by the physical grid part, calculated according to pure integration
The volume Vp of the Mathematical mesh unit and physical grid intersection.
Preferably, according to the geological structure in the mineral deposit and surface of stratum information, it is all that the acquisition mineral deposit includes
Block.
Preferably, the Mathematical mesh be shaped as it is following at least one：Tetrahedron, hexahedron, spheroid.
Preferably, the Mathematical mesh is shaped as identical or different, and the volume of the Mathematical mesh is identical or different.
Preferably, the border of the physical grid represent surface where mineral deposit, geological structure face, hand excavation face it is each
Discontinuous face is planted, and using joint, crack, tomography, the line of demarcation of material, drill hole line, hand excavation boundary line as each of representative
Plant discontinuous line.
Preferably, the border in the mineral deposit be represented as it is following at least one：Point, straight line, curve, plane, curved surface.
By above-mentioned technical proposal, set up physical grid and Mathematical mesh so that each of physical grid with it is described
Each of block is corresponded and both overlapping margins, and causes Mathematical mesh that physical grid is completely covered, and calculates number
The volume of grid is learned, and then the volume for quickly and accurately determining mineral deposit can be realized.
Other features and advantages of the present invention will be described in detail in subsequent embodiment part.
Brief description of the drawings
Accompanying drawing is, for providing a further understanding of the present invention, and to constitute a part for specification, with following tool
Body embodiment is used to explain the present invention together, but is not construed as limiting the invention.In the accompanying drawings：
Fig. 1 is the flow chart that mineral deposit volume provided by the present invention determines method；And
Fig. 2 is the flow chart for setting up physical grid and Mathematical mesh according to embodiment of the present invention.
Embodiment
The embodiment of the present invention is described in detail below in conjunction with accompanying drawing.It should be appreciated that this place is retouched
The embodiment stated is merely to illustrate and explain the present invention, and is not intended to limit the invention.
Fig. 1 is the flow chart that mineral deposit volume provided by the present invention determines method.As shown in figure 1, the ore deposit that the present invention is provided
Bed volume determines that method can include：At step 11, all blocks that the mineral deposit includes are obtained；At step 12, root
Physical grid is set up according to resulting block, wherein, each of each of described physical grid and described block one
One correspondence and both overlapping margins；At step 13, Mathematical mesh is set up, to cause the Mathematical mesh that the thing is completely covered
Manage grid；And at step 14, calculate the volume of the Mathematical mesh to determine the volume of all blocks.Pass through institute of the present invention
The method carried can quickly and accurately determine the volume in mineral deposit.
Wherein, above-mentioned physical grid represent intrinsic continuous, the discrete boundary line of nature, interface and some
Spatial point, for example, the various non-companies such as the surface of rock stratum, geological structure face (joint, fault plane etc.), hand excavation face where mineral deposit
Continuous face, it is each of representative that it, which can also represent joint, crack, tomography, the line of demarcation of material, drill hole line, hand excavation boundary line,
The discontinuous line of kind, the wherein line of demarcation of material, for example, the interface of coal seam and rock stratum, the mineral ore boundary line of metal mineral reserves, the boundary
Line can be it is following at least one：Point, straight line, plane, curve, curved surface.
Due to the complexity of architectonic complexity, and the block in composition mineral deposit, therefore, mineral deposit volume is accurately being calculated
There is very big difficulty in aspect.Therefore, applicant herein proposes the method for setting up Mathematical mesh, by the optional number of shape
Learn grid to cover to mineral deposit, and Mathematical mesh is covered more than existing physical grid, that is to say, that so that Mathematical mesh
Include physical grid completely, the boundary line of physical grid can be coincided with Mathematical mesh or physical grid side
Boundary line is fallen into Mathematical mesh.
The mistake of physical grid and Mathematical mesh is set up below with reference to Fig. 2 by specific embodiment to be described in detail
Journey, but the present invention is not restricted to this.
Fig. 2 is the flow chart for setting up physical grid and Mathematical mesh according to embodiment of the present invention.As shown in Fig. 2
Step 21 place, according to the geological structure in the mineral deposit and surface of stratum information, obtains all blocks that the mineral deposit includes, its
In, can be by the database of mineral deposit and its surrounding rock stratum at the structure such as inspecting hole data and geological structure data, and further obtain
Geological structure and bottom surface information are obtained, furthermore it is possible to all blocks that mineral deposit includes are obtained by block method for searching,
Wherein block method for searching can use existing any suitable method, will generally be cut in the range of mineral deposit by various discontinuous faces
The part to be formed is cut as block；At step 22, physical grid is set up so that physical grid and the natural boundary of the block
Coincide；At step 23, it is determined that each maximum and minimum value of the block on three-dimensional, and determine the scope in mineral deposit；
At step 24, based on identified each block, maximum, minimum value and the scope in mineral deposit set up number on three-dimensional
Learn grid so that the mineral deposit (physical grid) is completely covered in Mathematical mesh, and wherein Mathematical mesh can include Mathematical mesh list
Member；At step 25, judge whether each Mathematical mesh unit is completely covered by the physical grid, if Mathematical mesh unit
It is completely covered by the physical grid, then goes to step 26, retains the Mathematical mesh unit；If Mathematical mesh unit is not by thing
Reason grid is completely covered, then goes to step 27, judge whether the Mathematical mesh unit is covered by physical grid part, if mathematics
Grid cell is covered by physical grid part, then goes to step 26；If Mathematical mesh unit is not covered by physical grid part
(not covering completely), then delete the Mathematical mesh unit, and then can set up the inclusion relation of Mathematical mesh and physical grid (i.e.
Covering relation).
For at step 23, it is determined that maximum and minimum value that each block extends on three-dimensional, and determine base
The maximum and minimum value extended in all blocks on three-dimensional determines that the mineral deposit extends most on three-dimensional
Big value and minimum value (i.e. the scope in mineral deposit), its specific implementation method can be as described below.
Each block i MinXi, MaxXi, MinYi, MaxYi, MinZi are calculated first, and (i.e. block i is three by MaxZi
Tie up the maximum and minimum value on direction), MinX, the MaxX in the whole mineral deposit comprising all blocks can be so solved,
MinY, MaxY, MinZ, MaxZ (maximum and minimum value of the i.e. whole mineral deposit on three-dimensional)；Then, the mathematics is selected
The geomery (the Mathematical mesh unit that identical or different geomery can be selected according to actual needs) of grid cell, according to
The whole mineral deposit scope of above-mentioned calculating and the geomery of mathematics grid cell determine the quantity of Mathematical mesh unit, to determine
State the geometrical correspondence of physical grid and the Mathematical mesh unit；And according to the physical grid and the Mathematical mesh
The geometrical correspondence of unit is to determine whether the Mathematical mesh unit is completely covered by physical grid, part covering or complete
Do not cover entirely.
It is for instance possible to use the mode of simple decile is the Mathematical mesh sectionalization for mineral deposit, e.g., from cubical
Mathematical mesh unit Xa × Yb × Zc, wherein, Xa, Yb, Zc is the size on three directions of cuboid, then along three directions
On isodisperse can be tried to achieve using following formula,
Nx=(MaxX-MinX)/Xa
Ny=(MaxY-MinY)/Yb
Nz=(MaxZ-MinZ)/Zc
A shared N=Nx × Ny × Nz mathematics grid cell can be thus calculated, Mathematical mesh list is being determined
After the quantity of member, such as, according to the physical dimension of physical grid and position, each Mathematical mesh unit and thing can be calculated
The geometrical correspondence of grid is managed, such as the two intersecting point, line and face, and then can be according to physical grid and Mathematical mesh list
The geometrical correspondence of member is to determine whether Mathematical mesh unit is completely covered by physical grid, part covers or do not covered completely
Lid, so as to set up the inclusion relation (i.e. covering relation) of Mathematical mesh and physical grid, can equally calculate physics certainly
The inclusion relation of grid and mathematics grid cell.
The volume of determination Mathematical mesh will be described in detail how by embodiment below, but the present invention is not
Limited to this.
For each block in mineral deposit, judge what is overlapped for each Mathematical mesh unit of the block with Block Boundary
The physical grid is completely covered or part is covered；If the thing that the Mathematical mesh unit is overlapped with Block Boundary
Reason grid is completely covered, then calculates the volume Ve of the Mathematical mesh；If the Mathematical mesh unit is overlapped with Block Boundary
Physical grid part covering, then computational mathematics grid cell and the volume Vp of physical grid intersection；For each
Mathematical mesh unit calculates Ve or Vp, and all Ve and Vp calculated are summed, to obtain the mathematics net for each block
The volume of lattice；And the Mathematical mesh for being directed to the mineral deposit is obtained according to the volume of the Mathematical mesh for each block
Volume, to obtain the volume in mineral deposit.Wherein, can be according to simple for the Mathematical mesh unit covered by physical grid part
Property integration the pure of Numerical Manifold Method (for example, integration) come computational mathematics grid cell part and physical grid intersection
Volume Vp.
For the ease of calculating volume, the shape of above-mentioned Mathematical mesh unit can for it is following at least one：Tetrahedron, six faces
The regular shapes such as body, spheroid, and the shape of Mathematical mesh unit can be with identical or different, and its size can be identical or not
Together, this can be set according to actual mineral deposit border, joint, crack, tomography, boundary line etc..If the volume phase of Mathematical mesh unit
Together, then the quantity of Mathematical mesh unit being completely covered by physical grid can be determined, is multiplied by with the volume of Mathematical mesh unit
The quantity for the Mathematical mesh unit being completely covered by physical grid obtains all Mathematical meshs being completely covered by physical grid
The volume of unit.
Method of the present invention by setting up physical grid and Mathematical mesh double-grid, is tied according to mineral deposit border, tomography etc.
Structure face sets up the relation (be completely covered or part is covered) of physical grid and Mathematical mesh to be determined to represent mineral deposit volume
The volume of Mathematical mesh, simplifies the determination of mineral deposit volume, and can improve the determination precision of mineral deposit volume；And then can be
Follow-up large deformation numerical computations provide good initial data；Can also be steady in follow-up mining design, mining process
Qualitative analysis provides direct help；And, can be at any time according to dug up mine bed and bed of not digging up mine in follow-up mining process
The variation on border and obtained more accurate data boundary is actually measured, further can accurately describe geological structure, ore deposit
The position of bed and the regularity of distribution, for more reasonably mining design provides convenient from now on.
The preferred embodiment of the present invention is described in detail above in association with accompanying drawing, still, the present invention is not limited to above-mentioned reality
The detail in mode is applied, in the range of the technology design of the present invention, a variety of letters can be carried out to technical scheme
Monotropic type, these simple variants belong to protection scope of the present invention.
It is further to note that each particular technique feature described in above-mentioned embodiment, in not lance
In the case of shield, can be combined by any suitable means, in order to avoid unnecessary repetition, the present invention to it is various can
The combination of energy no longer separately illustrates.
In addition, various embodiments of the present invention can be combined randomly, as long as it is without prejudice to originally
The thought of invention, it should equally be considered as content disclosed in this invention.
Claims (9)
1. a kind of mineral deposit volume determines method, it is characterised in that the determination method includes：
Obtain all blocks that the mineral deposit includes；
Physical grid is set up according to resulting block, wherein each of described physical grid with it is each in the block
Person corresponds and both overlapping margins；
Mathematical mesh is set up, to cause the Mathematical mesh that the physical grid is completely covered；And
The volume of the Mathematical mesh is calculated to determine the volume of all blocks；
Wherein, the Mathematical mesh includes Mathematical mesh unit, and the determination method also includes：
It is determined that maximum and minimum value that each block extends on three-dimensional；
The maximum and minimum value extended based on all blocks on three-dimensional determines that the mineral deposit is prolonged on three-dimensional
The maximum and minimum value stretched；
Select the geomery of the Mathematical mesh unit；
The maximum and minimum value and the shape chi of the Mathematical mesh unit extended based on the mineral deposit on three-dimensional
The very little quantity for determining the Mathematical mesh unit, to determine physical grid pass corresponding with the geometry of the Mathematical mesh unit
System；And
Determine that the Mathematical mesh unit is with the geometrical correspondence of the Mathematical mesh unit according to the physical grid
It is no to be completely covered by the physical grid, partly cover or do not cover completely.
2. determination method according to claim 1, it is characterised in that the determination method also includes：
If the Mathematical mesh unit is completely covered by the physical grid or part is covered, retain the Mathematical mesh list
Member；
If the Mathematical mesh unit is not covered completely by the physical grid, the Mathematical mesh unit is deleted；And
Set up the Mathematical mesh unit and the covering relation of physical grid.
3. determination method according to claim 1, it is characterised in that the Mathematical mesh includes Mathematical mesh unit, meter
Calculating the volume of the Mathematical mesh includes：
Judge that each Mathematical mesh unit is completely covered by the physical grid overlapped with Block Boundary or part is covered；
If the Mathematical mesh unit is completely covered by the physical grid overlapped with Block Boundary, the mathematics net is calculated
The volume Ve of lattice；
If the physical grid part covering that the Mathematical mesh unit is overlapped with Block Boundary, computational mathematics grid
The volume Vp of unit and physical grid intersection；
Ve or Vp is calculated for each Mathematical mesh unit, and all Ve and Vp calculated are summed, to obtain for each
The volume of the Mathematical mesh of block；And
The volume of the Mathematical mesh for the mineral deposit is obtained according to the volume of the Mathematical mesh for each block.
4. determination method according to claim 3, it is characterised in that for the mathematics covered by the physical grid part
Grid cell, the Mathematical mesh unit and the volume Vp of physical grid intersection are calculated according to Simplex Integration.
5. determination method according to claim 1, it is characterised in that according to the geological structure and surface of stratum in the mineral deposit
Information, obtains all blocks that the mineral deposit includes.
6. the determination method according to any one of claim 1-4 claim, it is characterised in that the Mathematical mesh
Be shaped as it is following at least one：Tetrahedron, hexahedron, spheroid.
7. the determination method according to any one of claim 1-4 claim, it is characterised in that the Mathematical mesh
It is shaped as identical or different, the volume of the Mathematical mesh is identical or different.
8. the determination method according to any one of claim 1-4 claim, it is characterised in that the physical grid
Border represents the surface, geological structure face, the various discontinuous faces in hand excavation face of rock stratum where mineral deposit, and with joint, split
Gap, tomography, the line of demarcation of material, drill hole line, the various discontinuous lines that hand excavation boundary line is representative.
9. determination method according to claim 8, it is characterised in that the border in the mineral deposit represent it is following at least one：
Point, straight line, curve, plane, curved surface.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410805000.9A CN104573332B (en) | 2014-12-19 | 2014-12-19 | A kind of mineral deposit volume determines method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410805000.9A CN104573332B (en) | 2014-12-19 | 2014-12-19 | A kind of mineral deposit volume determines method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104573332A CN104573332A (en) | 2015-04-29 |
| CN104573332B true CN104573332B (en) | 2017-07-28 |
Family
ID=53089378
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410805000.9A Active CN104573332B (en) | 2014-12-19 | 2014-12-19 | A kind of mineral deposit volume determines method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104573332B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113077460B (en) * | 2021-04-28 | 2022-10-04 | 河北工业大学 | Method and device for determining deformation amount of test piece and storage medium |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101114276A (en) * | 2007-08-28 | 2008-01-30 | 中国地质大学（武汉） | Solid body mineral deposit three-dimensional visual reserves calculation system and computation method |
-
2014
- 2014-12-19 CN CN201410805000.9A patent/CN104573332B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101114276A (en) * | 2007-08-28 | 2008-01-30 | 中国地质大学（武汉） | Solid body mineral deposit three-dimensional visual reserves calculation system and computation method |
Non-Patent Citations (3)
| Title |
|---|
| 三维数值流形方法的理论研究;骆少明等;《应用数学和力学》;20050930;第26卷(第9期);1027-1032页 * |
| 岩体工程数值模拟新方法-数值流形法（NMM）;李学辉;《石家庄铁路职业技术学院学报》;20070930;第6卷(第3期);22-26页 * |
| 数值流形方法及其在岩石力学中的应用;李树忱等;《力学进展》;20041125;第34卷(第4期);446-454页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104573332A (en) | 2015-04-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103514630B (en) | Fault structure three-dimensional modeling method | |
| CN104076392A (en) | Microearthquake focus positioning combined inversion method based on grid search and Newton iteration | |
| Jaboyedoff et al. | Complex landslide behaviour and structural control: a three-dimensional conceptual model of Åknes rockslide, Norway | |
| CN105093352A (en) | Method for measuring and calculating rock mass fracture rate in field | |
| Xie et al. | Three-dimensional stability evaluation of landslides and a sliding process simulation using a new geographic information systems component | |
| CN104536054A (en) | Stripe-shaped iron mine assessment method and device | |
| Akiska et al. | Three-dimensional subsurface modeling of mineralization: a case study from the Handeresi (Çanakkale, NW Turkey) Pb-Zn-Cu deposit | |
| CN104570070B (en) | A kind of method and apparatus for setting up two-dimentional near surface geological model | |
| CN104573332B (en) | A kind of mineral deposit volume determines method | |
| CN105487106B (en) | A kind of benefit big gun method based on the illumination of Gaussian ray bundle target zone energy | |
| KR101348787B1 (en) | 3d magnetic inversion method based on the assemblage of 2d forward modeling and inversion as a technique to reflect the a-priori information | |
| Elkadi et al. | 3D-GSIS geotechnical modelling of tunnel intersection in soft ground: the Second Heinenoord Tunnel, Netherlands | |
| CN105279305B (en) | A kind of Three Dimensional Ground laser scanner technique survey station choosing method | |
| Cook et al. | Civil Engineering Optimisation Tool for the Study of CERN's Future Circular Colliders | |
| Rink et al. | Data visualisation and validation for hydrological models | |
| Yingbin et al. | An initial study of the risk of coal-roof water flow based on the “three maps two predictions” method | |
| Zahm | Improved Mechanical and Natural Fracture Model of the Lower Austin Chalk Using Outcrop, Core, and Well Logs, Central Texas | |
| Poláček et al. | Some applications of GIS ARC/INFO in the evaluation of urban radon risk | |
| Elorrieta | Geometry Model for the Follo Line Project | |
| CN208795577U (en) | A kind of Multifunctional outdoor test piece embedded device | |
| Exposito Marcano | Determination of an Optimal Arrangement of Blocks for Dimension Stone Mining | |
| Wang et al. | 3-D geological visual model and numerical model based on the Secondary Development Technology of CAD | |
| Kincal et al. | Large-scale land-suitability mapping in the GIS environment for the construction site of the University Olympic Village in Izmir (Turkey) | |
| Alionescu et al. | Mathematic interpolation methods-support for an efficient 3D modeling of landslides in the context of displacements monitoring | |
| Russell | Increasing geotechnical data confidence through the Integration of laser scanner face mapping data into the Sishen iron ore mine geotechnical database |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |