Method for determining lithological type of igneous rock, storage medium and computer device

Abstract

The present invention provides a method for determining rock lithology, a storage medium and a computer device. The method comprises: S100: obtaining the molecular formula of each of a plurality of pure minerals; S200: obtaining the element types and the proportional coefficients of each element type in the target igneous rock; S300: based on the molecular formula of each of the plurality of pure minerals, according to the element types and the proportional coefficients of each element type in the target igneous rock, determining the mineral types and the mass content of each mineral type in the target igneous rock by an inversion method; S400: determining the igneous rock types and the mass content of each igneous rock type in the target igneous rock; S500: determining the lithology type of the target igneous rock according to the igneous rock types and the mass content of each igneous rock type in the target igneous rock. The method can be used to effectively determine the quantitative relationship between the mineral content and the elements of the igneous rock.

Description

Technical Field

[0001] The present invention relates to the technical field of lithology type determination, and in particular to a method for determining the lithology type of igneous rocks, a storage medium and a computer device. Background Art

[0002] Jin Yunzhi, Gao Chuqiao and others established a relationship model between the stratigraphic elements of igneous rocks and their lithology based on the analysis of stratigraphic components. Using the optimization algorithm, starting from the contents of stratigraphic elements such as Al, Ca, Fe, H, Si, S, and Ti in the stratigraphic element logging data, they quantitatively calculated the contents of various lithological components of igneous rocks. Some studies, based on standard grayscale maps such as QAPF of igneous rocks, quickly realized the automatic classification and naming of igneous rocks by projecting points of igneous rock mineral components, obtaining grayscale values ​​of projected point positions, and quickly obtaining igneous rock names based on hash table sets. Some studies have comprehensively used magnetic data, seismic data, core data, and logging data to predict igneous reservoirs. This method can not only reasonably use magnetism to identify igneous rocks, but also use three-dimensional seismic frameworks to make up for the lack of spatial resolution. And under the constraints of drilling data, it can achieve the best prediction effect of igneous reservoirs. Some studies use the combined characteristics of gravity, magnetic and electric anomalies to identify the lithology of igneous rocks. By collecting and recording the physical properties of the core and using k-means clustering to perform cluster analysis on the density, magnetic susceptibility and resistivity of the rock, the core clustering code is obtained. The gravity and magnetic anomalies are stripped from the gravity and magnetic exploration results to obtain the gravity, magnetic and electric anomaly codes. The same clustering method is used for three-dimensional spatial anomaly coding, and the rock physical property clustering codes are combined with the gravity, magnetic and electric anomaly codes to reflect the lithology of each igneous rock in the exploration area.

[0003] At present, the interpretation methods for obtaining mineral mass content can be summarized into two types: one is the empirical interpretation method represented by Schlumberger, and the other is the optimization interpretation method represented by Halliburton. Schlumberger's interpretation method is to obtain the relative content of igneous rock formation elements and mineral mass content by empirical relationship; Halliburton's interpretation method is to directly invert the formation mineral mass content by using the optimization method from the formation element content. The interpretation method represented by Schlumberger relies on a large amount of regional experimental data. For some new formations, it is necessary to re-establish the model, and the applicability of the newly established model still needs further verification; the interpretation method represented by Halliburton has strong adaptability, but there is a problem of uncertainty in inverted minerals. Figure 1 This is the interpretation result of minerals in igneous rocks by a foreign software in the prior art, where WPYR represents pyrite, WCAR represents ash content, WQFM represents silica content, and WCLA represents mud content. It can be seen that the foreign software can only calculate the mud content, silica content, ash content and pyrite content, and cannot invert all mineral types and their contents in igneous rocks.

[0004] The mineral types and contents of igneous rock formations are complex. To accurately determine the mineral mass content of igneous rock formations, the igneous rock core experimental method is usually used, or it is determined through a multi-mineral interpretation model. The number of experimental cores is limited, and the experimental results are not obtained in a timely manner. The problems with the results of the multi-mineral interpretation model are: the igneous rock mineral types do not match, and the errors in the calculated mineral content are large. Summary of the invention

[0005] The main purpose of the present invention is to provide a method for determining rock lithology, a storage medium and a computer device to accurately determine the type and mass content of minerals in igneous rocks, and then determine the lithology of igneous rocks.

[0006] In a first aspect, the present application provides a method for determining the lithology type of igneous rock, comprising the following steps: S100: obtaining the molecular formula of each of a plurality of pure minerals; S200: obtaining the element types in the target igneous rock and the proportional coefficients of each element type; S300: based on the molecular formula of each of the plurality of pure minerals, according to the element types in the target igneous rock and the proportional coefficients of each element type, determining the mineral types contained in the target igneous rock and the mass content of each mineral type by an inversion method; S400: determining the igneous rock types contained in the target igneous rock according to the mineral types contained in the target igneous rock, and determining the mass content of each igneous rock type contained in the target igneous rock according to the mass content of each mineral type contained in the target igneous rock; S500: determining the lithology type of the target igneous rock according to the igneous rock types contained in the target igneous rock and the mass content of each igneous rock type.

[0007] In one embodiment, S100 includes: for each pure mineral, performing the following steps: S110: determining the relative molecular weight of the pure mineral according to the mass of the pure mineral and a preset molar volume; S120: determining the molecular formula model of the pure mineral according to the relative molecular weight of the pure mineral and the types of elements contained therein, wherein the molecular formula model includes acid and base radicals and at least one main element other than the acid and base radicals; S130: determining the residual relative molecular weight of the at least one main element other than the acid and base radicals according to the relative molecular weight of the pure mineral and the relative molecular weight of its acid and base radicals; S140: determining the distribution coefficient of each main element in the molecular formula model by an inversion method according to the residual relative molecular weight and the relative atomic mass of each main element, thereby obtaining the molecular formula of the pure mineral.

[0008] In one embodiment, the distribution coefficient of each main element in the molecular formula model is determined by an inversion method, including: constructing a relationship model among the residual relative molecular weight, the relative atomic mass of each main element and the distribution coefficient of each main element in the molecular formula model; based on the relationship model, using the least squares method to determine the distribution coefficient of each main element in the pure mineral in the molecular formula model.

[0009] In one embodiment, the relationship model is constructed using the following formula:

[0010] Y=ΣA i X i

[0011] Wherein, Y represents the residual relative molecular weight, A i represents the relative atomic mass of the ith principal element, and X represents the distribution coefficient of the ith principal element in the molecular formula model.

[0012] In one embodiment, after S140, the method further includes: S150: based on the molecular formula of the pure mineral, determining the distribution coefficient of each element type in the pure mineral, and determining a first error of the distribution coefficient of each element type, comparing the first error with a first error threshold, and when the first error is greater than or equal to the first error threshold, returning to execute S120 to redetermine the molecular formula model of the pure mineral.

[0013] In one embodiment, the first error of the partition coefficient of each element type is determined using the following formula:

[0014]

[0015] Wherein, F(x,a) represents the error objective function of the inversion, x represents the distribution coefficient of each element type, a represents the residual relative molecular weight, i represents the number of measurements, and a i represents the residual relative molecular weight of the ith measurement, m represents the number of element types in the molecular formula of the pure mineral, and f i (x,z) represents the mass content of the element measured for the i-th time, σ i represents the measurement error of the element mass content measured for the i-th time, τ i represents the response error of the residual relative molecular weight measured for the ith time, g j (x) is the jth constraint condition for the distribution coefficient x, τ j is the error of the jth constraint, and p is the number of constraints.

[0016] In one embodiment, in S300, the mineral types contained in the target igneous rock and the mass content of each mineral type are determined by an inversion method using the following formula:

[0017]

[0018] Among them, W n represents the mass content of the nth mineral, N represents the number of mineral types, C in represents the proportion coefficient of the element measured in the nth mineral, y i Represents the mass content of the element measured for the i-th time.

[0019] In one embodiment, after S300 and before S400, the method further includes: S600: determining a second error in the mass content of each mineral type in the target igneous rock, comparing the second error with a second error threshold, and when the second error is greater than or equal to the second error threshold, returning to execute S140 in S100 to redetermine the distribution coefficient of each main element in each pure mineral in the molecular formula model by an inversion method.

[0020] In one embodiment, the second error of the mass content of each mineral type in the target igneous rock is calculated using the following formula:

[0021]

[0022] Among them, F(s,b) represents the error objective function of inversion, s represents the mass content of pure minerals, b represents the mass content of elements, i represents the number of measurements, and b i represents the mass content of the element measured for the i-th time, q represents the number of pure mineral species, and f i (s, c) represents the mass content of the mineral measured for the i-th time, δ i represents the measurement error of the mass content of the element measured for the i-th time, ρ i represents the response error of the element mass measured for the i-th time, g j (s) is the jth constraint condition of s, ρ j is the error of the jth constraint, and p is the number of constraints.

[0023] In a second aspect, the present application provides a storage medium storing computer program code, which, when executed by a processor, implements the steps of the method for determining the lithology type of igneous rock as described above.

[0024] In a third aspect, the present application provides a computer device, including a processor and a storage medium storing program code, wherein when the program code is executed by the processor, the steps of the method for determining the lithology type of igneous rock as described above are implemented.

[0025] The present invention obtains the formation element content based on the original element logging data or XRF logging, but due to the uncertainty of the relationship between the mineral mass and the element content in the igneous rock formation, it is necessary to determine the mineral content coefficient of the element according to different types of pure minerals, and finally calculate the mineral mass content of the igneous rock formation using an optimization algorithm. The application of this method can effectively determine the quantitative relationship between the igneous rock mineral content and the elements, providing a theoretical basis for identifying the mineral type and mass content of the igneous rock formation, and playing a positive role in accelerating the exploration and development of domestic igneous rock formations. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] The drawings constituting a part of the present application are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0027] Figure 1 The existing technology is the interpretation result of mineral types in igneous rocks by a certain software outside China;

[0028] Figure 2 A flowchart of a method for determining lithology type of igneous rock according to an exemplary embodiment of the present application;

[0029] Figure 3 A flowchart of a method for determining the lithology type of igneous rock according to a specific embodiment of the present application;

[0030] Figure 4 is a coefficient matrix of each element type and its distribution coefficient in each pure mineral according to a specific embodiment of the present application;

[0031] FIG. 5A to FIG. 5C The left, middle and right parts of the schematic diagram of mineral content in an igneous rock formation according to a specific embodiment of the present application are shown in sequence;

[0032] FIG. 6A to FIG. 6C The left, middle and right parts of the schematic diagram of mineral content in an igneous rock formation according to another specific embodiment of the present application are shown in sequence;

[0033] 7A to 7C The left, middle and right parts are respectively a schematic diagram of the mineral content and igneous rock type in an igneous rock formation according to a specific embodiment of the present application. DETAILED DESCRIPTION

[0034] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present invention will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0035] Embodiment 1

[0036] This embodiment provides a method for determining the lithology type of igneous rock. Figure 2 FIG. 1 is a flow chart of a method for determining the lithology type of igneous rock according to an exemplary embodiment of the present application. Figure 2 As shown, the method comprises the following steps:

[0037] S100: Get the molecular formulas of various pure minerals.

[0038] Specifically, for each pure mineral, the molecular formula of the pure mineral can be determined by the following steps:

[0039] S110: Determine the relative molecular weight of the pure mineral according to the mass of the pure mineral and the preset molar volume. Specifically, the relative molecular weight of the pure mineral can be determined by dividing the mass of the pure mineral by the molar volume, wherein the molar volume can be preset as the molar volume of a material similar to the pure mineral.

[0040] S120: Determine a molecular formula model of the pure mineral according to the relative molecular weight of the pure mineral and the types of elements contained therein, wherein the molecular formula model includes acid and base radicals and at least one main element other than acid and base radicals.

