Gas content detection method, device, equipment and storage medium

Abstract

The application provides a gas content detection method, device, equipment and storage medium, which comprises the following steps: collecting pre-stack angle domain gathers, determining an effective incident angle range and a frequency band range of the pre-stack angle domain gathers; performing angle stacking according to the incident angle range to obtain n angle stacked seismic data and a center angle of the corresponding incident angle stacking range; calculating the amplitude frequency spectrum of each trace of the n angle stacked seismic data to obtain an amplitude frequency spectrum area and an energy maximum value in the amplitude frequency spectrum in the frequency band range of the seismic data; applying the amplitude frequency spectrum area and the center angle of the incident angle stacking range to obtain a variation rate data body of the seismic data spectrum area in the incident angle direction; applying the energy maximum value in the amplitude frequency spectrum and the center angle of the incident angle stacking range to obtain a variation rate data body of the energy maximum value in the amplitude frequency spectrum in the incident angle direction; calibrating the data body and the data body with a known gas well in an actual work area to complete gas content prediction.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of oil exploration, in particular to a gas-bearing property detection method, device, equipment and storage medium. BACKGROUND

[0002] At present, most of the gas-bearing property detection methods are qualitative prediction of post-stack attenuation attributes, and quantitative prediction of pre-stack inversion. For qualitative prediction of gas-bearing property by post-stack attenuation attribute, only post-stack seismic data is used, which lacks the rich information of pre-stack seismic data. The quantitative prediction technology of gas-bearing property by pre-stack inversion is based on the relationship between elastic parameters and gas saturation to calculate the gas saturation, or through neural network algorithm to learn. Since the change of gas saturation in elastic parameters is very small, the effect of quantitative prediction of gas-bearing property by pre-stack is not ideal. When quantitative prediction of gas-bearing property is performed, the stability of gas saturation interpretation result of well logging interpretation is poor, so if the gas saturation data of well logging interpretation is referred in the prediction of gas-bearing property, the accuracy of the prediction result will also be reduced. SUMMARY

[0003] In view of the above problems, the present application provides a gas-bearing property detection method, device, equipment and storage medium.

[0004] The present application provides a gas-bearing property detection method, comprising:

[0005] S1: collecting pre-stack angle domain gathers, determining an effective incident angle range and determining a frequency band range of the pre-stack angle domain gathers which are effective;

[0006] S2: performing angle stacking according to the incident angle range to obtain n angle-stacked seismic data, and the center angles of the corresponding incident angle stacking ranges are angle1, angle2, angle3…anglen respectively;

[0007] S3: calculating the amplitude frequency spectrum of each trace of the n angle-stacked seismic data to obtain the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data;

[0008] S4: applying the amplitude frequency spectrum area area and the center angle of the incident angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction;

[0009] S5: applying the maximum energy max_amp in the amplitude frequency spectrum and the center angle of the incident angle stacking range to obtain the variation rate data body max_amp_rate of the maximum energy in the amplitude frequency spectrum in the incident angle direction;

[0010] S6: calibrate the data body area_rate and data body max_amp_rate with the actual work area known gas well, complete the gas-bearing property prediction.

[0011] In some embodiments, before the determining the effective frequency band range of the pre-stack angle domain gather, the method further comprises:

[0012] Performing residual time correction.

[0013] In some embodiments, the n is not greater than 5.

[0014] In some embodiments, the specific method of applying the amplitude frequency spectrum area area and the center angle of the incident angle stack range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction comprises:

[0015] Let the Y axis be the n amplitude frequency spectrum areas area, and the X axis be the center angles angle1, angle2, angle3…anglen of the n incident angle stack ranges. Through curve fitting, the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction is obtained.

[0016] In some embodiments, the specific method of applying the maximum energy max_amp in the amplitude frequency spectrum and the center angle of the incident angle stack range to obtain the variation rate data body max_amp_rate of the maximum energy in the amplitude frequency spectrum in the incident angle direction comprises:

[0017] Let the Y axis be the n maximum energies max_amp in the amplitude frequency spectrum, and the X axis be the center angles angle1, angle2, angle3…anglen of the n incident angle stack ranges. Through curve fitting, the variation rate data body max_amp_rate of the maximum energy in the amplitude frequency spectrum in the incident angle direction is obtained.

[0018] In some embodiments, the effective incident angle range is 0-36.

[0019] In some embodiments, the effective frequency band range of the pre-stack angle domain gather is 7-75 Hz.

[0020] Embodiments of the present application provide a gas-bearing property detection device, comprising:

[0021] The pre-stack angle domain gather collection module, the center angle calculation module, the amplitude frequency spectrum data calculation module, the data body area_rate calculation module, the data body max_amp_rate calculation module, and the gas-bearing property prediction module.

[0022] A pre-stack angle domain gather collection module: collecting pre-stack angle domain gathers, determining an effective incident angle range, and determining a frequency band range in which the pre-stack angle domain gathers are effective;

[0023] A central angle calculation module: performing angle stacking according to the incident angle range, obtaining n angle-stacked seismic data, and corresponding central angles of the incident angle stacking range are angle1, angle2, angle3…anglen respectively;

[0024] An amplitude frequency spectrum data calculation module: calculating the amplitude frequency spectrum of each trace of the n angle-stacked seismic data, and obtaining an amplitude frequency spectrum area area and an energy maximum value max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data;

[0025] A data body area_rate calculation module: applying the amplitude frequency spectrum area area and the central angle of the incident angle stacking range, and obtaining a data body area_rate of the variation rate of the seismic data spectrum area in the incident angle direction;

[0026] A data body max_amp_rate calculation module: applying the energy maximum value max_amp in the amplitude frequency spectrum and the central angle of the incident angle stacking range, and obtaining a data body max_amp_rate of the variation rate of the energy maximum value in the amplitude frequency spectrum in the incident angle direction;

[0027] A gas-bearing property prediction module: calibrating the data body area_rate and the data body max_amp_rate with known gas production wells in an actual work area, and completing gas-bearing property prediction.

[0028] The embodiment of the application provides a gas-bearing property detection device, comprising a memory and a processor, and the memory stores a computer program.

[0029] The embodiment of the application provides a storage medium, which stores a computer program capable of being executed by one or more processors and capable of being used to implement the gas-bearing property detection method.

[0030] The gas-bearing property detection method, device, equipment and storage medium provided by the application have the following beneficial effects:

[0031] (1) The application is completely based on seismic data for calculation, and the obtained result has strong certainty;

[0032] (2) The application does not refer to well logging interpretation results in the calculation process, and is not affected by human factors. BRIEF DESCRIPTION OF DRAWINGS

[0033] The present application will be described in more detail below based on the embodiments and with reference to the drawings.

[0034] Figure 1 An implementation flowchart of a gas-bearing property detection method provided for an embodiment of the present application;

[0035] Figure 2 A pre-stack angle domain gather after residual moveout correction provided for an embodiment of the present application;

[0036] Figure 3 A pre-stack angle domain gather amplitude frequency spectrum provided for an embodiment of the present application;

[0037] Figure 4 An actual work area data body area_rate well tie profile provided for an embodiment of the present application;

[0038] Figure 5 An actual work area data body max_amp_rate well tie profile provided for an embodiment of the present application;

[0039] Figure 6 An implementation flowchart of a gas-bearing property detection method provided for an embodiment of the present application.

[0040] In the drawings, the same components are designated by the same reference numerals, and the drawings are not drawn according to the actual scale. DETAILED DESCRIPTION

[0041] In order to make the objects, technical solutions and advantages of the present application clearer, the present application will be described in further detail below with reference to the drawings, and the described embodiments should not be regarded as limiting the present application, and all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0042] In the following description, "some embodiments" are described, which describe a subset of all possible embodiments, but it can be understood that "some embodiments" can be the same subset or different subsets of all possible embodiments, and can be combined with each other without conflict.

[0043] If similar descriptions of "first\second\third" appear in the application file, the following description is added, in the following description, the terms "first\second\third" referred to only distinguish similar objects, and do not represent a specific order of the objects, and it can be understood that "first\second\third" can be interchanged in a specific order or sequence as allowed, so that the embodiments of the present application described herein can be implemented in an order other than that illustrated or described herein.

[0044] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used in the description herein is for describing particular embodiments only and is not intended to be limiting of the application.

[0045] Before introducing a gas-bearing property detection method provided by the embodiments of the present application, the problems in the related art are briefly introduced.

[0046] At present, most of the gas-bearing property detection methods are qualitative prediction of post-stack attenuation attributes and quantitative prediction of pre-stack inversion. For qualitative prediction of gas-bearing property by using post-stack attenuation attributes, the information of pre-stack seismic data is lacked. The quantitative prediction of gas-bearing property by using pre-stack inversion is based on the relationship between elastic parameters and gas saturation to calculate the gas saturation, or is based on the neural network algorithm to learn. Since the change of gas saturation in the elastic parameters is very small, the effect of quantitative prediction of gas-bearing property is not ideal. When the quantitative prediction of gas-bearing property is performed, the stability of the gas saturation interpretation result of well logging interpretation is poor. Therefore, if the gas saturation data of well logging interpretation is referred in the prediction of gas-bearing property, the accuracy of the prediction result is also reduced.

[0047] Based on the problems in the related art, the embodiments of the present application provide a gas-bearing property detection method. The method is applied to a gas-bearing property detection device. The gas-bearing property detection device can be an electronic device, such as a computer, a mobile terminal, etc. The function realized by the gas-bearing property detection method provided by the embodiments of the present application can be realized by calling program code by the processor of the electronic device. The program code can be saved in a computer storage medium.

[0048] Example One

[0049] The embodiments of the present application provide a gas-bearing property detection method, Figure 6 The implementation flowchart of the gas-bearing property detection method provided by the embodiments of the present application is shown in FIG. 1. Figure 6 As shown in FIG. 1, the method comprises the following steps.

[0050] S1: Collecting pre-stack angle domain gathers, determining an effective incident angle range and determining a frequency band range of the pre-stack angle domain gathers;

[0051] S2: Performing angle stacking according to the incident angle range to obtain n angle-stacked seismic data. The center angles of the corresponding incident angle stacking ranges are angle1, angle2, angle3…anglen, respectively.

[0052] S3: calculating the amplitude frequency spectrum of each trace of the seismic data after n-angle stacking, obtaining the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum within the effective frequency band range of the seismic data;

[0053] S4: applying the amplitude frequency spectrum area area and the central angle of the incident angle stacking range, obtaining the rate data body area_rate of the change of the seismic data spectrum area in the direction of the incident angle;

[0054] S5: applying the maximum energy max_amp in the amplitude frequency spectrum and the central angle of the incident angle stacking range, obtaining the rate data body max_amp_rate of the change of the maximum energy in the amplitude frequency spectrum in the direction of the incident angle;

[0055] The meaning represented by the data body area_rate: due to the influence of the gas-bearing reservoir, the energy of the seismic data decays, and the change of the energy decay is reflected in the amplitude frequency spectrum, that is, the stacking data spectrum area gradually becomes smaller with the increase of the incident angle, and the total energy of the seismic data decreases with the increase of the incident angle due to the influence of the gas-bearing property, and only the energy decay of the seismic data related to the gas-bearing property is considered without considering the reservoir thickness; The meaning represented by the data body max_amp_rate: due to the influence of the gas-bearing reservoir, the energy of the seismic data decays, and the change of the energy decay is reflected in the amplitude frequency spectrum, that is, the maximum energy in the amplitude frequency spectrum of the stacking data gradually becomes smaller with the increase of the incident angle, and the peak energy of the seismic data decreases with the increase of the incident angle due to the influence of the gas-bearing property;

[0056] S6: calibrating the data body area_rate and the data body max_amp_rate with the known gas production wells in the actual work area to complete the gas-bearing property prediction.

[0057] The application provides a gas-bearing property detection method,

[0058] (1) The application is completely based on seismic data for calculation, and the obtained result has strong certainty;

[0059] (2) The application does not refer to the logging interpretation result in the calculation process, and is not affected by human factors.

[0060] Example Two

[0061] Based on the foregoing embodiment, the application further provides a gas-bearing property detection method, comprising:

[0062] S21: collecting a pre-stack angle domain gather, determining an effective incident angle range and determining an effective frequency band range of the pre-stack angle domain gather;

[0063] In some embodiments, before the determining the effective frequency band range of the pre-stack angle domain gather, the method further comprises:

[0064] performing residual time correction;

[0065] S22: angle stacking is performed according to the incident angle range to obtain n angle stacked seismic data, and the center angles of the corresponding incident angle stacking ranges are angle1, angle2, angle3…anglen respectively;

[0066] S23: the amplitude frequency spectrum of each trace of the n angle stacked seismic data is calculated to obtain the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data;

[0067] S24: the amplitude frequency spectrum area area and the center angle of the incident angle stacking range are applied to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction;

[0068] S25: the maximum energy max_amp in the amplitude frequency spectrum and the center angle of the incident angle stacking range are applied to obtain the variation rate data body max_amp_rate of the maximum energy in the amplitude frequency spectrum in the incident angle direction;

[0069] The meaning represented by the data body area_rate: due to the influence of the reservoir gas content, the energy of the seismic data is attenuated, and the change of the energy attenuation is reflected in the amplitude frequency spectrum as the energy attenuation is stronger with the increase of the incident angle, that is, the spectrum area of the stacked data gradually decreases with the increase of the incident angle, and the total energy of the seismic data decreases with the increase of the incident angle due to the influence of the gas content, and only the energy attenuation of the seismic data related to the gas content is considered without considering the reservoir thickness; The meaning represented by the data body max_amp_rate: due to the influence of the reservoir gas content, the energy of the seismic data is attenuated, and the change of the energy attenuation is reflected in the amplitude frequency spectrum as the energy attenuation is stronger with the increase of the incident angle, that is, the maximum energy in the amplitude frequency spectrum of the stacked data gradually decreases with the increase of the incident angle, and the peak energy of the seismic data decreases with the increase of the incident angle due to the influence of the gas content;

[0070] S26: the data body area_rate and the data body max_amp_rate are calibrated with the known gas production wells in the actual work area to complete the gas content prediction.

[0071] The application provides a gas content detection method,

[0072] (1) The application is completely based on seismic data for calculation, and the obtained result has strong certainty;

[0073] (2) The application does not refer to well logging interpretation results in the calculation process and is not affected by human factors.

[0074] Example Three

[0075] Based on the foregoing embodiments, the embodiments of the application further provide a gas content detection method, comprising:

[0076] S31: Collecting pre-stack angle domain gathers, determining an effective incident angle range, and determining a frequency band range of the pre-stack angle domain gathers;

[0077] In some embodiments, before the determination of the frequency band range of the pre-stack angle domain gathers, the method further comprises:

[0078] Performing remaining time correction;

[0079] S32: Performing angle stacking according to the incident angle range to obtain n angle-stacked seismic data, and the center angles of the corresponding incident angle stacking ranges are angle1, angle2, angle3…anglen, respectively;

[0080] In some embodiments, since the properties of the angle-stacked seismic data need to be extracted later, the signal-to-noise ratio of the seismic data needs to be ensured not to be too low, and therefore the number of stacked data bodies is generally not greater than 5 when angle stacking is performed, and the n is not greater than 5;

[0081] S33: Calculating the amplitude frequency spectrum of each trace of the n angle-stacked seismic data to obtain an amplitude frequency spectrum area area and an energy maximum value max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data;

[0082] S34: Applying the amplitude frequency spectrum area area and the center angle of the incident angle stacking range to obtain a variation rate data body area_rate of the seismic data spectrum area in the incident angle direction;

[0083] S35: Applying the energy maximum value max_amp in the amplitude frequency spectrum and the center angle of the incident angle stacking range to obtain a variation rate data body max_amp_rate of the energy maximum value in the amplitude frequency spectrum in the incident angle direction;

[0084] Data body area_rate represents the meaning: due to the influence of reservoir gas, the energy of seismic data produces attenuation, the change of this energy attenuation is embodied in the amplitude frequency spectrum, that is, the energy attenuation is stronger with the increase of the incident angle, that is, the frequency spectrum area of the stacked data gradually becomes smaller with the increase of the incident angle, the total energy of the seismic data affected by the gas content decreases with the increase of the incident angle, only the energy attenuation of the seismic data related to the gas content is considered, and the reservoir thickness is not considered; Data body max_amp_rate represents the meaning: due to the influence of reservoir gas, the energy of seismic data produces attenuation, the change of this energy attenuation is embodied in the amplitude frequency spectrum, that is, the energy attenuation is stronger with the increase of the incident angle, that is, the maximum energy in the amplitude frequency spectrum of the stacked data gradually becomes smaller with the increase of the incident angle, the peak energy of the seismic data affected by the gas content decreases with the increase of the incident angle.

[0085] S36: calibrate data body area_rate and data body max_amp_rate with the known gas production wells in the actual work area, and complete the gas content prediction.

[0086] The application provides a gas content detection method,

[0087] (1) the application is completely based on seismic data for calculation, and the obtained result has strong certainty;

[0088] (2) the application does not refer to the logging interpretation result in the calculation process, and is not affected by human factors.

[0089] Example Four

[0090] Based on the foregoing embodiment, the application further provides a gas content detection method, which comprises the following steps:

[0091] S41: collecting prestack angle domain gathers, determining an effective incident angle range, and determining a frequency band range of the effective prestack angle domain gathers;

[0092] In some embodiments, before the step of determining the frequency band range of the effective prestack angle domain gathers, the method further comprises the following steps:

[0093] performing residual time correction;

[0094] S42: performing angle stacking according to the incident angle range to obtain n angle-stacked seismic data, and the center angles of the corresponding incident angle stacking ranges are angle1, angle2, angle3…anglen respectively;

[0095] In some embodiments, because the properties of the angle-stacked seismic data need to be extracted later, the signal-to-noise ratio of the seismic data cannot be too low, therefore, when performing angle stacking, the number of stacked data bodies is generally not greater than 5, and n is not greater than 5.

[0096] S43: calculating the amplitude frequency spectrum of each trace of the n-angle stacked seismic data to obtain the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum within the effective frequency band range of the seismic data;

[0097] S44: applying the amplitude frequency spectrum area area and the central angle of the incidence angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction;

