Method, system and device for calculating angular unconformity erosion thickness of sedimentary basins

Abstract

The present invention relates to a method, system and device for calculating the angular unconformity denudation thickness of a sedimentary basin, the method comprising obtaining N samples; testing the N samples with fluid inclusion homogenization temperature to obtain N fluid inclusion temperature data and surface temperature; calculating the N fluid inclusion temperature data and surface temperature with fluid inclusion pressure to obtain N pressure data and surface pressure; correcting the homogenization temperature of the N fluid inclusion temperature data and the N pressure data with layer-by-layer projection and linear fitting to obtain temperature gradient correction data and pressure gradient correction data; and calculating the denudation thickness of each sample in the sedimentary basin with an erosion thickness calculation formula. The method obtains the denudation thickness of each sample in the sedimentary basin with the temperature gradient correction data or the pressure gradient correction data, thereby obtaining the vertical denudation thickness of the underlying strata of the corrected angular unconformity.

Description

Technical Field

[0001] The present invention relates to the field of geological technology, and in particular to a method, system and equipment for calculating the angular unconformity erosion thickness of a sedimentary basin. Background Art

[0002] Restoration of the denudation thickness of sedimentary basins has always been an issue that the geological community and the energy and mineral exploration community urgently need to explore. Restoration of the denudation thickness of sedimentary basins plays a vital role in the restoration of geological history, hydrocarbon accumulation, and mineral preservation. There are many technical details involved, and the improvement of each technical detail will promote the further application of this problem in different regions. However, the existing data for calculating the denudation thickness of sedimentary basins with angular unconformity is inaccurate, which affects the restoration of geological history, hydrocarbon accumulation, and mineral preservation. Summary of the invention

[0003] The embodiment of the present invention provides a method, system and device for calculating the angular unconformity denudation thickness of a sedimentary basin, which are used to solve the technical problem of inaccurate data in calculating the denudation thickness of a sedimentary basin with an angular unconformity.

[0004] In order to achieve the above purpose, the embodiment of the present invention provides the following technical solutions:

[0005] A method for calculating the angular unconformity erosion thickness of a sedimentary basin comprises the following steps:

[0006] Obtain the position of the angular unconformity of the sedimentary basin to be calculated for the denudation thickness and obtain N samples along the direction of the vertical strata of the sedimentary basin at the position of the angular unconformity;

[0007] Performing a fluid inclusion uniform temperature test on N samples to obtain N fluid inclusion temperature data and a surface temperature; performing a fluid inclusion pressure calculation on the N fluid inclusion temperature data and a surface temperature to obtain N pressure data and a surface pressure;

[0008] The uniform temperature of the N fluid inclusion temperature data and the N pressure data are respectively processed by layer-wise projection correction and linear fitting to obtain temperature gradient correction data and pressure gradient correction data;

[0009] The uniform temperature of the fluid inclusion temperature data and the temperature gradient correction data are calculated by using an erosion thickness calculation formula, or the pressure data and the pressure gradient correction data are calculated to obtain the erosion thickness of each sample in the sedimentary basin;

[0010] The fluid inclusion temperature data include the uniform temperature, initial melting temperature, freezing point temperature and composition of each sample.

[0011] Preferably, the erosion thickness calculation formula is: E=ΔT / G or ΔT / G-Δy, E is the erosion thickness of each sample in the sedimentary basin, ΔT is the difference between the homogenization temperature of the fluid inclusion temperature data of the sample and the surface temperature or the difference between the pressure data and the surface pressure, G is the temperature gradient correction data or the pressure gradient correction data, and Δy is the vertical distance between the sample and the angular unconformity.

[0012] Preferably, before the fluid inclusion homogenization temperature test is performed on the N samples, the following conditions are included: the position of the angular unconformity is that the angular unconformity is located in the overlying strata of the sedimentary basin, or the position of the angular unconformity is that the angular unconformity is located in the underlying strata of the sedimentary basin and the former maximum burial depth of the underlying strata is greater than the current burial depth.

[0013] Preferably, the method for calculating the angular unconformity denudation thickness of the sedimentary basin comprises: obtaining the former maximum burial depth of the underlying strata of the sedimentary basin by comparing vitrinite reflectance and regional strata.

[0014] Preferably, the method for calculating the angular unconformity denudation thickness of the sedimentary basin comprises: sampling the sedimentary basin at the position of the angular unconformity along the direction of the vertical strata strike of the sedimentary basin using surface or shallow well sampling to obtain N samples.

[0015] The present invention also provides a system for calculating the angular unconformity erosion thickness of a sedimentary basin, comprising a parameter acquisition module, a parameter processing module, a data processing module and a calculation module;

[0016] The parameter acquisition module is used to obtain the position of the angular unconformity of the sedimentary basin to be calculated for the denudation thickness and to obtain N samples along the direction of the vertical stratum strike of the sedimentary basin at the position of the angular unconformity;

[0017] The parameter processing module is used to perform fluid inclusion uniform temperature testing on N samples to obtain N fluid inclusion temperature data and a surface temperature; and to perform fluid inclusion pressure calculation on the N fluid inclusion temperature data and one surface temperature to obtain N pressure data and a surface pressure;

[0018] The data processing module is used to process the uniform temperature of the N fluid inclusion temperature data and the N pressure data by layer-by-layer projection correction and linear fitting to obtain temperature gradient correction data and pressure gradient correction data;

[0019] The calculation module is used to calculate the homogenization temperature of the fluid inclusion temperature data and the temperature gradient correction data, or to calculate the pressure data and the pressure gradient correction data, by using the denudation thickness calculation formula, to obtain the denudation thickness of each sample in the sedimentary basin;

[0020] The fluid inclusion temperature data include the uniform temperature, initial melting temperature, freezing point temperature and composition of each sample.

[0021] Preferably, the erosion thickness calculation formula is: E=ΔT / G or ΔT / G-Δy, E is the erosion thickness of each sample in the sedimentary basin, ΔT is the difference between the homogenization temperature of the fluid inclusion temperature data of the sample and the surface temperature or the difference between the pressure data and the surface pressure, G is the temperature gradient correction data or the pressure gradient correction data, and Δy is the vertical distance between the sample and the angular unconformity.

[0022] Preferably, the parameter acquisition module is used to perform fluid inclusion homogenization temperature tests on N samples according to the position of the angular unconformity being the angular unconformity located in the overlying strata of the sedimentary basin, or the position of the angular unconformity being the angular unconformity located in the underlying strata of the sedimentary basin and the former maximum burial depth of the underlying strata being greater than the current burial depth.

[0023] Preferably, the parameter acquisition module is used to obtain N samples by surface or shallow well sampling of the sedimentary basin along the direction of the vertical strata strike of the sedimentary basin at the position of the angular unconformity.

[0024] The present invention also provides a device for calculating the angular unconformity erosion thickness of a sedimentary basin, comprising a processor and a memory;

[0025] The memory is used to store program codes and transmit the program codes to the processor;

[0026] The processor is used to execute the above-mentioned method for calculating the angular unconformity erosion thickness of a sedimentary basin according to the instructions in the program code.

[0027] It can be seen from the above technical scheme that the embodiments of the present invention have the following advantages: the method, system and equipment for calculating the angular unconformity erosion thickness of the sedimentary basin, the method comprising: obtaining the position of the angular unconformity of the sedimentary basin to be calculated for the erosion thickness and obtaining N samples at the position of the angular unconformity along the direction of the vertical strata of the sedimentary basin; using fluid inclusion homogenization temperature test on the N samples respectively to obtain N fluid inclusion temperature data and a surface temperature; using fluid inclusion pressure calculation on the N fluid inclusion temperature data and a surface temperature to obtain N pressure data and a surface pressure; using layer-by-layer projection correction and linear fitting processing on the homogenization temperature of the N fluid inclusion temperature data and the N pressure data to obtain temperature gradient correction data and pressure gradient correction data; calculating the homogenization temperature and temperature gradient correction data of the fluid inclusion temperature data or calculating the pressure data and pressure gradient correction data through the erosion thickness calculation formula to obtain the erosion thickness of each sample in the sedimentary basin. The method for calculating the angular unconformity erosion thickness of a sedimentary basin obtains the erosion thickness of each sample in the sedimentary basin through temperature gradient correction data or pressure gradient correction data, thereby achieving the vertical erosion thickness of the underlying strata of the corrected angular unconformity, and solving the technical problem of inaccurate data in calculating the erosion thickness of sedimentary basins with angular unconformity. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] In order to more clearly illustrate the embodiments of the present invention 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 invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0029] Figure 1 A flowchart of the steps of the method for calculating the angular unconformity erosion thickness of a sedimentary basin according to an embodiment of the present invention;

[0030] Figure 2 It is a schematic diagram of the angular unconformity in the method for calculating the angular unconformity erosion thickness of a sedimentary basin according to an embodiment of the present invention;

[0031] Figure 3 A schematic diagram of a layer-by-layer projection of a method for calculating the angular unconformity erosion thickness of a sedimentary basin according to another embodiment of the present invention;

[0032] Figure 4 A schematic diagram of the erosion amount restored by the vitrinite reflectivity of a vertical section of a method for calculating the angular unconformity erosion thickness of a sedimentary basin according to an embodiment of the present invention;

[0033] Figure 5It is a framework diagram of a system for calculating the angular unconformity erosion thickness of a sedimentary basin according to an embodiment of the present invention. DETAILED DESCRIPTION

[0034] In order to make the purpose, features and advantages of the present invention more obvious and easy to understand, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the embodiments described below are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0035] The embodiments of the present application provide a method, system and device for calculating the angular unconformity denudation thickness of a sedimentary basin, which are used to solve the technical problem of inaccurate data in calculating the denudation thickness of a sedimentary basin with an angular unconformity.

[0036] Embodiment 1:

[0037] Figure 1 This is a flowchart of the steps of the method for calculating the angular unconformity erosion thickness of a sedimentary basin according to an embodiment of the present invention. Figure 2 It is a schematic diagram of an angular unconformity in the method for calculating the angular unconformity erosion thickness of a sedimentary basin according to an embodiment of the present invention.

[0038] like Figure 1 As shown, an embodiment of the present invention provides a method for calculating the angular unconformity erosion thickness of a sedimentary basin, comprising the following steps:

[0039] S1. Obtain the position of the angular unconformity of the sedimentary basin to be calculated for the erosion thickness and obtain N samples along the direction of the vertical stratigraphic strike of the sedimentary basin at the position of the angular unconformity.

[0040] It should be noted that the main purpose is to obtain the angular unconformity position of the sedimentary basin where the erosion thickness is to be calculated and to obtain samples along the direction perpendicular to the stratigraphic strike of the sedimentary basin.

[0041] In the embodiment of the present invention, Figure 2As shown in the figure, angular unconformity refers to a contact relationship between strata. When the upper and lower strata (i.e., overlying strata and underlying strata) have different strikes and intersect at a certain angle, it means that the overlying strata and the underlying strata are discontinuous in age, and there is a sedimentary hiatus between the two. During the sedimentary hiatus, they were uplifted and eroded by strong tectonic movements, forming an erosion surface. Later, the erosion surface dropped and accepted sedimentation again, forming a new stratum; and the "crack" between the two is the angular unconformity. Angular unconformity can be discovered through field outcrop surveys, seismic profile interpretation and other means. Of course, the difference between angular unconformity and parallel unconformity is that there is an angular intersection between the overlying strata and the underlying strata, and the underlying strata are tilted, while parallel unconformity has no angular intersection, which means that the underlying strata did not tilt before the unconformity was formed.

[0042] In the embodiment of the present invention, N samples are obtained by surface or shallow well sampling in the sedimentary basin along the direction perpendicular to the strata strike of the sedimentary basin at the position of the angular unconformity.

[0043] It should be noted that when the angular unconformity of a sedimentary basin is directly exposed on the surface, sampling along the surface outcrop can be used. When there is an overlying stratum above the angular unconformity of a sedimentary basin, shallow well coring can be used. Of course, surface sampling and shallow well coring can also be combined, thereby obtaining a series of rock samples collected under the angular unconformity and along the underlying stratum in the direction perpendicular to the stratum strike, that is, N samples.

