An isochronal seismic slice acquisition method, reservoir prediction method and device

Abstract

The application discloses an isochronous seismic slice acquisition method, a reservoir prediction method and device. The isochronous seismic slice acquisition method comprises the following steps: acquiring position information and a drift amount of a reference point and a reference point in a target layer slice of a seismic data body; for any point in the target layer slice of the seismic data body except the reference point and the reference point, the drift amount of the point is determined according to the position information of the point and the position information and the drift amount of the reference point and the reference point; and the isochronous seismic slice is obtained by drifting each point in the target layer slice in the seismic data body along a depth direction according to the drift amount of the point. The isochronous seismic slice can be reasonably cut from the seismic data body, and data basis is provided for realizing reasonable reservoir prediction.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil exploration and development, and particularly relates to an isochronal seismic slice acquisition method, a reservoir prediction method and device. BACKGROUND

[0002] In recent years, the lithologic strata oil and gas reservoirs account for more than 50% of the newly added proven reserves in China each year, and the proportion is increasing year by year. Whether in the large depression basins (Songliao Basin and Ordos Basin) and the faulted basins (Bohaiwan Basin) in the east and west of China, or in the foreland basins (Junggar Basin) in the west, various types of lithologic strata oil and gas reservoirs with abundant reserves have been discovered. However, the single layer thickness of the oil and gas layers contained in these reservoirs is generally less than the seismic vertical resolution, which brings great difficulties to seismic prediction.

[0003] In view of the technical bottleneck of seismic data resolution, many scholars have carried out inversion algorithm research from the perspective of improving the vertical resolution, so as to predict the thickness and distribution of thin reservoirs, and various inversion methods such as sparse reflection coefficient inversion, spectral inversion, waveform indication inversion and phase-controlled random inversion have been formed. Based on the tuning theory, some scholars have studied the relationships between thin layer spectrum, thin layer resolution, thin layer peak frequency and thickness, and the relationships between thin interbed tuning and resolution, and thin layer prediction methods such as spectral decomposition, peak instantaneous frequency attribute and frequency-to-amplitude ratio attribute have been formed. Zeng Hongliu and Posamentier have respectively proposed the concepts of seismic sedimentology and seismic geomorphology from the perspective of using the lateral resolution of seismic data. The core idea is to detect thin layers through stratal slice scanning. The so-called stratal slice is a slice obtained by proportional interpolation between two isochronal marker layers. It can be seen from the principle that the stratal slice is mainly suitable for areas where the lateral sedimentation rate remains stable. However, the sedimentation rate often changes greatly in different regions of a continental basin. The traditional stratal slice is easy to be out of time, which leads to the failure of thin sand body prediction. SUMMARY

[0004] In view of the above problems, the present application is proposed to provide an isochronal seismic slice acquisition method, a reservoir prediction method and device which can overcome the above problems or at least partially solve the above problems, and can reasonably cut the isochronal seismic slice from the seismic data volume, thereby providing data basis for realizing reasonable reservoir prediction.

[0005] In a first aspect, an embodiment of the present application provides an isochronal seismic slice acquisition method, comprising:

[0006] obtaining position information and drift of a reference point and a datum point in a target layer slice of seismic data volume, the drift of the datum point being 0, and the drift of the reference point not being equal to 0;

[0007] According to the position information of the point and the position information and the drift amount of the reference point, the drift amount of the point is determined.

[0008] Each point in the target layer slice is drifted in the depth direction in the seismic data volume according to the drift amount of the point, so as to obtain an isochronous seismic slice.

[0009] In a second aspect, an embodiment of the present application provides a reservoir prediction method, comprising:

[0010] A second seismic attribute volume of a set attribute is extracted from the seismic data volume, and a second seismic attribute slice at an isochronous seismic slice is extracted from the second seismic attribute volume, the isochronous seismic slice being cut from the seismic data volume according to the method above;

[0011] According to the reservoir characteristics of the well target layer, points corresponding to the well in the second seismic attribute slice are marked with reservoir characteristics;

[0012] According to the second seismic attribute slice and the reservoir characteristic marking, reservoir characteristics are predicted.

[0013] In a third aspect, an embodiment of the present application provides an isochronous seismic slice acquisition device, comprising:

[0014] A reference data acquisition module is configured to acquire position information and drift amounts of a reference point and a reference point in a target layer slice of a seismic data volume, the drift amount of the reference point being 0, and the drift amount of the reference point not being equal to 0;

[0015] A drift amount determination module is configured to determine, for any point in the target layer slice of the seismic data volume except the reference point and the reference point, a drift amount of the point according to position information of the point and the position information and the drift amount of the reference point and the reference point;

[0016] A drift module is configured to drift each point in the target layer slice in the depth direction in the seismic data volume according to the drift amount of the point, so as to obtain an isochronous seismic slice.

[0017] In a fourth aspect, an embodiment of the present application provides a reservoir prediction device, comprising:

[0018] An extraction module is configured to extract a second seismic attribute volume of a set attribute from a seismic data volume, and extract a second seismic attribute slice at an isochronous seismic slice from the second seismic attribute volume, the isochronous seismic slice being cut from the seismic data volume according to the method above;

[0019] A marking module is configured to mark points corresponding to a well in the second seismic attribute slice with reservoir characteristics according to reservoir characteristics of a well target layer;

[0020] a prediction module configured to predict reservoir characteristics according to the second seismic attribute slice and the reservoir characteristic label.

[0021] In a fifth aspect, an embodiment of the present application provides a computer program product with reservoir prediction function, comprising computer programs / instructions, wherein the computer programs / instructions are executed by a processor to implement the above-mentioned isochronous seismic slice acquisition method or implement the above-mentioned reservoir prediction method.

[0022] The above-mentioned technical solution provided by the embodiments of the present application has at least the following beneficial effects:

[0023] The isochronous seismic slice acquisition method provided by the embodiments of the present application determines the drift amount of any point in the target layer slice of the seismic data volume except the reference point and the reference point according to the position information of the point and the position information and the drift amount of the reference point, drifts each point in the target layer slice in the seismic data volume along the depth direction according to the drift amount of the point, and obtains the isochronous seismic slice. The time of the reference point is taken as a reference, and the drift amount of the reference point is simultaneously referred to to reasonably determine the drift amount of other points in the target layer slice, so that each point in the target layer slice is reasonably drifted, the isochronous seismic slice is cut in the seismic data volume, the time-penetrating property of the isochronous seismic slice is greatly reduced, and thin-layer reservoir prediction through seismic data becomes possible.

[0024] The technical solution of the present application is further described in detail below by means of the accompanying drawings and embodiments. BRIEF DESCRIPTION OF DRAWINGS

[0025] The accompanying drawings are used to provide a further understanding of the present application, and constitute a part of the specification, and are used to explain the present application together with embodiments of the present application, and do not constitute a limitation on the present application. In the accompanying drawings:

[0026] Figure 1 is a schematic diagram of an extracted seismic wavelet waveform;

[0027] Figure 2 is a flowchart of the isochronous seismic slice acquisition method in the first embodiment of the present application;

[0028] Figure 3 is a specific implementation flowchart of the reference point determination method in the first embodiment of the present application;

[0029] Figure 4 is a specific implementation flowchart of the reference point determination method in the first embodiment of the present application;

[0030] Figure 5 is a specific implementation flowchart of the isochronous seismic slice acquisition method in the first embodiment of the present application;

[0031] Figure 6A is an example diagram of a seismic profile;

[0032] Figure 6B Seismic amplitude attribute map of the target layer obtained by the conventional slicing method;

[0033] Figure 6C Seismic amplitude attribute map of the target layer obtained by the method of Example 1;

[0034] Figure 7 Schematic view of a palaeotopographic map;

[0035] Figure 8 Structure schematic view of the medium-time seismic slicing acquisition device in the embodiment of the present application;

[0036] Figure 9 Structure schematic view of the reservoir prediction device in the embodiment of the present application. DETAILED DESCRIPTION

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

[0038] It should be understood that the terms used in the present application merely describe particular embodiments and are not intended to limit the present application. In addition, for numerical ranges in the present application, it should be understood that each intermediate value between the upper limit and the lower limit of the range is specifically disclosed. Each smaller range within the range of the intermediate values and any other stated value or intermediate value within the stated range is also included in the present application. The upper limit and the lower limit of these smaller ranges can be included or excluded independently from the range.

[0039] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present application, preferred methods and materials are described. All documents mentioned in this specification are incorporated by reference to disclose and describe the methods and / or materials in connection with which the documents are cited. In case of conflict between the content of the specification and that of any document incorporated herein by reference, the content of the specification prevails.

[0040] In the description of the present application, it should be noted that the terms "comprising", "including", "having", "containing", and the like, are open-ended terms, i.e., meaning "including, but not limited to". In addition, the terms "first", "second", and "third" and the like, are merely used to describe the objects and do not indicate or imply relative importance.

[0041] In order to solve the problem that the thin sand body prediction fails due to the time penetration of the seismic data body slicing in the prior art, the embodiment of the present application provides an isochronous seismic slicing acquisition method, a reservoir prediction method and device, which can reasonably slice the isochronous seismic slice from the seismic data body, and provides data basis for realizing the reasonable reservoir prediction.

[0042] The basic data required by the embodiment of the present application and the processing thereof include:

[0043] 1. The original post-stack seismic data body.

[0044] 2. The establishment of the time-depth correspondence.

[0045] The synthetic record is made for the well with the logging curve in the study area, the synthetic record is compared with the post-stack seismic data channel beside the well, and the correspondence between the seismic reflection time and the depth is established.

[0046] 3. The phase conversion of the original post-stack seismic data body.

[0047] The original post-stack seismic data body can be converted into a plurality of post-stack seismic data bodies with different phases through Hilbert transform. The logging lithology indicating curve (such as the natural potential, the natural gamma, etc.) and the post-stack seismic data body are matched through the time-depth correspondence, and a post-stack seismic data body with higher matching degree is selected.

[0048] The applicant finds that the matching degree of the post-stack seismic data body with-90 degrees phase and the logging lithology curve is higher through multiple experiments and comparisons, so the post-stack seismic data body with-90 degrees phase can be selected as the seismic data body used in the embodiment.

[0049] For example, taking the isochronous seismic slice acquisition of the target layer in a certain study area as an example, according to the wavelet extracted from the well logging data and the wavelet beside the well, Figure 1 It can be known that the original post-stack seismic data body is close to zero phase, the polarity is positive, and the wavelet length is about 70 ms. The synthetic seismic record of the typical well (well A) shows that the target layer depth is located at about 3540 meters, and the corresponding seismic data time is 3810 milliseconds. The original post-stack seismic data body is converted into the post-stack seismic data body with-90 degrees phase through Hilbert transform.

[0050] 4. The upper reference layer and the lower reference layer.

[0051] The isochronous marker seismic event is selected as the reference layer. The marker seismic event with the same geological time is selected as the reference layer in the post-stack seismic data body according to the logging calibration result, and in the usual case, the maximum flood (lake) flooding surface or the laterally stable distributed coal seam and the like reflection continuous layer can be selected. Finally, the upper reference layer and the lower reference layer of the target layer are selected.

[0052] 5. A target layer slice of the stacked seismic data volume.

[0053] According to the synthetic record of the well, the depth value of the target layer and the well trajectory data, a marker point of the target layer is determined in the seismic data volume, and a plurality of marker points are obtained according to a plurality of wells; under the constraint of the upper reference layer and the lower reference layer, the target layer slice is cut in the stacked seismic data volume according to the plurality of marker points by a thickness equal proportion interpolation method.

[0054] 6. Extraction of a seismic attribute volume of a set attribute.

[0055] The seismic attribute volume of the set attribute is extracted from the stacked seismic data volume.

[0056] Preferably, the set attribute can be a seismic amplitude; or can be a frequency or a phase or other attributes. The seismic attribute volume can be a seismic attribute volume of one seismic attribute, or can be a plurality of seismic attribute volumes of a plurality of seismic attributes.

[0057] 7. Generation of a palaeotopographic map before deposition of the target layer.

[0058] The thickness distribution data between the target layer slice and the reference layer are determined, and the palaeotopographic map before deposition of the target layer is obtained by referring to the upper reference layer or the lower reference layer. Here, the thickness can be a time thickness or an actual thickness.

[0059] Embodiment one

[0060] The embodiment one of the present application provides an isochronous seismic slice acquisition method, and a flowchart thereof is shown in Figure 2 The method comprises the following steps:

[0061] Step S21: acquiring position information and a drift amount of a reference point and a reference point in a target layer slice of a seismic data volume.

[0062] The drift amount of the reference point is 0, and the drift amount of the reference point is not equal to 0.

[0063] In some embodiments, referring to Figure 3 The reference point is determined by the following steps:

[0064] Step S2111: extracting a first seismic attribute volume of a set attribute from the seismic data volume, and extracting a first seismic attribute slice at the target layer slice from the first seismic attribute volume.

[0065] The set attribute can be one or more of a seismic amplitude, a frequency and a phase.

[0066] Step S2112: according to the well target layer reservoir characteristics and the corresponding marker point, marking the reservoir characteristics of the point corresponding to the marker point of the well in the first seismic attribute slice.

[0067] The reservoir feature can be a lithofacies feature, a sedimentary microfacies feature, or other features.

[0068] Step S2113: determining a reference point from the first seismic attribute slice according to the reservoir feature mark.

[0069] Generally, a point with obvious reservoir feature is selected as the reference point.

[0070] In some embodiments, referring to Figure 4 As shown, the reference point is determined by the following steps:

[0071] Step S2121: determining thickness distribution data between the target layer slice and the reference layer.

[0072] The parameter layer is the upper reference layer or the lower reference layer.

[0073] Step S2122: determining the reference point and the drift amount of the reference point from the target layer slice according to the thickness distribution data.

[0074] The thickness contour map obtained according to the thickness distribution data can represent a palaeotopographic map of the target layer before deposition. With the palaeotopographic thickness at the reference point as a reference, if the palaeotopographic thickness at the reference point is greater than the palaeotopographic thickness at the reference point, the reference point drifts downward, otherwise, the reference point drifts upward. According to the difference between the palaeotopographic thickness at the reference point and the palaeotopographic thickness at the reference point, the drift amplitude of the reference point is determined. The drift amount of the reference point is determined according to the drift direction and the drift amplitude of the reference point.

[0075] On the basis of the target layer slice, the moving direction is determined by the palaeotopographic indication slice, and the optimal drift amount is determined by the change of the seismic geomorphologic feature, so that the isochronous slice is obtained and the time-traveling property is avoided.

[0076] Step S22: for any point in the target layer slice of the seismic data volume except the reference point and the reference point, the drift amount of the point is determined according to the position information of the point and the position information and the drift amount of the reference point and the reference point.

[0077] Specifically, the drift amount of the point can be determined according to the depth value of the point and the depth value and the drift amount of the reference point and the reference point.

[0078] The depth value can be a time-depth value, an actual depth value, or a relative depth value.

[0079] In some embodiments, the drift amount of the point can be determined according to the following formula:

[0080]

[0081] Wherein, AutoShift is the drift amount of the point, T AutoShift is the depth value of the point, T RpDepth value of the reference point, T manual Depth value of the reference point, S manual Drift amount of the reference point.

[0082] Step S23: Drifting each point in the target layer slice in the seismic data volume along the depth direction according to the drift amount of the point to obtain an isochronous seismic slice.

[0083] Taking the thickness of the paleotopography at the reference point as a reference, if the thickness of the paleotopography in other regions is greater than the thickness of the paleotopography at the reference point, the target layer slice at the point is drifted downward, otherwise, the target layer slice is drifted upward; and the drift amount is different with the thickness of the paleotopography. Finally, the isochronous seismic slice is obtained.

[0084] The isochronous seismic slice acquisition method provided by the embodiment one of the present application is used for determining the drift amount of any point in the target layer slice of the seismic data volume except the reference point and the reference point according to the position information of the point and the position information and the drift amount of the reference point, drifting each point in the target layer slice in the seismic data volume along the depth direction according to the drift amount of the point to obtain an isochronous seismic slice. The time of the reference point is taken as a reference, and the drift amount of the other point in the target layer slice is reasonably determined by referring to the drift amount of the reference point, so that each point in the target layer slice is reasonably drifted, the isochronous seismic slice is cut in the seismic data volume, the time-penetration of the isochronous seismic slice is greatly reduced, and it is possible to predict the thin layer reservoir by using the seismic data.

[0085] In some embodiments, referring to FIG. 1, Figure 5 The isochronous seismic slice acquisition method can include:

[0086] 1. Generation of the paleotopography before the deposition of the target layer; 2. Extraction of the standard point of the isochronous slice; and 3. Generation of the isochronous slice. The generation of the paleotopography before the deposition of the target layer 1 includes: establishing the time-depth correspondence 101, generating the -90-degree phase seismic data volume 102, selecting the isochronous marker seismic event as the reference layer 103, generating the target layer stratum slice 104, and generating the paleotopography map before the deposition of the target layer 105. The extraction of the standard point of the isochronous slice 2 includes: extracting the seismic attribute value at the target layer stratum slice 201, and selecting the reference point of the isochronous slice 201. The generation of the isochronous slice 3 includes: determining the drift direction of the stratum slice 301, determining the drift amount of the stratum slice, and generating the isochronous slice 302.

[0087] Taking the isochronous seismic slice acquisition of the target layer in a certain research area as an example, referring to FIG. 2, Figures 6A-6C The isochronous seismic slice acquisition method can include: Figure 6A is Figure 6B the seismic profile at the profile position in FIG. 2, Figure 6At1 and t2 in the figure are respectively the upper reference layer and the lower reference layer of the target layer, isochronous marker seismic events t1 and t2 in the upper and lower of the target layer are selected as the reference layer, and the target layer slice s1 is extracted under the constraint of t1 and t2, and further, the seismic amplitude attribute corresponding to s1 is extracted, see Figure 6B . Figure 6B In the figure, the geological features marked according to well A are consistent with the seismic attribute, proving that the extracted seismic amplitude attribute has certain credibility. Figure 6B In the eastern part, a large-scale river channel with a north-east to south-west trend can be seen, the middle and southern part of the river channel is clearly displayed, but the northern part is fuzzy, and there is a river channel fuzzy area. In the seismic amplitude attribute, the middle point B of the river channel where the display is clear is selected as the isochronous slice standard point.

[0088] Figure 7 is the paleotopographic map of the sedimentary period according to the target layer slice s1 and the upper reference layer t1, the moving direction of the slice is indicated according to the paleotopographic map, the original target layer slice is subjected to drift processing, and the corresponding seismic amplitude value is extracted, Figure 6C is the amplitude attribute map corresponding to the isochronous slice after the drift processing, and Figure 6C It can be seen that the river channel in the northern part is clearly displayed.

[0089] It can be seen that the above method solves the problem that the traditional stratum slice technology is easy to be out of time due to the large difference in the sedimentation rate of different regions of the continental sedimentary stratum, and the thin layer prediction result is inaccurate, and a new thin layer prediction technology is provided, which is a major breakthrough of the seismic sedimentology in the thin interbed prediction.

[0090] Based on the inventive concept of the application, the embodiments of the application further provide a reservoir prediction method, comprising:

[0091] extracting a second seismic attribute volume of a set attribute from the seismic data volume, extracting a second seismic attribute slice at an isochronous seismic slice from the second seismic attribute volume, the isochronous seismic slice being cut from the seismic data volume according to the above method; marking the reservoir characteristics of the points corresponding to the well in the second seismic attribute slice according to the reservoir characteristics of the well target layer; and predicting the reservoir characteristics according to the second seismic attribute slice and the reservoir characteristic marking.

[0092] The reservoir characteristics can be lithofacies characteristics or sedimentary microfacies characteristics, and the reservoir prediction of the thin layer sand body can be realized through the above method.

[0093] Based on the inventive concept of the application, the embodiments of the application further provide an isochronous seismic slice acquisition device, the structure of the device is as shown in Figure 8 , comprising:

[0094] The reference data obtaining module 81 is configured to obtain position information and drift amounts of reference points and a reference point in a target layer slice of a seismic data volume, wherein the drift amount of the reference point is 0, and the drift amount of the reference point is not equal to 0.

[0095] The drift amount determining module 82 is configured to determine a drift amount of any point in the target layer slice of the seismic data volume other than the reference points and the reference point according to position information of the point and the position information and the drift amounts of the reference points and the reference point.

[0096] The drift module 83 is configured to drift each point in the target layer slice in a depth direction of the seismic data volume according to the drift amount of the point, to obtain an isochronous seismic slice.

[0097] In some embodiments, the drift amount determining module 82 is specifically configured to determine the drift amount of the point according to a depth value of the point and depth values and drift amounts of the reference points and the reference point.

[0098] In some embodiments, the drift amount determining module 82 is specifically configured to determine the drift amount of the point according to a depth value of the point and depth values and drift amounts of the reference points and the reference point.

[0099] In some embodiments, the drift amount determining module 82 is specifically configured to determine the drift amount of the point according to the following formula:

[0100]

[0101] wherein AutoShift is the drift amount of the point, T AutoShift is the depth value of the point, T Rp is the depth value of the reference point, T manual is the depth value of the reference point, and S manual is the drift amount of the reference point.

[0102] In some embodiments, the reference data obtaining module 81 is further configured to obtain the target layer slice by the following steps:

[0103] determining a mark point of the target layer in the seismic data volume according to a synthetic record of a well, a depth value of the target layer, and well trajectory data, and obtaining a plurality of mark points according to a plurality of wells;

[0104] under the constraint of an upper reference layer and a lower reference layer, cutting the target layer slice in the seismic data volume according to the plurality of mark points by a thickness-proportional interpolation method.

[0105] In some embodiments, the reference data obtaining module 81 is specifically configured to determine the reference point by the following steps:

[0106] extracting a first seismic attribute volume of a set attribute from the seismic data volume, extracting a first seismic attribute slice at the target layer slice from the first seismic attribute volume; marking reservoir features of points corresponding to the well's marked points in the first seismic attribute slice according to the well's target layer reservoir features and corresponding marked points; determining a reference point from the first seismic attribute slice according to the reservoir feature marking.

[0107] In some embodiments, the reference data obtaining module 81 is specifically configured to determine the reference point by the following steps:

[0108] determining thickness distribution data between the target layer slice and a reference layer, the parameter layer being the upper reference layer or the lower reference layer; determining a reference point and a drift amount of the reference point from the target layer slice according to the thickness distribution data.

[0109] Based on the inventive concept of the present application, the embodiments of the present application further provide a reservoir prediction device, the structure of which is shown in Figure 9 The device comprises:

[0110] The extraction module 91 is configured to extract a second seismic attribute volume of a set attribute from the seismic data volume, and extract a second seismic attribute slice at an isochronal seismic slice from the second seismic attribute volume, the isochronal seismic slice being cut from the seismic data volume according to the above method.

[0111] The marking module 92 is configured to mark reservoir features of points corresponding to the well in the second seismic attribute slice according to the well's target layer reservoir features.

[0112] The prediction module 93 is configured to perform reservoir feature prediction according to the second seismic attribute slice and the reservoir feature marking.

[0113] As to the device in the above embodiments, the specific manner in which each module performs operations has been described in detail in the embodiments of the method, and thus will not be described in detail here.

[0114] Based on the inventive concept of the present application, the embodiments of the present application further provide a computer program product with a reservoir prediction function, comprising computer programs / instructions, wherein the computer programs / instructions are executed by a processor to implement the above isochronal seismic slice obtaining method or the above reservoir prediction method.

[0115] It should be understood that the specific order or hierarchy of steps in the processes disclosed is an example. Based upon design preferences, it should be understood that the specific order or hierarchy of steps in the processes can be rearranged while remaining within the scope of the present disclosure. The accompanying method claims present elements of the various steps in a sample order, and as such claims should not necessarily be construed as depended on the particular order or hierarchy of the steps presented.

[0116] In the detailed description above, various features are grouped together in single embodiments for the purpose of streamlining the disclosure. This method of disclosure, however, is not to be interpreted as reflecting a necessity that the claimed subject matter requires more features than the claims do. Rather, aspects of the application are to be construed in accordance with the principles of patent law that section 112, sixth paragraph, of 35 U.S.C. expressly abolishes the requirement that a claimed disclosure be set forth in a single "specific" embodiment. Accordingly, the appended claims as they may come to be amended based on allowance in the USPTO, supersede this detailed description, and the complete scope of the claimed application is to be determined by references made to the claims as they may come to be amended, alone or in light of their legal equivalents.

[0117] The foregoing description includes example of one or more embodiments. Of course, not all possible combinations of components or method steps are described in order to describe the above mentioned embodiments, but one of ordinary skill in the art will recognize that further combinations and permutations of the above described elements are possible. Thus, the above described embodiments are intended to cover all such modifications and permutations of the disclosure incorporated herein. In addition, description of the term "comprising" throughout the specification is to be read as meanin the same as the term "including" as that term is used in the claims. Furthermore, any term "or" as used in the specification or claims is to be interpreted as an "inclusive OR".

Claims

1. A method for obtaining isochronous seismic slices, characterized in that, include: Based on the well's synthetic record, target layer depth, and well trajectory data, the target layer's annotation points are determined in the seismic data volume, and multiple annotation points are obtained from multiple wells. Based on the multiple annotation points, a slice of the target layer is cut from the seismic data volume using a thickness proportional interpolation method. Extract a first seismic attribute body with set attributes from the seismic data body, and extract a first seismic attribute slice at the target layer slice from the first seismic attribute body; mark the points in the first seismic attribute slice corresponding to the well's labeling points with reservoir features according to the reservoir characteristics of the target layer and the corresponding labeling points. Based on reservoir feature markers, a reference point is determined from the first seismic attribute slice; Determine the thickness distribution data between the target layer slice and the reference layer, and determine the reference point and the drift amount of the reference point from the target layer slice based on the thickness distribution data; The position information and drift amount of the reference point and the reference point in the target layer slice are obtained. The drift amount of the reference point is 0, and the drift amount of the reference point is not equal to 0. For any point in the target layer slice of the seismic data volume, excluding the reference point and the benchmark point, determine the drift of the point based on the location information of the point and the location information and drift of the reference point and the benchmark point; Each point in the target layer slice is drifted along the depth direction in the seismic data volume according to the drift amount of that point to obtain an isochronous seismic slice.

2. The method as described in claim 1, characterized in that, The step of determining the drift amount of a point based on its location information, the location information of the reference point, and the drift amount of the benchmark point and the reference point specifically includes: The drift amount of the point is determined based on the depth value of the point and the depth values ​​and drift amount of the reference and benchmark points.

3. The method as described in claim 2, characterized in that, The step of determining the drift amount of a point based on its depth value and the depth values ​​and drift amounts of the reference and benchmark points specifically includes determining the drift amount of the point according to the following formula: Where AutoShift is the drift amount at that point, and T AutoShift T is the depth value of that point. Rp T is the depth value of the reference point. manual S is the depth value of the reference point. manual The drift amount of the reference point.

4. The method as described in claim 1, characterized in that, The step of extracting a target layer slice from the seismic data volume based on the multiple annotation points using a thickness-proportional interpolation method includes: Under the constraints of the upper and lower reference layers, the target layer slices are cut from the seismic data volume based on the multiple annotation points using a thickness proportional interpolation method.

5. The method as described in claim 1, characterized in that, The reference layer is either an upper reference layer or a lower reference layer.

6. The method according to any one of claims 1 to 5, characterized in that, The seismic data volume is a -90 degree phase-stacked seismic data volume.

7. A reservoir prediction method, characterized in that, include: Extract a second seismic attribute body with set attributes from the seismic data body, and extract a second seismic attribute slice at the isochronous seismic slice from the second seismic attribute body, wherein the isochronous seismic slice is cut from the seismic data body according to the method of any one of claims 1 to 6; Based on the reservoir characteristics of the target well, the points in the second seismic attribute slice corresponding to the well are marked with reservoir characteristics. Based on the second seismic attribute slice and reservoir feature markers, reservoir feature prediction is performed.

8. A device for obtaining isochronous seismic slices, characterized in that, include: The reference data acquisition module is used to determine the annotation points of the target layer in the seismic data volume based on the well's synthetic record, the target layer depth value, and the well trajectory data, and to obtain multiple annotation points based on multiple wells; Based on the multiple annotation points, a target layer slice is cut from the seismic data volume using a thickness proportional interpolation method; a first seismic attribute volume with set attributes is extracted from the seismic data volume, and a first seismic attribute slice at the target layer slice is extracted from the first seismic attribute volume; based on the reservoir characteristics of the target layer of the well and the corresponding annotation points, the points in the first seismic attribute slice corresponding to the annotation points of the well are marked with reservoir characteristics. Based on reservoir feature markers, a reference point is determined from the first seismic attribute slice; thickness distribution data between the target layer slice and the reference layer is determined, and a reference point and its drift are determined from the target layer slice based on the thickness distribution data; the position information and drift of the reference point and the reference point in the target layer slice are obtained, wherein the drift of the reference point is 0, and the drift of the reference point is not equal to 0. The drift determination module is used to determine the drift amount of any point in the target layer slice of the seismic data volume, excluding the reference point and the benchmark point, based on the location information of the point and the location information of the reference point and the drift amount. The drift module is used to drift each point in the target layer slice along the depth direction in the seismic data volume according to the drift amount of that point, so as to obtain isochronous seismic slices.

9. A reservoir prediction device, characterized in that, include: An extraction module is used to extract a second seismic attribute volume with set attributes from a seismic data volume, and to extract a second seismic attribute slice at a time-seismic slice from the second seismic attribute volume, wherein the time-seismic slice is cut from the seismic data volume according to the method described in any one of claims 1 to 6; The marking module is used to mark the reservoir features of the points in the second seismic attribute slice corresponding to the well based on the reservoir features of the target well layer. The prediction module is used to predict reservoir features based on the second seismic attribute slice and reservoir feature markers.

10. A computer program product with reservoir prediction function, comprising a computer program / instructions, wherein, When the computer program / instruction is executed by the processor, it implements the isochronous seismic slice acquisition method according to any one of claims 1 to 6, or the reservoir prediction method according to claim 7.

Images