CN-122023695-A - Ore finding method of porphyry type tungsten-molybdenum ore based on comprehensive geological modeling

CN122023695ACN 122023695 ACN122023695 ACN 122023695ACN-122023695-A

Abstract

The application relates to the technical field of three-dimensional modeling and prospecting, in particular to a method for prospecting a porphyry tungsten-molybdenum ore based on comprehensive geologic modeling, which comprises the steps of calculating the original mineralisation favorability of each voxel in a three-dimensional geologic model of a new area to be predicted based on a known mineralization rule, and representing the possibility of the existence of the ore obtained by superposition of each voxel based on a single geologic factor; the gradient characteristics of the voxels are represented by constructing local change vectors, and the direction vectors reflecting the whole flow direction are generated by combining neighborhood weighted average, so that the spatial consistency of the two is analyzed, the original indexes are corrected, the original indexes are ensured to accord with geological rules, and mineralizable areas are calibrated in new areas to be predicted based on the optimized mineralization favorability. The method solves the problem that the existing three-dimensional modeling cannot be predicted quantitatively, and improves the prediction accuracy and the geological rationality by optimizing the ore formation favorability by combining the deposit rule and the space constraint.

Inventors

LIU XINWEI
XUE YUSHAN
YANG HAITAO

Assignees

西安西北有色地质研究院有限公司

Dates

Publication Date: 20260512
Application Date: 20260116

Claims (10)

1. The method for searching the porphyry tungsten-molybdenum ore based on comprehensive geological modeling is characterized by comprising the following steps of: constructing a three-dimensional geological model of a new area to be predicted and voxelizing the three-dimensional geological model; Calculating the original mineralization favor of each voxel in the three-dimensional geological model of the new area to be predicted based on the known mineralization rule, wherein the original mineralization favor is used for representing the possibility of the existence of the ore obtained by superposition of each voxel based on a single geological factor; The method comprises the steps of analyzing the original mineforming favorability difference between adjacent voxels of each voxel, constructing a local change vector of each voxel, representing the change rate and the maximum change direction of the original mineforming favorability of each voxel in a three-dimensional space, carrying out weighted average on unit vectors of all voxels in the neighborhood, which point to a central voxel, by utilizing the distance from each voxel to each voxel in the neighborhood and the original mineforming favorability, and generating a direction vector reflecting the overall attribute change flow direction in the neighborhood of each voxel; based on the optimized mineralisation benefits, the mineralisable zone is calibrated in the new zone to be predicted.
2. The method for prospecting a porphyry type tungsten molybdenum ore based on comprehensive geologic modeling as defined in claim 1, wherein constructing a three-dimensional geologic model of a new area to be predicted comprises: Collecting surface geological survey data, borehole logging data, geophysical logging data and geochemical analysis data of a new area to be predicted; all the acquired data are input into three-dimensional geological modeling software, and a three-dimensional geological model comprising a rock mass, a stratum, a fault, a contact zone and a mineral body is constructed through interpolation and geological interpretation.
3. The method for mining a porphyry tungsten molybdenum ore based on comprehensive geologic modeling as defined in claim 1, wherein the calculating of the original mining benefit of each voxel in the three-dimensional geologic model of the new area to be predicted based on the known mineralization rules comprises: Calculating the mineralization contribution weight of each geological attribute according to the numerical difference of each geological attribute between a mineralization area and a non-mineralization area in the known deposit model and the fluctuation degree of the numerical value of each geological attribute in the mineralization area; Carrying out normalized mapping on the numerical values of each geological attribute in the three-dimensional geological model to determine attribute favorability; and carrying out weighted summation on the attribute favorability of each geological attribute value in each voxel in the three-dimensional geological model and the corresponding mineralization favorability weight to obtain the original mineralization favorability of each voxel in the three-dimensional geological model.
4. The method for searching for ores of porphyry type tungsten molybdenum ore based on comprehensive geologic modeling as defined in claim 3, wherein the calculation process of the ore formation favorability weight is as follows: respectively calculating the numerical value average value of each geological attribute in all voxels in a mineralized area of each known deposit model and the numerical value average value of all voxels in a non-mineralized area, and calculating the absolute difference value between the two average values; calculating the numerical standard deviation of each geological attribute in all voxels in the mineralization area of each known deposit model; The ore-forming favorability weight respectively forms a positive correlation with the absolute difference value and forms a negative correlation with the numerical standard deviation.
5. A method for mining a porphyry tungsten molybdenum ore based on integrated geologic modeling as defined in claim 3, wherein said normalizing the mapping of the values of each geologic attribute in the three-dimensional geologic model to determine attribute benefit comprises: calculating the variation coefficient of the corresponding value of each geological attribute in all known deposit models; if the normalized value of the variation coefficient of the current geological attribute is smaller than a preset threshold, performing normalized mapping on each geological attribute value under the current geological attribute by using a Gaussian function to obtain attribute favorability of each geological attribute value; otherwise, performing normalized mapping on each geological attribute value under the current geological attribute by using a linear normalization method to obtain the attribute benefit of each geological attribute value.
6. The method for searching for the ore of the porphyry tungsten molybdenum ore based on comprehensive geologic modeling as defined in claim 1, wherein the construction process of the local change vector of each voxel is as follows: And calculating the finite differences of the original ore formation favorability among the upper and lower voxels, the left and right voxels and the front and rear voxels respectively by using a central difference method, and forming the local change vector of each voxel by using all the finite differences.
7. The method for searching for the ore of the porphyry type tungsten molybdenum ore based on comprehensive geological modeling as defined in claim 1, wherein the method for generating the direction vector is as follows: The method comprises the steps of respectively calculating a Gaussian function value of a distance from each voxel to each voxel in a neighborhood of each voxel and an original ore formation favorability difference between each voxel and each voxel in the neighborhood, taking a product of the Gaussian function value and the original ore formation favorability difference as a weight of a unit vector of each voxel, which points to each voxel, in the neighborhood of each voxel, so as to weight the unit vector to obtain a weighted vector, and averaging all the weighted vectors to obtain a direction vector.
8. The method for searching for ore of porphyry type tungsten molybdenum ore based on comprehensive geologic modeling as defined in claim 1, wherein optimizing the original ore formation availability comprises: and calculating an exponential function value taking a natural constant as a base and taking a normalized value of the spatial consistency as an independent variable, and taking the product of the original ore-forming favorability of each voxel and the corresponding exponential function value as the ore-forming favorability of each voxel after optimization.
9. The method for mining a porphyry tungsten molybdenum ore based on comprehensive geologic modeling as defined in claim 8, wherein the spatial consistency is an absolute value of projection of the local change vector of each voxel on the direction vector.
10. The method for mining a porphyry tungsten molybdenum ore based on comprehensive geologic modeling as defined in claim 1, wherein said calibrating the mineralizable zone in the new zone to be predicted comprises: Counting the cumulative probability distribution of the ore formation favorability of all voxels in the three-dimensional geological model after optimization, selecting a preset quantile as a threshold value, taking the voxels with the ore formation favorability greater than or equal to the threshold value after optimization as mineralized voxels, and taking the area formed by all connected mineralized voxels as a mineralizable area.

Description

Ore finding method of porphyry type tungsten-molybdenum ore based on comprehensive geological modeling Technical Field The application relates to the technical field of three-dimensional modeling and prospecting, in particular to a porphyry type tungsten-molybdenum ore prospecting method based on comprehensive geological modeling. Background With the development of three-dimensional geologic modeling technology, three-dimensional geologic modeling has become an important tool in the field of mineral exploration. Depending on Leapfrog, micromine and other professional software platforms, geologists can efficiently integrate multi-source heterogeneous data such as geology, geophysics, geochemistry and the like, and a high-precision digital three-dimensional geological model is constructed. However, in the prior art, three-dimensional modeling is mostly stopped at visual expression and geometric description of geological entities, namely, the three-dimensional modeling is mainly used as a static visual display tool of geological interpretation results, but is not a decision support system capable of actively carrying out quantitative analysis and intelligent prediction, and the lack of functions leads to insufficient excavation and effective utilization of deep mining information and space association rules contained in the model, and automatic circle determination of a target area cannot be realized directly based on the model, so that the accuracy and geological rationality of prospecting prediction are reduced. Disclosure of Invention In order to solve the technical problems, the application provides a method for searching for ores of porphyry type tungsten-molybdenum ores based on comprehensive geologic modeling, which aims to solve the existing problems. The application discloses a method for searching for ore of porphyry type tungsten-molybdenum ore based on comprehensive geological modeling, which adopts the following technical scheme: the application provides a method for searching for a porphyry tungsten-molybdenum ore based on comprehensive geological modeling, which comprises the following steps: constructing a three-dimensional geological model of a new area to be predicted and voxelizing the three-dimensional geological model; Calculating the original mineralization favor of each voxel in the three-dimensional geological model of the new area to be predicted based on the known mineralization rule, wherein the original mineralization favor is used for representing the possibility of the existence of the ore obtained by superposition of each voxel based on a single geological factor; The method comprises the steps of analyzing the original mineforming favorability difference between adjacent voxels of each voxel, constructing a local change vector of each voxel, representing the change rate and the maximum change direction of the original mineforming favorability of each voxel in a three-dimensional space, carrying out weighted average on unit vectors of all voxels in the neighborhood, which point to a central voxel, by utilizing the distance from each voxel to each voxel in the neighborhood and the original mineforming favorability, and generating a direction vector reflecting the overall attribute change flow direction in the neighborhood of each voxel; based on the optimized mineralisation benefits, the mineralisable zone is calibrated in the new zone to be predicted. Preferably, constructing a three-dimensional geological model of the new region to be predicted includes: Collecting surface geological survey data, borehole logging data, geophysical logging data and geochemical analysis data of a new area to be predicted; all the acquired data are input into three-dimensional geological modeling software, and a three-dimensional geological model comprising a rock mass, a stratum, a fault, a contact zone and a mineral body is constructed through interpolation and geological interpretation. Preferably, the calculating the original mineralisation benefit of each voxel in the three-dimensional geological model of the new region to be predicted based on the known mineralisation law includes: Calculating the mineralization contribution weight of each geological attribute according to the numerical difference of each geological attribute between a mineralization area and a non-mineralization area in the known deposit model and the fluctuation degree of the numerical value of each geological attribute in the mineralization area; Carrying out normalized mapping on the numerical values of each geological attribute in the three-dimensional geological model to determine attribute favorability; and carrying out weighted summation on the attribute favorability of each geological attribute value in each voxel in the three-dimensional geological model and the corresponding mineralization favorability weight to obtain the original mineralization favorability of each voxel in the three-dimensional geological model.