A method and related device for predicting the distribution range of gas-bearing sand bodies in narrow river channels

Abstract

The present application discloses a method and related device for predicting the distribution range of gas-bearing sand bodies in narrow channels, and relates to the field of petroleum geological exploration technology. The method comprises: obtaining seismic exploration data, calculating the instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer according to the seismic exploration data, and calculating the characteristic index of gas-bearing sand bodies in each layer accordingly, and finally determining the distribution range of gas-bearing sand bodies according to the characteristic index of gas-bearing sand bodies in each layer and the threshold value of characteristic index of gas-bearing sand bodies. When calculating the characteristic index of gas-bearing sand bodies, the present application takes into account the thin layer and mesh characteristics of narrow channels and the influence on seismic energy absorption and attenuation under gas-bearing conditions, and can more accurately predict the characteristics of gas-bearing sand bodies in narrow channels, and improve the accuracy of gas-bearing sand body distribution estimation.

Description

Technical Field

[0001] The present application relates to the technical field of petroleum geological exploration, and in particular to a method for predicting the distribution range of gas-bearing sand bodies in narrow channels and related devices. Background Art

[0002] Narrow channel oil and gas reservoirs have the characteristics of good storage and high production capacity. In recent years, they have received extensive attention from the industry. The distribution of narrow channel gas reservoirs has become a key target for oil and gas exploration. However, this type of oil and gas reservoir is a lithological oil and gas reservoir. Inaccurate identification will result in drilling failure or low drilling rate, and high geological risk of exploration. Radovich (1998) first proposed the concept of oil and gas sweet spots. When they studied the Gorgon oil field in Australia, they found that the gas-bearing sand bodies had strong instantaneous amplitude and low instantaneous frequency. Hart (2008) proposed to use the reflection intensity (also known as "instantaneous amplitude") divided by the square root of the instantaneous frequency to represent the seismic attributes of the sweet spot. However, these methods are more practical in high-porosity and permeability sandstones, but it is difficult to accurately determine the distribution range of gas-bearing sand bodies in medium- and low-porosity and permeability reservoirs with fast lateral changes in narrow channels, which increases the difficulty of geological identification and risk assessment of such oil and gas reservoirs, resulting in exploration well failures or inaccurate reserve estimates. In addition, although the current seismic data contains rich information on the reflection structure characteristics of river channel sand bodies and fluids, the use of information in seismic data in oil and gas exploration is far from enough, and the existing technical methods do not use seismic data to directly predict the distribution range of gas-bearing narrow channel sand bodies. Therefore, how to quickly and effectively identify and predict the boundaries of gas-bearing narrow channel sand bodies is an urgent problem that needs to be solved in exploration under similar geological conditions. Summary of the invention

[0003] The purpose of the present application is to provide a method and related device for predicting the distribution range of gas-bearing sand bodies in narrow river channels, which can quickly and effectively identify and predict the distribution range of gas-bearing sand bodies in narrow river channels.

[0004] To achieve the above objectives, this application provides the following solutions:

[0005] In a first aspect, the present application provides a method for predicting the distribution range of a narrow channel gas-bearing sand body, comprising the following steps:

[0006] Obtain seismic exploration data; seismic exploration data include well logging data of different layers, petroleum seismic data and their interpretation results data.

[0007] The instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer segment are calculated according to the seismic exploration data.

[0008] According to the instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer, the characteristic index of gas-bearing sand body of each layer is calculated respectively.

[0009] The distribution range of the gas-bearing sand body is determined according to the gas-bearing sand body characteristic index and the gas-bearing sand body characteristic index threshold of each layer.

[0010] Optionally, before determining the distribution range of the gas-bearing sand body according to the gas-bearing sand body characteristic index and the gas-bearing sand body characteristic index threshold of each layer, the distribution range prediction method of the narrow channel gas-bearing sand body further includes the following steps:

[0011] For the gas-bearing sand bodies that have been discovered through drilling in the region, the correlation between the characteristic parameters of the gas-bearing sand bodies and the characteristic index of the gas-bearing sand bodies is established to verify the accuracy of the characteristic index of the gas-bearing sand bodies.