[0041] The types of elements contained in pure minerals can be obtained by analyzing the ECS element logging data or XRF data.

[0042] S130: Determine the remaining relative molecular weight of the at least one main element other than the acid and base radicals according to the relative molecular weight of the pure mineral and the relative molecular weight of the acid and base radicals thereof.

[0043] S140: According to the residual relative molecular weight and the relative atomic mass of each main element, the distribution coefficient of each main element in the molecular formula model is determined by an inversion method, and then the molecular formula of the pure mineral is obtained.

[0044] Among them, the distribution coefficient of each main element in the molecular formula model is determined by an inversion method, including: constructing a relationship model among the residual relative molecular weight, the relative atomic mass of each main element and the distribution coefficient of each main element in the molecular formula model; based on the relationship model, using the least squares method to determine the distribution coefficient of each main element in the pure mineral in the molecular formula model.

[0045] The relationship model can be constructed using the following formula:

[0046] Y=ΣA i X i

[0047] Wherein, Y represents the residual relative molecular weight, A i represents the relative atomic mass of the ith principal element, and X represents the distribution coefficient of the ith principal element in the molecular formula model.

[0048] In some cases, the molecular formula obtained by only one inversion may be inaccurate. Therefore, after S140, the method further includes: S150: based on the molecular formula of the pure mineral, determining the distribution coefficient of each element type in the pure mineral, and determining a first error of the distribution coefficient of each element type, comparing the first error with a first error threshold, and when the first error is greater than or equal to the first error threshold, returning to execute S120 to redetermine the molecular formula model of the pure mineral.

[0049] The first error of the distribution coefficient for each element species can be determined using the following formula:

[0050]

[0051] Wherein, F(x,a) represents the error objective function of the inversion, x represents the distribution coefficient of each element type, a represents the residual relative molecular weight, i represents the number of measurements, and a i represents the residual relative molecular weight of the ith measurement, m represents the number of element types in the molecular formula of the pure mineral, and f i (x,z) represents the mass content of the element measured for the i-th time, σ i represents the measurement error of the element mass content measured for the i-th time, τ i represents the response error of the residual relative molecular weight measured for the ith time, g j (x) is the jth constraint condition for the distribution coefficient x, τ j is the error of the jth constraint, and p is the number of constraints.

[0052] After the molecular formulas of the plurality of pure minerals are determined, the process proceeds to S200 .

[0053] S200: Obtaining the element types in the target igneous rock and the proportional coefficients of the element types. Specifically, the element types in the target igneous rock can be obtained by deconstructing the ECS element logging data or XRF data.

[0054] S300: Based on the molecular formulas of the plurality of pure minerals, according to the element types in the target igneous rock and the proportional coefficients of the element types, the mineral types contained in the target igneous rock and the mass content of the mineral types are determined by an inversion method.

[0055] In S300, the mineral types contained in the target igneous rock and the mass content of each mineral type are determined by an inversion method using the following formula:

[0056]

[0057] Among them, W nrepresents the mass content of the nth mineral, N represents the number of mineral types, C in represents the proportion coefficient of the element measured in the nth mineral, y i Represents the mass content of the element measured for the i-th time.

[0058] After S300 and before S400, the method further includes: S600: determining a second error in the mass content of each mineral type in the target igneous rock, comparing the second error with a second error threshold, and when the second error is greater than or equal to the second error threshold, returning to execute S140 in S100 to redetermine the distribution coefficient of each main element in each pure mineral in the molecular formula model by an inversion method.

[0059] The second error of the mass content of each mineral type in the target igneous rock can be calculated using the following formula:

[0060]

