Method and device for acquiring geological layer biological reef information and storage medium

Abstract

The embodiment of the invention provides a method and device for acquiring geological stratum biological reef information and a storage medium. The method comprises the steps that the frequency range of a to-be-measured geological stratum is acquired; dividing the frequency range into a first frequency band, a second frequency band and a third frequency band from low to high; frequency band filtering is carried out on the different frequency bands to obtain energy focusing seismic data of the different frequency bands; respectively calculating energy focusing instantaneous phase data of the energy focusing seismic data of different frequency bands; fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data; and determining the biological reef information of the geological layer to be measured by using the fusion data. According to the method, detailed information of the outline and the internal structure of the biological reef is finely represented.

Description

Technical Field

[0001] The invention relates to the fields of oil and gas exploration and geophysics, and in particular to a method, a device and a storage medium for obtaining geological layer bioherm information. Background Art

[0002] Carbonate bioherms have good physical properties, large development scale, and are easy to form oil and gas accumulations, making them important oil and gas geophysical research objects. Because carbonate bioherms have undergone sedimentary construction, they have some special reflection structure characteristics on seismic reflection profiles.

[0003] The existing methods for obtaining bioherm characterization include planar seismic phase prediction methods based on reflection feature analysis, such as profile geological body carving methods based on amplitude and frequency seismic attributes. However, due to the low resolution of seismic data, the morphology of carbonate bioherms is difficult to accurately depict. Relying solely on a single seismic attribute can only reflect one aspect of the characteristics of carbonate bioherms, and it is difficult to simultaneously depict the specific details of the outer contour and internal structure of carbonate bioherms. Therefore, there is an urgent need for a description technology that can depict the mound-like outer contour and internal structure of carbonate bioherms. Summary of the invention

[0004] The purpose of the embodiments of the present invention is to provide a method, device and storage medium for obtaining bioherm information of a geological layer, wherein the method can simultaneously depict the detailed information of the outer contour and the internal structure of the bioherm.

[0005] In order to achieve the above object, an embodiment of the present invention provides a method for obtaining geological layer bioherm information, the method comprising:

[0006] Obtain the frequency range of the geological layer to be measured;

[0007] According to the frequency range, it is divided into a first frequency band, a second frequency band and a third frequency band from low to high;

[0008] Performing frequency band filtering on different frequency bands respectively to obtain energy focused seismic data of different frequency bands;

[0009] Calculate the energy focusing instantaneous phase data of energy focusing seismic data in different frequency bands respectively;

[0010] fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data;

[0011] The bioherm information of the geological layer to be measured is determined by using the fused data.

[0012] Optionally, the respectively calculating energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands includes:

[0013] Ph(t)=arctan[f(t) / g(t)]

[0014] Among them, Ph(t) is the instantaneous phase data of energy focusing,

[0015] f(t) is the earthquake signal,

[0016] g(t) is the Hilbert transform data of the seismic signal.

[0017] Optionally, fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data includes:

[0018] Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high

[0019] Among them, Phase_ mix To fuse the data,

[0020] Phase_ low is the energy focusing instantaneous phase data of the first frequency band,

[0021] Phase_ mid is the energy focusing instantaneous phase data of the second frequency band,

[0022] Phase_ high The instantaneous phase data of energy focusing in the third frequency band,

[0023] w A is the weighted coefficient of the energy focusing instantaneous phase data of the first frequency band,

[0024] w B is the weighted coefficient of the energy focusing instantaneous phase data of the second frequency band,

[0025] w C is the weighting coefficient of the energy focusing instantaneous phase data of the third frequency band.

[0026] Optionally, the energy focused seismic data in different frequency bands include: energy focused seismic data in a first frequency band, energy focused seismic data in a second frequency band, and energy focused seismic data in a third frequency band;

[0027] The bioherm information includes the outer contour and internal structure of the bioherm.

[0028] Optionally, performing frequency band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands includes:

[0029] Band filtering is performed on the first frequency band, the second frequency band and the third frequency band respectively through corresponding frequency division energy focusing filters;

[0030] The bandwidth of the frequency division energy focusing filter is 60 Hz;

[0031] The energy range of the frequency division energy focusing filter is 10-70 Hz;

[0032] The energy of the frequency division energy focusing filter is focused to the first frequency band, the second frequency band and the third frequency band;

[0033] The main frequency of the first frequency band is 20Hz;

[0034] The main frequency of the first frequency band is 40 Hz;

[0035] The main frequency of the first frequency band is 60 Hz.

[0036] The present invention also provides a device for obtaining bioherm information of geological layers, characterized in that the device comprises:

[0037] An acquisition module, used for acquiring the frequency range of the geological layer to be measured;

[0038] A first processing module, used for dividing the frequency range into a first frequency band, a second frequency band and a third frequency band from low to high;

[0039] The second processing module is used to perform frequency band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands;

[0040] A third processing module is used to respectively calculate energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands;

[0041] A fourth processing module, used for fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data;

[0042] The characterization module is used to determine the bioherm information of the geological layer to be measured by using the fused data.

[0043] Optionally, the respectively calculating energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands includes:

[0044] Ph(t)=arctan[f(t) / g(t)]

[0045] Among them, Ph(t) is the instantaneous phase data of energy focusing,

[0046] f(t) is the earthquake signal,

[0047] g(t) is the Hilbert transform data of the seismic signal.

[0048] Optionally, fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data includes:

[0049] Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high

[0050] Among them, Phase_ mix To fuse the data,

[0051] Phase_ low is the energy focusing instantaneous phase data of the first frequency band,

[0052] Phase_ mid is the energy focusing instantaneous phase data of the second frequency band,

[0053] Phase_ high The instantaneous phase data of energy focusing in the third frequency band,

[0054] w A is the weighted coefficient of the energy focusing instantaneous phase data of the first frequency band,

[0055] w B is the weighted coefficient of the energy focusing instantaneous phase data of the second frequency band,

[0056] w C is the weighting coefficient of the energy focusing instantaneous phase data of the third frequency band.

[0057] Optionally, the energy focused seismic data in different frequency bands include: energy focused seismic data in a first frequency band, energy focused seismic data in a second frequency band, and energy focused seismic data in a third frequency band;

[0058] The bioherm information includes the outer contour and internal structure of the bioherm.

[0059] Optionally, performing frequency band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands includes:

[0060] Band filtering is performed on the first frequency band, the second frequency band and the third frequency band respectively through corresponding frequency division energy focusing filters;

[0061] The bandwidth of the frequency division energy focusing filter is 60 Hz;

[0062] The energy range of the frequency division energy focusing filter is 10-70 Hz;

[0063] The energy of the frequency division energy focusing filter is focused to the first frequency band, the second frequency band and the third frequency band;

[0064] The main frequency of the first frequency band is 20Hz;

[0065] The main frequency of the first frequency band is 40 Hz;

[0066] The main frequency of the first frequency band is 60 Hz.

[0067] The present invention also proposes a machine-readable storage medium, on which instructions are stored, and the instructions are used to enable a machine to execute the above-mentioned method for obtaining bioherm characterization.

[0068] The method for obtaining bioherm information of the geological layer of the present invention includes: obtaining the frequency range of the geological layer to be measured; dividing the frequency range from low to high into a first frequency band, a second frequency band and a third frequency band; performing band filtering on different frequency bands to obtain energy focusing seismic data of different frequency bands; calculating energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands; fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data; and determining the bioherm information of the geological layer to be measured using the fused data. The present invention determines the frequency division energy focusing data bodies of different frequency bands by using seismic data of different frequency bands, and then obtains multi-scale instantaneous phase fusion data by performing RGB fusion on the frequency division energy focusing instantaneous phase data bodies of different frequency bands, thereby finely characterizing the outline and internal structure details of carbonate bioherms.

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

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

[0071] Figure 1 It is a flow chart of a method for obtaining bioherm information of a geological layer according to the present invention;

[0072] Figure 2 It is a schematic diagram of an embodiment of a method for obtaining bioherm information of a geological layer according to the present invention;

[0073] Figure 3 A schematic diagram of the spectrum obtained by performing spectrum scanning on the target layer segment of seismic data;

[0074] Figure 4 It is a schematic diagram of the frequency band characteristics of the low frequency, intermediate frequency and high frequency of the crossover filter;

[0075] Figure 5 It is a schematic diagram of focusing the main frequency of the low-frequency energy focusing filter to the low-frequency dominant frequency;

[0076] Figure 6 It is a schematic diagram of focusing the main frequency of the intermediate frequency energy focusing filter to the intermediate frequency dominant frequency;

[0077] Figure 7 It is a schematic diagram of focusing the main frequency of the high-frequency energy focusing filter to the high-frequency dominant frequency;

[0078] Figure 8 This is a schematic diagram of the original seismic data;

[0079] Fig. 9 This is a schematic diagram of low-frequency energy focusing seismic data;

[0080] Fig.10 This is a schematic diagram of medium frequency energy aggregation seismic data;

[0081] Fig.11 This is a schematic diagram of high-frequency energy focusing seismic data;

[0082] Fig.12 It is a schematic diagram of low-frequency energy focusing instantaneous phase data and bioherm interpretation results;

[0083] Fig.13 It is a schematic diagram of the medium frequency energy focusing instantaneous phase data and the interpretation results of the reef;

[0084] Fig.14 This is a schematic diagram of high-frequency energy focusing instantaneous phase data and bioherm interpretation results;

[0085] Fig.15 Schematic diagram of multi-scale instantaneous phase fusion data and reef interpretation results. DETAILED DESCRIPTION

[0086] The specific implementation of the embodiment of the present invention is described in detail below in conjunction with the accompanying drawings. It should be understood that the specific implementation described here is only used to illustrate and explain the embodiment of the present invention, and is not used to limit the embodiment of the present invention.

[0087] Embodiment 1:

[0088] Figure 1 It is a flow chart of a method for obtaining bioherm information of geological layers of the present invention, such as Figure 1As shown, the method for obtaining geological layer bioherm information of the present invention comprises: step S101 is to obtain the frequency range of the geological layer to be measured. The geological layer to be measured can be the target layer segment of seismic data, and the frequency spectrum of the target layer segment of seismic data is scanned to determine the frequency range.

[0089] Step S102 is to divide the frequency range from low to high into a first frequency band, a second frequency band and a third frequency band. For the convenience of description, the first frequency band is uniformly referred to as a low frequency band, the second frequency band is referred to as a medium frequency band, and the third frequency band is referred to as a high frequency band. The range of the low frequency band is 10 to 30 Hz, the range of the medium frequency band is 30 to 50 Hz, and the range of the high frequency band is 50 to 70 Hz.

