Method, device, equipment and readable storage medium for determining ablation amount

Abstract

The present application provides a method, device, equipment and readable storage medium for determining the amount of erosion. The present application can determine the temperature of each primary fluid inclusion of each sample collected from a geological vertical section where an unconformity exists, and convert the temperature data of each primary fluid inclusion in each sample into pressure data; based on the pressure data of each primary fluid inclusion in each sample, a linear fit is performed on the pressure data distribution of each sample of the upper geological layer and the lower geological layer located on the unconformity, and the formation pressure gradient is determined based on the fitting result; finally, the amount of erosion of the unconformity in the geology is determined based on the fitting result and the formation pressure gradient, thereby providing important reference parameters for geologists to conduct geological research and explore for oil.

Description

Technical Field

[0001] The present application relates to the field of geological exploration technology, and in particular to a method, device, equipment and readable storage medium for determining the amount of erosion. Background Art

[0002] The thickness of the denudation of sedimentary discontinuities in sedimentary basins has always been an important parameter that geologists have been trying to explore. It is one of the essential key factors in the fields of geological research and oil exploration. Determining the thickness of the stratum denudation is an important basic work for studying the evolution history of basins and conducting quantitative evaluation of oil and gas resources. By calculating the denudation of the stratum to determine the maximum burial depth of the stratum, geologists can determine the oil generation period and gas generation period of the source rock below the Tertiary System, and then accurately evaluate the potential of oil and gas resources in geology and select the optimal exploration target. This is particularly important for the exploration of oil and gas resources below the Tertiary System.

[0003] Therefore, how to determine the erosion thickness of sedimentary discontinuities in sedimentary basins has always been a concern. Summary of the invention

[0004] The present application aims to solve at least one of the above-mentioned technical defects. In view of this, the present application provides a method, device, equipment and readable storage medium for determining the amount of erosion, which are used to solve the technical problem of determining the amount of erosion.

[0005] A method for determining an amount of ablation, comprising:

[0006] Determine the temperature of each primary fluid inclusion in each sample collected from a geologic vertical section where the unconformity exists;

[0007] Convert the temperature data of each primary fluid inclusion in each sample into pressure data;

[0008] Based on the pressure data of each primary fluid inclusion in each sample, determine the pressure value and standard error of each sample;

[0009] Based on the pressure value and standard error of each sample, linear fitting is performed on the pressure data distribution of each sample of the upper geological layer and the lower geological layer located on the unconformity surface to obtain an upper pressure fitting line and a lower pressure fitting line;

[0010] Calculating the pressure value corresponding to the upper layer intersection point of the unconformity surface and the upper layer pressure fitting line, and calculating the pressure value corresponding to the lower layer intersection point of the unconformity surface and the lower layer pressure fitting line;

[0011] Determining a formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line;

[0012] The erosion amount of the unconformity surface is determined based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient.

[0013] Preferably, the pressure value and standard error of each sample are determined based on the pressure data of each primary fluid inclusion in each sample, including:

[0014] The pressure average and standard error of the pressure data of each primary fluid inclusion in each sample are calculated as the pressure value and standard error of each sample.

[0015] Preferably, before converting the temperature data of each primary fluid inclusion in each sample into pressure data, the method further comprises:

[0016] Determine the composition of individual primary fluid inclusions in each sample;

[0017] Determining a computational method for converting temperature data of each primary fluid inclusion in each sample into pressure data based on the composition of each primary fluid inclusion in each sample;

[0018] Based on this, the temperature data of each primary fluid inclusion in each sample is converted into pressure data, including:

[0019] Based on the calculation method corresponding to each primary fluid inclusion in each sample, the temperature data of each primary fluid inclusion in each sample is converted into pressure data.

[0020] Preferably, the formation pressure gradient includes an overlying formation pressure gradient and an underlying formation pressure gradient, and determining the formation pressure gradient based on the overlying formation pressure fitting line and the underlying formation pressure fitting line includes:

[0021] Calculating the inverse of the slope of the upper layer pressure fitting line as the overlying formation pressure gradient;

[0022] The inverse of the slope of the underlying formation pressure fitting line is calculated as the underlying formation pressure gradient.

[0023] Preferably, the formation pressure gradient includes an overlying formation pressure gradient and an underlying formation pressure gradient, and determining the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper layer intersection point and the pressure value corresponding to the lower layer intersection point and the formation pressure gradient includes:

[0024] Determine the magnitude of the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point;

[0025] If there is an overlying stratum on the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, then the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point divided by the pressure gradient of the underlying stratum; if the pressure value corresponding to the lower intersection point is less than the pressure value corresponding to the upper intersection point, then the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the lower intersection point and the surface pressure value divided by the pressure gradient of the overlying stratum.

[0026] Preferably, the determining the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient, further includes:

[0027] If there is no overlying stratum above the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, then the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the lower intersection point and the surface pressure value divided by the pressure gradient of the underlying stratum.

[0028] A device for determining an amount of ablation, comprising:

[0029] a temperature determination unit for determining the temperature of each primary fluid inclusion in each sample collected from a geological vertical section where an unconformity exists;

[0030] A conversion unit, used to convert the temperature data of each primary fluid inclusion in each sample into pressure data;

[0031] A first calculation unit is used to determine the pressure value and standard error of each sample based on the pressure data of each primary fluid inclusion in each sample;

[0032] A fitting unit is used to perform linear fitting on the pressure data distribution of each sample of the upper geological layer and the lower geological layer located on the unconformity surface based on the pressure value and standard error of each sample, so as to obtain an upper pressure fitting line and a lower pressure fitting line;

[0033] A second calculation unit is used to calculate the pressure value corresponding to the intersection point between the unconformity surface and the upper layer of the upper layer pressure fitting line, and calculate the pressure value corresponding to the intersection point between the unconformity surface and the lower layer of the lower layer pressure fitting line;

[0034] a pressure gradient determining unit, configured to determine a formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line;

[0035] The erosion amount determination unit is used to determine the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient.

[0036] Preferably, the execution process of the first computing unit includes:

[0037] The pressure average and standard error of the pressure data of each primary fluid inclusion in each sample are calculated as the pressure value and standard error of each sample.

[0038] An ablation amount determination device includes: one or more processors, and a memory;

[0039] The memory stores computer-readable instructions, and when the computer-readable instructions are executed by the one or more processors, the steps of the method for determining the amount of erosion as described in any one of the above descriptions are implemented.

[0040] A readable storage medium stores computer-readable instructions. When the computer-readable instructions are executed by one or more processors, the one or more processors implement the steps of the method for determining the amount of erosion as described in any one of the above descriptions.

[0041] It can be seen from the above technical scheme that the embodiment of the present application can determine the temperature of each primary fluid inclusion of each sample collected from the geological vertical section where the unconformity exists. Since the pressure data of each primary fluid inclusion in each sample can reflect the amount of erosion of the unconformity, the temperature data of each primary fluid inclusion in each sample can be converted into pressure data; based on the pressure data of each primary fluid inclusion in each sample, the pressure value and standard error of each sample are determined; based on the pressure value and standard error of each sample, the pressure value and standard error of each sample located at the upper geological layer and the lower geological layer of the unconformity are respectively determined. The pressure data distribution of the sample is linearly fitted to obtain the upper pressure fitting line and the lower pressure fitting line; after obtaining the upper pressure fitting line and the lower pressure fitting line, the pressure value corresponding to the upper intersection of the unconformity surface and the upper pressure fitting line can be calculated, and the pressure value corresponding to the lower intersection of the unconformity surface and the lower pressure fitting line can be calculated; the formation pressure gradient can be determined based on the upper pressure fitting line and the lower pressure fitting line; finally, the erosion amount of the unconformity surface can be determined based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient. The present application can determine the erosion amount of the unconformity surface based on the pressure data of each primary fluid inclusion of each sample collected on the unconformity surface, providing important reference parameters for geologists to conduct geological research and explore for oil. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0043] Figure 1 A flow chart of a method for determining the amount of erosion provided in an embodiment of the present application;

[0044] Figure 2 A schematic diagram of the fitting effect of pressure data of a primary fluid inclusion exemplified in an embodiment of the present application;

[0045] Figure 3 A schematic diagram of the fitting effect of pressure data of a primary fluid inclusion exemplified in an embodiment of the present application;

[0046] Figure 4 A schematic diagram of the structure of an apparatus for determining the amount of ablation as exemplified in an embodiment of the present application;

[0047] Figure 5 This is a hardware structure block diagram of an erosion amount determination device disclosed in an embodiment of the present application. DETAILED DESCRIPTION

[0048] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0049] Combine the following Figure 1 , introduces the process of the method for determining the amount of ablation given in the embodiment of the present application, which may include the following steps:

[0050] Step S101, determining the temperature of each primary fluid inclusion in each sample collected from a geological vertical section where an unconformity exists.

[0051] Specifically, if a set of rock layers is deposited in an area, and later the geology of the area rises to the surface and suffers erosion, resulting in a long period of sedimentation interruption, and then descends again to accept sedimentation, that is, a certain period of strata is missing between the strata deposited successively, resulting in discontinuity of the upper and lower strata. This contact relationship between the upper and lower strata is called unconformity contact. In geology, there is an erosion surface between the upper and lower strata of the geology with unconformity contact relationship, and this surface is called unconformity. Therefore, in order to determine the erosion of the unconformity, a vertical section can be selected in the geology with unconformity, and multiple samples can be collected at different heights of the overlying strata and underlying strata of the unconformity surface of the vertical section to analyze the erosion amount of the unconformity. There are multiple primary fluid inclusions in the geology with unconformity. It is a common phenomenon that primary fluid inclusions appear in mineral crystals. Primary fluid inclusions are almost formed at the same time as the main mineral and by the same substance. After the primary fluid is filled in the crystal defects, it can be immediately sealed by the main mineral that continues to grow, with basically no material leakage and the volume basically unchanged. Therefore, the primary fluid inclusion is a representative of the original mineralization, diagenetic solution or magma melt. As a sample of ore-forming fluid, the primary fluid inclusion is one of the most important typological characteristics of minerals. By studying the primary fluid inclusion, reliable data can be provided for solving some geological problems. For example, the temperature of the primary fluid inclusion reflects the erosion of the unconformity surface in geology.

[0052] Among them, the temperature of the primary fluid inclusion can be the uniform temperature, freezing point temperature, and initial melting temperature of the primary fluid inclusion.