[0061] Among them, F(s, b) represents the error objective function of the inversion, s represents the mass content of pure minerals, b represents the mass content of elements, i represents the number of measurements, and b i represents the mass content of the element measured for the i-th time, q represents the number of pure mineral species, and f i (s,c) represents the mass content of the mineral measured for the i-th time, δ i represents the measurement error of the mass content of the element measured for the i-th time, ρ i represents the response error of the element mass measured for the i-th time, g j (s) is the jth constraint condition of s, ρ j is the error of the jth constraint, and p is the number of constraints.

[0062] S400: Determine the type of igneous rock contained in the target igneous rock according to the mineral types contained in the target igneous rock, and determine the mass content of each igneous rock type contained in the target igneous rock according to the mass content of each mineral type contained in the target igneous rock.

[0063] S500: Determine the lithology type of the target igneous rock according to the igneous rock types contained in the target igneous rock and the mass content of each igneous rock type.

[0064] The present invention obtains the relative content of formation elements based on the original element logging data or XRF logging, but due to the uncertainty of the relationship between the mineral mass and the element content of the igneous rock formation, it is necessary to determine the mineral content coefficient of the element according to different types of pure minerals, and finally calculate the mineral mass content of the igneous rock formation using an optimization algorithm. The application of this method can effectively determine the quantitative relationship between the igneous rock mineral content and the elements, providing a theoretical basis for identifying the mineral type and mass content of the igneous rock formation, and playing a positive role in accelerating the exploration and development of domestic igneous rock formations.

[0065] Embodiment 2

[0066] This embodiment provides a specific embodiment of a method for determining the lithology type of igneous rock. Figure 3 FIG. 1 is a flow chart of a method for determining the lithology type of igneous rock according to a specific embodiment of the present application. Figure 3 As shown, the method may include the following steps:

[0067] (1) Select the types of pure minerals that may be contained in igneous rocks, and obtain the types of elements contained in each pure mineral type and the distribution coefficient of each element type in the molecular formula of the pure mineral through ECS logging data or XRF data analysis. The types of pure minerals may include: olivine, pyroxene, albite, calcium feldspar, amphibole, quartz, zircon, muscovite, biotite, apatite, zeolite, orthoclase, nepheline, etc.

[0068] (2) Measure the mass of the pure mineral and calculate the molar mass and relative molecular weight of the pure mineral based on the selected molar volume.

[0069] (3) Determine the molecular formula model of the pure mineral, and use A(X) to represent the main elements in the molecular formula model except acid and base radicals and their distribution coefficients in the molecular formula model, where A is the relative atomic mass of the main element contained in the mineral, and X is the distribution coefficient of the main element.

[0070] (4) According to the decomposition results in (1), the main element A contained in the pure mineral is determined. A may include but is not limited to elements such as Si, Al, Fe, Ca, K, Ti, Mg, Mn, Na, P, and Zr.

[0071] (5) Calculate the remaining relative molecular weight Y of all main elements except acid and base radicals in the relative molecular weight of pure minerals. Y is equal to the relative molecular weight of pure minerals minus the relative molecular weight of acid and base radicals.

[0072] (6) Establish the relationship between Y and the distribution coefficient X, Y = ΣAiXi.

[0073] (7) Calculate the distribution coefficient X of each main element using the least square method, and use the error objective function to determine whether the distribution coefficient is appropriate. If it is appropriate, determine the molecular formula of the pure mineral. If it is not appropriate, return to step (3) and redetermine the molecular formula model of the pure mineral.

[0074] (8) By repeatedly executing the above steps (1) to (7), the molecular formulas of various pure minerals are obtained, that is, the main elements and acid and base radicals contained in each pure mineral are determined.

[0075] (9) For each pure mineral, according to its molecular formula, calculate the distribution coefficient of each element in the molecular formula in the pure mineral. The distribution coefficient is equal to the product of the relative atomic mass of each element and the number of atoms of the element in the molecular formula divided by the relative molecular weight of the mineral. Finally, determine the distribution coefficient of each element in all pure minerals, and form a coefficient matrix of the distribution coefficients of multiple pure minerals and each element.

[0076] (11) Determine the element types and proportional coefficients of each element type in the igneous rock formation by interpreting the ECS logging data or XRF data.

[0077] (12) Based on the coefficient matrix of the distribution coefficients of various pure minerals and various elements, the mineral components contained in the igneous rock and the mass content of each mineral component are determined according to the proportional coefficients of various element types in the igneous rock formation, and the type of igneous rock contained in the igneous rock formation is determined according to the mineral composition.

[0078] (13) Determine the lithology type of the igneous rock in the igneous rock formation based on the type of igneous rock contained in the igneous rock formation.

[0079] This embodiment first determines the molecular formula of each of the multiple pure minerals, and then generates a coefficient matrix of the distribution coefficient of each element type in each pure mineral based on the molecular formula of each of the multiple pure minerals. Based on the coefficient matrix, the mineral types contained in the igneous rock formation and the mass content of each mineral are inverted according to the spectrum interpretation results of the ECS logging data or XRF data, and then the igneous rock type is determined. The whole process realizes modular operation, which plays a positive role in forming igneous rock interpretation software with its own intellectual property rights to adapt to the rapid and efficient development of domestic igneous rocks.

[0080] Embodiment 3

[0081] This embodiment provides a specific application example of a method for determining the lithology type of igneous rock, which is applied to igneous rock.

[0082] Olivine is a common mineral in igneous rocks. Based on the analysis of element logging data, we can first assume that the molecular formula model of olivine is (MgFeCa)SiO4, but we cannot determine the distribution coefficients between the main elements Mg, Fe, and Ca.

[0083] In order to determine the distribution coefficient of each main element, firstly, olivine can be taken from the target layer, its molar mass (mass / molar volume) can be calculated, and the relative molecular weight of olivine can be determined according to the molar mass; then the molecular formula model of olivine is set to ASiO4, A represents the main element, and the main element A is an element with a relatively high content of Mg, Fe, Ca, etc. in olivine except Si and O; the residual relative molecular weight Y of the main element is calculated according to the relative molecular weight of olivine and the relative molecular weight of acid and base radicals; the respective distribution coefficients (C1, C2, C3) of the main elements such as Mg, Fe, Ca, etc. are determined in combination with the spectrum analysis results of the element logging data of olivine, that is, the distribution coefficients of each main element are determined by the least square method according to the results of multiple measuring points by combining the expression Y=28*C1+56*C2+40*C3, and C1=0.6, C1=0.4, C3=0 are obtained by calculation, and finally the molecular formula of olivine is determined to be Mg3Fe2(SiO4)5.

[0084] Similarly, the molecular formulas of the various pure minerals are determined, thereby determining the coefficient matrix of the distribution coefficients of the various elements in the various pure minerals. Figure 4 It is a coefficient matrix of each element type and its distribution coefficient in each pure mineral according to a specific embodiment of the present application.

[0085] FIG. 5A to FIG. 5C The left, middle and right parts are respectively a schematic diagram of the mineral content in an igneous rock formation according to a specific embodiment of the present application. In the schematic diagram of the mineral content in the igneous rock formation, the first track is the formation depth; the second to tenth tracks are the relative contents of the logging elements, the reconstructed element response values ​​and the confidence intervals of the inversion results; the eleventh track is the error target function curve; the twelfth to seventeenth tracks are the inverted mineral contents, which are respectively the mud content, silica content, calcium content, iron content, gypsum mineral content and other mineral contents; the eighteenth track is the inverted mineral mass content combined profile.

[0086] FIG. 6A to FIG. 6C The left, middle and right parts of the mineral content diagram in the igneous rock formation according to another specific embodiment of the present application are shown in sequence. For the data in the mineral content diagram in the igneous rock formation, the processing depth is 5320-5350 meters. The meanings and FIG. 5A to FIG. 5C The same schematic diagram of the mineral content in the igneous rock formations is shown.

[0087] According to the interpretation results of the mineral content of the igneous rock formation, the mineral composition of the formation is obtained. Since the igneous rock minerals are formed in a specific intrusive environment, the igneous rock type of the formation is judged. By analyzing the minerals in the igneous rock and combining the igneous rock type and the relationship between the mineral composition and the intrusive environment, the igneous rock type contained in the igneous rock formation is determined, as shown in Table 1.

[0088] Table 1

[0089]

[0090] 7A to 7C The left, middle and right parts of the schematic diagram of the mineral content and igneous rock type in the igneous rock formation according to a specific embodiment of the present application are shown in sequence. The schematic diagram of the mineral content and igneous rock type in the igneous rock formation is the calculation result of 5400-5430m of a well. The inversion result shows that the formation contains 19 types of minerals. According to the inversion mineral result, it is further determined that the igneous rock intrusion environment is a basic or ultrabasic environment, mainly basalt.

[0091] The present invention uses pure minerals to establish a coefficient matrix between the type of formation minerals and the element content based on the relative content of igneous formation elements from original element logging / XRF data, and then uses an optimization algorithm to calculate the mineral mass content of the igneous formation to further determine the type of igneous formation.

[0092] Embodiment 4

[0093] This embodiment provides a storage medium storing computer program codes. When the computer program codes are executed by a processor, the steps of the method for determining the lithology type of igneous rock as described above are implemented.

[0094] These program codes can also be loaded onto a computer or other programmable data processing device so that a series of operation steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions executed on the computer or other programmable device for implementing the process. Figure 1 The steps of a specified function in a process or multiple processes.

[0095] Storage media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer-readable instructions, data structures, modules of programs or other data. Examples of storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic disk storage or other magnetic storage devices or any other non-transmission medium that can be used to store information that can be accessed by a computing device.

[0096] Embodiment 5

[0097] This embodiment provides a computer device, including a processor and a storage medium storing program code, and when the program code is executed by the processor, the steps of the method for determining the lithology type of igneous rock as described above are implemented.

[0098] In one embodiment, a computer device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0099] The memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash memory (FLASH RAM). The memory is an example of a computer-readable medium.

[0100] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. When the terms "include" and / or "include" are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or their combinations.

[0101] It is to be understood that the terms so used are interchangeable under appropriate circumstances such that the embodiments of the application described herein are, for example, capable of operation in sequences other than those illustrated or described herein.

[0102] It should be understood that the exemplary embodiments in this specification can be implemented in a variety of different forms and should not be interpreted as being limited to the embodiments described herein. These embodiments are provided to make the disclosure of this application thorough and complete, and to fully convey the concepts of these exemplary embodiments to those of ordinary skill in the art, and should not be interpreted as limiting the present invention.

Claims

1. A method for determining the lithology type of igneous rock, characterized in that: The following steps are involved: S100: Get the molecular formulas of various pure minerals; S200: Obtaining the element types and the proportion coefficients of the element types in the target igneous rock; S300: Based on the molecular formulas of the plurality of pure minerals, according to the element types in the target igneous rock and the proportional coefficients of the element types, the mineral types and the mass content of the mineral types in the target igneous rock are determined by an inversion method; S400: determining the type of igneous rock contained in the target igneous rock according to the mineral types contained in the target igneous rock, and determining the mass content of each igneous rock type contained in the target igneous rock according to the mass content of each mineral type contained in the target igneous rock; S5 00: Determine the lithology type of the target igneous rock based on the igneous rock types contained in the target igneous rock and the mass content of each igneous rock type.

2. The method for determining the lithology type of igneous rock according to claim 1, characterized in that: S100, comprising: for each pure mineral, performing the following steps: S110: determining the relative molecular weight of the pure mineral according to the mass of the pure mineral and the preset molar volume; S120: determining a molecular formula model of the pure mineral according to the relative molecular weight of the pure mineral and the types of elements contained therein, wherein the molecular formula model includes acid and base radicals and at least one main element other than the acid and base radicals; S130: Determine the remaining relative molecular weight of the at least one main element other than the acid and base radicals according to the relative molecular weight of the pure mineral and the relative molecular weight of the acid and base radicals thereof; S140: According to the residual relative molecular weight and the relative atomic mass of each main element, the distribution coefficient of each main element in the molecular formula model is determined by an inversion method, and then the molecular formula of the pure mineral is obtained.

