CN-116009094-B - Ultra-deep carbonate rock dome beach body describing method based on frequency division attribute

CN116009094BCN 116009094 BCN116009094 BCN 116009094BCN-116009094-B

Abstract

The invention discloses an ultra-deep carbonate beach body characterization method based on frequency division attributes, which is comprehensively applied to innovation in the prior art, utilizes the constraint of a geological prediction model to reduce the multi-resolution of a seismic prediction means by fully combining the geological model and the seismic prediction method, solves the limitation of single geology or the seismic prediction method, improves the beach plane distribution prediction precision, can be better applied to the exploration and development of a beach reservoir, and is characterized by focusing on the conventional attributes in a full-frequency data body, wherein the predicted beach body precision is limited, and the boundary range of the beach body can be finely characterized based on the frequency division attribute of high-frequency and low-frequency information of seismic data, and meanwhile, the depiction precision of beach body details can be effectively improved.

Inventors

ZHANG KUN
TAN LEI
YU TONG
WANG DONG
LIU HONG
LIANG FENG
HE BING
WANG LIEN
LIU QIANYU
Hu Luojia
YANG YING
YANG MENGXIANG

Assignees

西南石油大学

Dates

Publication Date: 20260508
Application Date: 20230104

Claims (7)

1. An ultra-deep carbonate hill body characterization method based on frequency division attribute, which is characterized by comprising the following steps: s1, carrying out single-well sequence division by combining drilling logging information loaded into an earthquake interpretation work area, rock core and sheet information loaded into the earthquake interpretation work area and field outcrop information, and establishing a sequence stratum grid; s2, well logging data are subjected to well earthquake fine calibration by using stratum layering data after synthetic earthquake records are calculated by using sound waves and density curves, a time-depth relation is obtained, and after a regional earthquake grid section is established by using the time-depth relation, a full-region three-dimensional earthquake horizon is accurately tracked and explained; S3, determining the scale level, the petrophysical parameters and the longitudinal and transverse distribution rules of the hilly beach body of the research area based on drilling logging data, field outcrop data, coring data and the established layer sequence stratum grillwork, and establishing a hilly beach body geological prediction model; s4, based on the established geological prediction model of the beach body, performing forward modeling of the earthquake to obtain a forward modeling result of the earthquake; s5, performing seismic frequency division processing on the forward modeling result of the earthquake, determining seismic reflection characteristics of the hilly beach bodies of different scales under different frequency bands, establishing a seismic response mode of the hilly beach body of the research area based on the seismic reflection characteristics, and matching optimal identification frequencies corresponding to the hilly beach bodies of different scales based on the seismic response mode; S6, performing seismic frequency division processing on the three-dimensional post-stack seismic data to obtain a plurality of single-frequency-band data volumes which are in reference to the optimal identification frequency obtained in the step S5, wherein the main frequency and the effective bandwidth of the seismic data in the research area are obtained; s7, based on the full-area three-dimensional seismic horizon, carrying out attribute extraction and algorithm optimization on each single-frequency band data body with different frequencies to obtain an optimized frequency division attribute and layer sequence stratum thickness; s8, after matching statistical correlation is carried out on the conditions of the real drilling drill meeting the hilly beach and the data related to the development of the hilly beach and the optimized frequency division attribute and the layer sequence stratum thickness, a conversion relation between the frequency division attribute and the hilly beach is established; S9, eliminating abnormal values in the frequency division attribute, and fusing the frequency division attribute to describe the plane spread of the hillside body of the research area based on the conversion relation between the frequency division attribute and the hillside body established in the step S8.
2. The method for characterizing an ultra-deep carbonate dome beach based on frequency division properties of claim 1, further comprising, prior to step S1: The method comprises the steps of obtaining basic data in a research area, wherein the basic data comprise field outcrop data, coring data, drilling logging data, stratum layering data and three-dimensional post-stack seismic data; and loading the basic data in the research area into the earthquake interpretation work area, and performing spectrum analysis on the three-dimensional post-stack earthquake data loaded into the earthquake interpretation work area to obtain the main frequency and the effective bandwidth of the earthquake data in the research area.
3. The method for describing the ultra-deep carbonate hill beach based on the frequency division attribute according to claim 1, wherein in the step S4, the excitation parameters and the three-dimensional seismic data volume parameters are ensured to be consistent in forward modeling, the method comprises main frequency and effective bandwidth information, and meanwhile, bandpass wavelets are required to be used in forward modeling, the wavelet main frequency and bandwidth are kept consistent with the seismic data volume parameters, a forward modeling result, namely a modeling section, is required to be compared with an over-real well seismic section, forward modeling effects are analyzed, and model parameters are optimized and adjusted.
4. The method of claim 1, wherein in steps S5 and S6, the seismic frequency division method uses a frequency division technique based on a matching pursuit algorithm, which decomposes any signal into a linear expansion of a set of basis functions as a technical means for processing highly unstable signals, iteratively searches the constructed ultra-complete wavelet library D for wavelet signals that most approximate residual signals, and gradually decomposes the original seismic signals into optimal linear superposition of a plurality of matching wavelets at a selected frequency As the center frequency, the frequency band range is used Reconstructing the seismic data to obtain a frequency division seismic data volume.
5. The method for describing the ultra-deep carbonate hill body based on the frequency division attribute according to claim 1 is characterized in that in the step S7, a low-frequency body is used for describing the phase band boundary of a large-scale hill body, a high-frequency body is used for carving small-scale hill bodies and internal details, root mean square amplitude attribute is used for describing the phase band boundary of the large-scale hill body, and waveform clustering attribute is used for carving the small-scale hill bodies and the internal details.
6. The method for describing the ultra-deep carbonate beach body based on the frequency division attribute according to claim 1 is characterized in that in the step S8, the sediment phase type of the beach body is determined to belong to an underdetermined problem according to the seismic attribute, constraint conditions are added in practical application, real drilling data are used as constraint conditions, the corresponding attribute types are calibrated by using geological research results and logging phases corresponding to the geological research results of the coring well of the beach body, and the geological research results are used for establishing the corresponding relationship between the geological research results and the logging phases and are recursively applied to a non-well area, so that the seismic phases are converted into sediment phases, and the purpose of finely describing the plane spread characteristics of the beach body in a region Neiqiu is achieved.
7. The method for describing the ultra-deep carbonate hill body based on the frequency division attribute according to claim 1, wherein in the step S9, the attribute abnormal value caused by the data reason is removed, the data quality is poor aiming at the three-dimensional data non-full coverage area or the broken belt, the extracted attribute value tends to have an abnormal point, the abnormal value is removed to be mapped by combining the actual situation, and finally the fine engraving of the hill body is completed.

