Shale deposition rhythm identification method and device and storage medium
阅读说明：本技术 一种页岩沉积韵律识别方法、装置及存储介质 (Shale deposition rhythm identification method and device and storage medium ) 是由 严德天 魏小松 张宝 刘紫璇 牛杏 李潼 于 2020-06-24 设计创作，主要内容包括：本发明属于页岩地层信息处理技术领域,公开了一种页岩沉积韵律识别方法、装置及存储介质,获取反映页岩沉积的韵律性变化的测试数据,存入数据库；获取替代指标,调用数据库数据,进行预处理；将预处理后的数据进行插值,将获得的等间距数据导入开源软件进行频谱分析,并识别出天文轨道周期；将识别出的天文轨道周期进行带通滤波,在目的层深度段标定顶底后将带通滤波,岩相以及数据曲线进行对比和调幅；依据岩心井或剖面识别的主控天文周期的波长识别和划分页岩沉积韵律性高频层序。本发明优化了泥页岩中的分析缺陷,首次在页岩地层中应用,解决了页岩地层高频韵律识别的难点,降低了分析成本。(The invention belongs to the technical field of shale stratum information processing, and discloses a method, a device and a storage medium for identifying shale deposition rhythm, which are used for acquiring test data reflecting rhythm change of shale deposition and storing the test data into a database; acquiring a substitution index, calling database data, and preprocessing; interpolating the preprocessed data, importing the obtained equidistant data into open source software for spectrum analysis, and identifying an astronomical orbit period; performing band-pass filtering on the identified astronomical orbit period, calibrating the top and the bottom of the target layer depth section, and then performing comparison and amplitude modulation on lithofacies and data curves; and identifying and dividing the shale deposition prosodic high-frequency sequence according to the wavelength of the main control astronomy period identified by the core well or the section. The method optimizes the analysis defects in the shale, is applied to the shale stratum for the first time, solves the difficulty of high-frequency rhythm identification of the shale stratum, and reduces the analysis cost.)
1. A shale deposition prosody recognition method is characterized by comprising the following steps:
acquiring test data reflecting prosodic changes of shale deposition, and storing the test data into a database; acquiring a substitution index, calling database data, and preprocessing;
interpolating the preprocessed data, importing the obtained equidistant data into open source software for spectrum analysis, and identifying an astronomical orbit period;
performing band-pass filtering on the identified astronomical orbit period, calibrating the top and the bottom of the target layer depth section, and then performing comparison and amplitude modulation on lithofacies and data curves;
and identifying the shale deposition prosody according to the wavelength of the master control astronomy period identified by the core well or the section.
2. The shale deposition prosody recognition method of claim 1, wherein the preprocessing step comprises: removing low-frequency background and high-frequency noise of the data, removing median and linear trend, and performing interpolation to enable data points to be equidistant;
the surrogate markers include: organic carbon content, stable oxygen isotopes, stable carbon isotopes, XRF element content, sediment color, acoustic waves, natural gamma, clay/dust content, magnetic susceptibility.
3. The shale deposition prosody recognition method of claim 1, wherein the call number spectral analysis comprises:
and calling 1-2 processed data of the substitution indexes for spectrum analysis, selecting open source software, performing spectrum analysis, evolving harmonic analysis, filtering and identifying the astronomical orbit period.
4. The shale deposition prosody recognition method of claim 1, wherein the shale deposition prosody recognition comprises:
comparing and amplitude modulating the astronomical orbit periodic filtering curve with the sedimentary stratum and surrogate index data curve, determining the top and bottom of the target layer section after the astronomical orbit periodic filtering curve is consistent with the sedimentary stratum and surrogate index data curve, and dividing the high-frequency rhythm.
5. The shale deposition prosody recognition method of claim 1, wherein the shale deposition prosody recognition further comprises: calibrating two or more wells or profiles, establishing a high-frequency stratum framework, carrying out primary top and bottom calibration on two or more wells or two profiles or one well and one profile, prolonging the identified sedimentary rhythm, and establishing the high-frequency sequence framework of the whole shale stratum.
6. The shale deposition prosody recognition method of claim 1, wherein after the shale deposition prosody is recognized, a filtering curve tuned by deposition rate is derived, and then accurate dating data is selected, a curve is calibrated, and a Laskar curve is selected for comparison;
the tuning method comprises the following steps: performing conversion of a depth domain and a time domain according to the deposition rate generated by software, then calibrating by using collected absolute annual data, comparing by using a theoretical Laskar curve, then obtaining an optimal astronomical convolution curve of a shale deposition target interval, and estimating the shale deposition duration;
and tuning the curve of the depth domain to the time domain, and identifying and dividing the shale deposition prosody in the time domain.
7. A shale deposition prosody recognition device, comprising:
the analysis end is used for calling data and preprocessing before spectrum analysis;
and the storage end is used for collecting various geological data and the data storage medium corresponding to various substitute indexes to be a database. The data corresponding to the substitute indexes can be directly used for calling after being preprocessed;
and the output end is used for outputting the required result.
8. The rock deposition prosody recognition device of claim 7, wherein the shale deposition prosody recognition device further comprises:
the processor 1 is used for performing spectrum analysis, classifying the called data, performing spectrum analysis and evolution harmonic analysis, identifying the astronomical orbit period, and performing band-pass filtering on the astronomical orbit period;
the processor 2 is used for high-frequency deposition prosody identification and division, determining the top and the bottom of the depth section of the target layer section, and then comparing and modulating the shale lithofacies with the called data curve and the called filter curve; when the amplitude modulation is consistent, dividing high-frequency rhythm according to wavelength (track period), and entering a processor 3 to build an astronomical cycle high-frequency sequence stratum framework; when the amplitude modulation is inconsistent, calling other index curves to pass through the steps of the processor 1 and the processor 2 again to enable the final results to be consistent;
and the processor 3 is used for building an astronomical cycle high-frequency sequence stratum framework, performing comparison work of multiple wells or multiple sections or sections and wells, and prolonging the divided high-frequency rhythm to realize high-frequency rhythm identification of the shale stratum and building the sequence stratum framework.
9. A rock deposition prosody recognition device as claimed in claim 7, wherein the output comprises:
the output end 1 is used for outputting a shale stratum high-frequency prosody graph;
the output end 2 is used for outputting an astronomical convolution high-frequency sequence stratigraphic framework diagram;
the output end 3 is used for outputting a time domain shale formation high-frequency prosody graph;
and the output end 4 is used for outputting the time domain shale formation high-frequency sequence trellis diagram.
10. A program storage medium storing a computer program for causing an electronic device to perform the method of rock deposition prosody recognition according to any one of claims 1 to 6, the method comprising the steps of:
step 1, obtaining test data reflecting prosodic changes of shale deposition and storing the test data into a database; calling database data, preprocessing, and obtaining a substitute index;
step 2, carrying out interpolation on the substitute indexes, importing the obtained equidistant data into open source software for spectrum analysis, and identifying the astronomical orbit period;
step 3, performing band-pass filtering on the identified astronomical orbit period, calibrating the top and the bottom of the depth section of the target layer, and then comparing and modulating the amplitude of lithofacies and data curves;
and 4, identifying the shale deposition prosody according to the wavelength of the master control astronomy period identified by the core well or the section.
Technical Field
The invention belongs to the technical field of shale stratum information processing, and particularly relates to a shale deposition rhythm identification method, a shale deposition rhythm identification device and a storage medium.
Background
Shale formations are different from sandstone and carbonate formations, and conventional sandstone or carbonate formations are not suitable for identification and division of sedimentary rhythm cycles of shale formations according to a sedimentary cycle division scheme established by field outcrop observation, log facies, earthquakes and the like. Although researchers constructively combine the Milnacidae Virginia theory and the sedimentary sequence to obtain a certain effect, obviously, for shale formations or strata which are stable in vertical variation and have no obvious logging response, the method for performing the gyratory partitioning by combining the conventional logging data with the Milnacidae Virginia gyratory theory has no effect, and even if other index data are used for replacing the conventional logging data, the partitioning method still catches the toggle. The scheme of prosodic identification and sequence segmentation of very paged shale layer systems by log data becomes almost impossible. This requires prosodic identification of the mechanisms by which the shale layer system is produced in combination with surrogate markers having variations.
Through the above analysis, the problems and defects of the prior art are as follows: (1) in the prior art, subjectivity exists when shale stratum deposition cycle division is specially aimed at, and the accuracy is low; in the prior art, the shale stratum high-frequency rhythm identification difficulty is large, and the analysis cost is high.
(2) The composition and the cause of the shale in different areas are different, the formation mechanism is different, but the formation cycle is divided by using the alternative indexes without being uniform, and the innovative technology is lacked.
(3) A more systematic, complete apparatus and a complete set of process technologies have not been developed.
The difficulty in solving the above problems and defects is: the shale stratum prosody identification needs large sampling density, high cost consumption, various substitution indexes and large workload. The significance of solving the problems and the defects is as follows: a set of special method, technology and device for identifying the shale deposition rhythm are formed, and great help is provided for theoretical innovation and guidance of unconventional oil and gas exploration.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a shale deposition rhythm identification method, a device and a storage medium.
The invention is realized in such a way that a shale deposition rhythm identification method comprises the following steps:
acquiring substitute index test data reflecting prosodic changes of shale deposition, and storing the substitute index test data into a database; and selecting a substitute index, calling database data, and preprocessing. The pretreatment steps are as follows: removing low-frequency background and high-frequency noise of the data, removing median and linear trend, and performing equidistant interpolation to enable data points to be equidistant;
importing the obtained equidistant data into open source software for spectrum analysis, and identifying an astronomical orbit period; the method comprises the steps that a plurality of peak values of spectrum analysis are obtained, the reciprocal of the frequency corresponding to each peak value is the wavelength, the ratio of the wavelength corresponding to the peak value is compared with the ratio of theoretical astronomical periods, and the astronomical period corresponding to a substitute index curve is identified;
performing band-pass filtering on the identified astronomical orbit period, wherein the band-pass filtering is filtering signals for filtering out eccentricity, slope and time lag periods, then calibrating the top and bottom of the target layer depth section, and comparing and modulating the amplitude of lithofacies and data curves;
and identifying the shale deposition prosody according to the wavelength of the master control astronomy period identified by the core well or the section.
Further, the method of pre-processing comprises: removing low-frequency background and high-frequency noise of data, removing median and linear trend, and performing interpolation to make the data points equidistant;
the surrogate markers include: organic carbon content, stable oxygen isotopes, stable carbon isotopes, XRF element content, sediment color, acoustic waves, natural gamma, clay/dust content, magnetic susceptibility.
Further, the invoking data spectrum analysis comprises:
calling 1-2 processed data capable of reflecting the substitute indexes of the shale prosody change for spectral analysis, selecting open source software, performing spectral analysis, evolving harmonic analysis, filtering and identifying the astronomical orbit period.
Further, the shale deposition prosody recognition comprises:
comparing and amplitude modulating the astronomical orbit periodic filtering curve with the sedimentary stratum and surrogate index data curve, determining the top and bottom of the target layer section after the astronomical orbit periodic filtering curve is consistent with the sedimentary stratum and surrogate index data curve, and dividing the high-frequency rhythm.
Further, the shale deposition prosody identification further comprises: calibrating two or more wells or profiles, establishing a high-frequency stratum framework, carrying out primary top and bottom calibration on two or more wells or two profiles or one well and one profile, prolonging the identified sedimentary rhythm, and establishing the high-frequency sequence framework of the whole shale stratum.
Further, after the shale deposition prosody is identified, a filtering curve tuned by the deposition rate is derived, then accurate dating data is selected, the curve is calibrated, and a Laskar curve is selected for comparison;
the tuning method comprises the following steps: performing conversion of a depth domain and a time domain according to the deposition rate generated by software, then calibrating by using collected absolute values, comparing by using a theoretical Laskar curve, then obtaining an optimal astronomical convolution curve of a shale deposition target interval, and estimating the shale deposition duration;
and tuning the curve of the depth domain to the time domain, and identifying and dividing the shale deposition prosody in the time domain.
Another object of the present invention is to provide a shale deposition prosody recognition apparatus, including:
the analysis end is used for calling data and preprocessing before spectrum analysis;
and the storage end is used for collecting various geological data and data corresponding to various substitute indexes and storing the data into a database. Data corresponding to the substitute indexes are preprocessed and directly used for calling;
and the output end is used for outputting the required result.
Further, the shale deposition prosody recognition device further comprises:
the
the
and the
Further, the output terminal includes:
the
the
the
and the output end 4 is used for outputting the time domain shale formation high-frequency sequence trellis diagram.
Another object of the present invention is to provide a program storage medium storing a computer program for causing an electronic device to execute the rock deposition prosody recognition method, including the steps of:
and 4, identifying the shale deposition prosody according to the wavelength of the master control astronomy period identified by the core well or the section.
By combining all the technical schemes, the invention has the advantages and positive effects that:
the embodiment of the invention collects a plurality of materials capable of reflecting the content change in the shale stratum in a certain research area, realizes a more reliable analysis result through an analysis library, substitution index selection, database calling, spectrum analysis and the like, is used for high-frequency prosody identification of the shale stratum and establishment of a high-frequency astronomical gyrus stratigraphic framework, and optimizes the analysis defects in the shale; the invention can make up the problem of difficult prosody recognition by a single index by mutually combining a plurality of alternative indexes.
In the embodiment of the invention, prosody recognition and division are respectively carried out on a large set of shale stratum by using two different indexes (GR and Fe elements). In large sets of shale formations where GR curves do not change significantly, it is possible to identify prosodic convolutions (fig. 4). For very short shale formations (below 20-40 m), it is almost impossible to identify high-precision prosodic gyre, and prosodic changes are well identified by using other substitution indexes Fe (figure 5).
The comparison of the two experimental effects provided by the embodiment of the invention shows that in the shale sedimentary stratum, different substitute indexes have different effects on prosody recognition, and because a plurality of index data are difficult to obtain in the sedimentary stratum, the substitute index data of a plurality of regions need to be collected as much as possible to establish a database, so that the database can be directly called in the subsequent prosody recognition, the workload and the difficulty are reduced, and meanwhile, the saving and effective effects are achieved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings needed to be used in the embodiments of the present application will be briefly described below, and it is obvious that the drawings described below are only some embodiments of the present application, and it is obvious for those skilled in the art that other drawings can be obtained from the drawings without creative efforts.
FIG. 1 is a flow chart of a shale deposition prosody recognition method according to an embodiment of the present invention
FIG. 2 is a schematic diagram of the natural gamma curve (GR) preprocessing for the alternative index according to the embodiment of the present invention.
Fig. 3 is a spectrum analysis diagram of a section of a shale in the Yunnan area according to an embodiment of the present invention.
Fig. 4 is a schematic view of an astronomical convolution high-frequency sequence of a shale formation depressed in the southwest of the north bay basin Wei according to an embodiment of the present invention.
Fig. 5 is a schematic diagram of prosody recognition of a shale formation in a certain interval of the north bay basin according to an embodiment of the present invention.
Fig. 6 is a schematic diagram of a device for identifying shale deposition prosody according to an embodiment of the present invention.
Fig. 7 is a schematic interface diagram of a shale deposition prosody recognition device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.