Frequency division phase control parameter model establishing method and device

Abstract

The invention discloses a frequency division phase control parameter model establishing method and device. The method comprises the following steps: acquiring a frequency division low-frequency seismic trace, a frequency division intermediate-frequency seismic trace, a frequency division high-frequency seismic trace, a frequency division low-frequency model trace, a frequency division intermediate-frequency model trace, a frequency division high-frequency model trace and a low-frequency parameter model trace of each trace to be calculated in a research area; respectively calculating root-mean-square amplitudes of the frequency-division low-frequency, intermediate-frequency and high-frequency seismic traces, and respectively determining frequency-division low-frequency, intermediate-frequency and high-frequency weight coefficients according to the obtained root-mean-square amplitudes; combining the frequency division low-frequency intermediate-frequency and high-frequency model channels and the low-frequency parameter model channel by utilizing each weight coefficient to obtain a frequency division phase control parameter model channel; and obtaining a frequency division phase control parameter model data body according to each frequency division phase control parameter model channel. According to the method, the response characteristics of different frequency-division seismic data to the reservoir sedimentary structure are fully excavated, the high-precision parameter model with the transverse change fully referring to the seismic characteristics is constructed, the seismic facies control degree of the model can be improved, and then the reliability of reservoir prediction is improved.

Description

Technical Field

[0001] The present invention relates to the field of geophysical exploration technology, and in particular to a method and device for establishing a frequency-division phase-controlled parameter model. Background Art

[0002] In the actual work process of reservoir quantitative prediction, it is often encountered that the wave impedance curve cannot effectively identify the reservoir, resulting in the application of wave impedance seismic inversion technology in the study area failing to achieve the purpose of reservoir quantitative identification. Based on the above situation, the relevant technology of using characteristic curves that can effectively identify reservoirs to build parameter models and conduct quantitative prediction of reservoirs has developed rapidly, which is an important development direction to make up for the shortcomings of wave impedance inversion technology.

[0003] At present, the technical method for constructing parameter models is mainly to interpolate well curves into parameter models based on the construction of grid models using mathematical algorithms such as local weighting, simple Kriging or variogram. The above methods can quickly form multi-type parameter models, but they mainly use well logging data in the modeling process, and their spatial changes are unrelated to seismic data. Therefore, there are major problems in the spatial predictability of parameter models.

[0004] With the further deepening of oil and gas exploration and development, the proportion of complex oil and gas reservoirs that cannot be identified by wave impedance is increasing. There is an urgent need for a parameter model establishment method that can effectively combine well seismic data to improve the accuracy of reservoir prediction and thus improve the success rate of exploration and development. Summary of the invention

[0005] In order to at least partially solve the technical problems existing in the prior art, the inventors have made the present invention, and through a specific implementation method, provide a method and device for establishing a frequency-fractionated phase-controlled parameter model, which fully explores the response characteristics of different frequency-fractionated seismic data to reservoir sedimentary structure, and constructs a high-precision parameter model with full reference to seismic characteristics in lateral changes, which can effectively improve the seismic phase control degree of the parameter model, thereby improving the reliability of reservoir prediction.

[0006] In a first aspect, an embodiment of the present invention provides a method for establishing a frequency division phase control parameter model, comprising:

[0007] Obtain the divided low-frequency seismic trace, divided medium-frequency seismic trace, divided high-frequency seismic trace, divided low-frequency model trace, divided medium-frequency model trace, divided high-frequency model trace and low-frequency parameter model trace of each trace to be calculated in the study area;

[0008] Calculating the root mean square amplitude of the divided-frequency low-frequency seismic trace, the root mean square amplitude of the divided-frequency intermediate-frequency seismic trace, and the root mean square amplitude of the divided-frequency high-frequency seismic trace, and determining the divided-frequency low-frequency weight coefficient, the divided-frequency intermediate-frequency weight coefficient, and the divided-frequency high-frequency weight coefficient according to the obtained root mean square amplitudes;

[0009] By using the frequency division low frequency weight coefficient, the frequency division intermediate frequency weight coefficient and the frequency division high frequency weight coefficient, the frequency division low frequency model channel, the frequency division intermediate frequency model channel, the frequency division high frequency model channel and the low frequency parameter model channel are merged to obtain the frequency division phase control parameter model channel of the channel to be calculated;

[0010] A frequency-division phase-control parameter model data body is obtained from the frequency-division phase-control parameter model channels of each channel to be calculated.

[0011] In a second aspect, an embodiment of the present invention provides a frequency division phase control parameter model establishment device, comprising:

[0012] The data acquisition module is used to obtain the divided low-frequency seismic traces, divided medium-frequency seismic traces, divided high-frequency seismic traces, divided low-frequency model traces, divided medium-frequency model traces, divided high-frequency model traces and low-frequency parameter model traces of each trace to be calculated in the study area;

[0013] A frequency division weight coefficient determination module is used to calculate the root mean square amplitude of the frequency division low-frequency seismic trace, the root mean square amplitude of the frequency division medium-frequency seismic trace and the root mean square amplitude of the frequency division high-frequency seismic trace, and determine the frequency division low-frequency weight coefficient, the frequency division medium-frequency weight coefficient and the frequency division high-frequency weight coefficient according to the obtained root mean square amplitudes;

[0014] A frequency division phase control parameter model channel establishment module is used to combine the frequency division low frequency model channel, the frequency division medium frequency model channel, the frequency division high frequency model channel and the low frequency parameter model channel by using the frequency division low frequency weight coefficient, the frequency division medium frequency weight coefficient and the frequency division high frequency weight coefficient to obtain the frequency division phase control parameter model channel of the channel to be calculated;

[0015] The frequency division phase control parameter model data body establishment module is used to obtain the frequency division phase control parameter model data body from the frequency division phase control parameter model channels of each channel to be calculated.

[0016] In a third aspect, an embodiment of the present invention provides a computer storage medium, wherein the computer storage medium stores computer executable instructions, and when the computer executable instructions are executed by a processor, the above-mentioned method for establishing a frequency division phase control parameter model is implemented.

[0017] In a fourth aspect, an embodiment of the present disclosure provides a server, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for establishing a frequency division phase control parameter model when executing the program.

[0018] The beneficial effects of the above technical solution provided by the embodiment of the present invention include at least:

[0019] The method for establishing a frequency-divided phase-controlled parameter model provided in an embodiment of the present invention extracts frequency-divided seismic traces, frequency-divided model traces and low-frequency parameter model traces, and then calculates the frequency-divided weight coefficients of different frequency band parameters with reference to the frequency-divided seismic attribute values, and finally uses the frequency-divided weight coefficients to merge the frequency-divided low-frequency model trace, the frequency-divided medium-frequency model trace, the frequency-divided high-frequency model trace and the low-frequency parameter model trace into a frequency-divided phase-controlled parameter model. This method fully exploits the response characteristics of different frequency-divided seismic data to reservoir sedimentary structures, constructs a high-precision parameter model with full reference to seismic characteristics for lateral changes, and can effectively improve the seismic phase control degree of the parameter model, thereby improving the reliability of reservoir prediction.

[0020] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures particularly pointed out in the written description, claims, and drawings.

[0021] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the present invention and do not constitute a limitation of the present invention. In the accompanying drawings:

[0023] Figure 1 It is a flow chart of a method for establishing a multi-frequency phase-controlled parameter model in an embodiment of the present invention;

[0024] Figure 2 It is the frequency-divided low-frequency seismic trace of a certain trace to be calculated in a certain block;

[0025] Figure 3 is the frequency-divided intermediate frequency seismic channel of the channel to be calculated;

[0026] Figure 4 is the frequency-divided high-frequency seismic channel of the channel to be calculated;

[0027] Figure 5 is the frequency-divided low-frequency model channel of the channel to be calculated;

[0028] Figure 6 is the divided intermediate frequency model channel of the channel to be calculated;

[0029] Figure 7 is the frequency-divided high-frequency model channel of the channel to be calculated;

[0030] Figure 8 is the low-frequency parameter model channel of the channel to be calculated;

[0031] Fig. 9A frequency division phase control parameter model established for an embodiment of the present invention;

[0032] Fig.10 is the simple Kriging parameter model profile;

[0033] Fig.11 A cross section of a frequency division phase control parameter model data volume established in an embodiment of the present invention;

[0034] Fig.12 It is a structural schematic diagram of a device for establishing a frequency division phase control parameter model in an embodiment of the present invention. DETAILED DESCRIPTION

