Method and device for calculating geochemical parameters of rock, medium and equipment

Abstract

The invention discloses a rock geochemical parameter calculation method and device, a medium and equipment. The method comprises the following steps: S10, acquiring TOC and pyrolysis data; s20, establishing a type dominance index classification function, and obtaining a type dominance index TPI; s30, dividing advantage type units based on a Fisher discriminant method; s40, according to a division result, establishing a calculation model of S2 in a segmented manner; s50, the TOC of the whole well section is calculated; and S60, substituting the calculation result of the TOC of the whole well section into the calculation model to obtain an S2 value corresponding to each measuring point so as to realize continuous calculation of S2. According to the method, advantage type units are divided through a Fisher discriminant method; establishing a calculation model of S2 in a segmented manner through a division result; the TOC of the whole well section is calculated through a well logging curve; according to the method, the calculation result of the TOC of the whole well section is substituted into the calculation model, and the S2 value corresponding to each measuring point is obtained, so that the continuous calculation of S2 is realized, the absolute error between the calculation parameter and the test parameter is small, and the precision can completely meet the technical requirements of deep oil and gas exploration and development.

Description

Technical Field

[0001] This invention relates to the field of well logging evaluation technology, and in particular to a method, apparatus, medium and equipment for calculating rock geochemical parameters. Background Technology

[0002] The geochemical parameter S2 represents the weight percentage (mg / g) of organic carbon obtained from pyrolysis of organic matter within the temperature range of 300℃ to 600℃. This parameter signifies the amount of organic matter that can be converted into petroleum under current maturity conditions, i.e., the oil generation potential of kerogen. At a given maturity level, differences in S2 values ​​can indicate different types of organic matter. Accurate calculation of S2 allows for a more comprehensive understanding of the characteristics of source rocks, providing a basis for organic carbon recovery and hydrocarbon generation / explosion calculations, thereby guiding exploration and development more efficiently. Therefore, research on well logging calculation methods for the geochemical parameter S2 of source rocks is of great significance.

[0003] However, the geochemical characteristics of rocks and well logging performance are inherently related, but due to measurement bias, unclear internal mechanisms and weak correlation between the two, the accurate calculation of geochemical parameters has always been a challenge. Summary of the Invention

[0004] The main objective of this invention is to provide a method and apparatus for calculating rock geochemical parameters, thereby solving the technical problem that rock geochemical parameters are difficult to calculate accurately in the prior art.

[0005] To achieve the above objectives, the present invention provides a method for calculating rock geochemical parameters, the method comprising the following steps: S10, obtaining TOC and pyrolysis data; S20, establishing a type dominance index classification function and obtaining the type dominance index (TPI); S30, classifying dominant type units based on Fisher's discriminant method; S40, Establish a calculation model for rock geochemical parameter S2 in segments according to the division results; S50, Calculate TOC for the entire well section; S60, Substitute the calculation result of TOC for the entire well section into the calculation model to obtain the S2 value corresponding to each measuring point, thereby realizing the continuous calculation of S2.

[0006] Optionally, step S10 includes the following steps: S110, conducting total organic carbon (TOC) experiments and rock pyrolysis experiments; S120, obtaining TOC and pyrolysis data based on the experimental results, wherein the pyrolysis data includes S1, S2 and HI.

[0007] Optionally, step S20 includes the following step: S210, establishing a type advantage index classification function based on the following formula: (1) In the formula: TPI is the type dominance index, dimensionless; S1 is the content of liquid hydrocarbons preserved in the rock sample evaporated at a certain heating rate to 300℃, mg / g; HI is the hydrogen index.

[0008] Optionally, step S30 includes the following steps: S310, arranging the acquired TPIs in ascending order to create a scatter plot, using the inflection points of the curves as segmentation points, and dividing the measurement points into 3 types; S320, for different types of measurement points, establishing a discriminant function using Fisher's discriminant method and well logging data to divide the dominant type units.

[0009] Optionally, the logging data includes lithology, formation pressure, temperature, and fluid properties.

[0010] Optionally, step S50 includes the following steps: S510, obtaining logging curves based on downhole measurements and laboratory analysis; S520, calculating the TOC for the entire well section based on the logging curves.

