CN113900141B - Oil-gas distribution prediction method and device - Google Patents

Oil-gas distribution prediction method and device Download PDF

Info

Publication number
CN113900141B
CN113900141B CN202010639251.XA CN202010639251A CN113900141B CN 113900141 B CN113900141 B CN 113900141B CN 202010639251 A CN202010639251 A CN 202010639251A CN 113900141 B CN113900141 B CN 113900141B
Authority
CN
China
Prior art keywords
stack
post
dimensional seismic
seismic data
data body
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010639251.XA
Other languages
Chinese (zh)
Other versions
CN113900141A (en
Inventor
于亮亮
刘宝鸿
刘志江
张高
庚琪
刘飞
毕二刚
苑晓娇
刘东旭
郭嘉丽
仇劲涛
倪唯一
陈星州
王一鸣
王乔
陈俊杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010639251.XA priority Critical patent/CN113900141B/en
Publication of CN113900141A publication Critical patent/CN113900141A/en
Application granted granted Critical
Publication of CN113900141B publication Critical patent/CN113900141B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/307Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/63Seismic attributes, e.g. amplitude, polarity, instant phase

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides an oil gas distribution prediction method and device, wherein the method comprises the following steps: according to the post-stack three-dimensional seismic data, using a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window, and establishing a post-stack three-dimensional seismic data volume; gray scale calculation is carried out on the three-dimensional seismic data body after the stacking to obtain a gray scale data body; gain amplification is carried out on the three-dimensional seismic data body after the stacking according to the gray data body, so that a gain data body is obtained; and obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body. The method can improve the accuracy of oil gas distribution prediction, and further provides effective reference for oil gas exploration and development well position deployment and decision optimization.

Description

Oil-gas distribution prediction method and device
Technical Field
The invention belongs to the field of oil and gas geophysical prospecting engineering, and particularly relates to an oil and gas distribution prediction method and device.
Background
The seismic attribute technology originates in the 70 th century of 20 and develops over 40 years, and has developed from a single bright point identification technology into a strong attribute group with multidisciplinary background and multiple technical means, and the whole flow of the seismic interpretation work is penetrated. The seismic attribute technology has better computing capacity and certain geological fitness in the aspects of fault identification, reflection horizon tracking, reservoir inversion, physical parameter prediction, oil gas detection and the like.
The oil gas detection technology is divided into two types of pre-stack and post-stack according to different adopted data. The prestack oil gas detection technology is related attribute calculation and analysis work based on AVO measurement results, but has large data volume and calculation amount, has higher requirements on various aspects such as hardware facilities, system environments and the like, and is limited by the technical characteristics of weak near offset reflection energy, easy occurrence of stretching distortion of far offset of prestack channels, so that in daily work, post-stack seismic data are often adopted to carry out oil gas detection work. The post-stack oil gas detection technology is generally established on a bidirectional medium theory (generalized), mainly comprises an effective medium theory, a wave equation, a self-adaptive theory and the like, and is mainly developed on the basis of a low-frequency resonance and high-frequency attenuation principle, but is influenced by a detection time window, the actual thickness of an oil and gas layer, oil and gas saturation, continuous low-frequency strong reflection information in seismic data and the like, and the predicted oil and gas distribution result is difficult to match with the actual oil and gas distribution characteristics.
Disclosure of Invention
The embodiment of the invention provides an oil gas distribution prediction method, which is used for improving the accuracy of oil gas distribution prediction, and comprises the following steps:
According to the post-stack three-dimensional seismic data, using a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window, and establishing a post-stack three-dimensional seismic data volume;
Gray scale calculation is carried out on the three-dimensional seismic data body after the stacking to obtain a gray scale data body;
gain amplification is carried out on the three-dimensional seismic data body after the stacking according to the gray data body, so that a gain data body is obtained;
and obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body.
The embodiment of the invention also provides an oil gas distribution prediction device, which is used for improving the accuracy of oil gas distribution prediction, and comprises the following steps:
The post-stack three-dimensional seismic data volume module is used for establishing a post-stack three-dimensional seismic data volume by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window according to the post-stack three-dimensional seismic data;
the gray data body module is used for carrying out gray calculation on the three-dimensional seismic data body after the stacking to obtain a gray data body;
the gain data body module is used for carrying out gain amplification on the three-dimensional seismic data body after the stacking according to the gray data body to obtain a gain data body;
And the oil-gas distribution prediction module is used for obtaining an oil-gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the method when executing the computer program.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the above method.
In the embodiment of the invention, according to the post-stack three-dimensional seismic data, a post-stack three-dimensional seismic data volume is established by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window; gray scale calculation is carried out on the three-dimensional seismic data body after the stacking to obtain a gray scale data body; gain amplification is carried out on the three-dimensional seismic data body after the stacking according to the gray data body, so that a gain data body is obtained; and obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body. Compared with the post-stack oil gas distribution prediction method based on spectrum analysis in the prior art, the method provided by the invention is based on the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body, so that the oil gas distribution prediction result is obtained, the result has higher coincidence degree with the actual oil gas distribution characteristics, the adverse effect of continuous low-frequency strong reflection information on oil gas prediction is reduced, the accuracy of oil gas distribution prediction is improved, and effective reference can be provided for well position deployment and decision optimization in oil gas exploration and development.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions of the prior art, the drawings which are used in the description of the embodiments or the prior art will be briefly described, it being obvious that the drawings in the description below are only some embodiments of the invention, and that other drawings can be obtained according to these drawings without inventive faculty for a person skilled in the art. In the drawings:
FIG. 1 is a flow chart of a method for predicting oil and gas distribution according to an embodiment of the present invention;
FIG. 2 is a schematic cross-sectional view of a seismic line of a gray scale data volume B according to an embodiment of the present invention;
FIG. 3 is a schematic cross-sectional view of a seismic line of a gain data volume P according to an embodiment of the invention;
FIG. 4 is a schematic diagram showing the planar characteristics of the predicted result E of the oil and gas distribution according to an embodiment of the present invention;
Fig. 5 is a schematic structural diagram of an oil-gas distribution prediction apparatus according to an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Seismic co-axial events corresponding to hydrocarbon formations often exhibit amplitude energy anomalies, which have been demonstrated and accepted in seismic mechanics studies. However, the amplitude variation characteristics caused by oil gas do not have the same standard or visual characteristics, and the oil gas enrichment beneficial zone cannot be defined simply by the amplitude variation on the seismic profile.
Therefore, the embodiment of the invention provides a method for predicting the oil gas distribution, which improves the accuracy of the oil gas distribution prediction by fully excavating the oil gas information carried by the seismic data. FIG. 1 is a flow chart of a method for predicting oil and gas distribution according to an embodiment of the present invention. As shown in fig. 1, the method for predicting the oil-gas distribution in the embodiment of the present invention may include:
step 101, building a post-stack three-dimensional seismic data volume according to the post-stack three-dimensional seismic data by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window;
102, carrying out gray calculation on the three-dimensional seismic data body after the stacking to obtain a gray data body;
Step 103, gain amplification is carried out on the three-dimensional seismic data body after the stacking according to the gray data body, so as to obtain a gain data body;
104, obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body;
As can be seen from the flow shown in fig. 1, in the embodiment of the present invention, according to the post-stack three-dimensional seismic data, a post-stack three-dimensional seismic data volume is established by using the seismic interpretation layer of the post-stack three-dimensional seismic data as a time window; gray scale calculation is carried out on the three-dimensional seismic data body after the stacking to obtain a gray scale data body; gain amplification is carried out on the three-dimensional seismic data body after the stacking according to the gray data body, so that a gain data body is obtained; and obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body. Compared with the post-stack oil gas distribution prediction method based on spectrum analysis in the prior art, the method provided by the invention is based on the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body, so that the oil gas distribution prediction result is obtained, the result has higher coincidence degree with the actual oil gas distribution characteristics, the adverse effect of continuous low-frequency strong reflection information on oil gas prediction is reduced, the accuracy of oil gas distribution prediction is improved, and effective reference can be provided for well position deployment and decision optimization in oil gas exploration and development.
When the method is implemented, firstly, a post-stack three-dimensional seismic data volume can be established according to the post-stack three-dimensional seismic data by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window.
In an embodiment, a seismic interpretation horizon of post-stack three-dimensional seismic data may be selected; interpolation is carried out on points in the seismic interpretation horizon of the post-stack three-dimensional seismic data, so that the seismic interpretation horizon of the post-stack three-dimensional seismic data is obtained; based on the post-stack three-dimensional seismic data, a post-stack three-dimensional seismic data volume is established by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window. By the steps, a post-stack three-dimensional seismic data volume taking a seismic interpretation layer as a time window (top and bottom) can be obtained.
In the implementation, after the three-dimensional seismic data volume after the stack is built, gray scale calculation can be performed on the three-dimensional seismic data volume after the stack to obtain a gray scale data volume.
In an embodiment, a maximum amplitude value, a minimum amplitude value and an amplitude value of each point in the post-stack three-dimensional seismic data body can be obtained first; and carrying out gray calculation on the post-stack three-dimensional seismic data body according to the maximum amplitude value, the minimum amplitude value and the amplitude value of each point in the post-stack three-dimensional seismic data body to obtain a gray data body.
In an embodiment, the gray data volume may be obtained as follows:
Wherein B is a gray data body, A value for each point within the gray data volume B; /(I)For the amplitude value of each point in the post-stack three-dimensional seismic data body, A max is the maximum amplitude value in the post-stack three-dimensional seismic data body, and A min is the minimum amplitude value in the post-stack three-dimensional seismic data body; c is a scaling factor, which is a constant different from 0; x i、yj is the horizontal and vertical coordinate values of each point, and z k is the time coordinate value of each point along the direction of the seismic trace, i, j, k epsilon N *.
In the implementation, the gain amplification can be performed on the post-stack three-dimensional seismic data body according to the gray data body to obtain a gain data body.
In an embodiment, a gray value for each point in the gray data volume may be obtained; and carrying out gain amplification on the amplitude value of each point in the three-dimensional seismic data body after the stacking according to the gray value of the point in the gray data body corresponding to the point to obtain a gain data body.
In an embodiment, the gain data volume may be obtained as follows:
Wherein P is the gain data volume, For the value of each point in the gain data volume,/>Amplitude value for each point in post-stack three-dimensional seismic data volume,/>For the gray value of each point in the gray data body, x i、yj is the horizontal and vertical coordinate values of each point, z k is the time coordinate value of each point along the direction of the seismic trace, i, j, k epsilon N *.
In the implementation, the oil gas distribution prediction result can be obtained according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body.
In an embodiment, the relative change of amplitude energy between the gain data volume and the post-stack three-dimensional seismic data volume can be obtained according to the amplitude value of each point in the gain data volume and the amplitude value of each point in the post-stack three-dimensional seismic data volume; and obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body.
In an embodiment, the oil and gas distribution prediction result may be obtained according to the following formula:
Wherein E is the oil gas distribution prediction result, Predicting values for oil and gas at each point in the plane; /(I)Is the amplitude value of any point along the direction of a seismic channel, namely the vertical direction, in the post-stack three-dimensional seismic data body,/>Amplitude values for each point along the direction of the seismic trace, i.e., in the vertical direction, within the gain data volume; x i、yj is the horizontal and vertical coordinate values of each point, and z k is the time coordinate value of each point along the direction of the seismic trace, i, j, k epsilon N *.
In an embodiment, the oil gas distribution prediction result may be drawn into a plan view, and the oil gas distribution may be analyzed through the plan view.
The following describes, by way of a specific embodiment, the application effect of the oil-gas distribution prediction method according to the embodiment of the present invention in a certain exploration area:
1. And acquiring post-stack three-dimensional seismic data of the region, and establishing a post-stack three-dimensional seismic data volume according to the post-stack three-dimensional seismic data by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window.
2. And carrying out gray level calculation on the three-dimensional seismic data body after the stacking to obtain a gray level data body. As shown in fig. 2, fig. 2 is a schematic cross-sectional view of a seismic line of a gray-scale data volume B.
3. And carrying out gain amplification on the post-stack three-dimensional seismic data body according to the gray data body to obtain a gain data body. As shown in FIG. 3, FIG. 3 is a schematic cross-sectional view of a seismic line of the gain data volume P.
4. And obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body.
As shown in fig. 4, fig. 4 is a schematic diagram of the oil gas distribution plane characteristics drawn according to the obtained oil gas distribution prediction result. By contrast, the oil and gas distribution plane characteristics of FIG. 4 have a high degree of agreement with the oil and gas distribution characteristics revealed by the well completion in the region.
Based on the same inventive concept, the embodiment of the invention also provides an oil gas distribution prediction device, as described in the following embodiment. Because the principle and the method for solving the problems of the oil gas distribution prediction device are similar, the implementation of the oil gas distribution prediction device can be referred to the implementation of the oil gas distribution prediction method, and the repetition is omitted.
Fig. 5 is a schematic structural diagram of an oil-gas distribution prediction apparatus according to an embodiment of the present invention. As shown in fig. 5, in an embodiment of the present invention, the oil gas distribution prediction apparatus may include:
A post-stack three-dimensional seismic data volume module 501, configured to establish a post-stack three-dimensional seismic data volume according to the post-stack three-dimensional seismic data with a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window;
The gray data body module 502 is configured to perform gray calculation on the post-stack three-dimensional seismic data body to obtain a gray data body;
A gain data body module 503, configured to gain-amplify the post-stack three-dimensional seismic data body according to the gray data body, to obtain a gain data body;
the hydrocarbon distribution prediction module 504 is configured to obtain a hydrocarbon distribution prediction result according to a relative change of amplitude energy between the gain data volume and the post-stack three-dimensional seismic data volume.
In an embodiment, post-stack three-dimensional seismic data volume module 501 may be specifically configured to:
Selecting a seismic interpretation horizon of the three-dimensional seismic data after the stacking;
Interpolation is carried out on points in the seismic interpretation horizon of the post-stack three-dimensional seismic data, so that the seismic interpretation horizon of the post-stack three-dimensional seismic data is obtained;
based on the post-stack three-dimensional seismic data, a post-stack three-dimensional seismic data volume is established by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window.
In an embodiment, the gray-scale data volume module 502 may be specifically configured to:
obtaining the maximum amplitude value, the minimum amplitude value and the amplitude value of each point in the three-dimensional seismic data body after the stack;
and carrying out gray calculation on the post-stack three-dimensional seismic data body according to the maximum amplitude value, the minimum amplitude value and the amplitude value of each point in the post-stack three-dimensional seismic data body to obtain a gray data body.
In an embodiment, the gain data body module 503 may specifically be used to:
obtaining a gray value of each point in the gray data body;
and carrying out gain amplification on the amplitude value of each point in the three-dimensional seismic data body after the stacking according to the gray value of the point in the gray data body corresponding to the point to obtain a gain data body.
In an embodiment, the hydrocarbon distribution prediction result module 504 may be specifically configured to:
Obtaining the relative change of amplitude energy between the gain data volume and the post-stack three-dimensional seismic data volume according to the amplitude value of each point in the gain data volume and the amplitude value of each point in the post-stack three-dimensional seismic data volume;
and obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body.
The embodiment of the invention also provides computer equipment, which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the oil gas distribution prediction method when executing the computer program.
The embodiment of the invention also provides a computer readable storage medium, and the computer readable storage medium stores a computer program for executing the oil gas distribution prediction method.
In summary, in the embodiment of the present invention, according to the post-stack three-dimensional seismic data, a post-stack three-dimensional seismic data volume is established by using the seismic interpretation layer of the post-stack three-dimensional seismic data as a time window; gray scale calculation is carried out on the three-dimensional seismic data body after the stacking to obtain a gray scale data body; gain amplification is carried out on the three-dimensional seismic data body after the stacking according to the gray data body, so that a gain data body is obtained; and obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body. Compared with the post-stack oil gas distribution prediction method based on spectrum analysis in the prior art, the method provided by the invention is based on the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body, so that the oil gas distribution prediction result is obtained, the result has higher coincidence degree with the actual oil gas distribution characteristics, the adverse effect of continuous low-frequency strong reflection information on oil gas prediction is reduced, the accuracy of oil gas distribution prediction is improved, and effective reference can be provided for well position deployment and decision optimization in oil gas exploration and development.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, or as a computer program product. Thus, the present invention may take the form of an entirely software embodiment. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, and various modifications and variations can be made to the embodiments of the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (9)

1. A method of predicting oil and gas distribution, comprising:
According to the post-stack three-dimensional seismic data, using a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window, and establishing a post-stack three-dimensional seismic data volume;
Gray scale calculation is carried out on the three-dimensional seismic data body after the stacking to obtain a gray scale data body;
gain amplification is carried out on the three-dimensional seismic data body after the stacking according to the gray data body, so that a gain data body is obtained;
Obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body;
gain amplification is carried out on the three-dimensional seismic data body after the stacking according to the gray data body, and the gain data body is obtained, and the method comprises the following steps:
obtaining a gray value of each point in the gray data body;
Gain amplification is carried out on the amplitude value of each point in the three-dimensional seismic data body after the stacking according to the gray value of the point in the gray data body corresponding to the point, so as to obtain a gain data body;
gain data volume is obtained according to the following formula:
Wherein P is the gain data volume, For the value of each point in the gain data volume,/>Amplitude value for each point in post-stack three-dimensional seismic data volume,/>For the gray value of each point in the gray data body, x i、yj is the horizontal coordinate value and the vertical coordinate value of each point, and z k is the time coordinate value of each point along the direction of the seismic trace, i, j, k epsilon N *;
obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data volume and the post-stack three-dimensional seismic data volume, comprising:
Obtaining the relative change of amplitude energy between the gain data volume and the post-stack three-dimensional seismic data volume according to the amplitude value of each point in the gain data volume and the amplitude value of each point in the post-stack three-dimensional seismic data volume;
Obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body;
Obtaining an oil gas distribution prediction result according to the following formula:
Wherein E is the oil gas distribution prediction result, Predicting values for oil and gas at each point in the plane; /(I)Is the amplitude value of any point along the direction of a seismic channel, namely the vertical direction, in the post-stack three-dimensional seismic data body,/>Amplitude values for each point along the direction of the seismic trace, i.e., in the vertical direction, within the gain data volume; x i、yj is the horizontal and vertical coordinate values of each point, and z k is the time coordinate value of each point along the direction of the seismic trace, i, j, k epsilon N *.
2. The method of claim 1, wherein creating a post-stack three-dimensional seismic data volume from the post-stack three-dimensional seismic data using a seismic interpretation level of the post-stack three-dimensional seismic data as a time window comprises:
Selecting a seismic interpretation horizon of the three-dimensional seismic data after the stacking;
Interpolation is carried out on points in the seismic interpretation horizon of the post-stack three-dimensional seismic data, so that the seismic interpretation horizon of the post-stack three-dimensional seismic data is obtained;
based on the post-stack three-dimensional seismic data, a post-stack three-dimensional seismic data volume is established by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window.
3. The method of claim 1, wherein performing a gray scale calculation on the post-stack three-dimensional seismic data volume to obtain a gray scale data volume comprises:
obtaining the maximum amplitude value, the minimum amplitude value and the amplitude value of each point in the three-dimensional seismic data body after the stack;
and carrying out gray calculation on the post-stack three-dimensional seismic data body according to the maximum amplitude value, the minimum amplitude value and the amplitude value of each point in the post-stack three-dimensional seismic data body to obtain a gray data body.
4. A method as claimed in claim 3, characterized in that the gray data volume is obtained according to the following formula:
Wherein B is a gray data body, A value for each point within the gray data volume B; /(I)For the amplitude value of each point in the post-stack three-dimensional seismic data body, A max is the maximum amplitude value in the post-stack three-dimensional seismic data body, and A min is the minimum amplitude value in the post-stack three-dimensional seismic data body; c is a scaling factor, which is a constant different from 0; x i、yj is the horizontal and vertical coordinate values of each point, and z k is the time coordinate value of each point along the direction of the seismic trace, i, j, k epsilon N *.
5. An oil and gas distribution prediction device, comprising:
The post-stack three-dimensional seismic data volume module is used for establishing a post-stack three-dimensional seismic data volume by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window according to the post-stack three-dimensional seismic data;
the gray data body module is used for carrying out gray calculation on the three-dimensional seismic data body after the stacking to obtain a gray data body;
the gain data body module is used for carrying out gain amplification on the three-dimensional seismic data body after the stacking according to the gray data body to obtain a gain data body;
The oil-gas distribution prediction module is used for obtaining an oil-gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body;
The gain data volume module is specifically configured to:
obtaining a gray value of each point in the gray data body;
Gain amplification is carried out on the amplitude value of each point in the three-dimensional seismic data body after the stacking according to the gray value of the point in the gray data body corresponding to the point, so as to obtain a gain data body;
gain data volume is obtained according to the following formula:
Wherein P is the gain data volume, For the value of each point in the gain data volume,/>Amplitude value for each point in post-stack three-dimensional seismic data volume,/>For the gray value of each point in the gray data body, x i、yj is the horizontal coordinate value and the vertical coordinate value of each point, and z k is the time coordinate value of each point along the direction of the seismic trace, i, j, k epsilon N *;
the oil gas distribution prediction module is specifically configured to:
Obtaining the relative change of amplitude energy between the gain data volume and the post-stack three-dimensional seismic data volume according to the amplitude value of each point in the gain data volume and the amplitude value of each point in the post-stack three-dimensional seismic data volume;
Obtaining an oil gas distribution prediction result according to the relative change of amplitude energy between the gain data body and the post-stack three-dimensional seismic data body;
Obtaining an oil gas distribution prediction result according to the following formula:
Wherein E is the oil gas distribution prediction result, Predicting values for oil and gas at each point in the plane; /(I)Is the amplitude value of any point along the direction of a seismic channel, namely the vertical direction, in the post-stack three-dimensional seismic data body,/>Amplitude values for each point along the direction of the seismic trace, i.e., in the vertical direction, within the gain data volume; x i、yj is the horizontal and vertical coordinate values of each point, and z k is the time coordinate value of each point along the direction of the seismic trace, i, j, k epsilon N *.
6. The apparatus of claim 5, wherein the post-stack three-dimensional seismic data volume module is configured to:
Selecting a seismic interpretation horizon of the three-dimensional seismic data after the stacking;
Interpolation is carried out on points in the seismic interpretation horizon of the post-stack three-dimensional seismic data, so that the seismic interpretation horizon of the post-stack three-dimensional seismic data is obtained;
based on the post-stack three-dimensional seismic data, a post-stack three-dimensional seismic data volume is established by taking a seismic interpretation layer of the post-stack three-dimensional seismic data as a time window.
7. The apparatus of claim 5, wherein the gray data volume module is specifically configured to:
obtaining the maximum amplitude value, the minimum amplitude value and the amplitude value of each point in the three-dimensional seismic data body after the stack;
and carrying out gray calculation on the post-stack three-dimensional seismic data body according to the maximum amplitude value, the minimum amplitude value and the amplitude value of each point in the post-stack three-dimensional seismic data body to obtain a gray data body.
8. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of claims 1 to 4 when executing the computer program.
9. A computer readable storage medium, characterized in that the computer readable storage medium stores a computer program for executing the method of any one of claims 1 to 4.
CN202010639251.XA 2020-07-06 2020-07-06 Oil-gas distribution prediction method and device Active CN113900141B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010639251.XA CN113900141B (en) 2020-07-06 2020-07-06 Oil-gas distribution prediction method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010639251.XA CN113900141B (en) 2020-07-06 2020-07-06 Oil-gas distribution prediction method and device

Publications (2)

Publication Number Publication Date
CN113900141A CN113900141A (en) 2022-01-07
CN113900141B true CN113900141B (en) 2024-05-28

Family

ID=79186537

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010639251.XA Active CN113900141B (en) 2020-07-06 2020-07-06 Oil-gas distribution prediction method and device

Country Status (1)

Country Link
CN (1) CN113900141B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117784225A (en) * 2022-09-22 2024-03-29 中国石油天然气股份有限公司 Synchronous detection method and device for river channel and fault

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103869358A (en) * 2014-02-18 2014-06-18 中国石油天然气股份有限公司 Histogram equalization-based fault identification method and equipment
CN105572740A (en) * 2016-02-23 2016-05-11 北京海思派克科技有限公司 Apparatus and method for controlling gain of earthquake data

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009039422A1 (en) * 2007-09-21 2009-03-26 Headwave, Inc. Seismic data processing and visualization

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103869358A (en) * 2014-02-18 2014-06-18 中国石油天然气股份有限公司 Histogram equalization-based fault identification method and equipment
CN105572740A (en) * 2016-02-23 2016-05-11 北京海思派克科技有限公司 Apparatus and method for controlling gain of earthquake data

Also Published As

Publication number Publication date
CN113900141A (en) 2022-01-07

Similar Documents

Publication Publication Date Title
CN101334483B (en) Method for attenuating rayleigh wave scattered noise in earthquake data-handling
CN104101902B (en) Seismic properties clustering method and device
CN106896406A (en) The method and apparatus that High-quality Reservoir is predicted based on impedance Domain Properties
CN106443774A (en) Method for improving pre-stack depth migration imaging precision of irregular earth surface
CN105425289A (en) Method and device of determining low frequency wave impedance
CN110879412A (en) Underground transverse wave velocity inversion method, device, computing equipment and storage medium
CN103869362A (en) Method and equipment for obtaining body curvature
CN110632652A (en) Method and device for determining oil-gas migration mode of fault block oil-gas reservoir
CN108873065B (en) Sandstone high-quality reservoir prediction method and device
CN113552621B (en) Shale gas ground stress determination method and device
CN111766635A (en) Sand body communication degree analysis method and system
CN105986819B (en) The method and apparatus with integrated interpretation are automatically processed for well-log information
CN109298464A (en) Tight sandstone reservoir Diagenetic Facies Logging Identification Method and device
CN108629459A (en) The detection method and device of the hydrocarbonaceous hole of reservoir
CN102253414B (en) Reservoir detecting method based on analysis of earthquake lines
CN113900141B (en) Oil-gas distribution prediction method and device
Cho et al. Quasi 3D transdimensional Markov-chain Monte Carlo for seismic impedance inversion and uncertainty analysis
CN109471162A (en) Interbed multiple processing method, system, electronic equipment and readable medium
CN114861515A (en) Method, device, equipment and medium for calculating layer speed data volume
CN112684498A (en) Reservoir fracture prediction method and system based on wide-azimuth seismic data
CN111381280A (en) Method and device for predicting hydrocarbon saturation of reservoir
CN115880455A (en) Three-dimensional intelligent interpolation method based on deep learning
CN109581498B (en) Reservoir thickness distribution determination method, system, device and readable medium
CN113050191A (en) Shale oil TOC prediction method and device based on double parameters
CN105467447A (en) Phased tendency energy matching seismic amplitude preservation evaluation method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant