A method and system for processing data of a material source area and a deposition area

Abstract

The application provides a kind of source area and sedimentary area data processing method and system, the processing method includes: the bottom surface is explained to the intersection of top surface and basement along the base surface to ensure the intersection of both;Respectively, the top surface and bottom surface of the overlying strata in the sedimentary area are converted into time-depth;The top surface and bottom surface of the overlying strata in the sedimentary area after time-depth conversion are respectively interpolated and encrypted;The negative value of the stratum thickness of the sedimentary area is calculated according to the interpolation and encryption data;Horizontal seismic interpretation is made, and the basement of the source area is interpreted;The basement of the source area and the horizontal interpretation layer are respectively converted into time-depth;The basement after time-depth conversion of the source area and the horizontal interpretation layer after time-depth conversion are respectively interpolated and encrypted;The stratum thickness of the source area above the horizontal plane during the stratum deposition period is calculated according to the interpolation and encryption data of the source area;Elevation data is calculated. It can make the integrated paleogeomorphology map more coordinated and accurate than the late splicing of the paleogeomorphology map of a single sedimentary area or source area.

Description

Technical Field

[0001] The present invention relates to the technical field of petroleum exploration, and in particular to a method and system for processing data of provenance areas and sedimentary areas. Background Art

[0002] In oil and gas exploration, except for the special study of water depth during the deposition period, paleo-geomorphological maps are generally drawn without considering water depth, converting the thickness of sedimentary strata into paleo-geomorphological maps during the deposition period, or directly using the current geomorphology of the basement strata in the provenance area as the paleo-geomorphology during the deposition period. Figure 1 It is generally a single expression of the sedimentary area or provenance area, which can only reflect the paleogeomorphic environment of the sedimentary system under the water surface during the sedimentary period or the paleogeomorphic characteristics of the provenance area that provided terrigenous debris for the sedimentary system. When analyzing geological information such as the transportation channels and sedimentary areas of terrigenous debris from the provenance area to the sedimentary area, geological researchers are usually required to manually combine the paleogeomorphic maps of the sedimentary area and the provenance area to conduct comprehensive research. On the one hand, this working method is not easy to accurately splice the boundaries of the sedimentary area and the provenance area. On the other hand, the sources of the paleogeomorphic maps of different regions are often different, and there are also difficulties in adjusting and unifying the vertical scales of the two. Summary of the Invention

[0003] In view of the above problems, the present invention is proposed to provide a method and system for processing provenance and deposition area data that overcomes the above problems or at least partially solves the above problems.

[0004] According to one aspect of the present invention, a method for processing data of a provenance area and a deposition area is provided, the processing method comprising:

[0005] Step 1: On a seismic section with overlapping sedimentary features, at the intersection of the bottom surface and the basement, interpret the bottom surface along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect;

[0006] Step 2, performing time-depth conversion on the top and bottom surfaces of the overlapping strata in the sedimentary area;

[0007] Step 3: interpolate and encrypt the top and bottom surfaces of the overlapping strata in the sedimentary area after time-depth conversion;

[0008] Step 4, calculating the negative value of the stratum thickness in the sedimentary area based on the interpolated encrypted data;

[0009] Step 5: Perform horizontal seismic interpretation using a sparse grid to interpret the basement of the provenance area;

[0010] Step 6: Perform time-depth conversion on the basement and horizontal interpretation layers of the provenance area respectively;

[0011] Step 7: interpolate and encrypt the basement and horizontal interpretation layers of the provenance area after time-depth conversion;

[0012] Step 8, calculating the thickness of the provenance strata above the horizontal plane during the deposition period based on the interpolated infill data of the provenance;

[0013] Step 9, calculating the normalized provenance thickness;

[0014] Step 10: Calculate continuous elevation data from the provenance area above the water surface to the sedimentation area below the water surface based on the negative value of the stratum thickness of the sedimentation area obtained in step 4 and the stratum thickness of the provenance area obtained in step 9.

[0015] Optionally, step 1, on a seismic profile having overlapping sedimentary features, at the intersection of the bottom surface and the basement, interpreting the bottom surface along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect, specifically includes:

[0016] In areas with overlapping sedimentary features, the top and bottom surfaces of the strata are interpreted using seismic workstations. At the intersection of the bottom surface and the basement, the bottom surface is interpreted along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect and avoid null values ​​during interpolation.

[0017] Manually connect the strata disconnected by the fault along the fault plane at the fault to ensure stratum continuity.

[0018] Optionally, the step 2 of performing time-depth conversion on the top and bottom surfaces of the overlying strata in the sedimentary area specifically includes:

[0019] At the seismic workstation, the interpreted top and bottom surfaces are converted into time and depth using the velocity field formula of the study area to obtain two depth layers.

[0020] Optionally, step 4, calculating the negative value of the formation thickness of the sedimentary area based on the interpolated encrypted data, specifically includes: subtracting the interpolated data of the formation bottom surface after time-depth conversion from the interpolated data of the formation top surface after time-depth conversion to obtain the negative value of the formation thickness of the sedimentary area.

[0021] Optionally, the step 5, performing horizontal seismic interpretation in a sparse grid manner and interpreting the basement of the provenance area specifically includes: performing horizontal seismic interpretation in a sparse grid manner from the intersection of the top surface of the sedimentary area stratum and the basement to the provenance area and interpreting the basement of the provenance area according to research needs or data range limitations.

[0022] Optionally, the step 6 of performing time-depth conversion on the basement and horizontal interpretation layer of the provenance area respectively includes: performing time-depth conversion on the basement and horizontal interpretation layer of the provenance area respectively using the velocity field formula of the study area.

[0023] Optionally, the step 8, calculating the stratigraphic thickness of the provenance area above the horizontal plane during the stratigraphic deposition period based on the interpolated encrypted data of the provenance area, specifically includes: subtracting the interpolated data of the basement after time-depth conversion from the interpolated data of the horizontal interpretation layer after time-depth conversion of the provenance area, to obtain the stratigraphic thickness of the provenance area above the horizontal plane during the stratigraphic deposition period under geological understanding.

[0024] Optionally, step 9, calculating the normalized provenance area stratum thickness specifically includes: dividing the maximum stratum thickness of the sedimentary area by the maximum stratum thickness of the provenance area to obtain a normalization coefficient, and then multiplying the provenance area stratum thickness data by the normalization coefficient to obtain the normalized provenance area stratum thickness.

[0025] Optionally, the step 10, based on the negative value of the stratum thickness of the sedimentary area obtained in step 4 and the stratum thickness of the provenance area obtained in step 9, calculates the continuous elevation data from the provenance area above the water surface to the sedimentary area below the water surface, specifically including: integrating the negative value of the stratum thickness of the sedimentary area obtained in step 4 and the normalized stratum thickness of the provenance area obtained in step 9 to obtain continuous elevation data from the provenance area above the water surface to the sedimentary area below the water surface.

[0026] The present invention further provides a provenance and sedimentation area data processing system, which applies the above-mentioned provenance and sedimentation area data processing method, and the processing system includes:

[0027] An interpretation module is used to interpret the bottom surface along the basement surface to the intersection of the top surface and the basement on a seismic section with overlapping sedimentary features to ensure that the two intersect;

[0028] Sedimentary area time-depth conversion module, used to perform time-depth conversion on the top and bottom surfaces of the overlapping strata in the sedimentary area;

[0029] Sedimentary area interpolation and encryption module, used for interpolation and encryption of the top and bottom surfaces of the overlapping strata in the sedimentary area after time-depth conversion;

[0030] A sedimentary area stratum thickness calculation module is used to calculate the negative value of the stratum thickness in the sedimentary area based on the interpolated encrypted data;

[0031] The provenance basement interpretation module is used to interpret horizontal seismic data in a sparse grid and to interpret the basement of the provenance area;

[0032] The provenance time-depth conversion module is used to perform time-depth conversion on the basement and horizontal interpretation layers of the provenance area respectively;

[0033] The interpolation and encryption module is used to interpolate and encrypt the basement and horizontal interpretation layers of the provenance area after time-depth conversion;

[0034] The provenance area stratigraphic thickness module is used to calculate the provenance area stratigraphic thickness above the horizontal plane during the stratigraphic deposition period based on the interpolated and intensified data of the provenance area;

[0035] Normalization processing module, used to calculate the normalized provenance thickness;

[0036] The elevation data calculation module is used to calculate continuous elevation data from the provenance area above the water surface to the sedimentation area below the water surface based on the negative value of the stratum thickness of the sedimentation area obtained in step 4 and the stratum thickness of the provenance area obtained in step 9.

[0037] The present invention provides a method and system for processing data of a provenance area and a sedimentary area, the processing method comprising: step 1, on a seismic profile having overlapping sedimentary characteristics, at the intersection of the bottom surface and the basement, interpreting the bottom surface along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect; step 2, performing time-depth conversion on the top surface and the bottom surface of the overlapping stratum in the sedimentary area respectively; step 3, interpolating and encrypting the top surface and the bottom surface of the overlapping stratum in the sedimentary area after time-depth conversion; step 4, calculating the negative value of the stratum thickness of the sedimentary area based on the interpolated encrypted data; step 5, performing horizontal seismic interpretation and solving in a sparse grid manner Release the provenance basement; Step 6: Perform time-depth conversion on the provenance basement and horizontal interpretation layers; Step 7: Interpolate and indent the provenance basement and horizontal interpretation layers after time-depth conversion; Step 8: Calculate the provenance stratum thickness above the horizontal plane during the depositional period based on the interpolated and indented provenance data; Step 9: Calculate the normalized provenance stratum thickness; Step 10: Calculate continuous elevation data from the provenance above the water surface to the sedimentary area below the water surface based on the negative value of the sedimentary area stratum thickness obtained in Step 4 and the provenance area stratum thickness obtained in Step 9. This makes the integrated paleogeomorphological map more coordinated and accurate than later splicing paleogeomorphological maps based on single sedimentary or provenance areas.

[0038] The above description is only an overview of the technical solution of the present invention. In order to more clearly understand the technical means of the present invention, it can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are specifically listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for use in the description of the embodiments. 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 any creative work.

[0040] Figure 1A flow chart of a specific embodiment of the data processing method of the provenance area and the sedimentary area provided by the embodiment of the present application;

[0041] Figure 2 A schematic diagram of the top and bottom of the stratum interpreted by the specific embodiment of the present application;

[0042] Figure 3 A schematic diagram of the stratum bottom surface interpreted to the top surface boundary in the specific embodiment of the present application;

[0043] Figure 4 A schematic diagram of the boundary points of the same stratum on the upper and lower plates of the fault connected in the specific embodiment of the present application;

[0044] Figure 5 A schematic diagram of the horizontal seismic interpretation to the provenance area in the specific embodiment of the present application;

[0045] Figure 6 A schematic diagram of the base of the provenance area interpreted in the specific embodiment of the present application. DETAILED DESCRIPTION

[0046] Exemplary embodiments of the present disclosure will be described more fully hereinafter with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it is to be understood that the present disclosure can be embodied in various forms without being limited by the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the present disclosure to those skilled in the art.

[0047] The terms "include" and "have" and any variations thereof in the specification and the claims of the present application and the accompanying drawings are intended to cover the non-exclusive inclusion, for example, the inclusion of a series of steps or units.

[0048] The technical solutions of the present application will be described in further detail below in combination with the drawings and embodiments.

[0049] As shown in Figure 1 , Figure 1 A flow chart of the data processing method of the provenance area and the sedimentary area based on the stratum thickness for drawing the paleogeomorphology map of the present application.

[0050] Step 101, the top and bottom surfaces of the stratum are interpreted in the seismic workstation, as shown in Figure 2 . At the intersection of the bottom surface and the base, the bottom surface is interpreted along the base to the intersection of the top surface and the base to ensure the intersection, as shown in Figure 3 . The boundary points of the same stratum on the upper and lower plates of the fault are manually connected to ensure the continuity of the stratum. As shown in Figure 4 .

[0051] Step 102, the top surface and bottom surface of the deposition area are respectively converted into time-depth in the seismic workstation by using the velocity field formula of the study area.

[0052] Step 103, the top surface and bottom surface of the overlying strata in the deposition area are respectively interpolated and encrypted after time-depth conversion.

[0053] Step 104, the interpolated top surface is subtracted from the interpolated bottom surface to obtain the negative value (0 to negative value) of the stratum thickness of the deposition area.

[0054] Step 105, according to the research needs or the range limitation of the data, the horizontal seismic interpretation is made in the source area from the intersection of the stratum top surface and the base surface of the deposition area to the source area in a sparse grid manner, as shown in Figure 5 The base interpretation is made in the source area from the intersection of the top surface and the base, as shown in Figure 6

[0055] Step 106, the base and the horizontal interpretation layer of the source area are respectively converted into time-depth in the seismic workstation by using the velocity field formula of the study area.

[0056] Step 107, the interpolated base and the interpolated horizontal interpretation layer of the source area are respectively interpolated and encrypted.

[0057] Step 108, the interpolated horizontal interpretation layer of the source area is subtracted from the interpolated base to obtain the stratum thickness (positive value to 0) of the source area above the horizontal surface of the deposition period under the geological understanding.

[0058] Step 109, the maximum value of the stratum thickness of the deposition area is divided by the maximum value of the stratum thickness of the source area to obtain a normalization coefficient, and then the stratum thickness data of the source area is multiplied by the normalization coefficient to obtain the normalized stratum thickness (positive value to 0) of the source area.

[0059] Step 110, the negative value (0 to negative value) of the stratum thickness of the deposition area obtained in step 104 and the normalized stratum thickness (positive value to 0) of the source area obtained in step 9 are integrated into one, to obtain the continuous elevation data from the source area above the water surface to the deposition area below the water surface.

[0060] The purpose of the present application is to provide a processing method for integrating the data of the source area and the deposition area in the process of compiling the paleogeomorphology map based on the stratum with overlying deposition characteristics by using the stratum thickness.

[0061] ​The present invention provides a data processing method for source areas and sedimentary areas for compiling paleo-geomorphological maps based on stratum thickness. The data processing method for source areas and sedimentary areas for compiling paleo-geomorphological maps based on stratum thickness includes: step 1, on a seismic profile with overlapping sedimentary characteristics, interpreting the top and bottom surfaces of the stratum in a sparse grid manner. At the intersection of the bottom surface and the basement, interpret the bottom surface along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect, and manually connect the boundary points of the same stratum on the upper and lower plates of the fault at the fault to ensure stratum continuity; step 2, respectively perform time-depth conversion on the top and bottom surfaces of the overlapping strata in the sedimentary area using the velocity field formula of the study area; step 3, interpolate and encrypt the top and bottom surfaces of the overlapping strata in the sedimentary area after time-depth conversion; step 4, subtract the interpolation data of the bottom surface of the stratum after time-depth conversion from the interpolation data of the top surface of the stratum after time-depth conversion to obtain a negative value (0 to negative value) of the stratum thickness in the sedimentary area; step 5, according to the research If research needs or data scope are limited, a sparse grid horizontal seismic interpretation is performed from the intersection of the top of the sedimentary strata and the overlap boundary toward the provenance area, and the basement of the provenance area is interpreted. Step 6: Time-depth conversion is performed on the basement and horizontal interpretation layers of the provenance area using the velocity field formula for the study area. Step 7: The time-depth-converted basement and the time-depth-converted horizontal interpretation layers of the provenance area are interpolated and densified. Step 8: The interpolated data of the time-depth-converted basement are subtracted from the interpolated data of the time-depth-converted horizontal interpretation layers of the provenance area to obtain the provenance area stratum thickness above the horizontal plane during the sedimentary period (positive values ​​to 0) based on geological understanding. Step 9: The maximum sedimentary stratum thickness is divided by the maximum provenance stratum thickness to obtain a normalization coefficient. The provenance area stratum thickness data is then multiplied by the normalization coefficient to obtain the normalized provenance area stratum thickness (positive values ​​to 0). In step 10, the negative value (0 to negative value) of the stratum thickness of the sedimentary area obtained in step 4 is integrated with the normalized stratum thickness (positive value to 0) of the provenance area obtained in step 9 to obtain continuous elevation data from the provenance area above the water surface to the sedimentary area below the water surface.

[0062] The purpose of the present invention can also be achieved by the following technical measures:

[0063] In step 1, on a seismic section with overlapping sedimentary features, the top and bottom surfaces of the strata are interpreted using a sparse grid. At the intersection of the bottom surface and the basement, the bottom surface is interpreted along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect. At the fault, the boundary points of the same stratum above and below the fault are manually connected to ensure stratigraphic continuity.

[0064] In step 2, the top and bottom surfaces interpreted in step 1 are converted into time-depth using the velocity field formula of the study area.

[0065] In step 3, the top surface and bottom surface of the deposition area that have undergone time-depth conversion in step 2 are interpolated and encrypted respectively.

[0066] In step 4, the interpolated data of the formation bottom after time-depth conversion in step 3 is subtracted from the interpolated data of the formation top after time-depth conversion to obtain a negative value (0 to negative value) of the formation thickness in the sedimentary area.

[0067] In step 5, based on research needs or data range limitations, horizontal seismic interpretation is performed from the intersection of the top surface and basement of the sedimentary area to the provenance area in a sparse grid manner, and the basement of the provenance area is interpreted;

[0068] In step 6, the basement and horizontal interpretation layers of the provenance interpretation in step 5 are converted in time and depth.

[0069] In step 7, the basement after time-depth conversion of the provenance area in step 6 and the horizontal interpretation layer after time-depth conversion are interpolated and encrypted respectively.

[0070] In step 8, the interpolated data of the horizontal interpretation layer after time-depth conversion in the provenance area in step 7 is subtracted from the interpolated data of the basement after time-depth conversion to obtain the stratigraphic thickness of the provenance area above the horizontal plane of the stratigraphic deposition period under geological understanding (positive value to 0).

[0071] In step 9, the maximum stratum thickness in the sedimentary area is divided by the maximum stratum thickness in the provenance area to obtain a normalization coefficient, and then the stratum thickness data in the provenance area is multiplied by the normalization coefficient to obtain the normalized stratum thickness in the provenance area (positive value to 0).

[0072] In step 10, the negative value (0 to negative value) of the stratigraphic thickness of the sedimentary area obtained in step 4 is integrated with the normalized stratigraphic thickness (positive value to 0) of the provenance area obtained in step 9 to obtain continuous elevation data from the provenance area above the water surface to the sedimentary area below the water surface.

[0073] Beneficial effects: The data processing method of provenance and sedimentary areas for compiling paleogeomorphological maps based on stratum thickness in the present invention can make the integrated paleogeomorphological map more coordinated and accurate than the later splicing of paleogeomorphological maps using single sedimentary areas or provenance areas.

[0074] The above specific implementation methods further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above are only specific implementation methods of the present invention and are not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for processing provenance and sedimentation area data, characterized in that: The processing method comprises: Step 1: On a seismic section with overlapping sedimentary features, at the intersection of the bottom surface and the basement, interpret the bottom surface along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect; Step 2, performing time-depth conversion on the top and bottom surfaces of the overlapping strata in the sedimentary area; Step 3: interpolate and encrypt the top and bottom surfaces of the overlapping strata in the sedimentary area after time-depth conversion; Step 4, calculating the negative value of the stratum thickness in the sedimentary area based on the interpolated encrypted data; Step 5: Perform horizontal seismic interpretation using a sparse grid to interpret the basement of the provenance area; Step 6: Perform time-depth conversion on the basement and horizontal interpretation layers of the provenance area respectively; Step 7: interpolate and encrypt the basement and horizontal interpretation layers of the provenance area after time-depth conversion; Step 8, calculating the thickness of the provenance strata above the horizontal plane during the deposition period based on the interpolated infill data of the provenance; Step 9, calculating the normalized provenance thickness; Step 10: Calculate continuous elevation data from the provenance area above the water surface to the sedimentation area below the water surface based on the negative value of the stratum thickness of the sedimentation area obtained in step 4 and the stratum thickness of the provenance area obtained in step 9.

2. A method for processing provenance and sedimentation area data according to claim 1, characterized in that: The step 1, on a seismic section with overlapping sedimentary features, at the intersection of the bottom surface and the basement, interpreting the bottom surface along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect, specifically includes: In areas with overlapping sedimentary features, the top and bottom surfaces of the strata are interpreted using seismic workstations. At the intersection of the bottom surface and the basement, the bottom surface is interpreted along the basement surface to the intersection of the top surface and the basement to ensure that the two intersect and avoid null values ​​during interpolation. Manually connect the strata disconnected by the fault along the fault plane at the fault to ensure stratum continuity.

3. The method for processing provenance and sedimentation area data according to claim 1, characterized in that: The step 2 of performing time-depth conversion on the top and bottom surfaces of the overlapping strata in the sedimentary area specifically includes: At the seismic workstation, the interpreted top and bottom surfaces are converted into time and depth using the velocity field formula of the study area to obtain two depth layers.

4. The method for processing provenance and sedimentation area data according to claim 1, characterized in that: The step 4, calculating the negative value of the formation thickness of the sedimentary area based on the interpolated encrypted data, specifically includes: subtracting the interpolated data of the formation bottom surface after time-depth conversion from the interpolated data of the formation top surface after time-depth conversion to obtain the negative value of the formation thickness of the sedimentary area.

5. The method for processing provenance and sedimentation area data according to claim 1, characterized in that: The step 5, performing horizontal seismic interpretation in a sparse grid manner and interpreting the basement of the provenance area, specifically includes: performing horizontal seismic interpretation in a sparse grid manner from the intersection of the top surface of the sedimentary area and the basement to the provenance area and interpreting the basement of the provenance area according to research needs or data range limitations.

6. The method for processing provenance and sedimentation area data according to claim 1, characterized in that: Said step 6, performing time-depth conversion on the basement and horizontal interpretation layer of the provenance area respectively, specifically comprises: performing time-depth conversion on the basement and horizontal interpretation layer of the provenance area respectively using the velocity field formula of the study area.

7. The method for processing provenance and sedimentation area data according to claim 1, characterized in that: Said step 8, calculating the stratigraphic thickness of the provenance area above the horizontal plane during the stratigraphic deposition period based on the interpolated encrypted data of the provenance area, specifically includes: subtracting the interpolated data of the basement after time-depth conversion from the interpolated data of the horizontal interpretation layer after time-depth conversion of the provenance area, to obtain the stratigraphic thickness of the provenance area above the horizontal plane during the stratigraphic deposition period under geological understanding.

8. The method for processing provenance and sedimentation area data according to claim 1, characterized in that: The step 9, calculating the normalized provenance area stratum thickness, specifically includes: dividing the maximum stratum thickness of the sedimentary area by the maximum stratum thickness of the provenance area to obtain a normalization coefficient, and then multiplying the provenance area stratum thickness data by the normalization coefficient to obtain the normalized provenance area stratum thickness.

9. The method for processing provenance and sedimentation area data according to claim 1, characterized in that: The step 10, based on the negative value of the stratum thickness of the sedimentary area obtained in step 4 and the stratum thickness of the provenance area obtained in step 9, calculates the continuous elevation data from the provenance area above the water surface to the sedimentary area below the water surface, specifically including: integrating the negative value of the stratum thickness of the sedimentary area obtained in step 4 and the normalized stratum thickness of the provenance area obtained in step 9 to obtain continuous elevation data from the provenance area above the water surface to the sedimentary area below the water surface.

10. A provenance and sedimentation area data processing system, applying a provenance and sedimentation area data processing method according to any one of claims 1 to 9, characterized in that: The processing system comprises: An interpretation module is used to interpret the bottom surface along the basement surface to the intersection of the top surface and the basement on a seismic section with overlapping sedimentary features to ensure that the two intersect; Sedimentary area time-depth conversion module, used to perform time-depth conversion on the top and bottom surfaces of the overlapping strata in the sedimentary area; Sedimentary area interpolation and encryption module, used for interpolation and encryption of the top and bottom surfaces of the overlapping strata in the sedimentary area after time-depth conversion; A sedimentary area stratum thickness calculation module is used to calculate the negative value of the stratum thickness in the sedimentary area based on the interpolated encrypted data; The provenance basement interpretation module is used to interpret horizontal seismic data in a sparse grid and to interpret the basement of the provenance area; The provenance time-depth conversion module is used to perform time-depth conversion on the basement and horizontal interpretation layers of the provenance area respectively; The interpolation and encryption module is used to interpolate and encrypt the basement and horizontal interpretation layers of the provenance area after time-depth conversion; The provenance area stratigraphic thickness module is used to calculate the provenance area stratigraphic thickness above the horizontal plane during the stratigraphic deposition period based on the interpolated and intensified data of the provenance area; Normalization processing module, used to calculate the normalized provenance thickness; The elevation data calculation module is used to calculate continuous elevation data from the provenance area above the water surface to the sedimentation area below the water surface based on the negative value of the stratum thickness of the sedimentation area obtained in step 4 and the stratum thickness of the provenance area obtained in step 9.

Images