CN111650670A

CN111650670A - Shale deposition rhythm identification method and device and storage medium

Info

Publication number: CN111650670A
Application number: CN202010587149.XA
Authority: CN
Inventors: 严德天; 魏小松; 张宝; 刘紫璇; 牛杏; 李潼
Original assignee: China University of Geosciences
Current assignee: China University of Geosciences
Priority date: 2020-06-24
Filing date: 2020-06-24
Publication date: 2020-09-11
Anticipated expiration: 2040-06-24
Also published as: CN111650670B

Abstract

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.

Description

Shale deposition rhythm identification method and device and storage medium

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 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 the wavelength, 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.

Further, the output terminal includes:

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.

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:

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.

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.

In the prior art, the method is specially aimed at the subjectivity existing in the shale stratum deposition cycle division, 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.

Aiming at the problems in the prior art, the invention provides a method, a device and a storage medium for identifying the prosody of shale deposition, and the invention is described in detail with reference to the accompanying drawings.

As shown in fig. 1, the method for identifying the prosody of shale deposition provided by the present invention includes:

s101: analyzing the aim of the shale deposition prosody identification method to be realized, then selecting a proper substitute index, preprocessing the data of the substitute index, and using the preprocessed data for the subsequent spectrum analysis function.

S102: the interpolated equidistant data were imported into open source software for spectral analysis and possible astronomical orbit periods were identified (eccentricities 405kyr and 100kyr, slopes 52kyr and 40kyr and ages 23kyr,22kyr and 18 kyr).

S103: and performing band-pass filtering on the identified astronomical orbit period, then calibrating the top and the bottom of the target layer depth section, and then comparing and modulating the amplitude of the lithofacies and the data curve.

S104: and identifying the shale deposition prosody according to the wavelength of the master astronomical period identified by the typical core well or the section.

S105: and (3) deriving a filtering curve after the deposition rate is tuned (the deposition rate and the tuning process are finished together, and only operation is needed), then selecting accurate dating data, calibrating the curve, and selecting a theoretical Laskar (04,10a,10b and 11) curve for comparison.

The invention is further described with reference to specific examples.

Examples

The shale deposition prosody identification method provided by the embodiment of the invention comprises the following steps:

firstly, for a field section, whether a collected sample exists, whether a section picture and a sample picture exist, whether analysis test data exist and the like are determined.

For well drilling, whether the well is a typical well with a rock core or not, whether a photo exists or not, whether test analysis data exists or not, whether well logging data exists or not and the like are determined.

The shale photo can directly observe the change of the deposition prosody of the shale, has different colors, represents the difference of internal sediment elements, and can obtain two indexes of gray scale and reflectivity from the colors.

Second, test data reflecting prosodic changes in shale deposits are collected or tested and included in the database.

Furthermore, according to the data of the analysis library, a suitable alternative index is selected, and the alternative index selection column has the following indexes: organic carbon content (TOC), stable oxygen isotopes, stable carbon isotopes, XRF elemental content, sediment color, sonic waves (DT), natural Gamma Rays (GR), clay/dust content, magnetic susceptibility, etc.

Moreover, the preprocessing of the data in the database mainly includes: the acquired various data are put into a database and preprocessed, and the preprocessing process comprises the following steps: and removing low-frequency background and high-frequency noise of the data, removing median and linear trend, and enabling data points to be equidistant through interpolation.

Furthermore, the data in the database is called for spectrum analysis, which mainly comprises: calling the data of 1-2 processed substitution indexes for spectrum analysis, selecting proper open source software, performing spectrum analysis (MTM), Evolution Harmonic Analysis (EHA), filtering (Bandbass) and the like, and identifying the astronomical orbit period.

Furthermore, the identified astronomical orbit period is used for shale deposition prosody identification, and the method mainly comprises the following steps: determining the top and the bottom of the target layer section, comparing and amplitude modulating the identified main astronomical orbit period filter curve with the sedimentary stratum and the substitution index data curve, and dividing the high-frequency rhythm after the comparison and amplitude modulation are consistent. Amplitude modulation here means that the amplitude of the surrogate index curve and the shale lithofacies variation is consistent with the amplitude of the filter curve.

