Method and apparatus for acquiring geological formation biogenic reef information, and storage medium

Abstract

A method for acquiring geological formation biogenic reef information, comprising: acquiring a frequency range of a geological formation to be detected (S101); on the basis of the frequency range, dividing a first frequency band, a second frequency band, and a third frequency band from low to high (S102); respectively performing frequency band filtering on different frequency bands to obtain energy focused seismic data of different frequency bands (S103); respectively calculating energy focused instantaneous phase data of the energy focused seismic data of different frequency bands (S104); fusing the energy focused instantaneous phase data of different frequency bands to obtain fused data (S105); and using the fused data to determine biogenic reef information of said geological formation (S106). Detailed information of the contour and the internal structure of biogenic reefs is finely characterized. Further provided are an apparatus for acquiring geological formation biogenic reef information and a storage medium.

Description

Method, device and storage medium for obtaining geological layer bioherm information

[0001] CROSS-REFERENCE TO RELATED APPLICATIONS

[0002] This application claims the benefit of Chinese patent application No. 202311220869.2 filed on September 20, 2023, the contents of which are incorporated herein by reference. Technical Field

[0003] The present 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

[0004] Carbonate bioherms have good physical properties, are large in scale, and are prone to forming oil and gas accumulations, making them important targets for oil and gas geophysical research. Because carbonate bioherms have undergone sedimentary construction, they exhibit some unique reflection structural characteristics on seismic reflection profiles.

[0005] Existing methods for characterizing reefs include planar seismic phase prediction methods based on reflection feature analysis, such as profile carving methods based on seismic attributes such as amplitude and frequency. However, due to the low resolution of seismic data, the morphology of carbonate reefs is difficult to accurately characterize. Relying solely on a single seismic attribute can only reflect a single aspect of a carbonate reef's characteristics, making it difficult to simultaneously depict both its external contours and the specific details of its internal structure. Therefore, a descriptive technology that can characterize both the mound-like contours and internal structure of carbonate reefs is urgently needed.

[0006] Summary of the Invention

[0007] 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 representing the outer contour and internal structure of the bioherm.

[0008] In order to achieve the above-mentioned purpose, an embodiment of the present invention provides a method for obtaining bioherm information of a geological layer, the method comprising: obtaining a 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-focused seismic data of different frequency bands; calculating energy-focused instantaneous phase data of the energy-focused seismic data of different frequency bands; fusing the energy-focused 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.

[0009] Optionally, the energy-focused instantaneous phase data of the energy-focused seismic data of different frequency bands are calculated separately, including: Ph(t) = arctan[g(t) / f(t)]

[0010] 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.

[0011] Optionally, 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

[0012] Among them, Phase _mix To fuse data, Phase _low The energy focusing instantaneous phase data of the first frequency band, Phase _mid 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 C The weighting coefficient for the energy focusing instantaneous phase data of the third frequency band.

[0013] Optionally, 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; and the reef information includes the outer contour and internal structure of the reef.

[0014] Optionally, the band filtering of different frequency bands to obtain energy-focused seismic data of different frequency bands includes: band filtering 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 60Hz; the energy range of the frequency-division energy focusing filter is 10-70Hz; 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 second frequency band is 40Hz; the main frequency of the third frequency band is 60Hz.

[0015] The present invention also proposes a device for obtaining bioherm information of a geological layer, characterized in that the device includes: an acquisition module for obtaining the frequency range of the 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-focused seismic data of different frequency bands; a third processing module for respectively calculating the energy-focused instantaneous phase data of the energy-focused seismic data of different frequency bands; a fourth processing module for fusing the energy-focused instantaneous phase data of the different frequency bands to obtain fused data; and a characterization module for using the fused data to determine the bioherm information of the geological layer to be measured.

[0016] Optionally, the energy-focused instantaneous phase data of the energy-focused seismic data of different frequency bands are calculated separately, including: Ph(t) = arctan[g(t) / f(t)]

[0017] 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.

[0018] Optionally, 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

[0019] Among them, Phase _mix To fuse data, Phase _low The energy focusing instantaneous phase data of the first frequency band, Phase _mid 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 C The weighting coefficient for the energy focusing instantaneous phase data of the third frequency band.

[0020] Optionally, 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; and the reef information includes the outer contour and internal structure of the reef.

[0021] Optionally, the band filtering of different frequency bands to obtain energy-focused seismic data of different frequency bands includes: band filtering 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 60Hz; the energy range of the frequency-division energy focusing filter is 10-70Hz; 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 second frequency band is 40Hz; the main frequency of the third frequency band is 60Hz.

[0022] The present invention also provides a machine-readable storage medium having instructions stored thereon, the instructions being used to enable a machine to execute the above-mentioned method for obtaining bioherm characterization.

[0023] The method of obtaining bioherm information of a geological layer of the present invention includes: obtaining a 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 according to the frequency range; performing frequency band filtering on different frequency bands to obtain energy-focused seismic data of different frequency bands; calculating energy-focused instantaneous phase data of the energy-focused seismic data of different frequency bands; fusing the energy-focused 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 frequency-divided energy-focused data volumes of different frequency bands using seismic data of different frequency bands, and then performs RGB fusion on the frequency-divided energy-focused instantaneous phase data volumes of different frequency bands to obtain multi-scale instantaneous phase fusion data, thereby finely characterizing the outline and internal structure details of carbonate bioherms.

[0024] 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

[0025] 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 detailed description, they are used to explain the embodiments of the present invention, but do not constitute a limitation of the embodiments of the present invention. In the accompanying drawings:

[0026] FIG1 is a schematic flow chart of a method for obtaining bioherm information of a geological layer according to the present invention;

[0027] FIG2 is a schematic diagram of an embodiment of a method for obtaining bioherm information of a geological layer according to the present invention;

[0028] FIG3 is a schematic diagram of a spectrum obtained by performing spectrum scanning on a target layer segment of seismic data;

[0029] FIG4 is a schematic diagram showing the frequency band characteristics of the low frequency, intermediate frequency and high frequency of the crossover filter;

[0030] FIG5 is a schematic diagram showing the main frequency of the low-frequency energy focusing filter focusing on the low-frequency dominant frequency;

[0031] FIG6 is a schematic diagram of an intermediate frequency energy focusing filter focusing the main frequency to the intermediate frequency dominant frequency;

[0032] FIG7 is a schematic diagram of a high-frequency energy focusing filter focusing the main frequency to a high-frequency dominant frequency;

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

[0034] Figure 9 is a schematic diagram of low-frequency energy focusing seismic data;

[0035] Figure 10 is a schematic diagram of medium frequency energy concentration seismic data;

[0036] Figure 11 is a schematic diagram of high-frequency energy focused seismic data;

[0037] Figure 12 is a schematic diagram of low-frequency energy focusing instantaneous phase data and reef interpretation results;

[0038] Figure 13 is a schematic diagram of the medium frequency energy focusing instantaneous phase data and the reef interpretation results;

[0039] Figure 14 is a schematic diagram of high-frequency energy focusing instantaneous phase data and reef interpretation results;

[0040] Figure 15 is a schematic diagram of multi-scale instantaneous phase fusion data and reef interpretation results. DETAILED DESCRIPTION

[0041] The following describes the specific implementation of the embodiment of the present invention in detail with reference to the accompanying drawings. It should be understood that the specific implementation described herein 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.

[0042] Example 1:

[0043] FIG1 is a flow chart of a method for obtaining bioherm information from a geological layer according to the present invention. As shown in FIG1 , the method for obtaining bioherm information from a geological layer according to the present invention includes: Step S101 of obtaining a frequency range of a geological layer to be measured. The geological layer to be measured may be a target layer segment of seismic data, and a frequency spectrum scan is performed on the target layer segment of the seismic data to determine the frequency range.

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

[0045] The frequency-division energy focusing filter is divided into low frequency, medium frequency and high frequency, maintaining the bandwidth (10-70Hz, i.e. 60Hz), focusing energy to 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 (as shown in Figure 5); the medium-frequency energy range is: 10-70Hz, bandwidth 60Hz, main frequency 40Hz (as shown in Figure 6); the high-frequency energy range is: 10-70Hz, bandwidth 60Hz, main frequency 60Hz (as shown in Figure 7). The bandwidth of the frequency-division energy focusing filter is 60Hz, which not only achieves the purpose of frequency division, but also each filter retains the full-band information of 60Hz.

[0046] Step S103 is to perform frequency band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands. Specifically, the energy-focused seismic data for different frequency bands includes low-frequency energy-focused seismic data, medium-frequency energy-focused seismic data, and high-frequency energy-focused seismic data; the reef information includes the outer contours and internal structure of the reef. Performing frequency band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands includes performing frequency band filtering on the different frequency bands using corresponding frequency-dividing energy-focusing filters.

[0047] 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-band energy-focused seismic data, second-band energy-focused seismic data and third-band energy-focused seismic data, namely, low-frequency energy-focused seismic data, medium-frequency energy-focused seismic data and high-frequency energy-focused seismic data.

[0048] Step S104 calculates the energy-focused instantaneous phase data for the energy-focused seismic data in different frequency bands. The energy-focused instantaneous phase data specifically includes low-frequency energy-focused seismic data, medium-frequency energy-focused seismic data, and high-frequency energy-focused seismic data. The following formula is used to calculate the resulting data volume: low-frequency energy-focused seismic instantaneous phase data, medium-frequency energy-focused seismic instantaneous phase data, and high-frequency energy-focused seismic instantaneous phase data.

[0049] Specifically, the energy-focused instantaneous phase data of the energy-focused seismic data of different frequency bands are calculated separately. Ph(t) = arctan[g(t) / f(t)]

[0050] 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.

[0051] Step S105 is to fuse the instantaneous phase data of energy focus of different frequency bands to obtain fused data. Specifically, the instantaneous phase data of energy focus of different frequency bands are fused to obtain fused data. The fusion method is preferably weighted average fusion, including: Phase _mix =w A *Phase _low +w B *Phase _mid +w C *Phase _high Among them, Phase _mix To fuse 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.

[0052] Step S106 uses the fused data to determine the bioherm information of the geological layer being measured. Specifically, the bioherm information includes the bioherm's outline and internal structure. The low-frequency energy focused instantaneous phase data volume reflects the outline of the carbonate bioherm, while the high-frequency energy focused instantaneous phase data volume reflects the internal structural details of the carbonate bioherm.

[0053] 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 segment 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 attributes are calculated based on the 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 structural 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 structural detail information of carbonate bioherms at the same time.

[0054] Example 2:

[0055] FIG2 is a schematic diagram of an embodiment of a method for obtaining bioherm information of a geological layer according to the present invention. As shown in FIG2 , step S1 is to perform spectrum scanning on a target layer segment of seismic data to determine a frequency range.

[0056] Specifically, a spectrum scan is performed on the target layer segment of the seismic data volume. The target layer segment is used to study the target top and bottom boundary layers. Typically, the spectrum scan segment is no less than two wavelet lengths, that is, no less than 200ms, to ensure the stability of the scanned spectrum. When the target layer segment is thin, it can be appropriately extended up and down to the top and bottom layers. The resulting spectrum is shown in Figure 3 (the abscissa is frequency in Hz, and the ordinate is normalized seismic amplitude, ranging from 0 to 1, dimensionless). With a 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 a normalized seismic amplitude of 0.5 as the boundary, the dominant frequency band range of the seismic data is obtained, namely the low-frequency dominant frequency Flow2 = 20 Hz, the high-frequency dominant frequency Fhigh2 = 60 Hz, and the mid-frequency dominant frequency range, i.e., the main seismic frequency (the frequency at the maximum normalized seismic amplitude), is Fmid = 40 Hz.

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

[0058] Specifically, conventional frequency division technology filters usually select a range of 10Hz around the dominant frequency for filtering, that is, a single-frequency body, as shown in Figure 4: 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. The bandwidth of the conventional filter is about 20Hz, which is relatively narrow overall. In order to overcome the defect of the narrow bandwidth of conventional frequency division filters, the present invention proposes a frequency division energy focusing filter, specifically: 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 shown in Figure 5, the low-cut, low-pass, high-pass, and high-cut frequencies of the mid-frequency energy focusing filter are [Flow1, Fmid, Fmid, Fhigh1], respectively. As shown in Figure 6, the low-cut, low-pass, high-pass, and high-cut frequencies of the high-frequency energy focusing filter are [Flow1, Fhigh2, Fhigh2, Fhigh1], respectively. As shown in Figure 7, the advantage of the frequency-splitting energy focusing filter is that the main frequency is focused on the low, mid, and high frequencies, achieving the purpose of frequency division. At the same time, each filter retains full-bandwidth information, overcoming the narrow bandwidth disadvantage of conventional frequency division technology.

[0059] The frequency band of the original earthquake is 10-70Hz, and the bandwidth is 60Hz. The bandwidth of the conventional frequency division filter is 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.

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

[0061] Specifically, the three frequency-division energy focusing filters designed by S2 are used to perform frequency band filtering on the seismic data respectively, and three frequency-division 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 rock reef, Figure 8 is the original earthquake number (Seis), and the triangle in the figure is the outline of the reef. Figure 9 is the low-frequency energy focusing seismic data (Seis_low) obtained by filtering the original seismic data using the low-frequency energy focusing filter in step 2, Figure 10 is the medium-frequency energy focusing seismic data (Seis_mid) obtained by filtering the original seismic data using the medium-frequency energy focusing filter of S2, and Figure 11 is the high-frequency energy focusing seismic data (Seis_high) obtained by filtering the original seismic data using the high-frequency energy focusing filter in step 2.

[0062] 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.

[0063] Specifically, the physical meaning of instantaneous phase is the angle between the phase vector and the real axis as a function of time. 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 is effective in depicting the lateral boundary characteristics of special carbonate reefs.

[0064] Complex seismic traces are used to calculate instantaneous phase attributes. First, the complex seismic traces are calculated, and then the instantaneous phase attributes are calculated from the complex seismic traces. The instantaneous phase attributes are instantaneous attributes derived from the complex seismic trace attributes. Complex seismic trace attributes are attributes obtained by performing a Hilbert transform on the original seismic traces. The seismic signal and its Hilbert transform are represented as f(t) and g(t), respectively. The complex seismic trace is represented as: F(t) = f(t) + ig(t) (Formula 1)

[0065] The expression of the instantaneous phase property is: Ph(t) = arctan[g(t) / f(t)] (Formula 2)

[0066] 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).

[0067] Figure 12 shows the low-frequency energy focusing instantaneous phase data (Phase_low) and the reef interpretation results. The right figure in Figure 12 is the reef interpretation result. The outer contour features of the reef body can be clearly identified from the figure. Due to the low resolution, the interpretation of the internal structure of the reef is not detailed enough.

[0068] Figure 13 shows the instantaneous phase data (Phase_mid) obtained using medium-frequency energy focusing and the interpretation results of the reefs. The right figure in Figure 13 shows the interpretation results of the reefs. As can be seen from the figure, the outline of the reefs is not as clear as that obtained using the instantaneous phase data obtained using low-frequency energy focusing, but the accuracy of the details of the internal structure of the reefs has been improved.

[0069] Figure 14 shows the high-frequency energy focusing instantaneous phase data (Phase_high) and the reef interpretation results. The right figure in Figure 14 shows the reef interpretation result. As can be seen from the figure, the outline of the reef becomes unclear, but the detailed information of the internal structure of the reef can be clearly portrayed.

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

[0071] This application prefers weighted average fusion. A reasonable data fusion strategy is key to obtaining high-quality fused images. Weighted average fusion is relatively simple to implement, not only has low system overhead, but also provides real-time results. It also suppresses noise in the original data. In many applications, this data fusion method often achieves satisfactory fusion results.

[0072] 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 the 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 weighted. In this example, all three weight coefficients are set to 0.33. After weighting, the 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 composite image. _mix =w A *Phase _low +w B *Phase _mid +w C *Phase _high (Formula 3)

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

[0074] As shown in Figure 15, the multi-scale instantaneous phase fusion data (Phase_mix) can clearly depict both the external contour features of the bioherm and the internal structure of the bioherm. It can also finely divide the bioherms into three phases, achieving the goal of simultaneously and clearly depicting the external contours and internal structural details of carbonate bioherms.

[0075] Example 3:

[0076] The present invention also proposes a device for obtaining bioherm information of geological layers, which includes: an acquisition module for obtaining the frequency range of the 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-focused seismic data of different frequency bands; a third processing module for respectively calculating energy-focused instantaneous phase data of energy-focused seismic data of different frequency bands; a fourth processing module for fusing the energy-focused 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.

[0077] Specifically, the energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands are calculated respectively, including: Ph(t) = arctan[g(t) / f(t)]

[0078] 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.

[0079] 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

[0080] Among them, Phase _mix To fuse data, Phase _low The energy focusing instantaneous phase data of the first frequency band, Phase _mid 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 coefficient of the energy-focused instantaneous phase data in the third frequency band. The energy-focused seismic data in different frequency bands includes low-frequency energy-focused seismic data, mid-frequency energy-focused seismic data, and high-frequency energy-focused seismic data. The bioherm information includes the external outline and internal structure of the bioherm. Band filtering the different frequency bands to obtain the energy-focused seismic data in the different frequency bands includes band filtering the different frequency bands using corresponding frequency-dividing energy-focusing filters.

[0081] The method of obtaining bioherm information of a geological layer of the present invention includes: obtaining a 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 according to the frequency range; performing frequency band filtering on different frequency bands to obtain energy-focused seismic data of different frequency bands; calculating energy-focused instantaneous phase data of the energy-focused seismic data of different frequency bands; fusing the energy-focused 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 frequency-divided energy-focused data volumes of different frequency bands using seismic data of different frequency bands, and then performs RGB fusion on the frequency-divided energy-focused instantaneous phase data volumes of different frequency bands to obtain multi-scale instantaneous phase fusion data, thereby finely characterizing the outline and internal structure details of carbonate bioherms.

[0082] 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.

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

[0084] 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.

[0085] 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.

[0086] 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.

[0087] An embodiment of the present invention provides a device comprising 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 performed: (method claim steps, sole claim + subordinate claim). The device herein may be a server, a PC, a PAD, a mobile phone, or the like.

[0088] 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 a frequency range of a 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-focused seismic data of different frequency bands; calculating energy-focused instantaneous phase data of the energy-focused seismic data of different frequency bands; fusing the energy-focused 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.

[0089] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can 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 can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0090] The present application is described with reference to the flow chart and / or block diagram of the method, device (system), and computer program product according to the embodiment of the present application. It should be understood that each flow process and / or box in the flow chart and / or block diagram and the combination of the flow process and / or box in the flow chart and / or block diagram can be realized 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 processing machine or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for realizing the function specified in one flow chart flow or multiple flows and / or one box or multiple boxes of the block diagram.

[0091] These computer program instructions may also be stored in a computer-readable memory that can direct 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 a product including an instruction device that implements the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

[0092] 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, so that the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

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

[0094] 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.

[0095] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. The 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 disc read-only memory (CD-ROM), digital versatile disc (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 transitory computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0096] It should also be noted that the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, commodity, or apparatus that includes a series of elements includes not only those elements but also other elements not explicitly listed, or includes elements inherent to such process, method, commodity, or apparatus. In the absence of further limitations, an element defined by the phrase "comprises a ..." does not exclude the presence of other identical elements in the process, method, commodity, or apparatus that includes the element.

[0097] The above are merely 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 modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present application should all 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; Determine the bioherm information of the geological layer to be measured by using the fused data; 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.

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: P(t)=arctan[g(t) / f(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: _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 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 20 Hz; The main frequency of the second frequency band is 40Hz; The main frequency of the third frequency band is 60 Hz.

5. A device for obtaining bioherm information of 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; A characterization module, used for determining the bioherm information of the geological layer to be measured by using the fused data; 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.

6. The device according to claim 5, characterized in that The energy focusing instantaneous phase data of energy focusing seismic data of different frequency bands are calculated respectively, including: P(t)=arctan[g(t) / f(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.

7. The device according to claim 5, 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.

8. The device according to claim 5, 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 second frequency band is 40Hz; The main frequency of the third frequency band is 60 Hz.

9. 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 4.

Images