Method for picking up bead string energy of carbonate reservoir
By constructing the root mean square amplitude attribute data volume and ratio under different window size conditions and setting threshold values, the automatic pickup of beaded energy in carbonate reservoirs is achieved, and the beaded amplitude energy and stratigraphic reflection in earthquake data is solved The problem of energy aliasing provides a detailed portrayal and identification of beads at small and medium scales.
Patent Information
- Application Number
- CN202310738529.2
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2023-06-20
- Publication Date
- 2024-12-20
AI Technical Summary
In marine carbonate reservoirs, the beaded amplitude energy in conventional seismic data and the seismic reflection energy of conventional formations are aliased, resulting in beaded anomalies of small and medium-scale beaded at easy to be ignored and difficult to effectively identify on the seismic profile.
By constructing the root mean square amplitude attribute data body under different window size conditions, calculating its ratio, setting the threshold value to identify the beading position and energy, the automatic pickup of beading energy is achieved. The specific steps include inputting post-stack seismic data, presetting the first and second processing windows, building a root mean square amplitude attribute data body, calculating the ratio value, setting the threshold value, building a beaded recognition ratio data body based on the ratio, and picking out the reservoir Beading position and energy.
Effectively identify and extract beaded anomalies in carbonate reservoirs, and finely portray the small-scale beading characteristics, solving the identification problem caused by the aliasing of beaded amplitude energy and the formation reflection energy in conventional seismic data, which is deep heterogeneity. The reservoir provides a fine data foundation.
Smart Images
Figure CN119165528A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to a seismic exploration data processing technology, and in particular to a method for automatically extracting the energy of beaded abnormal amplitude in a marine carbonate reservoir. Background Art
[0002] The statements in this section merely provide background information related to the present disclosure and may not constitute prior art.
[0003] In recent years, the exploration and development results of deep fractured-cave carbonate oil and gas reservoirs in the Tarim Basin, Sichuan Basin and other places, especially limestone karst oil reservoirs, have shown that most of the main oil and gas enrichment areas are closely related to deep and large faults. Oil and gas are enriched in deep and large fault fracture zones connected to oil sources. The nearby karst development areas have become an important field for exploring karst oil reservoirs. The formation mechanism is that underground karst water seeps along deep and large fault zones, dissolves and transforms the fracture fracture zones to form a karst fracture-cave system with a complex three-dimensional spatial structure, which is a typical fault-controlled karst; oil and gas migrate along the faults to the dissolution and fracture zones, and gather to form a special type of oil and gas reservoir under the sealing of the overlying marlstone or lateral dense limestone.
[0004] A large number of fault-karst oil reservoir exploration and development practices have shown that the development pattern of fault-karst reservoirs is characterized by a "five-part structure" with the oil source fault as the center and a laterally symmetrical structure, namely, the "surrounding rock-fracture zone-main oil source fault zone-fracture zone-surrounding rock" pattern. The reservoir types, logging and drilling characteristics in different zones are different. The main oil source fault zone is dominated by chaotic weak reflections, chaotic strong reflections and relatively regular "beaded" reflection characteristics.
[0005] There are two main types of geological origins of beads: 1) Fracture-cavity carbonate reservoirs are mainly formed by typical karst surface freshwater dissolution, with the main storage space being large caves, dissolution pores and cracks. The seismic anomalies are mainly "beads" distributed along the unconformity surface of karst residual hills and ancient water systems; 2) Fault-controlled fracture-cavity carbonate reservoirs are mainly formed by the combined action of multi-stage active strike-slip faults and superimposed later fluid transformation. The storage space is mainly caves (fault cavities), fracture zones and pores, with extremely strong heterogeneity. The seismic anomalies are mainly "beads", "chaotic reflections" and "linear weak reflections" distributed along high-steep strip-shaped strike-slip fault fracture zones of different scales.
[0006] The reservoirs encountered in the actual drilling are cave-type and fracture-hole-type. The drilling process shows emptying and leakage. The broken zone is mainly characterized by chaotic reflection and weak reflection. The reservoirs encountered in the drilling are mainly fracture-hole-type. The drilling emptying and leakage rate is lower than that in the main fault zone of the oil source. There is basically no abnormal reflection in the surrounding rock zone, and the actual drilling reservoir is not developed. This type of oil and gas reservoir is deeply buried and highly heterogeneous, and it is very difficult to identify it in detail through seismic interpretation.
[0007] Cave-type and fracture-hole-type reservoirs often show beaded abnormal wave field response characteristics on seismic sections. Their amplitude energy is stronger than the seismic reflection energy of adjacent strata. At the same time, the amplitude energy of beads in conventional seismic data is usually overlapped with the seismic reflection energy of conventional strata. Small and medium-scale beaded anomalies with weak energy are easily ignored. Therefore, they cannot be effectively identified from seismic sections.
[0008] It should be noted that the information disclosed in the above background technology section is only used to enhance the understanding of the background of the present disclosure, and therefore may contain information that does not constitute prior art. Summary of the invention
[0009] In order to refine and implement the spatial and planar distribution characteristics of beaded anomalies, the present invention provides a carbonate reservoir beaded energy picking method to solve the problem that the beaded amplitude energy of conventional seismic data is usually overlapped with the conventional formation seismic reflection energy, resulting in small and medium-scale, weaker energy beaded anomalies being easily ignored and unable to be effectively identified on seismic profiles.
[0010] To achieve the above-mentioned purpose of the invention, the carbonate reservoir bead energy pickup method comprises: Input the post-stack seismic data volume of the reservoir; Preset a first processing window and a second processing window; Respectively constructing root mean square amplitude attribute data volumes of the post-stack seismic data volume under the first processing window condition and the second processing window condition and dividing the two to obtain a ratio data volume of two root mean square amplitude attributes; Setting a threshold value, and constructing a bead recognition ratio data body according to a comparison result between the ratio data body and the threshold value; The bead position and energy of the reservoir are picked out according to the bead identification ratio data body.
[0011] In the present disclosure and possible embodiments, the method for presetting the first processing window and the second processing window includes: The preset is performed by determining the length in the time direction or the number of seismic sampling points and the number of seismic traces or the distance occupied in the space direction according to the prominent characteristics of the bead amplitude energy of conventional seismic data on the seismic section; the second processing window is not equal to the first processing window.
[0012] In the present disclosure and possible embodiments, the number of seismic sampling points in the time direction is less than half of the total number of seismic sampling points in the time direction; The number of seismic channels in the spatial direction is less than half of the maximum number of seismic channels in the spatial direction.
[0013] In the present disclosure and possible embodiments, the method for constructing the root mean square amplitude attribute data volume under the conditions of the first processing window and the second processing window includes: Select a sampling point in the post-stack seismic data volume, place the sampling point at the time window center of the first processing window or the second processing window, input the value of the sampling point, and input the values of other sampling points around the sampling point that are included in the first processing window or the second processing window, set the values of the sampling points that are not included in the first processing window or the second processing window to zero, and calculate the root mean square amplitude of the sampling point through the sampling point value and the values of other sampling points included in the first processing window or the second processing window; Similarly, the root mean square amplitude of all sampling point positions in the post-stack seismic data volume is completed to obtain the root mean square amplitude attribute data volume of the post-stack seismic data volume under the first processing window or the second processing window.
[0014] In the present disclosure and possible embodiments, the root mean square amplitude of the sampling points is calculated as follows: The sum of the squares of the values of all sampling points in the first processing window or the second processing window is divided by the number of sampling points and then the square root is taken to obtain the root mean square amplitude of the sampling points under the conditions of the first processing window or the second processing window.
[0015] In the present disclosure and possible embodiments, the threshold value is set by: According to whether the second processing window is larger than or smaller than the first processing window, setting the threshold value to be larger than or smaller than 1; The degree of proximity of the threshold value to 1 corresponds to the size of the identifiable bead energy, and the closer the threshold value is to 1, the smaller the bead energy is identified.
[0016] In the present disclosure and possible embodiments, the method of constructing the bead recognition ratio data body is: Based on the threshold value being greater than 1, the value of the position less than the threshold value in the ratio data body is set as the threshold value, so that the ratio value of the two root mean square amplitude attributes at the position greater than the threshold value in the ratio data body remains unchanged; all spaces and sampling point positions of the reservoir profile are judged and set according to the same judgment method and setting method to obtain the bead recognition ratio data body; conversely, if the threshold value is less than 1, the method for judging whether it contains bead energy is completely opposite.
[0017] In the present disclosure and possible embodiments, the bead identification ratio data volume is multiplied with the post-stack seismic data volume of the reservoir to enhance the energy of the bead position.
[0018] In the present disclosure and possible embodiments, the values of the positions where the values in the ratio data volume are equal to the threshold value are set to 0, so as to obtain a seismic data volume containing only beaded energy; The seismic data volume containing only beading energy is superimposed with the post-stack seismic data volume to highlight the position of weak beading energy in the seismic section.
[0019] In the present disclosure and possible embodiments, the post-stack seismic data volume is a time migration data volume or a depth migration data volume; The post-stack seismic data volume is subjected to random noise suppression processing.
[0020] The picking method disclosed in the present invention constructs root mean square amplitude attribute data bodies under two different window size conditions to obtain a ratio data body of the two, and compares the ratio data body with a set threshold value to judge the bead energy of all spaces and sampling point positions of the profile one by one, thereby realizing automatic picking of bead energy, solving the problem that the bead amplitude energy of conventional seismic data and the seismic reflection energy of the formation are overlapped with each other and cannot be effectively identified, and laying a data foundation for the fine characterization of small-scale beads in carbonate reservoirs. BRIEF DESCRIPTION OF THE DRAWINGS
[0021] The above and other objects, features and advantages of the present disclosure will become more apparent through the following description of the embodiments of the present disclosure with reference to the accompanying drawings, in which: Figure 1 is a post-stack seismic data volume containing beads according to an embodiment of the present disclosure; Figure 2 is a ratio attribute volume diagram of an embodiment of the present disclosure; Figure 3 is a post-stack seismic data volume after bead energy enhancement according to an embodiment of the present disclosure; Figure 4 It is a seismic data volume with only bead energy remaining in the embodiment of the present disclosure; Figure 5 The picked bead energy according to the embodiment of the present disclosure is superimposed on the post-stack seismic data volume. Implementation
[0022] The present disclosure is described below based on embodiments, but it is worth noting that the present disclosure is not limited to these embodiments. In the detailed description of the present disclosure below, some specific details are described in detail. However, for the parts that are not described in detail, those skilled in the art can also fully understand the present disclosure.
[0023] In addition, those skilled in the art should understand that the drawings are only provided to illustrate the purpose, features and advantages of the present disclosure, and the drawings are not actually drawn to scale.
[0024] At the same time, unless the context clearly requires otherwise, the words "include", "comprises" and similar words throughout the specification and claims should be interpreted as inclusive rather than exclusive or exhaustive; that is, the meaning is "including but not limited to".
[0025] The carbonate reservoir bead energy picking method of the disclosed embodiment includes: step S10: inputting a post-stack seismic data body of the reservoir; step S20: presetting a first processing window and a second processing window; step S30: respectively constructing a root mean square amplitude attribute data body of the post-stack seismic data body under the first processing window condition and the second processing window condition and dividing the two to obtain a ratio data body of two root mean square amplitude attributes; step S40: setting a threshold value, and constructing a bead recognition ratio data body according to a comparison result between the ratio data body and the threshold value; and picking up the bead position and energy of the reservoir according to the bead recognition ratio data body.
[0026] Specifically, each step of the carbonate reservoir bead energy pickup method according to the embodiment of the present disclosure is described in detail as follows: Step S10: Input the post-stack seismic data volume of the reservoir: A post-stack seismic data volume of a carbonate reservoir is obtained, which may be a time-shifted data volume or a depth-shifted data volume. After the post-stack seismic data volume is input, random noise may be appropriately suppressed based on the resultant data volume according to conventional technical means in the art to improve the signal-to-noise ratio processing.
[0027] Step S20: Preset the first processing window and the second processing window. The specific presetting method of the embodiment of the present disclosure is as follows: The preset principle is: the length in the time direction or the number of seismic sampling points and the number of seismic traces or the distance occupied in the spatial direction are determined according to the prominent characteristics of the bead amplitude energy of conventional seismic data on the seismic profile for the preset; wherein, the second processing window is required not to be equal to the first processing window.
[0028] The preset means of the embodiment of the present disclosure is: the number of seismic sampling points in the time direction is less than half of the total number of seismic sampling points in the time direction; the number of seismic traces in the space direction is less than half of the maximum number of seismic traces in the space direction.
[0029] Specifically, the first processing window is preset, including the number of seismic sampling points in the vertical direction and the number of seismic traces in the spatial direction, which needs to be an odd number. Among them, N1t is the number of seismic sampling points in the time direction, which is usually less than half of the total number of seismic sampling points in the time direction; for the two-dimensional case, N1x is the number of seismic traces in the spatial direction, and for the three-dimensional case, N1x and N1y are the number of seismic traces in the two spatial directions of x and y, respectively, which are also less than half of the maximum number of traces in the corresponding directions.
[0030] The preset second processing window also includes the number of seismic sampling points in the vertical direction and the number of seismic traces in the spatial direction, which needs to be an odd number. Among them, N2t is the number of seismic sampling points in the time direction, which is usually less than half of the total number of seismic sampling points in the time direction; for the two-dimensional case, N2x is the number of seismic traces in the spatial direction, and for the three-dimensional case, N2x and N2y are the number of seismic traces in the two spatial directions of x and y, respectively, which are also less than half of the maximum number of traces in the corresponding directions.
[0031] Step S30: constructing root mean square amplitude attribute data volumes of the post-stack seismic data volume under the first processing window condition and the second processing window condition respectively and dividing the two to obtain a ratio data volume of two root mean square amplitude attributes. The specific process is as follows: A sampling point in the post-stack seismic data volume input in step S10 is selected, the position of the sampling point is placed at the time window center of the first processing window and the value of the sampling point is input, and at the same time, the values of other sampling points around the sampling point position and included in the first processing window are input, and the values of the sampling points not included in the first processing window are set to zero, and the root mean square amplitude of the sampling point is calculated through the value of the sampling point and the values of other sampling points included in the first processing window.
[0032] The calculation method of the root mean square amplitude of the sampling point is: The sum of the squares of the values of all sampling points in the first processing window is divided by the number of sampling points and then squared to obtain the root mean square amplitude of the sampling point under the condition of the first processing window.
[0033] The root mean square amplitudes of all sampling points in the post-stack seismic data volume input in step S10 are completed in the same manner to obtain the root mean square amplitude attribute data volume of the post-stack seismic data volume under the first processing window condition.
[0034] Similarly, the root mean square amplitude attribute data volume of the post-stack seismic data volume under the second processing window condition is constructed in the same way as the above construction method.
[0035] Step S40: setting a threshold value, constructing a bead recognition ratio data body according to the comparison result between the ratio data body and the threshold value; picking up the bead position and energy of the reservoir according to the bead recognition ratio data body, the specific process is as follows: The method for setting the threshold value in the embodiment of the present disclosure is as follows: based on the premise that the second processing window is larger or smaller than the first processing window, the threshold value is correspondingly set to be larger or smaller than 1; The degree of proximity of the threshold value to 1 corresponds to the size of the identifiable bead energy, and the closer the threshold value is to 1, the smaller the bead energy is identified.
[0036] If the threshold value is greater than 1, the value of the position less than the threshold value in the ratio data body is set as the threshold value, so that the ratio value of the two root mean square amplitude attributes of the position greater than the threshold value in the ratio data body remains unchanged; all spaces and sampling point positions of the reservoir profile are judged and set according to the same judgment method and setting method to obtain the bead recognition ratio data body; conversely, if the threshold value is less than 1, the method for judging whether it contains bead energy is completely opposite.
[0037] Finally, the bead recognition ratio data body is used to pick out the bead position and energy of the reservoir, realizing the automatic picking of bead energy, laying a solid data foundation for the fine characterization of small-scale beads in carbonate reservoirs.
[0038] Furthermore, based on the bead recognition ratio data body, the bead recognition ratio data body is multiplied with the post-stack seismic data body input in step S10, so that the energy of the bead position can be enhanced to achieve an enhanced picking effect.
[0039] In addition, by setting the value of the position where the value in the ratio data body is equal to the threshold value to 0, a seismic data body containing only beaded energy can be obtained, and then the seismic data body containing only beaded energy is superimposed with the post-stack seismic data body input in step S10, so that the position of weak beaded energy in the seismic profile can be highlighted, achieving the effect of highlighting the picking.
[0040] Specifically, the bead energy extraction method disclosed in the present invention is applied to a carbonate reservoir to perform bead energy extraction, and the specific process is as follows: 1. Input the post-stack seismic data volume of the reservoir Specifically, a line of post-stack seismic data is obtained, which is a time migration data volume. On the basis of this data volume, random noise is appropriately suppressed to improve the signal-to-noise ratio, such as Figure 1 As shown in the figure, the bead amplitude energy of conventional seismic data is overlapped with the formation seismic reflection energy, making it impossible to effectively identify the bead amplitude energy.
[0041] 2. Preset the first processing window The embodiment of the present disclosure utilizes conventional technical means in the field to preset the first processing window, specifically based on the prominent features of the beaded amplitude energy of conventional seismic data on the seismic profile, which refers to the number of seismic traces or distance in the spatial direction and the length or number of sampling points in the time direction, for example, at least 2 consecutive traces in the spatial direction and at least 5 sampling points in the time direction.
[0042] The preset first processing window includes the preset number of seismic sampling points in the vertical direction (i.e., the time direction) and the number of seismic traces in the spatial direction, which need to be an odd number. N1t is the number of seismic sampling points in the time direction, which is usually less than half of the total number of seismic sampling points in the time direction; for the two-dimensional case, N1x is the number of seismic traces in the spatial direction, and for the three-dimensional case, N1x and N1y are the number of seismic traces in the two spatial directions of x and y, respectively, which are also less than half of the maximum number of seismic traces in the corresponding directions.
[0043] The first processing window preset in this embodiment is that the number of vertical seismic sampling points N1t=11 and the number of spatial seismic traces N1x=11.
[0044] The preset first processing window also needs to preset the time window center, which refers to half of the time window (window) size. Taking the number of seismic channels in the spatial direction as an example, that is, 11 / 2=5.5, the time window center is the 6th channel, with 5 channels on the left and 5 channels on the right; the number of seismic channels in the spatial direction needs to be an odd number. If it is an even number, it needs to be changed to an odd number by adding 1 or subtracting 1 from the value, because there is no unique time window center point in the even number case. Similarly, the time window center in the time direction is the 6th sampling point number, so the time window center of this embodiment is at the 6th sampling point number and the 6th seismic channel number point position in the first processing window.
[0045] 3. Preset the second processing window The same method as the preset first processing window is adopted, where N2t is the number of seismic sampling points in the time direction, which is usually less than half of the total number of seismic sampling points in the time direction; for the two-dimensional case, N2x is the number of seismic traces in the spatial direction, and for the three-dimensional case, N2x and N2y are the number of seismic traces in the two spatial directions of x and y, respectively, which are also less than half of the maximum number of seismic traces in the corresponding directions.
[0046] In this embodiment, the number of vertical seismic sampling points N2t of the second processing window is 11, the number of seismic traces in the spatial direction is N2x=21, and the time window center of the second processing window is located at the 6th sampling point and the 11th seismic trace point in the window.
[0047] The second processing window preset in the embodiment of the present disclosure is larger than the size of the first processing window preset.
[0048] For the preset sizes of the first processing window and the second processing window, either the second processing window is larger than the first processing window, or the second processing window is smaller than the first processing window. Preferably, the second processing window is larger than the first processing window.
[0049] 4. Construct the root mean square amplitude attribute data volume of the post-stack seismic data volume under the first processing window condition Select a sampling point in the post-stack seismic data volume, place the sampling point at the center of the time window of the first processing window and obtain the value of the sampling point, and at the same time obtain the values of other sampling points around the sampling point that are included in the first processing window, set the values of the sampling points that are not included in the first processing window to zero, and calculate the root mean square amplitude of the sampling point through the value of the sampling point and the values of other sampling points included in the first processing window. Specifically, the sum of the squares of all the sampling point values in the first processing window is divided by the number of sampling points and then squared to obtain the root mean square amplitude of the sampling point under the conditions of the first processing window.
[0050] According to the above steps, the root mean square amplitude of all sampling points in the post-stack seismic data volume is completed, thereby obtaining the root mean square amplitude attribute data volume of the post-stack seismic data volume under the first processing window condition.
[0051] 5. Construct the root mean square amplitude attribute data volume of the post-stack seismic data volume under the second processing window condition The same method as constructing the root mean square amplitude attribute data volume within the first processing window can be used.
[0052] 6. Get the ratio data of two RMS amplitude attributes The RMS amplitude attribute data volume under the first processing window condition is divided by the RMS amplitude attribute data volume under the second processing window condition to obtain the ratio data volume of the two RMS amplitude attributes, such as Figure 2 shown.
[0053] 7. Set the threshold Because the second processing window preset in this embodiment is larger than the first processing window, a threshold value greater than 1 indicates that there are beads, and a threshold value equal to 1 indicates that there are no beads. In theory, any threshold value greater than 1 is acceptable. A threshold value that is too large can only identify large beads, and a smaller threshold value can identify more (small and medium-sized) beads; but the closer the threshold value is to 1, the less likely it is to reliably identify beads, because at this time, basically all positions are identified as beads. In this embodiment, based on experience, the best effect is selected when 1.5.
[0054] 8. Construct a bead recognition ratio data body according to the comparison result between the ratio data body and the threshold value, and pick up the bead position and bead energy of the reservoir according to the bead recognition data body.
[0055] Because the second processing window of this embodiment is larger than the first processing window, the root mean square amplitude attribute obtained under the second processing window conditions represents the macro background, while the root mean square amplitude attribute obtained under the first processing window conditions represents the local anomaly. The root mean square amplitude attribute body under the second processing window conditions is divided by the root mean square amplitude attribute body under the first processing window conditions to obtain a ratio data body. When the ratio data body is greater than a preset threshold value, it is considered that beading energy exists at the position. If it is less than the threshold value, it is considered that beading energy does not exist.
[0056] For this reason, since the positions smaller than the threshold value do not contain beading energy, the values of these positions are set as the threshold value; while the positions larger than the threshold value contain beading energy, the ratio values of the two root mean square amplitude attributes of these positions remain unchanged. According to the same judgment method and setting method, all spaces and sampling point positions of the reservoir profile are judged and set, and finally the beading recognition ratio data body is obtained.
[0057] Then, based on the bead recognition ratio data body, all bead positions in the reservoir are automatically identified, and the automatic picking of the bead positions is realized; further, the corresponding bead energy can be identified and extracted based on the set value size at the bead position.
[0058] 9. Enhanced after bead energy is picked up: The disclosed embodiment multiplies the bead identification ratio data volume with the post-stack seismic data volume of the reservoir input in step 1. Since the ratio coefficient multiplied by the bead position is larger than the coefficient multiplied by the non-bead position, the energy of the bead position will become stronger and the wave field characteristics will be more prominent, thus realizing the enhancement of the bead energy after picking up. Figure 3 shown.
[0059] 10. Construct a seismic data volume with only beading energy remaining to highlight the small and medium-scale beading characteristics and the location of weak beads in the seismic profile: The position values where the values in the ratio data volume of the embodiment of the present disclosure are equal to the preset threshold value are set to 0, and the post-stack seismic data volume containing only beading energy is obtained, such as Figure 4 As shown in FIG. 1 , the post-stack seismic data volume containing only the beaded energy is superimposed with the post-stack seismic data volume input in step 1, as shown in FIG. Figure 5 As shown, the beading characteristics, especially the small and medium-scale beading characteristics and the location of weak beads in the seismic section are highlighted, while the location and characteristics of large beads remain unchanged.
[0060] The above-described embodiments are only embodiments of the present disclosure, and the descriptions thereof are relatively specific and detailed, but they cannot be construed as limiting the scope of the present disclosure. It should be noted that, for a person of ordinary skill in the art, without departing from the concept of the present disclosure, several variations, equivalent substitutions, improvements, etc. may be made, all of which fall within the scope of protection of the present disclosure. Therefore, the scope of protection of the present disclosure shall be subject to the attached claims.
Claims
1. A carbonate reservoir bead energy pickup method, It is characterized in that include: Input the post-stack seismic data volume of the reservoir; Preset a first processing window and a second processing window; Respectively constructing root mean square amplitude attribute data volumes of the post-stack seismic data volume under the first processing window condition and the second processing window condition and dividing the two to obtain a ratio data volume of two root mean square amplitude attributes; Setting a threshold value, and constructing a bead recognition ratio data body according to a comparison result between the ratio data body and the threshold value; The bead position and energy of the reservoir are picked out according to the bead identification ratio data body.
2. The carbonate reservoir bead energy pickup method according to claim 1, It is characterized in that The method for presetting the first processing window and the second processing window includes: The preset is performed by determining the length in the time direction or the number of seismic sampling points and the number of seismic traces or the distance occupied in the space direction according to the prominent characteristics of the bead amplitude energy of conventional seismic data on the seismic section; the second processing window is not equal to the first processing window.
3. The carbonate reservoir bead energy pickup method according to claim 2, Features: The number of seismic sampling points in the time direction is less than half of the total number of seismic sampling points in the time direction; The number of seismic channels in the spatial direction is less than half of the maximum number of seismic channels in the spatial direction.
4. The carbonate reservoir bead energy pickup method according to claim 3, It is characterized in that The method for constructing the root mean square amplitude attribute data body under the conditions of the first processing window and the second processing window includes: Select a sampling point in the post-stack seismic data volume, place the sampling point at the time window center of the first processing window or the second processing window, input the value of the sampling point, and input the values of other sampling points around the sampling point that are included in the first processing window or the second processing window, set the values of the sampling points that are not included in the first processing window or the second processing window to zero, and calculate the root mean square amplitude of the sampling point through the sampling point value and the values of other sampling points included in the first processing window or the second processing window; Similarly, the root mean square amplitude of all sampling point positions in the post-stack seismic data volume is completed to obtain the root mean square amplitude attribute data volume of the post-stack seismic data volume under the first processing window or the second processing window.
5. The carbonate reservoir bead energy pickup method according to claim 4, It is characterized in that The root mean square amplitude of the sampling points is calculated as: The sum of the squares of the values of all sampling points in the first processing window or the second processing window is divided by the number of sampling points and then the square root is taken to obtain the root mean square amplitude of the sampling points under the conditions of the first processing window or the second processing window.
6. The carbonate reservoir bead energy pickup method according to any one of claims 1 to 5, It is characterized in that The threshold value is set as follows: Based on the premise that the second processing window is greater than or less than the first processing window, the threshold value is correspondingly set to be greater than or less than 1; The degree of proximity of the threshold value to 1 corresponds to the size of the identifiable bead energy, and the closer the threshold value is to 1, the smaller the bead energy is identified.
7. The carbonate reservoir bead energy pickup method according to claim 6, It is characterized in that The method for constructing the bead recognition ratio data body is: Based on the threshold value being greater than 1, the value of the position less than the threshold value in the ratio data body is set as the threshold value, so that the value of the ratio of the two root mean square amplitude attributes at the position greater than the threshold value in the ratio data body remains unchanged; all spaces and sampling point positions of the reservoir section are judged and set according to the same judgment method and setting method to obtain the bead recognition ratio data body; On the contrary, when the threshold value is less than 1, the method for determining whether beaded energy is contained is completely opposite.
8. The carbonate reservoir bead energy pickup method according to claim 7, Features: The bead identification ratio data volume is multiplied with the post-stack seismic data volume of the reservoir to enhance the energy of the bead locations.
9. The carbonate reservoir bead energy pickup method according to claim 7, Features: The values of the positions where the values in the ratio data body are equal to the threshold value are set to 0, so as to obtain a seismic data body containing only beaded energy; The seismic data volume containing only beading energy is superimposed with the post-stack seismic data volume to highlight the position of weak beading energy in the seismic section.
10. The carbonate reservoir bead energy pickup method according to any one of claims 1-5 or 7-9, Features: The post-stack seismic data volume is a time migration data volume or a depth migration data volume; The post-stack seismic data volume is subjected to random noise suppression processing.