Three-dimensional modeling method and device for ore body

Abstract

The application provides a method and device for three-dimensional modeling of ore bodies, and relates to the field of earth exploration and information technology.The application introduces volume change evaluation indexes and reaction change evaluation indexes, realizes fine simulation of dynamic behaviors of ore bodies under different temperature and pressure conditions, improves the precision and dynamic response capability of the three-dimensional model, and makes the three-dimensional model more truly reflect the actual situation of the ore bodies.The evaluation indexes of volume change and reaction change provide more scientific and comprehensive decision support for managers, optimize the production process of the mining industry, and improve the overall operation efficiency.

Description

Technical Field

[0001] This invention relates to the field of Earth exploration and information technology, specifically to a method and apparatus for three-dimensional modeling of ore bodies. Background Technology

[0002] Three-dimensional modeling technology has wide applications in geological exploration and ore body research. Early ore body 3D modeling relied primarily on surface mapping and a limited amount of borehole data. While this method reflected the distribution of the ore body to some extent, the accuracy and reliability of the model were significantly limited due to insufficient data. With the continuous advancement of geophysical exploration, remote sensing, and computer technologies, ore body 3D modeling has gradually shifted towards more refined and comprehensive methods, such as incorporating high-resolution seismic imaging, geomagnetic measurements, and gravity measurements. These new technologies can provide richer subsurface information, thereby significantly improving the accuracy of the 3D model.

[0003] In the prior art, publication number CN102279980A, entitled "A Method and Apparatus for Three-Dimensional Modeling of Geological Exploration Ore Bodies," the method includes: Step 1, delineating the ore body for a single exploration project based on the exploration project to obtain the upper and lower boundaries of the ore body; Step 2, projecting the exploration project onto an exploration profile and connecting the upper and lower boundaries of the ore body on the exploration profile to obtain the geological body boundary of the ore body on the exploration profile; Step 3, connecting the geological body boundaries between exploration profiles using a three-dimensional curved surface to obtain the three-dimensional reconstructed geological body. This method can realize geological three-dimensional modeling from borehole to profile to curved surface to finally generate a solid entity.

[0004] Despite significant advancements in existing 3D modeling techniques, numerous challenges remain in handling the dynamic behavior of ore bodies under varying temperature and pressure. Current 3D models are primarily based on static data, lacking effective characterization of dynamic changes in mineral composition and ore body volume, and thus failing to accurately predict ore body behavior during mining. This limitation leads to low ore body resource utilization, frequent unforeseen problems during mining, and negatively impacts the economic efficiency and safety of mining production.

[0005] The information disclosed in the background section above is only intended to enhance the understanding of the background of this disclosure, and therefore may include information that does not constitute prior art known to those skilled in the art. Summary of the Invention

[0006] The purpose of this invention is to provide a method and apparatus for three-dimensional modeling of ore bodies to solve the problems mentioned in the background art.

[0007] To achieve the above objectives, the present invention provides the following technical solution:

[0008] The 3D modeling method for ore bodies includes the following specific steps:

[0009] Step S1: Determine the ore body study area and collect geological exploration data for the area. The geological exploration data includes mineral composition, initial volume of the ore body, and geological mapping data. Based on this data, construct a three-dimensional model of the current ore body study area.

[0010] Step S2: Collect core samples of the ore body from Step S1 and measure the volume change data of the core samples under increasing temperature or pressure.

[0011] Step S3: Calculate the volume change rate of the core sample in Step S2 for each increase in temperature or pressure, and analyze these volume change rates to generate a comprehensive volume change evaluation index.

[0012] Step S4: Determine the mineral composition reaction rate of the core sample under increasing temperature or pressure conditions;

[0013] Step S5: Calculate the rate of change of reaction rate of mineral composition of core sample in step S4 after each increase of temperature or pressure, and analyze the reaction rate change rate to generate a comprehensive reaction change evaluation index.

[0014] Step S6: Combine and analyze the volume change evaluation index and the reaction change evaluation index to generate a fine-tuning strategy for adjusting the three-dimensional model of the ore body study area.

[0015] Furthermore, a fine-tuning strategy is comprehensively generated for adjusting the 3D model of the ore body study area, specifically including:

[0016] The output value of the fine-tuning strategy is denoted as FineTuning, and the calculation formula is as follows:

[0017]

[0018] Among them, I V I is the index for evaluating volume change. R It serves as an index for evaluating changes in response;

[0019] The output value of the fine-tuning strategy, FineTuning, takes a value in the range of (0,1), and this range is divided into the following intervals;

[0020] The interval (0, 0.25) indicates that the 3D model of the current ore body study area needs to be adjusted at level one.

[0021] The first-level adjustment specifically includes:

[0022] More than 50% of the mineral composition samples, initial volume data of the ore body, and geological mapping data were re-acquired and analyzed; the grid resolution, number of sampling points, and sampling point spacing in the 3D model were adjusted to their current values. and Adjustments;

[0023] The interval (0.25, 0.5) indicates that the 3D model of the current ore body study area needs to be adjusted to a secondary level.

[0024] The second-level adjustments specifically include:

[0025] 25%–50% of the mineral composition samples, initial volume data of the ore body, and geological mapping data were remeasured and analyzed;

[0026] The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. and Adjustments;

[0027] The interval (0.5, 0.75) indicates that the 3D model of the current ore body study area needs to be adjusted to level three.

[0028] The Level 3 adjustment specifically includes:

[0029] The mineral composition samples, initial volume data of the ore body, and geological mapping data of 10% to 25% were remeasured and analyzed.

[0030] The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. and Adjustments;

[0031] Interval 4 (0.75, 1) indicates that the 3D model of the current ore body study area needs to be adjusted to level 4.

[0032] The Level 4 adjustment specifically includes:

[0033] The mineral composition samples, initial volume data of the ore body, and geological mapping data of 0% to 10% were remeasured and analyzed;

[0034] The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. and Adjustments.

[0035] A three-dimensional modeling device for ore bodies, the device being used to execute the aforementioned three-dimensional modeling method for ore bodies, comprising:

[0036] 3D model building module: used to determine the ore body study area and collect geological exploration data for the area, including mineral composition, initial volume of the ore body and geological mapping data, and based on this data, to build a 3D model of the current ore body study area;

[0037] The first measurement module is used to collect core samples of the ore body and measure the volume change data of the core sample under increasing temperature or pressure.

[0038] Volume change evaluation index generation module: used to calculate the volume change rate of core samples after each increase in temperature or pressure, and to analyze these volume change rates to generate a comprehensive volume change evaluation index.

[0039] The second measurement module is used to determine the reaction rate of mineral components in core samples under increasing temperature or pressure conditions.

[0040] The reaction change evaluation index generation module is used to calculate the rate of change of reaction rate of mineral composition of core samples with each increase of temperature or pressure, and to analyze these rate of change of reaction rate to generate a comprehensive reaction change evaluation index.

[0041] Fine-tuning module: This module combines the volume change evaluation index and the reaction change evaluation index to generate a fine-tuning strategy for adjusting the three-dimensional model of the ore body study area.

[0042] Compared with existing technologies, the beneficial effects of this invention are: by introducing volume change evaluation index and reaction change evaluation index, it realizes a fine simulation of the dynamic behavior of ore bodies under different temperature and pressure conditions; it improves the accuracy and dynamic response capability of the three-dimensional model, enabling it to more realistically reflect the actual situation of the ore body; and by integrating the evaluation indices of volume change and reaction change, it provides managers with more scientific and comprehensive decision support, optimizes the mining production process, and improves overall operational efficiency. Attached Figure Description

[0043] Figure 1 This is a schematic diagram of the overall method flow of the present invention;

[0044] Figure 2 This is a block diagram of the device of the present invention. Detailed Implementation

[0045] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to specific embodiments.

[0046] It should be noted that, unless otherwise defined, the technical or scientific terms used in this invention should have the ordinary meaning understood by one of ordinary skill in the art to which this invention pertains. The terms "first," "second," and similar terms used in this invention do not indicate any order, quantity, or importance, but are merely used to distinguish different components. Terms such as "comprising" or "including" mean that the element or object preceding the word encompasses the elements or objects listed following the word and their equivalents, without excluding other elements or objects. Terms such as "connected" or "linked" are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. Terms such as "upper," "lower," "left," and "right" are used only to indicate relative positional relationships; when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0047] Example 1:

[0048] Please see Figure 1 The present invention provides a technical solution:

[0049] The 3D modeling method for ore bodies includes the following specific steps:

[0050] Step S1: Determine the ore body study area and collect geological exploration data for the area. The geological exploration data includes mineral composition, initial volume of the ore body, and geological mapping data. Based on this data, construct a three-dimensional model of the current ore body study area.

[0051] Step S2: Collect representative core samples of the ore body from Step S1, and measure the volume change data of the core samples under increasing temperature or pressure.

[0052] Step S3: Calculate the volume change rate of the core sample in Step S2 for each increase in temperature or pressure, and analyze these volume change rates to generate a comprehensive volume change evaluation index.

[0053] Step S4: Determine the mineral composition reaction rate of the core sample under increasing temperature or pressure conditions;

[0054] Step S5: Calculate the rate of change of reaction rate of mineral composition of core sample in step S4 after each increase of temperature or pressure, and analyze the reaction rate change rate to generate a comprehensive reaction change evaluation index.

[0055] Step S6: Combine and analyze the volume change evaluation index and the reaction change evaluation index to generate a fine-tuning strategy for adjusting the three-dimensional model of the ore body study area.

[0056] To further explain, the study area of ​​the ore body was identified, and geological exploration data for that area was collected. This data included mineral composition, initial volume of the ore body, and geological mapping data. Based on this data, a three-dimensional model of the current ore body study area was constructed, specifically including:

[0057] Representative core samples of the ore body in step S1 were determined by an expert group.

[0058] The boundaries of the ore body study area were determined using a Geographic Information System (GIS); existing mining area maps and satellite imagery data were used to clarify the geographic coordinates and boundary range of the ore body study area.

[0059] After determining the study area, X-ray fluorescence spectrometry (XRF) and electron probe microanalysis (EPMA) were used to analyze core samples within the study area to obtain mineral composition data; the specific operating steps are as follows:

[0060] Drill holes were laid out within the study area to obtain core samples;

[0061] The collected core samples were prepared into thin sections for XRF and EPMA analysis.

[0062] Record the mineral composition data for each borehole location and mark it in a three-dimensional coordinate system; the mineral composition data includes the mineral composition ratio, mineral density, and mineral hardness;

[0063] The initial volume of the ore body was measured using a laser rangefinder and a ground laser scanner; the specific operating steps are as follows:

[0064] Scanning points were deployed within the ore body study area to conduct comprehensive laser scanning and establish initial three-dimensional point cloud data of the ore body.

[0065] Data processing: Using point cloud processing software, the initial point cloud data is processed into a three-dimensional model, and the initial volume of the ore body is calculated;

[0066] Geological mapping data of the ore body study area was obtained using geological mapping technology. This data includes stratum thickness, fault distribution, geological structure, and stratigraphic unit boundaries. The specific operational steps are as follows:

[0067] 1. Geological mapping: Conduct detailed geological mapping within the study area to record geological structural features such as faults and folds;

[0068] 2. Data processing: Integrate geological mapping data with mineral composition data and initial volume data to establish a comprehensive geological model;

[0069] Mineral composition data, initial volume data, and geological mapping data are imported into 3D modeling software. In this embodiment, Surpac, Leapfrog, or MICROMINE are selected for 3D modeling software to perform data fusion.