[0012] After the accuracy of the gas-bearing sand body characteristic index is verified, an intersection diagram of the gas-bearing sand body characteristic index and the characteristic parameters of the gas-bearing sand body is established, and the threshold value of the gas-bearing sand body characteristic index is determined.

[0013] Optionally, calculate the instantaneous amplitude property according to the following formula:

[0014] .

[0015] in, E ( t ) is the instantaneous amplitude property, T ( t ) is the original earthquake record, H ( t ) is the imaginary part obtained by Hilbert transforming the original seismic signal.

[0016] Optionally, the instantaneous frequency property is calculated according to the following formula:

[0017] .

[0018] in, f ( t ) is the instantaneous frequency attribute, is the instantaneous phase of the original seismic signal, is the time derivative of the instantaneous phase.

[0019] Optionally, a weighted variance attribute is calculated according to the following formula:

[0020] ;

[0021] in, is the weighted variance attribute, For the j The variance value attribute at each time point, , , ω j For the j The weight of the variance value attribute at each time point,l To obtain the time window length of the variance, j The variance value attribute of each time point It can be calculated by the following formula:

[0022] ;

[0023] in, x ij For the i Channels in j The earthquake amplitude value at time, w To find the number of data channels for the variance, For each channel j The average earthquake amplitude over time.

[0024] Optionally, the gas-bearing sand body characteristic index is calculated according to the following formula:

[0025] .

[0026] in, Nc is the characteristic index of gas-bearing sand body, E ( t ) is the instantaneous amplitude property, B is the regional adjustment coefficient, f ( t ) is the instantaneous frequency attribute, A is the instantaneous frequency deviation coefficient, is the weighted variance attribute.

[0027] In a second aspect, the present application provides a distribution range prediction device for narrow channel gas-bearing sand bodies, comprising the following modules:

[0028] The exploration data acquisition module is used to acquire seismic exploration data; seismic exploration data includes well logging data of different layers, petroleum seismic data and its interpretation results data.

[0029] The layer segment attribute calculation module is used to calculate the instantaneous amplitude attribute, instantaneous frequency attribute and variance value attribute of each layer segment according to the seismic exploration data.

[0030] The characteristic index calculation module is used to calculate the characteristic index of the gas-bearing sand body of each layer according to the instantaneous amplitude attribute, instantaneous frequency attribute and variance value attribute of each layer.

[0031] The gas-bearing sand body identification module is used to determine the distribution range of the gas-bearing sand body according to the gas-bearing sand body characteristic index and the gas-bearing sand body characteristic index threshold of each layer.

[0032] In a third aspect, the present application provides a computer device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the method for predicting the distribution range of narrow channel gas-bearing sand bodies described above.

[0033] In a fourth aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method for predicting the distribution range of narrow channel gas-bearing sand bodies as described above.

[0034] In a fifth aspect, the present application provides a computer program product, including a computer program, which, when executed by a processor, implements the steps of the method for predicting the distribution range of narrow channel gas-bearing sand bodies as described above.

[0035] According to the specific embodiments provided in this application, this application discloses the following technical effects:

[0036] The present application provides a method and related device for predicting the distribution range of gas-bearing sand bodies in narrow river channels. In the method, the instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer are calculated based on the acquired seismic exploration data, and the characteristic index of the gas-bearing sand bodies in each layer is calculated respectively based on this, and finally the distribution range of the gas-bearing sand bodies is determined according to the characteristic index of the gas-bearing sand bodies in each layer and the threshold value of the characteristic index of the gas-bearing sand bodies. The above scheme of the present application calculates the characteristic index of the gas-bearing sand bodies according to the instantaneous amplitude attribute, the instantaneous frequency attribute and the weighted variance value attribute, and takes into account the thin layer and mesh characteristics of the narrow river channel and the influence on the seismic energy absorption attenuation under gas-bearing conditions, so as to more accurately predict the characteristics of the gas-bearing sand bodies in the narrow river channel, and improve the accuracy of the estimation of the distribution of the gas-bearing sand bodies. The method can be implemented by an automated device, and can automatically complete data input, calculation and result output, thereby improving the work efficiency of the prediction of the distribution of gas-bearing sand bodies in narrow river channels. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0038] Figure 1 A schematic flow chart of a method for predicting the distribution range of gas-bearing sand bodies in narrow channels provided in one embodiment of the present application.

[0039] Figure 2A schematic diagram of interpreting the sedimentary subfacies characteristics of narrow channel sand bodies based on Ncg values ​​in a method for predicting the distribution range of narrow channel gas-bearing sand bodies provided in one embodiment of the present application.

[0040] Figure 3 A schematic diagram of the distribution range of narrow channel gas-bearing sand bodies generated in a method for predicting the distribution range of narrow channel gas-bearing sand bodies provided in one embodiment of the present application.

[0041] Figure 4 A flowchart of determining a threshold value of a characteristic index of a gas-bearing sand body in a method for predicting the distribution range of a narrow channel gas-bearing sand body provided in one embodiment of the present application.

[0042] Figure 5 An intersection diagram of a gas-bearing sand body characteristic index and water saturation in a method for predicting the distribution range of a narrow channel gas-bearing sand body provided in one embodiment of the present application.

[0043] Figure 6 A schematic diagram of the functional modules of a distribution range prediction device for narrow channel gas-bearing sand bodies provided in one embodiment of the present application.

[0044] Figure 7 A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0045] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0046] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0047] The present application provides a method for predicting the distribution range of a narrow channel gas-bearing sand body, such as Figure 1 As shown, the following steps are included:

[0048] A1. Obtain seismic exploration data; seismic exploration data include well logging data of different layers, petroleum seismic data and their interpretation results data. Generally speaking, well logging data include conventional well logging curves and layer data, petroleum seismic data include three-dimensional post-stack seismic data, and interpretation results data include layer interpretation and fault interpretation data.

[0049] A2. Calculate the instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer segment based on the seismic exploration data.

[0050] The instantaneous amplitude attribute is a key seismic attribute that describes the strength of seismic signals. It provides a measure of the overall energy of the signal at any given moment and is often used to highlight features such as changes in reflection intensity. In narrow channel sandstone formations, the lithology has little change in the reflection intensity of seismic waves. The instantaneous amplitude can help identify areas with weak reflection intensity but relatively stable distribution, which are usually related to dense rock formations or low permeability layers. That is, the higher the instantaneous amplitude intensity, the higher the sandstone content. The instantaneous amplitude attribute can be calculated according to the following formula:

[0051] .

[0052] in, E ( t ) is the instantaneous amplitude property, T ( t ) is the original earthquake record, H ( t ) is the imaginary part obtained by Hilbert transforming the original seismic signal.

[0053] In this embodiment, the instantaneous amplitude attribute is first obtained by calculation, and the instantaneous amplitude attribute data of the target layer segment is extracted using the target layer position, and the amplitude intensity data ranges from 0 to 35000.

[0054] The instantaneous frequency attribute is the frequency of the seismic signal at a certain moment. It is used to describe the change of the signal frequency over time. It is usually calculated by the analytical signal obtained by Hilbert transforming the seismic signal. The instantaneous frequency attribute can be calculated according to the following formula:

[0055] .

[0056] in, f ( t ) is the instantaneous frequency attribute, is the instantaneous phase of the original seismic signal, is the derivative of the instantaneous phase with respect to time, that is, the rate of change of phase with time.

[0057] In this embodiment, a time window of 33 ms is used, and the calculation result obtains the instantaneous frequency attribute, and the data range is -154 to 532 Hz.

[0058] The weighted variance attribute refers to the degree of spatial variation of adjacent seismic traces. The larger the weighted variance value, the worse the lateral continuity of the seismic wave, which can indicate the point of change caused by abnormal geological bodies. In terms of river channel identification, the weighted variance attribute can highlight the distribution characteristics and internal structure of the river channel. The weighted variance attribute can be calculated according to the following formula:

[0059] .

[0060] in, is the weighted variance attribute, For the j The variance value attribute at each time point, , , ω j For the j The weight of the variance value attribute at each time point, l To obtain the time window length of the variance, j The variance value attribute of each time point It can be calculated by the following formula:

[0061] .

[0062] in, x ij For the i Channels in j The earthquake amplitude value at time, w To find the number of data channels for the variance, For each channel j The average earthquake amplitude over time.