[0098] In some embodiments, the specific method of applying the amplitude frequency spectrum area area and the central angle of the incidence angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction comprises:

[0099] Let Y-axis be n amplitude frequency spectrum areas area, and X-axis be the central angles angle1, angle2, angle3…anglen of n incidence angle stacking ranges. Through curve fitting, the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction is obtained;

[0100] S45: applying the maximum energy max_amp in the amplitude frequency spectrum and the central angle of the incidence angle stacking range to obtain the variation rate data body max_amp_rate of the maximum energy in the amplitude frequency spectrum in the incidence angle direction;

[0101] The meaning represented by the data body area_rate: due to the influence of reservoir gas content, the energy of seismic data decays. The change of this energy decay is reflected in the amplitude frequency spectrum as the energy decay becomes stronger with the increase of the incidence angle, that is, the stacking data spectrum area gradually becomes smaller with the increase of the incidence angle. The total energy of the seismic data affected by the gas content decreases with the increase of the incidence angle. Here, only the energy decay of the seismic data related to the gas content is considered, and the reservoir thickness is not considered. The meaning represented by the data body max_amp_rate: due to the influence of reservoir gas content, the energy of seismic data decays. The change of this energy decay is reflected in the amplitude frequency spectrum as the energy decay becomes stronger with the increase of the incidence angle, that is, the maximum energy in the amplitude frequency spectrum of the stacking data gradually becomes smaller with the increase of the incidence angle. The peak energy of the seismic data affected by the gas content decreases with the increase of the incidence angle;

[0102] S46: calibrating the data body area_rate and the data body max_amp_rate with the known gas production wells in the actual work area to complete the gas content prediction.

[0103] The application provides a gas content detection method,

[0104] (1) the present application is completely based on seismic data to calculate, the obtained result has strong certainty;

[0105] (2) the present application does not refer to well logging interpretation result in the calculation process, is not influenced by artificial factors.

[0106] Example Five

[0107] Based on the foregoing embodiments, the embodiments of the application further provide a gas content detection method, comprising:

[0108] S51: collect prestack angle domain gathers, determine effective incidence angle range and determine effective frequency band range of prestack angle domain gathers;

[0109] In some embodiments, before the determination of the effective frequency band range of the prestack angle domain gathers, the method further comprises:

[0110] Performing remaining time correction;

[0111] S52: angle stacking according to the incidence angle range, obtaining n angle stacked seismic data, and the corresponding center angles of the incidence angle stacking range are angle1, angle2, angle3…anglen respectively;

[0112] In some embodiments, since the properties of the angle stacked seismic data need to be extracted afterwards, the signal-to-noise ratio of the seismic data cannot be too low, therefore, when angle stacking is performed, the number of stacked data bodies is generally not more than 5, and the n is not more than 5;

[0113] S53: calculating the amplitude frequency spectrum of each trace of the n angle stacked seismic data, obtaining the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data;

[0114] S54: applying the amplitude frequency spectrum area area and the center angle of the incidence angle stacking range, obtaining the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction;

[0115] In some embodiments, the specific method for applying the amplitude frequency spectrum area area and the center angle of the incidence angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction comprises:

[0116] Let Y axis be n amplitude frequency spectrum areas area, X axis be n center angles angle1, angle2, angle3…anglen of incidence angle stacking ranges, and through curve fitting, obtain the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction;

[0117] S55: applying the energy maximum value max_amp in the amplitude frequency spectrum and the central angle of the incident angle stacking range to obtain the variation rate data body max_amp_rate of the energy maximum value in the amplitude frequency spectrum in the direction of the incident angle;

[0118] In some embodiments, the specific method of applying the energy maximum value max_amp in the amplitude frequency spectrum and the central angle of the incident angle stacking range to obtain the variation rate data body max_amp_rate of the energy maximum value in the amplitude frequency spectrum in the direction of the incident angle includes:

[0119] Let Y axis be the energy maximum value max_amp in the n amplitude frequency spectrums, and X axis be the central angles angle1, angle2, angle3…anglen of the n incident angle stacking ranges. Through curve fitting, the variation rate data body max_amp_rate of the energy maximum value in the amplitude frequency spectrum in the direction of the incident angle is obtained.

[0120] The meaning represented by the data body area_rate: due to the influence of reservoir gas content, the energy of seismic data decays. The change of this energy decay is reflected in the amplitude frequency spectrum as the energy decay becomes stronger with the increase of the incident angle, that is, the stacking data spectrum area gradually becomes smaller with the increase of the incident angle. The total energy of the seismic data affected by the gas content decreases with the increase of the incident angle. Here, only the energy decay of the seismic data related to the gas content is considered, and the reservoir thickness is not considered. The meaning represented by the data body max_amp_rate: due to the influence of reservoir gas content, the energy of seismic data decays. The change of this energy decay is reflected in the amplitude frequency spectrum as the energy maximum value in the amplitude frequency spectrum gradually becomes smaller with the increase of the incident angle, that is, the peak energy of the seismic data affected by the gas content decreases with the increase of the incident angle.

[0121] S56: calibrating the data body area_rate and the data body max_amp_rate with the known gas production wells in the actual work area to complete the gas content prediction.

[0122] The application provides a gas content detection method,

[0123] (1) The application is completely based on seismic data for calculation, and the obtained result has strong certainty;

[0124] (2) The application does not refer to the logging interpretation result in the calculation process, and is not affected by human factors.

[0125] Example Six

[0126] Based on the foregoing embodiments, the application further provides a gas content detection method, which comprises:

[0127] S61: collect pre-stack angle domain gathers, determine the effective incident angle range and determine the effective frequency band range of the pre-stack angle domain gathers;

[0128] In some embodiments, the effective incident angle range is 0-36;

[0129] In some embodiments, before the determining the effective frequency band range of the pre-stack angle domain gathers, the method further comprises:

[0130] performing residual time correction;

[0131] S62: perform angle stacking according to the incident angle range, obtain n angle stacked seismic data, and the corresponding center angles of the incident angle stacking range are angle1, angle2, angle3…anglen respectively;

[0132] In some embodiments, since the properties of the angle stacked seismic data need to be extracted later, the signal-to-noise ratio of the seismic data cannot be too low, and therefore the number of stacked data bodies is generally not greater than 5 when angle stacking is performed, and the n is not greater than 5;

[0133] S63: calculate the amplitude frequency spectrum of each trace of the n angle stacked seismic data, and obtain the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data;

[0134] S64: apply the amplitude frequency spectrum area area and the center angle of the incident angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction;

[0135] In some embodiments, the specific method of applying the amplitude frequency spectrum area area and the center angle of the incident angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction comprises:

[0136] let the Y axis be the n amplitude frequency spectrum areas area, the X axis be the center angles angle1, angle2, angle3…anglen of the n incident angle stacking ranges, and the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction is obtained through curve fitting;

[0137] S65: apply the maximum energy max_amp in the amplitude frequency spectrum and the center angle of the incident angle stacking range to obtain the variation rate data body max_amp_rate of the maximum energy in the amplitude frequency spectrum in the incident angle direction;

[0138] In some embodiments, the specific method of superimposing the center angle of the incident angle range on the energy maximum max_amp in the amplitude frequency spectrum to obtain the variation rate data body max_amp_rate of the energy maximum in the amplitude frequency spectrum in the direction of the incident angle includes:

[0139] Let Y axis be the energy maximum max_amp in the n amplitude frequency spectrums, and X axis be the center angles angle1, angle2, angle3…anglen of the n incident angle superimposition ranges, and the variation rate data body max_amp_rate of the energy maximum in the amplitude frequency spectrum in the direction of the incident angle is obtained through curve fitting.

[0140] The meaning represented by the data body area_rate is that due to the influence of the reservoir gas content, the energy of the seismic data is attenuated, and the change of the energy attenuation is embodied in the amplitude frequency spectrum as the energy attenuation is stronger with the increase of the incident angle, that is, the superimposed data spectrum area gradually becomes smaller with the increase of the incident angle, and the total energy of the seismic data affected by the gas content decreases with the increase of the incident angle, and only the energy attenuation of the seismic data related to the gas content is considered without considering the reservoir thickness; The meaning represented by the data body max_amp_rate is that due to the influence of the reservoir gas content, the energy of the seismic data is attenuated, and the change of the energy attenuation is embodied in the amplitude frequency spectrum as the energy maximum in the amplitude frequency spectrum gradually becomes smaller with the increase of the incident angle, that is, the peak energy of the seismic data affected by the gas content decreases with the increase of the incident angle.

[0141] S66: Calibrate the data body area_rate and the data body max_amp_rate with the known gas production wells in the actual work area to complete the gas content prediction.

[0142] The application provides a gas content detection method,

[0143] (1) The application is completely based on the calculation of the seismic data, and the obtained result has strong certainty;

[0144] (2) The application does not refer to the logging interpretation result in the calculation process, and is not affected by human factors.

[0145] Example Seven

[0146] Based on the foregoing embodiments, the application further provides a gas content detection method, which includes:

[0147] S71: Collect the prestack angle domain gather, determine the effective incident angle range, and determine the effective frequency band range of the prestack angle domain gather;

[0148] In some embodiments, the effective incident angle range is 0-36.

[0149] In some embodiments, the effective frequency band range of the pre-stack angle domain gather is 7-75 Hz;

[0150] In some embodiments, before the determining the effective frequency band range of the pre-stack angle domain gather, the method further comprises:

[0151] performing remaining time correction;

[0152] S72: angle stacking according to the incident angle range to obtain n angle stacked seismic data, and the corresponding center angles of the incident angle stacking range are angle1, angle2, angle3…anglen respectively;

[0153] In some embodiments, since the angle stacked seismic data needs to be extracted for attributes later, the signal-to-noise ratio of the seismic data needs to be ensured not to be too low, and therefore the number of stacked data bodies is generally not greater than 5 when angle stacking is performed, and the n is not greater than 5;

[0154] S73: calculating the amplitude frequency spectrum of each trace of the n angle stacked seismic data to obtain the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data;

[0155] S74: applying the amplitude frequency spectrum area area and the center angle of the incident angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction;

[0156] In some embodiments, the specific method of applying the amplitude frequency spectrum area area and the center angle of the incident angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction comprises:

[0157] letting the Y axis be the n amplitude frequency spectrum areas area, the X axis be the center angles angle1, angle2, angle3…anglen of the n incident angle stacking ranges, and obtaining the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction through curve fitting;

[0158] S75: applying the maximum energy max_amp in the amplitude frequency spectrum and the center angle of the incident angle stacking range to obtain the variation rate data body max_amp_rate of the maximum energy in the amplitude frequency spectrum in the incident angle direction;

[0159] In some embodiments, the specific method of superimposing the center angle of the n incident angle superposition range on the amplitude frequency spectrum energy maximum max_amp to obtain the amplitude frequency spectrum energy maximum rate of change data body max_amp_rate in the direction of the incident angle includes:

[0160] Let Y axis be the n amplitude frequency spectrum energy maximum max_amp, and X axis be the n incident angle superposition range center angle angle1, angle2, angle3…anglen, and the amplitude frequency spectrum energy maximum rate of change data body max_amp_rate in the direction of the incident angle is obtained through curve fitting.

[0161] The meaning represented by the data body area_rate is that due to the influence of the reservoir gas content, the energy of the seismic data is attenuated, and the change of this energy attenuation is reflected in the amplitude frequency spectrum as the energy attenuation being stronger with the increase of the incident angle, that is, the superposition data spectrum area gradually becomes smaller with the increase of the incident angle, and the total energy of the seismic data affected by the gas content decreases with the increase of the incident angle, and here, only the seismic data energy attenuation related to the gas content is considered without considering the reservoir thickness; the meaning represented by the data body max_amp_rate is that due to the influence of the reservoir gas content, the energy of the seismic data is attenuated, and the change of this energy attenuation is reflected in the amplitude frequency spectrum as the amplitude frequency spectrum energy maximum gradually becoming smaller with the increase of the incident angle, that is, the peak value of the seismic data affected by the gas content decreases with the increase of the incident angle.

[0162] S76: Calibrate the data body area_rate and the data body max_amp_rate with the known gas production wells in the actual work area to complete the gas content prediction.

[0163] The application provides a gas content detection method,

[0164] (1) The application is completely based on seismic data for calculation, and the obtained result has strong certainty;

[0165] (2) The application does not refer to the logging interpretation result in the calculation process, and is not affected by human factors.

[0166] Example Eight

[0167] Based on the method of embodiment seven, as shown in Figure 1 and Figure 6 The application according to the real data gives the embodiments as shown in

[0168] S81: Collect the prestack angle domain gather, determine the effective incident angle range, and determine the effective frequency band range of the prestack angle domain gather;

[0169] In some embodiments, as shown in FIG. 6, the effective incidence angle range is 0-36; Figure 2

[0170] In some embodiments, as shown in FIG. 6, the effective incidence angle range is 0-36; Figure 3

[0171] In some embodiments, before determining the effective frequency band range of the pre-stack angle domain gather, the method further comprises:

[0172] performing residual time correction;

[0173] S82: angle stacking according to the incidence angle range to obtain n angle stacked seismic data, and the center angles of the corresponding incidence angle stacking ranges are angle1, angle2, angle3…anglen, respectively;

[0174] In some embodiments, because the angle stacked seismic data needs to be extracted for attributes later, the signal-to-noise ratio of the seismic data needs to be ensured not to be too low, and therefore the number of stacked data bodies is generally not greater than 5 when angle stacking is performed, and n is not greater than 5;

[0175] In some embodiments, angle stacking is performed according to the incidence angle range to obtain 3 angle stacked seismic data, and the incidence angle ranges of each angle stacked data are 0-12, 12-24, and 24-36, respectively. The 3 angle stacked seismic data can be expressed as seis_6, seis_18, and seis_30 according to the center angle.

[0176] S83: calculating the amplitude frequency spectrum of each trace of the n angle stacked seismic data to obtain the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data;

[0177] In some embodiments, the amplitude frequency spectrum of each trace of the 3 angle stacked seismic data is calculated to obtain the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum in the effective frequency band of the seismic data. The data bodies area_6, area_18, and area_30 represent the amplitude frequency spectrum area information of the seismic data in different incidence angle ranges, and the data bodies max_amp_6, max_amp_18, and max_amp_30 represent the maximum energy information of the amplitude frequency spectrum of the seismic data in different incidence angle ranges.

[0178] S84: applying the amplitude frequency spectrum area area and the center angle of the incidence angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction; ​​

[0179] In some embodiments, the specific method of applying the amplitude frequency spectrum area area and the center angle of the incident angle stack range to obtain the rate of change data volume area rate of the amplitude frequency spectrum area in the direction of the incident angle includes:

[0180] Let the Y axis be the n amplitude frequency spectrum areas area, and the X axis be the center angles angle1, angle2, angle3…anglen of the n incident angle stack ranges. Through curve fitting, the rate of change data volume area rate of the amplitude frequency spectrum area in the direction of the incident angle is obtained.

[0181] In some embodiments, let the Y axis be the amplitude frequency spectrum areas area_fre at the center angles 6, 18, and 30, and the X axis be 6, 18, and 30. Through linear regression, the corresponding slope values are calculated to form the change slope volume area_rate of the amplitude frequency spectrum area with the increase of the incident angle.

[0182] S85: Apply the maximum energy value max_amp in the amplitude frequency spectrum and the center angle of the incident angle stack range to obtain the rate of change data volume max_amp_rate of the maximum energy value in the amplitude frequency spectrum in the direction of the incident angle.

[0183] In some embodiments, the specific method of applying the maximum energy value max_amp in the amplitude frequency spectrum and the center angle of the incident angle stack range to obtain the rate of change data volume max_amp_rate of the maximum energy value in the amplitude frequency spectrum in the direction of the incident angle includes:

[0184] Let the Y axis be the n maximum energy values max_amp in the amplitude frequency spectrum, and the X axis be the center angles angle1, angle2, angle3…anglen of the n incident angle stack ranges. Through curve fitting, the rate of change data volume max_amp_rate of the maximum energy value in the amplitude frequency spectrum in the direction of the incident angle is obtained.

[0185] In some embodiments, according to step 5, let the Y axis be the maximum energy values max_fre in the amplitude frequency spectrum at the center angles 6, 18, and 30, and the X axis be 6, 18, and 30. Through linear regression, the corresponding slope values are calculated to form the change slope volume max_amp_rate of the maximum energy value in the amplitude frequency spectrum with the increase of the incident angle.

[0186] The data body area_rate represents the meaning: due to the influence of reservoir gas, the energy of seismic data produces attenuation, and the change of the energy attenuation is embodied in the amplitude frequency spectrum, that is, the energy attenuation is stronger with the increase of the incidence angle, that is, the frequency spectrum area of the stacked data gradually becomes smaller with the increase of the incidence angle, and the total energy of the seismic data affected by the gas saturation decreases with the increase of the incidence angle; the data body max_amp_rate represents the meaning: due to the influence of reservoir gas, the energy of seismic data produces attenuation, and the change of the energy attenuation is embodied in the amplitude frequency spectrum, that is, the energy maximum value in the amplitude frequency spectrum of the stacked data gradually becomes smaller with the increase of the incidence angle, and the peak energy of the seismic data affected by the gas saturation decreases with the increase of the incidence angle.

[0187] S86: calibrate the data body area_rate and the data body max_amp_rate with the known gas production wells in the actual work area, and complete the gas saturation prediction.

[0188] In some embodiments, the data body area_rate and the data body max_amp_rate are calibrated with the known gas production wells in the actual work area. There are three known wells in the work area, well A is a gas production well (gas saturation 70%), well B is a water production well (gas saturation 40%), well C is a super high gas production well, and well D is a poor gas production well. Wells A and B are used as calibration wells, and wells C and D are used as verification wells.

[0189] Figure 4 The data body area_rate is shown in the figure, and the values of the four wells A, B, C and D are-0.14, -0.04, -0.2 and-0.1 respectively. Well A is a gas production well (gas saturation 70%), and well B is a water production well (gas saturation 40%), which are used as calibration wells, so it can be determined that the gas saturation of well C is obviously higher than 70%, and the gas saturation of well D is between 40-70%, which is completely consistent with the actual yield of the well.

[0190] Figure 5 The data body max_amp_rate is shown in the figure, and the values of the four wells A, B, C and D are-0.1, -0.04, -0.14 and-0.08 respectively. Well A is a gas production well (gas saturation 70%), and well B is a water production well (gas saturation 40%), which are used as calibration wells, so it can be determined that the gas saturation of well C is obviously higher than 70%, and the gas saturation of well D is between 40-70%, which is completely consistent with the actual yield of the well. This method can effectively predict the gas saturation.

[0191] The application provides a gas saturation detection method,

[0192] (1) the present application is completely based on seismic data to calculate, the result is strong in certainty;

[0193] (2) the present application does not refer to the logging interpretation result in the calculation process, is not influenced by artificial factors.

[0194] Example Nine

[0195] Based on the foregoing embodiment, the embodiment of the application provides a gas content detection device, comprising:

[0196] The prestack angle domain gather collection module, the central angle calculation module, the amplitude frequency spectrum data calculation module, the data body area_rate calculation module, the data body max_amp_rate calculation module and the gas content prediction module are used to calculate the central angle of the corresponding incident angle stacking range angle1, angle2, angle3…anglen.

[0197] The prestack angle domain gather collection module: collecting prestack angle domain gathers, determining the effective incident angle range and determining the effective frequency band range of the prestack angle domain gathers;

[0198] The central angle calculation module: angle stacking is performed according to the incident angle range, n angle-stacked seismic data are obtained, and the central angles of the corresponding incident angle stacking ranges are angle1, angle2, angle3…anglen.

[0199] The amplitude frequency spectrum data calculation module: the amplitude frequency spectrum of each trace of the n angle-stacked seismic data is calculated, the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data are obtained; wherein, before determining the effective frequency band range of the prestack angle domain gathers, it further includes: performing residual time correction, and the n is not greater than 5.

[0200] The data body area_rate calculation module: the amplitude frequency spectrum area area and the central angle of the incident angle stacking range are applied to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction.

[0201] The specific method for applying the amplitude frequency spectrum area area and the central angle of the incident angle stacking range to obtain the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction includes:

[0202] Let Y axis be n amplitude frequency spectrum areas area, X axis be n central angles angle1, angle2, angle3…anglen of incident angle stacking ranges, the variation rate data body area_rate of the seismic data spectrum area in the incident angle direction is obtained through curve fitting.

[0203] The data body max_amp_rate is calculated by using the energy maximum value max_amp in the amplitude frequency spectrum and the central angle of the incident angle stacking range.

[0204] The specific method for obtaining the data body max_amp_rate of the variation rate of the energy maximum value in the amplitude frequency spectrum in the direction of the incident angle by using the energy maximum value max_amp in the amplitude frequency spectrum and the central angle of the incident angle stacking range comprises the following steps:

[0205] The Y axis is the energy maximum value max_amp in the n amplitude frequency spectrums, and the X axis is the central angle angle1, angle2, angle3…anglen of the n incident angle stacking ranges. The data body max_amp_rate of the variation rate of the energy maximum value in the amplitude frequency spectrum in the direction of the incident angle is obtained by curve fitting.

[0206] The gas-bearing property prediction module calibrates the data body area_rate and the data body max_amp_rate with the known gas production wells in the actual work area, and completes the gas-bearing property prediction.

[0207] The effective incident angle range is 0-36.

[0208] The effective frequency band range of the pre-stack angle domain gather is 7-75 Hz.

[0209] It should be noted that, in the embodiments of the present application, if the gas-bearing property detection method is realized in the form of a software function module and is sold or used as an independent product, it can also be stored in a computer-readable storage medium. Based on this understanding, the technical solutions of the embodiments of the present application can be embodied in the form of a software product, and the computer software product is stored in a storage medium, includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the method described in the embodiments of the present application. The foregoing storage medium includes a U disk, a mobile hard disk, a read-only memory (ROM, Read Only Memory), a magnetic disk or an optical disk, and various storage medium that can store program codes. Thus, the embodiments of the present application are not limited to any specific hardware and software combination.

[0210] Correspondingly, the embodiments of the present application provide a storage medium having a computer program stored thereon, characterized by comprising computer program instructions, wherein when the computer program instructions are executed by a processor, the steps of the gas-bearing property detection method provided in the above embodiments are implemented.

[0211] Embodiment ten

[0212] The embodiment of the present application provides a gas-containing detection device memory and a processor, the memory stores a computer program, the computer program is executed by the processor, the processor is configured to execute the program of the multi-scale electromagnetic field component denoising method stored in the memory, so that the steps in the multi-scale electromagnetic field component denoising method provided in the above embodiment are realized.

[0213] The description of the above display device and storage medium embodiments is similar to the description of the above method embodiments, and has similar beneficial effects to the method embodiments. For technical details not disclosed in the computer device and storage medium embodiments of the present application, please refer to the description of the method embodiments of the present application.

[0214] It should be pointed out here that the description of the above storage medium and device embodiments is similar to the description of the above method embodiments, and has similar beneficial effects to the method embodiments. For technical details not disclosed in the storage medium and device embodiments of the present application, please refer to the description of the method embodiments of the present application.

[0215] It should be understood that the "one embodiment" or "an embodiment" mentioned throughout the specification means that the specific features, structures or characteristics related to the embodiment are included in at least one embodiment of the present application. Therefore, "in one embodiment" or "in an embodiment" appearing throughout the specification does not necessarily refer to the same embodiment. In addition, these specific features, structures or characteristics can be combined in one or more embodiments in any suitable manner. It should be understood that in various embodiments of the present application, the size of the sequence number of each process does not mean the order of execution, and the execution order of each process should be determined according to its function and inherent logic, and should not constitute any limitation on the implementation process of the embodiments of the present application. The sequence number of the above embodiments of the present application is only for description, not representing the advantages and disadvantages of the embodiments.

[0216] It should be noted that in this paper, the term "include", "contain" or any other variant thereof is intended to cover non-exclusive inclusion, so that the process, method, article or device including a series of elements not only includes those elements, but also includes other elements not explicitly listed, or includes elements inherent to such process, method, article or device. Without more limitations, the element defined by the sentence "including a…" does not exclude the presence of other identical elements in the process, method, article or device including the element.