[0070] The 3D model is constructed using interpolation algorithms in 3D modeling software, including Inverse Distance Weighting (IDW) and Kriging. The 3D model includes model parameters such as mesh resolution, number of sampling points, and sampling point spacing.

[0071] The accuracy and reliability of the 3D model were verified through cross-validation and review by geological experts.

[0072] Export the final 3D model as a visualization file format. Visualization file formats include 3DPDF, OBJ, or STL.

[0073] To further explain, the increments of temperature and pressure are both set to n, forming an increasing sequence {1,2,…,i,…,n}, where i represents the i-th increment; and the initial values ​​of temperature and pressure are set to T0 and P0, respectively; and the increments of temperature and pressure for each increment are set to q1 and q2, respectively.

[0074] The initial volume of the core sample was measured using a 3D scanner as V0.

[0075] Based on the pressure value P0, the i-th incrementing temperature value is denoted as T. i And T i =T0+q1×i; The core sample is placed at temperature T i The volume below is denoted as

[0076] Based on the temperature value T0, the i-th increment of pressure value is labeled as P. i And P i =P0+q2×i; The core sample is subjected to pressure value P i The volume below is denoted as

[0077] The construction of the volume change evaluation index specifically includes:

[0078] It uses a high-temperature and high-pressure reactor, which has precise temperature and pressure control functions;

[0079] The rate of volume change of the core sample after each increase in temperature or pressure is expressed by the following formula:

[0080]

[0081] Where X∈{T,P}; Let be the rate of change of volume, and This represents the rate of volume change of the core sample after the i-th temperature increase. This represents the rate of volume change of the core sample after the i-th pressure increase.

[0082] Divide the increasing sequence {1,2,…,i,…,n} into the first sequence {1,2,…,k} and the second sequence {k+1,k+2,…,n}.

[0083] Calculate the average rate of volume change of the first segment of the sequence:

[0084]

[0085] in, Let k be the average volume change rate of the first segment of the sequence, and k be the dividing point, which is the median value of the current increasing sequence.

[0086] Calculate the average rate of volume change of the second segment of the sequence:

[0087]

[0088] in, The average volume change rate of the second sequence segment;

[0089] Define the volume change evaluation index as I V The calculation formula is as follows:

[0090]

[0091] Where η1 is the adjustment constant; 0.15≤η1≤0.68; the volume change evaluation index I is set. V The value of is in the range of (0,1);

[0092] When I V The closer it gets to 0, the lower the rate of volume change of the core sample in the first half of the temperature or pressure increase is compared to the rate of volume change in the second half. This means that the volume change of the core sample is relatively slow under lower temperature or pressure conditions, while the volume change becomes more significant as the temperature or pressure increases.

[0093] When I V The closer the value is to 1, the higher the rate of volume change of the core sample in the first half of the temperature or pressure increase phase compared to the second half. This indicates that the volume change of the core sample is more dramatic under lower temperature or pressure conditions, while the volume change tends to level off as temperature or pressure increases.

[0094] To further explain, the reaction rate change is analyzed by combining the changes in reaction rate after each increase in temperature and pressure, and a comprehensive reaction change evaluation index is generated; specifically including:

[0095] Determine the reaction rate of mineral components in core samples under increasing temperature or pressure conditions;

[0096] Mineral components include, but are not limited to: quartz Potassium feldspar Sodium feldspar muscovite calcite Diopside

[0097] The mineral components are divided into categories to form a classification sequence {1,2,…,j,…,m}, where j represents the index of the j-th category of mineral component and m represents the total number of mineral component categories;

[0098] Initial mineral composition analysis of core samples was performed using X-ray diffraction (XRD) and scanning electron microscopy (SEM) to record the initial mineral composition;

[0099] The formula for calculating the average reaction rate of mineral components after the i-th increase in temperature or pressure is as follows:

[0100]

