Method, device, electronic device and storage medium for obtaining dry weight of minerals

By iteratively updating the weight matrix, the problem of large errors in inversion of mineral dry weight is solved, and the calculation accuracy of mineral dry weight is improved, especially suitable for elements with less content in the strata. The results are more accurate.

CN114388074BActive Publication Date: 2025-01-14CHINA UNIV OF PETROLEUM (BEIJING)
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202210036836.1
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2022-01-13
Publication Date
2025-01-14
Estimated Expiration
2042-01-13

AI Technical Summary

Technical Problem

In the prior art, the weighted least squares method with fixed weight values ​​is used to invert the dry weight of minerals, and the measurement error of elements with less content in the formation is large, resulting in a large error of the dry weight of minerals.

Method used

By obtaining the dry weight value of the element in the formation to be measured, an inversion model from the dry weight of the element to the dry weight of the mineral is established, the initial weight matrix is ​​set to a diagonal matrix, and the weight matrix is ​​updated iteratively until the current dry weight of the mineral is the previous dry weight of the mineral is the dry weight of the mineral. The difference between weights meets the preset conditions and improves the inversion accuracy of the mineral dry weight.

Benefits of technology

The inversion accuracy of mineral dry weight is significantly improved, especially for elements with less content in the strata, which reduces errors. The obtained mineral dry weight is closer to the true value, and the core profile results are more in line with the actual strata.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN114388074B_ABST
    Figure CN114388074B_ABST
Patent Text Reader

Abstract

The present application provides a method, device, electronic device and storage medium for obtaining mineral dry weight, the method comprising: obtaining the element dry weight value of the stratum to be measured; establishing an inversion model from element dry weight to mineral dry weight according to a predetermined linear correlation between element dry weight and mineral dry weight and a random error vector; setting the initial value of the weight matrix corresponding to the inversion model; determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix; updating the weight matrix according to the currently estimated element dry weight value and the element dry weight value of the stratum to be measured, and executing again the step of determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix, until the difference between the current mineral dry weight and the previous mineral dry weight meets the preset conditions, and the mineral dry weight is obtained. The weight matrix is ​​continuously updated by the currently estimated element dry weight value and the element dry weight value obtained in the stratum to be measured, so as to improve the inversion accuracy of the mineral dry weight.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present application relates to the field of oil and gas exploration, and in particular to a method, device, electronic device and storage medium for obtaining the dry weight of a mineral. Background Art

[0002] In the field of oil and gas exploration, accurate acquisition of the dry weight of minerals in the reservoir is very important for understanding the reservoir. Since the elements in the formation are closely related to the formation of minerals and the chemical composition of minerals, the mineral content in the formation, i.e., the dry weight of minerals, is determined through the conversion relationship between element content and mineral content.

[0003] In the related art, when the element content of the stratum to be measured and the above-mentioned conversion relationship are known, the weighted least square method with a fixed weight value can usually be used to solve the problem and then invert the mineral content.

[0004] However, for elements with less content in the formation, the measured element content has a large error. When the weighted least squares method with fixed weight values ​​is used for solution, the inverted mineral dry weight has a large error. Summary of the invention

[0005] The method, device, electronic device and storage medium for obtaining the dry weight of a mineral provided in the present application are used to solve the problem of inaccurate error in the dry weight of a mineral inverted by the weighted least square method with fixed weight values ​​in the prior art.

[0006] In a first aspect, the present application provides a method for obtaining the dry weight of a mineral, the method comprising:

[0007] Obtain the dry weight value of the element of the formation to be tested;

[0008] According to the predetermined linear correlation between element dry weight and mineral dry weight and the random error vector, an inversion model from element dry weight to mineral dry weight is established;

[0009] Setting an initial value of a weight matrix corresponding to the inversion model, wherein the weight matrix is ​​a diagonal matrix;

[0010] Determine the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix;

[0011] According to the currently estimated element dry weight value and the element dry weight value of the stratum to be measured, the value of the weight matrix is ​​updated, and the step of determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix is ​​executed again until the difference between the current mineral dry weight and the previous mineral dry weight meets the preset conditions, thereby obtaining the mineral dry weight in the stratum to be measured, wherein the currently estimated element dry weight value is obtained based on the linear correlation and the current mineral dry weight value.

[0012] In a possible implementation, updating the value of the weight matrix based on the currently estimated element dry weight value and the element dry weight value of the formation to be measured includes:

[0013] Update the value of the weight matrix according to the following formula:

[0014]

[0015] Among them, the weight matrix W can be expressed as follows:

[0016]

[0017] j represents the number of updates of the weight matrix, i ranges from [1, n], n is a positive integer greater than or equal to 1, B j represents the j-th estimated element dry weight value, b represents the element dry weight value of the stratum to be measured, and the value range of p is (0, 1).

[0018] In a possible implementation, the preset condition is that the difference between the current mineral dry weight and the previous mineral dry weight converges to 0.

[0019] In a possible implementation, the method further includes:

[0020] Obtain the mineral composition of the stratum to be tested;

[0021] Establishing multiple content combinations of the mineral combination, wherein each content combination includes the content corresponding to each mineral in the mineral combination, and there is at least one distinguishing mineral between different content combinations, and the corresponding content of the distinguishing mineral in different content combinations is different;

[0022] Based on the plurality of content combinations and the relative molecular formula of each mineral in the mineral combination, a linear correlation between the element dry weight and the mineral dry weight is determined.

[0023] In a second aspect, the present application provides a device for obtaining a mineral dry weight, the device comprising: a first obtaining unit, for obtaining an element dry weight value of a formation to be measured;

[0024] A modeling unit, used for establishing an inversion model of element dry weight to mineral dry weight according to a predetermined linear correlation between element dry weight and mineral dry weight and a random error vector;

[0025] An initial unit, used to set an initial value of a weight matrix corresponding to the inversion model, wherein the weight matrix is ​​a diagonal matrix;

[0026] A first determination unit, configured to determine a current mineral dry weight value based on the element dry weight value and a current value of a weight matrix;

[0027] The second determination unit is used to update the value of the weight matrix according to the currently estimated element dry weight value and the element dry weight value of the stratum to be measured, and to again execute the step of determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix, until the difference between the current mineral dry weight and the previous mineral dry weight meets the preset conditions, thereby obtaining the mineral dry weight in the stratum to be measured, wherein the currently estimated element dry weight value is obtained based on the linear correlation and the current mineral dry weight value.

[0028] In a possible implementation manner, the first determining unit is specifically configured to update the value of the weight matrix according to the following formula:

[0029]

[0030] Among them, the weight matrix W can be expressed as follows:

[0031]

[0032] j represents the number of updates of the weight matrix, i ranges from [1, n], n is a positive integer greater than or equal to 1, B j represents the j-th estimated element dry weight value, b represents the element dry weight value of the stratum to be measured, and the value range of p is (0, 1).

[0033] In a possible implementation, the preset condition is that the difference between the current mineral dry weight and the previous mineral dry weight converges to 0.

[0034] In a possible implementation manner, the device further includes:

[0035] The second acquisition unit is used to acquire the mineral combination of the stratum to be measured;

[0036] A simulation unit, used for establishing a plurality of content combinations of the mineral combination, wherein each content combination includes the content corresponding to each mineral in the mineral combination, and there is at least one distinguishing mineral between different content combinations, and the corresponding content of the distinguishing mineral in different content combinations is different;

[0037] The third determination unit is used to establish a linear correlation between the element dry weight and the mineral dry weight based on the plurality of content combinations and the relative molecular formula of each mineral in the mineral combination.

[0038] In a third aspect, the present application provides an electronic device, including: a memory and a processor.

[0039] A memory for storing instructions executable by the processor;

[0040] Wherein, the processor is used to execute the method as described in any one of the first aspects according to the executable instructions.

[0041] In a fourth aspect, the present application provides a computer-readable storage medium, wherein the computer-readable storage medium stores computer-executable instructions, and when the computer-executable instructions are executed by a processor, they are used to implement the method as described in any one of the first aspects.

[0042] In a fifth aspect, the present application provides a computer program product, including a computer program, which, when executed by a processor, implements the method as described in any one of the first aspects.

[0043] The present application provides a method, device, electronic device and storage medium for obtaining mineral dry weight, the method comprising: obtaining the element dry weight value of the stratum to be measured; establishing an inversion model from element dry weight to mineral dry weight according to a predetermined linear correlation between element dry weight and mineral dry weight and a random error vector; setting the initial value of the weight matrix corresponding to the inversion model; determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix; updating the weight matrix according to the currently estimated element dry weight value and the element dry weight value of the stratum to be measured, and executing again the step of determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix, until the difference between the current mineral dry weight and the previous mineral dry weight meets the preset conditions, and the mineral dry weight is obtained. The weight matrix is ​​continuously updated by the currently estimated element dry weight value and the element dry weight value obtained in the stratum to be measured, so as to improve the inversion accuracy of the mineral dry weight. BRIEF DESCRIPTION OF THE DRAWINGS

[0044] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present application and, together with the description, serve to explain the principles of the present application.

[0045] Figure 1 A flow chart of a method for obtaining dry weight of a mineral provided in an embodiment of the present application;

[0046] Figure 2 A schematic diagram of the dry weight results of an element provided in an embodiment of the present application;

[0047] Figure 3 A schematic diagram for comparing the results of mineral dry weight inversion based on weighted least squares method and iterative least squares method provided in an embodiment of the present application;

[0048] Figure 4 A schematic diagram of a core profile result provided in an embodiment of the present application;

[0049] Figure 5 A schematic flow chart of a method for determining a linear correlation between an element dry weight and a mineral dry weight provided in an embodiment of the present application;

[0050] Figure 6A schematic diagram of the structure of a device for obtaining the dry weight of a mineral is provided for an embodiment of the present application;

[0051] Figure 7 A schematic diagram of the structure of another device for obtaining dry weight of minerals provided in an embodiment of the present application;

[0052] Figure 8 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application.

[0053] The above drawings have shown clear embodiments of the present application, which will be described in more detail later. These drawings and text descriptions are not intended to limit the scope of the present application in any way, but to illustrate the concept of the present application to those skilled in the art by referring to specific embodiments. DETAILED DESCRIPTION

[0054] Exemplary embodiments will be described in detail herein, examples of which are shown in the accompanying drawings. When the following description refers to the drawings, the same numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the following exemplary embodiments do not represent all implementations consistent with the present application. Instead, they are merely examples of devices and methods consistent with some aspects of the present application as detailed in the appended claims.

[0055] At present, geo-element logging is one of the important logging methods that cannot be ignored in reservoir evaluation. It can be used to evaluate the lithology of complex formations, study sedimentary environments, and determine the types and contents of minerals in reservoirs. In addition, significant progress has been made in the development of element logging instruments and theoretical research. The influence of complex formation mineral composition and wellbore environment on the relative yield of elements has also been corrected relatively accurately. Element logging has been widely used in the exploration and development of oil and gas resources. In related technologies, since there is a good linear correlation between element dry weight and mineral dry weight, the calculation of mineral dry weight is firstly to establish an element-mineral conversion matrix based on mineral type and element type, and then the matrix and element content are combined with weighted least squares method to obtain it.

[0056] In the related art, when the element content of the stratum to be measured and the above-mentioned conversion relationship are known, the weighted least square method with a fixed weight value can usually be used to solve the problem and then invert the mineral content.

[0057] However, for elements with less content in the formation, the measured element content has a large error. When the weighted least squares method with fixed weight values ​​is used for solution, the inverted mineral dry weight has a large error.

[0058] The method, device, electronic device and storage medium for obtaining the dry weight of a mineral provided in the present application are intended to solve the above technical problems of the prior art.

[0059] The technical solution of the present application and how the technical solution of the present application solves the above-mentioned technical problems are described in detail below with specific embodiments. The following specific embodiments can be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments. The embodiments of the present application will be described below in conjunction with the accompanying drawings.

[0060] Figure 1 A flowchart of a method for obtaining the dry weight of a mineral provided in an embodiment of the present application, the method comprising the following steps:

[0061] S101, obtaining the dry weight value of the element of the formation to be measured.

[0062] For example, in order to determine the dry weight of minerals in the stratum to be measured in this embodiment, the content of each element contained in the stratum to be measured, that is, the dry weight value of the element, is firstly obtained.

[0063] In one example, when determining the dry weight value of elements in the formation to be measured, the earth element logging method can be used. For example, the pulse seed source of the earth element logging instrument can be selected to emit fast neutrons into the formation to be measured. After the fast neutrons enter the formation, they successively undergo inelastic scattering and capture reactions with the elements in the formation, and a crystal detector is used to obtain the gamma energy spectrum when the formation to be measured undergoes inelastic scattering and capture reactions, and the dry weight value of each element contained in the formation to be measured is obtained from the gamma energy spectrum.

[0064] In one example, after obtaining the dry weight value of each element through geo-element logging, the obtained elements can be screened. For example, indicative elements can be screened out from the obtained elements. For example, Si, Al, Ca, Fe, K, Na, Mg and other elements can all be used as indicator elements of common minerals in the formation, such as Si as an indicator element of quartz, Ca as an indicator element of calcite, and Na and K as indicator elements of feldspar. At the same time, since the content of an element in a certain mineral is basically fixed, the content value can be obtained by the relative molecular formula of the mineral. However, since the formation basically contains O element except for special minerals such as sulfide and fluoride, O element cannot be used as an indicative element of minerals. Therefore, the dry weight value of O element can be ignored.

[0065] For example, Figure 2 A schematic diagram of the dry weight results of an element provided in an embodiment of the present application, in which the vertical axis represents the depth of the formation, and the figure includes dry weight result diagrams obtained by measuring the elements Ca, Si, Fe, Al, Mg, Na, Ti, S, Mn, TC, and T0C.

[0066] S102, establishing an inversion model from element dry weight to mineral dry weight based on a predetermined linear correlation between element dry weight and mineral dry weight and a random error vector.

[0067] For example, since there is a linear correlation between element dry weight and mineral dry weight, the generalized least squares regression model is used as the inversion model between element dry weight and mineral dry weight in this embodiment. Specifically, the following formula can be referred to:

[0068] b=Ax+e,E(e)=0,Var(e)=σ 2 Σ,Cov(e,x)=0

[0069] In the above formula, b represents the dry weight value of each element in the measured stratum; x represents the dry weight value corresponding to each mineral contained in the measured stratum; A represents the conversion matrix between the dry weight value of each element in the measured stratum and the dry weight value of each mineral. The linear relationship between the dry weight of the element and the dry weight of the mineral mentioned in this embodiment can be expressed by the following relationship: b = Ax. In the formula, e represents the random error vector, and Σ is a positive definite symmetric matrix.

[0070] When the dry weight of each element and the conversion matrix are known, the dry weight of the minerals in the formation to be measured can be obtained by solving the above inversion model. In addition, in the solution process, it is necessary to satisfy the minimum residual sum of squares RSS(x) in the following formula. The formula for the residual sum of squares is:

[0071] RSS(x)=(b-xA) T Σ -1 (b-xA)

[0072] Furthermore, according to the above residual square sum formula, the estimated value of the mineral dry weight in the measured stratum can be obtained: It can be calculated using the following formula:

[0073]

[0074] S103, setting an initial value of a weight matrix corresponding to the inversion model, where the weight matrix is ​​a diagonal matrix.

[0075] S104. Determine the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix.

[0076] For example, in this embodiment, when solving the above inversion model, the estimated value of the dry weight of the above mineral can be considered to be Solve. That is, for the matrix Σ, let the matrix be a diagonal matrix, that is, let the variance of the random vector satisfy the following formula:

[0077]

[0078] Specifically, when solving, let W = ∑ -1 , and W is used as the weight matrix corresponding to the inversion model to estimate the mineral dry weight value. Then, the estimation formula of the mineral dry weight value is:

[0079]

[0080] in, In this matrix, ω1 to ω n Their respective values ​​are non-zero numbers.

[0081] In addition, in this embodiment, in order to avoid the problem of large errors when using the fixed weight value method in the existing weighted least squares method for solution, in the solution method provided by the present application, first, an initial value is set for the weight matrix corresponding to the inversion model to obtain an initial weight matrix, and then the current mineral dry weight value can be obtained based on the initial weight matrix, the obtained element dry weight value and the above-mentioned mineral dry weight estimation formula.

[0082] In some examples, when setting the initial weight matrix, ω1 to ω in the weight matrix W can be set to n The initial value of is 1.

[0083] S105. Update the value of the weight matrix based on the currently estimated element dry weight value and the element dry weight value, and execute again the step of determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix, until the current mineral dry weight value and the previous mineral dry weight value meet the preset conditions, and obtain the mineral dry weight in the stratum to be measured, wherein the currently estimated element dry weight value is obtained based on the linear correlation and the current mineral dry weight value.

[0084] Exemplarily, after the current mineral dry weight value is obtained using the initial weight matrix, the current estimated element dry weight value can be obtained based on the current mineral dry weight value and the preset linear correlation between the mineral dry weight and the element dry weight. Afterwards, the weight matrix is ​​updated by the current estimated element dry weight value and the element dry weight value obtained by step S101. Afterwards, step S104 and step S105 are continuously executed in a loop until the current mineral dry weight value and the previous mineral dry weight value meet the preset conditions, then the above steps are stopped and the currently obtained mineral dry weight value is used as the mineral dry weight value in the measured stratum. For example, the preset condition can be set to that the difference between the current mineral dry weight value and the previous mineral dry weight value is less than or equal to the preset value.

[0085] In one example, the weight matrix of this time may be updated by the difference between the currently estimated element dry weight value and the element dry weight value obtained from the formation to be measured.

[0086] In this embodiment, when determining the dry weight of minerals in the stratum to be measured, in order to avoid the problem that the mineral dry weight value obtained by using the traditional weighted least squares method, i.e., using a fixed weight value method to solve the mineral dry weight is inaccurate, especially for some elements with less content in some strata, the inverted mineral dry weight value has a large error, in this embodiment, when inverting the mineral dry weight, the weight matrix is ​​continuously updated by the currently estimated element dry weight value and the element dry weight value obtained in the stratum to be measured, until the current mineral dry weight value and the previous mineral dry weight value meet the preset conditions, and the current mineral dry weight is used as the final determined mineral dry weight value of the stratum to be measured, and then the weight matrix is ​​updated through multiple iterations to improve the inversion accuracy of the mineral dry weight.

[0087] In some specific embodiments, when updating the values ​​of the weight matrix based on the currently estimated element dry weight values ​​and the element dry weight values ​​of the formation to be measured, the weight matrix may be updated by using an iterative reweighted least squares method, that is, the weight matrix is ​​updated according to the following formula:

[0088]

[0089] Where j represents the number of updates of the weight matrix, i ranges from [1, n], n is a positive integer greater than or equal to 1, B j represents the j-th estimated dry weight of the element, b represents the dry weight of each element in the measured stratum, and the value range of p is (0, 1).

[0090] In determining the B in the above formula j When , it can be obtained based on the mineral dry weight value calculated for the jth time, which can be expressed by the following formula:

[0091] B j =Ax j

[0092] In this embodiment, the weight matrix and the inverse are updated by iterative reweighted least squares, thereby avoiding the linear non-convexity problem when calculating the mineral dry weight value by the inversion model, and improving the accuracy of the final mineral dry weight. Especially for some elements with low content in some formations, the above-mentioned method for obtaining mineral dry weight improves the calculation accuracy of mineral dry weight, and the core profile finally obtained is more consistent with the real formation.

[0093] For example, Figure 3 A schematic diagram for comparing the results of mineral dry weight inversion based on weighted least squares method and iterative least squares method provided in an embodiment of the present application.

[0094] Figure 3 In the figure, the vertical axis represents the depth of the stratum, which mainly includes 7 kinds of minerals, namely minerals 1 to 7 in the figure. Figure 3 In the figure, the left box represents the mineral dry weight inversion result obtained by iterative reweighted least squares method, and the right box represents the mineral dry weight inversion result obtained by weighted least squares method, and the value of the weight matrix in the weighted least squares method here is fixed.

[0095] contrast Figure 3 After the iterative weighted least squares method of mineral dry weight obtained by two methods is applied to the mineral dry weight inversion, compared with the weighted least squares method with fixed weight value, the core profile curve obtained has no mutation value and the smoothness is significantly enhanced. For example, the profile curves of the areas pointed by the two arrows in the same frame can be compared. It can be seen that the left profile curve is smoother than the profile curve of the area pointed by the right arrow in the same frame. Compared with the weighted least squares method, the iterative reweighted least squares method eliminates the errors in the inversion process, making the core profile results more consistent with the actual formation conditions, which is conducive to the study of formation storage.

[0096] For example, Figure 4 A schematic diagram of a core profile result provided in an embodiment of the present application. Figure 4 It can be seen that by inferring the rock skeleton density through mineral combination, and then obtaining the skeleton porosity and comparing it with the core porosity calibration, it is found that each depth point has a high degree of consistency, indicating that the new method proposed in this application has a good application effect in obtaining mineral dry weight by earth element logging.

[0097] In some embodiments, in order to determine the linear correlation between the dry weight value of the element and the dry weight value of the mineral in the measured formation, Figure 5 The following is a flow chart of a method for determining the linear correlation between the dry weight of an element and the dry weight of a mineral provided in an embodiment of the present application. Figure 5 As shown, the following steps are included:

[0098] S201. Obtain the mineral composition of the stratum to be measured.

[0099] For example, when determining the mineral composition in the formation to be measured, it can be determined based on core analysis data or other logging data of the formation to be measured.

[0100] In one example, when determining the mineral combination, it is possible to consider selecting the main minerals with a larger content in the stratum to be tested.

[0101] S202, establishing multiple content combinations of the mineral combination, and determining that each content combination includes the content corresponding to each mineral in the mineral combination, and there is at least one distinguishing mineral between different content combinations, and the corresponding content of the distinguishing mineral in different content combinations is different.

[0102] For example, in this embodiment, after the mineral combination is determined, multiple content combinations are simulated, wherein each content combination includes the mineral content corresponding to each mineral in the mineral combination, and the content of at least one mineral is different between different content combinations.

[0103] S203. Determine the linear correlation between the element dry weight and the mineral dry weight based on the multiple content combinations and the relative molecular formula of each mineral in the mineral combination.

[0104] For example, in this embodiment, after simulating multiple groups of content combinations, since the percentage content of an element in a certain mineral is fixed, for example, when the mineral is quartz, the relative molecular formula contained in the mineral at this time is SiO2, and the weight percentage of Si in the molecular formula can be calculated by the molecular formula to be 46.7%. Furthermore, after the mineral combination and the content combination corresponding to the mineral combination are known, the content of each element corresponding to each content combination (i.e., the dry weight value of the element) can be determined. Afterwards, when the content combination and the content of each element corresponding to the content combination are known, the conversion matrix in the linear correlation relationship between the element dry weight value and the mineral dry weight value can be determined, and then the linear correlation between the element dry weight value and the mineral dry weight value can be established.

[0105] For example, suppose that the mineral assemblage includes n kinds of minerals, and each mineral includes m kinds of oxides or carbonate components, that is, m kinds of elements are selected. Then the corresponding relationship between the mineral dry weight values ​​of these n kinds of minerals and the dry weight values ​​of m kinds of elements can be expressed by the following relationship:

[0106]

[0107] Among them, A ij It is expressed as the response relationship between the jth mineral and the ith element, i = 1, 2, ..., n, j = 1, 2, ..., m,; b i is the dry weight of the ith element. ij is the value of the i-th row and j-th column in the conversion matrix A. That is, there exists a conversion matrix A that can satisfy the linear correlation Ax=b between the mineral dry weight value x and the element dry weight value b.

[0108] According to the above linear correlation, it can be deduced that the relationship satisfied by the transformation matrix A is:

[0109] A=(x T x) -1 x T b

[0110] Then, based on the multiple content combinations corresponding to the pre-simulated mineral combinations, the dry weight values ​​of each element corresponding to each content combination are inferred, and then the conversion matrix A is determined through the relationship satisfied by the above conversion matrix A.

[0111] For example, the following Table 1 is a schematic diagram of a conversion matrix provided in an embodiment of the present application.

[0112] Table 1

[0113]

[0114] In this embodiment, by simulating multiple groups of simulated strata with different content combinations under the same mineral combination, a conversion matrix in the linear correlation between the dry weight value of elements in the strata of the mineral combination and the dry weight value of minerals is deduced. In practical applications, the mineral combination can be continuously adjusted according to the lithology in the actual strata or the core analysis data in the strata to be suitable for different strata.

[0115] Figure 6 A schematic diagram of a device for obtaining dry weight of a mineral is provided in an embodiment of the present application, and the device comprises:

[0116] The first acquisition unit 41 is used to acquire the dry weight value of the element of the formation to be measured;

[0117] A modeling unit 42 is used to establish an inversion model of element dry weight to mineral dry weight according to a predetermined linear correlation between element dry weight and mineral dry weight and a random error vector;

[0118] An initialization unit 43 is used to set the initial value of the weight matrix corresponding to the inversion model, and the weight matrix is ​​a diagonal matrix;

[0119] A first determination unit 44 is used to determine the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix;

[0120] The second determination unit 45 is used to update the value of the weight matrix according to the currently estimated element dry weight value and the element dry weight value of the stratum to be measured, and to re-execute the step of determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix, until the difference between the current mineral dry weight and the previous mineral dry weight meets the preset conditions, thereby obtaining the mineral dry weight in the stratum to be measured, wherein the currently estimated element dry weight value is obtained based on the linear correlation and the current mineral dry weight value.

[0121] In a possible implementation, the first determining unit 44 is specifically configured to update the value of the weight matrix according to the following formula:

[0122]

[0123] Among them, the weight matrix W can be expressed as follows:

[0124]

[0125] j represents the number of updates of the weight matrix, i ranges from [1, n], n is a positive integer greater than or equal to 1, B j represents the j-th estimated element dry weight value, b represents the element dry weight value of the stratum to be measured, and the value range of p is (0, 1).

[0126] In a possible implementation, the preset condition is that the difference between the current mineral dry weight and the previous mineral dry weight converges to 0.

[0127] The device provided in this embodiment is used to implement the technical solution provided by the above method. Its implementation principle and technical effect are similar and will not be described in detail.

[0128] Figure 7 A schematic diagram of the structure of another mineral dry weight acquisition device provided in an embodiment of the present application, in Figure 6 Based on the device structure shown, the device also includes:

[0129] A second acquisition unit 46 is used to acquire the mineral composition of the stratum to be measured;

[0130] The simulation unit 47 is used to establish multiple content combinations of the mineral combination, wherein each content combination includes the content corresponding to each mineral in the mineral combination, and there is at least one distinguishing mineral between different content combinations, and the corresponding content of the distinguishing mineral in different content combinations is different;

[0131] The third determination unit 48 is used to establish a linear correlation between the element dry weight and the mineral dry weight based on the multiple content combinations and the relative molecular formula of each mineral in the mineral combination.

[0132] The device provided in this embodiment is used to implement the technical solution provided by the above method. Its implementation principle and technical effect are similar and will not be described in detail.

[0133] Figure 8 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application is shown in FIG. Figure 8 As shown, the electronic device includes:

[0134] The electronic device includes a processor 291 and a memory 292; it may also include a communication interface 293 and a bus 294. The processor 291, the memory 292, and the communication interface 293 may communicate with each other through the bus 294. The communication interface 293 may be used for information transmission. The processor 291 may call the logic instructions in the memory 294 to execute the method of the above embodiment.

[0135] In addition, the logic instructions in the above-mentioned memory 292 can be implemented in the form of software functional units and can be stored in a computer-readable storage medium when sold or used as an independent product.

[0136] The memory 292 is a computer-readable storage medium that can be used to store software programs and computer executable programs, such as program instructions / modules corresponding to the methods in the embodiments of the present application. The processor 291 executes functional applications and data processing by running the software programs, instructions, and modules stored in the memory 292, that is, implementing the methods in the above method embodiments.

[0137] The memory 292 may include a program storage area and a data storage area, wherein the program storage area may store an operating system and an application required for at least one function; the data storage area may store data created according to the use of the terminal device, etc. In addition, the memory 292 may include a high-speed random access memory and may also include a non-volatile memory.

[0138] An embodiment of the present application provides a computer-readable storage medium, in which computer-executable instructions are stored. When the computer-executable instructions are executed by a processor, they are used to implement the method provided in the above embodiment.

[0139] An embodiment of the present application provides a computer program product, including a computer program, which implements the method provided in the above embodiment when executed by a processor.

[0140] Those skilled in the art will readily appreciate other embodiments of the present application after considering the specification and practicing the invention disclosed herein. The present application is intended to cover any modification, use or adaptation of the present application, which follows the general principles of the present application and includes common knowledge or customary techniques in the art that are not disclosed in the present application. The specification and examples are intended to be exemplary only, and the true scope and spirit of the present application are indicated by the following claims.

[0141] It should be understood that the present application is not limited to the precise structures that have been described above and shown in the drawings, and that various modifications and changes may be made without departing from the scope thereof. The scope of the present application is limited only by the appended claims.

Claims

1. A method for obtaining the dry weight of a mineral, characterized in that: include: Obtain the dry weight value of the element of the formation to be tested; According to the predetermined linear correlation between element dry weight and mineral dry weight and the random error vector, an inversion model from element dry weight to mineral dry weight is established; Setting an initial value of a weight matrix corresponding to the inversion model, wherein the weight matrix is ​​a diagonal matrix; Determine the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix; According to the currently estimated element dry weight value and the element dry weight value of the stratum to be measured, the value of the weight matrix is ​​updated, and the step of determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix is ​​performed again, until the difference between the current mineral dry weight and the previous mineral dry weight meets the preset conditions, and the mineral dry weight in the stratum to be measured is obtained, wherein the currently estimated element dry weight value is obtained based on the linear correlation and the current mineral dry weight value; The updating of the value of the weight matrix based on the currently estimated element dry weight value and the element dry weight value of the formation to be measured includes: Update the value of the weight matrix according to the following formula: 1 / h i (j+1) =|B j -b| (p-2) / 2 Among them, the weight matrix W can be expressed as follows: j represents the number of updates of the weight matrix, i ranges from [1, n], n is a positive integer greater than or equal to 1, B j represents the j-th estimated element dry weight value, b represents the element dry weight value of the stratum to be measured, and the value range of p is (0, 1).

2. The method according to claim 1, characterized in that The preset condition is that the difference between the current mineral dry weight and the previous mineral dry weight converges to 0.

3. The method according to claim 1 or 2, characterized in that: The method further comprises: Obtain the mineral composition of the stratum to be tested; Establishing multiple content combinations of the mineral combination, wherein each content combination includes the content corresponding to each mineral in the mineral combination, and there is at least one distinguishing mineral between different content combinations, and the corresponding content of the distinguishing mineral in different content combinations is different; Based on the plurality of content combinations and the relative molecular formula of each mineral in the mineral combination, a linear correlation between the element dry weight and the mineral dry weight is determined.

4. A device for obtaining the dry weight of a mineral, characterized in that: include: A first acquisition unit is used to acquire the dry weight value of the element of the formation to be measured; A modeling unit, used for establishing an inversion model of element dry weight to mineral dry weight according to a predetermined linear correlation between element dry weight and mineral dry weight and a random error vector; An initial unit, used to set an initial value of a weight matrix corresponding to the inversion model, wherein the weight matrix is ​​a diagonal matrix; A first determination unit, configured to determine a current mineral dry weight value based on the element dry weight value and a current value of a weight matrix; A second determination unit is used to update the value of the weight matrix according to the currently estimated element dry weight value and the element dry weight value of the stratum to be measured, and to perform again the step of determining the current mineral dry weight value based on the element dry weight value and the current value of the weight matrix, until the difference between the current mineral dry weight and the previous mineral dry weight meets a preset condition, thereby obtaining the mineral dry weight in the stratum to be measured, wherein the currently estimated element dry weight value is obtained based on the linear correlation and the current mineral dry weight value; The first determining unit is specifically used to update the value of the weight matrix according to the following formula: 1 / oh i (j+1) =|B j -b| (p-2) / 2 Among them, the weight matrix W can be expressed as follows: j represents the number of updates of the weight matrix, i ranges from [1, n], n is a positive integer greater than or equal to 1, B j represents the j-th estimated element dry weight value, b represents the element dry weight value of the stratum to be measured, and the value range of p is (0, 1).

5. The device according to claim 4, characterized in that The preset condition is that the difference between the current mineral dry weight and the previous mineral dry weight converges to 0.

6. The device according to claim 4 or 5, characterized in that The device also includes: The second acquisition unit is used to acquire the mineral combination of the stratum to be measured; A simulation unit, used for establishing a plurality of content combinations of the mineral combination, wherein each content combination includes the content corresponding to each mineral in the mineral combination, and there is at least one distinguishing mineral between different content combinations, and the corresponding content of the distinguishing mineral in different content combinations is different; The third determination unit is used to establish a linear correlation between the element dry weight and the mineral dry weight based on the plurality of content combinations and the relative molecular formula of each mineral in the mineral combination.

7. An electronic device, characterized in that: include: Memory, processor; Memory; a memory for storing instructions executable by the processor; Wherein, the processor is used to execute the method according to any one of claims 1-3 according to the executable instructions.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions, which are used to implement the method according to any one of claims 1 to 3 when executed by a processor.

Citation Information

Patent Citations

  • Shale gas reservoir brittleness earthquake prediction method and device

    CN105954804A

  • Soil heavy metal bioavailability estimation method, system and apparatus

    CN108549620A