The basis of the sedimentary prosody is that an astronomical convolution period can be used as a convolution of a certain level of sequence.

In the invention, two or more wells or sections are calibrated to establish a high-frequency stratum framework, which mainly comprises the following steps:

and (3) carrying out primary top and bottom calibration on two or more wells (one of which has a rock core or a plurality of data and is subjected to spectral analysis), or two profiles (one of which is a research profile and the other of which can be a database profile), or one well and one profile, and prolonging the identified sedimentary prosody to establish a high-frequency sequence trellis of the whole shale stratum.

The invention can further tune the curve of the depth domain to the time domain, and the identification and division of the shale deposition prosody are carried out in the time domain. If the concept of time is involved, support for accurate dating data, such as isotopes, zircon dating, magnetic data, etc., is required.

The method used for tuning is to perform conversion of a depth domain and a time domain according to the deposition rate generated by software, then calibrate by using collected absolute annual data, and also compare with a theoretical Laskar (04,10a,10b and 11) curve, then obtain an optimal astronomical convolution curve of the shale deposition target interval, and estimate the duration of the shale deposition.

The invention also provides a device for identifying the shale deposition prosody, which comprises:

and the analysis end is used for data gating in the early stage and preparing before calling data and performing spectrum analysis. Specifically how many wells, several sections, how many samples, whether photographs are available, how many logging information, what test data are, how many sampling intervals can be made to correspond to or be less than the interval in the shale deposition color change? Suitable for the operation of the subsequent processor 1 with what data, etc.

And the storage end is used for collecting various geological data and data corresponding to various substitute indexes, and the data corresponding to the substitute indexes can be directly used for calling after being preprocessed. Specifically, data of the substitute indexes needs to be collected and imported into a database as much as possible, and then placed into a storage end.

The output mainly outputs the required result (graph).

And an output end 1 mainly outputs a shale stratum high-frequency prosody graph.

And the output end 2 is mainly used for outputting an astronomical convolution high-frequency sequence stratigraphic framework diagram.

And an output end 3 outputs a time domain shale formation high-frequency prosody graph.

And an output end 4 for outputting the time domain shale stratum high-frequency sequence trellis diagram.

The following are collection and processing schemes for surrogate markers: the first kind of data of organic carbon content (TOC), stable oxygen isotope, stable carbon isotope and clay/dust content should collect data of foreigners in a research area or an adjacent area as much as possible, and then coverage sampling and testing are carried out on horizons without data. Such indicators may be preferred if ancient climate studies are to be conducted after identifying the prosody and compartmentalization of the shale deposits.

The second type: XRF element content, magnetic susceptibility, etc.

The second type of index does not need sampling measurement, and is directly measured on the surface of a sample by an instrument, the sample indicates that smooth and fresh surfaces are required to be ensured, the sampling interval is as small as possible, and the sampling interval is ensured to be between 3cm and 5 cm.

In the third category: and (3) observing the change of the shale stratum, then photographing the shale stratum, and extracting the gray scale or color reflectivity through open source software. The shale deposition color is the visual reflection of the shale stratum change, and the prosody of the shale deposition can be rapidly and accurately identified and a plurality of gyrations can be divided by the gray scale or the color reflectivity of the index in the stratum with obvious change.

The fourth type: sonic (DT), natural Gamma (GR), sonic and natural gamma curves are commonly used in well data, and measurements are typically taken manually with instruments on the profile. The acoustic waves and the natural gamma rays can reflect the content of radioactive elements in the stratum, the content can reflect the cause of the stratum, and therefore the content can be used for identifying the sedimentary prosody.

The processor in the shale formation prosody recognition device provided by the invention has the function of processing all data in the database.

The processor 1 is mainly used for spectrum analysis. The method comprises the steps of firstly classifying called data, then clicking a spectrum analysis function to respectively perform spectrum analysis (MTM), Evolution Harmonic Analysis (EHA) and identification of an astronomical orbit period, and then respectively performing band-pass filtering (Bandpass) on the astronomical orbit period.

The processor 2 is mainly used for high-frequency deposition prosody recognition and division functions. The top and the bottom of a depth section of a target interval are mainly determined, and then shale lithofacies and the called data curve, the filter curve are compared and amplitude modulated. When the amplitude modulation is consistent, high-frequency rhythm is divided according to the wavelength, and the high-frequency rhythm enters the processor 3 to perform an astronomical cycle high-frequency sequence stratum framework building function. And when the amplitude modulation is inconsistent, calling other index curves to pass through the steps of the processor 1 and the processor 2 again, so that the final result is consistent and the reliability is higher.

The processor 3 mainly has the astronomical convolution high-frequency sequence stratum framework building function. And receiving the unfinished work of the processor 2, performing the comparison work of a plurality of wells or a plurality of sections or sections and wells, prolonging the divided high-frequency prosody, and finally realizing the high-frequency prosody recognition of the shale stratum and the construction of a stratum framework.

The invention is further described with reference to the accompanying drawings and detailed analysis.

In the present invention, fig. 2 is a schematic diagram of the natural gamma curve (GR) preprocessing of the substitute index provided by the embodiment of the present invention.

The main alternative indexes currently used are: organic carbon content (TOC), stable oxygen isotopes, stable carbon isotopes, XRF elemental content, sediment color, sonic waves (DT), natural Gamma Rays (GR), clay/dust content, magnetic susceptibility, etc.; wherein organic carbon content (TOC), stable oxygen isotopes and carbon isotopes, XRF element content, sediment color, and the like are preferred for shale formations. Generally, data acquisition of shale formations may have skip points or abnormal values, so before performing spectrum analysis, the abnormal values are removed. And removing trend, and after removing extreme values, interpolation is needed to enable data to be equidistant. For the well logging curve, because a lot of noises and other signals are contained, the noises are also reduced in advance, and low-frequency background and high-frequency noises are processed. Fig. 2 is a schematic diagram of noise reduction processing of a natural gamma curve in a substitute index provided in an embodiment of the present invention, and is implemented by matlab software.

In the present invention, fig. 3 is a diagram of spectrum analysis provided by the embodiment of the present invention.

In the example, the content of Ca in XRF elements of certain shale layers in Yunnan is selected for spectrum analysis. After spectral analysis, 4 peaks exceeding 95% of confidence were found, wherein the peaks are respectively 30cm,15.9cm,14.3cm and 7cm, and the ratio of the peaks to the peaks is 30: 15.9: 14.3: 1-8: 0.53: 0.477: 0.267, this ratio is similar to the astronomical orbit period ratio of 100kyr:52kyr:48kyr:23kyr, thus determining the main astronomical periods of 100kyr,52kyr,28kyr,23 kyr.

In the present invention, fig. 4 is a schematic view of an astronomical convolution high-frequency sequence of a mud shale formation depressed in the southwest of the north gulf basin Wei according to an embodiment of the present invention.

The example is a high-frequency sequence identification case which is carried out by determining an astronomical orbit period after spectral analysis by using a natural gamma curve and then filtering. According to the flow from the analysis library to the spectrum analysis library in the right, the natural gamma curve is preprocessed, and after calling, spectrum analysis is carried out to identify eccentricity, slope and age period. The method comprises the steps of selecting a 405kyr eccentricity period of a main control for filtering and amplitude modulation, comparing a GR curve with a lithofacies, preliminarily dividing the shale in a target interval into 14 five-level sedimentary sequences, calibrating the top and bottom of the example before implementation, and correcting the top and bottom of the example with a theoretical Laskar (2004, 2011) curve.

In the present invention, fig. 5 is a schematic diagram of prosodic recognition classification of a shale interval in the northern gulf basin according to an embodiment of the present invention.

The example is a case of determining an astronomical orbit period by using a tested Fe element after spectral analysis, and then filtering to perform high-precision prosody recognition on a shale sedimentary stratum. According to the flow from the analysis library to the spectrum analysis library in the patent application, after Fe element data are preprocessed, the eccentricity, the slope and the age period are mainly identified through spectrum analysis. The slope period 41kyr is selected for filtering and amplitude modulation to achieve the effect of identifying and demarcating prosodic loops in short shale intervals without log response.