3. The method for determining the lithology type of igneous rock according to claim 2, characterized in that: The distribution coefficient of each main element in the molecular formula model is determined by an inversion method, including: Constructing a relationship model among the residual relative molecular weight, the relative atomic mass of each main element and the distribution coefficient of each main element in the molecular formula model; Based on the relationship model, the distribution coefficient of each main element in the pure mineral in the molecular formula model is determined using the least square method.

4. The method for determining the lithology type of igneous rock according to claim 3, characterized in that: The relationship model is constructed using the following formula: Y=ΣA i X i Wherein, Y represents the residual relative molecular weight, A i represents the relative atomic mass of the ith principal element, and X represents the distribution coefficient of the ith principal element in the molecular formula model.

5. The method for determining the lithology type of igneous rock according to claim 2, characterized in that: After S140, the method further includes: S150: Based on the molecular formula of the pure mineral, determine the distribution coefficient of each element type in the pure mineral, and determine the first error of the distribution coefficient of each element type, compare the first error with the first error threshold, and when the first error is greater than or equal to the first error threshold, return to execute S120 to redetermine the molecular formula model of the pure mineral.

6. The method for determining the lithology type of igneous rock according to claim 5, characterized in that: The first error of the distribution coefficient for each element species is determined using the following formula: Wherein, F(x,a) represents the error objective function of the inversion, x represents the distribution coefficient of each element type, a represents the residual relative molecular weight, i represents the number of measurements, and a i represents the residual relative molecular weight of the ith measurement, m represents the number of element types in the molecular formula of the pure mineral, and f i (x,z) represents the mass content of the element measured for the i-th time, σ i represents the measurement error of the element mass content measured for the i-th time, τ i represents the response error of the residual relative molecular weight measured for the ith time, g j (x) is the jth constraint condition for the distribution coefficient x, τ j is the error of the jth constraint, and p is the number of constraints.

7. The method for determining the lithology type of igneous rock according to claim 1, characterized in that: In S300, the mineral types contained in the target igneous rock and the mass content of each mineral type are determined by an inversion method using the following formula: Among them, W n represents the mass content of the nth mineral, N represents the number of mineral types, C in represents the proportion coefficient of the element measured in the nth mineral, y i Represents the mass content of the element measured for the i-th time.

8. The method for determining the lithology type of igneous rock according to claim 2, characterized in that: After S300 and before S400, the method further includes: S600: Determine the second error of the mass content of each mineral type in the target igneous rock, compare the second error with the second error threshold, and when the second error is greater than or equal to the second error threshold, return to execute S140 in S100 to redetermine the distribution coefficient of each main element in each pure mineral in the molecular formula model through an inversion method.

9. The method for determining the lithology type of igneous rock according to claim 8, characterized in that: The second error of the mass content of each mineral type in the target igneous rock is calculated using the following formula: Among them, F(s,b) represents the error objective function of inversion, s represents the mass content of pure minerals, b represents the mass content of elements, i represents the number of measurements, and b i represents the mass content of the element measured for the i-th time, q represents the number of pure mineral species, and f i (s, c) represents the mass content of the mineral measured for the i-th time, δ i represents the measurement error of the mass content of the element measured for the i-th time, ρ i represents the response error of the element mass measured for the i-th time, g j (s) is the jth constraint condition of s, ρ j is the error of the jth constraint, and p is the number of constraints.

10. A storage medium storing computer program code, characterized in that: When the computer program code is executed by a processor, the steps of the method for determining the lithology type of igneous rock according to any one of claims 1 to 9 are implemented.

11. A computer device comprising a processor and a storage medium storing program codes, wherein when the program codes are executed by the processor, the steps of the method for determining the lithology type of igneous rocks according to any one of claims 1 to 9 are implemented.

Images