[0035] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0036] It should be understood that the terms described in the present invention are only for describing special embodiments and are not intended to limit the present invention. In addition, for the numerical range in the present invention, it should be understood that each intermediate value between the upper and lower limits of the scope is also specifically disclosed. Each smaller range between the intermediate value in any stated value or stated range and any other stated value or intermediate value in the described range is also included in the present invention. The upper and lower limits of these smaller ranges can be independently included or excluded in the scope.

[0037] Unless otherwise indicated, all technical and scientific terms used herein have the same meanings as commonly understood by those skilled in the art to which the invention belongs. Although the present invention describes only preferred methods and materials, any methods and materials similar or equivalent to those described herein may also be used in the implementation or testing of the present invention. All documents mentioned in this specification are incorporated by reference to disclose and describe the methods and / or materials related to the documents. In the event of a conflict with any incorporated document, the content of this specification shall prevail.

[0038] Example

[0039] The embodiment of the present invention provides a method for establishing a frequency division phase control parameter model, and the process thereof is as follows Figure 1 As shown, the following steps are included:

[0040] Step S11: Obtain the divided-frequency low-frequency seismic channel, divided-frequency medium-frequency seismic channel, divided-frequency high-frequency seismic channel, divided-frequency low-frequency model channel, divided-frequency medium-frequency model channel, divided-frequency high-frequency model channel and low-frequency parameter model channel of each channel to be calculated in the study area.

[0041] According to the structural grid and synthetic record method, the divided low-frequency seismic traces, divided medium-frequency seismic traces, divided high-frequency seismic traces, divided low-frequency model traces, divided medium-frequency model traces, divided high-frequency model traces and low-frequency parameter model traces of each trace to be calculated in the study area are extracted.

[0042] Structural framework: refers to the structural framework built by seismic horizon interpretation, which is used to constrain the lateral changes of model parameters. Synthetic records: used for time-depth calibration. The parameters obtained by well logging are in the depth domain, and the seismic parameters are in the time domain. The depth domain is converted to the depth domain through synthetic records, and the well and seismic are matched.

[0043] The above two technologies can be carried out using existing conventional technologies.

[0044] Figure 2 It is the frequency-divided low-frequency seismic trace of a certain block to be calculated. Figure 3 is the frequency-divided intermediate frequency seismic channel of the channel to be calculated, Figure 4 is the frequency-divided high-frequency seismic trace of the trace to be calculated, Figure 5 is the frequency-divided low-frequency model channel of the channel to be calculated, Figure 6 is the frequency-divided intermediate frequency model channel of the channel to be calculated, Figure 7 is the frequency-divided high-frequency model channel of the channel to be calculated, Figure 8 is the low-frequency parameter model channel of the channel to be calculated.

[0045] Step S12: Calculate the root mean square amplitude of the divided low-frequency seismic trace, the root mean square amplitude of the divided medium-frequency seismic trace, and the root mean square amplitude of the divided high-frequency seismic trace.

[0046] Calculate the sum of the squares of the sampling points in the divided-frequency low-frequency seismic trace, determine the ratio of the sum to the number of sampling points in the divided-frequency low-frequency seismic trace, and determine the square root of the ratio as the root mean square amplitude of the divided-frequency low-frequency seismic trace:

[0047]

[0048] In formula (1), Al is the root mean square amplitude of the low-frequency seismic trace, L Sl is the wavelength, that is, the number of sampling points of the wave, and also the number of sampling points of the frequency-divided low-frequency seismic trace. i is the value of the ith sampling point of the divided low-frequency seismic trace.

[0049] Calculate the sum of the squares of the sampling points in the divided-frequency IF seismic trace, determine the ratio of the sum to the number of sampling points in the divided-frequency IF seismic trace, and determine the square root of the ratio as the root mean square amplitude of the divided-frequency IF seismic trace:

[0050]

[0051] In formula (2), Am is the root mean square amplitude of the intermediate frequency seismic trace, L Smis the wavelength, that is, the number of sampling points of the wave, and also the number of sampling points of the frequency-divided intermediate frequency seismic trace. i is the value of the i-th sampling point of the divided intermediate frequency seismic trace.

[0052] Calculate the sum of the squares of the sampling points in the divided-frequency high-frequency seismic trace, determine the ratio of the sum to the number of sampling points in the divided-frequency high-frequency seismic trace, and determine the square root of the ratio as the root mean square amplitude of the divided-frequency high-frequency seismic trace:

[0053]

[0054] In formula (3), Ah is the root mean square amplitude of the frequency-divided high-frequency seismic trace, L sh is the wavelength, that is, the number of sampling points of the wave, and also the number of sampling points of the frequency-divided high-frequency seismic trace. Sh i is the value of the i-th sampling point of the divided high-frequency seismic trace.

[0055] The calculation results of the root mean square amplitude of the divided low-frequency seismic trace, the root mean square amplitude of the divided medium-frequency seismic trace and the root mean square amplitude of the divided high-frequency seismic trace are shown in Table 1 below:

[0056] Table 1 Calculation results of RMS amplitude of frequency-divided seismic traces

[0057] RMS amplitude of the divided low-frequency seismic trace RMS amplitude of frequency-divided high-frequency seismic traces RMS amplitude of the frequency-divided intermediate frequency seismic trace 2169.23 8549.32 3482.48

[0058] Step S13: determining the frequency division low-frequency weight coefficient, the frequency division middle-frequency weight coefficient and the frequency division high-frequency weight coefficient according to the obtained root mean square amplitudes.

[0059] The sum of the obtained root mean square amplitudes is determined, and the ratio of the root mean square amplitude of the divided low frequency seismic trace to the sum is determined as the divided low frequency weight coefficient Ql:

[0060]

[0061] The ratio of the RMS amplitude of the intermediate frequency seismic trace to the sum is determined as the intermediate frequency weight coefficient Qm:

[0062]

[0063] The ratio of the RMS amplitude of the divided high-frequency seismic trace to the sum is determined as the divided high-frequency weight coefficient Qh:

[0064]

[0065] The calculation results of the frequency division low frequency weight coefficient, frequency division middle frequency weight coefficient and frequency division high frequency weight coefficient of the above-mentioned channel to be calculated are shown in the following Table 2:

[0066] Table 2 Calculation results of frequency division weight coefficient

[0067] Crossover low frequency weight coefficient Crossover intermediate frequency weight coefficient Frequency division high frequency weight coefficient 0.24 0.39 0.38

[0068] Step S14: using the frequency-divided low-frequency weight coefficient, the frequency-divided intermediate-frequency weight coefficient and the frequency-divided high-frequency weight coefficient, the frequency-divided low-frequency model channel, the frequency-divided intermediate-frequency model channel, the frequency-divided high-frequency model channel and the low-frequency parameter model channel are merged to obtain the frequency-divided phase-controlled parameter model channel of the channel to be calculated.

[0069] The following formula (7) can be used to combine the divided-frequency low-frequency model channel, the divided-frequency intermediate-frequency model channel, the divided-frequency high-frequency model channel and the low-frequency parameter model channel to obtain the divided-frequency phase-controlled parameter model channel of the channel to be calculated:

[0070] Z=Cl+Ql*Ml+Qm*Mm+Qh*Mh (7)

[0071] In formula (7), Z is the frequency division phase control parameter model channel, Cl is the low-frequency parameter model channel, Ql is the frequency division low-frequency weight coefficient, Ml is the frequency division low-frequency model channel, Qm is the frequency division intermediate frequency weight coefficient, Mm is the frequency division intermediate frequency model channel, Qh is the frequency division high-frequency weight coefficient, and Mh is the frequency division high-frequency model channel.

[0072] The frequency division phase control parameter model of the above channel to be calculated is as follows: Fig. 9 shown.

[0073] Step S15: Obtain a frequency-division phase-control parameter model data body from the frequency-division phase-control parameter model channels of each channel to be calculated.

[0074] The method for establishing a frequency-divided phase-controlled parameter model provided in an embodiment of the present invention extracts frequency-divided seismic traces, frequency-divided model traces and low-frequency parameter model traces, and then calculates the frequency-divided weight coefficients of different frequency band parameters with reference to the frequency-divided seismic attribute values, and finally uses the frequency-divided weight coefficients to merge the frequency-divided low-frequency model trace, the frequency-divided medium-frequency model trace, the frequency-divided high-frequency model trace and the low-frequency parameter model trace into a frequency-divided phase-controlled parameter model. This method fully exploits the response characteristics of different frequency-divided seismic data to reservoir sedimentary structures, constructs a high-precision parameter model with full reference to seismic characteristics for lateral changes, and can effectively improve the seismic phase control degree of the parameter model, thereby improving the reliability of reservoir prediction.