Description

Ultra-deep carbonate rock dome beach body describing method based on frequency division attribute Technical Field The invention belongs to the field of petroleum exploration and development, and particularly relates to a method for describing ultra-deep carbonate hillock based on frequency division attribute. Background Carbonate rock building formed by microbial rocks formed by capturing and binding debris sediments by bacteria and algae microorganisms and particle rocks related to the microbial rocks are called a microbial hill complex (called hill body for short), are taken as important novel oil and gas resource carriers, also show great exploration potential in global oil and gas exploration practice, and become a hot spot of current research due to special sediment characteristics and storage processes. The existing carbonate rock beach earthquake prediction method is mainly based on conventional earthquake data, adopts methods such as earthquake reflection geometric structures, earthquake attributes, earthquake inversion and the like, but the methods are relatively dependent on the quality of the earthquake data, and meanwhile, due to the fact that the carbonate rock beach is generally buried deeply and poor in quality of the earthquake data and lacks proper geological model guidance, the reliability of explaining the distribution of the beach is low by utilizing the earthquake data, and the accuracy of the distribution prediction result of the beach is restricted. Therefore, the ultra-deep carbonate hill beach body describing method based on the frequency division attribute is invented, and the technical problems are solved. Disclosure of Invention The invention aims to solve the problems in the background art, and aims to provide an ultra-deep carbonate hill beach body describing method based on frequency division attribute, which improves the problems in the prior art. In order to solve the technical problems, the technical scheme of the invention is as follows: an ultra-deep carbonate hill body characterization method based on frequency division properties, the method comprising: s1, carrying out single-well sequence division by combining drilling logging information loaded into an earthquake interpretation work area, rock core and sheet information loaded into the earthquake interpretation work area and field outcrop information, and establishing a sequence stratum grid; s2, well logging data are subjected to well earthquake fine calibration by using stratum layering data after synthetic earthquake records are calculated by using sound waves and density curves, a time-depth relation is obtained, and after a regional earthquake grid section is established by using the time-depth relation, a full-region three-dimensional earthquake horizon is accurately tracked and explained; S3, determining the scale level, the petrophysical parameters and the longitudinal and transverse distribution rules of the hilly beach body of the research area based on drilling logging data, field outcrop data, coring data and the established layer sequence stratum grillwork, and establishing a hilly beach body geological prediction model; s4, based on the established geological prediction model of the beach body, performing forward modeling of the earthquake to obtain a forward modeling result of the earthquake; s5, performing seismic frequency division processing on the forward modeling result of the earthquake, determining seismic reflection characteristics of the hilly beach bodies of different scales under different frequency bands, establishing a seismic response mode of the hilly beach body of the research area based on the seismic reflection characteristics, and matching optimal identification frequencies corresponding to the hilly beach bodies of different scales based on the seismic response mode; S6, performing seismic frequency division processing on the three-dimensional post-stack seismic data to obtain a plurality of single-frequency-band data volumes which are in reference to the optimal identification frequency obtained in the step S5, wherein the main frequency and the effective bandwidth of the seismic data in the research area are obtained; s7, based on the full-area three-dimensional seismic horizon, carrying out attribute extraction and algorithm optimization on each single-frequency band data body with different frequencies to obtain an optimized frequency division attribute and layer sequence stratum thickness; s8, after matching statistical correlation is carried out on the conditions of the real drilling drill meeting the hilly beach and the data related to the development of the hilly beach and the optimized frequency division attribute and the layer sequence stratum thickness, a conversion relation between the frequency division attribute and the hilly beach is established; S9, eliminating abnormal values in the frequency division attribute, and fusing the frequency division attribute to describ