[0063] Specifically in this embodiment, the data range of the weighted variance value attribute obtained by the above calculation is between 0 and 1, and the smaller the value, the better the continuity, and vice versa.

[0064] A3. According to the instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer segment, the characteristic index of the gas-bearing sand body of each layer segment is calculated respectively.

[0065] Considering that seismic waves have an attenuation effect when they propagate underground, they are closely related to the underground lithofacies, porosity, and fluid properties. From the perspective of petroleum seismic exploration data research, the absorption of seismic waves by gas-bearing sand layers is related to the frequency of seismic waves, and increases with the increase of seismic wave frequency. Relatively high-porosity reservoirs have larger amplitudes than surrounding rocks, and their instantaneous frequencies are relatively small. When oil and gas are contained, the instantaneous amplitude properties are E ( t ) becomes larger, the instantaneous frequency property f ( t ) becomes smaller, by eliminating the negative value of instantaneous frequency and transforming the instantaneous frequency, the rate of change of low frequency is highlighted. After the two attributes are divided, the abnormal feature is further amplified and the gas-bearing feature is highlighted. In this embodiment, the gas-bearing sand body characteristic index is calculated according to the following formula:

[0066] .

[0067] in, Nc is the characteristic index of gas-bearing sand body, E ( t) is the instantaneous amplitude property, B is the regional adjustment coefficient, f ( t ) is the instantaneous frequency attribute, A is the instantaneous frequency deviation coefficient, is the weighted variance attribute.

[0068] Comprehensively analyze the instantaneous frequency attribute range, and set the instantaneous frequency offset coefficient to eliminate the influence of negative frequency. A is 200 Hz, the regional regulation coefficient B The gas-bearing sand body characteristic index is calculated to be 1E+6. Nc , which can be used for attribute extraction, time slicing and profile analysis. Its main function is to show which part of the plane is a high-value area and which part is a low-value area, such as Figure 2 As shown, the sedimentary subfacies characteristics (such as gas-bearing main channel, breach fan, floodplain, etc.) and gas-bearing distribution of narrow channel sand bodies are explained. The specific operation steps are as follows:

[0069] Layer interpolation: Perform layer interpolation on the layers to ensure the continuity and accuracy of the data.

[0070] Select the attribute extraction method: use one or two layers as the base surface, determine the time window by the top and bottom layers, or set the time window according to the earthquake travel time range.

[0071] Attribute extraction: Display gas-bearing sand body characteristic index on horizontal time slices or vertical sections Nc The value in the corresponding time window.

[0072] A4. According to the gas-bearing sand characteristic index and gas-bearing sand characteristic index threshold of each layer, the distribution range of the gas-bearing sand body is determined. Specifically, the calculated gas-bearing sand characteristic index within the range of a narrow river channel sediment to be predicted is Nc The data is compared with the threshold of gas-bearing sand body characteristic index. If it is greater than the threshold, the distribution range of gas-bearing sand body can be quantitatively determined, such as Figure 3 shown.

[0073] Specifically in this embodiment, before executing step A4, the method for predicting the distribution range of narrow channel gas-bearing sand bodies also includes a process for determining a threshold value of a characteristic index of a gas-bearing sand body, such as Figure 4 As shown, the process includes the following steps:

[0074] B1. For the gas-bearing sand bodies that have been discovered by drilling in the region, the correlation between the characteristic parameters of the gas-bearing sand bodies and the characteristic index of the gas-bearing sand bodies is constructed to verify the accuracy of the characteristic index of the gas-bearing sand bodies. In this embodiment, the porosity, permeability, water saturation, mud content, instantaneous amplitude, instantaneous frequency, variance characteristics, and characteristic index of the gas-bearing sand bodies in the narrow river channel that have been discovered by drilling in the region are statistically analyzed to establish an intersection relationship and obtain Table 1. Among the characteristic parameters of the narrow river channel gas-bearing sand bodies, the value 1 is used to subtract each characteristic parameter. It can be seen that the water saturation mainly reflects the gas-bearing property, while the other characteristic parameters reflect the lithology or physical properties. Therefore, the water saturation is selected to evaluate the characteristic index of the gas-bearing sand body. Nc The ability to identify gas-bearing sand body characteristics through the gas-bearing sand body characteristic index Nc The relationship between the gas-bearing sand body characteristic index and water saturation Nc It has the ability to distinguish gas-bearing sand bodies from non-gas-bearing sand bodies, verifying the gas-bearing sand body characteristic index calculated in this application Nc accuracy.