[0075] See also Fig.10 and Fig.11 As shown, they are respectively the profile of the simple Kriging parameter model and the profile of the frequency-divided phase-controlled parameter model data volume established in this application. Comparison between the two shows that the frequency-divided phase-controlled parameter model has richer spatial variations and higher vertical resolution, which can effectively improve the accuracy of parameter modeling.

[0076] Based on the inventive concept of the present invention, an embodiment of the present invention further provides a frequency division phase control parameter model establishment device, the structure of which is as follows: Fig.12 As shown, including:

[0077] The data acquisition module 121 is used to obtain the divided low-frequency seismic trace, divided medium-frequency seismic trace, divided high-frequency seismic trace, divided low-frequency model trace, divided medium-frequency model trace, divided high-frequency model trace and low-frequency parameter model trace of each trace to be calculated in the study area;

[0078] The frequency division weight coefficient determination module 122 is used to calculate the root mean square amplitude of the frequency division low-frequency seismic trace, the root mean square amplitude of the frequency division medium-frequency seismic trace and the root mean square amplitude of the frequency division high-frequency seismic trace, and determine the frequency division low-frequency weight coefficient, the frequency division medium-frequency weight coefficient and the frequency division high-frequency weight coefficient according to the obtained root mean square amplitudes;

[0079] The frequency division phase control parameter model channel establishment module 123 is used to combine the frequency division low frequency model channel, the frequency division medium frequency model channel, the frequency division high frequency model channel and the low frequency parameter model channel by using the frequency division low frequency weight coefficient, the frequency division medium frequency weight coefficient and the frequency division high frequency weight coefficient to obtain the frequency division phase control parameter model channel of the channel to be calculated;

[0080] The frequency division phase control parameter model data body establishing module 124 is used to obtain the frequency division phase control parameter model data body from the frequency division phase control parameter model channels of each channel to be calculated.

[0081] In some embodiments, the data acquisition module 121 is used to extract the divided-frequency low-frequency seismic traces, divided-frequency medium-frequency seismic traces, divided-frequency high-frequency seismic traces, divided-frequency low-frequency model traces, divided-frequency medium-frequency model traces, divided-frequency high-frequency model traces and low-frequency parameter model traces of each trace to be calculated in the study area according to the method of constructing a grid and synthesizing records.

[0082] In some embodiments, the frequency division weight coefficient determination module 122 is used to calculate the sum of the squares of the sampling points in the frequency division low-frequency seismic channel, determine the ratio of the square to the number of sampling points in the frequency division low-frequency seismic channel, and determine the square root of the ratio as the root mean square amplitude of the frequency division low-frequency seismic channel; calculate the sum of the squares of the sampling points in the frequency division medium-frequency seismic channel, determine the ratio of the square to the number of sampling points in the frequency division medium-frequency seismic channel, and determine the square root of the ratio as the root mean square amplitude of the frequency division medium-frequency seismic channel; calculate the sum of the squares of the sampling points in the frequency division high-frequency seismic channel, determine the ratio of the square to the number of sampling points in the frequency division high-frequency seismic channel, and determine the square root of the ratio as the root mean square amplitude of the frequency division high-frequency seismic channel.

[0083] In some embodiments, the frequency division weight coefficient determination module 122 is used to determine the sum of the obtained root mean square amplitudes, and determine the ratio of the root mean square amplitude of the frequency division low-frequency seismic channel to the sum as the frequency division low-frequency weight coefficient, determine the ratio of the root mean square amplitude of the frequency division medium-frequency seismic channel to the sum as the frequency division medium-frequency weight coefficient, and determine the ratio of the root mean square amplitude of the frequency division high-frequency seismic channel to the sum as the frequency division high-frequency weight coefficient.

[0084] In some embodiments, the frequency division phase control parameter model channel establishing module 123 is used to combine the frequency division low frequency model channel, the frequency division intermediate frequency model channel, the frequency division high frequency model channel and the low frequency parameter model channel to obtain the frequency division phase control parameter model channel of the channel to be calculated using the following formula:

[0085] Z=Cl+Ql*Ml+Qm*Mm+Qh*Mh