[0011] Furthermore, to achieve the above objectives, this application embodiment also provides a calculation device for rock geochemical parameters. The device includes: a data acquisition module for acquiring TOC and pyrolysis data; a function establishment module for establishing a type dominance index classification function and acquiring the type dominance index (TPI); a unit division module for dividing dominant type units based on Fisher's discriminant method; a model establishment module for establishing a segmented calculation model for S2 based on the division results; a TOC calculation module for calculating the TOC of the entire well section; and a calculation result substitution module for substituting the calculation results of the TOC of the entire well section into the calculation model to obtain the S2 value corresponding to each measuring point, thereby realizing continuous calculation of S2.

[0012] In addition, to achieve the above objectives, embodiments of this application also provide a computer-readable storage medium including instructions that, when executed on a computer, cause the computer to perform the method for calculating rock geochemical parameters as described in any embodiment of this application.

[0013] Furthermore, to achieve the above objectives, embodiments of this application also provide a computing device, which includes: at least one processor, a memory, and an input / output unit; wherein the memory is used to store a computer program, and the processor is used to call the computer program stored in the memory to execute the rock geochemical parameter calculation method described in any embodiment of this application.

[0014] Compared with the prior art, the beneficial effects of the present invention are as follows: The method for calculating rock geochemical parameters provided in this application involves dividing the dominant type units using the Fisher discriminant method; establishing a calculation model for S2 by segmenting the division results; calculating the TOC of the entire well section using well logging curves; and obtaining the S2 value corresponding to each measuring point by substituting the calculation result of the TOC of the entire well section into the calculation model, thereby achieving continuous calculation of S2. This results in a small absolute error between the calculated parameters and the test parameters, and the accuracy fully meets the technical requirements for deep oil and gas exploration and development. Attached Figure Description

[0015] Figure 1 A flowchart illustrating the method for calculating rock geochemical parameters provided in the embodiments of this application; Figure 2 A structural block diagram of the device for calculating rock geochemical parameters provided in the embodiments of this application; Figure 3 Different types of TPI measurement point division diagrams provided for embodiments of this application; Figure 4 Well logging curves provided for embodiments of this application; Figure 5 The S2 calculation error analysis diagram provided in the embodiments of this application; Figure 6 This is a schematic diagram of the structure of the medium provided in the embodiments of this application; Figure 7 A schematic diagram of the structure of a computing device provided in an embodiment of this application.

[0016] The realization of the purpose, functional features and advantages of this application will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0017] It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of the application. Rather, these embodiments are provided to make the disclosure more thorough and complete, and to fully convey the scope of the disclosure to those skilled in the art.

[0018] To address the aforementioned technical problems, embodiments of this application provide a method for calculating rock geochemical parameters, such as... Figure 1 , Figures 3-5 As shown, the method may include the following steps: S10, obtain TOC (Total Organic Carbon) and pyrolysis data.

[0019] Specifically, pyrolysis refers to the reaction process in which a substance decomposes when heated. Pyrolysis data includes S1, S2, and HI.

[0020] In an exemplary embodiment, step S10 may include the following steps: S110 was used for core total organic carbon experiments and rock pyrolysis experiments. S120, TOC and pyrolysis data are obtained based on experimental results, wherein the pyrolysis data includes S1, S2 and HI.

[0021] Specifically, the total organic carbon (TOC) test in the core was used to obtain TOC, and the rock pyrolysis test was used to obtain S1, S2, and HI.

[0022] Furthermore, S1 represents the content of liquid hydrocarbons preserved in the rock sample evaporated at a certain heating rate to 300℃, S2 represents the rock geochemical parameters, and HI represents the hydrogen index.

[0023] S20, Establish the Type Advantage Index (TPI) classification function and obtain the Type Advantage Index (TPI).

[0024] In an exemplary embodiment, step S20 may include the following steps: S210, the type dominance index classification function is established based on the following formula: (1) In the formula: TPI is the type dominance index, dimensionless; S1 is the content of liquid hydrocarbons preserved in the rock sample evaporated at a certain heating rate to 300℃, mg / g; HI is the hydrogen index; S220, Acquire Type Advantage Index (TPI).

