A method and device for predicting a fractured reservoir by all-around attribute optimization

Abstract

The application provides a full-azimuth attribute optimization fracture reservoir prediction method and device, a storage medium and an electronic equipment, and relates to the oil and gas field exploration and development technical field.The method comprises the following steps: performing full-azimuth migration imaging processing on pre-stack CMP gathers of a work area to obtain azimuth pre-stack gathers; performing pre-stack fracture anisotropy inversion based on the azimuth pre-stack gathers to obtain a predicted fracture development azimuth; performing attribute calculation on the azimuth pre-stack gathers to obtain an attribute calculation result; performing optimization on the attribute calculation result by using the predicted fracture development azimuth to obtain an optimized attribute; and predicting a reservoir of the work area based on the optimized attribute.The technical scheme provided by the application can more simply and accurately predict the reservoir.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil and gas field exploration and development, in particular to a fracture reservoir prediction method and device based on all-around attribute optimization. BACKGROUND

[0002] The existing reservoir prediction is based on post-stack seismic data. The reservoir prediction based on post-stack seismic data is relatively cumbersome in method and has large amount of calculation, because post-stack seismic data needs to be obtained first. Moreover, the existing reservoir prediction is only based on the conventional attribute in post-stack seismic data, and the fracture information in the reservoir is not considered, so the prediction result is not accurate enough. SUMMARY

[0003] In view of the above problems in the prior art, the present application provides a fracture reservoir prediction method and device based on all-around attribute optimization, which can more simply and accurately predict the reservoir.

[0004] To achieve the above object, the technical scheme of the present application is as follows:

[0005] In a first aspect, the present application provides a fracture reservoir prediction method based on all-around attribute optimization, which comprises the following steps:

[0006] Performing all-around migration imaging processing on the pre-stack CMP gather of the work area to obtain a pre-stack azimuth gather;

[0007] Performing pre-stack fracture anisotropy inversion based on the pre-stack azimuth gather to obtain a predicted fracture development azimuth;

[0008] Performing attribute calculation on the pre-stack azimuth gather to obtain an attribute calculation result;

[0009] Optimizing the attribute calculation result by using the predicted fracture development azimuth to obtain an optimized attribute;

[0010] Performing reservoir prediction on the work area based on the optimized attribute.

[0011] Preferably, the all-around migration imaging processing on the pre-stack CMP gather of the work area to obtain a pre-stack azimuth gather comprises:

[0012] Performing all-around migration imaging processing on the pre-stack CMP gather of the work area by using a local angle domain imaging method or an all-around OVT processing method or a split-azimuth migration imaging processing method to obtain a pre-stack azimuth gather.

[0013] Preferably, the attribute calculation on the pre-stack azimuth gather to obtain an attribute calculation result comprises:

[0014] perform single-attribute calculation on the azimuth pre-stack gathers to obtain azimuth pre-stack attribute gathers;

[0015] perform azimuth angle stacking on the azimuth pre-stack gathers to obtain azimuth stacking data volume;

[0016] perform three-dimensional attribute calculation on the azimuth stacking data volume to obtain three-dimensional azimuth stacking attribute volume;

[0017] obtain the azimuth pre-stack attribute gathers and the three-dimensional azimuth stacking attribute volume as the attribute calculation result.

[0018] Preferably, the single-attribute calculation on the azimuth pre-stack gathers to obtain azimuth pre-stack attribute gathers comprises:

[0019] perform attribute calculation of instantaneous amplitude, instantaneous phase, instantaneous frequency and root mean square amplitude on each seismic trace of the azimuth pre-stack gathers to obtain the azimuth pre-stack attribute gathers.

[0020] Preferably, the azimuth angle stacking on the azimuth pre-stack gathers to obtain azimuth stacking data volume comprises:

[0021] for the pre-stack gathers in each azimuth of the azimuth pre-stack gathers, perform the following operations:

[0022] based on a preset azimuth angle and a preset total offset, stack the pre-stack gathers in the azimuth to obtain a stacking data volume in the azimuth; wherein the preset azimuth angle is preset based on the number of the pre-stack gathers in the azimuth;

[0023] obtain the stacking data volume in each azimuth as the azimuth stacking data volume.

[0024] Preferably, the three-dimensional attribute calculation on the azimuth stacking data volume to obtain three-dimensional azimuth stacking attribute volume comprises:

[0025] perform calculation of coherent attribute, curvature attribute and ant attribute on the azimuth stacking data volume to obtain the three-dimensional azimuth stacking attribute volume.

[0026] Preferably, the adoption of the predicted fracture development azimuth to the attribute calculation result to obtain preferred attribute comprises:

[0027] adopt the predicted fracture development azimuth to the azimuth pre-stack attribute gathers to obtain preferred azimuth pre-stack attribute gathers;

[0028] adopt the predicted fracture development azimuth to the three-dimensional azimuth stacking attribute volume to obtain preferred three-dimensional azimuth stacking attribute volume;

[0029] The preferred azimuth pre-stack attribute gather and the preferred three-dimensional azimuth stacking attribute volume are obtained as the preferred attribute.

[0030] Preferably, the preferred azimuth pre-stack attribute gather is obtained by preferentially selecting the azimuth pre-stack attribute gather according to the predicted fracture development direction.

[0031] For each sampling point in each seismic trace of the azimuth pre-stack attribute gather, the following operations are performed:

[0032] An attribute gather of the sampling point in the seismic trace is obtained by selecting an attribute trace within a preset angle range from the azimuth pre-stack attribute gather, wherein the angle value of the predicted fracture development direction is taken as a center point of the preset angle range, and the preset angle range is determined based on the density of the sampling points.

[0033] The attribute gather of each sampling point in each seismic trace is obtained as the preferred azimuth pre-stack attribute gather.

[0034] Preferably, the preferred three-dimensional azimuth stacking attribute volume is obtained by preferentially selecting the three-dimensional azimuth stacking attribute volume according to the predicted fracture development direction.

[0035] For each sampling point in each seismic trace of the three-dimensional azimuth stacking attribute volume, the following operations are performed:

[0036] An attribute volume corresponding to an angle with the smallest absolute value of the difference value of the angle value of the predicted fracture development direction is obtained by selecting the attribute volume from the three-dimensional azimuth stacking attribute volume, wherein the attribute volume is corresponding to the angle of the sampling point in the seismic trace.

[0037] The attribute volume of each sampling point in each seismic trace is obtained as the preferred three-dimensional azimuth stacking attribute volume.

[0038] In a second aspect, an embodiment of the present application provides a full-azimuth attribute optimization fracture reservoir prediction device, which comprises:

[0039] The migration imaging processing module is configured to perform full-azimuth migration imaging processing on the pre-stack CMP gather of the work area to obtain an azimuth pre-stack gather.

[0040] The inversion module is configured to perform pre-stack fracture anisotropy inversion based on the azimuth pre-stack gather to obtain a predicted fracture development direction.

[0041] The attribute calculation module is configured to perform attribute calculation on the azimuth pre-stack gather to obtain an attribute calculation result.

[0042] preferably module, for preferentially selecting the attribute calculation result by using the predicted fracture development direction, to obtain a preferred attribute;

[0043] prediction module, for predicting the reservoir of the work area based on the preferred attribute.

[0044] In a third aspect, an embodiment of the present application provides a storage medium, which has program code stored thereon, and the program code, when executed by a processor, implements the fracture reservoir prediction method with omnibearing attribute optimization.

[0045] In a fourth aspect, an embodiment of the present application provides an electronic device, which comprises a memory and a processor, and the memory has program code stored thereon and executable on the processor, and the program code, when executed by the processor, implements the fracture reservoir prediction method with omnibearing attribute optimization.

[0046] The fracture reservoir prediction method with omnibearing attribute optimization, the device, the storage medium and the electronic device provided by the embodiments of the present application can obtain an azimuth prestack gather by performing omnibearing migration imaging processing on a prestack CMP gather of a work area, perform prestack fracture anisotropy inversion based on the azimuth prestack gather, obtain a predicted fracture development direction, perform attribute calculation on the azimuth prestack gather, obtain an attribute calculation result, preferentially select the attribute calculation result by using the predicted fracture development direction, obtain a preferred attribute, and predict the reservoir of the work area based on the preferred attribute, so that the prediction of the reservoir can be performed based on prestack seismic data, rather than poststack seismic data as in the prior art, and the prediction method is simplified. Meanwhile, the present application performs prestack fracture anisotropy inversion based on the azimuth prestack gather and preferentially selects the attribute calculation result by using the predicted fracture development direction, that is, the present application fully considers the fracture development information in the reservoir, and can obtain a more accurate prediction result compared with the prior art. It can be seen that the technical solution provided by the embodiments of the present application can more simply and accurately predict the reservoir. BRIEF DESCRIPTION OF DRAWINGS

[0047] The scope of the present application can be better understood by reading the following detailed description of exemplary embodiments in conjunction with the attached drawings, in which:

[0048] Figure 1 Method flow of the embodiments of the present application Figure One ;

[0049] Figure 2 Method flow of the embodiments of the present application Figure Two ;

[0050] Figure 3A schematic diagram of azimuth pre-stack gather in the embodiment of the present application;

[0051] Figure 4 A schematic diagram of fracture anisotropy strength in the embodiment of the present application;

[0052] Figure 5 A profile diagram of predicted fracture development result in the embodiment of the present application;

[0053] Figure 6 A schematic diagram of azimuth pre-stack attribute gather (instantaneous phase) in the embodiment of the present application;

[0054] Figure 7 A schematic diagram of azimuth stack data volume in the embodiment of the present application;

[0055] Figure 8 A schematic diagram of three-dimensional azimuth stack attribute volume in the embodiment of the present application;

[0056] Figure 9 A schematic diagram of preferred azimuth pre-stack attribute gather in the embodiment of the present application;

[0057] Figure 10 A result diagram of a certain fracture-type reservoir development area in a certain oilfield predicted by using the method described in the embodiment of the present application.

[0058] Figure 11 A device structure diagram of the embodiment of the present application. DETAILED DESCRIPTION

[0059] In order to make the purpose, technical scheme and advantages of the present application more clear, the implementation method of the present application will be described in detail below in combination with the drawings and embodiments, so that the implementation process of how to apply technical means to solve the technical problem and achieve the technical effect can be fully understood and implemented.

[0060] In the following description, many specific details are set forth in order to provide a thorough understanding of the present application, however, the present application can also be implemented in other ways different from those described herein, therefore, the protection scope of the present application is not limited by the specific embodiments disclosed below.

[0061] Example One

[0062] For conventional three-dimensional post-stack seismic data volume, using seismic attribute to carry out reservoir prediction is a direct and effective method, however, how to carry out attribute calculation on pre-stack gather, especially full-azimuth pre-stack gather, and make the attribute be able to be directly and effectively combined with geological target, at present, there is still no effective means and method.

[0063] The application faces all-around pre-stack gathers, uses all-around AVAZ inversion to obtain the strength and direction of azimuth anisotropy, carries out attribute calculation on the all-around gathers, carries out all-around attribute optimization under the direction constraint, and obtains the all-around optimized seismic attribute calculation result.

[0064] According to the embodiment of the application, a fracture reservoir prediction method for all-around attribute optimization is provided, as shown in the figure. Figure 1 The method comprises the steps of:

[0065] In step S101, all-around migration imaging processing is performed on the pre-stack CMP gathers of the working area to obtain azimuth pre-stack gathers.

[0066] In this embodiment, the all-around migration imaging processing on the pre-stack CMP (Common Middle Point) gathers of the working area to obtain azimuth pre-stack gathers comprises:

[0067] The local angle domain imaging method, the all-around OVT processing method or the split-azimuth migration imaging processing method is used to perform all-around migration imaging processing on the pre-stack CMP gathers of the working area to obtain azimuth pre-stack gathers.

[0068] Specifically, the local angle domain imaging method, the all-around OVT processing method or the split-azimuth migration imaging processing method is used to perform all-around migration imaging processing on the pre-stack CMP gathers of the working area to obtain split-azimuth split-offset gathers or split-azimuth split-angle gathers, i.e., to obtain the above-mentioned azimuth pre-stack gathers.

[0069] In step S102, pre-stack fracture anisotropy inversion is performed based on the azimuth pre-stack gathers to obtain a predicted fracture development direction.

[0070] In this embodiment, the velocity, the time difference or the amplitude of the seismic wave can be used to calculate the anisotropy strength and development direction of the fracture, and the calculated development direction of the fracture is the predicted fracture development direction.

[0071] In step S103, attribute calculation is performed on the azimuth pre-stack gathers to obtain an attribute calculation result.

[0072] Attribute calculation is a conventional method for reservoir prediction using seismic attributes, and can be divided into single-channel attribute calculation and three-dimensional attribute calculation according to different algorithms. The single-channel attribute calculation is an attribute algorithm that can be carried out on a single seismic channel, such as three-instant attributes (instantaneous amplitude, instantaneous phase, instantaneous frequency, etc.) and some other simple attributes. Some more complex attributes (such as coherence attribute, curvature attribute, ant attribute, etc.) need to be calculated on a three-dimensional data body. The embodiment can simultaneously perform all-around optimization for single-channel attribute calculation and three-dimensional attribute calculation.

[0073] In the embodiment, the attribute calculation on the azimuth pre-stack gather is performed to obtain an attribute calculation result, including:

[0074] The single-attribute calculation is performed on the azimuth pre-stack gather to obtain an azimuth pre-stack attribute gather;

[0075] The azimuth angle stacking is performed on the azimuth pre-stack gather to obtain an azimuth stacking data volume;

[0076] The three-dimensional attribute calculation is performed on the azimuth stacking data volume to obtain a three-dimensional azimuth stacking attribute volume;

[0077] The azimuth pre-stack attribute gather and the three-dimensional azimuth stacking attribute volume are obtained as the attribute calculation result.

[0078] In the embodiment, the single-attribute calculation is performed on the azimuth pre-stack gather to obtain an azimuth pre-stack attribute gather, including:

[0079] The attribute calculation of the instantaneous amplitude, the instantaneous phase, the instantaneous frequency and the root mean square amplitude is performed on each seismic trace of the azimuth pre-stack gather to obtain the azimuth pre-stack attribute gather.

[0080] For a more complex three-dimensional seismic attribute, the azimuth angle stacking needs to be performed on the azimuth pre-stack gather to obtain an azimuth stacking data volume, and then the three-dimensional attribute calculation is performed on the azimuth stacking data volume to obtain a three-dimensional azimuth stacking attribute volume.

[0081] In the embodiment, the azimuth angle stacking is performed on the azimuth pre-stack gather to obtain an azimuth stacking data volume, including:

[0082] For the pre-stack gather in each azimuth of the azimuth pre-stack gather, the following operation is performed: based on a preset azimuth angle and a preset total offset, the pre-stack gather in the azimuth is stacked to obtain a stacking data volume in the azimuth; wherein the preset azimuth angle is preset based on the number of the pre-stack gather in the azimuth;

[0083] The stacking data volume in each azimuth is obtained as the azimuth stacking data volume.

[0084] Specifically, for the pre-stack gathers of each azimuth in the azimuth pre-stack gather, an azimuth angle and a full offset are preset. The azimuth angle is determined based on the number of the pre-stack gathers of the azimuth, for example, the azimuth angle is set to 30° when the number of the pre-stack gathers is 6, and the azimuth angle is set to 45° when the number of the pre-stack gathers is 4. The value of the full offset ranges from 0 to the maximum value of the full offset. Based on the azimuth angle and the full offset, the pre-stack gathers of the azimuth are superimposed, that is, the pre-stack gathers of the azimuth are all added together, and then divided by the number of the pre-stack gathers of the azimuth, so that the superposition data volume of the azimuth can be obtained. The superposition data volumes of different azimuths, that is, the above-mentioned azimuth superposition data volumes, can be obtained by obtaining the superposition data volume of each azimuth. The data line number, the trace number, the number of seismic traces, and the sampling number of each seismic trace of the superposition data volumes of different azimuths are all the same, and only the seismic data of some azimuths may be different due to the superposition of different azimuths.

[0085] In the embodiment, the azimuth superposition data volume is subjected to three-dimensional attribute calculation to obtain a three-dimensional azimuth superposition attribute volume, including:

[0086] The azimuth superposition data volume is subjected to calculation of coherent attribute, curvature attribute and ant attribute to obtain the three-dimensional azimuth superposition attribute volume.

[0087] In step S104, the predicted fracture development azimuth is used to optimize the attribute calculation result to obtain an optimized attribute.

[0088] In the embodiment, the predicted fracture development azimuth is used to optimize the attribute calculation result to obtain an optimized attribute, including:

[0089] The predicted fracture development azimuth is used to optimize the azimuth pre-stack attribute gather to obtain an optimized azimuth pre-stack attribute gather.

[0090] The predicted fracture development azimuth is used to optimize the three-dimensional azimuth superposition attribute volume to obtain an optimized three-dimensional azimuth superposition attribute volume.

[0091] The optimized azimuth pre-stack attribute gather and the optimized three-dimensional azimuth superposition attribute volume are obtained as the optimized attribute.

[0092] The predicted fracture development azimuth is used to optimize the azimuth pre-stack attribute gather to obtain an optimized azimuth pre-stack attribute gather, including:

[0093] For each sampling point in each seismic trace of the azimuth pre-stack attribute gather, the following operation is performed: attribute traces within a preset angle range are selected from the azimuth pre-stack attribute gather to obtain an attribute gather of the sampling point in the seismic trace; wherein the preset angle range takes the angle value of the predicted fracture development azimuth as a center point; and the preset angle range is determined based on the density of the sampling points;

[0094] The attribute gather of each sampling point in each seismic trace is obtained as the preferred azimuth pre-stack attribute gather.

[0095] Specifically, since each seismic trace corresponds to a fracture azimuth, and the specific value of the azimuth is between 0° and 180°, for each sampling point in each seismic trace of the azimuth pre-stack attribute gather, there are M offset (incident angle) and N azimuth angle sampling points. According to the predicted fracture development azimuth obtained in step S102, attribute traces within a preset angle range are selected from the azimuth pre-stack attribute gather to obtain an attribute gather of the sampling point in the seismic trace. The preset angle range is generally not more than 30°. Of course, in actual application, the preset angle range can be appropriately reduced according to the number of gathers, but generally the minimum is not less than 5°.

[0096] For example, if the angle value of the predicted fracture development azimuth is 78°, the preset angle range can be set to 63°-93°. Of course, if the sampling points are sufficiently dense, the traces between 63° and 93° can also be selected between 73° and 83°.

[0097] By analogy, the attribute gather of each sampling point in each seismic trace is obtained as the preferred azimuth pre-stack attribute gather in the above manner. The preferred azimuth pre-stack attribute gather carries not only the reservoir properties represented by conventional attributes (such as instantaneous amplitude representing energy changes, instantaneous phase reflecting underground anomalies, and instantaneous frequency reflecting lithology), but also the direction information of the fractures, which is very advantageous for seismic prediction in fracture reservoirs.

[0098] In this embodiment, the predicted fracture development azimuth is used to optimize the three-dimensional azimuth stacking attribute volume to obtain a preferred three-dimensional azimuth stacking attribute volume, including:

[0099] For each sampling point in each seismic trace of the three-dimensional azimuth stacking attribute volume, the following operation is performed: an attribute volume corresponding to an angle with the smallest absolute value difference from the angle value of the predicted fracture development azimuth is selected from the three-dimensional azimuth stacking attribute volume to obtain an attribute volume of the sampling point in the seismic trace;

[0100] acquire the attribute body at each sampling point in each seismic trace as the preferred three-dimensional azimuthal prestack attribute body.

[0101] Specifically, since each seismic trace corresponds to a fracture azimuth, and the specific value of the azimuth is between 0° and 180°, for each sampling point in each seismic trace of the three-dimensional azimuthal prestack attribute body, there are N azimuthal angle (without offset and incidence angle information) sampling points. According to the predicted fracture development azimuth obtained in step S102, the attribute body corresponding to the angle with the smallest absolute value of the difference between the angle value of the predicted fracture development azimuth is selected from the three-dimensional azimuthal prestack attribute body, that is, the azimuth attribute body seismic data closest to the fracture azimuth is directly preferred, and the attribute body at the sampling point in the seismic trace is obtained.

[0102] For example, the angle value of the predicted fracture development azimuth is 78°, and the angle of the attribute body closest to the angle value is 90°, so the attribute body corresponding to 90° is selected as the attribute body at the sampling point in the seismic trace.

[0103] By analogy, the attribute body at each sampling point in each seismic trace is acquired as the preferred three-dimensional azimuthal prestack attribute body in the above manner. The preferred three-dimensional azimuthal prestack attribute body carries not only the reservoir properties represented by conventional attributes (for example, the discontinuity of a coherent event represents the development of a small fault, the curvature represents the bending degree of a coherent event, and the ant body more intuitively represents the distribution of faults / fractures), but also the direction information of the fractures, which is very advantageous for seismic prediction in a fracture reservoir.

[0104] In step S105, the reservoir of the work area is predicted based on the preferred attribute.

[0105] In this embodiment, predicting the reservoir of the work area based on the preferred attribute comprises:

[0106] Based on the preferred azimuthal prestack attribute gather and the preferred three-dimensional azimuthal prestack attribute body, the reservoir of the work area is jointly predicted to predict the overall development of the reservoir.

[0107] The fracture reservoir prediction method based on omnidirectional attribute optimization provided in this invention involves performing omnidirectional migration imaging on pre-stack CMP gathers in the work area to obtain azimuth pre-stack gathers. Based on these azimuth pre-stack gathers, pre-stack fracture anisotropy inversion is performed to obtain the predicted fracture development azimuth. Attribute calculations are then performed on these azimuth pre-stack gathers to obtain attribute calculation results. The predicted fracture development azimuth is then used to optimize the attribute calculation results to obtain preferred attributes. Based on these preferred attributes, reservoir prediction in the work area is performed. This allows reservoir prediction to be based on pre-stack seismic data, eliminating the need to calculate post-stack seismic data as in existing technologies, thus simplifying the prediction method. Furthermore, because this invention performs pre-stack fracture anisotropy inversion based on azimuth pre-stack gathers and optimizes the attribute calculation results using the predicted fracture development azimuth, it fully considers fracture development information in the reservoir, resulting in more accurate prediction results compared to existing technologies. Therefore, the technical solution provided in this invention can predict reservoirs more simply and accurately.

[0108] Example Two

[0109] This embodiment uses actual reservoir prediction as an example to further explain in detail the fracture reservoir prediction method with comprehensive attribute optimization.

[0110] like Figure 2 As shown, the method described in this embodiment includes:

[0111] Step S201: Perform omnidirectional migration imaging processing on the pre-stack CMP gathers of the work area to obtain the azimuth pre-stack gathers.

[0112] Specifically, the pre-stack CMP gathers in the work area are processed using either local angle domain imaging, omnidirectional OVT processing, or azimuth-based migration imaging, resulting in azimuth-based and offset-based gathers or azimuth-based and angle-based gathers, thus obtaining the aforementioned azimuth pre-stack gathers. A schematic diagram of this azimuth pre-stack gather is shown below. Figure 3 As shown.

[0113] Step S202: Perform pre-stack fracture anisotropy inversion based on the azimuth pre-stack gather to obtain the predicted fracture development azimuth.

[0114] In this embodiment, the anisotropy intensity and development direction of the crack can be calculated by using the velocity, time difference or amplitude of the seismic wave, that is, the predicted crack development result is obtained, and the development direction of the crack calculated above is the predicted crack development direction.

[0115] Figure 5 A cross-sectional view of the predicted crack development results is shown above. A schematic diagram of the crack anisotropy intensity is shown below. Figure 4 As shown.

[0116] Step S203, single-channel attribute calculation is performed on the azimuth pre-stack gathers to obtain azimuth pre-stack attribute gathers;

[0117] Specifically, attribute calculation of instantaneous amplitude, instantaneous phase, instantaneous frequency and root mean square amplitude is performed on each seismic trace of the azimuth pre-stack gathers to obtain the azimuth pre-stack attribute gathers.

[0118] The azimuth pre-stack attribute gathers obtained by performing attribute calculation of instantaneous phase are as shown in Figure 6 .

[0119] Step S204, azimuth angle stacking is performed on the azimuth pre-stack gathers to obtain azimuth stacking data volumes;

[0120] In the embodiment, the azimuth angle stacking is performed on the azimuth pre-stack gathers to obtain azimuth stacking data volumes, which includes:

[0121] For the pre-stack gathers on each azimuth of the azimuth pre-stack gathers, the following operation is performed: based on a preset azimuth angle and a preset total offset, the pre-stack gathers on the azimuth are stacked to obtain a stacking data volume on the azimuth; wherein the preset azimuth angle is preset based on the number of the pre-stack gathers on the azimuth;

[0122] The stacking data volume on each azimuth is obtained as the azimuth stacking data volume.

[0123] Specifically, for the pre-stack gathers on each azimuth of the azimuth pre-stack gathers, an azimuth angle and a total offset are preset. Wherein the azimuth angle is determined based on the number of the pre-stack gathers on the azimuth, for example, if the number of the pre-stack gathers is 6, the azimuth angle is set to 30°; if the number of the pre-stack gathers is 4, the azimuth angle is set to 45°. The numerical range of the total offset is 0 to the maximum value of the full offset. Based on the above-mentioned azimuth angle and total offset, the pre-stack gathers on the azimuth are stacked, that is, the pre-stack gathers on the azimuth are all added, and then divided by the number of the pre-stack gathers on the azimuth, to obtain the stacking data volume on the azimuth. The stacking data volume on each azimuth can obtain the seismic stacking data volumes of different azimuths, that is, the above-mentioned azimuth stacking data volumes. Wherein the data line number, trace number, seismic trace number and sampling number of each seismic trace of the stacking data volumes of different azimuths are all the same, only the seismic data of some azimuths may be different due to the stacking of different azimuths.

[0124] The schematic diagram of the azimuth stacking data volume obtained by the above-mentioned method is as shown in Figure 7 .

[0125] Step S205, performing three-dimensional attribute calculation on the azimuth stack data volume to obtain a three-dimensional azimuth stack attribute volume;

[0126] In this embodiment, the three-dimensional attribute calculation on the azimuth stack data volume to obtain a three-dimensional azimuth stack attribute volume includes:

[0127] The three-dimensional attribute calculation on the azimuth stack data volume includes calculation of coherent attribute, curvature attribute and ant attribute to obtain the three-dimensional azimuth stack attribute volume. A schematic diagram of the three-dimensional azimuth stack attribute volume is shown in FIG. 3. Figure 8

[0128] Step S206, performing optimization on the azimuth pre-stack attribute gather using the predicted fracture development azimuth to obtain an optimized azimuth pre-stack attribute gather.

[0129] In this embodiment, the optimization on the azimuth pre-stack attribute gather using the predicted fracture development azimuth to obtain an optimized azimuth pre-stack attribute gather includes:

[0130] For each sampling point in each seismic trace of the azimuth pre-stack attribute gather, the following operation is performed: selecting attribute traces within a preset angle range from the azimuth pre-stack attribute gather to obtain an attribute gather at the sampling point in the seismic trace; wherein the preset angle range has the angle value of the predicted fracture development azimuth as a center point; and the preset angle range is determined based on the density of the sampling points.

[0131] The attribute gather at each sampling point in each seismic trace is obtained as the optimized azimuth pre-stack attribute gather.

[0132] Specifically, since each seismic trace corresponds to an azimuth of a fracture, and the specific value of the azimuth is between 0° and 180°, for each sampling point in each seismic trace of the azimuth pre-stack attribute gather, there are M offset distances (incident angles) and N azimuth angles of sampling points. According to the predicted fracture development azimuth obtained in step S102, attribute traces within a preset angle range are selected from the azimuth pre-stack attribute gather to obtain an attribute gather at the sampling point in the seismic trace. The preset angle range is generally not more than 30°. Of course, in actual application, the preset angle range can be appropriately reduced according to the number of gathers, but generally the minimum is not less than 5°.

[0133] For example, if the angle value of the predicted fracture development azimuth is 78°, the preset angle range can be set to 63°-93°. Of course, if the sampling points are dense enough, the traces between 63° and 93° can also be selected between 73° and 83°.

[0134] ​Similarly, the attribute gathers at each sampling point in each seismic trace are obtained in the above manner as the preferred azimuth pre-stack attribute gathers. The preferred azimuth pre-stack attribute gathers not only carry the reservoir properties represented by conventional attributes (such as the instantaneous amplitude representing the energy change, the instantaneous phase representing the underground anomaly, the instantaneous frequency representing the lithology, etc.), but also carry the direction information of the fractures, which is very advantageous for seismic prediction in the fractured reservoir.

[0135] A schematic diagram of the preferred azimuth pre-stack attribute gathers is shown in Figure 9

[0136] Step S207, the three-dimensional azimuth stacking attribute volume is preferred using the predicted fracture development azimuth to obtain a preferred three-dimensional azimuth stacking attribute volume.

[0137] In this embodiment, the three-dimensional azimuth stacking attribute volume is preferred using the predicted fracture development azimuth to obtain a preferred three-dimensional azimuth stacking attribute volume, including:

[0138] For each sampling point in each seismic trace of the three-dimensional azimuth stacking attribute volume, the following operation is performed: from the three-dimensional azimuth stacking attribute volume, an attribute volume corresponding to an angle with the smallest absolute value of the difference between the angle value of the predicted fracture development azimuth is selected to obtain the attribute volume at the sampling point in the seismic trace.

[0139] The attribute volume at each sampling point in each seismic trace is obtained as the preferred three-dimensional azimuth stacking attribute volume.

[0140] Specifically, since each seismic trace corresponds to a fracture azimuth, and the specific value of the azimuth is between 0° and 180°, for each sampling point in each seismic trace of the three-dimensional azimuth stacking attribute volume, there are N azimuth angle (without offset and incidence angle information) sampling points. According to the predicted fracture development azimuth obtained in step S102, an attribute volume corresponding to an angle with the smallest absolute value of the difference between the angle value of the predicted fracture development azimuth is selected from the three-dimensional azimuth stacking attribute volume, that is, the azimuth attribute volume seismic data closest to the fracture azimuth is directly preferred to obtain the attribute volume at the sampling point in the seismic trace.

[0141] For example, the angle value of the predicted fracture development azimuth is 78°, and the angle of the attribute volume closest to the angle value is 90°, so the attribute volume corresponding to 90° is selected as the attribute volume at the sampling point in the seismic trace.

[0142] ​Similarly, the attribute volume at each sampling point in each seismic trace is obtained in the above manner as the preferred three-dimensional azimuthal stack attribute volume. The preferred three-dimensional azimuthal stack attribute volume carries not only the reservoir property represented by the conventional attribute (for example, the fault development represented by the discontinuity of the event represented by the coherence, the tortuosity represented by the curvature, and the fracture represented by the ant tracking), but also the direction information of the fracture, which is very advantageous for the seismic prediction in the fractured reservoir.

[0143] In step S208, the reservoir in the working area is predicted based on the preferred azimuthal pre-stack attribute gather and the preferred three-dimensional azimuthal stack attribute volume.

[0144] Specifically, the development of the whole reservoir can be predicted based on the preferred azimuthal pre-stack attribute gather and the preferred three-dimensional azimuthal stack attribute volume.

[0145] Figure 10 The result of the prediction of the development of a fractured reservoir in an oilfield is shown in the figure, from which it can be seen that the development of the fractured reservoir is clearly predicted, and the coincidence rate is good after the comparison and analysis with the FMI imaging in the well.

[0146] The conventional seismic attribute method is usually based on the stacked seismic data, and the technical solution provided in the present application is calculated on the seismic pre-stack gather, the attribute calculation involved makes full use of the offset and azimuth information obtained by the current wide-azimuth seismic acquisition gather, and compared with the coherence volume obtained by the stacked seismic data, the method has more abundant fracture development information, and through the optimization of the azimuthal pre-stack attribute gather and the azimuthal stack attribute, the development and distribution map of the fracture or the geological anomaly volume in the working area can be obtained.

[0147] The fracture reservoir prediction method provided by the embodiment of the present application can obtain the predicted fracture development direction based on the azimuth prestack gather, perform prestack fracture anisotropy inversion based on the azimuth prestack gather, obtain the predicted fracture development direction, perform attribute calculation on the azimuth prestack gather, obtain the attribute calculation result, perform optimization on the attribute calculation result by using the predicted fracture development direction, obtain the optimized attribute, and predict the reservoir of the work area based on the optimized attribute, so that the prediction of the reservoir can be performed based on the prestack seismic data, and it is not necessary to calculate the poststack seismic data as in the prior art, and the prediction method is simplified. Meanwhile, the present application performs prestack fracture anisotropy inversion based on the azimuth prestack gather, and performs optimization on the attribute calculation result by using the predicted fracture development direction, that is, the present application fully considers the fracture development information in the reservoir, and can obtain more accurate prediction result compared with the prior art. It can be seen that the technical solution provided by the embodiment of the present application can more simply and accurately predict the reservoir.

[0148] Example Three

[0149] Corresponding to the method embodiment, the present application further provides a fracture reservoir prediction device for azimuth attribute optimization, as shown in Figure 11 The device comprises:

[0150] The migration imaging processing module 301 is configured to perform azimuth migration imaging processing on the prestack CMP gather of the work area to obtain the azimuth prestack gather.

[0151] The inversion module 302 is configured to perform prestack fracture anisotropy inversion based on the azimuth prestack gather to obtain the predicted fracture development direction.

[0152] The attribute calculation module 303 is configured to perform attribute calculation on the azimuth prestack gather to obtain the attribute calculation result.

[0153] The optimization module 304 is configured to perform optimization on the attribute calculation result by using the predicted fracture development direction to obtain the optimized attribute.

[0154] The prediction module 305 is configured to predict the reservoir of the work area based on the optimized attribute.

[0155] In the embodiment, the migration imaging processing module 301 obtains the azimuth prestack gather in the following manner:

[0156] The azimuth migration imaging processing is performed on the prestack CMP gather of the work area by using the local angle domain imaging method or the azimuth OVT processing method or the azimuth-dependent migration imaging processing method to obtain the azimuth prestack gather.

[0157] In this embodiment, the attribute calculation module 303 comprises:

[0158] a single-channel attribute calculation unit, configured to perform single-channel attribute calculation on the azimuth pre-stack gather to obtain an azimuth pre-stack attribute gather;

[0159] an azimuth angle stacking unit, configured to perform azimuth angle stacking on the azimuth pre-stack gather to obtain an azimuth stacking data volume;

[0160] a three-dimensional attribute calculation unit, configured to perform three-dimensional attribute calculation on the azimuth stacking data volume to obtain a three-dimensional azimuth stacking attribute volume;

[0161] an attribute calculation result acquisition unit, configured to acquire the azimuth pre-stack attribute gather and the three-dimensional azimuth stacking attribute volume as the attribute calculation result.

[0162] In this embodiment, the single-channel attribute calculation unit obtains the azimuth pre-stack attribute gather in the following manner:

[0163] perform attribute calculation of instantaneous amplitude, instantaneous phase, instantaneous frequency and root mean square amplitude on each seismic trace of the azimuth pre-stack gather to obtain the azimuth pre-stack attribute gather.

[0164] In this embodiment, the azimuth angle stacking unit obtains the azimuth stacking data volume in the following manner:

[0165] for the pre-stack gather at each azimuth in the azimuth pre-stack gather, the following operations are performed:

[0166] based on a preset azimuth angle and a preset total offset, the pre-stack gather at the azimuth is stacked to obtain a stacking data volume at the azimuth; wherein the preset azimuth angle is preset based on the number of the pre-stack gather at the azimuth;

[0167] the stacking data volume at each azimuth is acquired as the azimuth stacking data volume.

[0168] In this embodiment, the three-dimensional attribute calculation unit obtains the three-dimensional azimuth stacking attribute volume in the following manner:

[0169] perform calculation of coherent attribute, curvature attribute and ant attribute on the azimuth stacking data volume to obtain the three-dimensional azimuth stacking attribute volume.

[0170] In this embodiment, the preferred module 304 comprises:

[0171] a first preferred unit, configured to perform preference on the azimuth pre-stack attribute gather using the predicted fracture development azimuth to obtain preferred azimuth pre-stack attribute gather;

[0172] a second preferred unit configured to prefer the three-dimensional azimuthal stack attribute volume using the predicted fracture development azimuth, to obtain a preferred three-dimensional azimuthal stack attribute volume;

[0173] a preferred attribute acquisition unit configured to acquire the preferred azimuthal pre-stack attribute gather and the preferred three-dimensional azimuthal stack attribute volume as the preferred attributes.

[0174] In this embodiment, the first preferred unit obtains the preferred azimuthal pre-stack attribute gather in the following manner:

[0175] For each sampling point in each seismic trace of the azimuthal pre-stack attribute gather, the following operations are performed:

[0176] attribute traces within a preset angle range are selected from the azimuthal pre-stack attribute gather, to obtain an attribute trace gather at the sampling point in the seismic trace; wherein the preset angle range takes the angle value of the predicted fracture development azimuth as a center point; and the preset angle range is determined based on the density of the sampling points;

[0177] the attribute trace gather at each sampling point in each seismic trace is acquired as the preferred azimuthal pre-stack attribute gather.

[0178] In this embodiment, the second preferred unit obtains the preferred three-dimensional azimuthal stack attribute volume in the following manner:

[0179] For each sampling point in each seismic trace of the three-dimensional azimuthal stack attribute volume, the following operations are performed:

[0180] an attribute volume corresponding to an angle having the smallest absolute value of difference from the angle value of the predicted fracture development azimuth is selected from the three-dimensional azimuthal stack attribute volume, to obtain an attribute volume at the sampling point in the seismic trace;

[0181] the attribute volume at each sampling point in each seismic trace is acquired as the preferred three-dimensional azimuthal stack attribute volume.

[0182] The working principle, working process, and other contents related to the specific embodiments of the above device can be seen from the specific embodiments of the fracture reservoir prediction method provided by the present application, and the same technical contents will not be described in detail here.

[0183] The fracture reservoir prediction device provided by the embodiment of the present application can obtain the predicted fracture development direction based on the azimuth prestack gather, perform prestack fracture anisotropy inversion based on the azimuth prestack gather, obtain the predicted fracture development direction, perform attribute calculation on the azimuth prestack gather, obtain the attribute calculation result, perform optimization on the attribute calculation result by using the predicted fracture development direction, obtain the optimized attribute, and predict the reservoir of the work area based on the optimized attribute, so that the prediction of the reservoir can be performed based on the prestack seismic data, and the prediction method is simplified. Meanwhile, the present application performs prestack fracture anisotropy inversion based on the azimuth prestack gather, and performs optimization on the attribute calculation result by using the predicted fracture development direction, that is, the present application fully considers the fracture development information in the reservoir, and can obtain more accurate prediction result compared with the prior art. It can be seen that the technical scheme provided by the embodiment of the present application can more simply and accurately predict the reservoir.

[0184] Example Four

[0185] According to the embodiment of the present application, a storage medium is also provided, and the storage medium stores program code, and the program code is executed by a processor to implement the fracture reservoir prediction method with azimuth attribute optimization according to any one of the above embodiments.

[0186] The method comprises:

[0187] Performing azimuth migration imaging processing on the prestack CMP gather of the work area to obtain an azimuth prestack gather;

[0188] Performing prestack fracture anisotropy inversion based on the azimuth prestack gather to obtain a predicted fracture development direction;

[0189] Performing attribute calculation on the azimuth prestack gather to obtain an attribute calculation result;

[0190] Performing optimization on the attribute calculation result by using the predicted fracture development direction to obtain an optimized attribute;

[0191] Performing prediction on the reservoir of the work area based on the optimized attribute.

[0192] In the embodiment, the azimuth migration imaging processing on the prestack CMP gather of the work area to obtain an azimuth prestack gather comprises:

[0193] Performing azimuth migration imaging processing on the prestack CMP gather of the work area by using a local angle domain imaging method or a full-azimuth OVT processing method or a split-azimuth migration imaging processing method to obtain an azimuth prestack gather.

[0194] In the embodiment, the attribute calculation on the azimuth pre-stack gather is performed to obtain an attribute calculation result, including:

[0195] The single-attribute calculation is performed on the azimuth pre-stack gather to obtain an azimuth pre-stack attribute gather.

[0196] The azimuth angle stacking is performed on the azimuth pre-stack gather to obtain an azimuth stacking data volume.

[0197] The three-dimensional attribute calculation is performed on the azimuth stacking data volume to obtain a three-dimensional azimuth stacking attribute volume.

[0198] The azimuth pre-stack attribute gather and the three-dimensional azimuth stacking attribute volume are obtained as the attribute calculation result.

[0199] In the embodiment, the single-attribute calculation on the azimuth pre-stack gather is performed to obtain an azimuth pre-stack attribute gather, including:

[0200] The attribute calculation of the instantaneous amplitude, the instantaneous phase, the instantaneous frequency and the root mean square amplitude is performed on each seismic trace of the azimuth pre-stack gather to obtain the azimuth pre-stack attribute gather.

[0201] In the embodiment, the azimuth angle stacking is performed on the azimuth pre-stack gather to obtain an azimuth stacking data volume, including:

[0202] For the pre-stack gather in each azimuth of the azimuth pre-stack gather, the following operations are performed:

[0203] Based on a preset azimuth angle and a preset total offset, the pre-stack gather in the azimuth is stacked to obtain a stacking data volume in the azimuth, wherein the preset azimuth angle is preset based on the number of the pre-stack gather in the azimuth.

[0204] The stacking data volume in each azimuth is obtained as the azimuth stacking data volume.

[0205] In the embodiment, the three-dimensional attribute calculation is performed on the azimuth stacking data volume to obtain a three-dimensional azimuth stacking attribute volume, including:

[0206] The calculation of the coherence attribute, the curvature attribute and the ant attribute is performed on the azimuth stacking data volume to obtain the three-dimensional azimuth stacking attribute volume.

[0207] In the embodiment, the predicted fracture development azimuth is used to optimize the attribute calculation result to obtain an optimized attribute, including:

[0208] The predicted fracture development azimuth is used to optimize the azimuth pre-stack attribute gather to obtain an optimized azimuth pre-stack attribute gather.

[0209] The predicted fracture development direction is used to optimize the three-dimensional azimuth stacking attribute volume, and an optimized three-dimensional azimuth stacking attribute volume is obtained.

[0210] The optimized azimuth pre-stack attribute gather and the optimized three-dimensional azimuth stacking attribute volume are obtained as the optimized attributes.

[0211] In this embodiment, the predicted fracture development direction is used to optimize the azimuth pre-stack attribute gather, and an optimized azimuth pre-stack attribute gather is obtained.

[0212] For each sampling point in each seismic trace of the azimuth pre-stack attribute gather, the following operation is performed:

[0213] An attribute trace within a preset angle range is selected from the azimuth pre-stack attribute gather, and an attribute trace gather at the sampling point in the seismic trace is obtained; the preset angle range has the angle value of the predicted fracture development direction as a center point, and is determined based on the density of the sampling points.

[0214] The attribute trace gather at each sampling point in each seismic trace is obtained as the optimized azimuth pre-stack attribute gather.

[0215] In this embodiment, the predicted fracture development direction is used to optimize the three-dimensional azimuth stacking attribute volume, and an optimized three-dimensional azimuth stacking attribute volume is obtained.

[0216] For each sampling point in each seismic trace of the three-dimensional azimuth stacking attribute volume, the following operation is performed:

[0217] An attribute volume corresponding to an angle having the minimum absolute value of the difference value with the angle value of the predicted fracture development direction is selected from the three-dimensional azimuth stacking attribute volume, and an attribute volume at the sampling point in the seismic trace is obtained.

[0218] The attribute volume at each sampling point in each seismic trace is obtained as the optimized three-dimensional azimuth stacking attribute volume.

[0219] Example Five

[0220] According to an embodiment of the present application, an electronic device is also provided, which comprises a memory and a processor, and the memory stores program codes executable on the processor, and the program codes are executed by the processor to implement the full-azimuth attribute optimization fracture reservoir prediction method according to any one of the above embodiments.

[0221] The fracture reservoir prediction method, device, storage medium and electronic equipment provided by the embodiment of the present application can obtain azimuth prestack gathers by performing azimuth migration imaging processing on the prestack CMP gathers of the work area, perform prestack fracture anisotropy inversion based on the azimuth prestack gathers, obtain a predicted fracture development azimuth, perform attribute calculation on the azimuth prestack gathers, obtain an attribute calculation result, perform optimization on the attribute calculation result by using the predicted fracture development azimuth, obtain an optimized attribute, and predict the reservoir of the work area based on the optimized attribute, so that the prediction of the reservoir can be performed based on prestack seismic data, and it is not necessary to calculate poststack seismic data as in the prior art, and the prediction method is simplified. Meanwhile, the present application performs prestack fracture anisotropy inversion based on azimuth prestack gathers, and performs optimization on the attribute calculation result by using the predicted fracture development azimuth, that is, the present application fully considers the fracture development information in the reservoir, and can obtain more accurate prediction results compared with the prior art. It can be seen that the technical solution provided by the embodiment of the present application can more simply and accurately predict the reservoir.

[0222] The present application aims at the azimuth anisotropy problem commonly existing in actual oil and gas field exploration, performs prestack fracture anisotropy inversion on azimuth prestack gathers, performs attribute calculation on the azimuth prestack gathers, and develops a fracture reservoir prediction method based on azimuth attribute optimization. The method increases the azimuth information of fractures on the basis of conventional attribute information, and enhances the accuracy of identifying geological abnormal bodies such as fractures. The present application takes the prediction of fracture attributes as an example, and can also be extended to the seismic prediction of other geological abnormal bodies (such as fracture-cavity, shale, tight sandstone and igneous rock) with strong heterogeneity.

[0223] In several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other manners. For example, the above-described device embodiments are merely illustrative, and the division of units is merely a logical function division, and actual implementation can be in another manner, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed.

[0224] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, i.e., they can be located in one place or distributed on a plurality of network units. Part or all of the units can be selected according to actual needs to achieve the purpose of the embodiment of the present application.

[0225] In addition, each function unit in each embodiment of the present application can be integrated in one processing unit, or each unit can be physically present separately, or two or more units can be integrated in one unit. The integrated unit can be realized in the form of hardware or in the form of a software function unit.

