Method and system for estimating organic carbon burial rate and amount in deep fault lakes

Abstract

The present invention belongs to the technical field of sedimentology, and discloses a method and system for estimating the burial rate and burial amount of organic carbon in deep fault lakes. The method obtains the volume of dark mud shale deposited per unit time based on the obtained distribution range data of the demarcated dark mud shale per unit time; the density of the dark mud shale is calculated using the dark mud shale density calculation formula; the mass obtained by the mass deposition rate formula of the dark mud shale is divided by the time to obtain the average mass deposition rate of the dark mud shale per unit time; the burial rate of organic carbon is calculated; and the organic carbon burial amount of the lake is calculated. The deposition rate obtained by the present invention using the astronomical cycle method is more accurate; the deposition rate is accurate, and the subsequent calculation results are more accurate.

Description

Technical Field

[0001] The present invention belongs to the technical field of sedimentology, and in particular relates to a method and system for estimating the burial rate and amount of organic carbon in deep fault lakes. Background Art

[0002] The study of organic carbon source-sink system has a long history. Now with climate change, the study of the cycle process and dynamic changes among the three major carbon pools in the world needs to be further strengthened. Lakes are an important part of the terrestrial carbon pool. With the rapid development of society, the destruction of terrestrial green ecosystems such as forests by human activities and the large-scale use of fossil fuels such as oil and coal have increased the carbon storage in the atmosphere. With the increase of CO2 in the atmosphere, terrestrial plants absorb CO2 through photosynthesis and consume it as terrestrial organic carbon and transport it to lakes through rivers. In the past decade, the deposition rate and burial amount of organic carbon in lake basins have received increasing attention. In general, the organic carbon buried in lakes around the world is nearly half of that in the ocean, and small lakes account for 60%-70% of the organic carbon buried in lakes around the world. The study of organic carbon burial in lakes is of great significance to the study of global climate and environmental change.

[0003] The organic carbon in lake basins is usually composed of plant and microbial remains around the lakes. These organic matter are gradually deposited and buried in the lake basins. Over time, these organic matter are gradually converted into resources such as oil and natural gas. About 20% of the world's currently produced oil and natural gas comes from ancient large lake basins. The organic carbon buried in lake basins is also one of the important sources of oil and gas resources on the earth. Therefore, it is very necessary to strengthen the research on the deposition and burial of organic carbon in lakes.

[0004] Through the above analysis, the problems and defects of the existing technology are as follows: the existing technology has low recognition accuracy for the organic carbon deposition rate and organic carbon burial volume of deep ancient fault lakes, and the calculation results of the deposition rate are inaccurate. Summary of the invention

[0005] In order to overcome the problems existing in the related art, the disclosed embodiments of the present invention provide a method and system for estimating the burial rate and amount of organic carbon in deep fault lakes.

[0006] The technical scheme is as follows: A method for estimating the burial rate and amount of organic carbon in deep fault lakes, comprising:

[0007] S1. According to the deposition rate, the deposition thickness of the dark shale per unit time is obtained, and the volume of the dark shale is obtained by combining the distribution range of the dark shale per unit time.

[0008] S2. Calculate the density of dark shale using the calculation formula for the density of dark shale;

[0009] S3. Based on the density data of dark mud shale and using the mass deposition rate formula of dark mud shale, the average mass deposition rate of dark mud shale per unit time is obtained;

[0010] S4. Calculate the burial rate of organic carbon using the average mass deposition rate of dark shale per unit time;

[0011] S5. Calculate the amount of organic carbon burial in the lake using the obtained organic carbon burial rate.

[0012] In step S1, before obtaining the distribution range data of the dark shale within the defined unit time, the following steps need to be performed:

[0013] To determine the stratigraphic age of a deep lake, divide a million-year unit of time into x units of deposition time; x is the number of million years;

[0014] Define the distribution range S of dark shale per unit time.

[0015] Further, a time unit of one million years is divided into x units of deposition time, including:

[0016] (1) Analyze drilling natural gamma curve data;

[0017] (2) The natural gamma ray data of a depression in the Yellow River Delta were preprocessed using the Milankovitch cycle time series analysis method;

[0018] (3) Through MTM spectrum analysis, a long eccentricity orbit period of 405 years was obtained;

[0019] (4) Select detrended data and estimate the most likely sedimentation rate;

[0020] (5) Identify the astronomical orbital period, and then perform astronomical tuning after successful identification. If the tuning is successful, the astronomical age scale is finally obtained, and the sedimentation rate SR of the target depression in x unit time during the target period is obtained, in m / Ma.

[0021] Furthermore, the distribution range S of dark shale in unit time is delineated, including:

[0022] Based on the three-dimensional seismic data and drilling data, the combination of well and seismic data is adopted to obtain the thickness of the dark mud shale deposited in unit time according to the deposition rate. Combined with the analysis of the seismic profile of the drilling, the reflection axis corresponding to the top surface of the dark mud shale deposited in unit time is determined, and the distribution of the top surface on a seismic profile is obtained. The seismic reflection isochronous axis is tracked in the plane using the three-dimensional seismic network; the distribution range of the dark mud shale deposited in unit time on the plane is obtained; and the distribution range S of the dark mud shale in unit time is delineated in combination with the plane thickness distribution map.

[0023] Furthermore, the calculation formula for the volume of dark shale deposited per unit time is:

[0024] V SED =SR×S×1 / 3

[0025] Among them, V SED is the volume of dark mud shale, which is in the shape of a cone. The cone formula is used to calculate the volume of dark mud shale. SR is the sedimentation rate, and S is the distribution range of dark mud shale per unit time.

[0026] In step S2, the calculation formula of dark shale density is:

[0027] ρ SED = (2×%OC×1.1) + (%BSi×2.1) + [(100-2×%OC-%BSi)×2.65] / 100

[0028] Among them, ρ SED is the density of dark shale, %OC and %BSi are mass percentages, 1.1 is the density of organic matter, 2.1 is the density of biogenic silicon, and 2.65 is the density of minerals, i.e. silicate and carbonate minerals, all in g.cm -3 .

[0029] In step S3, the mass deposition rate of dark shale is calculated as follows:

[0030] MAR SED = (1-φ) × V SED ×ρ SED

[0031] Among them, MAR SED is the mass deposition rate of dark shale, φ is the sediment porosity; m=ρv, m is the mass of dark shale deposited per unit time, and the deposition rate of dark shale is expressed in t / Ma.

[0032] In step S4, the calculation formula for the burial rate of organic carbon is:

[0033] MARoc=MAR SED × (% / OC / 100)

[0034] MARoc is the burial rate of organic carbon.

[0035] Furthermore, the formula for calculating the organic carbon storage in lakes is: MAR per unit time OC × sedimentation time and then add them together to get the lake organic carbon burial capacity.

[0036] Another object of the present invention is to provide a system for estimating the burial rate and amount of organic carbon in deep fault lakes, which implements the method for estimating the burial rate and amount of organic carbon in deep fault lakes, and the system comprises:

[0037] The stratigraphic age determination module of deep-time lakes is used to determine the stratigraphic age of deep-time lakes, dividing a million-year unit of time into x units of sedimentation time; x is the number of million years;

[0038] The dark mud shale distribution range delineation module is used to delineate the distribution range S of the dark mud shale within a unit time;

[0039] The dark mud shale volume calculation module is used to obtain the deposition thickness of the dark mud shale per unit time according to the deposition rate, and to obtain the volume of the dark mud shale in combination with the distribution range of the dark mud shale per unit time;

[0040] A dark mud shale density calculation module is used to calculate the density of dark mud shale using a dark mud shale density calculation formula;

[0041] The average mass deposition rate acquisition module is used to obtain the average mass deposition rate of the dark shale per unit time based on the density data of the dark shale and the mass deposition rate formula of the dark shale;