In the present invention, fig. 6 is a schematic diagram of an apparatus for identifying shale deposition prosody according to an embodiment of the present invention. The left part of fig. 5 is divided into 5 libraries of analysis library, surrogate index selection column, database, spectrum analysis library and depositional cycle partitioning-trellis building library, whose functions are consistent with the procedure in example 1. The right part of fig. 5 is a device for shale deposition rhythm recognition, which has an analysis end, a storage end, a processor 1, a processor 2, a processor 3, an output end 1, and an output end 2. The processor operates according to the sequence of the processor 1, the processor 2 and the processor 3. The content of the left part is essentially the core carrier of the right part of the device. An analysis library is arranged in the analysis end and is mainly used for distinguishing wells or sections, judging whether samples are collected or not, pictures are collected, analysis data are tested, whether the typical wells exist or not can represent a sedimentary basin, whether well logging curves exist or not and the like. The processed database is arranged in the storage end, and the processed database are connected into the processor together. The analysis end wire can select the corresponding substitute index in the substitute index selection column. The

processor terminals

1, 2, 3 have three main functions: the processor 1 determines an astronomical orbit period for spectral analysis and extracts a corresponding period to obtain filtering; the processor 2 performs astronomical convolution high-frequency division on the well or the section; the processor 3 is used for well-to-well, well-to-profile or profile-to-profile contrast calibration and high-frequency sequence stratum trellis construction. The output end 1 is an output shale deposition prosody graph, and the output end 2 is an output shale stratum high-frequency sequence stratum framework building graph.

In the present invention, fig. 7 is a schematic interface diagram of a shale deposition prosody recognition apparatus designed in the present invention, mainly for showing the function and application of the apparatus. The schematic diagram represents only the functional interface of the device of the present invention, and is not a real device interface. Various modifications and other effective device component interface optimizations may be made thereto without departing from the scope of the device interface of the present invention.

In the embodiments provided in the present application, the realized apparatus can also be realized by other means. The device of the embodiment of the invention is only a schematic diagram, and particularly needs to be optimized in practice. For example, when the analysis end module and the storage end module enter the processor, appropriate codes are needed, a substitute index bar is selected, and the execution code part of the database is called to be optimized. The

processors

1, 2, 3 respectively execute different operations, and the input and output functions in each processor and the combination function with other software need to be optimized, so that the required purpose can be finally achieved. The instructions executed by the functional modules in the device can be stored in a computer storage medium when the functional modules are sold or used as independent products. In essence, the inventive apparatus may ultimately be run in a computer as software.

In the description of the present invention, "a plurality" means two or more unless otherwise specified; the terms "upper", "lower", "left", "right", "inner", "outer", "front", "rear", "head", "tail", and the like, indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, are only for convenience in describing and simplifying the description, and do not indicate or imply that the device or element referred to must have a particular orientation, be constructed in a particular orientation, and be operated, and thus, should not be construed as limiting the invention. Furthermore, the terms "first," "second," "third," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.

It should be noted that the embodiments of the present invention can be realized by hardware, software, or a combination of software and hardware. The hardware portion may be implemented using dedicated logic; the software portions may be stored in a memory and executed by a suitable instruction execution system, such as a microprocessor or specially designed hardware. Those skilled in the art will appreciate that the apparatus and methods described above may be implemented using computer executable instructions and/or embodied in processor control code, such code being provided on a carrier medium such as a disk, CD-or DVD-ROM, programmable memory such as read only memory (firmware), or a data carrier such as an optical or electronic signal carrier, for example. The apparatus and its modules of the present invention may be implemented by hardware circuits such as very large scale integrated circuits or gate arrays, semiconductors such as logic chips, transistors, or programmable hardware devices such as field programmable gate arrays, programmable logic devices, etc., or by software executed by various types of processors, or by a combination of hardware circuits and software, e.g., firmware.

The above description is only for the purpose of illustrating the present invention and the appended claims are not to be construed as limiting the scope of the invention, which is intended to cover all modifications, equivalents and improvements that are within the spirit and scope of the invention as defined by the appended claims.

Claims

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;

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;

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:

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;

7. A shale deposition prosody recognition device, comprising:

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 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;

9. A rock deposition prosody recognition device as claimed in claim 7, wherein the output comprises:

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: