Multi-river confluence source-sink analysis method and analysis device

Abstract

The application discloses a kind of multi-river confluence material source source-sink analysis method and analysis device, and multi-river confluence material source source-sink analysis method includes: obtaining the content of each heavy mineral in each mineral sample in research area, and research area is the collection area of multiple rivers;According to the heavy mineral assignment formula, the sum of the heavy mineral content of each mineral sample is calculated, and the mineral samples in the research area are classified according to the calculation result, and the source direction of each type of mineral is analyzed;According to the classification result of mineral sample, the ratio method is used to cluster analysis on the mineral sample in each type of mineral, and a cluster diagram is obtained;Based on the cluster diagram, and combined with the source direction of each type of mineral and the geological characteristics of research area, the source direction of each mineral sample in each type of mineral in research area is determined by using fuzzy cluster analysis method;It can further finely divide the source direction, improve the accuracy of source source-sink analysis result, and the analysis operation convenience.

Description

Technical Field

[0001] This invention relates to the fields of geochemical exploration and sedimentology, specifically to a source-sink analysis method and apparatus for multi-river confluence sediments. Background Technology

[0002] In studies of provenance orientation at the confluence of multiple rivers, provenance analysis plays a certain indicative role in paleocurrents, paleogeography, and oil and gas exploration. Literature review indicates that current source-sink analysis primarily utilizes paleogeography, paleocurrent orientation, and single or heavy mineral assemblages for description or semi-quantitative analysis. However, factors such as data collection may hinder the reconstruction of the lacustrine basin paleogeography in the study area, making further analysis using the first two methods impossible. While single minerals and heavy mineral assemblages can roughly clarify provenance orientation, it is difficult to distinguish the specific provenance orientation of minerals at the confluence of similar mineral sources, and quantitative parameters and reproducible quantitative procedures are lacking. Summary of the Invention

[0003] The main objective of this invention is to propose a source-sink analysis method and apparatus for multiple river confluences, aiming to solve the aforementioned problems.

[0004] To achieve the above objectives, this invention proposes a source-sink analysis method for multiple river confluences, which includes the following steps:

[0005] The content of each heavy mineral in each mineral sample in the study area, which is the confluence area of ​​multiple rivers, was obtained.

[0006] The sum of heavy mineral content of each mineral sample is calculated according to the heavy mineral assignment formula, and the mineral samples in the study area are classified according to the calculation results, and the provenance direction of each type of mineral is analyzed.

[0007] Based on the classification results of the mineral samples, the ratio method was used to perform cluster analysis on the mineral samples of each type of mineral to obtain a cluster diagram;

[0008] Based on the cluster diagram and combined with the provenance orientation of each type of mineral and the geological characteristics of the study area, fuzzy clustering analysis was used to determine the provenance orientation of each mineral sample in each type of mineral within the study area.

[0009] Optionally, the step of obtaining the content of each heavy mineral in each mineral sample in the study area specifically includes:

[0010] Multiple mineral samples were obtained by sampling the fifth-order sequence stratigraphy in the study area.

[0011] The content of each heavy mineral in each mineral sample was calculated according to the heavy mineral content formula, which is:

[0012]

[0013] In the formula, h l The content of the l-th heavy mineral in a mineral sample;

[0014] m l This represents the mass of the l-th heavy mineral element in the mineral sample, in grams.

[0015] m represents the total mass of the mineral sample, expressed in grams.

[0016] Optionally, the steps of calculating the sum of heavy mineral contents of each mineral sample according to the heavy mineral assignment formula, classifying the mineral samples in the study area according to the calculation results, and analyzing the provenance direction of each type of mineral specifically include:

[0017] Based on the calculation results of the content of each heavy mineral in each mineral sample, the sum of the heavy mineral content of each mineral sample is calculated according to the heavy mineral assignment formula to obtain the total heavy mineral content value.

[0018] The total heavy mineral content of multiple mineral samples was calculated statistically, and mineral samples with the same integer part of the total heavy mineral content in the calculation results were classified into the same type of minerals.

[0019] Based on the decimal value of the total heavy mineral content of each mineral sample in each type of mineral, and in combination with the sampling location of the mineral sample, the source direction of the mineral samples in the same type of mineral is analyzed. The distribution direction of multiple mineral samples in the same type of mineral along the position where the decimal value increases is the source direction of the mineral of that type, and the distribution direction of multiple mineral samples in the same type of mineral along the position where the decimal value decreases is the aggregation direction of the mineral of that type.

[0020] Optionally, the formula for assigning values ​​to heavy minerals is:

[0021] H l =2 l-1 +h l ;

[0022]

[0023] In the formula, n is the number of heavy minerals contained in mineral sample k;

[0024] H l The numerical value assigned to the content of the l-th heavy mineral is dimensionless.

[0025] Z k denoted as the total content of heavy minerals contained in mineral sample k.

[0026] Optionally, the step of performing cluster analysis on mineral samples of each type using the ratio method based on the classification results of the mineral samples to obtain a cluster diagram specifically includes:

[0027] The ratio of the contents of any two heavy minerals in each mineral sample is calculated using the ratio formula, and a ratio series is obtained.

[0028] Based on the ratio series of multiple mineral samples in each type of mineral, the fuzzy similarity matrix of each type of mineral is constructed using the maximum-minimum method;

[0029] Calculate the transitive closure of each fuzzy similarity matrix, and obtain the similarity between any two mineral samples of each type of mineral based on the transitive closure to obtain a clustering graph.

[0030] Optionally, the ratio formula is:

[0031]

[0032] In the formula, h ki Let be the content of the i-th heavy mineral in mineral sample k;

[0033] h kj Let be the content of the j-th heavy mineral in mineral sample k;

[0034] h ij Let be the ratio of the content of the i-th heavy mineral to the content of the j-th heavy mineral in mineral sample k.

[0035] Optionally, the fuzzy similarity matrix is:

[0036]

[0037] In the formula, h a Let be the ratio of the content of the i-th heavy mineral to the content of the j-th heavy mineral in mineral sample a;

[0038] h b Let be the ratio of the content of the i-th heavy mineral to the content of the j-th heavy mineral in mineral sample b.

[0039] Optionally, the step of calculating the transitive closure of each fuzzy similarity matrix and obtaining the similarity between any two mineral samples of each type of mineral based on the transitive closure to obtain a clustering graph specifically includes:

[0040] The transitive closure R of each fuzzy similarity matrix is ​​calculated using the transitive closure formula;

[0041] Based on the transitive closure r, the similarity between any two mineral samples of each type of mineral is determined using the Euclidean formula.

[0042] The cluster diagram is obtained based on the calculation results.

[0043] Optionally, the transitive closure formula is:

[0044]

[0045] In the formula, m is the number of mineral samples in each type of mineral;

[0046] R is the transitive closure of each fuzzy similarity matrix;

[0047] The Euclidean formula is:

[0048]

[0049] In the formula, m is the number of mineral samples in each type of mineral;

[0050] X a This represents a one-dimensional vector observation of mineral sample a;

[0051] X b This represents a one-dimensional vector observation of mineral sample b.

[0052] d(X a X b ) represents the distance between mineral sample a and mineral sample b.

[0053] The present invention also provides a source-sink analysis device for multiple river confluences, the source-sink analysis device comprising:

[0054] The heavy mineral content acquisition module is used to acquire the content of each heavy mineral in each mineral sample in the study area;

[0055] The mineral classification module is used to calculate the sum of heavy mineral contents of each mineral sample according to the heavy mineral assignment formula, classify the mineral samples in the study area according to the calculation results, and analyze the provenance direction of each type of mineral.

[0056] The clustering analysis module is used to perform clustering analysis on mineral samples of each type of mineral based on the classification results of the mineral samples, using the ratio method, to obtain a clustering diagram; and

[0057] The provenance orientation determination module is used to determine the provenance orientation of each mineral sample in each type of mineral within the study area based on the cluster map and in combination with the provenance orientation of each type of mineral and the geological characteristics of the study area, using fuzzy clustering analysis.

[0058] In the technical solution of this invention, by assigning values ​​to the content of individual heavy minerals and the content of heavy mineral combinations, the direction of different types of provenances and the direction of paleocurrents are first roughly quantitatively distinguished. Then, based on this, cluster analysis is used to cluster the same type of provenances, further refining the division of provenance directions. This solves the problem that the division of provenance directions is controlled by human subjective factors, and improves the accuracy and convenience of provenance division. Combining sedimentology, geochemistry characteristics and cluster analysis to quantitatively analyze provenances has the characteristics of strong operability, easy replication and promotion, and relatively accurate results. Attached Figure Description

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

[0060] Figure 1 A flowchart illustrating an embodiment of the source-sink analysis method for multiple river confluences provided by the present invention;

[0061] Figure 2 for Figure 1 Flowchart of step S100;

[0062] Figure 3 for Figure 1 Flowchart of step S200;

[0063] Figure 4 for Figure 1 Flowchart of step S300 in the middle section;

[0064] Figure 5 for Figure 4 Flowchart of step S330 in the middle section;

[0065] Figure 6 This is a statistical table of heavy mineral content data in 34 mineral samples used in a specific embodiment of the present invention;

[0066] Figure 7 This is a source clustering diagram of the northwestern part of the study area in a specific embodiment of the present invention;

[0067] Figure 8 This is a source-sink analysis diagram of the northwestern part of the study area in a specific embodiment of the present invention;

[0068] Figure 9 This is a schematic diagram of an embodiment of the source-sink analysis device for multi-river confluences provided by the present invention.

[0069] Explanation of icon numbers:

[0070]

[0071] The realization of the objective, functional features and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0072] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present invention.

[0073] It should be noted that if the embodiments of the present invention involve directional indicators (such as up, down, left, right, front, back, etc.), the directional indicators are only used to explain the relative positional relationship and movement of the components in a certain specific posture (as shown in the figure). If the specific posture changes, the directional indicators will also change accordingly.

[0074] Furthermore, if the embodiments of this invention involve descriptions such as "first" or "second," these descriptions are for descriptive purposes only and should not be construed as indicating or implying their relative importance or implicitly specifying the number of technical features indicated. Therefore, a feature defined with "first" or "second" may explicitly or implicitly include at least one of those features. Additionally, the meaning of "and / or" throughout the text includes three parallel solutions; for example, "A and / or B" includes solution A, solution B, or a solution where both A and B are satisfied simultaneously. Furthermore, the technical solutions of the various embodiments can be combined with each other, but this must be based on the ability of those skilled in the art to implement them. When the combination of technical solutions is contradictory or impossible to implement, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection claimed by this invention.

[0075] In studies of provenance orientation at the confluence of multiple rivers, provenance analysis plays a certain indicative role in paleocurrents, paleogeography, and oil and gas exploration. Literature review indicates that current source-sink analysis primarily utilizes paleogeography, paleocurrent orientation, and single or heavy mineral assemblages for description or semi-quantitative analysis. However, factors such as data collection may hinder the reconstruction of the lacustrine basin paleogeography in the study area, making further analysis using the first two methods impossible. While single minerals and heavy mineral assemblages can roughly clarify provenance orientation, it is difficult to distinguish the specific provenance orientation of minerals at the confluence of similar mineral sources, and quantitative parameters and reproducible quantitative procedures are lacking.

[0076] In view of this, the present invention provides a source-sink analysis method for multiple river confluences. Figures 1 to 5 This is an embodiment of the source-sink analysis method for multi-river confluences provided by the present invention. In a specific embodiment of the present invention, the southwestern part of the Tianhuan Depression in the Ordos Basin is used as the study area for specific analysis.

[0077] Please see Figure 1 The source-sink analysis method for multiple river confluences includes the following steps:

[0078] Step S100: Obtain the content of each heavy mineral in each mineral sample in the study area, where the study area is the confluence of multiple rivers.

[0079] This step mainly involves statistical analysis of heavy mineral data from multiple mineral samples obtained by sampling in the study area.

[0080] Further, please refer to Figure 2 Step S100 specifically includes:

[0081] Step S110: Sampling is performed in the fifth-order sequence stratigraphy of the study area to obtain multiple mineral samples.

[0082] In this step, the time issue is resolved by sampling the fifth-order sequence stratigraphy of the study area, ensuring that the parent rocks of the multiple mineral samples obtained from the sampling are basically the same.

[0083] In a specific embodiment of the present invention, samples were taken from the Chang 63 section (fifth-order sequence) of the Yanchang Formation in the southwestern part of the Tianhuan Depression in the Ordos Basin to obtain 34 ore samples (such as...). Figure 6 W1-W shown 34 ).

[0084] Step S120: Calculate the content of each heavy mineral in each of the mineral samples according to the heavy mineral content formula, wherein the heavy mineral content formula is:

[0085]

[0086] In the formula, h l The content of the l-th heavy mineral in a mineral sample;

[0087] m l This represents the mass of the l-th heavy mineral element in the mineral sample, in grams.

[0088] m represents the total mass of the mineral sample, expressed in grams.

[0089] In a specific embodiment of the present invention, the 34 mineral samples, namely mineral samples W1-W 34Each element includes at least one of the four heavy mineral elements: zircon, rutile, leucite, and hematite. More specifically, according to... Figure 1 Statistical results show that mineral samples W1-W 25 All samples contain three heavy mineral elements: zircon, rutile, and leucobite, i.e., mineral samples W1-W. 25 The contents of zircon, rutile, and leucobite need to be calculated for all samples; mineral sample W 26 -W 31 All samples contain three heavy mineral elements: zircon, leucite, and hematite, i.e., mineral sample W. 26 -W 31 The contents of zircon, leucite, and hematite need to be calculated for all samples; mineral sample W 32 The mineral sample W contains two heavy mineral elements: zircon and hematite. 32 The contents of zircon and hematite need to be calculated; mineral sample W 33 -W 34 All samples contain two heavy mineral elements: leucobite and hematite, i.e., mineral sample W. 33 -W 34 All require calculation of the content of leucobrine and hematite (e.g., Figure 6 (As shown).

[0090] Step S200: Calculate the sum of heavy mineral content of each mineral sample according to the heavy mineral assignment formula, classify the mineral samples in the study area according to the calculation results, and analyze the provenance direction of each type of mineral.

[0091] Further, please refer to Figure 3 Step S200 specifically includes:

[0092] Step S210: Based on the calculation results of the content of each heavy mineral in each mineral sample, calculate the sum of the heavy mineral content of each mineral sample according to the heavy mineral assignment formula to obtain the total heavy mineral content value.

[0093] Step S220: Calculate the total heavy mineral content of multiple mineral samples, and classify mineral samples with the same integer part of the total heavy mineral content in the calculation results into the same type of mineral.

[0094] Step S230: Based on the decimal value of the total heavy mineral content of each mineral sample in each type of mineral, and in combination with the sampling location of the mineral sample, analyze the source direction of the mineral samples in the same type of mineral. The distribution direction of multiple mineral samples in the same type of mineral along the position where the decimal value increases is the source direction of the mineral of that type, and the distribution direction of multiple mineral samples in the same type of mineral along the position where the decimal value decreases is the aggregation direction of the mineral of that type.

[0095] More specifically, the formula for assigning values ​​to heavy minerals is:

[0096] H l =2 l-1 +h l (2)

[0097]

[0098] In the formula, n is the number of heavy minerals contained in mineral sample k;

[0099] H l The numerical value assigned to the content of the l-th heavy mineral is dimensionless.

[0100] Z k denoted as the total content of heavy minerals contained in mineral sample k.

[0101] This invention classifies mineral samples by their total heavy mineral content, that is, by dividing multiple mineral samples into intervals. To improve the accuracy of classification, a heavy mineral content assignment formula is used to assign values ​​to each heavy mineral content, resulting in a value for each heavy mineral content, i.e., H. l The values ​​fall within the intervals [1,2), [2,3), [4,5), [8,9), [16,17), ..., [2l-1,2l-1+1), respectively, thus allowing the total heavy mineral content of different types of mineral samples to fall within different intervals, achieving efficient classification of mineral samples.

[0102] In a specific embodiment of the present invention, the 34 mineral samples, namely mineral samples W1-W 34 Each of these elements includes at least one of four heavy mineral elements: zircon, rutile, leucite, and hematite. The formula for assigning values ​​to these heavy minerals will be explained in detail using zircon, rutile, leucite, and hematite as examples.

[0103] If zircon is treated as the first heavy mineral element, its assignment formula is H1 = 2. 1-1 +h1=1+h1;

[0104] If rutile is considered the second heavy mineral element, its assignment formula is H2 = 2. 2-1 +h2=2+h2;

[0105] If leucobite is considered the third heavy mineral element, its assignment formula is H3 = 2. 3-1 +h3=4+h3;

[0106] If hematite is classified as the fourth heavy mineral element, its assignment formula is H4 = 2. 4-1 +h4=8+h1;

[0107] Thus, combined Figure 6 Statistical data shows that the total content of heavy minerals in mineral sample W1 is: Mineral sample W 12 The total content of heavy minerals contained therein is Mineral sample W 25 The total content of heavy minerals contained therein is Mineral sample W 26 The total content of heavy minerals contained therein is Mineral sample W 31 The total content of heavy minerals contained therein is Mineral sample W 33 The total content of heavy minerals contained therein is Thus, for mineral samples W1-W 34 The total content of heavy minerals in each sample is calculated, and then the mineral samples are classified according to the integer part of the total heavy mineral content value. Specifically, mineral samples W1-W2 are classified. 25 Mineral samples W can be classified into a category if the total content of heavy minerals contains integer digits of 7. 26 -W 31 Mineral samples W can be classified into a category if the total content of heavy minerals is all in integers of 13. 33 -W 34 Mineral samples W can be classified into a category if the total content of heavy minerals is all in integers of 12. 32 Mineral samples W1-W are classified into one category if the total content of heavy minerals is an integer multiple of 9. 34 Minerals are classified into four types.

[0108] Then, based on the decimal values ​​of the total heavy mineral content of various mineral samples, the aggregation direction of this type of mineral is analyzed. In a specific embodiment of the present invention, after comprehensive analysis, it is found that the study area has a source of minerals from five directions: west, northwest, southwest, north, and northeast. That is, the source direction of minerals in the study area includes five directions: west, northwest, southwest, north, and northeast, which can reflect the basic direction of multiple rivers.

[0109] Step S300: Based on the classification results of the mineral samples, the ratio method is used to perform cluster analysis on the mineral samples of each type of mineral to obtain the cluster map of each type of mineral.

[0110] In this step, the proportions of the same type of minerals from different directions vary due to different environments. However, at the source, the proportions of heavy minerals remain relatively constant as they flow into the confluence area via a river. Therefore, based on these characteristic ratios, it can be determined which mineral samples of the same type flowed into the confluence area via the same river.

[0111] Further, please refer to Figure 4 Step S300 specifically includes:

[0112] Step S310: Calculate the content ratio of any two heavy minerals in each mineral sample using the ratio formula to obtain the ratio series.

[0113] Specifically, the ratio formula is as follows:

[0114]

[0115] In the formula, h ki Let be the content of the i-th heavy mineral in mineral sample k;

[0116] h kj Let be the content of the j-th heavy mineral in mineral sample k;

[0117] h ij Let be the ratio of the content of the i-th heavy mineral to the content of the j-th heavy mineral in mineral sample k.

[0118] It should be noted that the content ratio of any two heavy minerals is not in any particular order. Taking zircon and rutile as an example, if they have already been calculated, then the rutile / zircon ratio will not be calculated again.

[0119] Step S320: Based on the ratio series of multiple mineral samples in each type of mineral, construct the fuzzy similarity matrix of each type of mineral using the minimax method.

[0120] Furthermore, the fuzzy similarity matrix is:

[0121]

[0122] In the formula, h a Let be the ratio of the content of the i-th heavy mineral to the content of the j-th heavy mineral in mineral sample a;

[0123] h b Let be the ratio of the content of the i-th heavy mineral to the content of the j-th heavy mineral in mineral sample b.

[0124] In a specific implementation of this invention, the maximum-minimum method is used to construct mineral samples of the same type (mineral samples W1-W). 25 The fuzzy similarity matrix of is:

[0125]

[0126] Step S330: Calculate the transitive closure of each fuzzy similarity matrix, and obtain the similarity between any two mineral samples of each type of mineral based on the transitive closure to obtain a clustering graph.

[0127] Further, please refer to Figure 5 Step S330 specifically includes:

[0128] Step S331: Calculate the transitive closure R of each fuzzy similarity matrix using the transitive closure formula.

[0129] Furthermore, the transitive closure formula is:

[0130]

[0131] In the formula, m is the number of mineral samples in each type of mineral;

[0132] R is the transitive closure of each fuzzy similarity matrix;

[0133] In a specific embodiment of the present invention, minerals of the same type (mineral samples W1-W) 25 This includes 25 mineral samples, i.e., m = 25, corresponding to... but,

[0134]

[0135] Step S332: Based on the transitive closure r, use the Euclidean formula to determine the similarity between any two mineral samples of each type of mineral.

[0136] Furthermore, the Euclidean formula is:

[0137]

[0138] In the formula, m is the number of mineral samples in each type of mineral;

[0139] X a This represents a one-dimensional vector observation of mineral sample a;

[0140] X b This represents a one-dimensional vector observation of mineral sample b.

[0141] d(X a X b ) represents the distance between mineral sample a and mineral sample b.

[0142] Step S333: Obtain the cluster diagram based on the calculation results.

[0143] In a specific embodiment of the present invention, the clustering graph obtained based on the similarity calculation results is as follows: Figure 7 As shown. Based on the cluster diagram, mineral sample W... 15 W 21 W 13 W 10 W 17 W 16 W9 and W 12 Mineral samples W7 and W7, from the same river, 11 W8, W 18 W 19 and W 20 Mineral samples W4, W5, and W6 all originated from the same river, while mineral samples W1 and W6 also originated from the same river. 14 W2, W3, W 23 W 22 W 24 and W 25 They come from the same river.

[0144] Step S400: Based on the cluster diagram and combined with the provenance direction of each type of mineral and the geological characteristics of the study area, the provenance direction of each mineral sample in each type of mineral in the study area is determined by fuzzy clustering analysis.

[0145] In a specific embodiment of the present invention, based on the clustering diagram, combined with the analysis results of the provenance orientation of each type of mineral obtained in step S200 and the geological characteristics of the study area, the fuzzy similarity matrix constructed for each type of mineral is set to 0.34, and fuzzy clustering analysis is performed (e.g., Figure 8 As shown, mineral sample W... 15 W 21 W 13 W 10 W 17 W 16 W9 and W 12 ( Figure 8 (Marked by a triangle) Primarily sourced from southwestern rivers; mineral samples W7 and W 11 W8, W 18 W 19 and W 20 ( Figure 8 Mineral samples W4, W5, and W6 (marked in the middle square) mainly come from northwestern river sources; mineral samples W4, W5, and W6 mainly come from north-northwestern river sources. Figure 8 (marked with a pentagon in the middle); mineral samples W1 and W 14 W2, W3, W 23 W 22 W 24 and W 25 ( Figure 8 (Marked in a five-pointed star) mainly comes from the western river source.

[0146] It should be noted that, following step S400, the following is also included:

[0147] Step S500: Compare and analyze the provenance orientation results of each mineral sample obtained in step S400 with the numerical distribution results of heavy mineral content obtained in step S200.

[0148] This improves the accuracy of the source-sink analysis results of the multi-river confluence obtained in step S400.

[0149] In the technical solution of this invention, by assigning values ​​to the content of individual heavy minerals and the content of heavy mineral combinations, the source directions of different types of sediments and the direction of paleocurrents are first roughly quantitatively distinguished. Then, based on this, cluster analysis is used to cluster the same type of sediment sources, further refining the division of source directions. This solves the problem that the current source direction division is controlled by human subjective factors, and improves the accuracy and convenience of source division. Combining sedimentology, geochemistry characteristics and cluster analysis to quantitatively analyze sediment sources has the characteristics of strong operability, easy replication and promotion, and relatively accurate results.

[0150] This invention also provides a source-sink analysis device 100 for multi-river confluences, employing the source-sink analysis method for multi-river confluences described above. Please refer to [link to relevant documentation]. Figure 9 The multi-river confluence source-sink analysis device 100 includes a heavy mineral content acquisition module 1, a mineral classification module 2, a cluster analysis module 3, and a source direction determination module 4. The heavy mineral content acquisition module 1 is used to acquire the content of each heavy mineral in each mineral sample in the study area. The mineral classification module 2 is used to calculate the sum of the heavy mineral contents of each mineral sample according to the heavy mineral assignment formula, classify the mineral samples in the study area according to the calculation results, and analyze the source direction of each type of mineral. The cluster analysis module 3 is used to perform cluster analysis on the mineral samples in each type of mineral according to the classification results output by the mineral classification module 2, using the ratio method to obtain a cluster diagram. The source direction determination module 4 is used to determine the source direction of each mineral sample in each type of mineral in the study area based on the cluster diagram output by the cluster analysis module 3, combined with the source direction of each type of mineral output by the mineral classification module 2 and the geological characteristics of the study area, using fuzzy cluster analysis.

[0151] The above description is merely a preferred embodiment of the present invention and does not limit the patent scope of the present invention. Any equivalent structural transformations made using the contents of the present invention's specification and drawings under the inventive concept of the present invention, or direct / indirect applications in other related technical fields, are included within the patent protection scope of the present invention.

Claims

1. A source-sink analysis method for multiple river confluences, characterized in that, The source-sink analysis method for multiple river confluences includes the following steps: The content of each heavy mineral in each mineral sample in the study area, which is the confluence area of ​​multiple rivers, was obtained. The sum of heavy mineral content of each mineral sample is calculated according to the heavy mineral assignment formula, and the mineral samples in the study area are classified according to the calculation results, and the provenance direction of each type of mineral is analyzed. Based on the classification results of the mineral samples, the ratio method was used to perform cluster analysis on the mineral samples of each type of mineral to obtain a cluster diagram; Based on the cluster diagram, and combined with the provenance direction of each type of mineral and the geological characteristics of the study area, the provenance direction of each mineral sample in each type of mineral in the study area is determined by fuzzy clustering analysis. The steps of calculating the sum of heavy mineral contents of each mineral sample according to the heavy mineral assignment formula, classifying the mineral samples in the study area according to the calculation results, and analyzing the provenance direction of each type of mineral specifically include: Based on the calculation results of the content of each heavy mineral in each mineral sample, the sum of the heavy mineral content of each mineral sample is calculated according to the heavy mineral assignment formula to obtain the total heavy mineral content value. The total heavy mineral content of multiple mineral samples was calculated statistically, and mineral samples with the same integer part of the total heavy mineral content in the calculation results were classified into the same type of minerals. Based on the decimal value of the total heavy mineral content of each mineral sample in each type of mineral, and in combination with the sampling location of the mineral sample, the source direction of the mineral samples in the same type of mineral is analyzed. The distribution direction of multiple mineral samples in the same type of mineral along the position where the decimal value increases is the source direction of the mineral of that type, and the distribution direction of multiple mineral samples in the same type of mineral along the position where the decimal value decreases is the aggregation direction of the mineral of that type. The formula for assigning values ​​to heavy minerals is: ； In the formula, h l The content of the l-th heavy mineral in a mineral sample; H l The numerical value assigned to the content of the l-th heavy mineral is dimensionless.

2. The source-sink analysis method for multiple river confluences as described in claim 1, characterized in that, The steps for obtaining the content of each heavy mineral in each mineral sample in the study area specifically include: Multiple mineral samples were obtained by sampling the fifth-order sequence stratigraphy in the study area. The content of each heavy mineral in each mineral sample was calculated according to the heavy mineral content formula, which is: In the formula, h l The content of the l-th heavy mineral in a mineral sample; m l This represents the mass of the l-th heavy mineral element in the mineral sample, in grams. m represents the total mass of the mineral sample, expressed in grams.

3. The source-sink analysis method for multiple river confluences as described in claim 1, characterized in that, The step of performing cluster analysis on mineral samples of each type using the ratio method based on the classification results of the mineral samples to obtain a cluster diagram specifically includes: The ratio of the contents of any two heavy minerals in each mineral sample is calculated using the ratio formula, and a ratio series is obtained. Based on the ratio series of multiple mineral samples in each type of mineral, the fuzzy similarity matrix of each type of mineral is constructed using the maximum-minimum method; Calculate the transitive closure of each fuzzy similarity matrix, and obtain the similarity between any two mineral samples of each type of mineral based on the transitive closure to obtain a clustering graph.

4. The source-sink analysis method for multiple river confluences as described in claim 3, characterized in that, The formula for the ratio is: ； In the formula, h ki Let be the content of the i-th heavy mineral in mineral sample k; h kj Let be the content of the j-th heavy mineral in mineral sample k; h ij Let be the ratio of the content of the i-th heavy mineral to the content of the j-th heavy mineral in mineral sample k.

5. The source-sink analysis method for multiple river confluences as described in claim 3, characterized in that, The steps of calculating the transitive closure of each fuzzy similarity matrix and obtaining the similarity between any two mineral samples of each type of mineral based on the transitive closure to obtain the clustering graph specifically include: The transitive closure R of each fuzzy similarity matrix is ​​calculated using the transitive closure formula; Based on the transitive closure R, the similarity between any two mineral samples of each type of mineral is calculated using the Euclidean formula. The cluster diagram is obtained based on the calculation results.

6. A source-sink analysis device for multiple river confluences, applicable to the source-sink analysis method for multiple river confluences as described in any one of claims 1-5, characterized in that, The source-sink analysis device for multi-river confluences includes: The heavy mineral content acquisition module is used to acquire the content of each heavy mineral in each mineral sample in the study area; The mineral classification module is used to calculate the sum of heavy mineral contents of each mineral sample according to the heavy mineral assignment formula, classify the mineral samples in the study area according to the calculation results, and analyze the provenance direction of each type of mineral. The clustering analysis module is used to perform clustering analysis on mineral samples of each type of mineral based on the classification results of the mineral samples, using the ratio method, to obtain a clustering diagram; and The provenance orientation determination module is used to determine the provenance orientation of each mineral sample in each type of mineral within the study area based on the cluster map and in combination with the provenance orientation of each type of mineral and the geological characteristics of the study area, using fuzzy clustering analysis.

Images