Well-free delta sedimentary area paleo-water depth contour map compiling method

Abstract

The invention provides a well-free delta sedimentary area paleo-water depth contour map compilation method, and belongs to the technical field of oil and gas geological exploration. The method comprises the following steps: firstly, determining an object source direction and a delta lobes range through geological background research, selecting an input point as an original point, arranging a plurality of groups of seismic sections at equal angle intervals, and determining a main object source section; tracking the top and bottom boundary lines of the delta lobes on each section, taking a point with the maximum thickness as a reference point, and forming an initial water depth contour line by connecting the reference point of the main object source section with reference points of other sections; further selecting multiple points along the top boundary of the main section, calculating the paleo-water depth value of each point in combination with the thickness data, and constructing a complete contour map in a mode of parallel initial contour lines. According to the method, the limitation of dependence on drilling data traditionally is broken through, and rapid reconstruction of the paleo-water depth in the well-free area is achieved by establishing the quantitative relation between the sedimentary thickness and the paleo-water depth through the geometrical morphology characteristics of the seismic section and the sedimentary body.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas geological exploration, and particularly relates to a method for compiling an isobath map of paleo-water depth in a delta sedimentary area without wells. Background Art

[0002] In oil and gas exploration, paleo-water depth analysis is one of the important bases for predicting the distribution of source rocks and favorable reservoir sand bodies, and is of great significance for paleo-environment reconstruction, basin analysis, sequence stratigraphy research, paleogeomorphology restoration, and the evaluation of source, reservoir, and cap rock conditions.

[0003] Common methods for paleo-water depth analysis include paleontology method, sedimentology method, geophysical method, geochemical method, and digital simulation method, etc. The paleontology method requires drilling core or cuttings data. When the data is not detailed enough, the results vary greatly and are easily affected by sediment reworking. Methods such as sedimentology method, geophysical method, and geochemical method need to rely on drilling data to obtain cuttings data and logging curves. The digital simulation method is applicable to the basin scale and requires a large amount of data accumulation. The shoreline trajectory method believes that the delta plain area belongs to the compensation area, and the calculated paleo-water depth is equal to 0. However, there is a certain slope from the delta plain to the delta front, and the paleo-water depth should be greater than 0. At the same time, the size of the accommodation space is affected by the shoreline position, and the accommodation space is often calculated too large. The pre-stack thickness calculated by the seismic data method is the cumulative thickness, resulting in a greatly overestimated calculated paleo-water depth.

[0004] In recent years, various methods have been developed for paleo-water depth analysis. For example, in the patent with the application publication number: CN114624775A, titled: A Comprehensive Quantitative Restoration Method for Paleo-Water Depth in a Sedimentary Lake Basin, it is proposed to comprehensively determine the paleo-water depth by using the shoreline trajectory method and the trace element (Fe / Co) method in mudstone, but a large amount of Fe and Co element analysis and test data obtained from drilling are required. In the patent with the application publication number: CN113970796A, titled: A Method for Accurately Restoring the Paleo-Water Depth of a Sedimentary Basin, it is proposed to flatten the top surface of the delta foreset structure, restore the compaction amount and erosion amount, and then use logging data and sample analysis and test data to restore the paleo-water depth of the basin. Similarly, drilling data and analysis and test data are required. In the patent with the publication number: CN111475920B, titled: A Method, System, Electronic Device, and Storage Medium for Obtaining the Paleo-Water Depth of a Deep-Water Basin, a method for restoring the paleo-water depth of a deep-water basin is proposed, which requires the use of paleontological information in cores and cuttings obtained from drilling, the characteristics of the shelf break zone, etc. The three-level sequence sedimentary configuration analysis method (Shu Liangfeng et al., 2022), these two methods are applicable to the large-scale paleo-water depth restoration at the basin level.

[0005] The above methods either rely on drilling data or are only applicable to the large-scale conditions of the basin, and are not applicable to the small-scale delta sedimentary area without drilling. Summary of the Invention