[0086] Among them, Z is the frequency division phase control parameter model channel, Cl is the low-frequency parameter model channel, Ql is the frequency division low-frequency weight coefficient, Ml is the frequency division low-frequency model channel, Qm is the frequency division intermediate frequency weight coefficient, Mm is the frequency division intermediate frequency model channel, Qh is the frequency division high-frequency weight coefficient, and Mh is the frequency division high-frequency model channel.

[0087] In some embodiments, the above apparatus further includes a reservoir prediction module 125, which is used to:

[0088] Reservoir prediction is performed on the target layer in the study area according to the frequency-divided phase-controlled parameter model data body.

[0089] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0090] Based on the inventive concept of the present invention, an embodiment of the present invention further provides a computer storage medium, wherein the computer storage medium stores computer executable instructions, and when the computer executable instructions are executed by a processor, the above-mentioned method for establishing a frequency division phase control parameter model is implemented.

[0091] Based on the inventive concept of the present invention, an embodiment of the present invention further provides a server, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned method for establishing a frequency division phase control parameter model when executing the program.

[0092] Unless otherwise specifically stated, terms such as processing, computing, calculating, determining, displaying, etc. may refer to the actions and / or processes of one or more processing or computing systems, or similar devices, which operate and convert data represented as physical (e.g., electronic) quantities within registers or memories of a processing system into other data similarly represented as physical quantities within memories, registers, or other such information storage, transmission, or display devices of the processing system. Information and signals may be represented using any of a variety of different techniques and methods. For example, data, instructions, commands, information, signals, bits, symbols, and chips mentioned throughout the above description may be represented by voltages, currents, electromagnetic waves, magnetic fields or particles, light fields or particles, or any combination thereof.

[0093] It should be understood that the specific order or hierarchy of steps in the disclosed process is an example of an exemplary method. Based on design preferences, it should be understood that the specific order or hierarchy of steps in the process can be rearranged without departing from the scope of protection of the present disclosure. The attached method claims present the elements of the various steps in an exemplary order and are not intended to be limited to the specific order or hierarchy described.

[0094] In the above detailed description, various features are grouped together in a single embodiment to simplify the disclosure. This method of disclosure should not be interpreted as reflecting an intention that an embodiment of the claimed subject matter requires more features than those set forth in each claim. On the contrary, as reflected in the appended claims, the invention is in a state of having less than all the features of the disclosed individual embodiments. Therefore, the appended claims are hereby expressly incorporated into the detailed description, with each claim standing on its own as a separate preferred embodiment of the invention.

[0095] Those skilled in the art will also appreciate that the various illustrative logic blocks, modules, circuits, and algorithmic steps described in conjunction with the embodiments herein can all be implemented as electronic hardware, computer software, or a combination thereof. In order to clearly illustrate the interchangeability between hardware and software, various illustrative components, blocks, modules, circuits, and steps are generally described above around their functions. Whether such functions are implemented as hardware or software depends on specific applications and the design constraints imposed on the entire system. A skilled person can implement the described functions in an alternative manner for each specific application, but such implementation decisions should not be interpreted as departing from the scope of protection of the present disclosure.

[0096] The steps of the method or algorithm described in conjunction with the embodiments herein may be directly embodied as hardware, a software module executed by a processor, or a combination thereof. The software module may be located in a RAM memory, a flash memory, a ROM memory, an EPROM memory, an EEPROM memory, a register, a hard disk, a mobile disk, a CD-ROM, or any other form of storage medium well known in the art. An exemplary storage medium is connected to the processor so that the processor can read information from the storage medium and can write information to the storage medium. Of course, the storage medium may also be an integral part of the processor. The processor and the storage medium may be located in an ASIC. The ASIC may be located in a user terminal. Of course, the processor and the storage medium may also be present in a user terminal as discrete components.

[0097] For software implementation, the techniques described in this application can be implemented with modules (e.g., procedures, functions, etc.) that perform the functions described in this application. These software codes can be stored in a memory unit and executed by a processor. The memory unit can be implemented within the processor or outside the processor. In the latter case, it is coupled to the processor in a communication manner via various means, which are well known in the art.

[0098] The above description includes examples of one or more embodiments. Of course, it is impossible to describe all possible combinations of components or methods for the purpose of describing the above embodiments, but it should be recognized by those skilled in the art that the various embodiments may be further combined and arranged. Therefore, the embodiments described herein are intended to cover all such changes, modifications and variations that fall within the scope of protection of the appended claims. In addition, with respect to the term "comprising" used in the specification or claims, the word is covered in a manner similar to the term "including", just as "including," is explained as a transitional word in the claims. In addition, any term "or" used in the specification of the claims is intended to mean "non-exclusive or".