[0042] The organic carbon burial rate calculation module is used to calculate the organic carbon burial rate using the average mass deposition rate of dark mud shale per unit time;

[0043] The lake organic carbon burial volume calculation module is used to calculate the lake organic carbon burial volume using the obtained organic carbon burial rate.

[0044] Combining all the above technical solutions, the beneficial effects of the present invention are as follows: the organic carbon buried in the world's lakes is nearly half of that in the ocean, and small lakes account for 60%-70% of the global lake organic carbon burial. The study of lake organic carbon burial is of great significance to the study of global climate and environmental changes. About 20% of the oil and natural gas currently being produced in the world comes from ancient large lake basins, and the organic carbon buried in lake basins is also one of the important sources of oil and gas resources on the earth. Therefore, it is very necessary to strengthen the research on the deposition and burial of lake organic carbon.

[0045] The existing technology mainly calculates the organic carbon burial of modern lakes, while the method of the present invention calculates deep-time fault lakes. The existing technology calculates the sedimentation rate by dividing the sedimentation thickness by the sedimentation time, and the average sedimentation rate is calculated. The sedimentation rate obtained by the present invention using the astronomical cycle method is more accurate; the sedimentation rate is accurate, and the subsequent calculation results are more accurate.

[0046] The existing methods for estimating the deposition rate and burial of organic carbon in lakes in China are for modern lakes, while this invention studies the deposition rate of organic carbon in deep lakes and estimates the amount of organic carbon burial in deep lakes. (3) In the past, the deposition rate in organic carbon burial was calculated by dividing the sediment thickness by the deposition time to obtain only the average deposition rate, which also made the subsequent calculation results of organic carbon burial inaccurate. However, the advantages of the Milankovitch cycle in calculating the deposition rate are linked to organic carbon burial, which provides reference for the subsequent calculation of organic carbon burial and deposition rate in deep lake basins. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] The accompanying drawings herein are incorporated in and constitute a part of the specification, illustrate embodiments consistent with the present disclosure, and together with the description, serve to explain the principles of the present disclosure;

[0048] Figure 1 is a flow chart of a method for estimating the burial rate and amount of organic carbon in deep fault lakes provided by an embodiment of the present invention;

[0049] Figure 2 It is a schematic diagram of the demarcation range of the seismic network provided by an embodiment of the present invention;

[0050] Figure 3 It is a schematic diagram of deposition thickness per unit time;

[0051] Figure 4 This is a lake seismic profile of a certain area in the ancient Yellow River Delta;

[0052] Figure 5 This is a planar distribution map of dark mud shale in a certain area of ​​the Yellow River Delta. DETAILED DESCRIPTION

[0053] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the specific embodiments of the present invention are described in detail below in conjunction with the accompanying drawings. In the following description, many specific details are set forth to facilitate a full understanding of the present invention. However, the present invention can be implemented in many other ways different from those described herein, and those skilled in the art can make similar improvements without violating the connotation of the present invention, so the present invention is not limited by the specific implementation disclosed below.

[0054] The innovations of the present invention are: ① combining the seismic network and the Milankovitch cycle to more accurately estimate the volume of dark mud shale; ② for the first time, a density formula and a mass deposition rate formula for calculating dark mud shale are proposed, providing a new idea for estimating the organic carbon burial amount and the organic carbon mass deposition rate.

[0055] Example 1: The method for estimating the burial rate and amount of organic carbon in deep fault lakes provided by the embodiment of the present invention includes:

[0056] S1. According to the deposition rate, the deposition thickness of the dark shale per unit time is obtained, and the volume of the dark shale is obtained by combining the distribution range of the dark shale per unit time.

[0057] First, determine the stratigraphic age of the deep lake, and divide it into x (x is the number of million years) units of deposition time, with one million years as the unit time;

[0058] (1) Analyze the natural gamma ray curve data of drilling (gamma ray curve data are usually used as a proxy indicator of paleoclimate);

[0059] (2) The natural gamma ray data of a depression in the Yellow River Delta were preprocessed (detrended) using the Milankovitch cycle time series analysis method;

[0060] (3) Through MTM spectrum analysis, a long eccentricity orbital period of 405 years was obtained (405 years is considered to be the most stable astronomical orbital period in geological history and the most commonly used orbital period in establishing astronomical age scales);

[0061] (4) Select detrended data and estimate the most likely sedimentation rate;

[0062] (5) Identify the astronomical orbital period, and then perform astronomical tuning after successful identification. If the tuning is successful, the (floating) astronomical age scale is finally obtained, and the sedimentation rate SR (m / Ma) in the target period of the depression within x unit time is obtained.

[0063] S2. Calculate the density of dark shale using the calculation formula for the density of dark shale;

[0064] S3. Based on the density data of dark mud shale and using the mass deposition rate formula of dark mud shale, the average mass deposition rate of dark mud shale per unit time is obtained;

[0065] S4. Calculate the burial rate of organic carbon using the average mass deposition rate of dark shale per unit time;

[0066] S5. Calculate the amount of organic carbon burial in the lake using the obtained organic carbon burial rate.

[0067] The distribution range S of dark shale in unit time is delineated, specifically including: based on 3D seismic data and drilling data, the thickness of dark shale deposited in unit time can be known according to the deposition rate, and combined with the analysis of the seismic profile of the well (calibrated by the well), the reflection axis corresponding to the top surface of the dark shale deposited in unit time (such as Figure 2Schematic diagram of the delineation range of the seismic network. Taking the blue line in the figure as an example, it is the top surface of the dark mudstone shale calibrated by drilling), and then the distribution of the top surface on a seismic profile is obtained, and then the three-dimensional seismic network is used to track the seismic reflection phase axis in the plane (taking the blue line in the figure as an example, it will be extended on the plane after being calibrated by drilling. Several seismic east-west profiles in the area are combined with the southeast-northwest profiles to form a network, which can delineate a range); that is, the distribution range of the dark mudstone deposited in a unit time on the plane can be obtained; and then combined with its plane thickness distribution map to delineate the distribution range S of the dark mudstone in a unit time.

[0068] According to the distribution range S of dark shale in the above-mentioned unit time, V SED =SR×S×1 / 3(V SED is the volume of dark mud shale. At this time, its shape is approximately a cone. The volume of dark mud shale is calculated using the cone formula), that is, the volume of dark mud shale deposited per unit time is obtained.

[0069] Regarding the density of dark shale, the present invention innovatively proposes:

[0070] The density of dark shale is calculated using ρ SED = (2 × %OC × 1.1) + (%BSi × 2.1) + [ (100-2 × %OC-%BSi) × 2.65] / 100 formula calculation (%OC, %BSi are mass percentages, 1.1 is the density of organic matter, 2.1 is the density of biosilicon, 2.65 is the density of minerals, i.e. silicate and carbonate minerals, the unit is g.cm -3 );

[0071] The mass deposition rate of dark shale is innovatively proposed in this invention: the mass deposition rate of dark shale is calculated by using MAR SED = (1-φ) × V SED ×ρ SED It is concluded that φ is the porosity of sediment; the principle is m=ρv, m is the mass of dark mud shale deposited per unit time, that is, the deposition rate of dark mud shale (t / Ma); the obtained mass minus the time can be obtained to get the average mass deposition rate of dark mud shale per unit time.

[0072] The present invention innovatively proposes to calculate the mass deposition rate of organic carbon in dark shale using MARoc=MAR SED × (% / OC / 100).

[0073] The organic carbon storage in the lake is the MAR per unit time. OC × deposition time and then add them together to get the result.

[0074] Through the above embodiments, it can be known that the data formula and data involved in the present invention are as shown in Table 1;

[0075] Table 1 Data formula and data

[0076]

[0077] Embodiment 2, the embodiment of the present invention provides a system for estimating the burial rate and amount of organic carbon in deep fault lakes, the system comprising:

[0078] The stratigraphic age determination module of deep-time lakes is used to determine the stratigraphic age of deep-time lakes, dividing a million-year unit of time into x units of sedimentation time; x is the number of million years;

[0079] The dark mud shale distribution range delineation module is used to delineate the distribution range S of the dark mud shale within a unit time;

[0080] The dark mud shale volume calculation module is used to obtain the deposition thickness of the dark mud shale per unit time according to the deposition rate, and to obtain the volume of the dark mud shale in combination with the distribution range of the dark mud shale per unit time;

[0081] A dark mud shale density calculation module is used to calculate the density of dark mud shale using a dark mud shale density calculation formula;

[0082] The average mass deposition rate acquisition module is used to obtain the average mass deposition rate of the dark shale per unit time based on the density data of the dark shale and the mass deposition rate formula of the dark shale;

[0083] The organic carbon burial rate calculation module is used to calculate the organic carbon burial rate using the average mass deposition rate of dark mud shale per unit time;

[0084] The lake organic carbon burial volume calculation module is used to calculate the lake organic carbon burial volume using the obtained organic carbon burial rate.

[0085] Application example. As a practical application implementation, the method for estimating the burial rate and amount of organic carbon in deep fault lakes provided in the embodiment of the present invention includes:

[0086] a, the duration of the second section of Kong strata is about 3 million years, and the unit time in the following text is 1 million years. The natural gamma data of a depression in a certain area of ​​the Yellow River Delta were preprocessed (detrended), and the energy spectrum was analyzed by MTM to obtain a long eccentricity orbital period of 405yr. The detrended data were selected to estimate the most likely sedimentation rate, and then an astronomical drive test was performed. If the astronomical drive test was successful, astronomical tuning was performed. After successful tuning, the (floating) astronomical age scale was finally obtained, and the sedimentation rates of the three unit times were 195m / Ma, 180m / Ma and 248m / Ma respectively.

[0087] b. The distribution range of dark mud shale is calculated by combining drilling and seismic profile data. The deposition rate per unit time obtained from the drilling data can be used to obtain the deposition thickness of dark mud shale per unit time (1Ma). Combined with the analysis of the seismic profile of the well (calibrated by the drilling), the reflection axis corresponding to the top surface of the dark mud shale deposition per unit time is determined, and then the distribution of the top surface on a seismic profile is obtained. The three-dimensional seismic network is then used to track the seismic reflection isochronous axis in the plane, and the distribution range of the dark mud shale deposited per unit time on the plane can be obtained. Similarly, the distribution range of dark mud shale in different unit times is obtained. It is estimated that the distribution range of dark mud shale in the ancient lake depression of a certain area in the Yellow River Delta for three unit times is approximately 930km. 2 、367km 2 、88.5km 2 ; and the sediment thickness in the three unit times is 195m, 180m, and 248m respectively. Figure 3 Schematic diagram of deposition thickness per unit time, Figure 4 Seismic profile of a lake in a certain area of ​​the ancient Yellow River Delta. Figure 5 Planar distribution of dark mud shale in a certain area of ​​the Yellow River Delta.

[0088] c. Substitute the area and thickness values ​​obtained in the above steps into V 暗色 =SR×S×1 / 3 (approximately the volume of the cone), the volume of dark shale is about 139.36 km 3 、44.59km 3 、7.24km 3 .

[0089] d, dark mudstone composition is taken as an example of organic matter, biogenic silica and minerals (silicate, carbonate minerals). The mass percentage of organic matter is 5%, biogenic silica is 30%, and minerals are 60%. Substituting the values ​​into the formula ρ SED = (2 × %OC × 1.1) + (%BSi × 2.1) + [ (100-2 × %OC-%BSi) × 2.65] / 100, it can be concluded that the density of dark shale is 2.42 g / cm3 .

[0090] e, the porosity is 20%, and the above value is substituted into MAR 暗色 = (1-φ) × V 暗色 ×ρ 暗色 The mass deposition rates per unit time are 2.69×10 11 t / Ma, 8.63×10 10 t / Ma, 1.4×10 10 t / Ma.

[0091] f, after obtaining the mass deposition rate of dark shale, substitute it into MARoc=MAR 暗色dan × (% / OC / 100), we can get the organic carbon mass deposition rates in the three unit times to be approximately 13.4×10 9 t / Ma, 4.31×10 9 t / Ma, 0.7×10 9 t / Ma.

[0092] g, Final OC T1 =MARoc T1 ×T1 to obtain the amount of organic carbon buried in the lake in the first unit time. Similarly, the amount of organic carbon buried in T2 and T3 can be obtained respectively. The sum of the three is the total amount of organic carbon buried in the lake during this period, which is 18.41×10 9 t.

[0093] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0094] Since the information interaction, execution process, etc. between the above-mentioned devices / units are based on the same concept as the method embodiment of the present invention, their specific functions and technical effects can be found in the method embodiment part and will not be repeated here.

[0095] Those skilled in the art can clearly understand that for the convenience and simplicity of description, only the division of the above-mentioned functional units and modules is used as an example for illustration. In practical applications, the above-mentioned function allocation can be completed by different functional units and modules as needed, that is, the internal structure of the device can be divided into different functional units or modules to complete all or part of the functions described above. The functional units and modules in the embodiments can be integrated in one processing unit, or each unit can exist physically separately, or two or more units can be integrated in one unit. The above-mentioned integrated unit can be implemented in the form of hardware or in the form of software functional units. In addition, the specific names of the functional units and modules are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of the present invention. The specific working process of the units and modules in the above-mentioned system can refer to the corresponding process in the aforementioned method embodiment.

[0096] To further illustrate the relevant effects of the embodiments of the present invention, the following experiments were conducted: Organic carbon burial volume is closely related to oil and gas exploration. Estimating the organic carbon burial volume of a deep lake basin is of great reference significance for estimating the oil and gas equivalent of the lake basin. Traditional methods for estimating oil and gas equivalent, such as the fault block method, are slow in efficiency, have a long cycle, and the data may be inaccurate. However, our method of estimating the organic carbon burial volume can provide a more accurate estimate of the oil and gas equivalent.

[0097] The above description is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited thereto. Any modifications, equivalent substitutions and improvements made by any technician familiar with the technical field within the technical scope disclosed by the present invention and within the spirit and principles of the present invention should be covered within the protection scope of the present invention.

Claims

1. A method for estimating the burial rate and amount of organic carbon in deep fault lakes, characterized in that: The method includes: S1. According to the deposition rate, the deposition thickness of the dark shale per unit time is obtained, and the volume of the dark shale is obtained by combining the distribution range of the dark shale per unit time. S2. Calculate the density of dark shale using the calculation formula for the density of dark shale; S3. Based on the density data of dark mud shale and using the mass deposition rate formula of dark mud shale, the average mass deposition rate of dark mud shale per unit time is obtained; S4. Calculate the burial rate of organic carbon using the average mass deposition rate of dark shale per unit time; S5. Calculate the amount of organic carbon burial in the lake using the obtained organic carbon burial rate.

2. The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to claim 1, characterized in that: In step S1, before obtaining the distribution range data of the dark shale within the defined unit time, the following steps need to be performed: To determine the stratigraphic age of a deep lake, divide a million-year unit of time into x units of deposition time; x is the number of million years; Define the distribution range S of dark shale per unit time.

3. The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to claim 2, characterized in that: One million years is divided into x units of deposition time, including: (1) Analyze drilling natural gamma curve data; (2) The natural gamma ray data of a depression in the Yellow River Delta were preprocessed using the Milankovitch cycle time series analysis method; (3) Through MTM spectrum analysis, a long eccentricity orbit period of 405 years was obtained; (4) Select detrended data and estimate the most likely sedimentation rate; (5) Identify the astronomical orbital period, and then perform astronomical tuning after successful identification. If the tuning is successful, the astronomical age scale is finally obtained, and the sedimentation rate SR of the target depression in x unit time during the target period is obtained, in m / Ma.

4. The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to claim 2, characterized in that: The distribution range S of dark shale in unit time is defined, including: Based on the three-dimensional seismic data and drilling data, the combination of well and seismic data is adopted to obtain the thickness of the dark mud shale deposited in unit time according to the deposition rate. Combined with the analysis of the seismic profile of the drilling, the reflection axis corresponding to the top surface of the dark mud shale deposited in unit time is determined, and the distribution of the top surface on a seismic profile is obtained. The seismic reflection isochronous axis is tracked in the plane using the three-dimensional seismic network; the distribution range of the dark mud shale deposited in unit time on the plane is obtained; and the distribution range S of the dark mud shale in unit time is delineated in combination with the plane thickness distribution map.

5. The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to claim 4, characterized in that: The calculation formula for the volume of dark shale deposited per unit time is: In SED =SR×S×1 / 3 Among them, V SED is the volume of dark mud shale, which is in the shape of a cone. The cone formula is used to calculate the volume of dark mud shale. SR is the sedimentation rate, and S is the distribution range of dark mud shale per unit time.

6. The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to claim 1, characterized in that: In step S2, the calculation formula of dark shale density is: ρ SED =(2×%OC×1.1)+(%BSi×2.1)+[(100-2×%OC-%BSi)×2.65] / 100 Among them, ρ SED is the density of dark shale, %OC and %BSi are mass percentages, 1.1 is the density of organic matter, 2.1 is the density of biogenic silicon, and 2.65 is the density of minerals, i.e. silicate and carbonate minerals, all in g.cm -3 .

7. The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to claim 1, characterized in that: In step S3, the mass deposition rate of dark shale is calculated as follows: MAR SED = (1-φ) × V SED ×r SED Among them, MAR SED is the mass deposition rate of dark shale, φ is the sediment porosity; m=ρv, m is the mass of dark shale deposited per unit time, and the deposition rate of dark shale is expressed in t / Ma.

8. The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to claim 7, characterized in that: In step S4, the calculation formula for the burial rate of organic carbon is: MOROCCO=MAR SED ×(% / OC / 100) MARoc is the burial rate of organic carbon.

9. The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to claim 1, characterized in that: The formula for calculating the organic carbon storage in lakes is: MAR per unit time OC × sedimentation time and then add them together to get the lake organic carbon burial capacity.

10. A system for estimating the burial rate and amount of organic carbon in deep fault lakes, characterized in that: The method for estimating the burial rate and amount of organic carbon in deep fault lakes according to any one of claims 1 to 9 is implemented, the system comprising: The stratigraphic age determination module of deep-time lakes is used to determine the stratigraphic age of deep-time lakes, dividing a million-year unit of time into x units of sedimentation time; x is the number of million years; The dark mud shale distribution range delineation module is used to delineate the distribution range S of the dark mud shale within a unit time; The dark mud shale volume calculation module is used to obtain the deposition thickness of the dark mud shale per unit time according to the deposition rate, and to obtain the volume of the dark mud shale in combination with the distribution range of the dark mud shale per unit time; A dark mud shale density calculation module is used to calculate the density of dark mud shale using a dark mud shale density calculation formula; The average mass deposition rate acquisition module is used to obtain the average mass deposition rate of the dark shale per unit time based on the density data of the dark shale and using the mass deposition rate formula of the dark shale; The organic carbon burial rate calculation module is used to calculate the organic carbon burial rate using the average mass deposition rate of dark mud shale per unit time; The lake organic carbon burial volume calculation module is used to calculate the lake organic carbon burial volume using the obtained organic carbon burial rate.

Images