[0053] For example, the homogenization temperature of primary fluid inclusions is the instantaneous temperature when the two-phase or multi-phase inclusions at room temperature are transformed into the original uniform single-phase fluid when the temperature rises to a certain degree after artificial heating. It is generally believed that the homogenization temperature of primary fluid inclusions is the lower limit of the mineral formation temperature. After pressure correction, the approximate mineral formation temperature, i.e., the capture temperature of primary fluid inclusions, can be obtained.

[0054] The initial melting temperature and freezing point temperature of primary fluid inclusions are important parameters for determining the composition of the salt system of primary fluid inclusions.

[0055] Among them, a hot and cold stage can be used to perform homogenization measurement on the hydrocarbon and water components in the primary fluid inclusions and their phase change characteristics, and the initial melting temperature, freezing point temperature and homogenization temperature of each primary fluid inclusion of each sample can be obtained.

[0056] Therefore, the amount of erosion of the unconformity in the geology can be analyzed by determining the temperature of each sample collected in the vertical section of the geology where the unconformity exists, so as to analyze the erosion of the unconformity.

[0057] In particular, before determining the temperature of each primary fluid inclusion of each sample, each sample collected from different depths of the geology where the unconformity exists can be made into thin sections to determine the temperature of each primary fluid inclusion in each sample. In order to obtain a sufficient number of primary fluid inclusions for observation, a sample can be ground into multiple thin sections for observation.

[0058] Step S102, converting the temperature data of each primary fluid inclusion in each sample into pressure data.

[0059] Specifically, from the above steps, it can be seen that the formation of primary fluid inclusions runs through the entire process. The relevant data of primary fluid inclusions record and preserve the physical or chemical characteristics of different stages of geological actions. Therefore, the pressure data of primary fluid inclusions can also reflect the erosion of unconformity surfaces in geology.

[0060] The pressure of primary fluid inclusions is the fluid pressure when the primary fluid inclusions are captured. It is not easy to directly obtain the pressure data of primary fluid inclusions. Generally, the pressure data of primary fluid inclusions is obtained by using the temperature data of primary fluid inclusions.

[0061] Therefore, after determining the temperature data of each primary fluid inclusion in each sample, the temperature data of each primary fluid inclusion in each sample can be converted into pressure data in order to further analyze the erosion amount of the unconformity surface.

[0062] Step S103, based on the pressure data of each primary fluid inclusion in each sample, determine the pressure value and standard error of each sample.

[0063] It is understandable that in numerical analysis, standard error can generally be used to illustrate the relative deviation between the actual value and its estimated value, and is mainly used to measure the representativeness of the regression equation. The smaller the value of the estimated standard error, the smaller the approximate error between the estimated value and its true value.

[0064] Therefore, after determining the pressure value of each sample collected in the geological vertical section where the unconformity exists, the erosion amount of the unconformity can be further analyzed in combination with the standard error of the pressure value of each sample, so as to better analyze the erosion of the unconformity.

[0065] Specifically, after converting the temperature data of each sample into pressure data, the pressure value and standard error of each sample can be determined based on the pressure data of each primary fluid inclusion in each sample, so as to better analyze the erosion amount of the unconformity surface.

[0066] Step S104, based on the pressure value and standard error of each sample, linear fitting is performed on the pressure data distribution of each sample of the upper geological layer and the lower geological layer located on the unconformity surface to obtain an upper pressure fitting line and a lower pressure fitting line.

[0067] Specifically, it can be seen from the above steps that the pressure value and standard error of each sample can be determined based on the pressure data of each sample.

[0068] It is understandable that in various physical and statistical problems, some data sets are obtained by multiple observations or experiments on relevant variables. They are scattered, not only inconvenient to process, but also usually cannot accurately and fully reflect their inherent laws. In order to obtain the inherent laws between data or use current data to predict expected data, the functional relationship between the coordinates represented by discrete points on the plane can be approximately described or compared with continuous curves. The corresponding problems in high-dimensional space also belong to this category. In particular, when the standard error of the pressure value of each sample is relatively small, the least squares method can be used to linearly fit the distribution of the pressure values ​​of each sample.

[0069] Therefore, after determining the pressure value and standard error of each sample, the pressure data distribution of each sample of the unconformity surface and the upper and lower geological layers located on the unconformity surface can be linearly fitted based on the pressure value and standard error of each sample to obtain the upper pressure fitting line and the lower pressure fitting line.

[0070] For example, the pressure value and the stratigraphic depth of each sample collected from the geology with unconformity surface can be set as the observed variables, x is set as the pressure value of each sample, and y is the stratigraphic depth of each sample. y can be expressed as a function of x. Assuming that this functional relationship has been theoretically determined by the actual problem, y can be called a theoretical function, but it contains several unknown parameters x. In numerical analysis, curve fitting is the use of analytical expressions to approximate discrete data, that is, the formulation of discrete data. Linear fitting is a type of curve fitting. Therefore, linear fitting operations can be performed on the pressure values ​​of each sample collected from the geology with unconformity surface, and the upper pressure fitting line and the lower pressure fitting line in the geology can be obtained.

[0071] Step S105 , calculating the pressure value corresponding to the upper layer intersection point of the unconformity surface and the upper layer pressure fitting line, and calculating the pressure value corresponding to the lower layer intersection point of the unconformity surface and the lower layer pressure fitting line.

[0072] Specifically, it can be seen from the above steps that the pressure data distribution of each sample of the upper geological layer and the lower geological layer located on the unconformity surface can be linearly fitted to obtain the upper pressure fitting line and the lower pressure fitting line.

[0073] Based on this, the pressure value corresponding to the upper intersection point of the unconformity surface and the upper pressure fitting line, and the pressure value corresponding to the lower intersection point of the unconformity surface and the lower pressure fitting line can be further calculated, so as to better analyze the erosion amount of the unconformity surface.

[0074] Step S106, determining the formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line.

[0075] Specifically, pressure gradient refers to the change in pressure per unit length along the direction of fluid flow, and the commonly used value for unit depth is 100 meters. Formation pressure gradient refers to the increase in pressure per unit depth increase in the vertical depth of the formation from the ground. When applied to oilfield systems, it is often called pressure gradient, omitting the word "formation" in front, and the commonly used value for unit depth is 100 meters.

[0076] It is worth noting that the pressure gradient corresponds to the depth of the formation. The pressure gradient at different depths in the same geological structure is not the same. The formation pressure gradient can reflect the degree of change in the vertical pressure of the formation. The value of the formation pressure gradient of an unused reservoir is mainly affected by the density of the formation liquid. Therefore, the formation pressure gradient can be used to determine the oil-water interface. The pressure gradient of a used reservoir is the gradient generated by the injection-production pressure difference, so it can be determined whether it is affected by water injection.

[0077] Therefore, after obtaining the upper layer pressure fitting line and the lower layer pressure fitting line, the formation pressure gradient is determined based on the upper layer pressure fitting line and the lower layer pressure fitting line, so as to further analyze the erosion of the geological unconformity surface.

[0078] Step S107, determining the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient.

[0079] Specifically, from the above introduction, it can be seen that the present application can obtain the pressure value corresponding to the intersection of the unconformity surface and the upper pressure fitting line, and the pressure value corresponding to the intersection of the unconformity surface and the lower pressure fitting line. After obtaining the difference between the pressure values ​​of the two intersections, the erosion amount of the unconformity surface can be determined in combination with the formation pressure gradient.

[0080] It can be seen from the above technical scheme that the embodiment of the present application can determine the erosion amount of the unconformity surface based on the pressure values ​​of each sample collected from different geological depths where the unconformity surface exists, understand the erosion situation of the unconformity surface, and provide important reference parameters for geologists to conduct geological research and explore petroleum resources.

[0081] As can be seen from the introduction of the above technical solution, the present application can determine the pressure value and standard error of each sample based on the pressure data of each primary fluid inclusion of each sample. The following describes the process in detail, which may include:

[0082] The pressure average and standard error of the pressure data of each primary fluid inclusion in each sample are calculated as the pressure value and standard error of each sample.

[0083] Specifically, as can be seen from the above introduction, the pressure data of primary fluid inclusions in geology can reflect the erosion of unconformity surfaces in geology. Therefore, before analyzing the erosion amount of unconformity surfaces, the temperature of each primary fluid inclusion of each sample collected in the geological vertical section where the unconformity surface exists can be measured. After determining the temperature of each primary fluid inclusion of each sample, the temperature data of each sample can be converted into pressure data.

[0084] For example, the density method and the isochoric method can be used to convert the temperature data of each primary fluid inclusion in each sample into pressure data.

[0085] After determining the pressure data of each primary fluid inclusion in each sample, the pressure average value and standard error of the pressure data of each primary fluid inclusion in each sample can be further calculated as the pressure value and standard error of each sample.

[0086] It can be seen from the above technical solution that the embodiment of the present application can determine the pressure value and standard error of each sample based on the pressure data of each primary fluid inclusion of each sample, providing a reference for analyzing the erosion amount of the unconformity surface.

[0087] It can be seen from the above technical solution that the pressure data of each sample collected from the geology can reflect the erosion of the unconformity surface in the geology, but the pressure data of each sample is not easy to measure, and the embodiment of the present application can determine the temperature of each primary fluid inclusion of each sample collected from the geological vertical section where the unconformity surface exists, and thus, the temperature data of each primary fluid inclusion can be converted into pressure.

[0088] Next, the process of converting the temperature data of each primary fluid inclusion into pressure data is introduced. This process can include the following steps:

[0089] Step S11, determining the composition of each primary fluid inclusion in each sample.

[0090] From the above introduction, it can be seen that the temperature data of primary fluid inclusions can be converted into pressure data. For example, the temperature data of each primary fluid inclusion in each sample can be converted into pressure data using methods such as the density method and the isochoric method.

[0091] In particular, the conversion method for converting the temperature data of the primary fluid inclusion into the pressure data can be determined according to the composition of the primary fluid inclusion.

[0092] Specifically, the composition of each primary fluid inclusion in each sample can be determined so that a conversion method for converting temperature data of the primary fluid inclusions into pressure data can be determined.

[0093] Step S12, based on the composition of each primary fluid inclusion in each sample, determine a calculation method for converting the temperature data of each primary fluid inclusion in each sample into pressure data.

[0094] Specifically, from the above steps, it can be seen that the composition of each primary fluid inclusion in each sample can be determined. The calculation methods for converting primary fluid inclusions of different compositions into pressure data are different. Therefore, based on the composition of each primary fluid inclusion in each sample, the calculation method for converting the temperature data of each primary fluid inclusion in each sample into pressure data can be determined.

[0095] For example, if the composition of the primary fluid inclusion is determined to be a H2O-CO2 system, the pressure data of the primary fluid inclusion can be calculated based on the temperature of the primary fluid inclusion using the isochoric method. If the composition of the primary fluid inclusion is determined to be a H2O-NaCl system, the pressure data of the primary fluid inclusion can be calculated based on the temperature of the primary fluid inclusion using the density method or the isochoric method.

[0096] The composition of each primary fluid inclusion in each sample can be determined based on the temperature of each primary fluid inclusion in each sample when the temperature drops and a phase transition occurs, or the composition of the primary fluid inclusion can be directly measured by laser Raman spectroscopy.

[0097] It can be understood that, since the temperature data of the primary fluid inclusions can be converted into pressure data based on the composition of the primary fluid inclusions, the above step S102, converting the temperature data of each primary fluid inclusion in each sample into pressure data, can include:

[0098] Based on the calculation method corresponding to each primary fluid inclusion in each sample, the temperature data of each primary fluid inclusion in each sample is converted into pressure data.

[0099] For example, after determining the temperature data of each primary fluid inclusion in each sample, if it is determined that the composition of each primary fluid inclusion is a H2O-NaCl system, the following isochoric formula (1) can be used to calculate the pressure data of each primary fluid inclusion.

[0100] P = A1 + A2 × T1 (1);

[0101] in,

[0102] A1=6.1×10 -3 +(2.385×10 -1 -a1)×T h -(2.855×10 -1 -a2)×T h 2 -(a3×T+a4×T h 2 ) ×m;

[0103] A2=a1+a2×T h +9.888×10 -6 ×T h 2 +(a3+a4×T ) ×m;

[0104] In formula (1), P is pressure, unit is 10 -1 MPa;

[0105] T is the capture temperature in degrees Celsius;

[0106] a1, a2, a3, and a4 are constants. For the H2O-NaCl system of the primary fluid inclusion, a1 = 2.848 × 10 -1 ; a2=6.445×10 -2 ; a3=4.159×10 -1 ; a4=7.438×10 -3 ;

[0107] T h is the homogenization temperature of the primary fluid inclusion, in degrees Celsius;

[0108] m is the molar concentration of salts, in mol / kg. The conversion formula between m and salinity w is:

[0109]

[0110] It can be seen from the above technical solution that the embodiment of the present application can determine the calculation method for converting the temperature data of each primary fluid inclusion in each sample into pressure data according to the composition of each primary fluid inclusion in each sample, thereby realizing the conversion of the temperature data of each primary fluid inclusion in each sample into pressure data, so that the erosion of the unconformity surface can be determined according to the pressure data of each sample.

[0111] It can be seen from the above technical solution that the embodiment of the present application can determine the formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line.

[0112] It can be understood that the formation pressure gradient includes the overlying formation pressure gradient and the underlying formation pressure gradient.

[0113] Next, a process of determining the formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line is introduced, and the process may include the following steps:

[0114] Step S21, calculating the inverse of the slope of the upper layer pressure fitting line as the overlying formation pressure gradient.

[0115] Specifically, from the above introduction, it can be known that the upper layer pressure fitting line can be obtained based on the pressure data of each sample, and the inverse of the slope of the upper layer pressure fitting line can be further calculated as the overlying formation pressure gradient.

[0116] Step S22, calculating the inverse of the slope of the underlying formation pressure fitting line as the underlying formation pressure gradient.

[0117] From the above introduction, it can be known that the lower layer pressure fitting line can be obtained based on the pressure data of each sample, and the inverse of the slope of the lower layer pressure fitting line can be further calculated as the underlying formation pressure gradient.

[0118] It can be seen from the above technical solution that the embodiment of the present application can determine the formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line, so as to further determine the erosion amount of the unconformity surface.

[0119] From the above introduction, it can be seen that the embodiment of the present application can determine the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, as well as the formation pressure gradient. Therefore, the erosion amount of the unconformity surface can be determined based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient.

[0120] Next, a process of determining the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point and the formation pressure gradient is introduced. The process may include the following steps:

[0121] Step S31, determining the magnitude of the pressure value corresponding to the upper layer intersection point and the pressure value corresponding to the lower layer intersection point.

[0122] Specifically, after determining the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, as well as the formation pressure gradient, the magnitude of the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point can be further determined. If there is an overlying stratum on the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, step S32 can be executed; if the pressure value corresponding to the lower intersection point is less than the pressure value corresponding to the upper intersection point, step S33 can be executed; if there is no overlying stratum on the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, step S34 can be executed.

[0123] Step S32, determining the erosion thickness of the unconformity surface is equal to the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point divided by the pressure gradient of the underlying stratum.

[0124] Specifically, if there is an overlying stratum on the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, it means that the underlying stratum was once buried to a greater depth than the current burial depth of the underlying stratum, and therefore, the underlying stratum was eroded and buried until now. Assuming that the surface pressure is consistent from ancient times to the present, it can be determined that the erosion thickness of the unconformity surface is equal to the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point divided by the pressure gradient of the underlying stratum.

[0125] Step S33, determining the erosion thickness of the unconformity surface is equal to the difference between the pressure value corresponding to the lower layer intersection point and the surface pressure value divided by the overlying stratum pressure gradient.

[0126] Specifically, if the pressure value corresponding to the lower intersection point is less than the pressure value corresponding to the upper intersection point, it means that although the underlying stratum was buried very deeply before erosion, the burial depth of the underlying stratum before erosion did not reach a depth greater than the current burial depth of the underlying stratum, and the current burial depth of the underlying stratum is the maximum depth experienced by the underlying stratum. Then, assuming that the ancient stratum pressure gradient is equal to the current stratum pressure gradient, it can be determined that the erosion thickness of the unconformity surface is equal to the difference between the pressure value corresponding to the lower intersection point and the surface pressure value divided by the overlying stratum pressure gradient.

[0127] Among them, since the underlying stratum was further compacted during the later burial process, the thickness of the underlying stratum became smaller, and the fitted underlying stratum pressure gradient became larger than that before the actual erosion.

[0128] Step S34, determining the erosion thickness of the unconformity surface is equal to the difference between the pressure value corresponding to the lower layer intersection point and the surface pressure value divided by the pressure gradient of the underlying stratum.

[0129] Specifically, if there is no overlying formation above the unconformity surface and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, it is determined that the erosion thickness of the unconformity surface is equal to the difference between the pressure value corresponding to the lower intersection point and the surface pressure value divided by the pressure gradient of the underlying formation.

[0130] For example, there is an unconformity surface between the overlying formation and the underlying formation in a certain geology. A series of samples are collected in the upper geology and the lower geology located at the unconformity surface respectively, and the temperature values of the primary fluid inclusions of each obtained sample are measured. Then, the temperature data of the primary fluid inclusions of each sample are converted into pressure data. Next, linear fitting is performed on the pressure data of each sample, and an upper pressure fitting line and a lower pressure fitting line can be obtained.

[0131] As Figure 2 、 Figure 3 shown, Figure 2 、 Figure 3 is a fitting effect diagram of pressure data exemplified in an embodiment of the present application;

[0132] Among them, Figure 2 exemplifies the intersection situation between the upper pressure fitting line and the unconformity surface pressure fitting line, Figure 3 exemplifies the intersection situation between the lower pressure fitting line and the unconformity surface pressure fitting line.

[0133] After obtaining the upper pressure fitting line and the lower pressure fitting line, calculate the slope of each fitting line as the formation pressure gradient (unit: Pa / km).

[0134] Among them, the reciprocal of the slope of the upper pressure fitting line is the overlying formation pressure gradient G1, and the reciprocal of the slope of the lower pressure fitting line is the underlying formation pressure gradient G2. Further extend the upper pressure fitting line and the lower pressure fitting line to the unconformity surface, and two intersection points of the upper pressure fitting line, the lower pressure fitting line and the unconformity surface can be obtained. The intersection point coordinates of the upper pressure fitting line are (P1, h), and the intersection point coordinates of the lower pressure fitting line are (P2, h). The pressures of the upper intersection point and the lower intersection point correspond to P1 and P2 respectively. Among them, P S represents the surface pressure.

[0135] When P2 > P1, the erosion thickness of the unconformity surface = (P2 - Ps) / G2 (as Figure 2 shown).

[0136] When P2 < P1, assuming that the formation pressure gradient in the past is equal to the formation pressure gradient at the current moment, the erosion thickness of the unconformity surface = (P2 - Ps) / G1 (as Figure 3 shown).

[0137] In addition, the pressure gradient of the overlying strata and the pressure gradient of the underlying strata can be determined by determining the difference in buried depth and pressure value corresponding to any two points on the upper pressure fitting line or the lower pressure fitting line.

[0138] The calculation formula of formation pressure gradient is:

[0139]

[0140] Wherein, K is the formation pressure gradient, including the overlying formation pressure gradient and the underlying formation pressure gradient;

[0141] p1 can be the pressure value corresponding to any point on the upper layer pressure fitting line or the lower layer pressure fitting line, and p2 can be the pressure value corresponding to any other point on the upper layer pressure fitting line or the lower layer pressure fitting line;

[0142] h1 may be the burial depth corresponding to the pressure value p1 on the upper pressure fitting line or the lower pressure fitting line; h2 may be the burial depth corresponding to the pressure value p2 on the upper pressure fitting line or the lower pressure fitting line.

[0143] It can be seen from the above technical solution that the embodiment of the present application can determine the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient, thereby providing important reference parameters for geologists to conduct geological research and oil exploration.

[0144] The following is a description of an ablation amount determination device provided in an embodiment of the present application. The ablation amount determination device described below and the ablation amount determination method described above can refer to each other.

[0145] See also Figure 4 , Figure 4 This is a schematic diagram of the structure of an erosion amount determination device disclosed in an embodiment of the present application.

[0146] like Figure 4 As shown, the ablation amount determining device may include:

[0147] A temperature determination unit 101 is used to determine the temperature of each primary fluid inclusion in each sample collected from a geological vertical section where an unconformity exists;

[0148] A conversion unit 102, used to convert the temperature data of each primary fluid inclusion in each sample into pressure data;

[0149] A first calculation unit 103 is used to determine the pressure value and standard error of each sample based on the pressure data of each primary fluid inclusion in each sample;

[0150] The fitting unit 104 is used to perform linear fitting on the pressure data distribution of each sample of the upper geological layer and the lower geological layer located on the unconformity surface based on the pressure value and the standard error of each sample, so as to obtain an upper pressure fitting line and a lower pressure fitting line;

[0151] The second calculation unit 105 is used to calculate the pressure value corresponding to the intersection point between the unconformity surface and the upper layer of the upper layer pressure fitting line, and calculate the pressure value corresponding to the intersection point between the unconformity surface and the lower layer of the lower layer pressure fitting line;

[0152] A pressure gradient determination unit 106, configured to determine a formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line;

[0153] The erosion amount determination unit 107 is used to determine the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient.

[0154] It can be seen from the above technical solutions that the device of the embodiment of the present application can use the temperature determination unit 101 to determine the temperature of each primary fluid inclusion of each sample collected from the geological vertical section where the unconformity surface exists. Since the pressure data of each primary fluid inclusion in each sample can reflect the erosion amount of the unconformity surface, the temperature data of each primary fluid inclusion in each sample can be converted into pressure data by using the conversion unit 102; the first calculation unit 103 can be used to determine the pressure value and standard error of each sample based on the pressure data of each primary fluid inclusion in each sample; the fitting unit 104 can be used to respectively fit the upper geological layer and the lower geological layer located on the unconformity surface based on the pressure values ​​and standard errors of each sample. The pressure data distribution of each sample is linearly fitted to obtain an upper pressure fitting line and a lower pressure fitting line; after obtaining the upper pressure fitting line and the lower pressure fitting line, the second calculation unit 105 can be used to calculate the pressure value corresponding to the upper intersection of the unconformity surface and the upper pressure fitting line, and the pressure value corresponding to the lower intersection of the unconformity surface and the lower pressure fitting line; the pressure gradient determination unit 106 can be used to determine the formation pressure gradient based on the upper pressure fitting line and the lower pressure fitting line; finally, the erosion amount determination unit 107 can be used to determine the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection and the pressure value corresponding to the lower intersection, and the formation pressure gradient. The present application can determine the erosion amount of the unconformity surface based on the pressure data of each primary fluid inclusion of each sample collected on the unconformity surface, providing important reference parameters for geologists to carry out geological research and explore for oil.

[0155] Further optionally, the execution process of the first calculation unit 103 may include:

[0156] The pressure average and standard error of the pressure data of each primary fluid inclusion in each sample are calculated as the pressure value and standard error of each sample.

[0157] Further optionally, the above-mentioned device for determining the amount of ablation may further include:

[0158] A composition determination unit, used to determine the composition of each primary fluid inclusion in each sample;

[0159] A calculation method determination unit, for determining a calculation method for converting temperature data of each primary fluid inclusion in each sample into pressure data based on the composition of each primary fluid inclusion in each sample;

[0160] Based on this, the execution process of the above conversion unit 102 may include:

[0161] Based on the calculation method corresponding to each primary fluid inclusion in each sample, the temperature data of each primary fluid inclusion in each sample is converted into pressure data.

[0162] Further optionally, the execution process of the pressure gradient determining unit 106 may include:

[0163] Calculating the inverse of the slope of the upper layer pressure fitting line as the overlying formation pressure gradient;

[0164] The inverse of the slope of the underlying formation pressure fitting line is calculated as the underlying formation pressure gradient.

[0165] Further optionally, the execution process of the above-mentioned erosion amount determination unit 107 may include:

[0166] Determine the magnitude of the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point;

[0167] If there is an overlying stratum on the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, then the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point divided by the pressure gradient of the underlying stratum;

[0168] If the pressure value corresponding to the lower intersection point is less than the pressure value corresponding to the upper intersection point, the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the lower intersection point and the surface pressure value divided by the pressure gradient of the overlying stratum;

[0169] If there is no overlying stratum above the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, then the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the lower intersection point and the surface pressure value divided by the pressure gradient of the underlying stratum.

[0170] Among them, the specific processing flow of each unit included in the above-mentioned erosion amount determination device can refer to the relevant introduction of the erosion amount determination method part above, and will not be repeated here.

[0171] The apparatus for determining the amount of ablation provided in the embodiment of the present application can be applied to an apparatus for determining the amount of ablation, such as a terminal: a mobile phone, a computer, etc. Optionally, Figure 5 A hardware structure diagram of a device for determining the amount of erosion is shown. Figure 5 The hardware structure of the device for determining the amount of erosion may include: at least one processor 1, at least one communication interface 2, at least one memory 3 and at least one communication bus 4.

[0172] In the embodiment of the present application, the number of the processor 1 , the communication interface 2 , the memory 3 , and the communication bus 4 is at least one, and the processor 1 , the communication interface 2 , and the memory 3 communicate with each other through the communication bus 4 .

[0173] The processor 1 may be a central processing unit CPU, or an application-specific integrated circuit ASIC (Application Specific Integrated Circuit), or one or more integrated circuits configured to implement the embodiments of the present application, etc.;

[0174] The memory 3 may include a high-speed RAM memory, and may also include a non-volatile memory, such as at least one disk memory;

[0175] The memory stores a program, and the processor can call the program stored in the memory, wherein the program is used to implement various processing flows of the aforementioned terminal in the ablation amount determination scheme.

[0176] The embodiment of the present application further provides a readable storage medium, which can store a program suitable for execution by a processor, wherein the program is used to implement various processing flows of the aforementioned terminal in the ablation amount determination scheme.

[0177] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0178] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0179] The above description of the disclosed embodiments enables professionals and technicians in the field to implement or use the present application. Various modifications to these embodiments will be apparent to professionals and technicians in the field, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present application. The various embodiments can be combined with each other. Therefore, the present application will not be limited to the embodiments shown herein, but will conform to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for determining the amount of ablation, characterized in that: include: Determine the temperature of each primary fluid inclusion in each sample collected from a geologic vertical section where the unconformity exists; Convert the temperature data of each primary fluid inclusion in each sample into pressure data; Based on the pressure data of each primary fluid inclusion in each sample, determine the pressure value and standard error of each sample; Based on the pressure value and standard error of each sample, linear fitting is performed on the pressure data distribution of each sample of the upper geological layer and the lower geological layer located on the unconformity surface to obtain an upper pressure fitting line and a lower pressure fitting line; Calculating the pressure value corresponding to the upper layer intersection point of the unconformity surface and the upper layer pressure fitting line, and calculating the pressure value corresponding to the lower layer intersection point of the unconformity surface and the lower layer pressure fitting line; Determining a formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line; Determining the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient; Before converting the temperature data of each primary fluid inclusion in each sample into pressure data, the method further includes: Determine the composition of individual primary fluid inclusions in each sample; Determining a computational method for converting temperature data of each primary fluid inclusion in each sample into pressure data based on the composition of each primary fluid inclusion in each sample; Based on this, the temperature data of each primary fluid inclusion in each sample is converted into pressure data, including: Based on a calculation method corresponding to each primary fluid inclusion in each sample, temperature data of each primary fluid inclusion in each sample is converted into pressure data; The formation pressure gradient includes an overlying formation pressure gradient and an underlying formation pressure gradient, and determining the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper layer intersection point and the pressure value corresponding to the lower layer intersection point and the formation pressure gradient includes: Determine the magnitude of the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point; If there is an overlying stratum on the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, then the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point divided by the pressure gradient of the underlying stratum; If the pressure value corresponding to the lower intersection point is less than the pressure value corresponding to the upper intersection point, then the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the lower intersection point and the surface pressure value divided by the pressure gradient of the overlying stratum; The method of determining the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point and the formation pressure gradient further includes: If there is no overlying stratum above the unconformity surface, and the pressure value corresponding to the lower intersection point is greater than the pressure value corresponding to the upper intersection point, then the erosion thickness of the unconformity surface is determined to be equal to the difference between the pressure value corresponding to the lower intersection point and the surface pressure value divided by the pressure gradient of the underlying stratum.

2. The method for determining the amount of ablation according to claim 1, characterized in that: The method of determining the pressure value and standard error of each sample based on the pressure data of each primary fluid inclusion in each sample includes: The pressure average and standard error of the pressure data of each primary fluid inclusion in each sample are calculated as the pressure value and standard error of each sample.

3. The method for determining the amount of ablation according to claim 1, characterized in that: The formation pressure gradient includes an overlying formation pressure gradient and an underlying formation pressure gradient, and determining the formation pressure gradient based on the overlying formation pressure fitting line and the underlying formation pressure fitting line includes: Calculating the inverse of the slope of the upper layer pressure fitting line as the overlying formation pressure gradient; The inverse of the slope of the underlying formation pressure fitting line is calculated as the underlying formation pressure gradient.

4. A device for determining the amount of ablation, characterized in that: The method for determining the amount of ablation as claimed in claim 1 comprises: a temperature determination unit for determining the temperature of each primary fluid inclusion in each sample collected from a geological vertical section where an unconformity exists; A conversion unit, used to convert the temperature data of each primary fluid inclusion in each sample into pressure data; A first calculation unit is used to determine the pressure value and standard error of each sample based on the pressure data of each primary fluid inclusion in each sample; A fitting unit is used to perform linear fitting on the pressure data distribution of each sample of the upper geological layer and the lower geological layer located on the unconformity surface based on the pressure value and standard error of each sample, so as to obtain an upper pressure fitting line and a lower pressure fitting line; A second calculation unit is used to calculate the pressure value corresponding to the intersection point between the unconformity surface and the upper layer of the upper layer pressure fitting line, and calculate the pressure value corresponding to the intersection point between the unconformity surface and the lower layer of the lower layer pressure fitting line; a pressure gradient determining unit, configured to determine a formation pressure gradient based on the upper layer pressure fitting line and the lower layer pressure fitting line; The erosion amount determination unit is used to determine the erosion amount of the unconformity surface based on the difference between the pressure value corresponding to the upper intersection point and the pressure value corresponding to the lower intersection point, and the formation pressure gradient.

5. The device for determining the amount of ablation according to claim 4, characterized in that: The execution process of the first computing unit includes: The pressure average and standard error of the pressure data of each primary fluid inclusion in each sample are calculated as the pressure value and standard error of each sample.

6. An apparatus for determining the amount of ablation, characterized in that: include: one or more processors, and memory; The memory stores computer-readable instructions, which, when executed by the one or more processors, implement the steps of the method for determining the amount of ablation as claimed in any one of claims 1 to 3.

7. A readable storage medium, characterized in that: The readable storage medium stores computer-readable instructions, and when the computer-readable instructions are executed by one or more processors, the one or more processors implement the steps of the method for determining the amount of ablation as claimed in any one of claims 1 to 3.

Images