Claims

1. A method for establishing a frequency division phase control parameter model, characterized in that: include: Obtain the divided low-frequency seismic trace, divided medium-frequency seismic trace, divided high-frequency seismic trace, divided low-frequency model trace, divided medium-frequency model trace, divided high-frequency model trace and low-frequency parameter model trace of each trace to be calculated in the study area; Calculating the root mean square amplitude of the divided-frequency low-frequency seismic trace, the root mean square amplitude of the divided-frequency intermediate-frequency seismic trace, and the root mean square amplitude of the divided-frequency high-frequency seismic trace, and determining the divided-frequency low-frequency weight coefficient, the divided-frequency intermediate-frequency weight coefficient, and the divided-frequency high-frequency weight coefficient according to the obtained root mean square amplitudes; By using the frequency division low frequency weight coefficient, the frequency division intermediate frequency weight coefficient and the frequency division high frequency weight coefficient, the frequency division low frequency model channel, the frequency division intermediate frequency model channel, the frequency division high frequency model channel and the low frequency parameter model channel are merged to obtain the frequency division phase control parameter model channel of the channel to be calculated; A frequency-division phase-control parameter model data body is obtained from the frequency-division phase-control parameter model channels of each channel to be calculated.

2. The method according to claim 1, characterized in that The method of obtaining the divided-frequency low-frequency seismic trace, divided-frequency medium-frequency seismic trace, divided-frequency high-frequency seismic trace, divided-frequency low-frequency model trace, divided-frequency medium-frequency model trace, divided-frequency high-frequency model trace and low-frequency parameter model trace of each trace to be calculated in the study area includes: According to the structural grid and synthetic record method, the divided low-frequency seismic traces, divided medium-frequency seismic traces, divided high-frequency seismic traces, divided low-frequency model traces, divided medium-frequency model traces, divided high-frequency model traces and low-frequency parameter model traces of each trace to be calculated in the study area are extracted.

3. The method according to claim 1, characterized in that The calculating of the root mean square amplitude of the divided low-frequency seismic trace, the root mean square amplitude of the divided medium-frequency seismic trace and the root mean square amplitude of the divided high-frequency seismic trace comprises: Calculating the sum of squares of the sampling points in the frequency-divided low-frequency seismic trace, determining the ratio of the square to the number of sampling points in the frequency-divided low-frequency seismic trace, and determining the square root of the ratio as the root mean square amplitude of the frequency-divided low-frequency seismic trace; Calculating the sum of squares of the sampling points in the divided-frequency intermediate frequency seismic trace, determining the ratio of the square to the number of sampling points in the divided-frequency intermediate frequency seismic trace, and determining the square root of the ratio as the root mean square amplitude of the divided-frequency intermediate frequency seismic trace; The sum of squares of the sampling points in the frequency-divided high-frequency seismic trace is calculated, the ratio of the square to the number of sampling points in the frequency-divided high-frequency seismic trace is determined, and the square root of the ratio is determined as the root mean square amplitude of the frequency-divided high-frequency seismic trace.

4. The method according to claim 1, characterized in that: The step of determining the frequency division low frequency weight coefficient, the frequency division intermediate frequency weight coefficient and the frequency division high frequency weight coefficient according to the obtained root mean square amplitudes includes: Determine the sum of the obtained root mean square amplitudes, determine the ratio of the root mean square amplitude of the divided-frequency low-frequency seismic channel to the sum as the divided-frequency low-frequency weight coefficient, determine the ratio of the root mean square amplitude of the divided-frequency intermediate-frequency seismic channel to the sum as the divided-frequency intermediate-frequency weight coefficient, and determine the ratio of the root mean square amplitude of the divided-frequency high-frequency seismic channel to the sum as the divided-frequency high-frequency weight coefficient.

