A method and device for obtaining geological layer bio-reef information and a storage medium
By using frequency band filtering and phase data fusion, the problem of difficulty in depicting the external outline and internal structure of carbonate reefs was solved, and high-resolution reef information was obtained.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- CHINA NAT PETROLEUM CORP
- Filing Date
- 2023-09-20
- Publication Date
- 2026-05-01
AI Technical Summary
Existing technologies are insufficient to accurately characterize the external outline and internal structure of carbonate bioherms, and seismic data has low resolution, making it difficult for a single seismic attribute to simultaneously reflect its detailed features.
By obtaining the frequency range of the geological layer to be measured, it is divided into the first frequency band, the second frequency band and the third frequency band. Frequency band filtering is performed to obtain energy-focused seismic data. Energy-focused instantaneous phase data is calculated and fused. Multi-scale instantaneous phase fusion data is used to finely characterize bioherm information.
It enables the simultaneous and clear depiction of the external outline and internal structural details of carbonate reefs, improving the resolution and accuracy of seismic data.
Smart Images

Figure CN119667778B_ABST
Abstract
Description
A method, apparatus, and storage medium for obtaining information on bioherms in geological layers. Technical Field
[0001] This invention relates to the fields of oil and gas exploration and geophysics, and specifically to a method, apparatus and storage medium for obtaining information on geological bioherms. Background Technology
[0002] Carbonate bioherms are characterized by their favorable physical properties and large scale, making them ideal for hydrocarbon accumulation and important subjects in petroleum geophysical research. Due to their sedimentary formation, carbonate bioherms exhibit unique reflective structural features on seismic reflection profiles.
[0003] Existing methods for characterizing reefs include planar seismic facies prediction methods based on reflection feature analysis, such as profile geological body sculpting methods based on amplitude and frequency-based seismic attributes. However, due to the low resolution of seismic data, the morphology of carbonate reefs is difficult to accurately depict. Relying solely on a single seismic attribute can only reflect one aspect of the characteristics of carbonate reefs, making it difficult to simultaneously depict both the external outline and the specific details of the internal structure. Therefore, there is an urgent need for a descriptive technique capable of depicting both the mound-like external outline and the internal structure of carbonate reefs. Summary of the Invention
[0004] The purpose of this invention is to provide a method, apparatus, and storage medium for obtaining information about geological bioherms. This method enables the simultaneous depiction of the external outline and internal structural details of the bioherm.
[0005] To achieve the above objectives, embodiments of the present invention provide a method for obtaining information on geological layer bioherms, the method comprising:
[0006] Obtain the frequency range of the geological layer to be measured;
[0007] The frequency range is divided into three bands from low to high: the first band, the second band, and the third band.
[0008] Energy-focused seismic data for different frequency bands were obtained by performing band filtering on different frequency bands.
[0009] Calculate the instantaneous phase data of energy-focused seismic data in different frequency bands;
[0010] Fusion data is obtained by fusing the instantaneous phase data of energy focusing in the different frequency bands.
[0011] Bioherm information of the geological layer to be tested was determined using fused data.
[0012] Optionally, the step of calculating the instantaneous phase data of energy-focused seismic data in different frequency bands includes:
[0013] Ph(t) = arctan[f(t) / g(t)]
[0014] Where Ph(t) represents the instantaneous phase data of energy focusing.
[0015] f(t) is the seismic signal.
[0016] g(t) represents the Hilbert transform data of the seismic signal.
[0017] Optionally, fusing the instantaneous phase data of energy focusing in 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 integrate data,
[0020] Phase_ low This provides instantaneous phase data for energy focusing in the first frequency band.
[0021] Phase_ mid This is the instantaneous phase data for energy focusing in the second frequency band.
[0022] Phase_ high This is the instantaneous phase data for energy focusing in the third frequency band.
[0023] w A The weighting coefficients are the instantaneous phase data of energy focusing in the first frequency band.
[0024] w B The weighting coefficients are used for the instantaneous phase data of energy focusing in the second frequency band.
[0025] w C The weighting coefficients are used for the instantaneous phase data of energy focusing in the third frequency band.
[0026] Optionally, the energy-focused seismic data in different frequency bands includes: 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 bioherm's external outline and internal structure.
[0028] Optionally, the step of performing band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands includes:
[0029] Band filtering is performed on the first, second, and third frequency bands using appropriate frequency-division energy focusing filters;
[0030] The bandwidth of the frequency-division energy focusing filter is 60Hz;
[0031] The energy range of the frequency division energy focusing filter is 10-70Hz;
[0032] The energy of the frequency-division energy focusing filter is focused onto 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 40Hz;
[0035] The main frequency of the first frequency band is 60Hz.
[0036] The present invention also proposes a device for acquiring information on bioherms in geological layers, characterized in that the device comprises:
[0037] The acquisition module is used to acquire the frequency range of the geological layer to be measured;
[0038] The first processing module is used to divide 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 for different frequency bands;
[0040] The third processing module is used to calculate the instantaneous phase data of energy focusing seismic data in different frequency bands.
[0041] The fourth processing module is used to fuse the instantaneous phase data of energy focusing in the different frequency bands to obtain fused data;
[0042] The characterization module is used to determine bioherm information of the geological layer under test using fused data.
[0043] Optionally, the step of calculating the instantaneous phase data of energy-focused seismic data in different frequency bands includes:
[0044] Ph(t) = arctan[f(t) / g(t)]
[0045] Where Ph(t) represents the instantaneous phase data of energy focusing.
[0046] f(t) is the seismic signal.
[0047] g(t) represents the Hilbert transform data of the seismic signal.
[0048] Optionally, fusing the instantaneous phase data of energy focusing in 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 integrate data,
[0051] Phase_ low This provides instantaneous phase data for energy focusing in the first frequency band.
[0052] Phase_ mid This is the instantaneous phase data for energy focusing in the second frequency band.
[0053] Phase_ high This is the instantaneous phase data for energy focusing in the third frequency band.
[0054] w A The weighting coefficients are the instantaneous phase data of energy focusing in the first frequency band.
[0055] w B The weighting coefficients are used for the instantaneous phase data of energy focusing in the second frequency band.
[0056] w C The weighting coefficients are used for the instantaneous phase data of energy focusing in the third frequency band.
[0057] Optionally, the energy-focused seismic data in different frequency bands includes: 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 bioherm's external outline and internal structure.
[0059] Optionally, the step of performing band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands includes:
[0060] Band filtering is performed on the first, second, and third frequency bands using appropriate frequency-division energy focusing filters;
[0061] The bandwidth of the frequency-division energy focusing filter is 60Hz;
[0062] The energy range of the frequency division energy focusing filter is 10-70Hz;
[0063] The energy of the frequency-division energy focusing filter is focused onto 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 40Hz;
[0066] The main frequency of the first frequency band is 60Hz.
[0067] The present invention also proposes a machine-readable storage medium storing instructions for causing a machine to perform the above-described method for obtaining bioherm characterization.
[0068] The method for obtaining bioherm information of geological layers according to the present invention includes: 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 band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands; calculating the instantaneous phase data of the energy-focused seismic data for different frequency bands; fusing the instantaneous phase data of the energy-focused seismic data for different frequency bands to obtain fused data; and using the fused data to determine the bioherm information of the geological layer to be measured. The present invention determines the frequency-divided energy-focused data volumes for different frequency bands using seismic data from different frequency bands, and then performs RGB fusion on the instantaneous phase data volumes of the frequency-divided energy-focused seismic data for different frequency bands to obtain multi-scale instantaneous phase fused data, thereby finely characterizing the external outline and internal structural details of carbonate bioherms.
[0069] Other features and advantages of the embodiments of the present invention will be described in detail in the following detailed description section. Attached Figure Description
[0070] The accompanying drawings are provided to further illustrate embodiments of the present invention and form part of the specification. They are used together with the following detailed description to explain the embodiments of the present invention, but do not constitute a limitation thereof. In the drawings:
[0071] Figure 1 is a flowchart illustrating a method for obtaining bioherm information in geological strata according to the present invention;
[0072] Figure 2 is a schematic diagram of an embodiment of a method for obtaining bioherm information of geological layers according to the present invention;
[0073] Figure 3 is a schematic diagram of the spectrum obtained by performing a spectral scan on the target layer of the seismic data.
[0074] Figure 4 is a schematic diagram of the frequency band characteristics of the low-frequency, mid-frequency and high-frequency sections of the frequency divider filter;
[0075] Figure 5 is a schematic diagram of a low-frequency energy focusing filter focusing the main frequency to the low-frequency dominant frequency.
[0076] Figure 6 is a schematic diagram of the intermediate frequency energy focusing filter focusing the main frequency to the intermediate frequency dominant frequency;
[0077] Figure 7 is a schematic diagram of the high-frequency energy focusing filter focusing the main frequency to the high-frequency dominant frequency;
[0078] Figure 8 is a schematic diagram of the original seismic data;
[0079] Figure 9 is a schematic diagram of low-frequency energy focusing seismic data;
[0080] Figure 10 is a schematic diagram of mid-frequency energy-accumulating seismic data;
[0081] Figure 11 is a schematic diagram of high-frequency energy focusing seismic data;
[0082] Figure 12 is a schematic diagram of the instantaneous phase data of low-frequency energy focusing and the interpretation results of bioherms;
[0083] Figure 13 is a schematic diagram of the instantaneous phase data of mid-frequency energy focusing and the interpretation results of bioherms;
[0084] Figure 14 is a schematic diagram of the instantaneous phase data of high-frequency energy focusing and the interpretation results of bioherms;
[0085] Figure 15 is a schematic diagram of multi-scale instantaneous phase fusion data and bioherm interpretation results. Detailed Implementation
[0086] The specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are for illustration and explanation only and are not intended to limit the scope of the present invention.
[0087] Example 1:
[0088] Figure 1 is a flowchart illustrating a method for obtaining geological layer bioherm information according to the present invention. As shown in Figure 1, the method for obtaining geological layer bioherm information according to the present invention includes: 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 of seismic data. The target layer of seismic data is subjected to spectral scanning to determine the frequency range.
[0089] Step S102 involves dividing 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 will be referred to as the low-frequency band, the second frequency band as the mid-frequency band, and the third frequency band 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.
[0090] The frequency-division energy focusing filter is divided into low-frequency, mid-frequency, and high-frequency sections, maintaining a bandwidth of 10-70Hz (i.e., 60Hz). Energy is focused to low-frequency (main frequency 20Hz), mid-frequency (main frequency 40Hz), and high-frequency (main frequency 60Hz). Specifically, the low-frequency energy range is 10-70Hz, with a bandwidth of 60Hz and a main frequency of 20Hz (as shown in Figure 5); the mid-frequency energy range is 10-70Hz, with a bandwidth of 60Hz and a main frequency of 40Hz (as shown in Figure 6); and the high-frequency energy range is 10-70Hz, with a bandwidth of 60Hz and a main frequency of 60Hz (as shown in Figure 7). The frequency-division energy focusing filter has a bandwidth of 60Hz, which not only achieves the purpose of frequency division but also ensures that each filter retains 60Hz of full-band information.
[0091] Step S103 involves performing band filtering on different frequency bands to obtain energy-focused seismic data for each band. Specifically, the energy-focused seismic data for 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 bioherm's external outline and internal structure. The step of performing band filtering on different frequency bands to obtain energy-focused seismic data for each band includes: performing band filtering on the different frequency bands using corresponding frequency-division energy-focusing filters.
[0092] According to a specific implementation method, three frequency-division energy focusing filters are used to perform frequency band filtering on the seismic data to obtain first-band energy focusing seismic data, second-band energy focusing seismic data, and third-band energy focusing seismic data, namely low-frequency energy focusing seismic data, mid-frequency energy focusing seismic data, and high-frequency energy focusing seismic data.
[0093] Step S104 involves calculating the instantaneous phase data of energy-focused seismic data for different frequency bands. Specifically, the instantaneous phase data includes low-frequency, mid-frequency, and high-frequency energy-focused seismic data. The data volumes obtained by calculating using the following formulas are the instantaneous phase data for low-frequency, mid-frequency, and high-frequency energy-focused seismic data.
[0094] Specifically, the calculation of instantaneous phase data of energy-focused seismic data in different frequency bands includes:
[0095] Ph(t) = arctan[f(t) / g(t)]
[0096] Where Ph(t) is the instantaneous phase data of energy focusing, f(t) is the seismic signal, and g(t) is the Hilbert transform data of the seismic signal.
[0097] Step S105 involves fusing the instantaneous phase data of energy focusing from different frequency bands to obtain fused data. Specifically, the fusion of the instantaneous phase data of energy focusing from different frequency bands yields fused data. The preferred fusion method is weighted average fusion, including:
[0098] Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high in,
[0099] Phase_ mix To integrate the data, Phase_ low This is the instantaneous phase data for energy focusing in the first frequency band (i.e., the low frequency band), Phase_ mid This is the instantaneous phase data for energy focusing in the second frequency band (i.e., the mid-frequency band), Phase_ high For the energy focusing instantaneous phase data of the third frequency band (i.e., the high frequency band), w A The weighting coefficients for the instantaneous phase data of energy focusing in the first frequency band (i.e., low frequency) are w. B The weighting coefficients for the instantaneous phase data of energy focusing in the second frequency band (i.e., the mid-frequency band), w C These are the weighting coefficients for the instantaneous phase data of energy focusing in the third frequency band (i.e., high frequency). The fused data is multi-scale instantaneous phase fusion data.
[0100] Step S106 involves using fused data to determine the bioherm information of the geological layer to be measured. Specifically, the bioherm information includes the external outline and internal structure of the bioherm. Low-frequency energy focusing instantaneous phase data volume reflects the external outline of the carbonate bioherm, while high-frequency energy focusing instantaneous phase data volume reflects the details of the internal structure of the carbonate bioherm.
[0101] This invention proposes a method for acquiring information about bioherms in geological layers. Based on determining the effective frequency range of the target seismic data segment, three frequency-division energy focusing filters (low-frequency, mid-frequency, and high-frequency) are designed. These filters retain full-band information while focusing the main frequency to low, mid, and high frequencies respectively, overcoming the limitations of conventional frequency-division techniques (single-frequency volumes) which suffer from narrow bandwidth and unclear geological significance. Furthermore, instantaneous phase attributes are calculated based on the frequency-division energy focusing, enabling the low-frequency energy-focused instantaneous phase data volume to reflect the external contour of the carbonate bioherm, and the high-frequency energy-focused instantaneous phase data volume to reflect the internal structural details of the carbonate bioherm. High-resolution image processing is achieved through RGB fusion, effectively identifying image details and highlighting anomalies. Weighted average fusion yields multi-scale instantaneous phase fusion data, achieving the effect of simultaneously and clearly depicting both the external contour and internal structural details of the carbonate bioherm.
[0102] Example 2:
[0103] Figure 2 is a schematic diagram of an embodiment of a method for obtaining geological layer bioherm information according to the present invention. As shown in Figure 2, step S1 is to perform a spectral scan on the target layer of seismic data to determine the frequency range.
[0104] Specifically, a spectral scan is performed on the target layer of the seismic data volume. This target layer is used to study the top and bottom boundaries of the target. Typically, the scanned layer is no less than two wavelet lengths, i.e., no less than 200 ms, to ensure the stability of the scanned spectrum. When the target layer thickness is small, it can be appropriately extended upwards and downwards towards the top and bottom layers. The resulting spectrum is shown in Figure 3 (the horizontal axis represents frequency in Hz, and the vertical axis represents normalized seismic amplitude, ranging from 0 to 1, dimensionless). With a normalized seismic amplitude of 0.25 as the boundary, the effective frequency band is 10–70 Hz, i.e., 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 range of the seismic data is obtained, i.e., 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 seismic main frequency (the frequency at the maximum normalized seismic amplitude), is Fmid = 40 Hz.
[0105] Step S2 involves designing frequency-division energy focusing filters for three different frequency bands: low frequency, mid frequency, and high frequency.
[0106] Specifically, conventional frequency-divergent filters typically select a 10Hz range around the dominant frequency for filtering, i.e., a single-frequency filter, as shown in Figure 4: the low-frequency energy focusing filter (Filter_low) has low-cutoff, low-pass, high-pass, and high-cutoff frequencies of [Flow2-10Hz, Flow2, Flow2, Flow2+10Hz]; the mid-frequency energy focusing filter (Filter_mid) has low-cutoff, low-pass, high-pass, and high-cutoff frequencies of [Fmid-10Hz, Fmid, Fmid, Fmid+10Hz]; and the high-frequency energy focusing filter (Filter_high) has low-cutoff, low-pass, high-pass, and high-cutoff frequencies of [Fhigh2-10Hz, Fhigh2, Fhigh2, Fhigh2+10Hz]. The bandwidth of conventional filters is typically 20Hz, which is relatively narrow. To overcome the narrow bandwidth of conventional frequency-divergent filters, this invention proposes a frequency-divergent energy focusing filter, specifically: the low-frequency energy focusing filter has low-cutoff, low-pass, high-pass, and high-cutoff frequencies of [Flow1, Flow2, Flow2, Fhigh1]. As shown in Figure 5, the low-pass, low-pass, high-pass, and high-cutoff frequencies of the intermediate frequency energy focusing filter are [Flow1, Fmid, Fmid, Fhigh1]. As shown in Figure 6, the low-pass, low-pass, high-pass, and high-cutoff frequencies of the high-frequency energy focusing filter are [Flow1, Fhigh2, Fhigh2, Fhigh1]. As shown in Figure 7, the advantage of the frequency-division energy focusing filter is that the main frequency is focused onto the low, intermediate, and high frequencies respectively, achieving the purpose of frequency division. At the same time, each filter retains full-band information, overcoming the narrow bandwidth defect of conventional frequency division technology.
[0107] The original earthquake's frequency range is 10-70Hz, with a bandwidth of 60Hz. Conventional frequency-division filters only have a bandwidth of 20Hz, which is too narrow and lacks clear geological significance. The frequency-division energy focusing filter proposed in this invention also divides into low-frequency, mid-frequency, and high-frequency sections, but maintains a bandwidth (10-70Hz, i.e., 60Hz), focusing energy to low-frequency (dominant frequency 20Hz), mid-frequency (dominant frequency 40Hz), and high-frequency (dominant frequency 60Hz). Specifically, the low-frequency energy range is 10-70Hz, with a bandwidth of 60Hz and a dominant frequency of 20Hz; the mid-frequency energy range is 10-70Hz, with a bandwidth of 60Hz and a dominant frequency of 40Hz; and the high-frequency energy range is 10-70Hz, with a bandwidth of 60Hz and a dominant frequency of 60Hz. The 60Hz bandwidth of this invention not only achieves the purpose of frequency division but also ensures that each filter retains 60Hz of full-band information.
[0108] Step S3 involves using three frequency-division energy focusing filters to perform frequency band filtering on the seismic data.
[0109] Specifically, the seismic data is band-filtered using three frequency-division energy focusing filters designed in S2, resulting in three frequency-division energy focusing data volumes: low-frequency energy focusing seismic data (Seis_low), mid-frequency energy focusing seismic data (Seis_mid), and high-frequency energy focusing seismic data (Seis_high). This embodiment uses a carbonate reef. Figure 8 shows the original seismic data (Seis), with the triangles representing the reef's outline. Figure 9 shows 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 two. Figure 10 shows the mid-frequency energy focusing seismic data (Seis_mid) obtained by filtering the original seismic data using the mid-frequency energy focusing filter in S2. Figure 11 shows 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 two.
[0110] Step S4 involves calculating the instantaneous phase attributes of the three frequency-division energy focusing seismic data volumes to obtain the three frequency-division energy focusing instantaneous phase data volumes.
[0111] Specifically, the physical meaning of instantaneous phase is the angle between the phase vector, which is a function of time, and the real axis. Instantaneous phase attributes are independent of seismic energy, can enhance weak reflection signals, highlight the lateral variation characteristics of the same phase axis, and have a good effect on characterizing the lateral boundary features of special carbonate reefs.
[0112] Complex seismic traces are used to calculate instantaneous phase attributes. First, the complex seismic trace is calculated, and then the instantaneous phase attributes are calculated from it. The instantaneous phase attributes are a type of instantaneous attribute derived from the attributes of the complex seismic trace, which are obtained by performing a Hilbert transform on the original seismic trace. If the seismic signal and its Hilbert transform are represented as f(t) and g(t) respectively, then the complex seismic trace is represented as:
[0113] F(t)=f(t)+ig(t) (Formula 1) The expression for the instantaneous phase attribute is:
[0114] Ph(t) = arctan[f(t) / g(t)] (Formula 2)
[0115] The instantaneous phase attributes are calculated for the three frequency-division energy focusing data volumes respectively, resulting in three frequency-division energy focusing instantaneous phase data, namely, low-frequency energy focusing instantaneous phase data (Phase_low), mid-frequency energy focusing instantaneous phase data (Phase_mid), and high-frequency energy focusing instantaneous phase data (Phase_high).
[0116] Figure 12 shows the instantaneous phase data (Phase_low) of low-frequency energy focusing and the interpretation results of the bioherm. The right figure in Figure 12 shows the interpretation results of the bioherm. The outline features of the bioherm can be clearly identified from the figure. Due to the low resolution, the interpretation of the internal structural details of the bioherm is not detailed enough.
[0117] Figure 13 shows the instantaneous phase data (Phase_mid) of mid-frequency energy focusing and the interpretation results of the bioherm. The right figure in Figure 13 shows the interpretation results of the bioherm. It can be seen from the figure that the outline of the bioherm is not as clear as the instantaneous phase data of low-frequency energy focusing, but the accuracy of the details of the internal structure of the bioherm is improved.
[0118] Figure 14 shows the instantaneous phase data (Phase_high) of high-frequency energy focusing and the interpretation results of the bioherm. The right figure in Figure 14 shows the interpretation results of the bioherm. It can be seen from the figure that the outline of the bioherm becomes unclear, but the internal structural details of the bioherm can be clearly depicted.
[0119] Step S5 involves weighted averaging and fusing the three frequency-division energy focusing instantaneous phase data volumes to obtain multi-scale instantaneous phase fusion data.
[0120] This application preferentially uses weighted average fusion, as a reasonable data fusion strategy is key to obtaining high-quality fused images. Weighted average fusion is relatively simple to implement, has low system overhead, and can obtain results in real time. Furthermore, this method can suppress noise in the original data. In many applications, this data fusion method often achieves satisfactory fusion results.
[0121] In this invention, instantaneous phase data of energy focusing at different frequencies are obtained through frequency division energy focusing. The instantaneous phase data of energy focusing at different frequencies are then fused using a weighted average fusion method (Formula 3), where w... A w B w C The weighting coefficients for the instantaneous phase data of energy focusing at different frequencies are used; in this example, all three weighting coefficients are set to 0.33. After weighting, multi-scale instantaneous phase fusion data (Phase_mix) is obtained, thus preserving the significant features of the original image in different frequency domains in the synthesized image.
[0122] Phase_ mix =w A *Phase_ low +w B *Phase_ mid +w C *Phase_ high (Formula 3)
[0123] Step S6 involves using multi-scale instantaneous phase fusion data to characterize the external outline and internal structural details of carbonate bioherms.
[0124] As shown in Figure 15, multi-scale instantaneous phase fusion data (Phase_mix) can clearly depict both the external outline and internal structure of bioherms, and can finely divide the three phases of bioherms, thus achieving the goal of simultaneously and clearly depicting both the external outline and internal structural details of carbonate bioherms.
[0125] Example 3:
[0126] This invention also proposes a device for acquiring bioherm information of geological layers. The device includes: an acquisition module for acquiring 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 band filtering on the different frequency bands to obtain energy-focused seismic data for each band; a third processing module for calculating the instantaneous phase data of the energy-focused seismic data for each band; a fourth processing module for fusing the instantaneous phase data of the energy-focused seismic data for each band 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.
[0127] Specifically, the calculation of instantaneous phase data of energy-focused seismic data in different frequency bands includes:
[0128] Ph(t) = arctan[f(t) / g(t)]
[0129] Where Ph(t) is the instantaneous phase data of energy focusing, f(t) is the seismic signal, and g(t) is the Hilbert transform data of the seismic signal.
[0130] The process of fusing the instantaneous phase data of energy focusing in 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 integrate the data, Phase_ low For the energy focusing instantaneous phase data of the first frequency band, Phase_ mid For the energy focusing instantaneous phase data of the second frequency band, Phase_high For the energy focusing instantaneous phase data of the third frequency band, w A w represents the weighting coefficients for the instantaneous phase data of energy focusing in the first frequency band. B w represents the weighting coefficients for the instantaneous phase data of energy focusing in the second frequency band. C The weighting coefficients are for the instantaneous phase data of energy focusing in the third frequency band. The energy focusing seismic data in different frequency bands include: low-frequency energy focusing seismic data, mid-frequency energy focusing seismic data, and high-frequency energy focusing seismic data; the bioherm information includes the bioherm's external outline and internal structure. The step of obtaining energy focusing seismic data in different frequency bands by performing band filtering on each band includes: performing band filtering on each frequency band using corresponding frequency-division energy focusing filters.
[0133] The method for obtaining bioherm information of geological layers according to the present invention includes: 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 band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands; calculating the instantaneous phase data of the energy-focused seismic data for different frequency bands; fusing the instantaneous phase data of the energy-focused seismic data for different frequency bands to obtain fused data; and using the fused data to determine the bioherm information of the geological layer to be measured. The present invention determines the frequency-divided energy-focused data volumes for different frequency bands using seismic data from different frequency bands, and then performs RGB fusion on the instantaneous phase data volumes of the frequency-divided energy-focused seismic data for different frequency bands to obtain multi-scale instantaneous phase fused data, thereby finely characterizing the external outline and internal structural details of carbonate bioherms.
[0134] On the other hand, the present invention also provides an electronic device comprising: at least one processor; and a memory connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, and the at least one processor implements the method for obtaining geological layer bioherm information described above by executing the instructions stored in the memory.
[0135] The processor contains a kernel that retrieves corresponding program units from memory. One or more kernels can be configured, and by adjusting kernel parameters, detailed information characterizing the external outline and internal structure of the bioherm can be achieved.
[0136] The memory may include non-permanent memory in computer-readable media, such as 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] This invention provides a storage medium storing a program that, when executed by a processor, implements the oil extraction method.
[0138] This invention provides a processor for running a program, wherein the program executes the method for obtaining geological layer bioherm information.
[0139] This invention provides a device including a processor, a memory, and a program stored in the memory and executable on the processor. When the processor executes the program, it performs the following steps: (method claim steps, exclusive claim + subordinate claim). The device described herein can be a server, PC, PAD, mobile phone, etc.
[0140] This 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 band filtering on the different frequency bands to obtain energy-focused seismic data for the different frequency bands; calculating the instantaneous phase data of the energy-focused seismic data for the different frequency bands; fusing the instantaneous phase data of the energy-focused seismic data for the different frequency bands to obtain fused data; and using the fused data to determine the bioherm information of the geological layer to be measured.
[0141] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0142] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in one or more flowchart illustrations and / or one or more block diagrams.
[0143] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means that implement the functions specified in one or more flowcharts and / or one or more block diagrams.
[0144] These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process, such that the instructions, which execute on the computer or other programmable apparatus, provide steps for implementing the functions specified in one or more flowcharts and / or one or more block diagrams.
[0145] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.
[0146] Memory may include non-persistent memory in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.
[0147] Computer-readable media includes both permanent and non-permanent, removable and non-removable media that can store information using any method or technology. Information can be computer-readable instructions, data structures, modules of programs, 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 technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transferable medium that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transient computer-readable media, such as modulated data signals and carrier waves.
[0148] It should also be noted that the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus. Unless otherwise specified, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes that element.
[0149] The above are merely embodiments of this application and are not intended to limit the scope of this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of the claims of this application.
Claims
1. A method for obtaining information on bioherms in geological layers, characterized in that, The method includes: acquiring 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 band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands; calculating the instantaneous phase data of energy-focused seismic data for different frequency bands; fusing the instantaneous phase data of energy-focused seismic data for different frequency bands to obtain fused data; and using the fused data to determine the bioherm information of the geological layer to be measured. The step of performing band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands includes: performing band filtering on the first, second, and third frequency bands respectively using 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; and the energy of the frequency-division energy-focusing filter is focused onto the first, second, and third frequency bands. The step of calculating the instantaneous phase data of energy-focused seismic data for different frequency bands includes: in, For energy focusing instantaneous phase data, For earthquake signals, The data consists of Hilbert transform data of seismic signals; the fusion of energy-focused instantaneous phase data from different frequency bands to obtain fused data includes: Among them, Phase_ mix To integrate the data, Phase_ low For the energy focusing instantaneous phase data of the first frequency band, Phase_ mid For the energy focusing instantaneous phase data of the second frequency band, Phase_ high This is the instantaneous phase data for energy focusing in the third frequency band. The weighting coefficients are the instantaneous phase data of energy focusing in the first frequency band. The weighting coefficients are used for the instantaneous phase data of energy focusing in the second frequency band. The weighting coefficients are used for the instantaneous phase data of energy focusing in the third frequency band.
2. The method according to claim 1, characterized in that, The energy-focused seismic data in different frequency bands includes: energy-focused seismic data in the first frequency band, energy-focused seismic data in the second frequency band, and energy-focused seismic data in the third frequency band; the bioherm information includes the external outline and internal structure of the bioherm.
3. The method according to claim 1, characterized in that, The main frequency of the first frequency band is 20Hz; the main frequency of the second frequency band is 40Hz; and the main frequency of the third frequency band is 60Hz.
4. A device for acquiring information on bioherms in geological layers, characterized in that, The device includes: an acquisition module for acquiring 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 band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands; a third processing module for calculating the instantaneous phase data of energy-focused seismic data for different frequency bands; a fourth processing module for fusing the instantaneous phase data of energy-focused seismic data for 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. The step of performing band filtering on different frequency bands to obtain energy-focused seismic data for different frequency bands includes: performing band filtering on the first, second, and third frequency bands respectively using 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; and the energy of the frequency-division energy-focusing filter is focused onto the first, second, and third frequency bands. The step of calculating the instantaneous phase data of energy-focused seismic data for different frequency bands includes: in, For energy focusing instantaneous phase data, For earthquake signals, The data consists of Hilbert transform data of seismic signals; the fusion of energy-focused instantaneous phase data from different frequency bands to obtain fused data includes: Among them, Phase_ mix To integrate the data, Phase_ low For the energy focusing instantaneous phase data of the first frequency band, Phase_ mid For the energy focusing instantaneous phase data of the second frequency band, Phase_ high This is the instantaneous phase data for energy focusing in the third frequency band. The weighting coefficients are the instantaneous phase data of energy focusing in the first frequency band. The weighting coefficients are used for the instantaneous phase data of energy focusing in the second frequency band. The weighting coefficients are used for the instantaneous phase data of energy focusing in the third frequency band.
5. The apparatus according to claim 4, characterized in that, The energy-focused seismic data in different frequency bands includes: energy-focused seismic data in the first frequency band, energy-focused seismic data in the second frequency band, and energy-focused seismic data in the third frequency band; the bioherm information includes the external outline and internal structure of the bioherm.
6. The apparatus according to claim 4, characterized in that, The main frequency of the first frequency band is 20Hz; the main frequency of the second frequency band is 40Hz; and the main frequency of the third frequency band is 60Hz.
7. A machine-readable storage medium having instructions stored thereon for causing a machine to perform the method for obtaining geological layer bioherm information as described in any one of claims 1 to 3.
Citation Information
Patent Citations
Reef reservoir engraving method
CN104977611A
Method for predicting reservoir stratum by utilizing seismic energy difference attribute
CN104977612A