[0006] The present invention aims to solve at least one of the technical problems existing in the prior art. To this end, the present invention provides a method for compiling an isobath map of paleo-water depth in a delta sedimentation area without wells, aiming to solve the problem that the existing methods for analyzing paleo-water depth mostly rely on drilling data or are only applicable to large-scale basin conditions and cannot be applied to small-scale delta sedimentation areas without wells.

[0007] The present invention provides a method for compiling an isobath map of paleo-water depth in a delta sedimentation area without wells, including:

[0008] Determine the provenance direction and the range of delta lobes based on the research results of the geological background of the target area, and determine the upstream position and input point of the provenance input of the delta lobes;

[0009] Taking the input point as the origin, take multiple seismic profiles S passing through the delta lobes i and set the multiple seismic profiles S i at equal angular intervals, and determine the main provenance profile S0 of the seismic profile with the largest angle of the delta wedge;

[0010] On the seismic profile S i trace the top boundary and bottom boundary of the delta lobe, and determine the point corresponding to the position with the largest thickness of the delta lobe on the top boundary as point P i wherein, the point corresponding to the top boundary and the position with the largest thickness of the delta lobe on the main provenance profile S0 is point P0, the thickness of the delta lobe corresponding to point P0 is H0, and the thickness of the delta lobe is the vertical distance between the top boundary and the bottom boundary;

[0011] Connect point P0 and P in sequence i to obtain the isobath L0;

[0012] Take multiple points K along the top boundary of the delta lobe on the main provenance profile S0 j and measure the thickness value H of the delta lobe at point K j ; j ;

[0013] Based on the thickness values H0 and H of the delta lobe j calculate the corresponding water depth values W0 and W at point P0 and point K j ; j ;

[0014] On the main provenance profile S0, starting from point K j draw an isobath L parallel to the isobath L0 j, a paleo - water depth isoline map of the delta sedimentation area without wells is obtained. The water depth value corresponding to the water depth isoline L0 is W0, and the water depth isoline L j corresponds to a water depth value of W j .

[0015] According to the method for compiling the paleo - water depth isoline map of the delta sedimentation area without wells provided by the present invention, the angle between two adjacent seismic profiles is 10 degrees to 20 degrees.

[0016] According to the method for compiling the paleo - water depth isoline map of the delta sedimentation area without wells provided by the present invention, the angle between two adjacent seismic profiles is 15 degrees.

[0017] According to the method for compiling the paleo - water depth isoline map of the delta sedimentation area without wells provided by the present invention, multiple points K j are selected at equal intervals.

[0018] According to the method for compiling the paleo - water depth isoline map of the delta sedimentation area without wells provided by the present invention, multiple points K are taken along the top boundary of the delta lobe of the main provenance profile S0 j including:

[0019] Starting from the point P0, multiple points K are respectively selected at equal intervals upstream and downstream of the top boundary j .

[0020] According to the method for compiling the paleo - water depth isoline map of the delta sedimentation area without wells provided by the present invention, the water depth value corresponding to the point K located upstream of the point P0 j is calculated by the following formula:

[0021]

[0022] where W b is the water depth of the normal wave base obtained through literature research.

[0023] According to the method for compiling the paleo - water depth isoline map of the delta sedimentation area without wells provided by the present invention, the water depth value corresponding to the point K located downstream of the point P0 j is calculated by the following formula:

[0024]

[0025] where W b is the water depth of the normal wave base obtained through literature research. The present invention also provides a non - transitory computer - readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the method for compiling the paleo - water depth isoline map of the delta sedimentation area without wells as described above are implemented.

[0026] The present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the steps of the method for compiling the paleo-water depth isoline map of the well-free delta deposition area as described above are implemented.

[0027] The present invention also provides a computer program product, including a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the steps of the method for compiling the paleo-water depth isoline map of the well-free delta deposition area as described above.

[0028] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored on the memory and operable on the processor. When the processor executes the program, the steps of the method for compiling the paleo-water depth isoline map of the well-free delta deposition area as described above are implemented.

[0029] Due to the above technical solutions adopted by the present invention, it has the following advantages:

[0030] In the method for compiling the paleo-water depth isoline map of the well-free delta deposition area provided by the present invention, first, based on the research results of the geological background of the target area, the provenance direction and the range of the delta lobe are determined, and the upstream position and input point of the provenance input of the delta lobe are determined. Then, taking the input point as the origin, multiple seismic profiles S i , and the multiple seismic profiles are set at equal intervals of angles, and the seismic profile with the largest angle of the delta wedge is determined as the main provenance profile S0. Then, on the seismic profile S i , the top boundary and the bottom boundary of the delta lobe are traced, and the point corresponding to the position with the largest thickness of the delta lobe on the top boundary is determined as point P i . Among them, the point corresponding to the position with the largest thickness of the delta lobe on the main provenance profile S0 is point P0, the thickness of the delta lobe at point P0 is H0, and the thickness of the delta lobe is the vertical distance between the top boundary and the bottom boundary. Then, P0 and P are sequentially connected i to obtain the water depth isoline L0. Then, multiple points K are taken along the top boundary of the delta lobe of the main provenance profile S0 j , and the thickness value H of the delta lobe at point K j is measured. Then, based on the thickness values H0 and H of the delta lobe j , the corresponding water depth values W0 and W at point P0 and point K j are calculated. Finally, on the seismic profile S0, starting from point K j , a water depth isoline L parallel to the water depth isoline L0 is made j . j j, the paleo - water depth contour map of the delta sedimentary area without wells is obtained. The water depth value corresponding to the water depth contour L0 is W0, and the water depth value corresponding to the water depth contour L j is W j . The method for compiling the paleo - water depth contour map of the delta sedimentary area without wells provided by the present invention can determine the paleo - water depth of the delta sedimentary area without wells without relying on drilling to obtain data, and is not affected by the research scale. It can quickly calculate the paleo - water depth of the delta sedimentary area without wells. Under the condition of the delta sedimentary area without wells during the oil and gas exploration stage, the method has more accurate and efficient paleo - water depth restoration results, and plays an important role in the research of paleo - water depth. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0032] Figure 1 is a schematic flow chart of the method for compiling the paleo - water depth contour map of the delta sedimentary area without wells provided by an embodiment of the present invention;

[0033] Figure 2 is a schematic diagram of the thickness map of the delta sedimentary area, input points, seismic profile, and main provenance profile S0 provided by an embodiment of the present invention;

[0034] Figure 3 is a schematic diagram of the maximum thickness point and thickness value of the main provenance profile S0 provided by an embodiment of the present invention;

[0035] Figure 4 is a schematic diagram of the water depth contour L0 formed by connecting the maximum thickness value points of the seismic profile provided by an embodiment of the present invention;

[0036] Figure 5 is a schematic diagram of the profile thickness value points of the main provenance profile S0 and their corresponding water depth values provided by an embodiment of the present invention;

[0037] Figure 6 is a schematic diagram of the restoration of the paleo - water depth contour of the delta sedimentary area without wells provided by an embodiment of the present invention;

[0038] Figure 7 is a schematic diagram of the structure of an electronic device provided by an embodiment of the present invention.

[0039] Reference Signs:

[0040] 810: Processor; 820: Communication Interface; 830: Memory; 840: Communication Bus. Detailed implementation manners

[0041] To make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below with reference to the accompanying drawings in the present invention. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present invention without creative efforts shall fall within the protection scope of the present invention.

[0042] In the description of the present invention, it should be understood that the orientation or positional relationships indicated by the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", "axial", "radial", "circumferential", etc. are based on the orientation or positional relationships shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus should not be construed as limiting the present invention.

[0043] In addition, the terms "first" and "second" are only used for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include one or more of such features. In the description of the present invention, "a plurality" means two or more, unless otherwise specifically defined.

[0044] In the present invention, unless otherwise clearly defined and limited, the terms "mounted", "connected", "coupled", "fixed", etc. should be understood in a broad sense. For example, it may be a fixed connection, a detachable connection, or integrated; it may be a mechanical connection or an electrical connection; it may be directly connected or indirectly connected through an intermediate medium, and it may be the internal communication of two elements or the interaction relationship between two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0045] In the present invention, unless otherwise clearly specified or limited, the first feature being "on" or "under" the second feature may mean that the first and second features are in direct contact, or the first and second features are indirectly in contact through an intermediate medium. Moreover, the first feature being "above", "over" and "on top of" the second feature may mean that the first feature is directly above or obliquely above the second feature, or simply means that the horizontal height of the first feature is higher than that of the second feature. The first feature being "under", "beneath" and "underneath" the second feature may mean that the first feature is directly below or obliquely below the second feature, or simply means that the horizontal height of the first feature is less than that of the second feature.

[0046] In the description of this specification, the description with reference to terms such as "an embodiment", "some embodiments", "examples", "specific examples", or "some examples", etc. means that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described can be combined in a suitable manner in any one or more embodiments or examples. In addition, without contradiction, those skilled in the art can combine and combine the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0047] The method for compiling the paleo - water depth contour map in the non - well delta sedimentation area provided by the present invention first determines the provenance direction and the range of delta lobes based on the research results of the geological background of the target area, and determines the upstream position and input point of the provenance input of the delta lobes. Then, taking the input point as the origin, multiple seismic profiles S i passing through the delta lobes are taken, and the multiple seismic profiles are set at equal - interval angles, and the seismic profile with the largest angle of the delta wedge is determined as the main provenance profile S0. Then, on the seismic profile S i the top boundary and the bottom boundary of the delta lobe are traced, and the point corresponding to the position with the largest thickness of the delta lobe on the top boundary is determined as point P i Among them, the point corresponding to the position with the largest thickness of the delta lobe on the top boundary of the main provenance profile S0 is point P0, the thickness of the delta lobe at point P0 is H0, and the thickness of the delta lobe is the vertical distance between the top boundary and the bottom boundary. Then, P0 and P are connected in sequence i to obtain the water depth contour L0. Then, multiple points K j are taken along the top boundary of the delta lobe of the main provenance profile S0, and the thickness value H j of the delta lobe at point K j is measured. Then, based on the thickness values H0 and H j of the delta lobe, the points P0 and K jThe corresponding water depth values W0 and W at the position j . Finally, on the seismic profile S0, starting from point K j , draw a water depth contour line L parallel to the water depth contour line L0 j , and obtain the paleo-water depth contour map of the delta sedimentation area without wells. The water depth value corresponding to the water depth contour line L0 is W0, and the water depth value corresponding to the water depth contour line L j is W j . The method for compiling the paleo-water depth contour map of the delta sedimentation area without wells provided by the present invention can determine the paleo-water depth of the delta sedimentation area without wells without relying on drilling to obtain data, and is not affected by the research scale, and can quickly calculate the paleo-water depth of the delta sedimentation area without wells. This method has more accurate results and higher efficiency in restoring the paleo-water depth under the conditions of the delta sedimentation area without wells during the oil and gas exploration stage, and plays an important role in the research of paleo-water depth.

[0048] The following combines Figures 1 to 6 to describe the method for compiling the paleo-water depth contour map of the delta sedimentation area without wells provided by the present invention.

[0049] An embodiment of the present invention provides a method for compiling a paleo-water depth contour map of a delta sedimentation area without wells, including the following steps:

[0050] Step S100: Determine the provenance direction and the range of the delta lobe based on the research results of the geological background of the target area, and determine the upstream position and input point of the provenance input of the delta lobe;

[0051] Step S200: Take multiple seismic profiles S passing through the delta lobe with the input point as the origin i , and the multiple seismic profiles S i are set at equal angular intervals, and determine the main provenance profile S0 of the seismic profile with the largest angle of the delta wedge;

[0052] Step S300: Trace the top boundary and bottom boundary of the delta lobe on the seismic profile S i , and determine the point corresponding to the position with the largest thickness of the delta lobe on the top boundary as point P i , where the point corresponding to the position with the largest thickness of the delta lobe on the main provenance profile S0 is point P0, the thickness of the delta lobe corresponding to point P0 is H0, and the thickness of the delta lobe is the vertical distance between the top boundary and the bottom boundary;

[0053] Step S400: Connect point P0 and P in sequence i to obtain the water depth contour line L0;

[0054] Step S500: Take multiple points K along the top boundary of the delta lobe of the main provenance profile S0 j , and measure point Kj The thickness value H of the delta lobe at j ;