[0101] Among them, C j (T i ) represents the temperature value T in the i-th increment. i In the equation, the mass fraction of the j-th mineral component; C j (P i ) represents the pressure value P in the i-th increment. i In the equation, the mass fraction of the j-th mineral component, R(T) i R(P) represents the average reaction rate of the mineral composition of the core sample after the i-th temperature increase. i ) represents the average reaction rate of the mineral composition of the core sample after the i-th pressure increase;

[0102] Formula for the rate of change of reaction rate:

[0103]

[0104] Among them, QV1 i QV2 represents the rate of change in the reaction rate of the mineral composition of the core sample after the i-th temperature increase; i denoted as the rate of change of the reaction rate of the mineral composition of the core sample after the i-th pressure increase;

[0105] Obtain the first segment {1,2,...,i,...,n} and the second segment {k+1,k+2,...,n} of the increasing sequence {1,2,...,k}.

[0106] First stage average reaction rate change rate:

[0107]

[0108] Second stage average reaction rate change rate:

[0109]

[0110] Where k is the dividing point, and DV1 and DV2 are the average reaction rate changes of the first and second segments, respectively;

[0111] Define the reaction change evaluation index as I R The calculation formula is as follows;

[0112]

[0113] Where η2 is the adjustment constant; 0.01≤η2≤0.77; set I R The value of is in the range of (0,1);

[0114] I R The closer it gets to 0, the smaller the change in reaction rate of the mineral components in the first half of the increasing temperature or pressure, and the larger the change in reaction rate in the second half; this means that the mineral reaction is slower under lower temperature or pressure conditions, but the reaction rate is faster as the temperature or pressure increases.

[0115] I R The closer the value is to 1, the greater the change in reaction rate of the mineral components in the first half of the increasing temperature or pressure, and the smaller the change in reaction rate in the second half. This indicates that the mineral reaction is more vigorous under lower temperature or pressure conditions, but the reaction rate tends to level off as temperature or pressure increases.

[0116] To further explain, the comprehensive fine-tuning strategy for adjusting the 3D model of the ore body study area includes:

[0117] The output value of the fine-tuning strategy is denoted as FineTuning, and the calculation formula is as follows:

[0118]

[0119] Among them, I V I is the index for evaluating volume change. R It serves as an index for evaluating changes in response;

[0120] The output value of the fine-tuning strategy, FineTuning, takes a value in the range of (0,1), and this range is divided into the following intervals;

[0121] The interval (0, 0.25) indicates that the 3D model of the current ore body study area needs to be adjusted at level one.

[0122] The first-level adjustment specifically includes:

[0123] More than 50% of the mineral composition samples, initial volume data of ore bodies, and geological mapping data were re-acquired and analyzed;

[0124] The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated.

[0125] The interval (0.25, 0.5) indicates that the 3D model of the current ore body study area needs to be adjusted to a secondary level.

[0126] The second-level adjustments specifically include:

[0127] 25%–50% of the mineral composition samples, initial volume data of the ore body, and geological mapping data were remeasured and analyzed;

[0128] The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. and Adjustments;

[0129] The interval (0.5, 0.75) indicates that the 3D model of the current ore body study area needs to be adjusted to level three.

[0130] The Level 3 adjustment specifically includes:

[0131] Mineral composition samples, initial ore body volume data, and geological mapping data of 10%–25% were remeasured and analyzed; the grid resolution, number of sampling points, and sampling point spacing in the 3D model were adjusted to their current values. and Adjustments;

[0132] Interval 4 (0.75, 1) indicates that the 3D model of the current ore body study area needs to be adjusted to level 4.

[0133] The Level 4 adjustment specifically includes:

[0134] The mineral composition samples, initial volume data of the ore body, and geological mapping data of 0% to 10% were remeasured and analyzed;

[0135] The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. and Adjustments.

[0136] When I V When the value is increased, the output value of the fine-tuning strategy increases, and the corresponding model needs to be adjusted less.

[0137] A 10% increase in the rate of volume change results in a 5% increase in the FineTuning value;

[0138] The rate of change of the reaction rate remains constant or the change is less than 5%;