[0075] Table 1 Correlation between characteristic parameters of gas-bearing sand bodies and characteristic indexes of gas-bearing sand bodies

[0076]

[0077] B2. After the accuracy of the gas-bearing sand body characteristic index is verified, an intersection diagram of the gas-bearing sand body characteristic index and the characteristic parameters of the gas-bearing sand body is established, and a threshold value of the gas-bearing sand body characteristic index is determined. In this embodiment, the intersection diagram of the gas-bearing sand body characteristic index and the water saturation is as follows: Figure 5 As shown, in Figure 5 There are two main clusters: one is the gas-containing cluster, Nc >0.9, water saturation is between 16% and 65%, and the other is a cluster without gas. Nc <0.6, water saturation>56%, thus it can be determined that the threshold value of the gas sand characteristic index of the narrow channel gas sand should be 0.9.

[0078] The above embodiment of the present application provides a method for predicting the distribution range of gas-bearing sand bodies in narrow channels, which is suitable for geological conditions with low porosity and low permeability in narrow channels, and has high practical value and application prospects. Based on the same inventive concept, the embodiment of the present application also provides a device for implementing the above-mentioned method for predicting the distribution range of gas-bearing sand bodies in narrow channels. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme recorded in the above method, so the specific limitations in one or more device embodiments provided below can refer to the above limitations on the method for predicting the distribution range of gas-bearing sand bodies in narrow channels, and will not be repeated here.

[0079] In an exemplary embodiment, Figure 6As shown, a distribution range prediction device for narrow channel gas-bearing sand bodies is provided, comprising the following modules:

[0080] The exploration data acquisition module is used to acquire seismic exploration data; seismic exploration data includes well logging data of different layers, petroleum seismic data and its interpretation results data.

[0081] The layer segment attribute calculation module is used to calculate the instantaneous amplitude attribute, instantaneous frequency attribute and variance value attribute of each layer segment according to the seismic exploration data.

[0082] The characteristic index calculation module is used to calculate the characteristic index of the gas-bearing sand body of each layer according to the instantaneous amplitude attribute, instantaneous frequency attribute and variance value attribute of each layer.

[0083] The gas-bearing sand body identification module is used to determine the distribution range of the gas-bearing sand body according to the gas-bearing sand body characteristic index and the gas-bearing sand body characteristic index threshold of each layer.

[0084] certainly, Figure 6 The architecture shown is only exemplary and can be omitted according to actual needs when implementing different functions. Figure 6 One or at least two components of the system shown.

[0085] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Figure 7 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. Among them, the processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store seismic exploration data and gas sand body characteristic index thresholds. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, the steps of the distribution range prediction method of the narrow channel gas sand body or the functions of the distribution range prediction device of the narrow channel gas sand body provided in the above embodiment can be implemented.

[0086] Those skilled in the art will understand that Figure 7The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0087] In an exemplary embodiment, a computer device is further provided, including a memory and a processor, wherein a computer program is stored in the memory, and the processor implements the steps in the above-mentioned method embodiments when executing the computer program.

[0088] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0089] In an exemplary embodiment, a computer program product is provided, including a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0090] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.

[0091] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0092] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but is not limited thereto.

[0093] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0094] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A method for predicting the distribution range of gas-bearing sand bodies in narrow channels, characterized in that: include: Acquiring seismic exploration data; the seismic exploration data includes well logging data of different layers, petroleum seismic data and interpretation results data; Calculating instantaneous amplitude attributes, instantaneous frequency attributes and weighted variance value attributes of each layer segment according to the seismic exploration data; According to the instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer, the characteristic index of gas-bearing sand body of each layer is calculated respectively; According to the gas-bearing sand body characteristic index and gas-bearing sand body characteristic index threshold of each layer, the distribution range of the gas-bearing sand body is determined; The gas-bearing sand body characteristic index is calculated according to the following formula: ； in, Nc is the characteristic index of gas-bearing sand body, E ( t ) is the instantaneous amplitude property, B is the regional adjustment coefficient, f ( t ) is the instantaneous frequency attribute, A is the instantaneous frequency deviation coefficient, is the weighted variance attribute; The instantaneous amplitude property is calculated according to the following formula: ； in, E ( t ) is the instantaneous amplitude property, T ( t ) is the original earthquake record, H ( t ) is the imaginary part obtained by Hilbert transforming the original seismic signal; The instantaneous frequency property is calculated according to the following formula: ； in, f ( t ) is the instantaneous frequency attribute, is the instantaneous phase of the original seismic signal, is the derivative of the instantaneous phase with respect to time; The weighted variance value attribute is calculated according to the following formula: ； in, is the weighted variance attribute, For the j The variance value attribute at each time point, , , ω j For the j The weight of the variance value attribute at each time point, l To obtain the time window length of the variance, j The variance value attribute of each time point Calculated by the following formula: ； in, x ij For the i Channels in j The earthquake amplitude value at time, w To find the number of data channels for the variance, For each channel j The average earthquake amplitude over time.

2. The method for predicting the distribution range of narrow channel gas-bearing sand bodies according to claim 1, characterized in that: Before determining the distribution range of the gas-bearing sand body according to the gas-bearing sand body characteristic index and the gas-bearing sand body characteristic index threshold of each layer, the distribution range prediction method of the narrow channel gas-bearing sand body further includes: For the gas-bearing sand bodies that have been discovered through drilling in the region, the correlation between the characteristic parameters of the gas-bearing sand bodies and the characteristic index of the gas-bearing sand bodies is constructed to verify the accuracy of the characteristic index of the gas-bearing sand bodies; After the accuracy of the gas-bearing sand body characteristic index is verified, an intersection diagram of the gas-bearing sand body characteristic index and the characteristic parameters of the gas-bearing sand body is established, and a threshold value of the gas-bearing sand body characteristic index is determined.

3. A distribution range prediction device for narrow channel gas-bearing sand bodies, characterized in that: include: An exploration data acquisition module is used to acquire seismic exploration data; the seismic exploration data includes well logging data of different layers, petroleum seismic data and interpretation results data; A layer segment attribute calculation module, used for calculating the instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer segment according to the seismic exploration data; The characteristic index calculation module is used to calculate the characteristic index of the gas-bearing sand body of each layer according to the instantaneous amplitude attribute, instantaneous frequency attribute and weighted variance value attribute of each layer; A gas-bearing sand body identification module is used to determine the distribution range of the gas-bearing sand body according to the gas-bearing sand body characteristic index and the gas-bearing sand body characteristic index threshold of each layer; The gas-bearing sand body characteristic index is calculated according to the following formula: ； in, Nc is the characteristic index of gas-bearing sand body, E ( t ) is the instantaneous amplitude property, B is the regional adjustment coefficient, f ( t ) is the instantaneous frequency attribute, A is the instantaneous frequency deviation coefficient, is the weighted variance attribute; The instantaneous amplitude property is calculated according to the following formula: ； in, E ( t ) is the instantaneous amplitude property, T ( t ) is the original earthquake record, H ( t ) is the imaginary part obtained by Hilbert transforming the original seismic signal; The instantaneous frequency property is calculated according to the following formula: ； in, f ( t ) is the instantaneous frequency attribute, is the instantaneous phase of the original seismic signal, is the derivative of the instantaneous phase with respect to time; The weighted variance value attribute is calculated according to the following formula: ； in, is the weighted variance attribute, For the j The variance value attribute at each time point, , , ω j For the j The weight of the variance value attribute at each time point, l To obtain the time window length of the variance, j The variance value attribute of each time point Calculated by the following formula: ； in, x ij For the i Channels in j The earthquake amplitude value at time, w To find the number of data channels for the variance, For each channel j The average earthquake amplitude over time.

4. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method for predicting the distribution range of narrow channel gas-bearing sand bodies according to any one of claims 1 to 2.

5. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the distribution range prediction method of the narrow channel gas-bearing sand body according to any one of claims 1 to 2 is implemented.

6. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the distribution range prediction method of the narrow channel gas-bearing sand body according to any one of claims 1 to 2 is implemented.

Images