Evaluation method for element enrichment and loss in geological structure based on anisotropic singularity

Abstract

The invention discloses a geologic structure element enrichment and loss evaluation method based on anisotropic singularity, and the method comprises the steps: building a space coupling relation between a mineral product enrichment degree and an element singularity index in a linear structure direction through calculating the element singularity index in each direction of the linear structure; and then evaluating the metallogenic enrichment degree based on the anisotropic singularity calculation model of the exploration geochemical region to obtain an evaluation report of the metallogenic enrichment degree in each direction of the linear structure.

Description

Technical Field

[0001] The present invention relates to the field of geology, and particularly to a method for evaluating element enrichment and depletion in geological structures based on anisotropy. Background Art

[0002] The calculation of regional anisotropy in exploration geochemistry is an important mathematical method for geologists to process geological data and find out the distribution law of mineral resources. The existing methods calculate point by point within a neighborhood or in several directions within a neighborhood, often without calculating the anisotropy in all directions of all sampling points and without fully considering the coupling relationship with linear structures, which will have an adverse impact on the evaluation of the final ore-forming enrichment degree. Summary of the Invention

[0003] In view of the above deficiencies in the prior art, the present invention provides a method for evaluating element enrichment and depletion in geological structures based on anisotropy.

[0004] In order to achieve the above invention object, the technical solution adopted by the present invention is as follows:

[0005] A method for evaluating element enrichment and depletion in geological structures based on anisotropy, comprising the following steps:

[0006] S1. Obtain the spatial distribution data of metal minerals and the spatial distribution data of geological linear structures in the study area;

[0007] S2. Calculate the anisotropy index of geochemical elements and the direction of geological linear structures based on the data obtained in S1;

[0008] S3. Establish a spatial coupling relationship between the ore enrichment degree and the anisotropy index in the direction of the linear structure according to the anisotropy index calculated in S2 and the direction of the geological linear structure;

[0009] S4. Evaluate the ore-forming enrichment degree based on the calculation model of regional anisotropy in exploration geochemistry to realize the evaluation of the direction of the linear structure and the ore-forming enrichment degree.

[0010] Further, the specific calculation method of the direction of the geological linear structure in S2 is as follows:

[0011] S201. Set calculation parameters, including step size, number of steps, and Euclidean space dimension;

[0012] S202. Calculate the direction of the geological linear structure by using the start-end connection method or the segmented calculation method.

[0013] Further, the specific calculation method of the start-end connection method is as follows:

[0014] A1. Obtain the starting point and ending point of the linear structure line data, and connect them into a straight line;

[0015] A2. Define the due east direction as 0°, the due north direction as 90°, and calculate the mathematical expression of the start-stop connection line by rotating counterclockwise;

[0016] A3. Take the counterclockwise angle between the geological linear structure line and the due east direction as the direction of the geological linear structure line.

[0017] Furthermore, the specific calculation method of the segmented calculation method is as follows:

[0018] B1. Obtain the non-smooth curve of the linear structure line data, and starting from its first point, find its adjacent point as the second point and connect them into a straight line as the start-stop connection line of this segment;

[0019] B2. Define the due east direction as 0°, the due north direction as 90°, and calculate the mathematical expression of the start-stop connection line in B1 by rotating counterclockwise. Take the counterclockwise angle between this straight line and the due east direction as the direction of this straight line segment;

[0020] B3. Repeat steps B1 - B2 to calculate the direction of all linear structure lines in the area segment by segment. Group and sort the directions according to the occurrence frequency, using the direction value and direction interval. The angle with the highest occurrence frequency is the direction of the calculated geological linear structure line.

[0021] Furthermore, the specific calculation method of the anisotropy index of each geochemical element in S2 is as follows:

[0022] S211. After selecting the direction of the geological linear structure line, set the anisotropy index calculation parameters, including step size, number of steps, Euclidean space dimension, and calculation element;

[0023] S212. Starting from the selected direction, draw a rectangular area according to the step size and number of steps, and determine whether there are exploration geochemical data points in the rectangular area. If so, calculate the total content of the elements in the area based on the total content data of the selected elements;

[0024] S213. Taking the current point of the selected geochemical exploration data point as the center, calculate the total amount of elements and the corresponding area at different distance scales in the selected direction to obtain the element content density; perform a linear fit after double logarithmic transformation with the element content density and different distance scales as the bivariate variables, and the obtained exponent of the linear fit is the anisotropy index in the specified direction;

[0025] S214. Execute the above steps S211 - S213 for all points in the study area to obtain the anisotropy of all points in the study area in the specified direction.

[0026] Further, the specific calculation method of the anisotropy index of geochemical elements in S213 is as follows:

[0027] logC[A(εi)] = c + (α - E)logεi

[0028] In the formula, C[A(εi)] is the total content of elements in the calculation area, εi is the total content data of the selected element; E is the dimension of Euclidean space; α is the anisotropy index of geochemical elements, and c is a constant.

[0029] The present invention has the following beneficial effects:

[0030] After simple parameter settings, the present invention can obtain the calculation result within 1 minute, greatly improving the calculation efficiency. After calculating in the direction of introducing linear structure data, multiple directions can be more accurately selected for batch calculation. The methods for calculating linear structure data are further divided into the start-stop connection method and the segmented calculation method. The calculation results of different methods in each direction can also be compared. After calculation, a spatial coupling relationship between the degree of mineral enrichment and the anisotropy index of elements in the linear structure direction is established, realizing an evaluation method for calculating the anisotropy index model of the linear structure direction and the degree of ore-forming enrichment. The linear structure direction calculated by the software is grouped and sorted by statistical methods and intuitively displayed to geological researchers. After selecting one or more directions, the anisotropy of all points in the entire area in the selected direction is calculated, covering the calculation in multiple directions without increasing the workload, greatly improving the calculation coverage rate and reducing the risk of omission; Description of the Drawings

[0031] Figure 1 It is a schematic flow chart of the method for evaluating element enrichment and depletion in geological structures based on anisotropy of the present invention.

[0032] Figure 2 It is the spatial coupling relationship between the degree of mineral enrichment and the anisotropy index of elements in the linear structure direction in the embodiment of the present invention. Detailed Embodiments

[0033] The following describes the specific embodiments of the present invention to facilitate those skilled in the art of the present technology to understand the present invention. However, it should be clear that the present invention is not limited to the scope of the specific embodiments. For those of ordinary skill in the art of the present technology, as long as various changes are within the spirit and scope of the present invention defined and determined by the appended claims, these changes are obvious, and all inventions made using the concept of the present invention are within the scope of protection.

[0034] A method for evaluating element enrichment and depletion in geological structures based on anisotropy, as Figure 1 shown, includes the following steps:

[0035] S1. Obtain the spatial distribution data of metal minerals and the spatial distribution data of geological linear structures within the research area;

[0036] Calibrate the gold mine location on the GIS software according to the data, and store the exploration geochemical data of the gold mine location in the shp file.

[0037] The description of the shp file is as follows:

[0038] Shapefile is a vector graphics format that can store the location of geometric figures and related attributes.

[0039] The main file of the Shapefile format contains georeferenced data. This file consists of a fixed-length file header and one or several variable-length record data. Each variable-length data record contains a record header and some record content.

[0040] The main file header contains 17 fields, a total of 100 bytes, including nine 4-byte (32-bit signed integer, int32) integer fields, followed by eight 8-byte (double-precision floating-point number) signed floating-point fields.

[0041] Then this file contains an indefinite number of variable-length data records. Each data record starts with an 8-byte record header, and the actual record follows the record header.

[0042] The content of the variable-length record is determined by the type of the figure. Shapefile supports multiple figure types, and two data types, points and polylines, will be used in the present invention.

[0043] When extracting the spatial data of the gold mine, the path needs to be specified, that is, the location of the shp file of the exploration geochemical points, and the figure type is point.

[0044] Shapefile usually uses X and Y to represent geographical coordinates, where X corresponds to longitude and Y corresponds to latitude. Relevant information required for subsequent calculations, especially spatial coordinate information, represented by X and Y, can be obtained by parsing the file.

[0045] The spatial coordinate system is the spatial coordinate system used for the geographical coordinates stored in the shp file. The purpose of obtaining the spatial coordinate system is to accurately perform spatial transformation and position calculation during visual data presentation.