[0139] When I R When the value is increased, the output value of the fine-tuning strategy increases, and the corresponding model needs to be adjusted less.

[0140] A 10% increase in the rate of change of the reaction rate results in a 5% increase in the FineTuning value;

[0141] The rate of volume change remains constant or the change is less than 5%.

[0142] Improving model accuracy and dynamic response: By introducing volume change evaluation indices and reaction change evaluation indices, the dynamic behavior of ore bodies under different temperature and pressure conditions can be simulated more accurately. This improves the accuracy of the 3D model, enabling it to more realistically reflect the actual ore body conditions.

[0143] Enhancing production safety: Dynamic models can predict potential hazards during mining, such as ore body instability and collapse. By providing early warnings and adjusting mining strategies, production risks can be effectively reduced, ensuring worker safety.

[0144] Reduce unforeseen problems: Existing static models may encounter unforeseen problems during the mining process, while dynamic models can more comprehensively consider the performance of the ore body under different environmental conditions, thereby reducing the frequency of these problems and improving the continuity and stability of production.

[0145] Scientific Decision Support: By integrating volume change evaluation indices and response change evaluation indices, managers can obtain more scientific and comprehensive decision support. This helps optimize mining production processes and improve overall operational efficiency.

[0146] Example 2:

[0147] Please see Figure 2 A three-dimensional modeling device for ore bodies, the device being used to execute the aforementioned three-dimensional modeling method for ore bodies, comprising:

[0148] 3D model building module: used to determine the ore body study area and collect geological exploration data for the area, including mineral composition, initial volume of the ore body and geological mapping data, and based on this data, to build a 3D model of the current ore body study area;

[0149] The first measurement module is used to collect representative core samples of the ore body and measure the volume change data of the core samples under increasing temperature or pressure.

[0150] Volume change evaluation index generation module: used to calculate the volume change rate of core samples after each increase in temperature or pressure, and to analyze these volume change rates to generate a comprehensive volume change evaluation index.

[0151] The second measurement module is used to determine the reaction rate of mineral components in core samples under increasing temperature or pressure conditions.

[0152] The reaction change evaluation index generation module is used to calculate the rate of change of reaction rate of mineral composition of core samples with each increase of temperature or pressure, and to analyze these rate of change of reaction rate to generate a comprehensive reaction change evaluation index.

[0153] Fine-tuning module: This module combines the volume change evaluation index and the reaction change evaluation index to generate a fine-tuning strategy for adjusting the three-dimensional model of the ore body study area.

[0154] The above formulas are all dimensionless calculations. The formulas are derived from software simulations based on a large amount of collected data to obtain the most recent real-world results. The preset parameters in the formulas are set by those skilled in the art according to the actual situation.

[0155] The above embodiments can be implemented, in whole or in part, by software, hardware, firmware, or any other combination thereof. When implemented in software, the above embodiments can be implemented, in whole or in part, as a computer program product. Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution.

[0156] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment, depending on actual needs.

[0157] The above description is merely a specific embodiment of this application, but the scope of protection of this application is not limited thereto. Any changes or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application.

Claims

1. A three-dimensional modeling method for ore bodies, characterized in that, The specific steps include: Step S1: Determine the ore body study area and collect geological exploration data for the area. The geological exploration data includes mineral composition, initial volume of the ore body, and geological mapping data. Based on this data, construct a three-dimensional model of the current ore body study area. Step S2: Collect core samples of the ore body from Step S1 and measure the volume change data of the core samples under increasing temperature and pressure. Step S3: Calculate the volume change rate of the core sample in Step S2 for each increase in temperature and pressure, and analyze these volume change rates to generate a comprehensive volume change evaluation index. Step S4: Determine the mineral composition reaction rate of the core sample under increasing temperature and pressure conditions; Step S5: Calculate the rate of change of reaction rate of mineral composition of core sample in step S4 with each increase of temperature and pressure, and analyze the reaction rate change rate to generate a comprehensive reaction change evaluation index. Step S6: Combine and analyze the volume change evaluation index and the reaction change evaluation index to generate a fine-tuning strategy for adjusting the three-dimensional model of the ore body study area.

