The occurrence state of rubidium and the methods, equipment and media for analyzing the genetic origin of the Shimen rubidium deposit.

Abstract

This invention discloses a method, equipment, and medium for analyzing the occurrence state and genetic genesis of rubidium in the Shimen rubidium deposit, belonging to the field of ore genetic analysis technology. The method includes: identifying structurally anomalous samples by calculating the deviation between theoretical and measured structural parameters; distinguishing between first-strain mode samples and second-strain mode samples based on the divergence characteristics of strain spatial distribution; determining the first type of rubidium occurrence layer and the first type of enrichment mechanism by utilizing the spatial lag relationship between isotope fractionation peak positions and strain peak positions; determining the second type of rubidium occurrence layer and the second type of enrichment mechanism by identifying the abrupt changes in rubidium concentration gradient around fluid inclusions; and delineating the mineralized region by coupling the enrichment direction at the sample scale with fracture structures or permeable strata. This invention achieves quantitative determination of the lattice position of rubidium occurrence, differentiation of the kinetic processes of mineralization mechanisms, and delineation of the geological structural constraints of mineralized regions.

Description

Technical Field

[0001] This invention relates to the field of ore deposit genetic analysis technology, specifically to the occurrence state of rubidium in the Shimen rubidium deposit and the methods, equipment and media for ore deposit genetic analysis. Background Technology

[0002] Rubidium exists in mica lattices in two ways: interlayer sites and octahedral sites. The extraction difficulty varies significantly between these two site types. Current techniques use electron probe microanalysis to determine rubidium content, but cannot distinguish the specific site type. X-ray diffraction analysis assumes equilibrium crystallization and does not consider non-equilibrium effects in explosive mineralization, thus failing to extract site information from lattice structure deviations.

[0003] Rubidium deposits exist in two types: explosive mineralization and infiltration mineralization. Current genetic analyses rely on inferences from geological background and lack the technical means to quantitatively distinguish mineralization mechanisms.

[0004] Mineralization zone delineation is based on interpolation of known orebody boundaries, without considering the spatial propagation of mineralization dynamics. Explosive mineralization is controlled by fault structures, while permeable mineralization is controlled by stratigraphic permeability. Current technologies lack methods to couple sample-scale enrichment characteristics with regional-scale geological structures, resulting in insufficient accuracy in mineral exploration prediction. Summary of the Invention

[0005] In view of the above-mentioned problems, the present invention provides a method, equipment and medium for analyzing the occurrence state of rubidium and the genetic formation of the Shimen rubidium deposit.

[0006] Therefore, the technical problem solved by this invention is: how to identify non-equilibrium crystal samples by means of mica lattice structure deviation, to determine the occurrence layer of rubidium in different lattice positions and to quantitatively distinguish between explosive and permeation enrichment mechanisms, and at the same time to associate the spatial propagation characteristics of rubidium enrichment with the direction of tectonic stress field or the direction of formation permeation in order to delineate the mineralization area.

[0007] To solve the above-mentioned technical problems, the present invention provides the following technical solution: a method for analyzing the occurrence state and genetic genesis of rubidium in the Shimen rubidium deposit, comprising,

[0008] Obtain rubidium-containing mica samples, calculate theoretical structural parameters based on the elemental composition of the rubidium-containing mica samples, obtain measured structural parameters through diffraction testing, and identify structurally abnormal samples based on the deviation between theoretical and measured structural parameters.

[0009] The spatial distribution characteristics of strain are obtained by analyzing the diffraction signal morphology of structurally anomalous samples. Based on the dominant direction of strain spatial distribution, the structurally anomalous samples are divided into first strain mode samples and second strain mode samples.

[0010] The isotopic composition of the first strain mode sample was spatially sampled and determined along the dominant direction. Based on the spatial variation characteristics of the isotopic composition, the first type of rubidium occurrence layer and the first type of enrichment mechanism were determined.

[0011] Based on the type I rubidium occurrence layer and enrichment mechanism, the spatial distribution of the samples was correlated with the direction of the tectonic stress field to delineate the type I mineralization region;

[0012] The directional characteristics of fluid inclusions in the second strain mode sample were observed and obtained. Based on the matching of the directional characteristics with the anisotropy of the formation, the second type of rubidium occurrence layer and the second type of enrichment mechanism were determined.

[0013] Based on the second type of rubidium occurrence horizon and enrichment mechanism, the spatial distribution of the samples was correlated with the formation permeability direction to delineate the second type of mineralization region.

[0014] As a preferred embodiment of the method for analyzing the occurrence state and genetic genesis of rubidium in the Shimen rubidium deposit described in this invention, the calculation of theoretical structural parameters based on the elemental composition of the rubidium-bearing mica sample includes:

[0015] The rubidium and potassium content of the rubidium-containing mica sample was determined.

[0016] A database of rubidium-potassium mica standard samples in equilibrium was established, and the relationship curve between rubidium mole fraction and interplanar spacing was extracted.

[0017] A dataset of experimental relationships between cooling rate and interlayer position occupancy ratio was established through high temperature and high pressure tests. The cooling rate of ore-forming fluid was estimated based on the uniform temperature spatial variation gradient of fluid inclusions.

[0018] The theoretical structural parameters are determined based on the rubidium content, potassium content, interlayer site occupancy ratio, and relationship curves.

[0019] As a preferred embodiment of the method for analyzing the occurrence state and genetic genesis of the Shimen rubidium deposit described in this invention, the method for obtaining measured structural parameters through diffraction testing includes:

[0020] X-ray diffraction tests were performed on rubidium-containing mica samples to identify the first, second, and third diffraction order secondary peaks of the mica's layered structure.

[0021] The deviation of the interplanar spacing corresponding to the three diffraction order sub-peaks is used to identify lattice periodic inhomogeneity. In response to the presence of lattice periodic inhomogeneity, the measured structural parameters are determined according to the inhomogeneous lattice correction method. In response to the absence of lattice periodic inhomogeneity, the measured structural parameters are determined from the first diffraction order sub-peak.

[0022] As a preferred embodiment of the method for analyzing the occurrence state and genetic genesis of the Shimen rubidium deposit described in this invention, the method of dividing the structurally abnormal samples into first strain mode samples and second strain mode samples according to the dominant direction of strain spatial distribution includes:

[0023] Spatially distributed thin sections were prepared for samples with structural anomalies and observed using a polarizing microscope. The interference color order of each observation point in the observation point array was recorded and converted into birefringence values.

[0024] Determine the dominant gradient direction and intensity at each observation point, and construct a gradient vector field distribution map;

[0025] The gradient vector field distribution map is integrated in the divergence space. If the absolute value of the divergence space integral is greater than the radial mode divergence threshold, the structurally anomalous sample is identified as a first strain mode sample. If the absolute value of the divergence space integral is less than or equal to the radial mode divergence threshold, the structurally anomalous sample is identified as a second strain mode sample.

[0026] As a preferred embodiment of the method for analyzing the occurrence state and genetic genesis of the Shimen rubidium deposit described in this invention, the method for determining the first type of rubidium occurrence stratigraphy and the first type of enrichment mechanism based on the spatial variation characteristics of isotopic composition includes:

[0027] A polar coordinate system was established for the first strain mode sample and multiple radiation lines were set. The rubidium 87 abundance and rubidium 85 abundance at each measurement position on each radiation line were measured.

[0028] The locations of strain peaks and isotope fractionation peaks were identified for each radiation survey line, and the first type of rubidium-bearing horizons were determined based on the statistical distribution of the radial offset of the peak locations.

[0029] The azimuth distribution characteristics of isotope attenuation amplitude were extracted from each radiation line, and the dominant propagation direction of rubidium enrichment and the distribution range of isotope attenuation amplitude were determined as the first type of enrichment mechanism.

[0030] As a preferred embodiment of the method for analyzing the occurrence state and genetic genesis of the Shimen rubidium deposit described in this invention, the step of delineating the first type of mineralization region by correlating the spatial distribution of the sample with the direction of the tectonic stress field based on the first type of rubidium occurrence strata and enrichment mechanism includes:

[0031] Identify samples with the first strain mode, compare the rubidium enrichment-dominant propagation direction with the fracture structure orientation to identify fracture-related samples, project the fracture-related samples onto the fracture structure and cluster them to identify the active fluid transport segments of the fracture structure;