[0055] Step S600: Based on the thickness values H0 and H of the delta lobe j , calculate the water depth values W0 and W corresponding to point P0 and point K j ; j ;

[0056] Step S700: On the main provenance profile S0, starting from point K j , draw a water depth contour line L parallel to the water depth contour line L0 to obtain the paleo-water depth contour map of the well-free delta sedimentation area. The water depth value corresponding to the water depth contour line L0 is W0, and the water depth value corresponding to the water depth contour line L j is W j ; j .

[0057] Specifically, first determine the provenance direction and the delta lobe range according to the regional geological background research results, and find the upstream position and input point of the delta lobe provenance input.

[0058] Then, taking the input point as the origin, at a certain angular interval, take a series of seismic profiles S i (i = 0, 1, 2,...) of the delta lobe. On the seismic profile, according to the angle of the delta wedge, determine the seismic profile with the largest angle of the delta wedge as the main provenance profile, denoted as the main provenance profile S0. The remaining seismic profiles can be sequentially denoted as S1, S2, S3, S4 in the clockwise direction starting from the main provenance profile, and marked as S5, S6, S7 in the counterclockwise direction.

[0059] Then, on the main provenance profile S0, trace the top boundary and bottom boundary of the delta lobe, calculate the height difference between the top boundary and bottom boundary at each position. This height difference is the thickness of the delta lobe. Denote the point at the maximum thickness position on the top boundary as point P0, and the maximum thickness value here is H0.

[0060] Similarly, on the profile Si (i = 1, 2,...), denote the point at the maximum thickness position on the top boundary as point P i , and connect point P i and point P0 in sequence to obtain the water depth contour line L0.

[0061] Then, take a series of points K j (j = 1, 2,...) on the main provenance profile S0. The thickness value of the delta lobe corresponding to point K j is H j (j = 1, 2,...). Using the relationship between the thickness value and the water depth value, calculate point P0 and point K jThe corresponding water depth values W0 and W at the location j (j = 1, 2,...).

[0062] Finally, on the main provenance profile S0, starting from the point K taken in the previous step j travel, make a water depth contour line L parallel to the water depth contour line L0 j (j =, 1, 2,...), L j is the paleo - water depth contour line during the delta lobe deposition period. The water depth value corresponding to the water depth contour line L0 is W0, and the water depth value corresponding to the water depth contour line L j is set as W j , and the paleo - water depth contour map of the delta sedimentation area without wells is obtained.

[0063] In some embodiments, the angle between two adjacent seismic profiles is 10 degrees to 20 degrees, specifically it can be 15 degrees.

[0064] In some embodiments, the point K j is equally spaced on the top boundary line of the delta lobe of the main provenance profile S0.

[0065] Furthermore, the method of taking multiple points K j along the top boundary line of the delta lobe of the main provenance profile S0 is to equally spacedly select multiple points K j upstream and downstream of the top boundary line starting from the point P0.

[0066] The water depth value corresponding to the point K j is calculated by the following formula:

[0067]

[0068] where W b is the normal wave base water depth obtained from literature research. Formula (1) represents the points located upstream of the point P0, and this formula is used to calculate the water depth value; formula (2) represents the points located downstream of the point P0, and this formula is used to calculate the water depth value.

[0069] In a specific embodiment, taking the area where the delta lobe is expected to develop in the Middle Miocene Sanya Formation of the Haikou 29A structure as an example, the method for compiling the paleo - water depth contour map of the delta sedimentation area without wells includes the following steps:

[0070] As Figure 2 shown, according to the results of regional geological background research, determine that the provenance direction is from the northeast direction. According to the thickness value between the top surface T52 and the bottom surface T60, find the upstream position and input point of the delta lobe provenance input.

[0071] In Figure 2In it, taking the input point as the origin, a series of seismic profiles S passing through the delta lobe are extracted in the 3D seismic work area at an interval of 15° included angle. i (i = 0, 1, 2,...). In the above seismic profiles, according to the included angle of the delta wedge, the seismic profile with the largest included angle of the delta wedge is determined as the main provenance profile S0.