[0025] S30, based on Fisher's discriminant method to classify dominant type units.

[0026] In an exemplary embodiment, step S30 may specifically include the following steps: S310, Reference Figure 3 The acquired TPIs are arranged in ascending order to form a scatter plot. The inflection points of the curves are used as segmentation points to divide the measurement points into three types. S320: For different types of measuring points, a discriminant function is established using Fisher's discriminant method and well logging data to classify dominant type units. The well logging data includes lithology, formation pressure, temperature, and fluid properties.

[0027] In this exemplary embodiment, Fisher's discriminant method is a linear discriminant method that finds a linear combination that maximizes the differences between different categories while minimizing the differences within the same category. This method does not require any assumptions about the distribution of the population and has good applicability.

[0028] Furthermore, the specific steps of Fisher's discriminant method include: Data preprocessing: Standardize the data to ensure that the dimensions of each feature are consistent.

[0029] Constructing the discriminant function: Based on the idea of ​​analysis of variance, construct a discriminant function to separate data of different categories as much as possible in the projected space.

[0030] Projection: Projecting the new sample into this low-dimensional space.

[0031] Classification: Determine the category of a new sample based on its projected position.

[0032] S40, based on the division results, establish a calculation model for the rock geochemical parameter S2 in segments.

[0033] S50, calculate TOC for the entire well section, and the logging curve is as follows: Figure 4 As shown.

[0034] In an exemplary embodiment, step S50 includes the following steps: S510, based on downhole measurements and laboratory analysis, obtains logging curves; S520 calculates the total cost of the entire well section based on well logging curves.

[0035] Specifically, well logging curves are geological data obtained through well logging instruments in underground boreholes. They describe the regular changes in various physical properties of rock formations with depth. These curves can reflect the composition, porosity, permeability, density, and other characteristics of rock formations, and are important evidence for geological exploration and oil and gas resource development.

[0036] Furthermore, obtaining well logging curves typically involves two main methods: downhole measurement and laboratory analysis. Downhole measurement utilizes various logging instruments to continuously measure formation parameters during drilling, obtaining real-time logging data. Laboratory analysis, on the other hand, involves obtaining core or cuttings samples and using various analytical testing methods to perform detailed physicochemical property tests on the samples, yielding logging data. These two methods complement each other, providing comprehensive information for geological interpretation.

[0037] S60, Substitute the calculation result of the TOC of the entire well section into the calculation model to obtain the S2 value corresponding to each measuring point, and realize the continuous calculation of S2.

[0038] Furthermore, the accuracy analysis table of the S2 calculation results is shown in Table 1 below, and the error analysis graph of the S2 calculation is shown below. Figure 5 As shown. Table 1

[0039] In one or more of these exemplary embodiments, the dominant type units are divided using the Fisher discriminant method; a calculation model for S2 is established by segmenting the division results; the TOC of the entire well section is calculated using well logging curves; and the S2 value corresponding to each measuring point is obtained by substituting the calculation result of the TOC of the entire well section into the calculation model, thereby realizing the continuous calculation of S2, making the absolute error between the calculation parameters and the test parameters small, and the accuracy fully meeting the technical requirements of deep oil and gas exploration and development.

[0040] Based on the above embodiments, refer to Figure 2 Another embodiment of this application also provides a calculation device for rock geochemical parameters. The calculation device 200 for rock geochemical parameters may include the following modules: Data acquisition module 210 is used to acquire TOC and pyrolysis data; Function creation module 220 is used to create a type advantage index classification function and obtain the type advantage index (TPI). Unit partitioning module 230 is used to partition dominant type units based on Fisher's discriminant method; The model building module 240 is used to build the calculation model of S2 segment by segment based on the partitioning results; TOC Calculation Module 250 is used to calculate TOC for the entire well section. The calculation result substitution module 260 is used to substitute the calculation result of the TOC of the entire well section into the calculation model to obtain the S2 value corresponding to each measuring point, so as to realize the continuous calculation of S2.