[0217] In several embodiments provided in the present application, it should be understood that the disclosed devices and methods can be implemented in other manners. The described device embodiments are merely schematic. For example, the division of the units is only a logical function division. There can be another division manner for the actual implementation, for example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the displayed or discussed coupling, or direct coupling or communication connection between the components can be indirect coupling or communication connection through some interfaces, devices, or units, and can be electrical, mechanical, or in other forms.

[0218] The units described as separate components can or can not be physically separate, and the components shown as units can or can not be physical units; they can be located in one place, or distributed on multiple network units; and some or all of the units can be selected according to actual needs to achieve the purpose of the embodiment.

[0219] In addition, each functional unit in each embodiment of the present application can be integrated into one processing unit, or each unit can be a separate unit, or two or more units can be integrated into one unit; the integrated unit can be realized in the form of hardware, or in the form of hardware plus software functional units.

[0220] Those of ordinary skill in the art can understand that all or part of the steps of the above-described method embodiments can be completed by a program instructing related hardware, and the foregoing program can be stored in a computer readable storage medium, and when the program is executed, the steps of the method embodiments are executed.

[0221] Alternatively, when the integrated units of the present application are implemented in the form of software functional modules and sold or used as independent products, they can also be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the embodiments of the present application can be embodied in the form of a software product, and the computer software product is stored in a storage medium, and includes several instructions for causing a computer to execute all or part of the methods described in the embodiments of the present application. The storage medium includes mobile storage devices, ROM, magnetic disks, or optical disks, and various other media that can store program codes.

[0222] The above merely provides the implementation of the present application, but the protection scope of the present application is not limited thereto, any person skilled in the art can easily think of the change or replacement within the technical range disclosed by the present application, which should be covered in the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A method for detecting gas content, characterized by, The method comprises the following steps: S1: collecting prestack angle domain gathers, determining an effective incidence angle range, and determining a frequency band range in which the prestack angle domain gathers are effective; S2: performing angle stacking according to the incidence angle range to obtain n angle-stacked seismic data, and the central angles of the corresponding incidence angle stacking ranges are angle1, angle2, angle3…anglen, respectively; S3: calculating the amplitude frequency spectrum of each trace of the n angle-stacked seismic data to obtain the amplitude frequency spectrum area area and the maximum energy value max_amp in the amplitude frequency spectrum in the effective frequency band range of the seismic data; S4: applying the amplitude frequency spectrum area area and the central angles of the incidence angle stacking ranges to obtain a variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction; S5: applying the maximum energy value max_amp in the amplitude frequency spectrum and the central angles of the incidence angle stacking ranges to obtain a variation rate data body max_amp_rate of the maximum energy value in the amplitude frequency spectrum in the incidence angle direction; S6: calibrating the data body area_rate and the data body max_amp_rate with known gas production wells in an actual work area to complete gas-bearing property prediction.

2. The method of claim 1, wherein, Before the step of determining the frequency band range in which the prestack angle domain gathers are effective, the method further comprises the following step: performing residual time correction.

3. The method of claim 1, wherein, The n is not greater than 5.

4. The method of claim 1, wherein, The specific method of applying the amplitude frequency spectrum area area and the central angles of the incidence angle stacking ranges to obtain the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction comprises the following steps: letting the Y axis be the n amplitude frequency spectrum areas area, the X axis be the central angles angle1, angle2, angle3…anglen of the n incidence angle stacking ranges, and obtaining the variation rate data body area_rate of the seismic data spectrum area in the incidence angle direction through curve fitting.

5. The method of claim 1, wherein, The specific method of applying the maximum energy value max_amp in the amplitude frequency spectrum and the central angles of the incidence angle stacking ranges to obtain the variation rate data body max_amp_rate of the maximum energy value in the amplitude frequency spectrum in the incidence angle direction comprises the following steps: letting the Y axis be the n maximum energy values max_amp in the amplitude frequency spectrum, the X axis be the central angles angle1, angle2, angle3…anglen of the n incidence angle stacking ranges, and obtaining the variation rate data body max_amp_rate of the maximum energy value in the amplitude frequency spectrum in the incidence angle direction through curve fitting.

6. The method of claim 1, wherein, The effective incidence angle range is 0-36.

7. The method of claim 1, wherein, The frequency band range in which the prestack angle domain gathers are effective is 7-75 Hz.

8. A gas content detecting device characterized by comprising: The method comprises the following steps: a prestack angle domain gather collection module, a central angle calculation module, an amplitude frequency spectrum data calculation module, a data body area_rate calculation module, a data body max_amp_rate calculation module, and a gas-bearing property prediction module; the prestack angle domain gather collection module collects prestack angle domain gathers, determines an effective incidence angle range, and determines a frequency band range in which the prestack angle domain gathers are effective; The center angle calculation module: according to the incident angle range, the angle superposition is carried out, n angle superposition seismic data are obtained, and the corresponding incident angle superposition range center angles are angle1, angle2, angle3…anglen respectively; The amplitude frequency spectrum data calculation module: the amplitude frequency spectrum of each trace of n angle superposition seismic data is calculated, the amplitude frequency spectrum area area and the maximum energy max_amp in the amplitude frequency spectrum within the effective frequency band range of the seismic data are obtained; The data body area_rate calculation module: the amplitude frequency spectrum area area and the center angle of the incident angle superposition range are applied to obtain the data body area_rate of the change rate of the seismic data spectrum area in the incident angle direction; The data body max_amp_rate calculation module: the maximum energy max_amp in the amplitude frequency spectrum and the center angle of the incident angle superposition range are applied to obtain the data body max_amp_rate of the change rate of the maximum energy in the amplitude frequency spectrum in the incident angle direction; The gas-bearing property prediction module: the data body area_rate and the data body max_amp_rate are calibrated with the known gas production wells in the actual work area to complete the gas-bearing property prediction.

9. A gas content detecting apparatus characterized by comprising: The memory and the processor are included, the computer program is stored on the memory, and the computer program is executed by the processor to execute the gas-bearing property detection method in any one of claims 1 to 7.

10. A storage medium, characterized by The computer program stored in the storage medium can be executed by one or more processors and can be used to implement the gas-bearing property detection method in any one of claims 1 to 7.

Images