[0072] As Figure 3 shown, on the main provenance profile S0, trace the top boundary and the bottom boundary of the delta lobe. According to the difference between the top boundary and the bottom boundary, find the position with the largest difference, and determine the point where the top boundary is located at this position as point P0. The maximum thickness value corresponding to point P0 is H0 = 176 m.

[0073] Similarly, in Figure 2 the seismic profile S i (i = 1, 2,...), trace the top boundary and the bottom boundary of the delta lobe, and find the position point P at the top boundary where the difference between the top boundary and the bottom boundary is the largest i . On the plane, connect point P0 and point P i in sequence to obtain the water depth isoline L0, as Figure 4 shown.

[0074] As Figure 5 shown, on the above main provenance profile S0, with point P0 as the center, set a series of points K j (j = 1, 2,...) at equal intervals upstream and downstream. The thickness values of the delta lobe corresponding to these points are H j (j = 1, 2,...). Use the calculation formula of the thickness value H j and the water depth value W j to calculate the corresponding water depth values at each point. The calculation formula is as follows:

[0075]

[0076] In the above formula, W b is the normal wave base water depth obtained from literature research, W b = 20 m. Formula (1) represents the points located upstream of point P0, and this formula is used to calculate the water depth value; formula (2) represents the points located downstream of point P0, and this formula is used to calculate the water depth value.

[0077] Through the above formula, calculate the water depth values W j of the delta lobe at point K j (j = 1, 2,...). Among them, the thickness value H0 corresponding to point P0 is 176 m, and the corresponding water depth value W0 is 10 m. The water depth value corresponding to the water depth isoline L0 is set to be 10 m.

[0078] On a plane, starting from point K set on the main source profile S0 j draw a water depth contour line L parallel to the water depth contour line L0 j and set the water depth value corresponding to the water depth contour line L j as W j L j is the paleo - water depth contour line during the delta lobe deposition period, and the water depth value corresponding to each curve is W j thus obtaining the paleo - water depth contour map of the delta sedimentation area without wells, as shown in Figure 6 shown

[0079] Due to the adoption of the above - mentioned technical solutions, the present invention has the following advantages:

[0080] The present invention avoids relying on data obtained from drilling wells

[0081] The present invention avoids calculation errors of the paleo - water depth of large - scale basins

[0082] Based on the relationship between the delta lobe thickness and the paleo - water depth, the present invention can restore the paleo - water depth of the delta sedimentation area without drilling wells

[0083] In summary, the present invention can be widely applied to the method and system for compiling the paleo - water depth contour map of the delta sedimentation area without drilling wells

[0084] Figure 7 Illustrates a schematic diagram of the physical structure of an electronic device, as shown in Figure 7 shown. The electronic device may include: a processor 810, a communication interface 820, a memory 830, and a communication bus 840. Among them, the processor 810, the communication interface 820, and the memory 830 complete communication with each other through the communication bus 840. The processor 810 can call the logical instructions in the memory 830 to execute the method for compiling the paleo - water depth contour map of the delta sedimentation area without wells, and this method includes:

[0085] Step S100: Determine the source direction and the range of the delta lobe based on the research results of the geological background of the target area, and determine the upstream position and input point of the delta lobe source input

[0086] Step S200: Take multiple seismic profiles S passing through the delta lobe with the input point as the origin i and the multiple seismic profiles S i are set at equal - interval angles, and determine the main source profile S0 of the seismic profile with the largest angle of the delta wedge

[0087] Step S300: On the seismic profile S iTrack the top and bottom boundaries of the delta lobe, and determine the point corresponding to the position with the maximum thickness of the delta lobe on the top boundary as point P. i Among them, the point corresponding to the position with the maximum thickness of the delta lobe on the main provenance profile S0 is point P0, the thickness of the delta lobe corresponding to point P0 is H0, and the thickness of the delta lobe is the vertical distance between the top and bottom boundaries.

[0088] Step S400: Connect point P0 and P in sequence. i Obtain the water depth isoline L0.

[0089] Step S500: Take multiple points K along the top boundary of the delta lobe on the main provenance profile S0. j And measure the thickness value H of the delta lobe at point K. j j ;

[0090] Step S600: Based on the thickness values H0 and H of the delta lobe. j Calculate the corresponding water depth values W0 and W at point P0 and point K. j j ;

[0091] Step S700: On the main provenance profile S0, starting from point K. j Make a water depth isoline L parallel to the water depth isoline L0. j Obtain the paleo-water depth isoline map of the delta sedimentation area without wells. The water depth value corresponding to the water depth isoline L0 is W0, and the water depth value corresponding to the water depth isoline L. j j ;

[0092] In addition, when the logical instructions in the above-mentioned memory 830 are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or a part of this technical solution can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0093] On the other hand, the present invention also provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the method for compiling the paleo-water depth isoline map in the non-well delta deposition area provided by the above-mentioned various methods. The method includes:

[0094] Step S100: Determine the provenance direction and the range of the delta lobe based on the research results of the geological background of the target area, and determine the upstream position and the input point of the provenance input of the delta lobe.

[0095] Step S200: Take multiple seismic profiles S passing through the delta lobe with the input point as the origin, i and the multiple seismic profiles S i are set at equal interval angles, and determine the main provenance profile S0 of the seismic profile with the largest angle of the delta wedge.

[0096] Step S300: Trace the top boundary and the bottom boundary of the delta lobe on the seismic profile S i , and determine the point P corresponding to the position with the largest thickness of the delta lobe on the top boundary. i Among them, the point corresponding to the position with the largest thickness of the delta lobe on the top boundary of the main provenance profile S0 is the point P0, the thickness of the delta lobe corresponding to the point P0 is H0, and the thickness of the delta lobe is the vertical distance between the top boundary and the bottom boundary.

[0097] Step S400: Connect the points P0 and P in sequence i to obtain the water depth isoline L0.

[0098] Step S500: Take multiple points K along the top boundary of the delta lobe of the main provenance profile S0 j , and measure the thickness value H of the delta lobe at the point K j ; j

[0099] Step S600: Based on the thickness values H0 and H of the delta lobe j , calculate the water depth values W0 and W corresponding to the points P0 and K j ; j

[0100] Step S700: On the main provenance profile S0, starting from the point K j , draw a water depth isoline L j parallel to the water depth isoline L0 to obtain the paleo-water depth isoline map of the non-well delta deposition area. The water depth value corresponding to the water depth isoline L0 is W0, and the water depth value corresponding to the water depth isoline L j is W j .

[0101] In another aspect, the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it is configured to execute the method for compiling the paleo-water depth contour map of the wellless delta deposition area provided above, and the method includes:

[0102] Step S100: Determine the provenance direction and the range of the delta lobe based on the research results of the geological background of the target area, and determine the upstream position and the input point of the provenance input of the delta lobe;

[0103] Step S200: Take multiple seismic profiles S passing through the delta lobe with the input point as the origin i and the multiple seismic profiles S i are set at equal interval angles, and determine the main provenance profile S0 as the seismic profile with the largest angle of the delta wedge;

[0104] Step S300: Trace the top boundary and the bottom boundary of the delta lobe on the seismic profile S i , and determine the point corresponding to the position with the largest thickness of the delta lobe on the top boundary as point P i . Among them, the point corresponding to the position with the largest thickness of the delta lobe on the top boundary of the main provenance profile S0 is point P0, the thickness of the delta lobe corresponding to point P0 is H0, and the thickness of the delta lobe is the vertical distance between the top boundary and the bottom boundary;

[0105] Step S400: Connect point P0 and P in sequence i to obtain the water depth contour line L0;

[0106] Step S500: Take multiple points K along the top boundary of the delta lobe of the main provenance profile S0 j , and measure the thickness value H of the delta lobe at point K j ; j

[0107] Step S600: Based on the thickness values H0 and H of the delta lobe j , calculate the water depth values W0 and W corresponding to point P0 and point K j ; j

[0108] Step S700: On the main provenance profile S0, starting from point K j , draw a water depth contour line L j parallel to the water depth contour line L0 to obtain the paleo-water depth contour map of the wellless delta deposition area. The water depth value corresponding to the water depth contour line L0 is W0, and the water depth value corresponding to the water depth contour line L j is W j . ​​

[0109] The device embodiments described above are merely illustrative. 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 to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art can understand and implement it without creative efforts.

[0110] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the essence of the above technical solution, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0111] Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present application, and are not intended to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements on some or all of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of each embodiment of the present application, and they should all be covered by the scope of the claims and the specification of the present application. In particular, as long as there is no structural conflict, the technical features mentioned in each embodiment can be combined in any way. The present application is not limited to the specific embodiments disclosed in the text, but includes all technical solutions falling within the scope of the claims.

Claims

1. A method for compiling paleo - water depth contour maps in a non - well - delta sedimentation area, characterized in that Including: Determine the provenance direction and the range of delta lobes based on the research results of the geological background of the target area, and determine the upstream position and input points of the provenance input of the delta lobes; Taking the input point as the origin, multiple seismic profiles S passing through the delta lobe are taken i , and multiple seismic profiles S i are set at equal interval angles, and the seismic profile with the largest angle of the delta wedge is determined as the main source profile S0; On the seismic profile S i Trace the top boundary and the bottom boundary of the delta lobe on it, and determine the point corresponding to the position where the thickness of the delta lobe is the largest as point P i , where the point corresponding to the position where the top boundary and the thickness of the delta lobe are the largest on the main provenance profile S0 is point P0, the thickness of the delta lobe corresponding to point P0 is H0, and the thickness of the delta lobe is the vertical distance between the top boundary and the bottom boundary; Connect the point P0 and the point P in sequence i to obtain an isobath L0; Take a plurality of points K along the top boundary of the delta lobe along the main source profile S0 j , and measure the thickness value H of the delta lobe at the point K j ; j ; Based on the thickness values H0 and H of the delta lobe j , calculate the water depth values W0 and W corresponding to the point P0 and the point K j ; j ; On the main source profile S0, starting from the point K j draw a water depth contour line L parallel to the water depth contour line L0 j to obtain the paleo-water depth contour map of the delta sedimentation area without wells. The water depth value corresponding to the water depth contour line L0 is W0, and the water depth value corresponding to the water depth contour line L j is W j .

2. The method for compiling the paleo-water depth contour map of the subaqueous delta sedimentation area according to claim 1, wherein The angle between two adjacent seismic profiles is 10 degrees to 20 degrees.

3. The method for compiling the paleo-water depth contour map of the non-well delta sedimentation area according to claim 2, wherein, The angle between two adjacent seismic profiles is 15 degrees.

4. The method for compiling the paleo-water depth isoline map of the well-free delta sedimentation area according to claim 1, wherein Multiple points K j Selected at equal intervals.

5. The method for compiling the paleo-water depth contour map of the non-well delta sedimentation area according to claim 4, characterized in that, Take a plurality of points K along the top boundary of the delta lobe along the main sediment source profile S0 j including: Starting from the point P0, a plurality of the points K are selected at equal intervals upstream and downstream of the top boundary respectively. j .

6. The method for compiling the paleo-water depth contour map of the non-well delta sedimentation area according to claim 5, characterized in that, The point K upstream of the point P0 j The corresponding water depth value is calculated by the following formula: Among them, W b is the water depth of the normal wave base obtained from literature research.

7. The method for compiling the paleo - water depth isoline map of the delta sedimentation area without wells according to claim 5, characterized in that, The point K downstream of the point P0 j The corresponding water depth value is calculated by the following formula: Among them, W b is the water depth of the normal wave base obtained from literature research.

8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method for compiling the paleo-water depth contour map of the well-free delta sedimentary area according to any one of claims 1 to 7.

9. A computer program product, characterized in that, Including a computer program stored on a non-transitory computer-readable storage medium, the computer program includes program instructions, and when the program instructions are executed by a computer, the computer can execute the steps of the method for compiling the paleo-water depth contour map of the well-free delta sedimentary area according to any one of claims 1 to 7.

10. An electronic device, characterized in that, Including a memory, a processor, and a computer program stored on the memory and executable on the processor, and when the processor executes the program, it implements the steps of the method for compiling the paleo-water depth contour map of the well-free delta sedimentary area according to any one of claims 1 to 7.

Images