[0046] The plane coordinate system is used for subsequent calculations, and the WGS84 spatial coordinate system is used for visual presentation. This coordinate system uses longitude and latitude to represent the position of points. Therefore, the process of calculation and visualization involves two transformations. The first is to convert the source data into plane coordinates, and the second is to convert the plane coordinates into WGS84 coordinates.

[0047] Draw the linear structure on the GIS software according to the data, and store the point set data included in the linear structure in the shp file. When extracting the spatial data of the linear structure, the path needs to be specified, that is, the location of the shp file of the linear structure, and the graphic type is a polyline.

[0048] S2. Calculate the anisotropy index of each geochemical element and the direction of the geological linear structure based on the data obtained in S1;

[0049] In this embodiment, the specific calculation method of the geological linear structure direction is as follows:

[0050] S201. Set the calculation parameters, including the step size, the number of steps, and the Euclidean space dimension;

[0051] S202. Calculate the direction of the geological linear structure by using the start-end connection method or the segmented calculation method. Among them,

[0052] Description of the linear structure direction calculation method - start-end connection method:

[0053] The spatial data of the linear structure line is expressed as a non-smooth curve, which is represented as a point set in the data structure (see the description of the shp file).

[0054] The start-end connection method finds the start point and the end point of the linear structure line data, that is, the first point and the last point of the point set, and connects them into a straight line;

[0055] Calculate the slope of this straight line: Define the due east direction as 0 degrees, and the mathematical expression is y = kx + b, where k = 0 and the slope is 0. The due north is 90 degrees, and the mathematical expression is x = b. The mathematical expression of the start-end connection line can be calculated by rotating counterclockwise. According to the corresponding expression, the direction of this straight line can be calculated; the value range of the direction is between 0 and 360 degrees, and the direction depends on the positions of the start and end points of the linear structure.

[0056] The direction of the linear structure line is the counterclockwise angle with the due east direction after connecting the start and end points.

[0057] Description of the linear structure direction calculation method - segmented calculation method:

[0058] The spatial data of the linear structure line is expressed as a non-smooth curve, which is represented as a point set in the data structure (see the description of the shp file).

[0059] Starting from the first point of the point set, that is, the start point, find its adjacent point, that is, the second point of the point set, and connect them into a straight line;

[0060] Calculate the slope of the straight line: Define the due east direction as 0 degrees, and the mathematical expression is y = kx + b, where k = 0 and the slope is 0. The due north direction is 90 degrees, and the mathematical expression is x = b. By calculating the rotation in the counterclockwise direction, the mathematical expression of the start-stop connection line can be obtained, and the direction of the straight line can be calculated based on the corresponding expression. The value range of the direction is 0 - 360 degrees, and the direction depends on the positions of the start and end points of the linear structure.

[0061] The direction of the linear structure line is the counterclockwise angle between the connection of the start and end points and the due east direction.

[0062] Using the above method, calculate the direction of each segment of all the linear structure lines in the area one by one, and group and sort the directions according to the frequency of occurrence, using the direction values and direction intervals (such as 0 - 30 degrees, etc.).

[0063] The direction angle field is selected according to the angle with the highest frequency of occurrence after calculating the direction of the linear structure line.

[0064] The calculation and statistical results of the direction of the linear structure are the pre - steps for calculating the singularity. The purpose is to select the direction. The specific calculation method of the anisotropy index of each geochemical element is as follows:

[0065] S211. After selecting the direction of the geological linear structure line, set the singularity index calculation parameters, including the step size, number of steps, Euclidean space dimension, and calculation element. Among them,

[0066] Step size: The side length of the rectangular box when calculating the singularity

[0067] Number of steps: The number of rectangular boxes when calculating the singularity

[0068] Euclidean space dimension: When calculating on a plane, the Euclidean space dimension is 2

[0069] Element: Gold, silver, copper, etc. The data is from the attribute words of the specified shp file.

[0070] S212. Starting from the selected direction, draw rectangular areas according to the step size and number of steps, and judge whether the rectangular areas contain the detection points of exploration geochemistry. If so, calculate the total content of the elements in the area based on the total content data of the selected elements;

[0071] The initial rectangular box faces the due east direction. First, calculate the positions of the vertices of the initial rectangular box, and then rotate counterclockwise according to the selected angle to obtain the positions of the vertices of the rotated rectangular box. Given that the plane coordinates of the starting point are (x, y), the initial vertex coordinate calculation formula of the starting rectangular box is:

[0072]

[0073] wherein, lt is the left-top vertex of the rectangular frame, rt is the right-top vertex of the rectangular frame, lb is the left-bottom vertex of the rectangular frame, rb is the right-bottom vertex of the rectangular frame, and the calculation formulas for the initial vertex coordinates of the remaining rectangular frames are as follows:

[0074]

[0075]

[0076] The calculated and selected angle is angle, r is the angle value used when substituting into the formula, w is the width of the rectangular frame when calculating singularity, h is the height of the rectangular frame, the four vertices of the rectangular frame constructed according to the step size are A, B, C, D, and their coordinates on the plane are newAX, newAY, and so on; the center point coordinates of the rectangular frame are described as centerX and CenterY.

[0077]

[0078] w = dis + dis;

[0079] h = dis;

[0080] newAX = centerX + (x - centerX) * cos(r) - (y - centerY) * sin(r);

[0081] newAY = centerY + (x - centerY) * sin(r) + (y - centerY) * cos(r);

[0082] newBX = centerX + (x + w - centerX) * cos(r) - (y - centerY) * sin(r);

[0083] newBY = centerY + (x + w - centerX) * sin(r) + (y - centerY) * cos(r);

[0084] newCX = centerX + (x + w - centerX) * cos(r) - (y + h - centerY) * sin(r);

[0085] newCY = centerY + (x + w + centerX) * sin(r) + (y + h - centerY) * cos(r);

[0086] newDX = centerX + (x - centerX) * cos(r) - (y + h - centerY) * sin(r);

[0087] newDY = centerY + (x - centerX) * sin(r) + (y + h - centerY) * cos(r);

[0088] The coordinates of the 4 pairs of vertices after rotation are obtained through the above formula.

[0089] S213. Calculate the corresponding data between the total content data of a set of selected elements and the content of the selected elements at the current point of the geochemical exploration detection point on the earth, and perform linear fitting on this corresponding data to obtain the anisotropy index of each geochemical element;

[0090] Determine whether the detection points of exploration geochemistry are included in the rectangular box (the determination method is not described here).

[0091] Then calculate the singularity according to the following formula

[0092] logC[A(εi)] = c + (α - E)logεi

[0093] For the total content data of the selected elements of the geochemical exploration detection points in the rectangular box, calculate the total content of the regional elements, that is, C[A(εi)];

[0094] The content data of the selected elements of the geochemical exploration detection point at the current point is εi;

[0095] The dimension of the Euclidean space, that is, E;

[0096] c is a constant;

[0097] A set of corresponding data of εi and C[A(εi)] can be calculated, and linear fitting on this set of data can obtain α, that is, the singularity.

[0098] S214. Execute the above steps S211 - S213 for all points in the study area to obtain the singularities of all points in the study area in the specified direction.

[0099] Execute the above calculation process for all points in the area, and the singularities of all points in the area in the specified direction can be obtained.

[0100] S3. Establish a spatial coupling relationship between the degree of mineral enrichment and the singularity index in the direction of the linear structure according to the singularity index calculated in S2 and the direction of the geological linear structure;

[0101] Perform mean division on the calculation results, bind the spatial positions of the sampling points, and mark them with different colors on the map, as Figure 2 shown, to obtain the spatial coupling model of the singularity index.

[0102] The model is described as follows: According to the selected direction, calculate the singularity of sampling points in the whole area in the selected direction, and classify the singularity values of all points into three levels by mean value, namely high, medium and low, which are represented by red, yellow and blue respectively. The position of the colored points represents the singularity in the linear structure direction at the spatial position of the sampling points, that is, the spatial coupling of the singularity index, as shown in Table 1.

[0103] Table 1 Representation of singularity in the linear structure direction at the spatial position

[0104] Singularity value representation Singularity segmented interval Color representation Numerical value < 1.1046 High Red 1.1046 ≤ value < 1.1905 Medium Yellow Value ≥ 1.1905 Low Blue

[0105] S4. Evaluate the degree of ore-forming enrichment based on the anisotropy calculation model of the exploration geochemistry region, and obtain an evaluation report on the linear structure direction and the degree of ore-forming enrichment.

[0106] Evaluate the degree of ore-forming enrichment based on the anisotropy calculation model of the exploration geochemistry region, and obtain an evaluation report on the linear structure direction and the degree of ore-forming enrichment. The report is a quantitative evaluation of the degree of ore-forming enrichment, which is reflected in the following aspects:

[0107] Relationship between singularity index and enrichment degree

[0108] Through the calculated singularity index, combined with the geological background and ore-forming theory, quantify the degree of ore-forming enrichment. When the singularity index is in the range of (0-2) and closer to 2, it indicates that the element enrichment degree in this area is higher and the ore-forming potential is greater. The enrichment degree can be divided into different levels (such as high, medium and low) according to the numerical range of the singularity index, so as to more intuitively evaluate the ore-forming potential.

[0109] Spatial distribution of enrichment degree

[0110] Combine the calculation results of the singularity index with the spatial distribution map to generate the spatial distribution map of the ore-forming enrichment degree. Through visualization means, geological workers can intuitively identify high-enrichment areas, which may be the key points for further exploration.

[0111] Comprehensive analysis of multi-direction singularity index

[0112] Comparative analysis of multi-direction singularity index: Since the present invention introduces the calculation of the direction of the linear structure, the singularity index can be calculated in multiple directions at the same time. By comparing and analyzing the singularity indices in different directions, it is possible to identify which directions have higher element enrichment degrees, thus providing a more comprehensive basis for ore-forming prediction.

[0113] Coupling relationship between direction and enrichment degree: By establishing the spatial coupling relationship between the linear structure direction and the singularity index, analyze the element enrichment characteristics in different directions. For example, the singularity index in certain directions may be significantly higher than that in other directions, indicating that the geological structures in these directions have a stronger control on ore-forming enrichment.

[0114] Evaluation report on ore-forming enrichment degree

[0115] Visual display of the report: Visualize data such as singularity index, linear structure direction, and ore-forming enrichment degree in the form of maps, charts, etc. The visualization results can help geological workers more intuitively understand the spatial distribution characteristics of ore-forming enrichment degree and provide decision-making support for subsequent exploration work.

[0116] Verification and optimization of ore-forming enrichment degree

[0117] Field verification: According to the high-enrichment areas in the evaluation report, conduct field verification work. Through on-site sampling and laboratory analysis, verify whether the relationship between the singularity index and the ore-forming enrichment degree is accurate. The field verification results can be fed back into the model to further optimize the calculation parameters and evaluation methods of the singularity index.

[0118] Model optimization: According to the field verification results and actual exploration experience, optimize the singularity calculation model. For example, adjust parameters such as step size and number of steps, or introduce more geological data (such as lithology, structure, etc.) to improve the prediction accuracy and reliability of the model.

[0119] Dynamic evaluation and update of ore-forming enrichment degree

[0120] Dynamic update and evaluation: As the exploration work progresses and new data is obtained, the calculation results of the singularity index can be dynamically updated, and the ore-forming enrichment degree can be re-evaluated. Through continuous iterative optimization, the prediction ability of the model will gradually improve, providing a more reliable basis for the exploration and development of mineral resources.

[0121] This invention is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems), and computer program products according to embodiments of the invention. It should be understood that each process and / or block in the flowchart and / or block diagram, and the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate a device for realizing the specified functions in Figure One one process or multiple processes and / or blocks Figure One one block or multiple blocks.

[0122] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to operate in a particular manner, such that the instructions stored in the computer-readable memory produce a manufacture including an instruction means that implements the functions specified in one or more of the processes and / or blocks Figure One of one or more of the processes and / or blocks Figure One specified in the one or more of the processes and / or blocks.

[0123] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus, such that a series of operational steps are performed on the computer or other programmable apparatus to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable apparatus provide steps for implementing the functions specified in one or more of the processes and / or blocks Figure One of one or more of the processes and / or blocks Figure One specified in the one or more of the processes and / or blocks.

[0124] Specific embodiments are used in the present invention to illustrate the principles and implementation manners of the present invention. The description of the above embodiments is only used to help understand the method and its core idea of the present invention; at the same time, for those of ordinary skill in the art, according to the idea of the present invention, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to the present invention.

[0125] Those of ordinary skill in the art will realize that the embodiments described herein are for helping the reader understand the principles of the present invention, and it should be understood that the protection scope of the present invention is not limited to such specific statements and embodiments. Those of ordinary skill in the art can make various other specific deformations and combinations that do not depart from the essence of the present invention based on the technical revelations disclosed in the present invention, and these deformations and combinations are still within the protection scope of the present invention.

Claims

1. A method for evaluating element enrichment and depletion in geological structures based on anisotropy, characterized in that, It includes the following steps: S1. Obtain the spatial distribution data of metallic minerals and the spatial distribution data of geological linear structures within the research area; S2. Calculate the anisotropy index of each geochemical element and the direction of geological linear structures based on the data obtained in S1; S3. Establish the spatial coupling relationship between the degree of mineral enrichment and the anisotropy index in the direction of linear structures according to the anisotropy index and the direction of geological linear structures calculated in S2; S4. Evaluate the degree of ore-forming enrichment based on the regional anisotropy calculation model of exploration geochemistry, and realize the evaluation of the direction of linear structures and the degree of ore-forming enrichment.

2. The method for evaluating element enrichment and depletion in a geological structure based on anisotropy according to claim 1, wherein The specific calculation method of the direction of geological linear structures in S2 is as follows: S201. Set calculation parameters, including step size, number of steps, and Euclidean space dimension; S202. Calculate the direction of geological linear structures by using the start-end connection method or the segmented calculation method.

3. The method for evaluating element enrichment and depletion in a geological structure based on anisotropy according to claim 2, characterized in that, The specific calculation method of the start-end connection method is as follows: A1. Obtain the starting point and ending point of the linear structure line data and connect them into a straight line; A2. Define the due east direction as 0°, the due north direction as 90°, and calculate the mathematical expression of the start-end connection line by rotating counterclockwise; A3. Take the counterclockwise angle between the geological linear structure line and the due east direction as the direction of the geological linear structure line.

4. The method for evaluating element enrichment and depletion in a geological structure based on anisotropy according to claim 2, wherein The specific calculation method of the segmented calculation method is as follows: B1. Obtain the non-smooth curve of the linear structure line data, and starting from its first point, find its adjacent point as the second point and connect them into a straight line as the start-end connection line of this segment; B2. Define the due east direction as 0°, the due north direction as 90°, calculate the mathematical expression of the start-end connection line in B1 by rotating counterclockwise, and take the counterclockwise angle between this straight line and the due east direction as the direction of this straight line segment; B3. Repeat steps B1 - B2 to calculate the direction of all linear structure lines in the area segment by segment, and statistically group and sort the directions according to the frequency of occurrence, using the direction value and direction interval. The angle with the highest frequency of occurrence is the direction of the calculated geological linear structure line.

5. The method for evaluating element enrichment and depletion in a geological structure based on anisotropy according to claim 1, wherein The specific calculation method of the anisotropy index of each geochemical element in S2 is as follows: S211. After selecting the direction of the geological linear structure line, set the anisotropy index calculation parameters, including step size, number of steps, Euclidean space latitude, and calculation element; S212. Starting from the selected direction, draw a rectangular area according to the step size and the number of steps, and judge whether the exploration geochemistry data points are included in the rectangular area. If so, calculate the total content of the elements in the area based on the total content data of the selected elements; S213. Taking the current point of the selected exploration geochemistry data point as the center, calculate the total amount of elements and the corresponding area at different distance scales in the selected direction to obtain the element content density; perform a linear fit after double logarithmic transformation with the element content density and different distance scales as the two variables, and the obtained exponent of the linear fit is the anisotropy index in the specified direction; S214. Execute the above steps S211 - S213 for all points in the research area to obtain the anisotropy of all points in the research area in the specified direction.

6. The method for evaluating element enrichment and depletion in a geological structure based on anisotropy according to claim 1, characterized in that, The specific calculation method of the anisotropy index of each geochemical element in S213 is as follows: logC[A(εi)] = c + (α - E)logεi Where C[A(εi)] is the total content of elements in the calculation area, εi is the total content data of the selected element; E is the dimension of Euclidean space; α is the anisotropy index of geochemical elements, and c is a constant.

Images