[0032] A spatial prediction domain for fluid migration is established along the active fluid migration segment of the fracture structure. The three-dimensional extension range of the spatial prediction domain for fluid migration is determined based on the distribution range of isotope attenuation amplitude, the geometric parameters of the active fluid migration segment of the fracture structure, and the dip angle of the fracture structure, and a first-type spatial envelope is formed.

[0033] Within the fluid transport space prediction domain defined by the first type of spatial envelope, rubidium preferential enrichment layers are identified and first type mineralization regions are delineated based on the first type of rubidium occurrence layers.

[0034] As a preferred embodiment of the method for analyzing the occurrence state and genetic genesis of the Shimen rubidium deposit described in this invention, the determination of the second type of rubidium occurrence horizon and the second type of enrichment mechanism based on directional characteristics and stratigraphic anisotropy matching includes:

[0035] Transmission electron microscopy was used to identify the lattice distortion region by imaging the mica lattice around the fluid inclusions in the sample of the second strain mode, and the rubidium concentration distribution within the radius of influence of the fluid inclusions was determined and the location of the abrupt change in the rubidium concentration gradient was identified.

[0036] The distribution pattern of fluid inclusions with abrupt changes in rubidium concentration gradients within mica crystals was used to determine the second type of rubidium occurrence sites.

[0037] The distribution of vector directions pointing from the center of the fluid inclusions to the location of abrupt changes in the rubidium concentration gradient is used to obtain the dominant orientation of rubidium migration in the region, and the spatial distribution characteristics of the locations of abrupt changes in the rubidium concentration gradient and the geometric characteristics of the lattice distortion region are extracted.

[0038] In response to the matching relationship between the dominant rubidium migration orientation and the stratigraphic strike satisfying the stratigraphic permeability consistency condition, the dominant rubidium migration orientation and fluid inclusion interface reaction characteristic parameters are identified as the second type of enrichment mechanism.

[0039] As a preferred embodiment of the method for analyzing the occurrence state and genetic genesis of the Shimen rubidium deposit described in this invention, the step of delineating the second type of mineralization region by correlating the spatial distribution of samples with the stratigraphic permeability direction based on the second type of rubidium occurrence horizon and enrichment mechanism includes:

[0040] Identify samples with the second strain mode, compare the dominant orientation of regional rubidium migration with the strike of the high-permeability stratigraphic unit to identify stratigraphic permeability-related samples, project the stratigraphic permeability-related samples onto the high-permeability stratigraphic unit and continuously analyze to identify active fluid permeable stratigraphic segments;

[0041] Spatial interpolation is performed on the fluid inclusion interface reaction characteristic parameters of permeability-associated samples in active fluid-permeable strata to identify active interface reaction regions and establish a spatial prediction domain for fluid permeability.

[0042] Based on the spatial geometric characteristics of the active interface reaction region and the geometric parameters of the active fluid-permeable stratum, the three-dimensional extension range of the fluid permeability spatial prediction domain is determined and a second type of spatial envelope is formed.

[0043] Within the fluid permeation space prediction domain defined by the second type of spatial envelope, rubidium preferential enrichment layers are identified and second type mineralization regions are delineated based on the second type of rubidium occurrence layers.

[0044] The present invention provides a computer device, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps of the method for analyzing the occurrence state of rubidium and the genetic formation of the Shimen rubidium deposit.

[0045] The present invention provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the steps of the method for analyzing the occurrence state of rubidium in the Shimen rubidium deposit and the genetic formation of the deposit are implemented.

[0046] The beneficial effects of this invention are as follows: This invention identifies structurally anomalous samples by the deviation between theoretical and measured structural parameters; it distinguishes between radial divergence and parallel flow modes based on the divergence topological characteristics of strain spatial distribution; it determines the interlayer or octahedral occurrence layer of rubidium by utilizing the spatial lag relationship between the peak position of isotope fractionation and the peak position of strain; it identifies the interfacial reaction enrichment mechanism by identifying the abrupt change in the rubidium concentration gradient around fluid inclusions; and it establishes a spatial prediction domain by coupling the enrichment direction information at the sample scale with the fracture structure or permeable strata at the regional scale. This invention achieves quantitative determination of the lattice position of rubidium occurrence, differentiation of the dynamic process of mineralization mechanism, and delineation of the geological structural constraints of mineralized areas, providing a quantitative technical method for rubidium deposit prospecting and prediction. Attached Figure Description

[0047] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0048] Figure 1 This is a flowchart illustrating the overall process of analyzing the occurrence state of rubidium and the genetic origin of the Shimen rubidium deposit, as provided in an embodiment of the present invention.

[0049] Figure 2 A flowchart illustrating the structural anomaly sample classification method for analyzing the occurrence state of rubidium and the genetic origin of the Shimen rubidium deposit, as provided in an embodiment of the present invention.

[0050] Figure 3This is a flowchart illustrating the first type of mineralization region delineation process for the analysis of the occurrence state of rubidium and the genetic formation of the Shimen rubidium deposit, as provided in an embodiment of the present invention. Detailed Implementation

[0051] To make the present invention more apparent and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the protection scope of the present invention.

[0052] Example 1, referring to Figures 1-3 This is one embodiment of the present invention, which provides a method for analyzing the occurrence state of rubidium and the genetic origin of the Shimen rubidium deposit, including:

[0053] S1: Obtain rubidium-containing mica samples, calculate theoretical structural parameters based on the elemental composition of the rubidium-containing mica samples, obtain measured structural parameters through diffraction testing, and identify structurally abnormal samples based on the deviation between theoretical and measured structural parameters.

[0054] S2: Obtain the strain spatial distribution characteristics by analyzing the diffraction signal morphology of the structurally anomalous sample, and divide the structurally anomalous sample into a first strain mode sample and a second strain mode sample according to the dominant direction of the strain spatial distribution.

[0055] S3: Spatial sampling and isotopic composition determination of the first strain mode sample along the dominant direction, and determination of the first type of rubidium occurrence layer and the first type of enrichment mechanism based on the spatial variation characteristics of isotopic composition.

[0056] S4: Based on the first type of rubidium occurrence layer and enrichment mechanism, the spatial distribution of the sample is correlated with the direction of the tectonic stress field to delineate the first type of mineralization region;

[0057] S5: Observe the fluid inclusions of the second strain mode sample to obtain directional characteristics, and determine the second type of rubidium occurrence layer and the second type of enrichment mechanism based on the directional characteristics and the anisotropy of the formation.

[0058] S6: Based on the second type of rubidium occurrence horizon and enrichment mechanism, the spatial distribution of the samples is correlated with the formation permeability direction to delineate the second type of mineralization area.

[0059] In some embodiments, step S1 specifically includes the following steps S11-S19:

[0060] S11: Rubidium-bearing mica samples were collected from different mineralized parts of the Shimen Rubidium Deposit. After crushing the rubidium-bearing mica samples, mica particles were selected and made into polished thin sheets.

[0061] S12: Elemental composition analysis of polished thin sections was performed using an electron probe microanalysis device. Multiple measurement points were selected on each mica particle to determine the rubidium and potassium content. The rubidium and potassium content of the mica particles were calculated based on the measurement results of multiple measurement points. The rubidium mole fraction of the rubidium-containing mica sample was calculated based on the rubidium and potassium content.

[0062] S13: Establish a database of rubidium-potassium mica standard samples in equilibrium. The standard samples in the database are grown under constant temperature and pressure conditions to achieve thermodynamic equilibrium in chemical composition and crystal structure. Extract the relationship curve between the rubidium mole fraction and interplanar spacing of the standard samples from the database.

[0063] S14: The explosive mineralization event in the rubidium deposit led to the rapid growth of mica crystals during the cooling process of the ore-forming fluid. The distribution of rubidium in the interlayer and octahedral positions of the mica lattice did not reach the equilibrium distribution ratio. Mica samples were synthesized under different cooling rates through high-temperature and high-pressure experiments, and the occupancy ratio of rubidium in the interlayer positions of the synthesized mica samples was measured. An experimental relationship dataset between cooling rate and interlayer position occupancy ratio was established. The cooling rate of the ore-forming fluid was estimated based on the spatial variation gradient of the homogenization temperature of the fluid inclusions in the rubidium deposit. The interlayer position occupancy ratio of the cooling rate of the ore-forming fluid was extracted from the experimental relationship dataset between cooling rate and interlayer position occupancy ratio.

[0064] The cooling rate is positively correlated with the proportion of rubidium in the interlayer.

[0065] S15: Analyze a series of mica samples with increasing rubidium content using crystal structure refinement methods and determine the critical rubidium mole fraction at which rubidium begins to enter the octahedral position. If the rubidium mole fraction of the rubidium-containing mica sample exceeds the critical rubidium mole fraction, calculate the effective rubidium mole fraction between layers based on the rubidium mole fraction of the rubidium-containing mica sample and the interlayer position occupancy ratio. If the rubidium mole fraction of the rubidium-containing mica sample does not exceed the critical rubidium mole fraction, the rubidium mole fraction of the rubidium-containing mica sample is the effective rubidium mole fraction between layers. Find the interplanar spacing value on the relationship curve between the rubidium mole fraction and the interplanar spacing based on the effective rubidium mole fraction between layers as a theoretical structural parameter.

[0066] The effective molar fraction of rubidium between layers determines the interplanar spacing of mica in equilibrium state.

[0067] S16: Grind the rubidium-containing mica sample into powder and sieve it to obtain powder components with different particle size ranges. Perform X-ray diffraction tests on the powder components and identify the first, second, and third diffraction order peaks of the mica layered structure in the diffraction pattern. Calculate the first interplanar spacing from the peak position of the first diffraction order peak, the second interplanar spacing from the peak position of the second diffraction order peak, and the third interplanar spacing from the peak position of the third diffraction order peak.

[0068] The first diffraction order peak of the mica layered structure corresponds to the (001) crystal plane family, the second diffraction order peak corresponds to the (002) crystal plane family, and the third diffraction order peak corresponds to the (003) crystal plane family. The interplanar spacing is calculated from the diffraction angle using Bragg's law.

[0069] S17: During the non-equilibrium crystallization process of rubidium-bearing mica, the interplanar spacing changes at different stages of crystal growth, resulting in systematic differences between the first, second, and third interplanar spacings. The second diffraction order peak reflects the cumulative effect of two lattice periods, and the third diffraction order peak reflects the cumulative effect of three lattice periods. The lattice period inhomogeneity is identified by comparing the deviations between the first, second, and third interplanar spacings. If the deviations between the first, second, and third interplanar spacings exceed the lattice period deviation threshold, the rubidium-bearing mica sample is determined to have lattice period inhomogeneity. The first, second, and third interplanar spacings are then combined according to the inhomogeneous lattice correction method to obtain the measured structural parameters. If the deviations between the first, second, and third interplanar spacings do not exceed the lattice period deviation threshold, the rubidium-bearing mica sample is determined to not have lattice period inhomogeneity, and the first interplanar spacing is used as the measured structural parameter.

[0070] Theoretically, the second interplanar spacing is twice that of the first interplanar spacing, and the third interplanar spacing is three times that of the first interplanar spacing. The deviation of the actual measured value reflects the degree of non-uniformity of the lattice period.

[0071] S18: The normalized structural deviation index is obtained by comparing the measured structural parameters with the theoretical structural parameters. The normalized structural deviation index of all rubidium-containing mica samples in the batch is statistically analyzed and the distribution characteristics are obtained. The structural anomaly judgment boundary is determined based on the distribution characteristics. In response to the normalized structural deviation index of the rubidium-containing mica sample exceeding the structural anomaly judgment boundary, the rubidium-containing mica sample is identified as a structural anomaly sample.

[0072] The normalized structural deviation index eliminates the influence of the difference in the absolute value of interplanar spacing between samples with different rubidium contents.

[0073] Through the above step S1, the present invention calculates theoretical structural parameters from elemental composition, obtains measured structural parameters through diffraction testing, and identifies structurally abnormal samples of non-equilibrium crystallization based on the deviation between the two. The sample identification success rate is about 50% higher than the existing single-crystal structure refinement method, providing a sample basis for subsequent mineralization mechanism identification and occurrence layer determination.

[0074] In some embodiments, step S2 specifically includes the following steps S21-S25:

[0075] S21: Spatial distribution thin films are prepared by cutting the structurally abnormal sample along the direction perpendicular to the mica cleavage plane. The spatial distribution thin films are observed using a polarizing microscope in orthogonal polarization mode. A two-dimensional coordinate system is established on the surface of the spatial distribution thin films. An observation point array is set in the two-dimensional coordinate system. The observation point array includes multiple observation point positions. The interference color order of each observation point position in the observation point array is recorded.

[0076] S22: Based on the relationship between the interference color order and the optical path difference, convert the interference color order of each observation point in the observation point array into the optical path difference value, and calculate the birefringence value of each observation point in the observation point array based on the optical path difference value and the thickness of the spatially distributed thin film.

[0077] Birefringence is equal to the optical path difference divided by the thickness of the sheet.

[0078] S23: Calculate the dominant gradient direction and dominant gradient intensity for each observation point in the observation point array. For any observation point in the array, extract the birefringence value of that observation point and the birefringence values ​​of its neighboring observation points (the observation points adjacent to the observation point in the two-dimensional coordinate system). Calculate the difference between the birefringence value of the observation point and the birefringence values ​​of each neighboring observation point. Based on the difference in birefringence values ​​and the spatial distance between the observation point and its neighboring observation points, calculate the spatial rate of change of birefringence in each spatial direction. Identify the spatial direction with the largest spatial rate of change of birefringence as the dominant gradient direction of the observation point, and take the maximum value of the spatial rate of change of birefringence as the dominant gradient intensity of the observation point. Obtain the dominant gradient direction and dominant gradient intensity for each observation point in the array by calculating them one by one.

[0079] S24: Construct a gradient vector field distribution map based on the dominant gradient direction and dominant gradient intensity of each observation point in the observation point array, calculate the vector field divergence distribution of the gradient vector field distribution map, and perform spatial integration of the vector field divergence distribution within the measurement area of ​​the two-dimensional coordinate system to obtain the divergence spatial integral value.

[0080] The divergence calculation employs a discrete point vector field divergence algorithm. For any observation point in a two-dimensional coordinate system, the dominant gradient direction and intensity at that observation point are extracted to form a gradient vector. The gradient vectors of neighboring observation points are also extracted. The rate of change of the gradient vector in the spatial direction is calculated. The rate of change is then linearly integrated over a closed path around the observation point and divided by the area enclosed by the closed path to obtain the divergence value at that observation point. The divergence values ​​for all observation points are summed over the measurement region to obtain the spatial integral value of the divergence. A positive divergence value in the radially divergent field corresponds to the strain field source characteristics caused by the outward migration of explosive fluid from the center. A negative divergence value in the radially converging field corresponds to the strain field sink characteristics caused by the convergence of fluid towards the center. A divergence value close to zero in the parallel flow field corresponds to the parallel propagation characteristics of the strain field caused by the directional migration of permeable fluid along the strata.

[0081] S25: Calculate the absolute value of the divergence space integral. If the absolute value of the divergence space integral is greater than the radial mode divergence threshold, determine that the dominant direction of the spatial distribution characteristics of strain is the radial divergence mode and identify the structurally anomalous sample as the first strain mode sample. If the absolute value of the divergence space integral is less than or equal to the radial mode divergence threshold, determine that the dominant direction of the spatial distribution characteristics of strain is the parallel flow mode and identify the structurally anomalous sample as the second strain mode sample.

[0082] Through the above step S2, the present invention quantitatively distinguishes between explosive mineralization and permeation mineralization mechanisms based on the divergence topological characteristics of strain spatial distribution, transforming the mineralization mechanism discrimination from qualitative description to quantitative calculation, and providing mineralization mechanism constraints for subsequent determination of rubidium occurrence strata and delineation of mineralized areas.

[0083] In some embodiments, step S3 specifically includes the following steps S31-S38:

[0084] S31: Determine the spatial distribution of rubidium isotopes on the first strain mode sample. Establish a polar coordinate system on the surface of the first strain mode sample. The origin of the polar coordinate system is located at the radial divergence center identified in step S2. Set multiple radiation measurement lines in the polar coordinate system. The radiation measurement lines extend outward from the origin of the polar coordinate system. The multiple radiation measurement lines are evenly distributed in the azimuth direction. Set multiple measurement positions on each radiation measurement line. The measurement positions are distributed sequentially outward from the origin of the polar coordinate system along the radiation measurement lines.

[0085] S32: Measure the rubidium 87 abundance and rubidium 85 abundance at each measurement location along each radiation line, and calculate the ratio of rubidium 87 abundance to rubidium 85 abundance at each measurement location as the isotopic abundance ratio.

[0086] S33: Extract the birefringence values ​​corresponding to each measurement position on each radiation measurement line from step S2, calculate the difference in birefringence values ​​between adjacent measurement positions on each radiation measurement line, and calculate the radial change rate of birefringence based on the difference in birefringence values ​​and the radial distance between adjacent measurement positions.

[0087] S34: Calculate the difference in isotopic abundance ratio between adjacent measurement locations on each radiation measurement line, and calculate the radial change rate of the isotopic abundance ratio based on the difference in isotopic abundance ratio and the radial distance between adjacent measurement locations.

[0088] S35: The extreme values ​​of the distribution of the radial change rate of birefringence along the radial distance on each radiation measurement line are identified, and the measurement position corresponding to the maximum value of the radial change rate of birefringence on each radiation measurement line is extracted as the strain peak position of each radiation measurement line. The extreme values ​​of the distribution of the radial change rate of isotope abundance ratio along the radial distance on each radiation measurement line are identified, and the measurement position corresponding to the maximum value of the radial change rate of isotope abundance ratio on each radiation measurement line is extracted as the isotope fractionation peak position of each radiation measurement line.

[0089] S36: Calculate the radial distance difference between the peak position of isotope fractionation and the peak position of strain on each radiation line as the radial offset of the peak position of each radiation line. Statistically analyze the radial offsets of the peak positions of multiple radiation lines and calculate the median of the radial offsets. If the median of the radial offsets of the peak positions is greater than the spatial lag determination threshold, the interlayer position is determined as a type I rubidium occurrence layer. If the median of the radial offsets of the peak positions is less than or equal to the spatial lag determination threshold, the octahedral position is determined as a type I rubidium occurrence layer.

[0090] When rubidium occupies interlayer positions, it lies between mica crystal layers. Strain propagation causes relative displacement of the crystal layers. The van der Waals forces at the interlayer positions are insufficient to constrain the movement of rubidium atoms, which then diffuse and migrate along the strain propagation direction. During diffusion, the mass difference between rubidium-87 and rubidium-85 leads to different migration activation energies. Due to its smaller mass, rubidium-85 preferentially migrates towards the diffusion front, while rubidium-87 becomes enriched behind the diffusion path. The most dramatic change in the isotopic abundance ratio occurs at the diffusion front, forming the isotopic fractionation peak position. The strain peak position corresponds to the spatial position with the largest relative displacement of the crystal layers. The time required for diffusion leads to the isotopic fractionation peak. The position lags behind the strain peak position. The spatial lag distance is related to the diffusion coefficient and strain propagation rate. When rubidium occupies the octahedral position, it is located at the oxygen-coordinated octahedral center of the mica lattice framework. The coordination bond energy of the octahedral position is higher than that of the van der Waals force at the interlayer position. During strain propagation, rubidium atoms are bound by the octahedral coordination bonds and undergo coordinated deformation with the lattice framework. Rubidium atoms cannot diffuse and migrate. Isotope fractionation occurs synchronously with strain at the strain peak position. The radial offset of the peak position is close to zero. The statistical distribution of the radial offset of the peak position can be used to distinguish between the interlayer position diffusion mechanism and the octahedral position binding mechanism.

[0091] S37: Calculate the isotopic abundance ratio at the measurement positions near the origin of the polar coordinate system on each radiation measurement line as the isotopic abundance ratio at the center of the radiation measurement line, calculate the isotopic abundance ratio at the measurement positions far from the origin of the polar coordinate system on each radiation measurement line as the isotopic abundance ratio at the edge of the radiation measurement line, and calculate the difference between the isotopic abundance ratio at the center of the radiation measurement line and the isotopic abundance ratio at the edge of the radiation measurement line as the isotopic attenuation amplitude of the radiation measurement line.

[0092] S38: Arrange the isotope attenuation amplitudes of multiple radiometric lines in azimuth order to form an isotope attenuation amplitude azimuth distribution sequence. Extract the maximum value of the isotope attenuation amplitude of the radiometric lines and its corresponding azimuth from the isotope attenuation amplitude azimuth distribution sequence as the dominant propagation direction of rubidium enrichment. Extract the minimum value of the isotope attenuation amplitude of the radiometric lines from the isotope attenuation amplitude azimuth distribution sequence. Calculate the difference between the maximum value and the minimum value of the isotope attenuation amplitude of the radiometric lines as the distribution range of the isotope attenuation amplitude of the radiometric lines. The dominant propagation direction of rubidium enrichment and the distribution range of the isotope attenuation amplitude of the radiometric lines are jointly determined as the first type of enrichment mechanism.

[0093] Through the above step S3, the present invention uses the spatial lag relationship between the peak position of isotope fractionation and the peak position of strain to determine the occurrence site of rubidium in the mica lattice, and the applicable sample range is significantly expanded compared with the existing single crystal structure refinement method.

[0094] In some embodiments, step S4 specifically includes the following steps S41-S46:

[0095] S41: Extract the spatial coordinates of multiple rubidium-bearing mica samples collected within the mining area, identify the rubidium-bearing mica samples determined to be first strain mode samples as first strain mode samples, and extract the first type of rubidium occurrence layer, the dominant propagation direction of rubidium enrichment, and the distribution range of isotope attenuation amplitude of radiometric lines corresponding to each first strain mode sample.

[0096] S42: Obtain structural geological data within the mining area, including the spatial location and occurrence parameters of fault structures. Extract fault structures from the structural geological data as potential fluid channels. Calculate the angle between the dominant rubidium enrichment propagation direction of each first strain mode sample and the strike of the nearest fault structure as the angle between the sample propagation direction and the fault strike.

[0097] S43: The first strain mode sample whose angle between the sample propagation direction and the fracture direction is less than the structural consistency angle judgment value is identified as the fracture-related sample. The spatial coordinate position of each fracture-related sample is vertically projected onto the nearest fracture structure to obtain the projection position of each fracture-related sample on the fracture structure. The multiple projection positions on the fracture structure are clustered according to the spatial distribution along the fracture direction. The fracture segments with spatial clustering of projection positions are identified as the active fluid transport segments of the fracture structure.

[0098] S44: Establish a fluid migration spatial prediction domain along the active fluid migration segment of the fracture structure. Extract the isotope attenuation amplitude distribution range values ​​of the radiometric isotope attenuation amplitude distribution range of the fracture-related samples. Calculate the maximum value of the isotope attenuation amplitude distribution range values ​​of the radiometric isotope attenuation amplitude distribution range of the fracture-related samples as the maximum attenuation amplitude distribution range. Extend the maximum attenuation amplitude distribution range from the active fluid migration segment of the fracture structure to both sides in a direction perpendicular to the fracture structure strike to form the lateral extension range of the fluid migration spatial prediction domain. Extend the active fluid migration segment of the fracture structure outward from both ends in a direction along the fracture structure strike to form the longitudinal extension range of the fluid migration spatial prediction domain. The extension length of the longitudinal extension range is determined according to the length of the active fluid migration segment of the fracture structure. Extend the active fluid migration segment of the fracture structure towards the hanging wall and footwall of the fracture structure in a direction along the fracture structure dip to form the dip extension range of the fluid migration spatial prediction domain. The dip extension range is determined according to the dip angle of the fracture structure.

[0099] S45: The lateral extension range, longitudinal extension range and directional extension range of the fluid transport spatial prediction domain are combined in three-dimensional space to form a first type of spatial envelope, which defines the boundary of the fluid transport spatial prediction domain in three-dimensional space.

[0100] S46: Stratification is performed within the fluid migration spatial prediction domain defined by the first type of spatial envelope. In response to the first type of rubidium occurrence layer being an interlayer location, the mica interlayer location is identified as a rubidium-preferred enrichment layer within the fluid migration spatial prediction domain. In response to the first type of rubidium occurrence layer being an octahedral location, the mica octahedral location is identified as a rubidium-preferred enrichment layer within the fluid migration spatial prediction domain. The fluid migration spatial prediction domain is thus delineated as the first type of mineralization region.

[0101] Through the above step S4, the present invention couples the rubidium enrichment characteristics at the sample scale with the fracture structure at the regional scale, establishes a spatial prediction domain for fluid migration and delineates the first type of mineralization region, thereby achieving the geological structural constraint delineation of the mineralization region.

[0102] In some embodiments, step S5 specifically includes the following steps S51-S59:

[0103] S51: Prepare transparent thin films for the second strain mode sample, use a microscope to locate the fluid inclusions inside the rubidium-containing mica crystal in the transparent thin film, record the spatial coordinates of each fluid inclusion in the mica crystal, and perform transmission electron microscopy imaging on the mica lattice around each fluid inclusion to obtain the atomic-scale structural information of the mica lattice.

[0104] S52: Identify lattice distortion regions from transmission electron microscopy images of mica lattices around each fluid inclusion. Lattice distortion regions are areas where the lattice arrangement deviates from the ideal periodic structure. Measure the distance from the boundary of the fluid inclusion to the boundary of the lattice distortion region as the fluid inclusion influence radius. The fluid inclusion influence radius reflects the spatial range of interfacial reactions between the fluid and mica.

[0105] S53: Perform elemental line scan analysis on mica within the influence radius of the fluid inclusion. The elemental line scan analysis is performed along a radial path from the center of the fluid inclusion away from the fluid inclusion. Measure the rubidium concentration at each spatial location along the radial path and plot the distribution curve of the rubidium concentration along the radial distance.

[0106] S54: Perform gradient analysis on the rubidium concentration distribution curve, calculate the spatial rate of change of rubidium concentration along the radial distance, identify the extreme position of the spatial rate of change as the abrupt change position of the rubidium concentration gradient, and mark the spatial position where the fluid-mica interface reaction is most active.

[0107] S55: The proportion of fluid inclusions with abrupt changes in rubidium concentration gradient in the second strain mode sample is used as the proportion of interfacial reaction inclusions. If the proportion of interfacial reaction inclusions is greater than the interfacial reaction determination threshold, it is determined that rubidium mainly originates from the fluid and enters the mica lattice through the interfacial reaction between the fluid and mica.

[0108] S56: Perform spatial location statistics on fluid inclusions with abrupt changes in rubidium concentration gradients, analyze the distribution pattern of these fluid inclusions within the mica crystal, and determine that rubidium preferentially accumulates at interlayer locations through interfacial reactions, based on the fact that fluid inclusions with abrupt changes in rubidium concentration gradients are concentrated in the interlayer locations of the mica crystal. The interlayer interfacial reaction enrichment mode is identified as a type II rubidium occurrence site. In response to the fact that fluid inclusions with abrupt changes in rubidium concentration gradients are uniformly distributed within the mica crystal, determine that rubidium undergoes interfacial reaction enrichment at all locations in the mica crystal, and the whole-lattice interfacial reaction enrichment mode is identified as a type II rubidium occurrence site.

[0109] S57: Perform directional analysis on fluid inclusions with abrupt changes in rubidium concentration gradient in the second strain mode sample. Calculate the vector direction from the center of each fluid inclusion to the abrupt change in rubidium concentration gradient as the dominant rubidium migration direction of the fluid inclusion. Statistically analyze the spatial distribution of the dominant rubidium migration directions of multiple fluid inclusions. Divide the dominant rubidium migration directions into multiple azimuth intervals according to azimuth angle. Calculate the number of dominant rubidium migration directions in each azimuth interval. Identify the azimuth interval with the most dominant rubidium migration directions as the concentrated azimuth interval. Calculate the central azimuth angle of the concentrated azimuth interval as the dominant rubidium migration azimuth in the region.

[0110] S58: Statistically analyze the distances from the locations of abrupt changes in the rubidium concentration gradient of multiple fluid inclusions to the center of the fluid inclusions, calculate the median distance as the characteristic response distance, measure the area of ​​the lattice distortion region around multiple fluid inclusions, and calculate the average area of ​​the lattice distortion region as the average distortion area.

[0111] S59: Obtain stratigraphic attitude data within the mining area, extract the strike of stratigraphic bedding planes from the stratigraphic attitude data, calculate the angle between the dominant rubidium migration orientation and the strike of stratigraphic bedding planes as the angle between rubidium migration and stratigraphic strike, and determine that the rubidium migration direction is controlled by stratigraphic permeability anisotropy when the angle between rubidium migration and stratigraphic strike is less than the stratigraphic permeability consistency angle threshold. The dominant rubidium migration orientation, characteristic response distance, and average distortion area of ​​the region are jointly determined as the second type of enrichment mechanism.

[0112] Through the above step S5, the present invention identifies the interfacial reaction enrichment mechanism by the location of abrupt changes in rubidium concentration gradient around the fluid inclusions, and determines the formation ore-controlling characteristics based on the matching relationship between the rubidium migration direction and the formation permeation direction.

[0113] In some embodiments, step S6 specifically includes the following steps S61-S69:

[0114] S61: Extract the spatial coordinates of multiple rubidium-bearing mica samples collected within the mining area, identify the rubidium-bearing mica samples that are determined to be second strain mode samples as second strain mode samples, and extract the second type of rubidium occurrence layer, the dominant orientation of regional rubidium migration, the characteristic response distance and the average distortion area corresponding to each second strain mode sample.

[0115] S62: Obtain stratigraphic attitude data and stratigraphic permeability data within the mining area. The stratigraphic attitude data includes the strike, dip, and dip angle of the strata. The stratigraphic permeability data includes the permeability test results of different lithological units. Identify high-permeability stratigraphic units from the stratigraphic attitude data and stratigraphic permeability data. High-permeability stratigraphic units are stratigraphic units whose permeability exceeds the lower permeability threshold.

[0116] S63: Perform geometric relationship analysis between the dominant rubidium migration orientation of each second strain mode sample and the bedding strike of the high-permeability stratigraphic unit. Calculate the angle between the dominant rubidium migration orientation of the region and the bedding strike of the high-permeability stratigraphic unit as the angle between the rubidium migration direction of the sample and the stratigraphic strike. Identify second strain mode samples whose angle between the rubidium migration direction of the sample and the stratigraphic strike is less than the stratigraphic permeability consistency angle judgment value as stratigraphic permeability associated samples.

[0117] S64: Vertically project the spatial coordinates of the permeability-associated samples from each formation onto the nearest high-permeability formation unit to obtain the projected positions of the permeability-associated samples on the high-permeability formation unit. Perform continuity analysis on the spatial distribution of multiple projected positions on the high-permeability formation unit along the strike of the formation. Identify formation segments with continuous distribution of projected positions along the strike direction as active fluid permeable formation segments. Continuity is determined by calculating the strike distance between adjacent projected positions. When the strike distance between adjacent projected positions is less than the continuity determination distance, it is determined to be a continuous distribution.

[0118] S65: Extract the characteristic reaction distances of the permeability-associated samples of each formation in the active fluid-permeable formation section. Perform spatial interpolation on the characteristic reaction distances of the permeability-associated samples of each formation in the active fluid-permeable formation section. The spatial interpolation generates a spatial distribution field of characteristic reaction distances within the active fluid-permeable formation section. Perform statistical analysis on the values ​​of the spatial distribution field of characteristic reaction distances. Calculate the upper quartile of the values ​​of the spatial distribution field of characteristic reaction distances as the threshold for determining active reaction. Identify the spatial regions in the spatial distribution field of characteristic reaction distances whose values ​​exceed the threshold for determining active reaction as the active interface reaction regions.

[0119] S66: The active interfacial reaction region indicates the spatial location where the interfacial reaction intensity between the fluid and mica is high. The high interfacial reaction intensity leads to high transfer efficiency of rubidium from the fluid to the mica. The spatial prediction domain of fluid permeation is established along the spatial distribution of the active interfacial reaction region.

[0120] S67: The lateral extension range of the fluid permeability prediction domain extends from the active fluid-permeable stratum segment to both sides in the direction perpendicular to the strike of the high-permeability stratum unit. The lateral extension range is determined based on the spatial width of the active interface reaction area in the direction perpendicular to the strike of the stratum. The longitudinal extension range extends outward from both ends of the active fluid-permeable stratum segment in the direction along the strike of the high-permeability stratum unit. The extension length of the longitudinal extension range is determined based on the length of the active fluid-permeable stratum segment. The dip extension range extends from the active fluid-permeable stratum segment to the hanging wall and footwall in the direction along the dip of the high-permeability stratum unit. The dip extension range is determined based on the average distortion area.

[0121] S68: The lateral extension range, longitudinal extension range and directional extension range of the fluid permeation spatial prediction domain are combined in three-dimensional space to form a second type of spatial envelope, which defines the boundary of the fluid permeation spatial prediction domain in three-dimensional space.

[0122] S69: Stratification is performed within the fluid penetration space prediction domain defined by the second type of spatial envelope. In response to the second type of rubidium occurrence layer being an interlayer interface reaction enrichment mode, the mica interlayer position is identified as a rubidium preferential enrichment layer within the fluid penetration space prediction domain. In response to the second type of rubidium occurrence layer being a whole lattice interface reaction enrichment mode, the entire mica is identified as a rubidium enrichment position within the fluid penetration space prediction domain, thus delineating the fluid penetration space prediction domain as the second type of mineralization region.

[0123] Through the above step S6, the present invention couples the fluid inclusion interface reaction characteristics at the sample scale with the high-permeability stratigraphic unit at the regional scale to establish a fluid permeability spatial prediction domain and delineate the second type of mineralization region.

[0124] Example 2, an embodiment of the present invention, provides a method for analyzing the occurrence state and genetic origin of rubidium in the Shimen rubidium deposit based on the previous embodiment, including:

[0125] In S14, the specific method for estimating the cooling rate of the ore-forming fluid is as follows: the homogenization temperature of fluid inclusions in multiple rubidium-bearing mica samples within the mining area is measured, and the ratio of the homogenization temperature difference between adjacent sample collection locations to the spatial distance is calculated as the temperature spatial gradient. The cooling rate of the ore-forming fluid is proportional to the temperature spatial gradient, and the proportionality coefficient is determined based on the estimated fluid transport velocity and the rock thermal diffusivity.

[0126] In S15, the effective rubidium mole fraction between layers is calculated as follows: when the rubidium mole fraction of the rubidium-containing mica sample exceeds the critical rubidium mole fraction, the effective rubidium mole fraction between layers is equal to the rubidium mole fraction of the rubidium-containing mica sample multiplied by the interlayer position occupancy ratio; when the rubidium mole fraction of the rubidium-containing mica sample does not exceed the critical rubidium mole fraction, the effective rubidium mole fraction between layers is equal to the rubidium mole fraction of the rubidium-containing mica sample.

[0127] In S17, the implementation steps of the non-uniform lattice correction method are as follows: divide the second interplanar spacing by 2 to obtain the single-period interplanar spacing corresponding to the second diffraction order, divide the third interplanar spacing by 3 to obtain the single-period interplanar spacing corresponding to the third diffraction order, and arithmetically average the first interplanar spacing, the single-period interplanar spacing corresponding to the second diffraction order, and the single-period interplanar spacing corresponding to the third diffraction order to obtain the measured structural parameters.

[0128] The lattice period deviation threshold is determined as follows: calculate the standard deviation of the single-period lattice interplanar spacing corresponding to the first lattice interplanar spacing, the second diffraction order, and the third diffraction order; calculate the average of the three lattice interplanar spacings; divide the standard deviation by the average to obtain the relative standard deviation; and use the relative standard deviation as the lattice period deviation threshold.

[0129] In S18, the normalized structural deviation index is calculated as follows: the normalized structural deviation index equals the measured structural parameter minus the theoretical structural parameter, and then divided by the theoretical structural parameter.

[0130] The method for determining the boundary of structural anomaly is as follows: calculate the average and standard deviation of the normalized structural deviation index of all rubidium-containing mica samples in the batch. The upper limit of the boundary of structural anomaly is equal to the average plus twice the standard deviation, and the lower limit of the boundary of structural anomaly is equal to the average minus twice the standard deviation.

[0131] In S25, the radial mode divergence threshold is determined as follows: For rubidium-bearing mica samples collected within the mining area, the absolute value of the divergence space integral is calculated according to steps S21 to S24. The distribution of the absolute value of the divergence space integral is statistically analyzed, and a bimodal distribution is identified. The two peaks of the bimodal distribution correspond to the radial divergence mode sample and the parallel flow mode sample, respectively. The appearance of the bimodal distribution reflects the existence of two different dynamic processes in the deposit. The rapid fluid pulse event of explosive mineralization generates a radial strain field centered on the point source in the sample, and the absolute value of the divergence space integral of the strain field is large. The slow fluid infiltration process of permeable mineralization generates a parallel strain field along the permeation channel in the sample, and the absolute value of the divergence space integral of the strain field is small. The two types of samples form two separate peaks in the distribution of the absolute value of the divergence space integral. The absolute value of the divergence space integral corresponding to the valley position between the two peaks is calculated as the radial mode divergence threshold. This threshold realizes the quantitative discrimination boundary for the two types of dynamic processes.

[0132] In S36, the spatial lag threshold is determined as follows: At least five mica samples with known interlayer occupancy and five mica samples with known occupancy at octahedral positions are identified as calibration samples using crystal structure refinement methods. Mica samples with known interlayer occupancy are those where the rubidium atom electron density is concentrated at interlayer positions as shown in the single-crystal structure refinement results. Mica samples with known octahedral occupancy are those where the rubidium atom electron density is concentrated at octahedral positions as shown in the single-crystal structure refinement results. The median radial offset of the peak position for each calibration sample is calculated according to steps S31 to S36. The interlayer occupancy calibration sample... The spatial lag caused by rubidium atom diffusion in the sample results in a median distribution of the radial offset of the peak position within a larger range. In contrast, the octahedral position occupancy calibration sample, due to the lack of diffusion caused by rubidium atom confinement, results in a median distribution of the radial offset of the peak position within a smaller range. The median distributions of the radial offsets of the peak positions of the two types of calibration samples form two separate ranges. The average of the lower limit of the median distribution of the radial offset of the peak positions of the interlayer position occupancy calibration sample and the upper limit of the median distribution of the radial offset of the peak positions of the octahedral position occupancy calibration sample is calculated as the spatial lag determination threshold. This threshold enables the discrimination of the rubidium occupancy mode of unknown samples.

[0133] In S43, the clustering implementation method is as follows: calculate the distance between adjacent projected positions along the fracture direction on the fracture structure, statistically analyze the distance distribution and calculate the median as the distance threshold. When the distance between adjacent projected positions along the fracture direction is less than the distance threshold, they are determined to be in the same cluster. The fracture segment corresponding to the cluster containing no less than 3 projected positions is identified as the active fluid transport segment of the fracture structure.

[0134] The method for determining the structural consistency angle is as follows: collect more than 10 rubidium-bearing mica samples near the fault structure within the mining area, calculate the dominant propagation direction of rubidium enrichment for each sample according to the method in steps S31 to S38, calculate the angle between the dominant propagation direction of rubidium enrichment for each sample and the strike of the nearest fault structure, statistically analyze the angle distribution and calculate the upper quartile, and use the upper quartile as the structural consistency angle determination value.

[0135] In S44, the longitudinal extension range is determined as follows: statistically analyze the positional distribution of fracture-associated samples along the fracture direction, calculate the distribution density of fracture-associated samples on the active fluid migration segment of the fracture structure, statistically analyze the spatial gradient of sample distribution density along the fracture direction from both ends of the active fluid migration segment of the fracture structure, identify the location where the spatial gradient of sample distribution density changes abruptly as the transition boundary from dense to sparse sample distribution, and use the distance from both ends of the active fluid migration segment of the fracture structure to the transition boundary as the longitudinal extension range.

[0136] The method for determining the dip extension range is as follows: obtain the average width and dip angle of the active fluid migration segment of the fracture structure, calculate the average width divided by the cosine of the dip angle to obtain the projected length of the fracture surface on the dip, and use this projected length as the dip extension range.

[0137] In S55, the method for determining the interface reaction threshold is as follows: calculate the proportion of interface reaction inclusions for the second strain mode samples collected within the mining area according to the methods in steps S51 to S55, statistically analyze the distribution of the proportion of interface reaction inclusions, and calculate the median of the distribution of the proportion of interface reaction inclusions as the interface reaction threshold.

[0138] In S59 and S63, the formation permeability consistency angle threshold and formation permeability consistency angle judgment value are determined as follows: For more than 10 rubidium-bearing mica samples collected near the stratigraphic unit within the mining area, the regional rubidium migration dominant orientation of each sample is calculated according to the method in step S57. The angle between the regional rubidium migration dominant orientation of each sample and the strike of the nearest high-permeability stratigraphic unit is calculated. The angle distribution is statistically analyzed and the upper quartile is calculated. The upper quartile is used as the formation permeability consistency angle threshold and formation permeability consistency angle judgment value.

[0139] In S62, the lower permeability threshold is determined as follows: permeability tests are conducted on different lithological units within the mining area to identify ore-bearing and non-ore-bearing strata units. The permeability distribution of ore-bearing strata units is statistically analyzed and the lower quartile is calculated. This lower quartile is used as the lower permeability threshold.

[0140] In S64, the continuity determination distance is determined by: calculating the strike distance between adjacent projected positions on a high-permeability stratigraphic unit, statistically analyzing the strike distance distribution and calculating the upper quartile, which is then used as the continuity determination distance.

[0141] In S65, the upper quartile is used as the basis for determining the active reaction: the characteristic reaction distance is spatially non-uniformly distributed in the active fluid-permeable strata. The characteristic reaction distance values ​​of the fluid permeable channel development locations differ from those of the non-developed fluid permeable channel locations. The upper quartile divides the spatial distribution field of the characteristic reaction distance into the top 25% of the numerical ranking region and the bottom 75% of the numerical ranking region. The top 25% of the numerical ranking region corresponds to the fluid permeable channel development locations.

[0142] In S67, the longitudinal extension range is determined as follows: the distribution of the projected positions of the formation permeability associated samples along the formation strike is statistically analyzed. The distances along the strike between adjacent projected positions are statistically analyzed from both ends of the active fluid permeable formation segment. The position where the distance along the strike first exceeds the continuity determination distance is identified as the boundary of the continuity interruption of the projected position. The distance from both ends of the active fluid permeable formation segment to the continuity interruption boundary is taken as the longitudinal extension range.

[0143] The method for determining the tendency extension range is as follows: calculate the square root of the average distortion area to obtain the characteristic length of the lattice distortion region, obtain the formation permeability through formation permeability testing and obtain the formation porosity through pore structure analysis, calculate the square root of the formation permeability and divide it by the formation porosity to obtain the formation pore connectivity coefficient, and multiply the characteristic length of the lattice distortion region by the formation pore connectivity coefficient to obtain the tendency extension range.

[0144] Example 3 is an embodiment of the present invention, which provides a method for analyzing the occurrence state of rubidium and the genesis of the Shimen rubidium deposit. In order to verify the beneficial effects of the present invention, scientific demonstration is carried out through experiments.

[0145] Twenty-four rubidium-bearing mica samples from the mineralization zone of the Shimen rubidium deposit were selected for analysis. The sample collection locations covered the known ore body and its periphery. The methods of this invention and existing technologies were used to conduct comparative experiments on rubidium occurrence stratigraphic determination, mineralization mechanism identification, and mineralization area delineation.

[0146] Existing techniques use electron probe microanalysis to determine the rubidium content of samples. For samples with high rubidium content, single-crystal structure refinement is attempted to determine rubidium occupancy. Of the 24 samples tested, only 3 yielded single-crystal particles suitable for structure refinement; the remaining 21 samples could not be refined due to small grain size or severe lattice distortion. X-ray diffraction was used to determine the mineral type of all samples; 18 samples were identified as muscovite, and 6 samples as lepidolite. Homogenization temperature of fluid inclusions was used to determine the mineralization temperature. Kriging interpolation was performed based on known orebody boundaries to delineate the mineralization range.

[0147] The method of this invention calculates theoretical and measured structural parameters for 24 samples. The normalized structural deviation index of 16 samples exceeds the boundary for structural anomaly judgment. Strain spatial distribution analysis of the 16 structurally anomalous samples reveals that 9 samples exhibit radial divergence and 7 exhibit parallel flow. Isotopic composition spatial variation analysis of the 9 radial divergence samples shows that the median radial offset of the peak position in 5 samples exceeds the spatial hysteresis judgment threshold, indicating interlayer occupancy, and 4 samples exhibit octahedral occupancy. Internal observation of fluid inclusions in the 7 parallel flow samples reveals that fluid inclusions with abrupt changes in rubidium concentration gradients in 4 samples are concentrated in the mica interlayer, indicating an interlayer interface reaction enrichment mode, and 3 samples exhibit a full-lattice interface reaction enrichment mode.

[0148] The dominant propagation direction of rubidium enrichment in nine radially divergent model samples was compared with the fault structures in the mining area. The enrichment direction of six samples aligned with the northeast-trending fault, and the enrichment direction of three samples aligned with the northwest-trending fault. Spatial prediction domains for fluid migration were established along the northeast-trending and northwest-trending faults, and the first-type mineralization areas were delineated. The dominant regional rubidium migration orientation of seven parallel flow model samples was compared with the stratigraphic attitude. The migration orientation of five samples aligned with the strike of the sandstone layer, and the migration orientation of two samples aligned with the strike of the shale layer. Spatial prediction domains for fluid permeability were established along the sandstone and shale layers, and the second-type mineralization areas were delineated.

[0149] In the first type of mineralization area delineated by the method of this invention, 12 verification boreholes were drilled, and 10 boreholes encountered mineralization. In the second type of mineralization area, 8 verification boreholes were drilled, and 7 boreholes encountered mineralization. Among the 20 verification boreholes, the ore-bearing strata of 17 boreholes were consistent with the rubidium preferential enrichment strata predicted by the method of this invention.

[0150] Ten historical boreholes with locations similar to the verification boreholes of this invention were selected within the mineralized area delineated by existing methods for comparison; six of these boreholes encountered mineralization. Four boreholes were drilled outside the boundary of the mineralized area delineated by existing methods at the predicted ore body extension location; three of these boreholes encountered mineralization.

[0151] Existing methods determine rubidium occupancy through single-crystal structure refinement, with 3 out of 24 samples successfully identified. The method of this invention determines rubidium occurrence sites through structural deviation and strain analysis, with 16 out of 24 samples successfully identified. Existing methods determine mineralization temperature through fluid inclusion homogenization temperature but cannot distinguish mineralization mechanism types. The method of this invention distinguishes between explosive and infiltrative mineralization through the divergence topological characteristics of strain spatial distribution. The method of this invention verifies a borehole mineralization rate of 17 mineralized boreholes out of 20 total boreholes, while the historical borehole mineralization rate of existing methods is 6 mineralized boreholes out of 10 total boreholes. The method of this invention achieves a higher mineralization rate than existing methods. Furthermore, the method of this invention found mineralization in 3 out of 4 borehole locations predicted outside the mineralization zone boundary delineated by existing methods, confirming that existing methods may miss mineralization extents. The method of this invention enables quantitative determination of the lattice position of rubidium occurrence layers, differentiation of the dynamic process of mineralization mechanism, and delineation of geological structural constraints of mineralized areas, thereby improving the sample coverage and accuracy of mineralization area delineation for determining the occurrence state of rubidium.

[0152] This embodiment also provides an electronic device applicable to the analysis of the occurrence state of rubidium and the genetic analysis method of the Shimen rubidium deposit, comprising: a memory and a processor; the memory is used to store computer-executable instructions, and the processor is used to execute the computer-executable instructions to realize the analysis of the occurrence state of rubidium and the genetic analysis method of the Shimen rubidium deposit as proposed in the above embodiment.

[0153] This embodiment also provides a storage medium storing a computer program that, when executed by a processor, implements the rubidium occurrence state and ore genesis analysis method of the Shimen rubidium deposit as proposed in the above embodiments.

[0154] The storage medium proposed in this embodiment belongs to the same inventive concept as the method for realizing the occurrence state of rubidium and the genetic analysis of the rubidium deposit in the Shimen rubidium deposit proposed in the above embodiments. Technical details not described in detail in this embodiment can be found in the above embodiments, and this embodiment has the same beneficial effects as the above embodiments.

[0155] Based on the above description of the implementation methods, those skilled in the art can clearly understand that the present invention can be implemented using software and necessary general-purpose hardware, and of course, it can also be implemented using hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as a computer floppy disk, read-only memory (ROM), random access memory (RAM), flash memory, hard disk, or optical disk, etc., including several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods of the various embodiments of the present invention.

[0156] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A method for analyzing the occurrence state of rubidium and the genesis of a Shimen rubidium deposit, characterized in that, include: Obtain rubidium-containing mica samples, calculate theoretical structural parameters based on the elemental composition of the rubidium-containing mica samples, obtain measured structural parameters through diffraction testing, and identify structurally abnormal samples based on the deviation between theoretical and measured structural parameters. The spatial distribution characteristics of strain are obtained by analyzing the diffraction signal morphology of structurally anomalous samples. Based on the dominant direction of strain spatial distribution, the structurally anomalous samples are divided into first strain mode samples and second strain mode samples. The isotopic composition of the first strain mode sample was spatially sampled and determined along the dominant direction. Based on the spatial variation characteristics of the isotopic composition, the first type of rubidium occurrence layer and the first type of enrichment mechanism were determined. Based on the type I rubidium occurrence layer and enrichment mechanism, the spatial distribution of the samples was correlated with the direction of the tectonic stress field to delineate the type I mineralization region; The directional characteristics of fluid inclusions in the second strain mode sample were observed and obtained. Based on the matching of the directional characteristics with the anisotropy of the formation, the second type of rubidium occurrence layer and the second type of enrichment mechanism were determined. Based on the second type of rubidium occurrence horizon and enrichment mechanism, the spatial distribution of the samples was correlated with the formation permeability direction to delineate the second type of mineralization region.

2. The method for analyzing the occurrence state of rubidium and the genesis of a rubidium deposit in a stone gate rubidium deposit according to claim 1, characterized by, The calculation of theoretical structural parameters based on the elemental composition of the rubidium-containing mica sample includes: The rubidium and potassium content of the rubidium-containing mica sample was determined. A database of rubidium-potassium mica standard samples in equilibrium was established, and the relationship curve between rubidium mole fraction and interplanar spacing was extracted. A dataset of experimental relationships between cooling rate and interlayer position occupancy ratio was established through high temperature and high pressure tests. The cooling rate of ore-forming fluid was estimated based on the uniform temperature spatial variation gradient of fluid inclusions. The theoretical structural parameters are determined based on the rubidium content, potassium content, interlayer site occupancy ratio, and relationship curves.

3. The method for analyzing the occurrence state of rubidium and the genesis of a rubidium deposit in a stone gate rubidium deposit according to claim 2, characterized by, The method of obtaining the measured structural parameters through diffraction testing includes: X-ray diffraction tests were performed on rubidium-containing mica samples to identify the first, second, and third diffraction order secondary peaks of the mica's layered structure. The deviation of the interplanar spacing corresponding to the three diffraction order sub-peaks is used to identify lattice periodic inhomogeneity. In response to the presence of lattice periodic inhomogeneity, the measured structural parameters are determined according to the inhomogeneous lattice correction method. In response to the absence of lattice periodic inhomogeneity, the measured structural parameters are determined from the first diffraction order sub-peak.

4. The method for analyzing the occurrence state of rubidium and the genesis of a rubidium deposit in a stone gate according to claim 3, characterized in that, The method of classifying structurally abnormal samples into first strain mode samples and second strain mode samples based on the dominant direction of strain spatial distribution includes: Spatially distributed thin sections were prepared for samples with structural anomalies and observed using a polarizing microscope. The interference color order of each observation point in the observation point array was recorded and converted into birefringence values. Determine the dominant gradient direction and intensity at each observation point, and construct a gradient vector field distribution map; The gradient vector field distribution map is integrated in the divergence space. If the absolute value of the divergence space integral is greater than the radial mode divergence threshold, the structurally anomalous sample is identified as a first strain mode sample. If the absolute value of the divergence space integral is less than or equal to the radial mode divergence threshold, the structurally anomalous sample is identified as a second strain mode sample.

5. The method for analyzing the occurrence state and genetic origin of rubidium in the Shimen rubidium deposit as described in claim 4, characterized in that, The determination of the first type of rubidium occurrence sites and the first type of enrichment mechanism based on the spatial variation characteristics of isotopic composition includes: A polar coordinate system was established for the first strain mode sample and multiple radiation lines were set. The rubidium 87 abundance and rubidium 85 abundance at each measurement position on each radiation line were measured. The locations of strain peaks and isotope fractionation peaks were identified for each radiation survey line, and the first type of rubidium-bearing horizons were determined based on the statistical distribution of the radial offset of the peak locations. The azimuth distribution characteristics of isotope attenuation amplitude were extracted from each radiation line, and the dominant propagation direction of rubidium enrichment and the distribution range of isotope attenuation amplitude were determined as the first type of enrichment mechanism.

6. The method for analyzing the occurrence state and genetic origin of rubidium in the Shimen rubidium deposit as described in claim 5, characterized in that, The process of delineating the first type of mineralization region by correlating the spatial distribution of samples with the direction of tectonic stress field based on the first type of rubidium occurrence horizon and enrichment mechanism includes: Identify samples with the first strain mode, compare the rubidium enrichment-dominant propagation direction with the fracture structure orientation to identify fracture-related samples, project the fracture-related samples onto the fracture structure and cluster them to identify the active fluid transport segments of the fracture structure; A spatial prediction domain for fluid migration is established along the active fluid migration segment of the fracture structure. The three-dimensional extension range of the spatial prediction domain for fluid migration is determined based on the distribution range of isotope attenuation amplitude, the geometric parameters of the active fluid migration segment of the fracture structure, and the dip angle of the fracture structure, and a first-type spatial envelope is formed. Within the fluid transport space prediction domain defined by the first type of spatial envelope, rubidium preferential enrichment layers are identified and first type mineralization regions are delineated based on the first type of rubidium occurrence layers.

7. The method for analyzing the occurrence state and genetic origin of rubidium in the Shimen rubidium deposit as described in claim 6, characterized in that, The determination of the second type of rubidium-bearing strata and the second type of enrichment mechanism based on the matching of directional characteristics and stratigraphic anisotropy includes: Transmission electron microscopy was used to identify the lattice distortion region by imaging the mica lattice around the fluid inclusions in the sample of the second strain mode, and the rubidium concentration distribution within the radius of influence of the fluid inclusions was determined and the location of the abrupt change in the rubidium concentration gradient was identified. The distribution pattern of fluid inclusions with abrupt changes in rubidium concentration gradients within mica crystals was used to determine the second type of rubidium occurrence sites. The distribution of vector directions pointing from the center of the fluid inclusions to the location of abrupt changes in the rubidium concentration gradient is used to obtain the dominant orientation of rubidium migration in the region, and the spatial distribution characteristics of the locations of abrupt changes in the rubidium concentration gradient and the geometric characteristics of the lattice distortion region are extracted. In response to the matching relationship between the dominant rubidium migration orientation and the stratigraphic strike satisfying the stratigraphic permeability consistency condition, the dominant rubidium migration orientation and fluid inclusion interface reaction characteristic parameters are identified as the second type of enrichment mechanism.

8. The method for analyzing the occurrence state and genetic origin of rubidium in the Shimen rubidium deposit as described in claim 7, characterized in that, The method of delineating the second type of mineralization region by correlating the spatial distribution of samples with the direction of formation permeability based on the second type of rubidium occurrence horizon and enrichment mechanism includes: Identify samples with the second strain mode, compare the dominant orientation of regional rubidium migration with the strike of the high-permeability stratigraphic unit to identify stratigraphic permeability-related samples, project the stratigraphic permeability-related samples onto the high-permeability stratigraphic unit and continuously analyze to identify active fluid permeable stratigraphic segments; Spatial interpolation is performed on the fluid inclusion interface reaction characteristic parameters of permeability-associated samples in active fluid-permeable strata to identify active interface reaction regions and establish a spatial prediction domain for fluid permeability. Based on the spatial geometric characteristics of the active interface reaction region and the geometric parameters of the active fluid-permeable stratum, the three-dimensional extension range of the fluid permeability spatial prediction domain is determined and a second type of spatial envelope is formed. Within the fluid permeation space prediction domain defined by the second type of spatial envelope, rubidium preferential enrichment layers are identified and second type mineralization regions are delineated based on the second type of rubidium occurrence layers.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of analyzing the occurrence state of rubidium in the Shimen rubidium deposit and the genetic analysis method of the deposit as described in any one of claims 1 to 8.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for analyzing the occurrence state of rubidium in the Shimen rubidium deposit and the genetic analysis of the deposit as described in any one of claims 1 to 8.

Images