Mineral inclusions state analysis method
By employing layered grinding and mathematical model correction, the problems of lag and one-sidedness in the detection of ore embedment state in existing technologies have been solved. This enables three-dimensional continuous observation and quantitative characterization of the target mineral embedment state, thereby improving the accuracy of grinding process design and the reliability of data.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- CHANGCHUN GOLD RES INST
- Filing Date
- 2026-04-03
- Publication Date
- 2026-06-16
Smart Images

Figure CN121954594B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the fields of mineral processing and geological testing and analysis technology, specifically to a method for analyzing the mineral embedding state, and more specifically to a quantitative analysis method for the embedding state of complex polymetallic minerals. Background Technology
[0002] In the field of mineral processing and beneficiation technology, the embedding state of target minerals in ore is a core process mineralogical parameter. This parameter encompasses key information such as the grain size, occurrence morphology, embedding characteristics, and inter-mineral relationships of the target minerals. It directly determines the selection and design of beneficiation process flow, the control of grinding fineness, the optimization and screening of separation methods, and the accurate prediction of beneficiation recovery rate. It is a key factor affecting beneficiation technical indicators, resource recovery efficiency, and technical and economic benefits.
[0003] For the analysis and detection of the embedding state of target minerals in ores, existing technologies mainly employ two conventional methods, as follows: The first type is the artificial heavy sand method, or automated detection methods based on automatic analysis systems for process mineralogical parameters (such as MLA, AMICS, etc.). These methods require selecting representative ore samples from a single crushing operation or at different grinding finenesses. After sample preparation, observation, and statistical analysis, parameters such as the degree of mineral liberation, intergrowth content, and type within a specific particle size range are obtained. However, this type of method has inherent technical limitations: it can only characterize the instantaneous embedding state of the ore at a specific crushing and grinding stage, and cannot systematically and continuously reflect the dynamic evolution of the target mineral embedding state throughout the entire process from coarse crushing to fine grinding. Furthermore, this method is a post-processing detection method, and the results lag behind the process design requirements, making it difficult to achieve forward-looking prediction and precise guidance of key process parameters. The second type is particle-by-particle crushing analysis and microscopic in-situ tracking methods. These methods attempt to restore the original embedding characteristics of target minerals by breaking single-particle ore through stepwise crushing and microscopic observation. However, they suffer from problems such as cumbersome operation procedures, extremely low detection efficiency, and long testing cycles. Furthermore, they are limited by the number of test samples and the observation range, making it impossible to obtain sufficient statistically representative test data. This makes them difficult to adapt to the actual needs of industrial mineral processing design and severely limits their engineering applications.
[0004] In view of this, existing embedding state analysis techniques generally suffer from technical problems such as limited detection, low efficiency, and insufficient data representativeness. There is a lack of an analytical method that can efficiently, systematically, and quantitatively characterize the complete embedding characteristics of target minerals in their original state. Constrained by this technical bottleneck, current mineral processing technology design and parameter optimization largely rely on engineering experience or limited detection data, failing to achieve precise matching and refined control of the entire crushing, grinding, and separation process. This severely restricts the upgrading and development of mineral processing technology towards intelligence, efficiency, and low consumption. Summary of the Invention
[0005] In view of the technical problems existing in the background art, the present invention provides a method for analyzing the embedding state of target polymetallic complex minerals based on layered grinding and data correction. This method simulates the grinding process, acquires embedding state data at different grinding depths in layers, and corrects errors through a mathematical model, thereby constructing a refined embedding state analysis process for target minerals under macroscopic ore conditions.
[0006] This invention provides a method for analyzing the mineral embedding state, which includes the following steps:
[0007] S1, Regular Sample Preparation: Select representative ore samples, pre-process them to prepare regular samples; determine the initial observation surface as the 0th layer measurement surface;
[0008] S2, Reference Surface Establishment and Measurement: Mineralogical analysis is performed on the 0th layer measurement surface to obtain the exposed area, embedding characteristic parameters, and reference values of the target mineral content in the state of fine particle encapsulation.
[0009] S3, Layered grinding and measurement: The initial observation surface is ground in layers n times, with a grinding depth of Δh each time, to obtain the measurement surfaces from the 1st layer to the nth layer, where n≥2; mineralogical analysis is performed on each new measurement surface under the same conditions as the 0th layer measurement surface to obtain the exposed area and embedding characteristic parameters of the target mineral in each layer, until the target mineral particles tend to be encapsulated by gangue.
[0010] S4, Gangue Intercalation Correction Coefficient Calculation: Based on the measured content of target minerals in the fine-grained state in the actual grinding product, a gangue intercalation correction coefficient is established to correct the gangue intercalation data obtained by layer measurement.
[0011] S5, Calculation of embedding state data: Based on the exposed area and embedding characteristic parameters of each layer of the measurement surface, combined with the gangue crystallization correction coefficient, the embedding state data of the target mineral in three-dimensional space is obtained through weighted calculation.
[0012] As a further improvement of the present invention, the embedded feature parameter is the crystallization ratio of the target mineral to other minerals or pores, including the crystallization ratio of the target mineral to metal sulfides, metal oxides, gangue minerals and pores.
[0013] As a further improvement of the present invention, in step S5, the weighted calculation adopts a weighted average algorithm based on the proportion of exposed area of each layer of measurement surface, and the formula is: Ce=1 / (n+1)×Σ(Sn / S0×Cn.e); or Ce=T / (n+1)×Σ(Sn / S0×Cn.e);
[0014] Where Ce is the crystallization ratio of the target mineral to the e-th type of mineral phase or pores; e=1 represents crystallization with metal sulfides; e=2 represents crystallization with metal oxides; e=3 represents crystallization with gangue minerals; e=4 represents crystallization with pores; Sn is the total exposed area of the target mineral on the n-th measurement surface; S0 is the total exposed area of the target mineral on the 0-th measurement surface; Cn.e is the crystallization ratio of the target mineral to the e-th type of mineral phase or pores on the n-th measurement surface; T is the gangue crystallization correction coefficient, which is introduced when e represents gangue minerals, and T=1 in other cases.
[0015] As a further improvement of the present invention, when e=1 or 2, Ce=1 / (n+1)×Σ(Sn / S0×Cn.e);
[0016] When e=3, Ce=T / (n+1)×Σ(Sn / S0×Cn.e);
[0017] When e=4, C4=1-C1-C2-C3.
[0018] As a further improvement of the present invention, the target mineral in the fine-particle encapsulated state is a target mineral particle with a particle size of 0.010~0.005mm; the correction coefficient T is the ratio of the content P' of the target mineral in the particle size range encapsulated by gangue in the actual grinding product of the representative ore sample to the content P of the target mineral in the particle size range encapsulated by gangue in the 0th layer measurement surface; P is the benchmark value of the target mineral content in the fine-particle encapsulated state, and the calculation formula is as follows: T=P' / P.
[0019] As a further improvement of the present invention, the formula for calculating the crystallization ratio Cn.e of the target mineral and the e-th type of mineral phase or pores on the nth measurement surface is as follows:
[0020] ;
[0021] Where n is the layer number of the measurement surface; e is the mineral phase or pore type of the intergrowth, e=1 represents intergrowth with metal sulfides; e=2 represents intergrowth with metal oxides; e=3 represents intergrowth with gangue minerals; e=4 represents intergrowth with pores; m is the number of the target mineral particle, m=1, 2, 3..., representing the 1st, 2nd, 3rd... target mineral particles detected in the automated mineralogical analysis; Sm.n is the exposed area of the mth target mineral particle on the nth measurement surface, in μm²; Cm.en is the proportion of the side length of the mth target mineral particle in contact with the eth type of mineral phase or pore on the nth measurement surface.
[0022] As a further improvement of the present invention, in step S3, the grinding depth Δh is 1 / 6 to 1 / 2 of the average particle size of the target mineral.
[0023] As a further improvement of the present invention, in step S3, the determination condition for the tendency to be encased by gangue is: the exposed area of the target mineral particles on the current measuring surface is zero, and / or the equivalent particle size is not greater than 10% of the initial particle size.
[0024] As a further improvement of the present invention, the mineral embedding state analysis method further includes an identification step before step S2, specifically: setting an identifier on a regular sample, placing it in a direction perpendicular to the measurement surface, with the width and height plane parallel to the measurement surface to establish a spatial reference system; the identifier is a metal sheet, made of copper or copper alloy, with dimensions of 0.1~0.5mm in length, 0.1~0.5mm in width, and 1.0~5.0mm in height.
[0025] As a further improvement of the present invention, the target mineral is a metal sulfide, including one or more combinations of chalcopyrite, pyrite, galena, sphalerite, and molybdenite; the ore to be tested is a polymetallic complex ore, including copper-molybdenum polymetallic ore, copper-lead-zinc polymetallic ore, precious metal ore, or rare metal ore.
[0026] As a further improvement of the present invention, the representative ore sample is a representative large ore block or drill core with a minimum particle size of not less than 4.0 cm; the measuring surface is ground, polished and carbon sprayed before mineralogical analysis, with polishing using 0.5~5.0 μm abrasive and carbon spraying thickness of 10~20 nm.
[0027] Beneficial effects:
[0028] The method provided by this invention is standardized in operation, with strong data comparability. It specifies a complete standardized process from representative sample selection, marker positioning, precision layered grinding to automated mineralogical analysis, ensuring the consistency and repeatability of test results across different samples, batches, and operators. This provides reliable technical support for mineral processing research and industrial production. Its technical advantages include:
[0029] 1. Improved accuracy of embedding state characterization: By creating multi-layer (at least three layers) continuous observation in three-dimensional space and introducing a weighted average algorithm based on the actual grinding depth and mineral exposure area changes, the limitations of traditional two-dimensional single-surface analysis in representing three-dimensional distribution are overcome, and the actual embedding relationship of the target mineral in the ore is reflected more realistically.
[0030] 2. Corrected the statistical bias of target minerals in the state of fine-grained inclusions: For the first time, a correction coefficient T for gangue crystal linkage was proposed. By comparing the content of target minerals in the state of fine-grained inclusions in the layered grinding statistics with the measured content in the final grinding product, the problem of undercounting fine-grained gangue inclusions caused by the limitation of the observation surface resolution was effectively corrected, making the evaluation data of "crystallization with gangue" more reliable.
[0031] 3. Achieved quantitative and normalized expression of the embedding state: The complex mineral embedding relationship (with metal sulfides, metal oxides, gangue, and pore intergrowth) is quantified into four proportional parameters (C1, C2, C3, C4), and the sum is 1. This provides a quantitative data foundation that can be directly used for modeling in subsequent grinding process simulation, selectivity prediction, and metallurgical response analysis.
[0032] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, and in order to make the above and other objects, features and advantages of the present invention more apparent and understandable, specific embodiments of the present invention are described below. Attached Figure Description
[0033] To more clearly illustrate the technical solution of the present invention, the accompanying drawings used in the present invention will be briefly described below. Obviously, the drawings described below are merely some embodiments of the present invention, and those skilled in the art can obtain other drawings based on these drawings without any creative effort.
[0034] Figure 1 This is a flowchart illustrating the mineral embedding state analysis method provided in this embodiment of the invention. Detailed Implementation
[0035] The embodiments of the technical solution of the present invention will now be described in detail with reference to the accompanying drawings. These embodiments are merely illustrative of the technical solution of the present invention and are therefore intended to limit the scope of protection of the present invention.
[0036] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains; the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the invention; the terms “comprising” and “having”, and any variations thereof, in the specification, claims, and foregoing description of the invention, are intended to cover non-exclusive inclusion.
[0037] In the description of the embodiments of this invention, technical terms such as "first" and "second" are used only to distinguish different objects and should not be construed as indicating or implying relative importance or implicitly specifying the number, specific order, or primary and secondary relationship of the indicated technical features. In the description of the embodiments of this invention, "multiple" means two or more, unless otherwise explicitly defined.
[0038] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of the invention. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.
[0039] In the description of the embodiments of this invention, the term "and / or" is merely a description of the relationship between related objects, indicating that three relationships can exist. For example, A and / or B can represent: A existing alone, A and B existing simultaneously, and B existing alone. Additionally, the character " / " in this document generally indicates that the preceding and following related objects have an "or" relationship.
[0040] In the description of the embodiments of the present invention, the term "multiple" refers to two or more (including two), similarly, "multiple sets" refers to two or more (including two sets), and "multiple pieces" refers to two or more (including two pieces).
[0041] In the description of the embodiments of the present invention, the technical terms "center," "longitudinal," "lateral," "length," "width," "thickness," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," "clockwise," "counterclockwise," "axial," "radial," and "circumferential" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are only for the convenience of describing the embodiments of the present invention and simplifying the description, and are not intended to indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on the embodiments of the present invention.
[0042] In the description of the embodiments of the present invention, unless otherwise explicitly specified and limited, the technical terms such as "installation," "connection," "joining," and "fixing" should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral part; they can refer to a mechanical connection or an electrical connection; they can refer to a direct connection or an indirect connection through an intermediate medium; they can refer to the internal communication of two components or the interaction between two components. Those skilled in the art can understand the specific meaning of the above terms in the embodiments of the present invention according to the specific circumstances.
[0043] To address the shortcomings of existing embedding state analysis techniques, such as limited detection, low efficiency, insufficient data representativeness, and inability to systematically and quantitatively describe the embedding state of target minerals in their original ore state, this invention provides a mineral embedding state analysis method. It proposes a layered grinding-reconstruction method for three-dimensional embedding states, extending traditional two-dimensional cross-sectional analysis to continuous three-dimensional spatial observation. By performing layer-by-layer controlled grinding on the same sample, multi-layer parallel measurement surface sequence data is obtained. An area-weighted algorithm is used to reconstruct multi-layer two-dimensional information into a three-dimensional embedding structure, overcoming the limitations of instantaneous observation on a single cross-section and systematically characterizing the continuous evolution of mineral embedding states during grinding. A layered statistical-overall correction mechanism for target minerals in the fine-grained encapsulation state is established, effectively solving the statistical bias caused by insufficient observation surface resolution. By statistically analyzing the gangue encapsulation content of target minerals in the fine-grained encapsulation state layer by layer and comparing it with the measured values of actual grinding products, correction coefficients are constructed to mathematically correct the three-dimensional reconstruction results, achieving quantitative correction of the mineral embedding state from observed values to true values. A quaternary normalized quantitative expression system for embedding states is constructed, transforming complex mineral intergrowth relationships into standardized calculable parameters. By defining four types of crystallization states and satisfying normalization constraints, a quantitative data foundation is formed that can be directly used for grinding simulation, ore dressing prediction, and metallurgical response analysis, providing a unified quantitative basis for modeling mineral processing processes.
[0044] Please refer to Figure 1 As shown, the present invention provides a method for analyzing the mineral embedding state, comprising the following steps:
[0045] S1, Regular Sample Preparation: Select representative ore samples and preprocess them to prepare regular samples; determine the initial observation surface as the 0th layer measurement surface; details are as follows:
[0046] Select representative large samples of the ore to be tested or drill cores (minimum grain size greater than or equal to 4.0 cm), cut them (after cutting, the maximum grain size should be smaller than the inner diameter of the sample preparation mold), and embed them into regular samples (such as cylinders or cubes). Specifically, place the cut sample in the sample preparation mold, add epoxy resin and its corresponding curing agent at a volume ratio of 1.0:0.2~1.0, control the resin layer height at 1.0~1.5 cm, treat with ultrasonic vibration for 10.0~20.0 min to remove air bubbles, and then place in a constant temperature environment of 45~65℃ to accelerate curing. Select the end with more visible metal sulfides as the initial observation surface (0th layer measurement surface).
[0047] In some implementations, an identification step is also included, as detailed below:
[0048] The sample is cut along the direction perpendicular to the measurement surface. A metal sheet (preferably made of pure copper, with a length, width and height of 0.1~0.5mm×0.1~0.5mm×1.0~5.0mm) is placed tangentially towards the centroid of the measurement surface, with its width and height planes parallel to the measurement surface.
[0049] S2, Reference Surface Establishment and Measurement: Mineralogical analysis is performed on the 0th layer measurement surface to obtain the exposed area, embedding characteristic parameters, and reference values of the target mineral content in the state of fine-grained encapsulation; specifically as follows:
[0050] S21. The end face to be measured (initial observation surface (0th layer measurement surface)) is initially ground and polished (polishing with 1.0~5.0 micrometer abrasive) to make a flat surface without obvious scratches. Epoxy resin and its corresponding curing agent are added again, with a volume ratio of 1.0:0.2~1.0. The height of the resin layer is controlled at about 0.1~0.5cm. Ultrasonic vibration treatment is performed for 10.0~20min.
[0051] The second mixed surface is then polished (using 0.5~1.0 micrometer abrasive) to create a smooth initial measurement surface, and then carbon is sprayed at 10~20nm, denoted as Ai (where i represents a different measurement surface number).
[0052] S22, perform automatic mineralogical analysis on i initial measurement surfaces to obtain the target mineral exposure area Sm.n (n=0, 1, 2, representing the initial measurement surface (layer 0), the first and second measurement surfaces after grinding, respectively, m is 1, 2, 3... representing the number of target mineral particles, and a total of m target mineral particles are measured on i measurement surfaces), the crystallization characteristics with other minerals Cm.en (e is 1, 2, 3, 4 representing the ratio of the edge length of the m-th target mineral particle in contact with metal sulfides, metal oxides, gangue minerals and pores (the ratio of the edge length of this particle in contact with metal sulfides, metal oxides, gangue minerals and pores to the total edge length of this particle, the total edge length of this particle is the sum of the edge lengths in contact with metal sulfides, metal oxides, gangue minerals and pores)), and the grain size characteristic parameter Dm (equivalent area circle diameter), as the baseline data for the "layer 0" analysis.
[0053] Sn is the sum of the measured areas of the target mineral on a certain measuring surface, i.e.
[0054] S0 = ΣSm.0;
[0055] S1=ΣSm.1;
[0056] S2 = ΣSm.2;
[0057] The formula for calculating the crystallization ratio Cn.e of the target mineral and the e-th type of mineral phase or pores on the nth measurement surface is:
[0058] ;
[0059] Where n is the layer number of the measurement surface; e is the mineral phase or pore type of the intergrowth, e=1 represents intergrowth with metal sulfides; e=2 represents intergrowth with metal oxides; e=3 represents intergrowth with gangue minerals; e=4 represents intergrowth with pores; m is the number of the target mineral particle, m=1, 2, 3..., representing the 1st, 2nd, 3rd... target mineral particles detected in the automated mineralogical analysis; Sm.n is the exposed area of the mth target mineral particle on the nth measurement surface, in μm²; Cm.en is the proportion of the side length of the mth target mineral particle in contact with the eth type of mineral phase or pore on the nth measurement surface.
[0060] Where C0.e represents the embedding state of the target mineral on the initial measurement surface, 0 indicates the measurement surface number, and e indicates the embedding state, i.e.:
[0061] The initial measurement surface shows a crystallization ratio of target minerals to metal sulfides of C0.1:
[0062] C0.1=Σ(Sm.0×Cm.1.0) / ΣSm.0;
[0063] The initial measurement surface shows a crystallization ratio of target minerals to metal oxides of C0.2.
[0064] C0.2=Σ(Sm.0×Cm.2.0) / ΣSm.0;
[0065] The initial measurement surface shows a crystallization ratio of target mineral to gangue minerals of C0.3.
[0066] C0.3=Σ(Sm.0×Cm.3.0) / ΣSm.0;
[0067] The initial measurement surface shows a target mineral to fracture intergrowth ratio of C0.4:
[0068] C0.4=Σ(Sm.0×Cm.4.0) / ΣSm.0.
[0069] S23, measure the content of target minerals encapsulated in gangue in the state of 0.010~0.005mm fine particles, and record it as P;
[0070] Where, P=ΣSm 细 .0 / ΣSm.0(Sm 细 .0 represents the area of target mineral particles with a particle size of 0.010~0.005mm in the 0th layer;
[0071] S3, Layered Grinding and Measurement: The initial observation surface is subjected to n layers of grinding, with a grinding depth of Δh each time, to obtain the measurement surfaces from layer 1 to layer n, where n≥2; each new measurement surface is subjected to mineralogical analysis under the same conditions as the 0th layer measurement surface to obtain the exposed area and embedding characteristic parameters of the target mineral in each layer, until the target mineral particles tend to be encapsulated by gangue; specifically as follows:
[0072] The initial observation surface is subjected to deep grinding of a fixed thickness (Δh) using precision grinding or polishing equipment with controllable depth (preferably a grinding machine with micron-level feed accuracy).
[0073] The thickness Δh is preferably 1 / 6 to 1 / 2 of the average grain size of the target mineral.
[0074] After grinding, polishing, and carbon spraying, a new measurement surface is obtained. The same automated mineralogical analysis as in step S2 is performed on this new measurement surface, and the marking position is consistent with the measurement of layer 0 (i.e., the marking position remains unchanged when the sample is placed).
[0075] Obtain the following embedding state parameters of the target mineral in this layer (the nth layer): mineral exposure area Sm.n (m is 1, 2, 3... representing the number of target mineral particles, and a total of m chalcopyrite particles are measured on i measurement surfaces), and the intercrystallization characteristics with other minerals Cm.en (e is 1, 2, 3 representing the proportion of the edge length of the contact between the mth target mineral particle and metal sulfides, metal oxides, gangue minerals, and pores). Repeat this step, performing n=2 separate grinding and testing, testing each target mineral particle until it tends to be completely gangue intercrystallized on the observation surface, and the particle size is less than or equal to 0.1×Dm.
[0076] Gangue intergrowth state, meaning the target mineral is tightly wrapped by gangue or embedded between gangue grains or in fissures.
[0077] C1.e represents the embedding state of the target mineral on the first measurement surface, where 1 indicates the measurement surface number and e indicates the embedding state, i.e.:
[0078] The crystallization ratio of the target mineral and metal sulfide on the first measuring surface is C1.1:
[0079] C1.1=Σ(Sm.1×Cm.1.1) / ΣSm.1;
[0080] The ratio of the intercrystallization of the target mineral and metal oxide on the first measuring surface is C1.2:
[0081] C1.2=Σ(Sm.1×Cm.2.1) / ΣSm.1;
[0082] The crystallization ratio of the target mineral to gangue minerals on the first measuring surface is C1.3:
[0083] C1.3=Σ(Sm.1×Cm.3.1) / ΣSm.1;
[0084] The ratio of the target mineral to the intercrystallization of the fracture on the first measuring surface is C1.4:
[0085] C1.4=Σ(Sm.1×Cm.4.1) / ΣSm.1.
[0086] C2.e represents the embedding state of the target mineral on the second measurement surface, where 2 indicates the measurement surface number and e indicates the embedding state, i.e.:
[0087] The crystallization ratio of the target mineral and metal sulfide on the second measuring surface is C2.1:
[0088] C2.1=Σ(Sm.2×Cm.1.2) / ΣSm.2;
[0089] The ratio of the intercrystallization of the target mineral and metal oxide on the second measuring surface is C2.2:
[0090] C2.2=Σ(Sm.2×Cm.2.2) / ΣSm.2;
[0091] The crystallization ratio of the target mineral to gangue minerals on the second measuring surface is C2.3:
[0092] C2.3=Σ(Sm.2×Cm.3.2) / ΣSm.2;
[0093] The ratio of the target mineral to the fractured intergrowth on the second measuring surface is C2.4:
[0094] C2.4=Σ(Sm.2×Cm.4.2) / ΣSm.2.
[0095] S4, Gangue Intercalation Correction Coefficient Calculation: Based on the measured content of target minerals in the fine-grained encapsulated state in actual grinding products, a gangue intercalation correction coefficient is established to correct the gangue intercalation data obtained from layered measurements; details are as follows:
[0096] Instead of sampling, finely grind 45-55% of the particles to -0.074mm, take 2-3g of sample, prepare an automated mineralogical analysis sample, perform automated mineralogical analysis, measure the content of target minerals encapsulated in gangue in the state of 0.010-0.005mm fine particles, and record it as P'. Then the gangue correction number T = P' / P.
[0097] S5, Calculation of embedding state data: Based on the exposed area and embedding characteristic parameters of each layer of the measurement surface, combined with the gangue crystallization correction coefficient, the embedding state data of the target mineral in three-dimensional space is obtained through weighted calculation.
[0098] The weighted calculation adopts a weighted average algorithm based on the proportion of exposed area of each layer of measurement surface, and the formula is: Ce=1 / (n+1)×Σ(Sn / S0×Cn.e); or Ce=T / (n+1)×Σ(Sn / S0×Cn.e);
[0099] Where Ce is the crystallization ratio of the target mineral to the e-th type of mineral phase or pores; Sn is the total exposed area of the target mineral on the n-th measurement surface; S0 is the total exposed area of the target mineral on the 0-th measurement surface; Cn.e is the crystallization ratio of the target mineral to the e-th type of mineral phase or pores on the n-th measurement surface; T is the gangue crystallization correction coefficient, which is introduced when e represents gangue minerals, and T=1 in other cases.
[0100] When e=1, 2, Ce=1 / (n+1)×Σ(Sn / S0×Cn.e);
[0101] When e=3, Ce=T / (n+1)×Σ(Sn / S0×Cn.e);
[0102] When e=4, C4=1-C1-C2-C3.
[0103] For example, when n is 2, the specific calculation is as follows:
[0104] C1 represents the crystallization ratio of the target mineral to the metal sulfide;
[0105] C1=1 / 3×(C0.1+S1 / S0×C1.1+S2 / S0×C2.1);
[0106] C2 is the ratio of the target mineral to the metal oxide crystallization;
[0107] C2=1 / 3×(C0.2+S1 / S0×C1.2+S2 / S0×C2.2);
[0108] C3 represents the crystallization ratio of the target mineral to the gangue mineral;
[0109] C3=T / 3×(C0.3+S1 / S0×C1.3+S2 / S0×C2.3);
[0110] C4 represents the ratio of the target mineral to the intercrystalline structure of the pores;
[0111] C4 = 1 - C1 - C2 - C3.
[0112] Example 1
[0113] Embodiment 1 of this invention provides a method for analyzing the mineral embedding state, specifically a method for analyzing the embedding state of chalcopyrite in a copper-molybdenum polymetallic ore in Inner Mongolia. The copper-molybdenum polymetallic ore in Inner Mongolia contains pyrite as the main metal sulfide, followed by chalcopyrite, with smaller amounts of other minerals such as chalcocite, molybdenite, and lapis lazuli. Metal oxides are present in trace amounts. The gangue minerals are mainly quartz, followed by potassium feldspar, mica-like minerals, etc., and contain small amounts of plagioclase, illite, and kaolinite, etc. The ore is a low-sulfide copper ore.
[0114] Specifically, the copper minerals in this ore are mainly chalcopyrite, followed by a small amount of secondary copper minerals, with very few free and bound copper oxide minerals. Chalcopyrite grains in the ore are predominantly in the range of 0.100–0.010 mm, accounting for over 80%. The chalcopyrite is mainly found intercalated with gangue minerals, with other states of chalcopyrite being less common.
[0115] The mineral embedding state analysis method includes the following steps:
[0116] S1, Sample Preparation
[0117] Select representative large samples of copper-molybdenum polymetallic ore (minimum particle size greater than or equal to 2.0 cm), cut them (maximum particle size less than the inner diameter of the sample preparation mold 3.0 cm), and embed them into regular samples (this sample is a cylinder, add epoxy resin and its corresponding curing agent, volume ratio 1.0:0.5, control the resin layer height to about 1.0 cm, ultrasonically vibrate for 10.0 min to remove air bubbles, and then place it in a constant temperature environment of 45℃ to accelerate curing).
[0118] Specifically, a cylindrical sample with a diameter of 3.0cm × 1.0cm was prepared, and the end with more metal sulfides observed by the naked eye was used as the initial observation surface (layer 0).
[0119] Adding markings: The sample is cut along the direction perpendicular to the measurement surface, tangential to the center of the cylinder. The length, width, and height of the cut are greater than 0.1mm×0.1mm×1.0mm respectively. Place a metal sheet (made of copper, with a length, width, and height of 0.1mm×0.1mm×1.0mm respectively) with the width and height planes of the metal sheet parallel to the measurement surface.
[0120] S2, Establishment and Measurement of the Reference Surface
[0121] S21, this end face (initial observation surface (layer 0)) is initially ground and polished (polishing with 3.0 micron abrasive) to make a flat surface without obvious scratches. Epoxy resin and its corresponding curing agent are added again at a volume ratio of 1.0:0.5. The height of the resin layer is controlled at about 0.2 cm. Ultrasonic vibration treatment is performed for 10.0 min.
[0122] The second mixed surface was then ground and polished (ground flat, and finally polished with 0.5-micron abrasive) to create a smooth initial measurement surface, and then carbon was sprayed at 10nm, denoted as Ai (where i represents the different measurement surface number).
[0123] S22, perform automatic mineralogical analysis on the initial measurement surface Ai to obtain the target mineral exposure area Sm.0 (m is 1, 2, 3... representing the number of target mineral particles respectively, and a total of m chalcopyrite particles are measured on the i measurement surfaces), the crystallization characteristics with other minerals Cm.e.0 (e is 1, 2, 3, 4 representing the proportion of the edge length of the m-th chalcopyrite particle in contact with metal sulfides, metal oxides, gangue minerals and pores respectively), and the grain size characteristic parameter Dm, as detailed in Table 1.
[0124] Table 1. Automated Mineralogical Analysis Data of Initial Measurement Surface (n=0)
[0125]
[0126] Example calculation:
[0127] S0=ΣSm.0
[0128] =796.24+704.33+4942.87+436.74……+365.32+98.51+47.53
[0129] =1.659×10 7 μm 2 ;
[0130] C0.1 = Σ(Sm.0 × Cm.1.0) / ΣSm.0
[0131] =(796.24×23.56%+704.33×42.88%+436.74×7.61%+……+67.92×64.97%+365.32×100.00%) / 1.659×10 7 =1.513×10 6 / 1.659×10 7
[0132] =9.12%;
[0133] C0.2 = Σ(Sm.0 × Cm.2.0) / ΣSm.0
[0134] =(436.74×14.23%+42.18×16.49%+394.62×56.77%+……+47.53×4.25%) / 1.659×10 7
[0135] =1.309×106 / 1.659×10 7
[0136] =7.89%;
[0137] C0.3 = Σ(Sm.0 × Cm.3.0) / ΣSm.0
[0138] =(796.24×76.44%+704.33×57.12%+……+98.51×90.15%+47.53×80.97%) / 1.659×10 7
[0139] =1.160×10 7 / 1.659×10 7
[0140] =69.92%;
[0141] C0.4 = Σ(Sm.0 × Cm.4.0) / ΣSm.0
[0142] =(436.74×12.83%+42.18×38.29%+394.62×2.31%+……+98.51×9.85%) / 1.659×10 7
[0143] =2.168×10 6 / 1.659×10 7
[0144] =13.07%.
[0145] S23, the content of chalcopyrite encapsulated in gangue in the state of 0.010~0.005mm fine particles is measured and recorded as P;
[0146] P=ΣSm 细 .0 / ΣSm.0
[0147] =(42.18+70.95+……+42.97+67.92+47.53) / 1.659×10 7
[0148] =8.693×10 5 / 1.659×10 7
[0149] =5.24%.
[0150] S3, Layered Grinding and Measurement
[0151] The initial observation surface was subjected to deep grinding with a fixed thickness (Δh) using a grinding machine.
[0152] The preferred thickness Δh is 1 / 3 of the average grain size of the target mineral, Δh = 35.23 × 1 / 3 = 11.74 μm.
[0153] After grinding, polishing, and carbon spraying to 10 nm, a new measurement surface is obtained. The same automatic mineralogical analysis as in step S3 is performed on this new measurement surface (with the marker position unchanged) to obtain the following embedding state parameters of the target mineral within this layer (the n=1th layer): the mineral exposure area Sm.n (m represents the number of target mineral particles, 1, 2, 3… respectively, with a total of m chalcopyrite particles measured across i measurement surfaces), and the intercrystallization characteristics with other minerals Cm.en (e represents the proportion of the edge length of the mth chalcopyrite particle in contact with metal sulfides, metal oxides, gangue minerals, and pores, 1, 2, 3 respectively). This step is repeated for the n=2nd grinding and testing, with each target mineral particle tested until it tends to be completely intercrystallized with gangue on the observation surface, and the particle size is less than or equal to 0.1 × Dm.e.0, as detailed in Tables 2 and 3.
[0154] Specifically, when Sm.1 = 0 μm 2 At that time, Cm.1.1=0, Cm.2.1=0, Cm.3.1=100, Cm.4.1=0;
[0155] Specifically, when Sm.2 = 0 μm 2 At that time, Cm.1.2=0, Cm.2.2=0, Cm.3.2=100, Cm.4.2=0;
[0156] That is, at this time, the measuring surface is in a state where the tested particle is wrapped by gangue.
[0157] Table 2. Automated Mineralogical Analysis Data for Measurement Surfaces (n=1)
[0158]
[0159] S1=ΣSm.1
[0160] =705.24+731.52+4270.46+……+316.97+54.17
[0161] =1.503×10 7 μm 2 .
[0162] Table 3. Automated mineralogical analysis data for measurement surfaces (n=2)
[0163]
[0164] S2=ΣSm.2
[0165] =684.97+712.94+3820.74+……+224.94+26.94
[0166] =1.677×10 7 μm 2 .
[0167] S4, Data Correction Coefficient
[0168] Instead of using a sample for analysis, finely grind 0.074mm particles (50%) and take 2-3g of the sample to prepare an automated mineralogical analysis sample. Perform automated mineralogical analysis to measure the content of chalcopyrite encapsulated in gangue with 0.010-0.005mm fine particles, denoted as P'=5.35%.
[0169] The correction number for the vein is T = P' / P = 5.35% / 5.24% = 1.02.
[0170] S5, Embedded state data calculation
[0171] C1=1 / 3×(C0.1+S1 / S0×C1.1+S2 / S0×C2.1)
[0172] =1 / 3×(9.12%+1.503 / 1.659×10.36%+1.677 / 1.659×8.26%)
[0173] =8.95%;
[0174] C2=1 / 3×(C0.2+S1 / S0×C1.2+S2 / S0×C2.2)
[0175] =1 / 3×(7.89%+1.503 / 1.659×4.27%+1.677 / 1.659×6.33%)
[0176] =6.05%;
[0177] C3=T / 3×(C0.3+S1 / S0×C1.3+S2 / S0×C2.3)
[0178] =1.02 / 3×(69.92%+1.503 / 1.659×71.68%+1.677 / 1.659×73.10%)
[0179] =70.98%;
[0180] C4 = 1 - C1 - C2 - C3
[0181] =1 - 8.95% - 6.05% - 70.98%
[0182] =14.02%.
[0183] Table 4. Results of Embedded State Analysis
[0184]
[0185] The above data shows that in Example 1, the proportion of chalcopyrite linked to gangue minerals was the highest, reaching 70.98%, followed by linkage with pores, while linkage with metal sulfides and metal oxides was less common. Under normal grinding conditions, chalcopyrite tends to be linked with gangue minerals, so it is necessary to strengthen grinding to improve its selectivity.
[0186] Example 2
[0187] Embodiment 2 of this invention provides a method for analyzing the mineral embedding state, specifically a method for analyzing the embedding state of galena in a lead-zinc polymetallic ore from Inner Mongolia. The ore contains approximately 5.00% metal sulfides, mainly pyrrhotite and galena, with minor amounts of sphalerite, chalcopyrite, and arsenopyrite. The main metal oxide is magnetite, followed by hematite and limonite. The gangue minerals are mainly quartz and mica-like minerals, followed by feldspar minerals, with other gangue minerals present in smaller quantities. The ore's oxidation rate, calculated based on lead, is approximately 3.50%, classifying it as a primary ore. This ore is classified as a medium-sulfide lead-zinc polymetallic ore.
[0188] The mineral embedding state analysis method includes the following steps:
[0189] S1, Sample Preparation
[0190] A representative large sample block of the lead-zinc polymetallic ore (minimum particle size greater than or equal to 2.0 cm) was selected, cut (maximum particle size less than the inner diameter of the sample preparation mold 3.0 cm), and inlaid into a regular sample (this sample is a cylinder, epoxy resin and its corresponding curing agent are added, the volume ratio is 1.0:0.5, the resin layer height is controlled at about 1.0 cm, ultrasonic vibration is treated for 10.0 min to remove air bubbles, and then placed in a constant temperature environment of 45℃ to accelerate curing).
[0191] Specifically, a cylindrical sample with a diameter of 3.0cm × 1.0cm was prepared, and the end with more metal sulfides observed by the naked eye was used as the initial observation surface (layer 0).
[0192] Adding markings: The sample is cut along the direction perpendicular to the measurement surface, tangential to the center of the cylinder. The length, width, and height of the cut are greater than 0.1mm×0.1mm×1.0mm respectively. Place a metal sheet (made of copper, with a length, width, and height of 0.1mm×0.1mm×1.0mm respectively) with the width and height planes of the metal sheet parallel to the measurement surface.
[0193] S2, Establishment and Measurement of the Reference Surface
[0194] S21, this end face (initial observation surface (layer 0)) is initially ground and polished (polishing with 3.0 micron abrasive) to make a flat surface without obvious scratches. Epoxy resin and its corresponding curing agent are added again at a volume ratio of 1.0:0.5. The height of the resin layer is controlled at about 0.2 cm. Ultrasonic vibration treatment is performed for 15.0 min.
[0195] The second mixed surface was then ground and polished (ground flat, and finally polished with 0.5-micron abrasive) to create a smooth initial measurement surface, and then carbon was sprayed at 10nm, denoted as Ai (where i represents the different measurement surface number).
[0196] S22, perform automatic mineralogical analysis on the initial measurement surface Ai to obtain the target mineral exposure area Sm.0 (m is 1, 2, 3... representing the number of target mineral particles respectively, and a total of m galena particles are measured on the i measurement surfaces), the crystallization characteristics with other minerals Cm.e.0 (e is 1, 2, 3, 4 representing the proportion of the edge length of the m-th galena particle in contact with metal sulfides, metal oxides, gangue minerals and pores respectively), and the grain size characteristic parameter Dm, as detailed in Table 5.
[0197] Table 5. Automated mineralogical analysis data for the initial measurement surface (n=0)
[0198]
[0199] Example calculation:
[0200] S0=ΣSm.0
[0201] =572.56+135.32+1743.27+……+43.96+172.57+652.76+75.98
[0202] =5.270×10 8 μm 2 ;
[0203] C0.1 = Σ(Sm.0 × Cm.1.0) / ΣSm.0
[0204] =(7572.56×54.32%+7135.32×24.47%+……+652.76×100.00%) / 5.27×10 8
[0205] =8.291×10 7 / / 5.270×10 8
[0206] =15.73%.
[0207] S23, the content of galena encapsulated in gangue in the state of 0.010~0.005mm fine particles is measured and recorded as P;
[0208] P=ΣSm 细 .0 / ΣSm.0=2.44%.
[0209] S3, Layered Grinding and Measurement
[0210] The initial observation surface was subjected to deep grinding with a fixed thickness (Δh) using a grinding machine.
[0211] The preferred thickness Δh is 1 / 3 of the average grain size of the target mineral, Δh = 52.15 × 1 / 3 = 17.38 μm.
[0212] After grinding, polishing, and carbon spraying to 10 nm, a new measurement surface is obtained. The same automatic mineralogical analysis as in step S3 is performed on this new measurement surface (with the marker position unchanged) to obtain the following embedding state parameters of the target mineral within this layer (the n=1th layer): the mineral exposure area Sm.n (m represents the number of target mineral particles, 1, 2, 3… respectively, with a total of m galena particles measured across i measurement surfaces), and the intercrystallization characteristics with other minerals Cm.en (e represents the proportion of the edge length of the mth galena particle in contact with metal sulfides, metal oxides, gangue minerals, and pores, 1, 2, 3 respectively). This step is repeated for the n=2nd grinding and testing, with each target mineral particle tested until it tends to be completely intercrystallized with gangue on the observation surface, and the particle size is less than or equal to 0.1 × Dm.e.0, as detailed in Tables 6 and 7.
[0213] Table 6. Automated mineralogical analysis data for measurement surfaces (n=1)
[0214]
[0215] S1=ΣSm.1
[0216] =473.28+138.22+1647.51+……+164.88+526.73+82.16
[0217] =5.152×10 8 μm 2 .
[0218] Table 7 Automated mineralogical analysis data for measurement surfaces (n=2)
[0219]
[0220] S2=ΣSm.2
[0221] =373.28+138.22+1247.51+……+426.46+72.16
[0222] =4.832×10 8 μm 2 .
[0223] S4, Data Correction Coefficient
[0224] Instead of using a sample tester, finely grind 50% of the material to a thickness of -0.074 mm, take 2-3 g of the sample, prepare an automated mineralogical analysis sample, perform automated mineralogical analysis, and measure the content of 0.010-0.005 mm galena encapsulated by gangue, denoted as P'=2.63%;
[0225] The correction number for the vein is T = P' / P = 2.63% / 2.44% = 1.08.
[0226] S5, Embedded state data calculation
[0227] C1=1 / 3×(C0.1+S1 / S0×C1.1+S2 / S0×C2.1)
[0228] =1 / 3×(15.73%+5.152 / 5.270×12.74%+4.832 / 5.270×10.64%)
[0229] =12.65%;
[0230] C2=1 / 3×(C0.2+S1 / S0×C1.2+S2 / S0×C2.2)
[0231] =1 / 3×(8.92%+5.152 / 5.270×6.75%+4.832 / 5.270×4.87%)
[0232] =6.66%;
[0233] C3=T / 3×(C0.3+S1 / S0×C1.3+S2 / S0×C2.3)
[0234] =1.08 / 3×(70.53%+5.152 / 5.270×74.65%+4.832 / 5.270×70.57%)
[0235] =74.96%;
[0236] C4==5.73%.
[0237] Table 8. Results of Embedded State Analysis
[0238]
[0239] The above data shows that in Example 2, the proportion of galena linked with gangue minerals is the highest, reaching 74.96%, followed by the proportion linked with metal sulfides, with less linkage with metal oxides, and the least linkage porosity. Under normal grinding conditions, galena tends to be linked with gangue, so it is necessary to strengthen grinding to improve its selectivity.
[0240] In summary, this invention provides a method for analyzing the embedment state of target polymetallic complex minerals based on layered grinding and data correction, relating to the fields of mineral processing and geological testing and analysis. This method simulates the grinding process, acquiring embedment state data at different grinding depths in layers, and corrects errors using mathematical models, thereby constructing a refined embedment state analysis process for target minerals under macroscopic ore conditions. This method is standardized in operation, with strong data comparability. It specifies a complete standardized process from representative sample selection, marker positioning, precise layered grinding to automated mineralogical analysis, ensuring the consistency and repeatability of test results across different samples, batches, and operators, providing reliable technical support for mineral processing research and industrial production.
[0241] It should be noted that the present invention is not limited to the above-described embodiments. The above embodiments are merely examples, and any embodiments that have the same structure and perform the same effects as the technical concept within the scope of the present invention are included within the scope of the present invention. Furthermore, various modifications that can be conceived by those skilled in the art to the embodiments, and other ways of constructing by combining some of the constituent elements of the embodiments, without departing from the spirit of the present invention, are also included within the scope of the present invention.
Claims
1. A method for analyzing the mineral embedding state, characterized in that, Includes the following steps: S1, Regular Sample Preparation: Select representative ore samples, pre-process them to prepare regular samples; determine the initial observation surface as the 0th layer measurement surface; S2, Reference Surface Establishment and Measurement: Mineralogical analysis is performed on the 0th layer measurement surface to obtain the exposed area, embedding characteristic parameters, and reference values of the target mineral content in the state of fine particle encapsulation. S3, Layered grinding and measurement: The initial observation surface is ground in layers n times, with a grinding depth of Δh each time, to obtain the measurement surfaces from the 1st layer to the nth layer, where n≥2; mineralogical analysis is performed on each new measurement surface under the same conditions as the 0th layer measurement surface to obtain the exposed area and embedding characteristic parameters of the target mineral in each layer, until the target mineral particles tend to be encapsulated by gangue. S4, Gangue Intercalation Correction Coefficient Calculation: Based on the measured content of target minerals in the fine-particle encapsulated state in the actual grinding product, a gangue intercalation correction coefficient is established to correct the gangue intercalation data obtained by layered measurement. S5, Calculation of embedding state data: Based on the exposed area and embedding characteristic parameters of each layer of the measurement surface, combined with the gangue crystallization correction coefficient, the embedding state data of the target mineral in three-dimensional space is obtained through weighted calculation.
2. The mineral embedding state analysis method according to claim 1, characterized in that, The embedded characteristic parameter is the intercrystallization ratio of the target mineral with other minerals or pores, including the intercrystallization ratio of the target mineral with metal sulfides, metal oxides, gangue minerals and pores.
3. The mineral embedding state analysis method according to claim 2, characterized in that, In step S5, the weighted calculation adopts a weighted average algorithm based on the proportion of exposed area of each layer measurement surface, where Ce is the crystallization ratio of the target mineral to the e-th type of mineral phase or pores; where e=1 represents crystallization with metal sulfides; e=2 represents crystallization with metal oxides; e=3 represents crystallization with gangue minerals; and e=4 represents crystallization with pores. When e=1, 2, Ce=1 / (n+1)×Σ(Sn / S0×Cn.e); When e=3, Ce=T / (n+1)×Σ(Sn / S0×Cn.e); When e=4, C4=1-C1-C2-C3; Where Sn is the total exposed area of the target mineral on the nth measurement surface; S0 is the total exposed area of the target mineral on the 0th measurement surface; Cn.e is the crystallization ratio of the target mineral on the nth measurement surface to the eth type of mineral phase or pores; and T is the gangue crystallization correction coefficient.
4. The mineral embedding state analysis method according to claim 3, characterized in that, The target mineral in the fine-particle encapsulated state is a target mineral particle with a particle size range of 0.010~0.005mm; the correction coefficient T is the ratio of the content P' of the target mineral in the particle size range encapsulated by gangue in the actual grinding product of the representative ore sample to the content P of the target mineral in the particle size range encapsulated by gangue in the 0th layer measurement surface; P is the benchmark value of the target mineral content in the fine-particle encapsulated state, and the calculation formula is as follows: T=P' / P.
5. The mineral embedding state analysis method according to claim 3, characterized in that, The formula for calculating the crystallization ratio Cn.e of the target mineral and the e-th type of mineral phase or pores on the nth measurement surface is: ; Where n is the layer number of the measurement surface; m is the number of the target mineral particle, m=1, 2, 3..., representing the 1st, 2nd, 3rd... target mineral particles detected in the automated mineralogical analysis; Sm.n is the exposed area of the mth target mineral particle on the nth layer measurement surface; Cm.en is the proportion of the side length of the mth target mineral particle in contact with the eth type of mineral phase or pore on the nth layer measurement surface.
6. The mineral embedding state analysis method according to claim 1, characterized in that, In step S3, the grinding depth Δh is 1 / 6 to 1 / 2 of the average particle size of the target mineral.
7. The mineral embedding state analysis method according to claim 1, characterized in that, In step S3, the criteria for determining the state of being encased in gangue are: the exposed area of the target mineral particles on the current measurement surface is zero, and / or the equivalent particle size is not greater than 10% of the initial particle size.
8. The mineral embedding state analysis method according to claim 1, characterized in that, The mineral embedding state analysis method further includes an identification step before step S2, specifically: setting an identifier on a regular sample, placing it in a direction perpendicular to the measurement surface, with its width and height plane parallel to the measurement surface to establish a spatial reference system; the identifier is a metal sheet, made of copper or copper alloy, with dimensions of 0.1~0.5mm in length, 0.1~0.5mm in width, and 1.0~5.0mm in height.
9. The mineral embedding state analysis method according to claim 1, characterized in that, The target mineral is a metal sulfide, including one or more combinations of chalcopyrite, pyrite, galena, sphalerite, and molybdenite; the ore to be tested is a polymetallic complex ore, including copper-molybdenum polymetallic ore, copper-lead-zinc polymetallic ore, precious metal ore, or rare metal ore.
10. The mineral embedding state analysis method according to claim 1, characterized in that, The representative ore sample is a representative large ore block or drill core with a minimum particle size of not less than 4.0 cm. The measurement surface is ground, polished and carbonized before mineralogical analysis. Polishing uses 0.5~5.0 μm abrasive and carbonizing thickness is 10~20 nm.