5. The method according to claim 1, characterized in that The step of merging the divided-frequency low-frequency model channel, the divided-frequency intermediate-frequency model channel, the divided-frequency high-frequency model channel and the low-frequency parameter model channel to obtain the divided-frequency phase-controlled parameter model channel of the channel to be calculated includes: The frequency-divided low-frequency model channel, the frequency-divided intermediate-frequency model channel, the frequency-divided high-frequency model channel and the low-frequency parameter model channel are combined using the following formula to obtain the frequency-divided phase-controlled parameter model channel of the channel to be calculated: Z=Cl+Ql*Ml+Qm*Mm+Qh*Mh Among them, Z is the frequency division phase control parameter model channel, Cl is the low-frequency parameter model channel, Ql is the frequency division low-frequency weight coefficient, Ml is the frequency division low-frequency model channel, Qm is the frequency division intermediate frequency weight coefficient, Mm is the frequency division intermediate frequency model channel, Qh is the frequency division high-frequency weight coefficient, and Mh is the frequency division high-frequency model channel.

6. The method according to any one of claims 1 to 5, characterized in that: Also includes: Reservoir prediction is performed on the target layer in the study area according to the frequency-divided phase-controlled parameter model data body.

7. A device for establishing a frequency division phase control parameter model, characterized in that: The device comprises: The data acquisition module is used to obtain the divided low-frequency seismic traces, divided medium-frequency seismic traces, divided high-frequency seismic traces, divided low-frequency model traces, divided medium-frequency model traces, divided high-frequency model traces and low-frequency parameter model traces of each trace to be calculated in the study area; A frequency division weight coefficient determination module is used to calculate the root mean square amplitude of the frequency division low-frequency seismic trace, the root mean square amplitude of the frequency division medium-frequency seismic trace and the root mean square amplitude of the frequency division high-frequency seismic trace, and determine the frequency division low-frequency weight coefficient, the frequency division medium-frequency weight coefficient and the frequency division high-frequency weight coefficient according to the obtained root mean square amplitudes; A frequency division phase control parameter model channel establishment module is used to combine the frequency division low frequency model channel, the frequency division medium frequency model channel, the frequency division high frequency model channel and the low frequency parameter model channel by using the frequency division low frequency weight coefficient, the frequency division medium frequency weight coefficient and the frequency division high frequency weight coefficient to obtain the frequency division phase control parameter model channel of the channel to be calculated; The frequency division phase control parameter model data body establishment module is used to obtain the frequency division phase control parameter model data body from the frequency division phase control parameter model channels of each channel to be calculated.

8. The device according to claim 7, characterized in that The frequency division weight coefficient determination module determines the frequency division low frequency weight coefficient, the frequency division intermediate frequency weight coefficient and the frequency division high frequency weight coefficient according to the obtained root mean square amplitudes, and is used for: Determine the sum of the obtained root mean square amplitudes, determine the ratio of the root mean square amplitude of the divided-frequency low-frequency seismic channel to the sum as the divided-frequency low-frequency weight coefficient, determine the ratio of the root mean square amplitude of the divided-frequency intermediate-frequency seismic channel to the sum as the divided-frequency intermediate-frequency weight coefficient, and determine the ratio of the root mean square amplitude of the divided-frequency high-frequency seismic channel to the sum as the divided-frequency high-frequency weight coefficient.

9. The device according to claim 7, characterized in that The frequency division phase control parameter model channel establishment module combines the frequency division low frequency model channel, the frequency division intermediate frequency model channel, the frequency division high frequency model channel and the low frequency parameter model channel to obtain the frequency division phase control parameter model channel of the channel to be calculated, which is used for: The frequency-divided low-frequency model channel, the frequency-divided intermediate-frequency model channel, the frequency-divided high-frequency model channel and the low-frequency parameter model channel are combined using the following formula to obtain the frequency-divided phase-controlled parameter model channel of the channel to be calculated: Z=Cl+Ql*Ml+Qm*Mm+Qh*Mh Among them, Z is the frequency division phase control parameter model channel, Cl is the low-frequency parameter model channel, Ql is the frequency division low-frequency weight coefficient, Ml is the frequency division low-frequency model channel, Qm is the frequency division intermediate frequency weight coefficient, Mm is the frequency division intermediate frequency model channel, Qh is the frequency division high-frequency weight coefficient, and Mh is the frequency division high-frequency model channel.

10. A computer storage medium, characterized in that: The computer storage medium stores computer executable instructions, and when the computer executable instructions are executed by a processor, the method for establishing a frequency division phase control parameter model according to any one of claims 1 to 6 is implemented.

11. A server, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the method for establishing a frequency division phase control parameter model as claimed in any one of claims 1 to 6 when executing the program.

Images