2. The three-dimensional modeling method for ore bodies according to claim 1, characterized in that: Constructing a three-dimensional model of the current ore body research area, specifically including: The boundaries of the ore body study area were determined using a geographic information system. After the study area was determined, X-ray fluorescence spectrometry and electron probe microanalysis equipment were used to analyze the core samples in the study area to obtain mineral composition data. Initial volume measurement of the ore body was performed using a laser rangefinder and a ground laser scanner. Geological mapping technology is used to obtain geological mapping data of the ore body study area. The geological mapping data includes the thickness of the stratum, fault distribution, geological structure and stratigraphic unit boundaries. Import mineral composition data, initial volume data, and geological mapping data into 3D modeling software; mineral composition data includes mineral composition ratio, mineral density, and mineral hardness. The 3D model is constructed using interpolation algorithms in 3D modeling software. The final 3D model is then exported as a visualization file, which includes model parameters such as mesh resolution, number of sampling points, and sampling point spacing.

3. The three-dimensional modeling method for ore bodies according to claim 2, characterized in that: Set the increment number of temperature and pressure to n, forming an increasing sequence {1,2,…,i,…,n}, where i represents the i-th increment; and set the initial values ​​of temperature and pressure to T0 and P0, respectively; and set the increment of temperature and pressure to q1 and q2, respectively. The initial volume of the core sample was measured using a 3D scanner. ; Based on the pressure value P0, the i-th incrementing temperature value is marked as... ,and ; core samples were subjected to temperature values The volume below is denoted as ; Based on the temperature value T0, the i-th increment of pressure value is marked as... ,and ; The core sample was subjected to pressure value The volume below is denoted as .

4. The three-dimensional modeling method for ore bodies according to claim 3, characterized in that: The construction of the volume change evaluation index specifically includes: The rate of volume change of the core sample after each increase in temperature and pressure is expressed by the following formula: Where, X∈{ , }; Let be the rate of change of volume, and This represents the rate of volume change of the core sample after the i-th temperature increase. This represents the rate of volume change of the core sample after the i-th pressure increase. Divide the increasing sequence {1,2,…,i,…,n} into the first sequence {1,2,…,k} and the second sequence {k+1,k+2,…,n}. Calculate the average rate of volume change of the first segment of the sequence: in, Let k be the average volume change rate of the first segment of the sequence, and k be the dividing point, which is the median value of the current increasing sequence. Calculate the average rate of volume change of the second segment of the sequence: in, The average volume change rate of the second sequence segment; Define the volume change evaluation index as: The calculation formula is as follows: in, To adjust the constant term; Set a volume change evaluation index The value of is in the range of (0,1); when The closer it gets to 0, the lower the rate of volume change of the core sample in the first half of the temperature and pressure increase is compared to the rate of volume change in the second half. when The closer it is to 1, the higher the rate of volume change of the core sample in the first half of the temperature and pressure increase is compared to the rate of volume change in the second half.