[0090] The frequency division energy focusing filter is divided into low frequency, medium frequency and high frequency, and the bandwidth (10-70Hz, i.e. 60Hz) is maintained, and the energy is focused on low frequency (main frequency 20Hz), medium frequency (main frequency 40Hz) and high frequency (main frequency 60Hz). Specifically, the low frequency energy range is: 10-70Hz, the bandwidth is 60Hz, the main frequency is 20Hz (such as Figure 5 The intermediate frequency energy range is 10-70Hz, the bandwidth is 60Hz, and the main frequency is 40Hz (as shown in Figure 6 As shown); high frequency energy range is: 10-70Hz, bandwidth 60Hz, main frequency 60Hz (as shown Figure 7 The frequency division energy focusing filter has a bandwidth of 60 Hz, which not only achieves the purpose of frequency division, but also each filter retains the full-band information of 60 Hz.

[0091] Step S103 is to perform band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands. Specifically, the energy focused seismic data of different frequency bands include: low-frequency energy focused seismic data, medium-frequency energy focused seismic data and high-frequency energy focused seismic data; the bioherm information includes the outer contour and internal structure of the bioherm. The band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands includes: band filtering on the different frequency bands respectively through corresponding frequency division energy focused filters.

[0092] According to a specific implementation method, three frequency-dividing energy focusing filters are used to perform frequency band filtering on the seismic data respectively to obtain first frequency band energy focusing seismic data, second frequency band energy focusing seismic data and third frequency band energy focusing seismic data, namely low-frequency energy focusing seismic data, medium-frequency energy focusing seismic data and high-frequency energy focusing seismic data.

[0093] Step S104 is to calculate the energy focusing instantaneous phase data of the energy focusing seismic data of different frequency bands respectively. The energy focusing instantaneous phase data specifically includes low-frequency energy focusing seismic data, medium-frequency energy focusing seismic data and high-frequency energy focusing seismic data. The following formula is used for calculation, and the data body obtained is the low-frequency energy focusing seismic instantaneous phase data, the medium-frequency energy focusing seismic instantaneous phase data and the high-frequency energy focusing seismic instantaneous phase data.

[0094] Specifically, the energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands are calculated respectively.

[0095] Ph(t)=arctan[f(t) / g(t)]

[0096] Among them, Ph(t) is the energy focusing instantaneous phase data, f(t) is the seismic signal, and g(t) is the Hilbert transform data of the seismic signal.

[0097] Step S105 is to fuse the energy focusing instantaneous phase data of the different frequency bands to obtain fused data. Specifically, the energy focusing instantaneous phase data of the different frequency bands are fused to obtain fused data. The fusion method is preferably weighted average fusion, including:

[0098] Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high in,

[0099] Phase_ mix To fuse the data, Phase_ low is the energy focusing instantaneous phase data of the first frequency band (i.e., low frequency band), Phase_ mid is the energy focusing instantaneous phase data of the second frequency band (i.e., the middle frequency band), Phase_ high is the energy focusing instantaneous phase data of the third frequency band (i.e., high frequency band), w A is the weighted coefficient of the energy focusing instantaneous phase data of the first frequency band (i.e. low frequency), w B is the weighted coefficient of the energy focusing instantaneous phase data of the second frequency band (i.e., the middle frequency), w C is the weighted coefficient of the energy focusing instantaneous phase data of the third frequency band (ie, high frequency). The fused data is multi-scale instantaneous phase fusion data.

[0100] Step S106 is to determine the bioherm information of the geological layer to be tested by using the fused data. Specifically, the bioherm information includes the outer contour and internal structure of the bioherm. The low-frequency energy focused instantaneous phase data volume reflects the outer contour of the carbonate bioherm, and the high-frequency energy focused instantaneous phase data volume reflects the internal structure details of the carbonate bioherm.

[0101] The present invention proposes a method for obtaining geological layer bioherm information. On the basis of determining the effective frequency range of the target layer of seismic data, frequency division energy focusing filters of three different frequency bands, low frequency, medium frequency and high frequency, are designed. The characteristic of the frequency division energy focusing filter is that the main frequency is focused to low frequency, medium frequency and high frequency respectively while retaining the full-band information, overcoming the defects of conventional frequency division technology (single frequency body) with narrow frequency band and unclear geological significance. In addition, the instantaneous phase attribute is calculated on the basis of frequency division energy focusing, so that the low-frequency energy focusing instantaneous phase data body reflects the outer contour of carbonate bioherms, and the high-frequency energy focusing instantaneous phase data body reflects the internal structure details of carbonate bioherms. High-resolution image processing is achieved through the RGB fusion method, which effectively identifies image detail information and highlights anomalies. Multi-scale instantaneous phase fusion data is obtained through weighted average fusion, which achieves the effect of clearly depicting the outer contour and internal structure detail information of carbonate bioherms at the same time.

[0102] Embodiment 2:

[0103] Figure 2 FIG. 1 is a schematic diagram of an embodiment of a method for obtaining bioherm information of a geological layer according to the present invention, such as Figure 2 As shown, step S1 is to perform spectrum scanning on the target layer segment of the seismic data to determine the frequency range.

[0104] Specifically, the target layer segment of the seismic data volume is spectrally scanned, and the target layer segment is used to study the target top and bottom boundary layers. Usually, the layer segment of the spectrum scan is not less than two sub-wavelengths, that is, not less than 200ms, so as to ensure the stability of the scanned spectrum. When the thickness of the target layer segment is small, it can be appropriately extended up and down to the top and bottom layers. The spectrum obtained by scanning is as follows Figure 3 As shown in the figure (the horizontal axis is frequency, in Hz, and the vertical axis is normalized seismic amplitude, with a value of 0 to 1, dimensionless), with the normalized seismic amplitude of 0.25 as the boundary, the effective frequency band unit is 10 to 70 Hz, that is, the low cutoff frequency is Flow1 = 10 Hz, and the high cutoff frequency is Fhigh1 = 70 Hz. With the normalized seismic amplitude of 0.5 as the boundary, the dominant frequency band range of seismic data is obtained, that is, the low frequency dominant frequency Flow2 = 20 Hz, the high frequency dominant frequency is Fhigh2 = 60 Hz, and the mid-frequency dominant frequency range of seismic data, that is, the main frequency of seismic data (the frequency at the maximum normalized seismic amplitude) is Fmid = 40 Hz.

[0105] Step S2 is to design frequency division energy focusing filters for three different frequency bands: low frequency, medium frequency and high frequency.

[0106] Specifically, conventional frequency division technology filters usually select a range of 10 Hz around the dominant frequency for filtering, that is, a single-frequency body, such as Figure 4 As shown: the low-cut, low-pass, high-pass and high-cut frequencies of the low-frequency energy focusing filter (Filter_low) are [Flow2-10Hz, Flow2, Flow2, Flow2+10Hz] respectively; the low-cut, low-pass, high-pass and high-cut frequencies of the medium-frequency energy focusing filter (Filter_mid) are [Fmid-10Hz, Fmid, Fmid, Fmid+10Hz] respectively; the low-cut, low-pass, high-pass and high-cut frequencies of the high-frequency energy focusing filter (Filter_high) are [Fhigh2-10Hz, Fhigh2, Fhigh2, Fhigh2+10Hz] respectively. Among them, the bandwidth of the conventional filter is 20Hz, which is relatively narrow overall. In order to overcome the defect of the narrow bandwidth of the conventional frequency division filter, the present invention proposes a frequency division energy focusing filter, and the specific method is: the low-cut, low-pass, high-pass and high-cut frequencies of the low-frequency energy focusing filter are [Flow1, Flow2, Flow2, Fhigh1] respectively. As Figure 5 As shown, the low-cut, low-pass, high-pass and high-cut frequencies of the intermediate frequency energy focusing filter are [Flow1, Fmid, Fmid, Fhigh1] respectively. Figure 6 As shown; the low-cut, low-pass, high-pass and high-cut frequencies of the high-frequency energy focusing filter are [Flow1, Fhigh2, Fhigh2, Fhigh1] respectively. Figure 7 As shown, the advantage of the frequency division energy focusing filter is that the main frequency is focused to the low frequency, medium frequency and high frequency respectively, thus achieving the purpose of frequency division, but at the same time each filter retains the full-band information, thus overcoming the defect of narrow frequency band of conventional frequency division technology.

[0107] The frequency band of the original earthquake is 10-70Hz, and the bandwidth is 60Hz. The conventional frequency division filter has a bandwidth of only 20Hz, the frequency band is too narrow, and the geological significance is unclear. The frequency division energy focusing filter proposed in the present invention is also divided into low frequency, medium frequency and high frequency, but it maintains the bandwidth (10-70Hz, i.e. 60Hz), and the energy is focused on low frequency (main frequency 20Hz), medium frequency (main frequency 40Hz) and high frequency (main frequency 60Hz). Specifically, the low-frequency energy range is: 10-70Hz, bandwidth 60Hz, main frequency 20Hz; the medium-frequency energy range is: 10-70Hz, bandwidth 60Hz, main frequency 40Hz; the high-frequency energy range is: 10-70Hz, bandwidth 60Hz, main frequency 60Hz. The bandwidth of the present invention is 60Hz, which not only achieves the purpose of frequency division, but also each filter retains the full-band information of 60Hz.

[0108] Step S3 is to use three frequency division energy focusing filters to perform frequency band filtering on the seismic data respectively.

[0109] Specifically, the three frequency-divided energy focusing filters designed by S2 are used to perform frequency band filtering on the seismic data respectively, and three frequency-divided energy focusing data bodies are obtained, namely, low-frequency energy focusing seismic data (Seis_low), medium-frequency energy focusing seismic data (Seis_mid) and high-frequency energy focusing seismic data (Seis_high). This embodiment is a carbonate bioherm, Figure 8 is the original earthquake number (Seis), and the triangle in the figure is the outline of the reef. Fig. 9 The low-frequency energy focused seismic data (Seis_low) obtained by filtering the original seismic data with the low-frequency energy focusing filter in step 2, Fig.10 The intermediate frequency energy focusing seismic data (Seis_mid) is obtained by filtering the original seismic data using the S2 intermediate frequency energy focusing filter. Fig.11 The high-frequency energy focused seismic data (Seis_high) is obtained by filtering the original seismic data using the high-frequency energy focusing filter in step 2.

[0110] Step S4 is to calculate the instantaneous phase attributes of the three frequency-division energy focusing seismic data volumes respectively, to obtain three frequency-division energy focusing instantaneous phase data volumes.

[0111] Specifically, the physical meaning of the instantaneous phase is the angle between the phase vector as a function of time and the real axis. The instantaneous phase attribute is independent of seismic energy. It can enhance weak reflection signals and highlight the lateral variation characteristics of the event axis. It has a good effect on characterizing the lateral boundary characteristics of special carbonate reefs.

[0112] Complex seismic traces are used to calculate instantaneous phase attributes, that is, first calculate the complex seismic traces, and then calculate the instantaneous phase attributes from the complex seismic traces. The instantaneous phase attribute is an instantaneous attribute derived from the complex seismic trace attributes. The complex seismic trace attributes are the attributes obtained by Hilbert transforming the original seismic traces. The seismic signal and its Hilbert transform are expressed as f(t) and g(t) respectively, and the complex seismic trace is expressed as:

[0113] F(t) = f(t) + ig(t) (Formula 1) The expression of the instantaneous phase property is:

[0114] Ph(t)=arctan[f(t) / g(t)] (Formula 2)

[0115] The instantaneous phase attributes of the three frequency-divided energy focusing data bodies are calculated respectively to obtain three frequency-divided energy focusing instantaneous phase data, namely low-frequency energy focusing instantaneous phase data (Phase_low), medium-frequency energy focusing instantaneous phase data (Phase_mid) and high-frequency energy focusing instantaneous phase data (Phase_high).

[0116] Fig.12 It is the low-frequency energy focusing instantaneous phase data (Phase_low) and the interpretation results of bioherms. Fig.12 The middle right picture is the result of the reef interpretation. The outline features of the reef body can be clearly identified from the picture. Due to the low resolution, the details of the internal structure of the reef are not explained in detail.

[0117] Fig.13 The intermediate frequency energy focusing instantaneous phase data (Phase_mid) and the reef interpretation results are shown in Figure 2. Fig.13 The middle right picture shows the interpretation result of the reef. It can be seen from the picture that the outline of the reef is not as clear as the low-frequency energy focusing instantaneous phase data, but the accuracy of the internal structure details of the reef has been improved.

[0118] Fig.14 It is the high-frequency energy focusing instantaneous phase data (Phase_high) and the interpretation results of bioherms. Fig.14 The middle right picture shows the interpretation result of the reef. From the picture, it can be seen that the outline of the reef becomes unclear, but the detailed information of the internal structure of the reef can be clearly depicted.

[0119] Step S5 is to perform weighted average fusion on the three frequency-divided energy focusing instantaneous phase data volumes to obtain multi-scale instantaneous phase fusion data.

[0120] This application prefers weighted average fusion. Reasonable data fusion strategy is the key to obtain high-quality fused images. Weighted average fusion is relatively simple to execute. It not only has low system overhead, but also can obtain results in real time. At the same time, this method can suppress noise in the original data. In many applications, this data fusion method can often achieve satisfactory fusion effects.

[0121] In the present invention, the energy focusing instantaneous phase data of different frequencies are obtained by frequency division energy focusing, and the energy focusing instantaneous phase data of different frequencies are fused by weighted average fusion method (Formula 3), where w A 、w B 、w C The weighted coefficients of the energy focusing instantaneous phase data of different frequencies are respectively used. In this example, all three weight coefficients are 0.33. After weighting, multi-scale instantaneous phase fusion data (Phase_mix) is obtained, thereby retaining the significant features of the original image in different frequency domains in the synthetic image.

[0122] Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high (Formula 3)

[0123] Step S6 is to use multi-scale instantaneous phase fusion data to depict the outline and internal structure details of the carbonate reef.

[0124] like Fig.15 As shown in the figure, the multi-scale instantaneous phase fusion data (Phase_mix) can clearly depict the external contour features of the bioherm, as well as the internal structure of the bioherm, and can finely divide the three-stage bioherms, thus achieving the purpose of clearly depicting the external contour and internal structure details of the carbonate bioherm at the same time.

[0125] Embodiment three:

[0126] The present invention also proposes a device for obtaining bioherm information of a geological layer, the device comprising: an acquisition module for obtaining the frequency range of a geological layer to be measured; a first processing module for dividing the frequency range into a first frequency band, a second frequency band and a third frequency band from low to high; a second processing module for performing frequency band filtering on different frequency bands respectively to obtain energy focusing seismic data of different frequency bands; a third processing module for respectively calculating energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands; a fourth processing module for fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data; and a characterization module for determining the bioherm information of the geological layer to be measured using the fused data.

[0127] Specifically, the energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands are calculated respectively, including:

[0128] Ph(t)=arctan[f(t) / g(t)]

[0129] Among them, Ph(t) is the energy focusing instantaneous phase data, f(t) is the seismic signal, and g(t) is the Hilbert transform data of the seismic signal.

[0130] The fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data includes:

[0131] Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high

[0132] Among them, Phase_ mix To fuse the data, Phase_ low is the energy focusing instantaneous phase data of the first frequency band, Phase_ mid is the instantaneous phase data of energy focusing in the second frequency band, Phase_ high is the energy focusing instantaneous phase data of the third frequency band, w A is the weighted coefficient of the energy focusing instantaneous phase data of the first frequency band, w B is the weighted coefficient of the energy focusing instantaneous phase data of the second frequency band, w Cis the weighted weight coefficient of the energy focusing instantaneous phase data of the third frequency band. The energy focusing seismic data of different frequency bands include: low-frequency energy focusing seismic data, medium-frequency energy focusing seismic data and high-frequency energy focusing seismic data; the bioherm information includes the outer contour and internal structure of the bioherm. The step of performing frequency band filtering on different frequency bands to obtain energy focusing seismic data of different frequency bands includes: performing frequency band filtering on the different frequency bands respectively through corresponding frequency division energy focusing filters.

[0133] The method for obtaining bioherm information of the geological layer of the present invention includes: obtaining the frequency range of the geological layer to be measured; dividing the frequency range from low to high into a first frequency band, a second frequency band and a third frequency band; performing band filtering on different frequency bands to obtain energy focusing seismic data of different frequency bands; calculating energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands; fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data; and determining the bioherm information of the geological layer to be measured using the fused data. The present invention determines the frequency division energy focusing data bodies of different frequency bands by using seismic data of different frequency bands, and then obtains multi-scale instantaneous phase fusion data by performing RGB fusion on the frequency division energy focusing instantaneous phase data bodies of different frequency bands, thereby finely characterizing the outline and internal structure details of carbonate bioherms.

[0134] On the other hand, the present invention also provides an electronic device, which includes: at least one processor; a memory connected to the at least one processor; wherein the memory stores instructions that can be executed by the at least one processor, and the at least one processor implements the above-mentioned method of obtaining geological layer bioherm information by executing the instructions stored in the memory.

[0135] The processor includes a kernel, which retrieves the corresponding program unit from the memory. One or more kernels can be set, and the detailed information of the outer contour and internal structure of the reef can be finely characterized by adjusting the kernel parameters.

[0136] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM, and the memory includes at least one memory chip.

[0137] An embodiment of the present invention provides a storage medium on which a program is stored. When the program is executed by a processor, the oil production method is implemented.

[0138] An embodiment of the present invention provides a processor, which is used to run a program, wherein the method for obtaining geological layer bioherm information is executed when the program is run.

[0139] The embodiment of the present invention provides a device, which includes a processor, a memory, and a program stored in the memory and executable on the processor. When the processor executes the program, the following steps are implemented: (method claim steps, sole right + subordinate right). The device in this article can be a server, a PC, a PAD, a mobile phone, etc.

[0140] The present application also provides a computer program product, which, when executed on a data processing device, is suitable for executing an initialization program having the following method steps: obtaining the frequency range of the geological layer to be measured; dividing the frequency range into a first frequency band, a second frequency band and a third frequency band from low to high; performing frequency band filtering on different frequency bands to obtain energy focusing seismic data of different frequency bands; calculating energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands; fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data; and determining the bioherm information of the geological layer to be measured using the fused data.

[0141] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0142] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0143] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0144] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0145] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0146] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0147] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0148] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0149] The above are only embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should be included within the scope of the claims of the present application.

Claims

1. A method for obtaining bioherm information of geological layers, characterized in that: The method includes: Obtain the frequency range of the geological layer to be measured; According to the frequency range, it is divided into a first frequency band, a second frequency band and a third frequency band from low to high; Performing frequency band filtering on different frequency bands respectively to obtain energy focused seismic data of different frequency bands; Calculate the energy focusing instantaneous phase data of energy focusing seismic data in different frequency bands respectively; fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data; The bioherm information of the geological layer to be measured is determined by using the fused data.

2. The method according to claim 1, characterized in that: The energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands are calculated respectively, including: Ph(t)=arctan[f(t) / g(t)] Among them, Ph(t) is the instantaneous phase data of energy focusing, f(t) is the earthquake signal, g(t) is the Hilbert transform data of the seismic signal.

3. The method according to claim 1, characterized in that The energy focusing instantaneous phase data of the different frequency bands are fused to obtain fused data, including: Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high Among them, Phase_ mix To fuse the data, Phase_ low is the energy focusing instantaneous phase data of the first frequency band, Phase_ mid is the energy focusing instantaneous phase data of the second frequency band, Phase_ high The instantaneous phase data of energy focusing in the third frequency band, w A is the weighted coefficient of the energy focusing instantaneous phase data of the first frequency band, w B is the weighted coefficient of the energy focusing instantaneous phase data of the second frequency band, w C is the weighting coefficient of the energy focusing instantaneous phase data of the third frequency band.

4. The method according to claim 1, characterized in that: The energy-focused seismic data of different frequency bands include: energy-focused seismic data of a first frequency band, energy-focused seismic data of a second frequency band, and energy-focused seismic data of a third frequency band; The bioherm information includes the outer contour and internal structure of the bioherm.

5. The method according to claim 1, characterized in that The method of performing frequency band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands includes: Band filtering is performed on the first frequency band, the second frequency band, and the third frequency band respectively through corresponding frequency division energy focusing filters; The bandwidth of the frequency division energy focusing filter is 60 Hz; The energy range of the frequency division energy focusing filter is 10-70 Hz; The energy of the frequency division energy focusing filter is focused to the first frequency band, the second frequency band and the third frequency band; The main frequency of the first frequency band is 20Hz; The main frequency of the first frequency band is 40 Hz; The main frequency of the first frequency band is 60 Hz.

6. A device for obtaining information about biological reefs in geological layers, characterized in that: The device includes: An acquisition module, used for acquiring the frequency range of the geological layer to be measured; A first processing module, configured to divide the frequency range into a first frequency band, a second frequency band and a third frequency band from low to high; The second processing module is used to perform frequency band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands; A third processing module is used to respectively calculate energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands; A fourth processing module, used for fusing the energy focusing instantaneous phase data of different frequency bands to obtain fused data; The characterization module is used to determine the bioherm information of the geological layer to be measured by using the fused data.

7. The device according to claim 6, characterized in that The energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands are calculated respectively, including: Ph(t)=arctan[f(t) / g(t)] Among them, Ph(t) is the instantaneous phase data of energy focusing, f(t) is the earthquake signal, g(t) is the Hilbert transform data of the seismic signal.

8. The device according to claim 6, characterized in that The fusing the energy focusing instantaneous phase data of the different frequency bands to obtain fused data includes: Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high Among them, Phase_ mix To fuse the data, Phase_ low is the energy focusing instantaneous phase data of the first frequency band, Phase_ mid is the energy focusing instantaneous phase data of the second frequency band, Phase_ high The instantaneous phase data of energy focusing in the third frequency band, w A is the weighted coefficient of the energy focusing instantaneous phase data of the first frequency band, w B is the weighted coefficient of the energy focusing instantaneous phase data of the second frequency band, w C is the weighting coefficient of the energy focusing instantaneous phase data of the third frequency band.

9. The device according to claim 6, characterized in that The energy-focused seismic data of different frequency bands include: energy-focused seismic data of a first frequency band, energy-focused seismic data of a second frequency band, and energy-focused seismic data of a third frequency band; The bioherm information includes the outer contour and internal structure of the bioherm.

10. The device according to claim 6, characterized in that The method of performing frequency band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands includes: Band filtering is performed on the first frequency band, the second frequency band, and the third frequency band respectively through corresponding frequency division energy focusing filters; The bandwidth of the frequency division energy focusing filter is 60 Hz; The energy range of the frequency division energy focusing filter is 10-70 Hz; The energy of the frequency division energy focusing filter is focused to the first frequency band, the second frequency band and the third frequency band; The main frequency of the first frequency band is 20Hz; The main frequency of the first frequency band is 40 Hz; The main frequency of the first frequency band is 60 Hz.

11. A machine-readable storage medium having instructions stored thereon, the instructions being used to enable a machine to execute the method for obtaining geological layer bioherm information as described in any one of claims 1 to 5.

Images