Method and system for quantitative restoration of provenance area of ​​source-sink system

Abstract

The present invention provides a quantitative restoration method and system for the provenance area of ​​a source-sink system, and obtains a plurality of first overshoot points and a plurality of farthest overshoot points in the denudation area of ​​the source-sink system to be evaluated; determines the first envelope line and the second envelope line corresponding to the first overshoot point and the farthest overshoot point, respectively, and determines the maximum denudation area area and the minimum denudation area area of ​​the denudation area of ​​the source-sink system to be evaluated according to the first envelope line and the second envelope line; performs a line extension and intersection operation on the first overshoot point of the wing and the farthest overshoot point of the wing located on the two wings of the denudation area of ​​the source-sink system to be evaluated, so as to determine the intersection point of the extended line of the overshoot point connection line; and obtains the minimum denudation height and the maximum denudation height respectively. The present invention does not require drilling core samples for U-Pb and (U-Th) / He dating analysis, and can be applied to no-well / few-well working conditions based on the accurate prediction of the area and distribution of favorable oil and gas reservoirs, and does not require high experimental testing costs.

Description

Technical Field

[0001] The present invention belongs to the technical field of erosion analysis of provenance areas in source-sink systems, and in particular, relates to a quantitative restoration method and system for provenance areas in source-sink systems. Background Art

[0002] Oil and natural gas (Oil and Gas) are important strategic resources for the development of the national economy. Finding favorable oil and gas reservoirs (reservoirs) is one of the keys to the success of oil and gas exploration. "How to accurately restore the amount of sediment erosion in the denudation area of ​​the source-sink system, so as to accurately predict the area and distribution of favorable oil and gas reservoirs" is crucial in oil and gas exploration.

[0003] The existing "quantitative restoration method of the amount of erosion in the provenance area of ​​the source-sink system" mostly relies on the dual dating method of detrital zircon U-Pb and (U-Th) / He, combined with petrological analysis, to determine the possible provenance areas of igneous rocks, metamorphic rocks and sedimentary rocks. On this basis, the range of each provenance area is imported into the software to calculate the drainage area of ​​the source-sink system in different periods. This traditional "quantitative restoration method of the provenance area of ​​the source-sink system based on detrital zircon dating" is heavily dependent on detrital zircon samples, and requires the use of core samples extracted from drilling to perform U-Pb and (U-Th) / He dating analysis.

[0004] Due to the high cost of drilling operations, the available core samples are often very limited or even non-existent. In addition, U-Pb and (U-Th) / He dating analysis often require high experimental testing costs. Therefore, it is urgent to creatively invent a "source-sink system provenance area erosion quantitative recovery method" that does not rely on core sample dating. Summary of the invention

[0005] In view of the above-mentioned defects or shortcomings of the prior art, the present invention provides a method and system for quantitative restoration of the provenance area of ​​a source-sink system to solve the technical problem of high cost in restoring the amount of sediment erosion in the erosion area of ​​a source-sink system in the prior art.

[0006] In order to achieve the above object, the present invention provides a method for quantitatively restoring a provenance area of ​​a source-sink system, the method comprising:

[0007] Step S1, obtaining a plurality of first overrun points and a plurality of farthest overrun points in the denudation area of ​​the source-sink system to be evaluated;

[0008] Step S2, determining a first envelope and a second envelope corresponding to the first overrun point and the farthest overrun point respectively, and determining the maximum erosion area and the minimum erosion area of ​​the erosion area of ​​the source-sink system to be evaluated according to the first envelope and the second envelope;

[0009] Step S3, extending and intersecting the line connecting the first overshoot point of the wing and the far overshoot point of the wing located at the two wings in the denudation area of ​​the source-sink system to be evaluated, so as to determine the intersection point of the extended line connecting the overshoot points;

[0010] Step S4, connecting two far-up points of the wing to obtain a first line, connecting two front-up points of the wing to obtain a second line, and drawing perpendicular lines from the intersection points to the first line and the second line to obtain the minimum erosion height and the maximum erosion height respectively;

[0011] Step S5, quantitatively calculating the maximum erosion zone volume, the residual erosion zone volume and the sediment erosion volume in the erosion zone of the source-sink system to be evaluated according to the maximum erosion zone area, the minimum erosion zone area, the minimum erosion height and the maximum erosion height.

[0012] In the embodiment of the present invention, the step S1 includes:

[0013] Step S11, collecting three-dimensional seismic data of the denudation area of ​​the source-sink system to be evaluated, and obtaining a seismic profile according to the three-dimensional seismic data;

[0014] Step S12, obtaining a preset color depth threshold and a preset continuity threshold of the seismic profile, and determining a plurality of end points of the stratigraphic segment in the seismic profile according to the preset color depth threshold and the preset continuity threshold;

[0015] Step S13, determining the first overtaking point and the farthest overtaking point among the plurality of termination points from the bottom to the top of the seismic profile.

[0016] In the embodiment of the present invention, the step S1 further includes:

[0017] Step S14, based on the first overrun point and the farthest overrun point in the seismic profile, using a seismic reflection event similarity comparison method, identifying and comparing the remaining first overrun points and the farthest overrun points with similar seismic reflection event termination characteristics in the three-dimensional seismic data body.

[0018] In the embodiment of the present invention, the step S5 includes:

[0019] Step S6, acquiring a plurality of equally spaced areas to be evaluated of the seismic profile with a preset step length;

[0020] Step S7, obtaining the sediment erosion volume of the area to be evaluated, and obtaining the maximum sediment erosion volume, the minimum sediment erosion volume and the average sediment erosion volume according to the multiple sediment erosion volumes.

[0021] In the embodiment of the present invention, the number of the areas to be evaluated is greater than 25.

[0022] In this embodiment of the present invention, step S4 includes:

[0023] Step S41, determining a regional marker layer on the seismic profile according to a preset formation reflection feature, wherein the preset formation reflection feature includes formation read-through length information and color depth information;

[0024] Step S42, performing a layer flattening operation on the seismic profile along the regional marker layer, connecting two wing far-up points on the flattened seismic profile to obtain a first line, and connecting two wing first-up points to obtain a second line.

[0025] In the embodiment of the present invention, the maximum denudation zone volume, the residual denudation zone volume and the sediment denudation volume are quantitatively calculated using the following formula:

[0026] V max =(A max ×H max ) / 3;

[0027] V min =(A min ×H min ) / 3;

[0028] V ero =V max -V min ;

[0029] Among them, V max is the maximum erosion volume, V min is the residual erosion volume, V ero is the volume of sediment erosion, A max is the maximum denudation area, A min is the minimum denudation area, H max is the maximum erosion height, H min is the minimum erosion height.

[0030] In this embodiment of the present invention, step S2 includes:

[0031] Step S21, projecting the first overtaking point and the farthest overtaking point on a horizontal plane;

[0032] Step S22, sequentially connecting the plurality of first overshoot points after projection to obtain a first envelope, and sequentially connecting the plurality of the farthest overshoot points after projection to obtain a second envelope;

[0033] Step S23, obtaining the maximum erosion area and the minimum erosion area according to the plane ranges defined by the first envelope and the second envelope respectively.

[0034] The present invention also provides a source-sink system provenance quantitative recovery system, the source-sink system provenance quantitative recovery system comprising:

[0035] An image processing module is used to obtain a plurality of first overshoot points and a plurality of farthest overshoot points of the erosion zone of the source-sink system to be evaluated; determine a first envelope line and a second envelope line corresponding to the first overshoot point and the farthest overshoot point, respectively, and determine the maximum erosion zone area and the minimum erosion zone area of ​​the erosion zone of the source-sink system to be evaluated according to the first envelope line and the second envelope line; perform a line extension and intersection operation on the first overshoot points of the wing and the far overshoot points of the wing located at the two wings of the erosion zone of the source-sink system to be evaluated, so as to determine the intersection point of the extended line of the overshoot point connection line; connect the two far overshoot points of the wing to obtain a first line, connect the two first overshoot points of the wing to obtain a second line, draw perpendicular lines along the intersection points to the first line and the second line, respectively, and obtain the minimum erosion height and the maximum erosion height respectively;

[0036] The data processing module is used to quantitatively calculate the maximum erosion zone volume, the residual erosion zone volume and the sediment erosion volume in the erosion zone of the source-sink system to be evaluated based on the maximum erosion zone area, the minimum erosion zone area, the minimum erosion height and the maximum erosion height.

[0037] In an embodiment of the present invention, the image processing module is also used to: collect three-dimensional seismic data of the erosion area of ​​the source-sink system to be evaluated, and obtain a seismic profile based on the three-dimensional seismic data; obtain a preset color depth threshold and a preset continuity threshold of the seismic profile, and determine multiple termination points of the stratigraphic segment in the seismic profile according to the preset color depth threshold and the preset continuity threshold; determine the first overtaking point and the farthest overtaking point among the multiple termination points from the bottom to the top of the seismic profile.

[0038] Through the above technical solution, the quantitative recovery method of the source-sink system provenance area provided by the embodiment of the present invention has the following beneficial effects:

[0039] By determining the first envelope and the second envelope through multiple first overshoot points and multiple farthest overshoot points in the denudation area of ​​the source-sink system to be evaluated, the denudation area of ​​the source-sink system to be evaluated can be accurately obtained. The minimum denudation height and the maximum denudation height can also be quantitatively calculated based on the first overshoot points and farthest overshoot points on the wings located on the two wings of the denudation area of ​​the source-sink system to be evaluated in the length direction; it can be helpful to find favorable oil and gas reservoirs (reservoirs), so as to accurately predict the area and distribution of favorable oil and gas reservoirs. On the basis of quantitatively calculating the maximum denudation area, the minimum denudation area, the minimum denudation height and the maximum denudation height, the maximum denudation area volume, the residual denudation area volume and the sediment denudation volume in the denudation area of ​​the source-sink system to be evaluated can also be quantitatively calculated. According to the sediment denudation volume range generated by the denudation of the source-sink system denudation area, the area and distribution of favorable oil and gas reservoirs can be further accurately predicted, which has important practical significance for favorable oil and gas reservoirs (reservoirs) and oil and gas exploration. In addition, compared with the "quantitative restoration method of source-sink system provenance area based on detrital zircon dating" in the prior art, the present invention does not require drilling core samples for U-Pb and (U-Th) / He dating analysis. On the basis of accurately predicting the area and distribution of favorable oil and gas reservoirs, it can also be applied to conditions with no wells or few wells, and does not require high experimental testing costs, thereby reducing the cost of exploration and prediction.

[0040] Other features and advantages of the present invention will be described in detail in the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] The accompanying drawings are used to provide an understanding of the present invention and constitute a part of the specification. Together with the following specific embodiments, they are used to explain the present invention, but do not constitute a limitation of the present invention. In the accompanying drawings:

[0042] Figure 1 is a schematic flow chart of a method for quantitatively restoring a source-sink system provenance area according to a first embodiment of the present invention;

[0043] Figure 2 It is a schematic diagram for determining the first super-point and the farthest super-point according to an embodiment of the present invention;

[0044] Figure 3 is a schematic diagram of a maximum erosion area and a minimum erosion area according to an embodiment of the present invention;

[0045] Figure 4 is a schematic diagram of the minimum erosion height and the maximum erosion height according to an embodiment of the present invention;

[0046] Figure 5 It is a schematic diagram of seismic section division according to one embodiment of the present invention. DETAILED DESCRIPTION

[0047] The specific embodiments of the present invention are described in detail below in conjunction with the accompanying drawings. It should be understood that the specific embodiments described herein are only used to illustrate and explain the present invention, and are not used to limit the present invention.

[0048] The method for quantitatively restoring the provenance area of ​​a source-sink system according to the present invention is described below with reference to the accompanying drawings.

[0049] like Figure 1 As shown, in an embodiment of the present invention, a method for quantitatively restoring a provenance area of ​​a source-sink system is provided, and the method for quantitatively restoring a provenance area of ​​a source-sink system comprises:

[0050] Step S1, obtaining a plurality of first overrun points and a plurality of farthest overrun points in the denudation area of ​​the source-sink system to be evaluated;

[0051] Step S2, determining a first envelope and a second envelope corresponding to the first overrun point and the farthest overrun point respectively, and determining the maximum and minimum erosion areas of the erosion area of ​​the source-sink system to be evaluated according to the first and second envelopes;

[0052] Step S3, extending and intersecting the line connecting the first overshoot point of the wing and the far overshoot point of the wing located at the two wings in the denudation area of ​​the source-sink system to be evaluated, so as to determine the intersection point of the extended line connecting the overshoot points;

[0053] Step S4, connecting two far-up points of the wing to obtain a first line, connecting two front-up points of the wing to obtain a second line, and drawing perpendicular lines from the intersection points to the first line and the second line, respectively, to obtain the minimum erosion height and the maximum erosion height respectively;

[0054] Step S5, quantitatively calculating the maximum denudation zone volume, the residual denudation zone volume and the sediment denudation volume in the denudation zone of the source-sink system to be evaluated according to the maximum denudation zone area, the minimum denudation zone area, the minimum denudation height and the maximum denudation height.

[0055] In this embodiment, the first envelope and the second envelope are determined by multiple first overshoot points and multiple farthest overshoot points in the erosion area of ​​the source-sink system to be evaluated, so that the erosion area of ​​the source-sink system to be evaluated can be accurately obtained, and the minimum erosion height and the maximum erosion height can be quantitatively calculated based on the first overshoot points and the farthest overshoot points of the wings located on the two wings of the erosion area of ​​the source-sink system to be evaluated in the length direction; it is helpful to find favorable oil and gas reservoirs (reservoirs), so as to accurately predict the area and distribution of favorable oil and gas reservoirs. On the basis of quantitatively calculating the maximum erosion area, the minimum erosion area, the minimum erosion height and the maximum erosion height, the maximum erosion area volume, the residual erosion area volume and the sediment erosion volume in the erosion area of ​​the source-sink system to be evaluated can also be quantitatively calculated. According to the sediment erosion volume range generated by the erosion of the source-sink system erosion area, the area and distribution of favorable oil and gas reservoirs can be further accurately predicted, which has important practical significance for favorable oil and gas reservoirs (reservoirs) and oil and gas exploration. In addition, compared with the "quantitative restoration method of source-sink system provenance area based on detrital zircon dating" in the prior art, this embodiment does not require drilling core samples for U-Pb and (U-Th) / He dating analysis. On the basis of accurately predicting the area and distribution of favorable oil and gas reservoirs, it can also be applied to conditions with no wells / few wells, and does not require high experimental testing costs, thereby reducing the cost of exploration and prediction.

[0056] like Figure 3 As shown, in one embodiment, based on the use of the first envelope line and the second envelope line, the "maximum erosion zone area and its corresponding maximum erosion zone height" and the "minimum erosion zone area and its corresponding minimum erosion zone height" are determined.

[0057] In one embodiment, step S1 includes:

[0058] Step S11, collecting three-dimensional seismic data of the denudation area of ​​the source-sink system to be evaluated, and obtaining a seismic profile according to the three-dimensional seismic data;

[0059] Step S12, obtaining a preset color depth threshold and a preset continuity threshold of the seismic profile, and determining a plurality of end points of the stratigraphic segment in the seismic profile according to the preset color depth threshold and the preset continuity threshold;

[0060] Step S13, determining the first overrun point and the farthest overrun point among the multiple termination points from the bottom to the top of the seismic profile. Figure 2 As shown, in one embodiment, it is a schematic diagram of identifying the first on-override point and the farthest on-override point in the erosion area (Bonan low uplift) of the second section source-sink system of the Dongying Formation in the Bozhong Sag of the Bohai Bay Basin, wherein the solid triangle icon in the figure may represent the farthest on-override point, and the solid positive direction icon may represent the first on-override point.

[0061] In this embodiment, the lateral continuity of the seismic reflection event is determined specifically by the preset color depth threshold and the preset continuity threshold on the seismic profile, so that the "termination point" where the seismic reflection event terminates one by one from the bottom position to the high position, that is, the overrun point, can be quickly identified and determined. Looking from bottom to top along the erosion area of ​​the source-sink system to be evaluated, the overrun point at the bottom is the "first overrun point", and the overrun point at the top is the "farthest overrun point". The specific preset color depth threshold and the preset continuity threshold can be set according to specific usage requirements, and the present invention does not limit them. In this embodiment, based on the accurate acquisition of multiple overrun points by the preset color depth threshold and the preset continuity threshold, the first overrun point and the farthest overrun point can be quickly compared and acquired, thereby improving the recognition accuracy of the first overrun point and the farthest overrun point.

[0062] In another embodiment, step S1 further includes:

[0063] Step S14, based on the first overrun point and the farthest overrun point in the seismic profile, a seismic reflection event similarity comparison method is used to identify and compare other first overrun points and farthest overrun points with similar seismic reflection event termination characteristics in the three-dimensional seismic data body.

[0064] In this embodiment, the seismic reflection event similarity comparison method can be used to identify all first overrun points and farthest overrun points in the three-dimensional space in the entire three-dimensional seismic data body around the erosion area of ​​the source-sink system to be evaluated, and the identification can be carried out quickly and accurately.

[0065] In the embodiment of the present invention, step S5 includes:

[0066] Step S6, obtaining a plurality of equally spaced areas to be evaluated of the seismic profile with a preset step length;

[0067] Step S7, obtaining the sediment erosion volume of the area to be evaluated, and obtaining the maximum sediment erosion volume, the minimum sediment erosion volume and the average sediment erosion volume according to the multiple sediment erosion volumes.

[0068] Specifically, a series of equally spaced seismic profiles S across the denuded area of ​​the source-sink system to be evaluated are obtained with the same measurement step length, such as 2.5 km, 5 km, 7.5 km, etc. i (i=1,2,3,…n), the number of seismic profiles obtained must be no less than 25 in the erosion area of ​​the source-sink system to be evaluated, and the geological requirements for sample points of statistical data must be met. According to the above steps S1 to S5, the sediment erosion volume (Vero i ,i=1,2,3,…n) On this basis, the maximum sediment erosion volume (Vero max ), minimum sediment erosion volume (Vero min) and the average sediment erosion volume (Vero mean ) to accurately obtain the denuded sediment volume of the source-sink system provenance area to be evaluated. In one embodiment, the number of areas to be evaluated is greater than 25. In other embodiments, the number of areas to be evaluated can be set according to actual usage requirements.

[0069] In this embodiment of the present invention, step S4 includes:

[0070] Step S41, determining a regional marker layer on a seismic profile according to a preset formation reflection feature, wherein the preset formation reflection feature includes formation read-through length information and color depth information;

[0071] Step S42, performing a layer flattening operation on the seismic profile along the regional marker layer, connecting two wing far-up points on the flattened seismic profile to obtain a first line, and connecting two wing first-up points to obtain a second line.

[0072] In this embodiment, the regional marker layer on the seismic profile is first obtained through the formation reflection characteristics as the basis for the layer pulling operation, which can provide an accurate graphic basis for the subsequent first connection line and second continuity determination, and further optimize the accuracy of the area and distribution prediction of favorable oil and gas reservoirs.

[0073] In another embodiment, the maximum denudation zone volume, the residual denudation zone volume and the sediment denudation volume are quantitatively calculated using the following formula:

[0074] V max =(A max ×H max ) / 3;

[0075] V min =(A min ×H min ) / 3;

[0076] V ero =V max -V min ;

[0077] Among them, V max is the maximum erosion volume, V min is the residual erosion volume, V ero is the volume of sediment erosion, A max is the maximum denudation area, A min is the minimum denudation area, H max is the maximum erosion height, H min In this embodiment, the maximum denudation area volume, the residual denudation area volume and the sediment denudation volume can be quantitatively obtained through the above formula, which is helpful for predicting the area and distribution of favorable oil and gas reservoirs.

[0078] In one embodiment, step S2 includes:

[0079] Step S21, projecting the first overtaking point and the farthest overtaking point on a horizontal plane;

[0080] Step S22, sequentially connecting the multiple first-overlapping points after projection to obtain a first envelope, and sequentially connecting the multiple farthest-overlapping points after projection to obtain a second envelope;

[0081] Step S23, obtaining the maximum erosion area and the minimum erosion area according to the plane ranges defined by the first envelope and the second envelope respectively.

[0082] In this embodiment, the first overtaking point and the farthest overtaking point are firstly projected on the horizontal plane, and then the first envelope line and the second envelope line are obtained in sequence by connecting them. This not only facilitates the connection operation of the first overtaking point and the farthest overtaking point, but also makes the determination of the first envelope line and the second envelope line more accurate, thereby accurately and scientifically obtaining the erosion area of ​​the source area of ​​the source-sink system to be evaluated.

[0083] The present invention also proposes a source-sink system provenance quantitative recovery system, which includes an image processing module and a data processing module:

[0084] The image processing module is used to obtain multiple first overrun points and multiple farthest overrun points in the erosion area of ​​the source-sink system to be evaluated; determine the first envelope line and the second envelope line corresponding to the first overrun point and the farthest overrun point respectively, and determine the maximum erosion area and the minimum erosion area of ​​the erosion area of ​​the source-sink system to be evaluated according to the first envelope line and the second envelope line; extend and intersect the wing head overrun points and the wing far overrun points located on the two wings of the erosion area of ​​the source-sink system to be evaluated, so as to determine the intersection point of the extended line of the overrun point connection line; connect the two wing far overrun points to obtain the first line, connect the two wing head overrun points to obtain the second line, draw perpendicular lines to the first line and the second line along the intersection point, and obtain the minimum erosion height and the maximum erosion height respectively; the data processing module is used to quantitatively calculate the maximum erosion area volume, the residual erosion area volume and the sediment erosion volume in the erosion area of ​​the source-sink system to be evaluated according to the maximum erosion area, the minimum erosion area, the minimum erosion height and the maximum erosion height.

[0085] In this embodiment, the image processing module determines the first envelope and the second envelope through multiple first overshoot points and multiple farthest overshoot points in the denudation area of ​​the source-sink system to be evaluated, and can accurately obtain the denudation area of ​​the source-sink system provenance area to be evaluated, and can also quantitatively calculate the minimum denudation height and the maximum denudation height according to the first overshoot points and farthest overshoot points of the wings located on the two wings of the denudation area of ​​the source-sink system to be evaluated in the length direction; it is helpful to find favorable oil and gas reservoirs (reservoirs), so as to accurately predict the area and distribution of favorable oil and gas reservoirs. On the basis of quantitatively calculating the maximum denudation area, the minimum denudation area, the minimum denudation height and the maximum denudation height, the maximum denudation area volume, the residual denudation area volume and the sediment denudation volume in the denudation area of ​​the source-sink system to be evaluated can also be quantitatively calculated by the data processing module. According to the sediment denudation volume range generated by the denudation of the source-sink system denudation area, the area and distribution of favorable oil and gas reservoirs can be further accurately predicted, which has important practical significance for favorable oil and gas reservoirs (reservoirs) and oil and gas exploration. In addition, compared with the "quantitative restoration method of source-sink system provenance area based on detrital zircon dating" in the prior art, this embodiment does not require drilling core samples for U-Pb and (U-Th) / He dating analysis. On the basis of accurately predicting the area and distribution of favorable oil and gas reservoirs, it can also be applied to conditions with no wells / few wells, and does not require high experimental testing costs, thereby reducing the cost of exploration and prediction.

[0086] In an embodiment of the present invention, the image processing module is also used to: collect three-dimensional seismic data of the erosion area of ​​the source-sink system to be evaluated, and obtain a seismic profile based on the three-dimensional seismic data; obtain a preset color depth threshold and a preset continuity threshold of the seismic profile, and determine multiple termination points of the stratigraphic segment in the seismic profile based on the preset color depth threshold and the preset continuity threshold; determine the first overrun point and the farthest overrun point among the multiple termination points from the bottom to the top of the seismic profile.

[0087] like Figure 2 As shown, the source-sink system provenance quantitative recovery method in this embodiment is used to perform quantitative recovery calculation on the “source-sink system of the second segment of the Dongying Formation in the Oligocene in the Bozhong Sag of the Bohai Bay Basin”:

[0088] Obtain seismic profiles across the denuded area of ​​the source-sink system to be evaluated, identify and regionally trace and interpret the first and farthest onlap points;

[0089] Using 3D seismic data, we can obtain Figure 2 Shown is a seismic profile across the denuded area of ​​the source-sink system to be evaluated (Bonan Low Uplift).

[0090] According to the lateral continuity of the seismic reflection event on the obtained seismic profile, the "termination points" where the seismic reflection event terminates one by one from the low position to the high position are identified.

[0091] The end points at the bottom and top of the second member of the Dongying Formation defined by the dot-dash line (top) and the double-dot-dash line (bottom) are the first onlap points of the Bonan low uplift of the second member of the Dongying Formation in the Oligocene in the Bozhong Sag of the Bohai Bay Basin ( Figure 2 Points A and A` in the Figure 2 Points B and B` in the figure.

[0092] The entire 3D seismic data volume is encircled by the "Oligocene Dongying Formation 2nd Section Source-Sink System Denudation Area (Bonan Low Uplift) in Bozhong Sag" to identify and interpret all the first onset points and the farthest onset points in the 3D space. All the first onset points and the farthest onset points are projected onto the horizontal plane, and they are connected into lines on the horizontal plane to obtain the first envelope and the second envelope. The first envelope and the second envelope. The encircled plane ranges are the area of ​​the maximum denudation zone (A max =1596km 2 )( Figure 3 The area of ​​the midpoint dash line) and the area of ​​the minimum erosion zone (A min =462km 2 )( Figure 3 The area within the double-dotted line).

[0093] The seismic profile across the Bonan low uplift was flattened by performing layer-pulling operations along the regional marker layer of the second segment of the Oligocene Dongying Formation in the Bozhong Sag to obtain the flattened seismic profile ( Figure 4 On the flattened seismic profile, the first overburden point (first overburden point on the wing) and the farthest overburden point (far overburden point on the wing) of the outermost wing of the denudation zone of the source-sink system of the second segment of the Dongying Formation in the Oligocene in the Bozhong Sag were connected into lines and extended to intersect, and the intersection point of the extended lines of the overburden point connection was determined ( Figure 4 midpoint C).

[0094] Get the minimum erosion height (H min =571m), maximum erosion height (H max =854m).

[0095] The maximum erosion volume (V max ), residual erosion volume (V min ) and the volume of sediment erosion (V ero ):

[0096] V max =(1596km×0.571km) / 3=304km 3

[0097] V min =(462km×0.854km) / 3=132km 3

[0098] V ero =304–132=172km 3

[0099] A series of equally spaced seismic profiles S were obtained with a measurement step length of 5 km across the eroded area of ​​the second segment of the Dongying Formation in the Oligocene in the Bozhong Sag (Bonan Low Uplift). i (i=1,2,3,…n), flatten the top interface of each area to be evaluated and calculate the sediment erosion volume (V ero i ,i=1,2,3,…n) On this basis, the maximum sediment erosion volume (V ero max ), minimum sediment erosion volume (V ero min ) and the average sediment erosion volume (V ero mean ) are 247km3, 87km3 and 159km3 respectively. It can be seen that the volume of sediments produced by erosion in the Bonan Low Uplift of the Bohai Bay Basin during the deposition of the second member of the Dongying Formation in the Oligocene ranges from 87km3 to 247km3, with an average of 159km3.

[0100] In the description of the present invention, it should be understood that the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of the present invention, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0101] In the present invention, unless otherwise clearly specified and limited, the terms "installed", "connected", "connected", "fixed" and the like should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral one; it can be a mechanical connection, an electrical connection, or communication with each other; it can be a direct connection, or an indirect connection through an intermediate medium, it can be the internal connection of two elements or the interaction relationship between two elements, unless otherwise clearly defined. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

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

[0103] Although the embodiments of the present invention have been shown and described above, it is to be understood that the above embodiments are exemplary and are not to be construed as limitations of the present invention. A person skilled in the art may change, modify, replace and vary the above embodiments within the scope of the present invention.

Claims

1. A quantitative restoration method for the provenance of a source-sink system. It is characterized in that The method for quantitatively restoring the provenance area of ​​the source-sink system comprises: Step S1, obtaining a plurality of first overrun points and a plurality of farthest overrun points in the denudation area of ​​the source-sink system to be evaluated; Step S2, determining a first envelope and a second envelope corresponding to the first overrun point and the farthest overrun point respectively, and determining the maximum erosion area and the minimum erosion area of ​​the erosion area of ​​the source-sink system to be evaluated according to the first envelope and the second envelope; Step S3, extending and intersecting the line connecting the first overshoot point of the wing and the far overshoot point of the wing located at the two wings in the denudation area of ​​the source-sink system to be evaluated, so as to determine the intersection point of the extended line connecting the overshoot points; Step S4, connecting two far-up points of the wing to obtain a first line, connecting two front-up points of the wing to obtain a second line, and drawing perpendicular lines from the intersection points to the first line and the second line to obtain the minimum erosion height and the maximum erosion height respectively; Step S5, quantitatively calculating the maximum erosion zone volume, the residual erosion zone volume and the sediment erosion volume in the erosion zone of the source-sink system to be evaluated according to the maximum erosion zone area, the minimum erosion zone area, the minimum erosion height and the maximum erosion height.

2. The method for quantitative restoration of source-sink system provenance according to claim 1, It is characterized in that The step S1 comprises: Step S11, collecting three-dimensional seismic data of the denudation area of ​​the source-sink system to be evaluated, and obtaining a seismic profile according to the three-dimensional seismic data; Step S12, obtaining a preset color depth threshold and a preset continuity threshold of the seismic profile, and determining a plurality of end points of the stratigraphic segment in the seismic profile according to the preset color depth threshold and the preset continuity threshold; Step S13, determining the first overtaking point and the farthest overtaking point among the plurality of termination points from the bottom to the top of the seismic profile.

3. The method for quantitative restoration of source-sink system provenance area according to claim 2, It is characterized in that The step S1 further comprises: Step S14, based on the first overrun point and the farthest overrun point in the seismic profile, using a seismic reflection event similarity comparison method, identifying and comparing the remaining first overrun points and the farthest overrun points with similar seismic reflection event termination characteristics in the three-dimensional seismic data body.

4. The method for quantitative restoration of source-sink system provenance area according to claim 2, It is characterized in that The step S5 then includes: Step S6, obtaining a plurality of equally spaced areas to be evaluated of the seismic profile with a preset step length; Step S7, obtaining the sediment erosion volume of the area to be evaluated, and obtaining the maximum sediment erosion volume, the minimum sediment erosion volume and the average sediment erosion volume according to the multiple sediment erosion volumes.

5. The method for quantitative restoration of source-sink system provenance area according to claim 4, It is characterized in that The number of the areas to be evaluated is greater than 25.

6. The method for quantitative restoration of source-sink system provenance according to any one of claims 2 to 5, It is characterized in that The step S4 comprises: Step S41, determining a regional marker layer on the seismic profile according to a preset formation reflection feature, wherein the preset formation reflection feature includes formation read-through length information and color depth information; Step S42, performing a layer flattening operation on the seismic profile along the regional marker layer, connecting two wing far-up points on the flattened seismic profile to obtain a first line, and connecting two wing first-up points to obtain a second line.

7. The method for quantitative restoration of source-sink system provenance according to any one of claims 1 to 5, It is characterized in that The maximum denudation area volume, the residual denudation area volume and the sediment denudation volume are quantitatively calculated using the following formula: In max =(A max ×H max ) / 3; In min =(A min ×H min ) / 3; V ero =V max -V min ; Among them, V max is the maximum erosion volume, V min is the residual erosion volume, V ero is the volume of sediment erosion, A max is the maximum denudation area, A min is the minimum denudation area, H max is the maximum erosion height, H min is the minimum erosion height.

8. The method for quantitative restoration of source-sink system provenance according to any one of claims 1 to 5, It is characterized in that The step S2 comprises: Step S21, projecting the first overtaking point and the farthest overtaking point on a horizontal plane; Step S22, sequentially connecting the plurality of first overshoot points after projection to obtain a first envelope, and sequentially connecting the plurality of the farthest overshoot points after projection to obtain a second envelope; Step S23, obtaining the maximum erosion area and the minimum erosion area according to the plane ranges defined by the first envelope and the second envelope respectively.

9. A quantitative restoration system for the provenance of source-sink systems. It is characterized in that The source-sink system provenance area quantitative recovery system comprises: An image processing module is used to obtain a plurality of first overshoot points and a plurality of farthest overshoot points of the erosion zone of the source-sink system to be evaluated; determine a first envelope line and a second envelope line corresponding to the first overshoot point and the farthest overshoot point, respectively, and determine the maximum erosion zone area and the minimum erosion zone area of ​​the erosion zone of the source-sink system to be evaluated according to the first envelope line and the second envelope line; perform a line extension and intersection operation on the first overshoot points of the wing and the far overshoot points of the wing located at the two wings of the erosion zone of the source-sink system to be evaluated, so as to determine the intersection point of the extended line of the overshoot point connection line; connect the two far overshoot points of the wing to obtain a first line, connect the two first overshoot points of the wing to obtain a second line, draw perpendicular lines along the intersection points to the first line and the second line, respectively, and obtain the minimum erosion height and the maximum erosion height respectively; The data processing module is used to quantitatively calculate the maximum erosion zone volume, the residual erosion zone volume and the sediment erosion volume in the erosion zone of the source-sink system to be evaluated based on the maximum erosion zone area, the minimum erosion zone area, the minimum erosion height and the maximum erosion height.

10. The source-sink system provenance quantitative recovery system according to claim 9, It is characterized in that The image processing module is also used to: collect three-dimensional seismic data of the erosion area of ​​the source-sink system to be evaluated, and obtain a seismic profile based on the three-dimensional seismic data; obtain a preset color depth threshold and a preset continuity threshold of the seismic profile, and determine multiple termination points of the stratigraphic segment in the seismic profile based on the preset color depth threshold and the preset continuity threshold; determine the first overtaking point and the farthest overtaking point among the multiple termination points from the bottom to the top of the seismic profile.

Images