[0041] Based on the above embodiments, this application also provides a computer-readable storage medium, see reference. Figure 6 The computer-readable storage medium shown is an optical disc 50, on which a computer program (i.e., a program product) is stored. When the computer program is run by a processor, it implements the steps described in the above-described method implementation, such as: acquiring TOC and pyrolysis data; establishing a type dominance index classification function to obtain the type dominance index TPI; dividing the dominant type units based on Fisher's discriminant method; establishing a segmented calculation model of S2 based on the division results; calculating the TOC of the entire well section; and substituting the calculation result of the TOC of the entire well section into the calculation model to obtain the S2 value corresponding to each measuring point, thereby realizing the continuous calculation of S2. The specific implementation methods of each step will not be repeated here.

[0042] It should be noted that examples of the computer-readable storage medium may also 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 optical and magnetic storage media, which will not be elaborated here.

[0043] In addition to the above embodiments, this application also provides a computing device. Figure 7 A block diagram is shown of an exemplary computing device 60 suitable for implementing embodiments of the present application. The computing device 60 may be a computer system or a server. Figure 7 The computing device 60 shown is merely an example and should not impose any limitation on the functionality and scope of use of the embodiments of this application.

[0044] like Figure 7 As shown, the components of computing device 60 may include, but are not limited to: one or more processors or processing units 601, system memory 602, and bus 603 connecting different system components (including system memory 602 and processing unit 601).

[0045] The computing device 60 typically includes a variety of computer system readable media. These media can be any available media that can be accessed by the computing device 60, including volatile and non-volatile media, removable and non-removable media.

[0046] System memory 602 may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 6021 and / or cache memory 6022. Computing device 60 may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, ROM 6023 may be used to read and write non-removable, non-volatile magnetic media (…). Figure 7 (Not shown in the image, usually referred to as "hard drive"). Although not shown in... Figure 7 The diagram illustrates that a disk drive for reading and writing to removable non-volatile disks (e.g., "floppy disks") and an optical disk drive for reading and writing to removable non-volatile optical disks (e.g., CD-ROMs, DVD-ROMs, or other optical media) can be provided. In these cases, each drive can be connected to a bus 603 connecting different system components via one or more data media interfaces. The system memory 602 may include at least one program product having a set (e.g., at least one) of program modules configured to perform the functions of the embodiments of this application.

[0047] A program / utility 6025 having a set (at least one) of program modules 6024 may be stored, for example, in system memory 602, and such program modules 6024 include, but are not limited to, an operating system, one or more application programs, other program modules, and program data, each or some combination of these examples may include an implementation of a network environment. Program modules 6024 typically perform the functions and / or methods described in the embodiments of this application.

[0048] The computing device 60 can also communicate with one or more external devices 604 (such as a keyboard, pointing device, display, etc.). This communication can be performed via input / output (I / O) interface 605. Furthermore, the computing device 60 can also communicate with one or more networks (e.g., local area networks (LANs), wide area networks (WANs), and / or public networks, such as the Internet) via network adapter 606. Figure 7 As shown, network adapter 606 communicates with other modules of computing device 60 (such as processing unit 601, etc.) via bus 603, which connects different system components. It should be understood that, although... Figure 7 Other hardware and / or software modules may be used in conjunction with computing device 60, as not shown in the diagram.

[0049] The processing unit 601 executes various functional applications and data processing by running programs stored in the system memory 602. For example, it acquires TOC and pyrolysis data; establishes a type dominance index classification function to obtain the type dominance index (TPI); divides dominant type units based on the Fisher discriminant method; establishes a segmented S2 calculation model based on the division results; calculates the TOC for the entire well section; and substitutes the TOC calculation results for the entire well section into the calculation model to obtain the S2 value corresponding to each measuring point, thus achieving continuous calculation of S2. The specific implementation methods of each step will not be repeated here. It should be noted that although several units / modules or sub-units / sub-modules of the rock geochemical parameter calculation device are mentioned in the detailed description above, this division is merely exemplary and not mandatory. In fact, according to the embodiments of this application, the features and functions of two or more units / modules described above can be embodied in one unit / module. Conversely, the features and functions of one unit / module described above can be further divided into multiple units / modules for embodiment.

[0050] In the description of this application, it should be noted that the terms "first", "second", and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance.