5. The three-dimensional modeling method for ore bodies according to claim 4, characterized in that: By analyzing the rate of change of reaction rate corresponding to each increase in temperature and pressure, a comprehensive evaluation index for reaction change is generated; specifically including: Determine the mineral composition reaction rate of core samples under increasing temperature and pressure conditions; The mineral components are divided into categories to form a classification sequence {1,2,…,j,…,m}, where j represents the index of the j-th category of mineral component and m represents the total number of mineral component categories; The formula for calculating the average reaction rate of mineral components after the i-th increase in temperature and pressure is as follows: in, This represents the temperature value during the i-th increment. In the text, the mass fraction of the j-th mineral component; This represents the pressure value at the i-th increment. In the text, the mass fraction of the j-th mineral component is... Let be the average reaction rate of the mineral composition of the core sample after the i-th temperature increase. denoted as the average reaction rate of the mineral composition of the core sample after the i-th pressure increase; Formula for the rate of change of reaction rate: in, denoted as the rate of change of the reaction rate of the mineral composition of the core sample after the i-th temperature increase; denoted as the rate of change of the reaction rate of the mineral composition of the core sample after the i-th pressure increase; Obtain the first segment {1,2,...,i,...,n} and the second segment {k+1,k+2,...,n} of the increasing sequence {1,2,...,k}. Calculate the rate of change of the average reaction rate in the first stage: Calculate the rate of change of the average reaction rate in the second stage: Where k is the dividing point, and DV1 and DV2 are the average reaction rate changes of the first and second segments, respectively; Define the reaction change evaluation index as: The calculation formula is as follows; in, To adjust the constant term; ;set up The value of is in the range of (0,1); The closer it gets to 0, the smaller the change in reaction rate of the mineral composition in the first half of the increasing temperature and pressure, and the greater the change in reaction rate in the second half. The closer it is to 1, the greater the change in reaction rate of the mineral components in the first half of the increasing temperature and pressure, and the smaller the change in reaction rate in the second half.

6. The three-dimensional modeling method for ore bodies according to claim 5, characterized in that: A comprehensive fine-tuning strategy is generated for adjusting the 3D model of the ore body study area, specifically including: The output value of the fine-tuning strategy is calibrated as The calculation formula is as follows: in, As an index for evaluating volume change, It serves as an index for evaluating changes in response; Fine-tuning strategy output value The value range of is (0,1), and this range is divided into the following intervals; The interval (0, 0.25) indicates that the 3D model of the current ore body study area needs to be adjusted at level one. The first-level adjustment specifically includes: More than 50% of the mineral composition samples, initial volume data of ore bodies, and geological mapping data were re-acquired and analyzed; The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. , and Adjustments; The interval (0.25, 0.5) indicates that the 3D model of the current ore body study area needs to be adjusted to a secondary level. The second-level adjustments specifically include: The mineral composition samples, initial volume data of the ore body, and geological mapping data of 25% to 50% were remeasured and analyzed. The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. , and Adjustments; The interval (0.5, 0.75) indicates that the 3D model of the current ore body study area needs to be adjusted to level three. The Level 3 adjustment specifically includes: The mineral composition samples, initial volume data of the ore body, and geological mapping data of 10% to 25% were remeasured and analyzed; The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. , and Adjustments; Interval 4 (0.75, 1) indicates that the 3D model of the current ore body study area needs to be adjusted to level 4. The Level 4 adjustment specifically includes: The mineral composition samples, initial volume data of the ore body, and geological mapping data of 0% to 10% were remeasured and analyzed; The current values ​​of mesh resolution, number of sampling points, and sampling point spacing in the 3D model are respectively calculated. , and Adjustments.

7. A three-dimensional modeling device for ore bodies, characterized in that: The apparatus is used to perform the ore body three-dimensional modeling method according to any one of claims 1-6, including: 3D model building module: used to determine the ore body study area and collect geological exploration data for the area, including mineral composition, initial volume of the ore body and geological mapping data, and based on this data, to build a 3D model of the current ore body study area; The first measurement module is used to collect core samples from the ore body and measure the volume change data of the core sample under increasing temperature and pressure. Volume change evaluation index generation module: used to calculate the volume change rate of core samples after each increase in temperature and pressure, and to analyze these volume change rates to generate a comprehensive volume change evaluation index. The second measurement module is used to determine the reaction rate of mineral components in core samples under increasing temperature and pressure conditions. The reaction change evaluation index generation module is used to calculate the rate of change of reaction rate of mineral composition of core samples with each increase of temperature and pressure, and to analyze these rate of change of reaction rate to generate a comprehensive reaction change evaluation index. Fine-tuning module: This module combines the volume change evaluation index and the reaction change evaluation index to generate a fine-tuning strategy for adjusting the three-dimensional model of the ore body study area.

Images