[0226] When the integrated unit is realized in the form of a software function unit and sold or used as an independent product, it can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the part of the prior art that contributes to the present application, or all or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes a number of instructions for making an electronic device (which can be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the embodiments of the present application. The foregoing storage medium includes a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

[0227] Although the disclosed embodiments of the present application are as above, the content described is only for the purpose of facilitating understanding of the embodiments adopted by the present application, and is not intended to limit the present application. Any person skilled in the art of the present application can make any modification and change in the form of implementation and details without departing from the spirit and scope of the present application, but the protection scope of the present application shall be subject to the scope defined by the appended claims.

Claims

1. A method for azimuthally anisotropic attribute optimization for fracture reservoir prediction, characterized in that, The method comprises: full azimuth migration imaging processing is performed on the pre-stack CMP gathers of the work area to obtain azimuth pre-stack gathers; based on the azimuth pre-stack gathers, pre-stack fracture anisotropy inversion is performed to obtain a predicted fracture development azimuth; attribute calculation is performed on the azimuth pre-stack gathers to obtain attribute calculation results; the attribute calculation results are optimized using the predicted fracture development azimuth to obtain optimized attributes; based on the optimized attributes, reservoirs in the work area are predicted; the attribute calculation on the azimuth pre-stack gathers to obtain attribute calculation results comprises: single-attribute calculation is performed on the azimuth pre-stack gathers to obtain azimuth pre-stack attribute gathers; azimuth angle stacking is performed on the azimuth pre-stack gathers to obtain an azimuth stacking data volume; three-dimensional attribute calculation is performed on the azimuth stacking data volume to obtain a three-dimensional azimuth stacking attribute volume; the azimuth pre-stack attribute gathers and the three-dimensional azimuth stacking attribute volume are obtained as the attribute calculation results; the azimuth angle stacking on the azimuth pre-stack gathers to obtain an azimuth stacking data volume comprises: for the pre-stack gathers in each azimuth of the azimuth pre-stack gathers, the following operations are performed: based on a preset azimuth angle and a preset full offset, the pre-stack gathers in the azimuth are stacked to obtain a stacking data volume in the azimuth; wherein the preset azimuth angle is preset based on the number of the pre-stack gathers in the azimuth; the stacking data volume in each azimuth is obtained as the azimuth stacking data volume; the optimization of the attribute calculation results using the predicted fracture development azimuth to obtain optimized attributes comprises: the azimuth pre-stack attribute gathers are optimized using the predicted fracture development azimuth to obtain optimized azimuth pre-stack attribute gathers; the three-dimensional azimuth stacking attribute volume is optimized using the predicted fracture development azimuth to obtain optimized three-dimensional azimuth stacking attribute volume; the optimized azimuth pre-stack attribute gathers and the optimized three-dimensional azimuth stacking attribute volume are obtained as the optimized attributes.

2. The all azimuth attribute optimization fracture reservoir prediction method of claim 1, wherein, the full azimuth migration imaging processing on the pre-stack CMP gathers of the work area to obtain azimuth pre-stack gathers comprises: full azimuth migration imaging processing is performed on the pre-stack CMP gathers of the work area using a local angle domain imaging method, a full azimuth OVT processing method or a split-azimuth migration imaging processing method to obtain azimuth pre-stack gathers.

3. The all azimuth attribute optimization fracture reservoir prediction method of claim 1, wherein, the single-attribute calculation on the azimuth pre-stack gathers to obtain azimuth pre-stack attribute gathers comprises: instantaneous amplitude, instantaneous phase, instantaneous frequency and root mean square amplitude attribute calculation is performed on each seismic trace of the azimuth pre-stack gathers to obtain the azimuth pre-stack attribute gathers.

4. The all azimuth attribute optimization fracture reservoir prediction method of claim 1, wherein, the three-dimensional attribute calculation on the azimuth stacking data volume to obtain a three-dimensional azimuth stacking attribute volume comprises: coherent attribute, curvature attribute and ant attribute calculation is performed on the azimuth stacking data volume to obtain the three-dimensional azimuth stacking attribute volume.

5. The all azimuth attribute optimization fracture reservoir prediction method of claim 1, wherein, the optimization of the azimuth pre-stack attribute gathers using the predicted fracture development azimuth to obtain optimized azimuth pre-stack attribute gathers comprises: for each sampling point in each seismic trace of the azimuth pre-stack attribute gathers, the following operations are performed: select attribute traces within a preset angle range from the azimuth pre-stack attribute trace set, to obtain the attribute trace set at the sampling point in the seismic trace; wherein the preset angle range takes the angle value of the predicted fracture development azimuth as a center point; and the preset angle range is determined based on the density of the sampling points; obtain the attribute trace set at each sampling point in each seismic trace as the preferred azimuth pre-stack attribute trace set.

6. The all azimuth attribute optimization fracture reservoir prediction method of claim 1, wherein, perform preference on the three-dimensional azimuth stacking attribute volume using the predicted fracture development azimuth, to obtain a preferred three-dimensional azimuth stacking attribute volume, including: for each sampling point in each seismic trace of the three-dimensional azimuth stacking attribute volume, perform the following operations: select the attribute volume corresponding to the angle with the smallest absolute value of the angle value difference from the predicted fracture development azimuth from the three-dimensional azimuth stacking attribute volume, to obtain the attribute volume at the sampling point in the seismic trace; obtain the attribute volume at each sampling point in each seismic trace as the preferred three-dimensional azimuth stacking attribute volume.

7. An all azimuth attribute optimization fracture reservoir prediction apparatus, characterized by, The device comprises: an offset imaging processing module, configured to perform full-azimuth offset imaging processing on a pre-stack CMP trace set of a work area, to obtain an azimuth pre-stack trace set; an inversion module, configured to perform pre-stack fracture anisotropy inversion based on the azimuth pre-stack trace set, to obtain a predicted fracture development azimuth; an attribute calculation module, configured to perform attribute calculation on the azimuth pre-stack trace set, to obtain an attribute calculation result; a preference module, configured to perform preference on the attribute calculation result using the predicted fracture development azimuth, to obtain a preferred attribute; a prediction module, configured to perform prediction on a reservoir of the work area based on the preferred attribute; The attribute calculation module comprises: a single-trace attribute calculation unit, configured to perform single-trace attribute calculation on the azimuth pre-stack trace set, to obtain an azimuth pre-stack attribute trace set; an azimuth angle stacking unit, configured to perform azimuth angle stacking on the azimuth pre-stack trace set, to obtain an azimuth stacking data volume; a three-dimensional attribute calculation unit, configured to perform three-dimensional attribute calculation on the azimuth stacking data volume, to obtain a three-dimensional azimuth stacking attribute volume; an attribute calculation result acquisition unit, configured to obtain the azimuth pre-stack attribute trace set and the three-dimensional azimuth stacking attribute volume as the attribute calculation result; The azimuth angle stacking unit is configured to: for the pre-stack trace set at each azimuth in the azimuth pre-stack trace set, perform the following operations: based on a preset azimuth angle and a preset full offset, stack the pre-stack trace set at the azimuth, to obtain a stacking data volume at the azimuth; wherein the preset azimuth angle is preset based on the number of the pre-stack trace set at the azimuth; obtain the stacking data volume at each azimuth as the azimuth stacking data volume; The preference module comprises: a first preference unit, configured to perform preference on the azimuth pre-stack attribute trace set using the predicted fracture development azimuth, to obtain a preferred azimuth pre-stack attribute trace set; a second preference unit, configured to perform preference on the three-dimensional azimuth stacking attribute volume using the predicted fracture development azimuth, to obtain a preferred three-dimensional azimuth stacking attribute volume; Preferable attribute acquisition unit is used for acquiring the preferred azimuth pre-stack attribute gather and the preferred three-dimensional azimuth pre-stack attribute volume as the preferable attribute.

8. A storage medium having stored thereon a program code, characterized in that, The program code is executed by the processor to implement the full-azimuth attribute optimization fracture reservoir prediction method in any one of claims 1 to 6.

9. An electronic device, comprising: The electronic device includes a memory and a processor, the memory stores program code executable on the processor, and the program code is executed by the processor to implement the full-azimuth attribute optimization fracture reservoir prediction method in any one of claims 1 to 6.

Images