[0051] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0052] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. The apparatus embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the coupling or direct coupling or communication connection shown or discussed may be through some communication interface; the indirect coupling or communication connection between apparatuses or units may be electrical, mechanical, or other forms.

[0053] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0054] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0055] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0056] Finally, it should be noted that the above-described embodiments are merely specific implementations of this application, used to illustrate the technical solutions of this application, and not to limit them. The protection scope of this application is not limited thereto. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments, or make equivalent substitutions for some of the technical features, within the technical scope disclosed in this application. Such modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be covered within the protection scope of this application. Therefore, the protection scope of this application should be determined by the protection scope of the claims.

[0057] Furthermore, although the operations of the method of this application are described in a specific order in the accompanying drawings, this does not require or imply that these operations must be performed in that specific order, or that all the operations shown must be performed to achieve the desired result. Additionally or alternatively, certain steps may be omitted, multiple steps may be combined into one step, and / or one step may be broken down into multiple steps.

Claims

1. A method for calculating rock geochemical parameters, characterized in that, The method for calculating the rock geochemical parameters includes the following steps: S10, obtain TOC and pyrolysis data; S20, Establish the Type Advantage Index (TPI) classification function and obtain the Type Advantage Index (TPI); S30, dominant type units are divided based on Fisher's discriminant method; S40, based on the division results, segmented calculation models for rock geochemical parameters S2 are established; S50, calculate the TOC for the entire well section; S60, Substitute the calculation result of the TOC of the entire well section into the calculation model to obtain the S2 value corresponding to each measuring point, and realize the continuous calculation of S2.

2. The method for calculating rock geochemical parameters according to claim 1, characterized in that, Step S10 includes the following steps: S110 was used for core total organic carbon experiments and rock pyrolysis experiments. S120, TOC and pyrolysis data are obtained based on experimental results, wherein the pyrolysis data includes S1, S2 and HI.

3. The method for calculating rock geochemical parameters according to claim 1, characterized in that, Step S20 includes the following steps: S210, The type dominance index classification function is established based on the following formula: （1） In the formula: TPI is the type dominance index, dimensionless; S1 is the content of liquid hydrocarbons preserved in the rock sample evaporated at a certain heating rate to 300℃, mg / g; HI is the hydrogen index; S220, obtain Type Advantage Index (TPI).

4. The method for calculating rock geochemical parameters according to claim 1, characterized in that, Step S30 includes the following steps: S310: Arrange the acquired TPIs in ascending order to create a scatter plot, and use the inflection points of the curves as segmentation points to divide the measurement points into 3 types. S320 uses Fisher's discrimination method and well logging data to establish a discrimination function for different types of measuring points, and divides them into dominant type units.

5. The method for calculating rock geochemical parameters according to claim 4, characterized in that, The logging data includes lithology, formation pressure, temperature, and fluid properties.

6. The method for calculating rock geochemical parameters according to claim 1, characterized in that, Step S50 includes the following steps: S510, based on downhole measurements and laboratory analysis, obtains logging curves; S520 calculates the total cost of the entire well section based on well logging curves.

7. A device for calculating rock geochemical parameters, characterized in that, include: The data acquisition module is used to acquire TOC and pyrolysis data; The function creation module is used to create a type advantage index classification function and obtain the type advantage index (TPI). The unit partitioning module is used to partition dominant type units based on Fisher's discriminant method. The model building module is used to build the calculation model of S2 segment by segment based on the partitioning results; The TOC calculation module is used to calculate the TOC for the entire well section. The calculation result substitution module is used to substitute the calculation results of TOC for the entire well section into the calculation model to obtain the S2 value corresponding to each measuring point, thereby realizing the continuous calculation of S2.

8. A computer-readable storage medium, characterized in that, It includes instructions that, when run on a computer, cause the computer to perform the method for calculating the rock geochemical parameters as described in any one of claims 1 to 6.

9. A computing device, characterized in that, The computing device includes: At least one processor, memory, and input / output unit; The memory is used to store computer programs, and the processor is used to call the computer programs stored in the memory to execute the method for calculating rock geochemical parameters according to any one of claims 1 to 6.

Images