[0044] S2. Perform fluid inclusion homogenization temperature test on N samples respectively to obtain N fluid inclusion temperature data and a surface temperature; perform fluid inclusion pressure calculation on the N fluid inclusion temperature data and a surface temperature to obtain N pressure data and a surface pressure. The fluid inclusion temperature data include the homogenization temperature, initial melting temperature, freezing point temperature and composition of each sample.

[0045] It should be noted that the fluid inclusion homogenization temperature test is mainly used for each sample to obtain N fluid inclusion temperature data corresponding to N samples, and the fluid inclusion pressure is calculated for each fluid inclusion temperature data and a surface temperature to obtain N pressure data and a surface pressure corresponding to N samples. In this embodiment, regardless of whether the fluid inclusion homogenization temperature test or the fluid inclusion pressure calculation is selected for each sample, N fluid inclusion temperature data and a surface temperature or N pressure data and a surface pressure required for step S3 are obtained. In the test process, the rock (i.e., the sample) is first ground into thin slices, then observed under a microscope, the primary fluid inclusions are selected for marking, and then the components are tested by Raman analysis, and then the primary fluid inclusions are subjected to a homogenization temperature test on a hot and cold table or the fluid inclusion pressure is calculated based on the fluid inclusion temperature data obtained from the inclusions. In the test, the primary fluid inclusions are subjected to a homogenization temperature test or the fluid inclusion pressure is calculated based on the fluid inclusion temperature data obtained from the inclusions. It is a common test technology for geological testing and will not be described in detail here. In the test, the components of each sample may include carbon dioxide, sodium chloride, etc.

[0046] S3. The uniform temperature of N fluid inclusion temperature data and N pressure data are processed by layer-wise projection correction and linear fitting to obtain temperature gradient correction data and pressure gradient correction data.

[0047] It should be noted that the uniform temperature of the N fluid inclusion temperature data and the N pressure data obtained in step S2 are mainly processed to obtain the temperature gradient correction data and pressure gradient correction data required for step S4.

[0048] In the embodiment of the present invention, the uniform temperature of N fluid inclusion temperature data and N pressure data are corrected by layer projection and linear fitting to obtain temperature gradient correction data and pressure gradient correction data. N is a natural number greater than 1.

[0049] It should be noted that bedding projection correction refers to transforming the acquired data through projection to ensure that the data is obtained under the same projection, reduce data errors, and achieve projection correction. Linear fitting processing refers to a form of curve fitting. Assume that x and Δy are both observed quantities, and Δy is a function of x: Δy = f(x; b). Curve fitting is to find the best estimate of parameter b through the observed values ​​of x and Δy, and to find the best theoretical curve Δy = f(x; b). Bedding projection correction and linear fitting processing are relatively mature technologies. How to use bedding projection correction and linear fitting processing to obtain temperature gradient correction data and pressure gradient correction data for the uniform temperature of N fluid inclusion temperature data and N pressure data will not be described in detail.

[0050] S4. The uniform temperature of the fluid inclusion temperature data and the temperature gradient correction data are calculated by the erosion thickness calculation formula, or the pressure data and the pressure gradient correction data are calculated to obtain the erosion thickness of each sample in the sedimentary basin.

[0051] It should be noted that the temperature or pressure is mainly selected to calculate the erosion thickness of each sample in the sedimentary basin through the erosion thickness calculation formula.

[0052] In an embodiment of the present invention, the calculation formula for the erosion thickness is: E=ΔT / G or ΔT / G-Δy, E is the erosion thickness of each sample in the sedimentary basin, ΔT is the difference between the homogenization temperature of the fluid inclusion temperature data of the sample and the surface temperature or the difference between the pressure data and the surface pressure, G is the temperature gradient correction data or the pressure gradient correction data, and Δy is the vertical distance between the sample and the angular unconformity.

[0053] The present invention provides a method for calculating the angular unconformity denudation thickness of a sedimentary basin, comprising: obtaining the position of the angular unconformity of the sedimentary basin whose erosion thickness is to be calculated and obtaining N samples along the direction of the vertical strata strike of the sedimentary basin at the position of the angular unconformity; performing fluid inclusion homogenization temperature test on the N samples respectively to obtain N fluid inclusion temperature data and a surface temperature; performing fluid inclusion pressure calculation on the N fluid inclusion temperature data and a surface temperature to obtain N pressure data and a surface pressure; performing layer-parallel projection correction and linear fitting processing on the homogenization temperature of the N fluid inclusion temperature data and the N pressure data to obtain temperature gradient correction data and pressure gradient correction data; calculating the homogenization temperature and temperature gradient correction data of the fluid inclusion temperature data or calculating the pressure data and the pressure gradient correction data by using an erosion thickness calculation formula to obtain the denudation thickness of each sample in the sedimentary basin. The method for calculating the angular unconformity erosion thickness of a sedimentary basin obtains the erosion thickness of each sample in the sedimentary basin through temperature gradient correction data or pressure gradient correction data, thereby achieving the vertical erosion thickness of the underlying strata of the corrected angular unconformity, and solving the technical problem of inaccurate data in calculating the erosion thickness of sedimentary basins with angular unconformity.

[0054] It should be noted that the method for calculating the angular unconformity erosion thickness of the sedimentary basin reduces the difficulty of sample sampling by obtaining N samples through surface or shallow well sampling in the sedimentary basin. Compared with the existing method that requires obtaining a complete drilling vertical profile at each point for the vertical erosion thickness, the method for calculating the angular unconformity erosion thickness of the sedimentary basin can complete the erosion thickness recovery work by scattered surface or shallow well sampling, thereby improving work efficiency.

[0055] When the underlying strata are not monoclinic strata but flexural folds, it is impossible to use surface or shallow well sampling in the sedimentary basin along the direction of the vertical strata of the sedimentary basin at the position of the angular unconformity to obtain N samples. Since the flexural folds may be irregular in shape and unevenly distributed after rock deformation, the amount of denudation cannot be deduced along the layer. The angular unconformity denudation thickness calculation method of the sedimentary basin is only for reference. However, if the sampling area is locally single-tilted, the angular unconformity denudation thickness calculation method of the sedimentary basin is still suitable in this area. The angular unconformity denudation thickness calculation method of the sedimentary basin can obtain the sedimentary basin with angular unconformity. Through the correction of the depth along the layer, the differential denudation thickness of each point on the profile or even the plane can be restored.

[0056] In one embodiment of the present invention, before N samples are subjected to fluid inclusion homogenization temperature test or pressure test, the following conditions are included: the position of the angular unconformity is that the angular unconformity is located in the overlying strata of the sedimentary basin, or the position of the angular unconformity is that the angular unconformity is located in the underlying strata of the sedimentary basin and the former maximum burial depth of the underlying strata is greater than the current burial depth.

[0057] It should be noted that the above is mainly to ensure that the sample is suitable for the angular unconformity erosion thickness calculation method of the sedimentary basin, and to ensure the accuracy of the erosion thickness of each sample in the sedimentary basin.

[0058] In one embodiment of the present invention, the method for calculating the angular unconformity denudation thickness of the sedimentary basin includes: obtaining the former maximum burial depth of the underlying strata of the sedimentary basin by comparing vitrinite reflectance and regional strata.

[0059] It should be noted that if there is an overlying stratum above the position of the angular unconformity, then it can be determined that the maximum burial depth of the underlying stratum was greater than the present burial depth by testing the homogenized temperature (or pressure) of the fluid inclusions in rock samples collected from the overlying strata and the underlying strata adjacent to each other above and below the angular unconformity. When the homogenized temperature (or pressure) of the overlying strata is greater than, or relatively close to, or almost equal to the homogenized temperature (or pressure) of the middle and underlying strata under the angular unconformity, it can be determined that the maximum burial depth of the underlying strata is most likely not greater than the present burial depth. It is also possible to indirectly determine whether the maximum burial depth was greater than the present burial depth through other methods such as vitrinite reflectance and regional stratigraphic comparison.

[0060] Figure 3 It is a schematic diagram of the layer-by-layer projection of a method for calculating the angular unconformity erosion thickness of a sedimentary basin according to another embodiment of the present invention.

[0061] like Figure 3As shown, it should be noted that the parameters obtained during projection are as follows: the horizontal distance of the layer projection = the vertical distance from the sample to the intersection of the ground layer and the surface, and the vertical distance of the layer projection = the vertical distance from the sample to the intersection of the ground layer and the surface * the tangent value of the inclination angle α. Figure 3 ). The horizontal distance of the layer projection Δd = the vertical distance from the sample to the intersection of the stratigraphic layer and the surface, and the vertical distance of the layer projection Δh = the vertical distance from the sample to the intersection of the stratigraphic layer. Δh = Δd*tgα. Stars A and B are sampling points; circle C is the projection point of sample B directly below sample A, and α is the inclination of the underlying stratum.

[0062] Figure 4 It is a schematic diagram of the erosion amount restored by the vertical section vitrinite reflectivity of the method for calculating the angular unconformity erosion thickness of a sedimentary basin according to an embodiment of the present invention.

[0063] In the embodiment of the present invention, Figure 4 As shown, Figure 4 The star on the left represents the sample. Figure 4 The circle on the left represents the projection point. Figure 4 The circle on the right in the middle is the homogenization temperature test result of the fluid inclusion; T represents the homogenization temperature of the fluid inclusion temperature data, Ts represents the surface temperature, and (T, Δy) is the temperature and depth coordinates. Figure 4 The two figures above show that the angular unconformity surface is directly exposed on the surface, and sampling can be carried out on the surface outcrop, which is relatively easy. Then the sample at the top of the underlying stratum is taken as the top sample ( Figure 4 The leftmost star in the upper left figure), other samples are projected to the same vertical direction along the formation dip, and then the uniform temperature or pressure data of the fluid inclusion temperature data is obtained by homogenizing the fluid inclusions for linear fitting, and the temperature gradient correction data or pressure gradient correction data (that is, correction data) G is obtained, then the denudation thickness of each sample point = (uniform temperature of fluid inclusion temperature data - surface temperature) / G. Similarly, the fluid inclusion temperature data is replaced by pressure data. Figure 4 The two figures below show that there is a shallow layer covering the angular unconformity, which can be sampled through shallow wells ( Figure 4 The difficulty of the work is still lower than that of coring in a vertical deep well system. Then the sample at the top of the underlying stratum is taken as the top sample ( Figure 4The leftmost star in the lower left figure), other samples are projected to the same vertical direction along the formation dip; then the uniform temperature or pressure data of the fluid inclusion temperature data is obtained by homogenizing the fluid inclusions, and linear fitting is performed to obtain the temperature gradient correction data or pressure gradient correction data (that is, correction data) G, then the denudation thickness of each sample point = (uniform temperature of fluid inclusion temperature data - surface temperature) / G-Δy. Similarly, the uniform fluid inclusion temperature data is replaced with pressure data.

[0064] Embodiment 2:

[0065] Figure 5 It is a framework diagram of a system for calculating the angular unconformity erosion thickness of a sedimentary basin according to an embodiment of the present invention.

[0066] like Figure 5 As shown, an embodiment of the present invention further provides a system for calculating the angular unconformity erosion thickness of a sedimentary basin, comprising a parameter acquisition module 10, a parameter processing module 20, a data processing module 30 and a calculation module 40;

[0067] The parameter acquisition module 10 is used to obtain the position of the angular unconformity of the sedimentary basin to be calculated for the denudation thickness and to obtain N samples along the vertical stratum strike direction of the sedimentary basin at the position of the angular unconformity;

[0068] The parameter processing module 20 is used to perform fluid inclusion homogenization temperature tests on N samples respectively to obtain N fluid inclusion temperature data and a surface temperature; and to perform fluid inclusion pressure calculation on the N fluid inclusion temperature data and a surface temperature to obtain N pressure data and a surface pressure;

[0069] The data processing module 30 is used to process the uniform temperature of the N fluid inclusion temperature data and the N pressure data by layer-by-layer projection correction and linear fitting to obtain temperature gradient correction data and pressure gradient correction data;

[0070] The calculation module 40 is used to calculate the uniform temperature and temperature gradient correction data of the fluid inclusion temperature data, or to calculate the pressure data and pressure gradient correction data through the erosion thickness calculation formula to obtain the erosion thickness of each sample in the sedimentary basin.

[0071] In an embodiment of the present invention, the calculation formula for the erosion thickness is: E=ΔT / G or ΔT / G-Δy, E is the erosion thickness of each sample in the sedimentary basin, ΔT is the difference between the homogenization temperature of the fluid inclusion temperature data of the sample and the surface temperature or the difference between the pressure data and the surface pressure, G is the temperature gradient correction data or the pressure gradient correction data, and Δy is the vertical distance between the sample and the angular unconformity.

[0072] In an embodiment of the present invention, the parameter acquisition module is used to perform fluid inclusion uniform temperature tests on N samples according to the position of the angular unconformity being the overlying strata in the sedimentary basin, or the position of the angular unconformity being the underlying strata in the sedimentary basin and the former maximum burial depth of the underlying strata being greater than the current burial depth.

[0073] In the embodiment of the present invention, the parameter acquisition module 10 is used to obtain N samples by sampling the sedimentary basin at the position of the angular unconformity along the direction of the vertical strata of the sedimentary basin using surface or shallow well sampling.

[0074] It should be noted that the modules in the system of the second embodiment correspond to the steps in the method of the first embodiment. The steps of the method of the first embodiment have been described in detail in the first embodiment, and the contents of the modules will not be described in detail in the second embodiment.

[0075] Embodiment three:

[0076] An embodiment of the present invention provides a device for calculating the angular unconformity denudation thickness of a sedimentary basin, including a processor and a memory;

[0077] A memory, used for storing program codes and transmitting the program codes to a processor;

[0078] The processor is used to execute the above-mentioned method for calculating the angular unconformity denudation thickness of the sedimentary basin according to the instructions in the program code.

[0079] It should be noted that the processor is used to execute the steps in the above-mentioned method for calculating the angular unconformity denudation thickness of a sedimentary basin according to the instructions in the program code. Alternatively, the processor implements the functions of each module / unit in the above-mentioned system / device embodiments when executing the computer program.

[0080] Exemplarily, the computer program may be divided into one or more modules / units, one or more modules / units are stored in a memory and executed by a processor to complete the present application. One or more modules / units may be a series of computer program instruction segments capable of completing a specific function, and the instruction segments are used to describe the execution process of the computer program in a terminal device.

[0081] The terminal device may be a computing device such as a desktop computer, a notebook, a PDA, or a cloud server. The terminal device may include, but is not limited to, a processor and a memory. Those skilled in the art will appreciate that this does not constitute a limitation on the terminal device, and may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the terminal device may also include input and output devices, network access devices, buses, etc.

[0082] The processor may be a central processing unit (CPU), other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

[0083] The memory may be an internal storage unit of the terminal device, such as a hard disk or memory of the terminal device. The memory may also be an external storage device of the terminal device, such as a plug-in hard disk equipped on the terminal device, a smart memory card (SmartMedia Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), etc. Further, the memory may also include both an internal storage unit of the terminal device and an external storage device. The memory is used to store computer programs and other programs and data required by the terminal device. The memory may also be used to temporarily store data that has been output or is to be output.

[0084] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0085] In the several embodiments provided in the present application, it should be understood that the disclosed systems, devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be an indirect coupling or communication connection through some interfaces, devices or units, which can be electrical, mechanical or other forms.

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

[0087] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0088] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk and other media that can store program codes.

[0089] As described above, the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features thereof may be replaced by equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for calculating the angular unconformity erosion thickness of a sedimentary basin, characterized in that: The following steps are involved: Obtain the position of the angular unconformity of the sedimentary basin to be calculated for the denudation thickness and obtain N samples along the direction of the vertical strata of the sedimentary basin at the position of the angular unconformity; Performing a fluid inclusion uniform temperature test on N samples to obtain N fluid inclusion temperature data and a surface temperature; performing a fluid inclusion pressure calculation on the N fluid inclusion temperature data and a surface temperature to obtain N pressure data and a surface pressure; The uniform temperature of the N fluid inclusion temperature data and the N pressure data are respectively processed by layer-wise projection correction and linear fitting to obtain temperature gradient correction data and pressure gradient correction data; The uniform temperature of the fluid inclusion temperature data and the temperature gradient correction data are calculated by using an erosion thickness calculation formula, or the pressure data and the pressure gradient correction data are calculated to obtain the erosion thickness of each sample in the sedimentary basin; Wherein, the fluid inclusion temperature data includes the uniform temperature, initial melting temperature, freezing point temperature and composition of each sample; The calculation formula for the erosion thickness is: E=ΔT / G-Δy, where E is the erosion thickness of each sample in the sedimentary basin, ΔT is the difference between the homogenization temperature of the fluid inclusion temperature data of the sample and the surface temperature or the difference between the pressure data and the surface pressure, G is the temperature gradient correction data or the pressure gradient correction data, and Δy is the vertical distance between the sample and the angular unconformity.

2. The method for calculating the angular unconformity erosion thickness of a sedimentary basin according to claim 1, characterized in that: Before the fluid inclusion homogenization temperature test is performed on the N samples, the following conditions are included: the position of the angular unconformity is the overlying strata of the sedimentary basin, or the position of the angular unconformity is the underlying strata of the sedimentary basin and the former maximum burial depth of the underlying strata is greater than the current burial depth.

3. The method for calculating the angular unconformity erosion thickness of a sedimentary basin according to claim 2, characterized in that: include: The maximum burial depth of the underlying strata of the sedimentary basin was obtained by comparing vitrinite reflectance and regional stratigraphy.

4. The method for calculating the angular unconformity erosion thickness of a sedimentary basin according to claim 1, characterized in that: include: At the location of the angular unconformity, the sedimentary basin is sampled by surface or shallow well sampling along the direction perpendicular to the stratigraphic strike of the sedimentary basin to obtain N samples.

5. A system for calculating the angular unconformity erosion thickness of a sedimentary basin, characterized in that: include: Parameter acquisition module, parameter processing module, data processing module and calculation module; The parameter acquisition module is used to obtain the position of the angular unconformity of the sedimentary basin to be calculated for the denudation thickness and to obtain N samples along the direction of the vertical stratum strike of the sedimentary basin at the position of the angular unconformity; The parameter processing module is used to perform fluid inclusion uniform temperature testing on N samples to obtain N fluid inclusion temperature data and a surface temperature; and to perform fluid inclusion pressure calculation on the N fluid inclusion temperature data and one surface temperature to obtain N pressure data and a surface pressure; The data processing module is used to process the uniform temperature of the N fluid inclusion temperature data and the N pressure data by layer-by-layer projection correction and linear fitting to obtain temperature gradient correction data and pressure gradient correction data; The calculation module is used to calculate the homogenization temperature of the fluid inclusion temperature data and the temperature gradient correction data, or to calculate the pressure data and the pressure gradient correction data, by using the denudation thickness calculation formula, to obtain the denudation thickness of each sample in the sedimentary basin; Wherein, the fluid inclusion temperature data includes the uniform temperature, initial melting temperature, freezing point temperature and composition of each sample; The calculation formula for the erosion thickness is: E=ΔT / G-Δy, where E is the erosion thickness of each sample in the sedimentary basin, ΔT is the difference between the homogenization temperature of the fluid inclusion temperature data of the sample and the surface temperature or the difference between the pressure data and the surface pressure, G is the temperature gradient correction data or the pressure gradient correction data, and Δy is the vertical distance between the sample and the angular unconformity.

6. The angular unconformity erosion thickness calculation system of a sedimentary basin according to claim 5, characterized in that: The parameter acquisition module is used to perform fluid inclusion homogenization temperature tests on N samples according to the position of the angular unconformity being the overlying strata in the sedimentary basin, or the position of the angular unconformity being the underlying strata in the sedimentary basin and the former maximum burial depth of the underlying strata being greater than the current burial depth.

7. The system for calculating the angular unconformity erosion thickness of a sedimentary basin according to claim 5, characterized in that: The parameter acquisition module is used to obtain N samples by using surface or shallow well sampling in the sedimentary basin along the direction of the vertical stratum strike of the sedimentary basin at the position of the angular unconformity.

8. A device for calculating the angular unconformity erosion thickness of a sedimentary basin, characterized in that: including a processor and a memory; The memory is used to store program codes and transmit the program codes to the processor; The processor is used to execute the method for calculating the angular unconformity erosion thickness of a sedimentary basin according to any one of claims 1 to 4 according to the instructions in the program code.

Images