CN111077570B - Method, device and storage medium for determining invasion position of diabase in oil reservoir - Google Patents
Method, device and storage medium for determining invasion position of diabase in oil reservoir Download PDFInfo
- Publication number
- CN111077570B CN111077570B CN201811231734.5A CN201811231734A CN111077570B CN 111077570 B CN111077570 B CN 111077570B CN 201811231734 A CN201811231734 A CN 201811231734A CN 111077570 B CN111077570 B CN 111077570B
- Authority
- CN
- China
- Prior art keywords
- seismic
- diabase
- oil reservoir
- determining
- invasion
- 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
Links
- 230000009545 invasion Effects 0.000 title claims abstract description 199
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005553 drilling Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 9
- 230000002093 peripheral effect Effects 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 9
- 239000004576 sand Substances 0.000 description 9
- 238000004891 communication Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 230000000644 propagated effect Effects 0.000 description 6
- 230000033001 locomotion Effects 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052655 plagioclase feldspar Inorganic materials 0.000 description 1
- 229910052611 pyroxene Inorganic materials 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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 discloses a method, a device and a storage medium for determining an invasion position of a diabase body in an oil reservoir, and belongs to the technical field of oil field geophysical prospecting. The method comprises the following steps: the method comprises the steps of obtaining amplitude data of each of a plurality of seismic channels in a target oil reservoir, determining an amplitude ratio value set of each seismic channel according to the amplitude data of each of the plurality of seismic channels, and determining an invasion position of diabase in the target oil reservoir according to the amplitude data and the amplitude ratio value set of each of the plurality of seismic channels. The invasion position of the diabase in the target oil reservoir is determined through the amplitude data and the amplitude ratio value set, and the problem that the result error of the determined invasion position of the diabase interval is large due to the fact that the difference between the amplitude of the diabase and the amplitude of the oil, gas or water-containing sandstone at the top-bottom interface of the diabase is small when a person skilled in the relevant art determines the invasion position of the diabase according to an oil reservoir seismic amplitude profile display picture is solved.
Description
Technical Field
The invention relates to the technical field of oil field geophysical prospecting, in particular to a method, a device and a storage medium for determining an invasion position of a diabase body in an oil reservoir.
Background
Diabase is a basic shallow invasion rock, often invaded in sand mudstone, whose main mineral components are pyroxene and basic plagioclase. Because diabase is developed more and has more perfect fracture-cavity units and larger porosity, and is a good reservoir body of oil and gas, the invasion position of diabase in an oil reservoir needs to be determined in the geological exploration process so as to be used as a basis for searching the oil and gas reservoir in the oil reservoir.
In the related art, since the propagation velocity of the seismic wave in diabase is higher than the propagation velocity in oil-, gas-or water-free sand mudstone at the top-bottom interface of diabase, the amplitude of the seismic wave in diabase is larger than the amplitude of the seismic wave in sand mudstone, based on this, when determining the invasion position of diabase in an oil reservoir by using seismic data, a technician may generally find an interval in which the amplitude of the seismic wave is larger than the amplitude of the seismic wave at an adjacent position in an oil reservoir seismic amplitude profile display graph drawn by seismic amplitude, determine the interval as a diabase interval, and determine the position of the diabase interval as the invasion position of diabase.
However, when the sandstone at the top-bottom interface of diabase is oil, gas or water, since the propagation velocity of the seismic wave in the oil, gas or water-containing sandstone is greater than that of the oil, gas or water-free sandstone, the amplitude of the seismic wave in the oil, gas or water-containing sandstone is also increased, so that the difference between the amplitudes of the seismic wave in the diabase and the oil, gas or water-containing sandstone is reduced, the intrusion position of the diabase cannot be accurately determined, and the determination result error is large.
Disclosure of Invention
The embodiment of the invention provides a method, a device and a storage medium for determining an invasion position of diabase in an oil reservoir, which can be used for solving the problem of larger result error of the invasion position of diabase determined according to the amplitude of seismic waves. The technical scheme is as follows:
in a first aspect, there is provided a method of determining the location of invasion of diabase in a reservoir, the method comprising:
acquiring amplitude data of each seismic channel in a plurality of seismic channels in a target oil reservoir, wherein the plurality of seismic channels are obtained by detecting seismic waves transmitted in the target oil reservoir at a plurality of positions of the target oil reservoir;
determining an amplitude ratio value set of each seismic channel according to the amplitude data of each seismic channel in the plurality of seismic channels, wherein the amplitude ratio value set of each seismic channel comprises the amplitude ratio of a plurality of seismic wave periods of each seismic channel, and the amplitude ratio of each seismic wave period refers to the ratio of the amplitude value of the peak position and the amplitude value of the trough position of each seismic wave period;
and determining the invasion position of diabase in the target oil reservoir according to the amplitude data and the amplitude ratio value set of each seismic channel in the plurality of seismic channels.
Optionally, the determining, according to the amplitude data and the amplitude ratio value set of each of the plurality of seismic traces, an invasion position of diabase in the target oil reservoir includes:
determining a positive polarity seismic display section of the target oil reservoir according to the amplitude data of each seismic channel in the plurality of seismic channels;
determining the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil deposit according to the positive polarity earthquake display section and the amplitude ratio value set of each seismic channel in the plurality of seismic channels;
and determining the invasion position of the diabase in the target oil reservoir according to the top boundary position of the invasion position of the diabase in the target oil reservoir and the bottom boundary position of the invasion position of the diabase in the target oil reservoir.
Optionally, the determining, according to the positive polarity seismic display section and the amplitude ratio value set of each of the plurality of seismic traces, a top boundary position and a bottom boundary position of an invasion position of diabase in the target oil reservoir includes:
for each seismic trace in the plurality of seismic traces, if the absolute value of the difference between the amplitude ratio of a first seismic wave period starting from the top of the target oil reservoir in each seismic trace and 1 is smaller than a first threshold, determining the first seismic wave period as the top boundary position period of each seismic trace;
if the amplitude ratio of the k +1 th seismic wave period starting from the top of the target oil deposit in each seismic channel is smaller than the amplitude ratio of the k +1 th seismic wave period, and the absolute value of the difference between the amplitude ratio of the k +1 th seismic wave period and 1 is smaller than the first threshold, determining the k +1 th seismic wave period as the top boundary position period of each seismic channel, wherein k is a positive integer;
and determining the top boundary position and the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the top boundary position periods of the plurality of seismic channels.
Optionally, the determining a top boundary position and a bottom boundary position of an invasion position of diabase in the target oil reservoir according to the top boundary position cycles of the plurality of seismic traces includes:
marking the peak positions of the top boundary position cycles of the seismic channels in the positive polarity seismic display section to obtain a plurality of top boundary position points of the invasion position of diabase;
determining the top boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of top boundary position points of the invasion position of the diabase;
marking the last 0 phase position of the top boundary position period of the seismic channels in the positive polarity seismic display section to obtain a plurality of bottom boundary position points of the invasion position of the diabase;
and determining the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of bottom boundary position points of the invasion position of the diabase.
Optionally, after determining the invasion position of diabase in the target reservoir according to the amplitude data of each of the plurality of seismic traces and the plurality of amplitude ratio values, the method further includes:
acquiring a wellhead position of a new well to be deployed in the target oil reservoir;
determining a target point of a new well to be deployed in the target oil reservoir according to the invasion position of diabase in the target oil reservoir, wherein the target point of the new well to be deployed refers to the position of the new well to be deployed, which needs to be drilled to meet diabase in the drilling process;
and determining the borehole trajectory of the new well to be deployed in the target oil reservoir according to the wellhead position and the target point of the new well to be deployed in the target oil reservoir.
Optionally, the determining a target point of a new well to be deployed in the target oil reservoir according to the invasion position of diabase in the target oil reservoir includes:
determining the oil-gas containing position in the diabase according to the invasion position of the diabase in the target oil reservoir and the amplitude data of each seismic channel in the plurality of seismic channels;
and determining a target point of a new well to be deployed in the target oil reservoir according to the oil-gas-containing position in the diabase.
In a second aspect, there is provided an apparatus for determining the location of invasion of diabase in a reservoir, the apparatus comprising:
the device comprises a first acquisition module, a second acquisition module and a third acquisition module, wherein the first acquisition module is used for acquiring amplitude data of each seismic channel in a plurality of seismic channels in a target oil reservoir, and the plurality of seismic channels are obtained by detecting seismic waves transmitted in the target oil reservoir at a plurality of positions of the target oil reservoir;
a first determining module, configured to determine an amplitude ratio set of each seismic trace according to amplitude data of each seismic trace in the multiple seismic traces, where the amplitude ratio set of each seismic trace includes amplitude ratios of multiple seismic wave periods of each seismic trace, and the amplitude ratio of each seismic wave period refers to a ratio of an amplitude value of a peak position to an amplitude value of a trough position of each seismic wave period;
and the second determination module is used for determining the invasion position of diabase in the target oil reservoir according to the amplitude data and the amplitude ratio value set of each seismic channel in the plurality of seismic channels.
Optionally, the second determining module includes:
the first determining unit is used for determining a positive polarity earthquake display section of the target oil deposit according to the amplitude data of each seismic channel in the plurality of seismic channels;
the second determining unit is used for determining the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil deposit according to the positive polarity earthquake display section and the amplitude ratio value set of each seismic channel in the plurality of seismic channels;
and the third determining unit is used for determining the invasion position of the diabase in the target oil reservoir according to the top boundary position of the invasion position of the diabase in the target oil reservoir and the bottom boundary position of the invasion position of the diabase in the target oil reservoir.
Optionally, the second determining unit includes:
the first determining subunit is configured to determine, for each of the plurality of seismic traces, a first seismic wave period as a top boundary position period of each of the seismic traces if an absolute value of a difference between an amplitude ratio of the first seismic wave period starting from the top of the target oil reservoir in each of the seismic traces and 1 is smaller than a first threshold;
a second determining subunit, configured to determine, if an amplitude ratio of a k +1 th seismic wave period from the top of the target oil reservoir in each seismic trace is smaller than the amplitude ratio of the k +1 th seismic wave period, and an absolute value of a difference between the amplitude ratio of the k +1 th seismic wave period and 1 is smaller than the first threshold, the k +1 th seismic wave period as a top boundary position period of each seismic trace, where k is a positive integer;
and the third determining subunit is used for determining the top boundary position and the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the top boundary position periods of the plurality of seismic channels.
Optionally, the second determining unit includes:
the first marking subunit is used for marking the peak positions of the top boundary position cycles of the seismic channels in the positive polarity seismic display section so as to obtain a plurality of top boundary position points of the invasion positions of diabase;
the fourth determining subunit is used for determining the top boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of top boundary position points of the invasion position of the diabase;
the second marking subunit is used for marking the last 0 phase position of the top boundary position period of the seismic channels in the positive polarity seismic display section so as to obtain a plurality of bottom boundary position points of the invasion position of the diabase;
and the fifth determining subunit is used for determining the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of bottom boundary position points of the invasion position of the diabase.
Optionally, the apparatus further comprises:
the second acquisition module is used for acquiring the well mouth position of a new well to be deployed in the target oil reservoir;
the third determining module is used for determining a target point of a new well to be deployed in the target oil reservoir according to the invasion position of diabase in the target oil reservoir, wherein the target point of the new well to be deployed refers to the position of the new well to be deployed, which needs to be drilled to meet diabase in the drilling process;
and the fourth determining module is used for determining the well track of the new well to be deployed in the target oil reservoir according to the well mouth position and the target point of the new well to be deployed in the target oil reservoir.
Optionally, the determining a third module comprises:
the fourth determining unit is used for determining the oil-gas-containing position in the diabase according to the invasion position of the diabase in the target oil reservoir and the amplitude data of each seismic channel in the plurality of seismic channels;
and the fifth determining unit is used for determining a target point of a new well to be deployed in the target oil reservoir according to the oil-gas containing position in the diabase.
In a third aspect, there is provided an apparatus for determining the location of invasion of diabase in a reservoir, the apparatus comprising:
a processor and a memory for storing processor-executable instructions;
wherein the processor is configured to perform any of the methods of the first aspect above.
In a fourth aspect, a computer-readable storage medium is provided, in which a computer program is stored, which, when executed by a processor, implements any of the methods provided in the first aspect above.
The technical scheme provided by the embodiment of the invention at least has the following beneficial effects: in the embodiment of the invention, a plurality of amplitude ratios can be determined according to the amplitude data of each of a plurality of seismic channels when the seismic waves in the target oil reservoir are transmitted, and the invasion position of diabase in the target oil reservoir can be determined according to the amplitude data of each of the plurality of seismic channels and the plurality of amplitude ratios. Because the difference of the change of the amplitude ratio value is more obvious than the difference of the change of the amplitude in the seismic data, the boundary position between the diabase and other rocks determined according to the amplitude ratio value is clearer than the boundary position between the diabase and other rocks determined according to the amplitude, on the basis, the accuracy of the determination result can be improved when the intrusion position of the diabase is determined by adopting the amplitude ratio method, and the problem that the result error of the determined intrusion position of the diabase section is larger because the amplitude difference of the diabase and the oil, gas or water-containing sand mudstone at the top and bottom interfaces of the diabase is smaller when the technical personnel in the related art determine the intrusion position of the diabase only according to a plurality of amplitudes in each seismic channel is solved. In addition, because the embodiment of the invention calculates and plots the amplitude data of each seismic channel in the plurality of seismic channels, and then determines the invasion position of the diabase interval according to the calculation result and the plotting result, the method is simple and easy to operate, and the working efficiency of determining the invasion position of the diabase in the oil reservoir is improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a flow chart of a method for determining an invasion location of diabase in a reservoir according to an embodiment of the present invention;
FIG. 2 is a flow chart of another method for determining the invasion location of diabase in a reservoir according to an embodiment of the present invention;
FIG. 3 is a positive polarity seismic display section in a target reservoir provided by an embodiment of the invention;
FIG. 4 is a reverse polarity seismic display profile in a target reservoir provided by an embodiment of the invention;
FIG. 5 is a schematic diagram of an apparatus for determining the invasion position of diabase in an oil reservoir according to an embodiment of the present invention;
fig. 6 is a schematic structural diagram of a terminal 600 according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
Before explaining the embodiments of the present invention in detail, terms related to the embodiments of the present invention will be explained.
Positive polarity seismic display section
The positive polarity earthquake display section refers to the earthquake display section of which the positive reflection coefficient interface corresponds to a wave crest and the negative reflection coefficient interface corresponds to a wave trough.
Phase 0
The 0 phase is a position where the amplitude value is 0 in one seismic wave cycle.
Well trajectory
The wellbore trajectory refers to the specific location of the wellbore axis in the target reservoir during drilling.
Fig. 1 is a flowchart of a method for determining an invasion position of diabase in an oil reservoir according to an embodiment of the present invention. Referring to fig. 1, the method comprises the steps of:
step 101: the method comprises the steps of obtaining amplitude data of each seismic channel in a plurality of seismic channels in a target oil reservoir, wherein the plurality of seismic channels are obtained by detecting seismic waves transmitted in the target oil reservoir at a plurality of positions of the target oil reservoir.
Step 102: and determining an amplitude ratio set of each seismic channel according to the amplitude data of each seismic channel in the plurality of seismic channels, wherein the amplitude ratio set of each seismic channel comprises the amplitude ratios of a plurality of seismic wave periods of each seismic channel, and the amplitude ratio of each seismic wave period refers to the ratio of the amplitude value of the peak position to the amplitude value of the trough position of each seismic wave period.
Step 103: and determining the invasion position of diabase in the target oil reservoir according to the amplitude data and the amplitude ratio value set of each seismic channel in the plurality of seismic channels.
In the embodiment of the invention, a plurality of amplitude ratios can be determined according to the amplitude data of each of a plurality of seismic channels when the seismic waves in the target oil reservoir are propagated, and the invasion position of diabase in the target oil reservoir is determined according to the amplitude data of each of the plurality of seismic channels and the plurality of amplitude ratios, so that the problem that the result error of the determined invasion position of the diabase interval is large due to the fact that the amplitude difference between the amplitude of the diabase and the amplitude of the oil, gas or water-containing sandstone at the top-bottom interface of the diabase is small when the invasion position of the diabase is determined according to the seismic amplitude profile display image of the oil reservoir by a person skilled in the relevant technology is solved.
Optionally, determining an invasion position of diabase in the target oil reservoir according to the amplitude data and the amplitude ratio value set of each of the plurality of seismic channels, including:
determining a positive polarity earthquake display section of the target oil reservoir according to the amplitude data of each of the plurality of seismic channels;
determining the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil deposit according to the positive earthquake display section and the amplitude ratio value set of each seismic channel in the plurality of seismic channels;
and determining the invasion position of the diabase in the target oil reservoir according to the top boundary position of the invasion position of the diabase in the target oil reservoir and the bottom boundary position of the invasion position of the diabase in the target oil reservoir.
Optionally, determining a top boundary position and a bottom boundary position of an invasion position of diabase in the target oil reservoir according to the positive polarity seismic display section and the amplitude ratio value set of each seismic channel in the plurality of seismic channels, including:
for each seismic trace in the plurality of seismic traces, if the absolute value of the difference value between the amplitude ratio of the first seismic wave period starting from the top of the target oil reservoir in each seismic trace and 1 is smaller than a first threshold value, determining the first seismic wave period as the top boundary position period of each seismic trace;
if the amplitude ratio of the k +1 seismic wave period starting from the top of the target oil reservoir in each seismic channel is smaller than the amplitude ratio of the k seismic wave period, and the absolute value of the difference between the amplitude ratio of the k +1 seismic wave period and 1 is smaller than a first threshold, determining the k +1 seismic wave period as the top boundary position period of each seismic channel, wherein k is a positive integer;
and determining the top boundary position and the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the top boundary position periods of the plurality of seismic channels.
Optionally, determining a top boundary position and a bottom boundary position of an invasion position of diabase in the target oil reservoir according to the top boundary position cycle of the plurality of seismic channels, including:
marking the peak positions of the top boundary position cycles of a plurality of seismic channels in the positive polarity seismic display section to obtain a plurality of top boundary position points of the invasion position of diabase;
determining the top boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of top boundary position points of the invasion position of the diabase;
marking the last 0 phase position of the top boundary position period of a plurality of seismic channels in the positive polarity seismic display section to obtain a plurality of bottom boundary position points of the invasion position of diabase;
and determining the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of bottom boundary position points of the invasion position of the diabase.
Optionally, after determining the invasion position of the diabase in the target oil reservoir according to the amplitude data and the amplitude ratios of each of the plurality of seismic channels, the method further includes:
acquiring a wellhead position of a new well to be deployed in a target oil reservoir;
determining a target point of a new well to be deployed in the target oil reservoir according to the invasion position of diabase in the target oil reservoir, wherein the target point of the new well to be deployed refers to the position of the new well to be deployed, which needs to be drilled to meet diabase in the drilling process;
and determining the borehole trajectory of the new well to be deployed in the target oil reservoir according to the wellhead position and the target point of the new well to be deployed in the target oil reservoir.
Optionally, determining a target point of a new well to be deployed in the target oil reservoir according to the invasion position of diabase in the target oil reservoir, including:
determining the oil-gas-containing position in the diabase according to the invasion position of the diabase in the target oil reservoir and the amplitude data of each seismic channel in the plurality of seismic channels;
and determining a target point of a new well to be deployed in the target oil reservoir according to the oil-gas-containing position in the diabase.
All the above optional technical solutions can be combined arbitrarily to form an optional embodiment of the present invention, which is not described in detail herein.
Fig. 2 is a flowchart of another method for determining an invasion position of diabase in an oil reservoir according to an embodiment of the present invention, where the method may be applied to a terminal, which may be a mobile phone, a tablet computer, a computer, or the like. Referring to fig. 2, the method comprises the steps of:
step 201: and acquiring the amplitude data of each seismic channel in a plurality of seismic channels in the target oil reservoir.
The multiple seismic channels are obtained by detecting seismic waves propagated in the target oil reservoir at multiple positions of the target oil reservoir, that is, the multiple seismic channels are seismic wave data detected at different positions of the target oil reservoir. Furthermore, since seismic waves are typically periodically propagated in the target reservoir, each seismic trace includes a plurality of seismic wave periods, and the amplitude data of each seismic trace is composed of the amplitude data of the plurality of seismic wave periods.
In practical application, the amplitude data of each seismic channel in the multiple seismic channels can be obtained by user input, can also be obtained by sending of other equipment, and can also be obtained by time domain seismic data volume analysis of a target oil reservoir. For example, a time domain seismic data volume of a target reservoir is obtained, where the time domain seismic data volume includes original seismic wave amplitude data of each of a plurality of seismic traces for describing a geological structure of the target reservoir and time data when the original seismic wave amplitude data is received, then the time domain seismic data volume is converted into a depth domain seismic data volume through a time domain-depth domain seismic data volume conversion formula, where the depth domain seismic data volume includes seismic wave amplitude data of each of the plurality of seismic traces for characterizing the target reservoir and formation depth data corresponding to each of the original seismic wave amplitude data, and finally, the amplitude data of each of the plurality of seismic traces in the obtained depth domain seismic data volume is determined as the amplitude data of each of the plurality of seismic traces in the target reservoir.
It should be noted that the time domain-depth domain seismic data volume conversion formula required for converting the time domain seismic data volume into the depth domain seismic data volume is common knowledge in the art, and the specific form of the time domain-depth domain seismic data volume conversion formula may be selected by a skilled person according to the actual characteristics of the target reservoir and the processing capability of the terminal, which is not specifically limited in this embodiment of the present invention.
In addition, the reservoir depth in the embodiment of the present invention refers to the shortest distance from the position of any point in the target reservoir to the top surface of the target reservoir, which is measured by using the top surface of the target reservoir as a reference surface. Since the distance is positive, the reservoir depth is also positive.
Step 202: and determining an amplitude ratio value set of each seismic channel according to the amplitude data of each seismic channel in the plurality of seismic channels.
The amplitude ratio value set of each seismic channel comprises the amplitude ratio values of a plurality of seismic wave periods of each seismic channel, and the amplitude ratio value of each seismic wave period refers to the ratio of the amplitude value of the peak position to the amplitude value of the trough position of each seismic wave period.
Specifically, after obtaining the amplitude data of each of the plurality of seismic channels in the target oil reservoir, the terminal may analyze the amplitude data of each of the plurality of seismic channels in the target oil reservoir to obtain a plurality of seismic wave periods of each seismic channel, an amplitude value at a peak position and an amplitude value at a trough position in each of the plurality of seismic wave periods, and an oil reservoir depth corresponding to positions of three 0 phases arranged in order from the top to the bottom of the oil reservoir in each seismic wave period.
After the amplitude value at the peak position and the amplitude value at the trough position in each seismic wave period in the multiple seismic wave periods of each seismic channel are obtained, determining the ratio of the amplitude value at the peak position and the amplitude value at the trough position in each seismic wave period in the multiple seismic wave periods of each seismic channel as the amplitude ratio of the seismic wave period, and determining the amplitude ratios of the multiple seismic wave periods of each seismic channel as the amplitude ratio set of the seismic channel.
In one possible embodiment, Table 1 shows amplitude data on a seismic trace of a target reservoir, where the first column shows the seismic wave period, the second column shows the amplitude values at the peak position, the third column shows the amplitude values at the trough position, and the fourth column shows the amplitude ratio.
TABLE 1
As can be seen from table 1, three seismic wave periods exist on one seismic trace of the target reservoir, where the amplitude ratio of the first period is 2.41, the amplitude ratio of the second period is 1.06, and the amplitude ratio of the third period is 2.39, and then the set {2.41,1.06,2.39} of the amplitude ratios of the three periods is determined as the amplitude ratio set of the seismic trace.
It should be noted that the amplitude data on the seismic trace of the target reservoir shown in table 1 above are all exemplary data given in the embodiment of the present invention, and do not constitute a specific limitation on the amplitude data on the seismic trace of the target reservoir.
Step 203: and determining the invasion position of diabase in the target oil reservoir according to the amplitude data and the amplitude ratio value set of each seismic channel in the plurality of seismic channels.
Specifically, the invasion position of diabase in the target oil reservoir can be determined through the following steps 2031 to 2033 according to the amplitude ratio value set of each seismic channel in the plurality of seismic channels.
Step 2031: and determining a positive polarity earthquake display section of the target oil deposit according to the amplitude data of each of the plurality of seismic channels.
Specifically, referring to fig. 3, when determining a positive polarity seismic display section of the target oil reservoir, a positive reflection coefficient interface and a negative reflection coefficient interface may be determined, where the positive reflection coefficient interface is a surface composed of a point at a position corresponding to a maximum amplitude in each of the plurality of seismic wave cycles, the negative reflection coefficient interface is a surface composed of a point at a position corresponding to a minimum amplitude in each of the plurality of seismic wave cycles, the positive reflection coefficient interface in the amplitude data of each of the plurality of seismic traces may be mapped to a peak of the positive polarity seismic display section, and the negative reflection coefficient interface is mapped to a trough of the positive polarity seismic display section, so as to obtain the positive polarity seismic display section of the target oil reservoir. In addition, in the positive polarity earthquake display section, the wave crest in each seismic channel in the plurality of seismic channels can be marked to be black, so that a plurality of coaxial continuous wave crests on the same axis in the positive polarity earthquake display section of the target oil deposit can be determined according to the black mark of the wave crest in each seismic channel in the plurality of seismic channels, and the lithologic development rule of the target oil deposit can be observed more clearly and more visually.
Step 2032: and determining the top boundary position and the bottom boundary position of the invasion position of the diabase in the target oil deposit according to the positive earthquake display section and the amplitude ratio value set of each seismic channel in the plurality of seismic channels.
Before determining the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil reservoir according to the positive seismic display section and the amplitude ratio value set of each seismic channel in the plurality of seismic channels, the top boundary position period of each seismic channel can be determined according to the amplitude ratio value set of each seismic channel in the plurality of seismic channels, and the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil reservoir can be determined according to the positive seismic display section and the top boundary position period of each seismic channel.
Optionally, the top boundary position period of each seismic trace may be determined according to the following steps:
(1) and for each seismic channel in the plurality of seismic channels, if the absolute value of the difference value between the amplitude ratio of the first seismic wave period starting from the top of the target oil reservoir in each seismic channel and 1 is smaller than a first threshold value, determining the first seismic wave period as the top boundary position period of each seismic channel.
The first preset threshold may be preset, specifically may be set by a technician according to actual conditions of different oil reservoirs, or may be set by a terminal default, which is not limited in the embodiment of the present invention. For example, the first preset threshold may be set to 0.2.
In a possible embodiment, assuming that the first preset threshold is 0.2, if the amplitude ratio of the first seismic wave period starting from the top of the target oil reservoir in each seismic trace is 1.1, and the absolute value of the difference between the amplitude ratio of the seismic wave period and 1 is 0.1, at this time, the absolute value of the difference between the amplitude ratio of the seismic wave period and 1 is less than the first threshold 0.2, the first seismic wave period starting from the top of the target oil reservoir in each seismic trace may be determined as the top boundary position period of each seismic trace; if the amplitude ratio of the first seismic wave period starting from the top of the target oil reservoir in each seismic channel is 2.5, the absolute value of the difference between the amplitude ratio of the seismic wave period and 1 is 1.5, at this time, the absolute value of the difference between the amplitude ratio of the seismic wave period and 1 is greater than a first threshold value 0.2, and the first seismic wave period starting from the top of the target oil reservoir in each seismic channel is not the top boundary position period of each seismic channel. The first threshold and the amplitude ratio are not particularly limited in the embodiments of the present invention.
(2) And if the amplitude ratio of the k +1 th seismic wave period starting from the top of the target oil reservoir in each seismic channel is smaller than the amplitude ratio of the k +1 th seismic wave period, and the absolute value of the difference between the amplitude ratio of the k +1 th seismic wave period and 1 is smaller than a first threshold, determining the k +1 th seismic wave period as the top boundary position period of each seismic channel, wherein k is a positive integer.
It should be noted that, in the target oil reservoir, the propagation speed of the seismic wave in the diabase is higher than the propagation speed of the seismic wave in the sand mudstone, and therefore, the amplitude value at the peak position and the amplitude value at the trough position in the seismic wave period of the diabase are both greater than the amplitude value at the peak position and the amplitude value at the trough position in the seismic wave period of the sand mudstone, and since the amplitude value at the peak position and the amplitude value at the trough position in the seismic wave period of the diabase are relatively stable, the amplitude ratio of the diabase is close to 1, and due to the layering effect of the sand mudstone, the amplitude ratio of the sand mudstone at the peak position and the amplitude value at the trough position in the seismic wave period of the sand mudstone is greater than 1, and therefore, it is possible to determine that the amplitude ratio of the k +1 seismic wave periods from the top of the target oil reservoir in each seismic channel is greater than 1 And determining the top boundary position period of each seismic channel by using the amplitude ratio smaller than the kth seismic wave period and the absolute value of the difference between the amplitude ratio of the kth +1 seismic wave period and 1 smaller than a first threshold.
In addition, when the amplitude ratio of at least two seismic wave periods existing in one seismic channel meets the judgment condition, the top boundary position period of at least two diabase existing in one seismic channel can be determined.
In one possible embodiment, where table 2 shows amplitude ratios of multiple seismic wave cycles in a seismic trace in a target reservoir, the first threshold may be 0.2. Where the first column shows the number of seismic wave cycles and the second column shows the amplitude ratio.
TABLE 2
Numbering of seismic wave periods | Amplitude ratio |
1 | 2.41 |
2 | 2.2 |
3 | 1.52 |
4 | 1.06 |
5 | 1.22 |
6 | 3.29 |
7 | 3.53 |
As can be seen from Table 2, the amplitude ratio of the 4 th seismic wave period from the top of the target reservoir is smaller than the amplitude ratio of the 3 rd seismic wave period, and the difference between the amplitude ratio of the 4 th seismic wave period and 1 is smaller than 0.2, then the 4 th seismic wave period can be determined as the top boundary position period of the seismic trace.
It should be noted that the amplitude ratio of a plurality of seismic wave periods on a seismic trace of a target reservoir and the first threshold shown in table 2 are exemplary data given in the embodiment of the present invention, and do not constitute a specific limitation on the amplitude ratio of a plurality of seismic wave periods and the first threshold on a seismic trace of a target reservoir.
Step 2033: and determining the top boundary position and the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the top boundary position periods of the plurality of seismic channels.
When determining the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil reservoir, the top boundary position of the invasion position of diabase in the target oil reservoir may be determined, and then the bottom boundary position of the invasion position of diabase in the target oil reservoir may be determined according to the top boundary position of the invasion position of diabase in the target oil reservoir.
Alternatively, the top boundary position of the invasion position of diabase in the target reservoir may be determined according to the following steps:
(1) in the positive polarity earthquake display section, the peak positions of the top boundary position cycles of a plurality of earthquake channels are marked so as to obtain a plurality of top boundary position points of the invasion position of diabase.
(2) And determining the top boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of top boundary position points of the invasion position of the diabase.
It should be noted that, because the positive polarity seismic display section has the peak corresponding to the positive reflection coefficient interface, that is, the peak position in the top boundary position period in the positive polarity seismic display section can be determined as the top boundary position of the invasion position of diabase, after the top boundary position period of each seismic channel is determined, the position of the peak in the top boundary position period of each seismic channel can be marked, so as to identify the top boundary position point on each seismic channel in the positive polarity seismic display section, and the position corresponding to the top boundary position point on each seismic channel is the top boundary position of the invasion position of diabase on each seismic channel.
After the top boundary position points of the invasion position of the diabase are determined, smoothly connecting the top boundary position points to obtain a connecting line of the top boundary position of the invasion position of the diabase, and determining the top boundary position of the invasion position of the diabase on each seismic channel according to the specific position of the connecting line of the top boundary position of the invasion position of the diabase in the positive polarity seismic display section.
After the top boundary position of the invasion position of the diabase on each seismic channel is determined, the bottom boundary position of the invasion position of the diabase on each seismic channel can be determined according to the top boundary position period in the positive polarity seismic display section. And when the bottom boundary position of the invasion position of the diabase on each seismic channel is determined according to the top boundary position period in the positive polarity seismic display section, the bottom boundary position of the invasion position of the diabase on each seismic channel can be directly determined through the top boundary position period in the positive polarity seismic display section, or a reverse polarity seismic display section can be determined according to the positive polarity seismic display section, and then the bottom boundary position of the invasion position of the diabase on each seismic channel can be determined according to the top boundary position period in the reverse polarity seismic display section.
Optionally, in the positive polarity seismic display section, marking the last 0 phase position of the top boundary position cycle of the multiple seismic channels to obtain multiple bottom boundary position points of the invasion position of the diabase; and determining the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of bottom boundary position points of the invasion position of the diabase.
It should be noted that, when diabase develops in the oil reservoir, the bottom boundary position of the invasion position of diabase and the corresponding top boundary position thereof are all in one seismic wave period, that is, the thickness of diabase when invading the oil reservoir is less than one seismic wave period, therefore, after the top boundary position period of a plurality of seismic channels is determined, the last 0 phase position point of the top boundary position period of the plurality of seismic channels can be determined as the bottom boundary position point of the invasion position of diabase of the plurality of seismic channels, the plurality of bottom boundary position points are smoothly connected, that is, the connecting line of the bottom boundary position of the invasion position of diabase can be obtained, and the bottom boundary position of the invasion position of diabase on each seismic channel can be determined according to the specific position of the connecting line of the bottom boundary position of diabase in the positive polarity seismic display section.
Optionally, referring to fig. 4, in order to more clearly and more intuitively show the bottom boundary position of the invasion position of the diabase on each seismic channel, when determining the bottom boundary position of the invasion position of the diabase on each seismic channel, the polarity of the seismic wave period on each seismic channel in the positive polarity seismic display section may be reversed by 180 degrees, so that the peak of the seismic wave period on each seismic channel in the positive polarity seismic display section becomes the trough in the reverse polarity seismic display section, and the trough of the seismic wave period on each seismic channel in the positive polarity seismic display section becomes the peak in the reverse polarity seismic display section
In the reverse polarity seismic display section, although the peak and the trough in the top boundary position period are subjected to polarity conversion by 180 degrees, the relative position of the top boundary position period in the reverse polarity seismic display section is not changed, so that the last 0-phase position point of the top boundary position period of a plurality of seismic channels in the reverse polarity seismic display section can be determined as the bottom boundary position point of the invasion positions of the diabase of the plurality of seismic channels, the bottom boundary position points are smoothly connected, the connecting line of the bottom boundary position of the invasion position of the diabase can be obtained, and the bottom boundary position of the invasion position of the diabase on each seismic channel can be determined according to the specific position of the connecting line of the bottom boundary position of the invasion position of the diabase in the reverse polarity seismic display section. In addition, in order to enable the bottom boundary position of the invasion position of the diabase in the target oil reservoir to be clearer in the reversed polarity earthquake display section, wave crests in all earthquake wave periods in a plurality of earthquake channels in the reversed polarity earthquake display section can be marked to be black, so that the bottom boundary position of the invasion position of the diabase in the target oil reservoir can be observed more clearly and more visually according to the black marks of the wave crests in all earthquake wave periods in the plurality of earthquake channels.
Step 2034: and determining the invasion position of the diabase in the target oil reservoir according to the top boundary position of the invasion position of the diabase in the target oil reservoir and the bottom boundary position of the invasion position of the diabase in the target oil reservoir.
After the top boundary position according to the invasion position of the diabase in the target oil reservoir and the bottom boundary position according to the invasion position of the diabase in the target oil reservoir are determined, the position between the top boundary position of the invasion position of the diabase in the target oil reservoir and the bottom boundary position of the invasion position of the diabase in the target oil reservoir can be determined as the invasion position of the diabase in the target oil reservoir.
Step 204: and determining a target point of a new well to be deployed in the target oil reservoir according to the invasion position of diabase in the target oil reservoir, wherein the target point of the new well to be deployed refers to the position of the new well to be deployed, which needs to be drilled to meet diabase in the drilling process.
Because diabase is developed more and has more perfect fracture-cavity units and porosity is larger, the diabase is a good reservoir body of oil gas and is a geological basis for determining the target point of the new well to be deployed, the target point of the new well to be deployed can be further determined after the invasion position of the diabase in the target oil deposit is determined, and the drilling success rate of the new well to be deployed in the target oil deposit is improved.
When determining a target point of a new well to be deployed in a target oil reservoir, determining the oil-gas-containing position in diabase according to the invasion position of diabase in the target oil reservoir and the amplitude data of each seismic channel in a plurality of seismic channels; and determining a target point of a new well to be deployed in the target oil reservoir according to the oil-gas-containing position in the diabase.
Specifically, because the oil-gas storage probability of the upper layer position in the invasion position of the diabase in the target oil reservoir is greater than that of the lower layer position, after the invasion position of the diabase in the target oil reservoir is determined, the invasion position of the diabase can be further layered according to the amplitude data of each seismic channel in a plurality of seismic channels in the target oil reservoir to obtain the invasion position of the upper layer diabase and the invasion position of the lower layer diabase, wherein the invasion position of the upper layer diabase is the oil-gas position in the diabase, and any position in the invasion position of the upper layer diabase can be determined as a target point of a new well to be deployed in the target oil reservoir.
In a possible embodiment, if 7 drilled positions in a target oil reservoir are assumed to be drilled all at diabase, diabase reservoir yield is high, oil test daily yield is up to 22t/d, and initial daily yield after production is up to 33t/d, therefore, the method provided by the embodiment of the invention can be used for determining the diabase invasion position of the target oil reservoir, and after the diabase invasion position in the target oil reservoir is determined, new wells are drilled in the target oil reservoir, and the drilling result is shown in table 3. Where the first column of table 3 shows the well number, the second column shows the amplitude value at the peak position of the diabase, the third column shows the amplitude value at the trough position of the diabase, the fourth column shows the amplitude ratio, the fifth column shows the target depth, the sixth column shows the actual drilling depth, and the seventh column shows the error rate.
TABLE 3
As can be seen from table 3, after the invasion position of diabase is determined according to the method provided by the embodiment of the present invention, the error rate between the target point of the new well to be deployed and the depth of the target point encountered by the actual drill is further determined to be not greater than 1%, which indicates that the accuracy of determining the target point of the new well to be deployed according to the invasion position of diabase determined by the method provided by the embodiment of the present invention is higher.
Step 205: and acquiring the well mouth position of the new well to be deployed in the target oil reservoir, and determining the well bore track of the new well to be deployed in the target oil reservoir according to the well mouth position and the target point of the new well to be deployed in the target oil reservoir.
It should be noted that, because of the influence of factors such as geological environment and government land use management policy, the well mouth position of the new well to be deployed cannot be directly deployed at the ground position vertically corresponding to the determined target point of the new well to be deployed, so that after the target point of the new well to be deployed in the target oil reservoir is determined, the well bore track of the new well to be deployed in the target oil reservoir can be determined according to the obtained well mouth position and the obtained target point position of the new well to be deployed in the target oil reservoir.
For example, when the method is not limited by geological environment and government land use management policies, the wellhead position of a new well to be deployed in the target oil reservoir can be selected on the ground position vertically corresponding to the target point, and the well track at the moment is a vertical line; when the method is limited by geological environment and government land use management policies, after the wellhead position of the new well to be deployed in the target oil reservoir is determined, the optimal borehole trajectory can be determined according to the wellhead position and the target point of the new well to be deployed in the target oil reservoir, so that the borehole trajectory is scientific and reasonable, and the drilling efficiency is also improved.
In the embodiment of the invention, a plurality of amplitude ratios can be determined according to the amplitude data of each of a plurality of seismic channels when the seismic waves in the target oil reservoir are propagated, and the invasion position of diabase in the target oil reservoir is determined according to the amplitude data of each of the plurality of seismic channels and the plurality of amplitude ratios, so that the problem that the result error of the determined invasion position of the diabase interval is large due to the fact that the amplitude difference between the amplitude of the diabase and the amplitude of the oil, gas or water-containing sandstone at the top-bottom interface of the diabase is small when the invasion position of the diabase is determined according to the seismic amplitude profile display image of the oil reservoir by a person skilled in the relevant technology is solved. In addition, because the embodiment of the invention calculates and draws the amplitude data of each seismic channel in the plurality of seismic channels, and determines the invasion position of the diabase interval according to the obtained calculation result and drawing result, the method is simple and easy to operate, and improves the working efficiency of determining the invasion position of the diabase in the oil reservoir.
In addition, after the invasion position of the diabase in the target oil reservoir is determined, the target point of the new well to be deployed in the target oil reservoir can be determined according to the invasion position of the diabase in the target oil reservoir, the determined target point is high in accuracy, the success rate of drilling is improved, and the determined well track can improve the drilling efficiency according to the well mouth position and the target point of the new well to be deployed in the target oil reservoir.
Fig. 5 is a schematic structural diagram of an apparatus for determining an invasion position of diabase in an oil reservoir according to an embodiment of the present invention. Referring to fig. 5, the apparatus may include:
the first obtaining module 501 is configured to obtain amplitude data of each of a plurality of seismic traces in a target oil reservoir, where the plurality of seismic traces are obtained by detecting seismic waves propagated in the target oil reservoir at a plurality of positions of the target oil reservoir.
The first determining module 502 is configured to determine an amplitude ratio set of each seismic trace according to amplitude data of each seismic trace in the plurality of seismic traces, where the amplitude ratio set of each seismic trace includes amplitude ratios of a plurality of seismic wave periods of each seismic trace, and the amplitude ratio of each seismic wave period refers to a ratio of an amplitude value of a peak position to an amplitude value of a trough position of each seismic wave period.
And a second determining module 503, configured to determine an invasion position of diabase in the target oil reservoir according to the amplitude data and the amplitude ratio value set of each of the plurality of seismic channels.
Optionally, the second determining module includes:
the first determining unit is used for determining a positive polarity earthquake display section of the target oil deposit according to the amplitude data of each seismic channel in the plurality of seismic channels;
the second determining unit is used for determining the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil deposit according to the positive earthquake display section and the amplitude ratio value set of each seismic channel in the plurality of seismic channels;
and the third determining unit is used for determining the invasion position of the diabase in the target oil reservoir according to the top boundary position of the invasion position of the diabase in the target oil reservoir and the bottom boundary position of the invasion position of the diabase in the target oil reservoir.
Optionally, the second determination unit includes:
the first determining subunit is used for determining a first seismic wave period as the top boundary position period of each seismic channel if the absolute value of the difference between the amplitude ratio of the first seismic wave period starting from the top of the target oil reservoir in each seismic channel and 1 is smaller than a first threshold value for each seismic channel;
the second determining subunit is used for determining the k +1 th seismic wave period as the top boundary position period of each seismic channel if the amplitude ratio of the k +1 th seismic wave period starting from the top of the target oil reservoir in each seismic channel is smaller than the amplitude ratio of the k +1 th seismic wave period and the absolute value of the difference between the amplitude ratio of the k +1 th seismic wave period and 1 is smaller than a first threshold, wherein k is a positive integer;
and the third determining subunit is used for determining the top boundary position and the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the top boundary position periods of the plurality of seismic channels.
Optionally, the second determination unit includes:
the first marking subunit is used for marking the peak positions of the top boundary position cycles of a plurality of seismic channels in the positive polarity seismic display section so as to obtain a plurality of top boundary position points of the invasion positions of diabase;
the fourth determining subunit is used for determining the top boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of top boundary position points of the invasion position of the diabase;
the second marking subunit is used for marking the last 0 phase position of the top boundary position period of the plurality of seismic channels in the positive polarity seismic display section so as to obtain a plurality of bottom boundary position points of the invasion position of the diabase;
and the fifth determining subunit is used for determining the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of bottom boundary position points of the invasion position of the diabase.
Optionally, the apparatus further comprises:
the second acquisition module is used for acquiring the well mouth position of a new well to be deployed in the target oil reservoir;
the third determining module is used for determining a target point of a new well to be deployed in the target oil reservoir according to the invasion position of diabase in the target oil reservoir, wherein the target point of the new well to be deployed refers to the position of the new well to be deployed, which needs to be drilled to meet diabase in the drilling process;
and the fourth determining module is used for determining the well track of the new well to be deployed in the target oil deposit according to the well mouth position and the target point of the new well to be deployed in the target oil deposit.
Optionally, determining the third module comprises:
the fourth determining unit is used for determining the oil-gas containing position in the diabase according to the invasion position of the diabase in the target oil reservoir and the amplitude data of each seismic channel in the plurality of seismic channels;
and the fifth determining unit is used for determining a target point of a new well to be deployed in the target oil reservoir according to the oil-gas containing position in the diabase.
In the embodiment of the invention, a plurality of amplitude ratios can be determined according to the amplitude data of each of a plurality of seismic channels when the seismic waves in the target oil reservoir are propagated, and the invasion position of diabase in the target oil reservoir is determined according to the amplitude data of each of the plurality of seismic channels and the plurality of amplitude ratios, so that the problem that the result error of the determined invasion position of the diabase interval is large due to the fact that the amplitude difference between the amplitude of the diabase and the amplitude of the oil, gas or water-containing sandstone at the top-bottom interface of the diabase is small when the invasion position of the diabase is determined according to the seismic amplitude profile display image of the oil reservoir by a person skilled in the relevant technology is solved. In addition, because the embodiment of the invention calculates and draws the amplitude data of each seismic channel in the plurality of seismic channels, and determines the invasion position of the diabase interval according to the obtained calculation result and drawing result, the method is simple and easy to operate, and improves the working efficiency of determining the invasion position of the diabase in the oil reservoir. After the invasion position of the diabase in the target oil reservoir is determined, the target point of the new well to be deployed in the target oil reservoir can be determined according to the invasion position of the diabase in the target oil reservoir, the determined target point is high in accuracy, the success rate of well drilling is improved, and the well track can be determined according to the well mouth position and the target point of the new well to be deployed in the target oil reservoir, so that the well drilling efficiency can be improved.
It should be noted that: the device for determining an invasion position of diabase in an oil reservoir provided in the above embodiment is illustrated by only the division of the above functional modules when determining an invasion position of diabase in an oil reservoir, and in practical applications, the above function allocation may be completed by different functional modules according to needs, that is, the internal structure of the device may be divided into different functional modules to complete all or part of the above described functions. In addition, the device for determining the invasion position of the diabase in the oil reservoir and the method embodiment for determining the invasion position of the diabase in the oil reservoir provided by the embodiments belong to the same concept, and specific implementation processes are detailed in the method embodiment and are not described again.
Fig. 6 is a schematic structural diagram of a terminal 600 according to an embodiment of the present invention. The terminal 600 may be: a smart phone, a tablet computer, an MP3 player (Moving Picture Experts Group Audio Layer III, motion video Experts compression standard Audio Layer 3), an MP4 player (Moving Picture Experts Group Audio Layer IV, motion video Experts compression standard Audio Layer 4), a notebook computer, or a desktop computer. The terminal 600 may also be referred to by other names such as user equipment, portable terminal, laptop terminal, desktop terminal, etc.
In general, the terminal 600 includes: a processor 601 and a memory 602.
The processor 601 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and so on. The processor 601 may be implemented in at least one hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), and a PLA (Programmable Logic Array). The processor 601 may also include a main processor and a coprocessor, where the main processor is a processor for Processing data in an awake state, and is also called a Central Processing Unit (CPU); a coprocessor is a low power processor for processing data in a standby state. In some embodiments, the processor 601 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content required to be displayed on the display screen. In some embodiments, processor 601 may also include an AI (Artificial Intelligence) processor for processing computational operations related to machine learning.
The memory 602 may include one or more computer-readable storage media, which may be non-transitory. The memory 602 may also include high-speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 602 is used to store at least one instruction for execution by processor 601 to implement the method of determining the invasion location of diabase in a reservoir provided by the method embodiments herein.
In some embodiments, the terminal 600 may further optionally include: a peripheral interface 603 and at least one peripheral. The processor 601, memory 602, and peripheral interface 603 may be connected by buses or signal lines. Various peripheral devices may be connected to the peripheral interface 603 via a bus, signal line, or circuit board. Specifically, the peripheral device includes: at least one of a radio frequency circuit 604, a touch screen display 604, a camera 606, an audio circuit 607, a positioning component 608, and a power supply 609.
The peripheral interface 603 may be used to connect at least one peripheral related to I/O (Input/Output) to the processor 601 and the memory 602. In some embodiments, the processor 601, memory 602, and peripheral interface 603 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 601, the memory 602, and the peripheral interface 603 may be implemented on a separate chip or circuit board, which is not limited in this embodiment.
The Radio Frequency circuit 604 is used for receiving and transmitting RF (Radio Frequency) signals, also called electromagnetic signals. The radio frequency circuitry 604 communicates with communication networks and other communication devices via electromagnetic signals. The rf circuit 604 converts an electrical signal into an electromagnetic signal to transmit, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 604 comprises: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and so forth. The radio frequency circuitry 604 may communicate with other terminals via at least one wireless communication protocol. The wireless communication protocols include, but are not limited to: metropolitan area networks, various generation mobile communication networks (2G, 3G, 4G, and 4G), Wireless local area networks, and/or WiFi (Wireless Fidelity) networks. In some embodiments, the rf circuit 604 may further include NFC (Near Field Communication) related circuits, which are not limited in this application.
The display screen 604 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display screen 605 is a touch display screen, the display screen 605 also has the ability to capture touch signals on or over the surface of the display screen 605. The touch signal may be input to the processor 601 as a control signal for processing. At this point, the display 605 may also be used to provide virtual buttons and/or a virtual keyboard, also referred to as soft buttons and/or a soft keyboard. In some embodiments, the display 605 may be one, providing the front panel of the terminal 600; in other embodiments, the display 605 may be at least two, respectively disposed on different surfaces of the terminal 600 or in a folded design; in still other embodiments, the display 605 may be a flexible display disposed on a curved surface or on a folded surface of the terminal 600. Even more, the display 605 may be arranged in a non-rectangular irregular pattern, i.e., a shaped screen. The Display 605 may be made of LCD (Liquid Crystal Display), OLED (Organic Light-Emitting Diode), and the like.
The camera assembly 606 is used to capture images or video. Optionally, camera assembly 606 includes a front camera and a rear camera. Generally, a front camera is disposed at a front panel of the terminal, and a rear camera is disposed at a rear surface of the terminal. In some embodiments, the number of the rear cameras is at least two, and each rear camera is any one of a main camera, a depth-of-field camera, a wide-angle camera and a telephoto camera, so that the main camera and the depth-of-field camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize panoramic shooting and VR (Virtual Reality) shooting functions or other fusion shooting functions. In some embodiments, camera assembly 606 may also include a flash. The flash lamp can be a monochrome temperature flash lamp or a bicolor temperature flash lamp. The double-color-temperature flash lamp is a combination of a warm-light flash lamp and a cold-light flash lamp, and can be used for light compensation at different color temperatures.
The positioning component 608 is used for positioning the current geographic Location of the terminal 600 to implement navigation or LBS (Location Based Service). The Positioning component 608 can be a Positioning component based on the united states GPS (Global Positioning System), the chinese beidou System, the russian graves System, or the european union's galileo System.
In some embodiments, the terminal 600 also includes one or more sensors 610. The one or more sensors 610 include, but are not limited to: acceleration sensor 611, gyro sensor 612, pressure sensor 613, fingerprint sensor 614, optical sensor 615, and proximity sensor 616.
The acceleration sensor 611 may detect the magnitude of acceleration in three coordinate axes of the coordinate system established with the terminal 600. For example, the acceleration sensor 611 may be used to detect components of the gravitational acceleration in three coordinate axes. The processor 601 may control the touch screen display 605 to display the user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 611. The acceleration sensor 611 may also be used for acquisition of motion data of a game or a user.
The gyro sensor 612 may detect a body direction and a rotation angle of the terminal 600, and the gyro sensor 612 and the acceleration sensor 611 may cooperate to acquire a 3D motion of the user on the terminal 600. The processor 601 may implement the following functions according to the data collected by the gyro sensor 612: motion sensing (such as changing the UI according to a user's tilting operation), image stabilization at the time of photographing, game control, and inertial navigation.
The pressure sensor 613 may be disposed on a side frame of the terminal 600 and/or on a lower layer of the touch display screen 605. When the pressure sensor 613 is disposed on the side frame of the terminal 600, a user's holding signal of the terminal 600 can be detected, and the processor 601 performs left-right hand recognition or shortcut operation according to the holding signal collected by the pressure sensor 613. When the pressure sensor 613 is disposed at the lower layer of the touch display screen 605, the processor 601 controls the operability control on the UI interface according to the pressure operation of the user on the touch display screen 605. The operability control comprises at least one of a button control, a scroll bar control, an icon control and a menu control.
The fingerprint sensor 614 is used for collecting a fingerprint of a user, and the processor 601 identifies the identity of the user according to the fingerprint collected by the fingerprint sensor 614, or the fingerprint sensor 614 identifies the identity of the user according to the collected fingerprint. Upon identifying that the user's identity is a trusted identity, the processor 601 authorizes the user to perform relevant sensitive operations including unlocking the screen, viewing encrypted information, downloading software, paying, and changing settings, etc. The fingerprint sensor 614 may be disposed on the front, back, or side of the terminal 600. When a physical button or vendor Logo is provided on the terminal 600, the fingerprint sensor 614 may be integrated with the physical button or vendor Logo.
The optical sensor 615 is used to collect the ambient light intensity. In one embodiment, processor 601 may control the display brightness of touch display 605 based on the ambient light intensity collected by optical sensor 615. Specifically, when the ambient light intensity is high, the display brightness of the touch display screen 605 is increased; when the ambient light intensity is low, the display brightness of the touch display screen 605 is turned down. In another embodiment, the processor 601 may also dynamically adjust the shooting parameters of the camera assembly 606 according to the ambient light intensity collected by the optical sensor 615.
A proximity sensor 616, also known as a distance sensor, is typically disposed on the front panel of the terminal 600. The proximity sensor 616 is used to collect the distance between the user and the front surface of the terminal 600. In one embodiment, when the proximity sensor 616 detects that the distance between the user and the front surface of the terminal 600 gradually decreases, the processor 601 controls the touch display 605 to switch from the bright screen state to the dark screen state; when the proximity sensor 616 detects that the distance between the user and the front surface of the terminal 600 gradually becomes larger, the processor 601 controls the touch display 605 to switch from the breath screen state to the bright screen state.
That is, not only is an embodiment of the present invention provide a terminal including a processor and a memory for storing processor-executable instructions, where the processor is configured to execute the method in the embodiment shown in fig. 1 or fig. 2, but also an embodiment of the present invention provides a computer-readable storage medium having a computer program stored therein, where the computer program, when executed by the processor, can implement the method for determining an invasion position of diabase in an oil reservoir in the embodiment shown in fig. 1 or fig. 2.
Those skilled in the art will appreciate that the configuration shown in fig. 6 is not intended to be limiting of terminal 600 and may include more or fewer components than those shown, or some components may be combined, or a different arrangement of components may be used.
It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program instructing relevant hardware, where the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (7)
1. A method of determining the invasion location of diabase in a reservoir, the method comprising:
acquiring amplitude data of each seismic channel in a plurality of seismic channels in a target oil reservoir, wherein the plurality of seismic channels are obtained by detecting seismic waves transmitted in the target oil reservoir at a plurality of positions of the target oil reservoir;
determining an amplitude ratio value set of each seismic channel according to the amplitude data of each seismic channel in the plurality of seismic channels, wherein the amplitude ratio value set of each seismic channel comprises the amplitude ratio of a plurality of seismic wave periods of each seismic channel, and the amplitude ratio of each seismic wave period refers to the ratio of the amplitude value of the peak position and the amplitude value of the trough position of each seismic wave period;
determining a positive polarity seismic display section of the target oil reservoir according to the amplitude data of each seismic channel in the plurality of seismic channels; for each seismic trace in the plurality of seismic traces, if the absolute value of the difference between the amplitude ratio of a first seismic wave period starting from the top of the target oil reservoir in each seismic trace and 1 is smaller than a first threshold, determining the first seismic wave period as the top boundary position period of each seismic trace; if the amplitude ratio of the k +1 th seismic wave period starting from the top of the target oil deposit in each seismic channel is smaller than the amplitude ratio of the k +1 th seismic wave period, and the absolute value of the difference between the amplitude ratio of the k +1 th seismic wave period and 1 is smaller than the first threshold, determining the k +1 th seismic wave period as the top boundary position period of each seismic channel, wherein k is a positive integer; determining the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil reservoir according to the top boundary position periods of the plurality of seismic channels; and determining the invasion position of the diabase in the target oil reservoir according to the top boundary position of the invasion position of the diabase in the target oil reservoir and the bottom boundary position of the invasion position of the diabase in the target oil reservoir.
2. The method of claim 1, wherein determining a top boundary position and a bottom boundary position of an invasion position of diabase in the target reservoir from a top boundary position period of the plurality of seismic traces comprises:
marking the peak positions of the top boundary position cycles of the seismic channels in the positive polarity variable density seismic display section to obtain a plurality of top boundary position points of the invasion position of diabase;
determining the top boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of top boundary position points of the invasion position of the diabase;
marking the last 0 phase position of the top boundary position period of the seismic channels in the positive polarity variable density seismic display section to obtain a plurality of bottom boundary position points of the invasion position of the diabase;
and determining the bottom boundary position of the invasion position of the diabase in the target oil reservoir according to the plurality of bottom boundary position points of the invasion position of the diabase.
3. The method of any one of claims 1-2, wherein after determining the invasion location of diabase in the target reservoir based on the amplitude data and the set of amplitude ratios for each of the plurality of seismic traces, further comprising:
acquiring a wellhead position of a new well to be deployed in the target oil reservoir;
determining a target point of a new well to be deployed in the target oil reservoir according to the invasion position of diabase in the target oil reservoir, wherein the target point of the new well to be deployed refers to the position of the new well to be deployed, which needs to be drilled to meet diabase in the drilling process;
and determining the borehole trajectory of the new well to be deployed in the target oil reservoir according to the wellhead position and the target point of the new well to be deployed in the target oil reservoir.
4. The method of claim 3, wherein determining the target point of a new well to be deployed in the target reservoir based on the invasion location of diabase in the target reservoir comprises:
determining the oil-gas containing position in the diabase according to the invasion position of the diabase in the target oil reservoir and the amplitude data of each seismic channel in the plurality of seismic channels;
and determining a target point of a new well to be deployed in the target oil reservoir according to the oil-gas-containing position in the diabase.
5. An apparatus for determining the location of diabase invasion in a reservoir, the apparatus comprising:
the device comprises a first acquisition module, a second acquisition module and a third acquisition module, wherein the first acquisition module is used for acquiring amplitude data of each seismic channel in a plurality of seismic channels in a target oil reservoir, and the plurality of seismic channels are obtained by detecting seismic waves transmitted in the target oil reservoir at a plurality of positions of the target oil reservoir;
a first determining module, configured to determine an amplitude ratio set of each seismic trace according to amplitude data of each seismic trace in the multiple seismic traces, where the amplitude ratio set of each seismic trace includes amplitude ratios of multiple seismic wave periods of each seismic trace, and the amplitude ratio of each seismic wave period refers to a ratio of an amplitude value of a peak position to an amplitude value of a trough position of each seismic wave period;
the first determining unit is used for determining a positive polarity earthquake display section of the target oil deposit according to the amplitude data of each seismic channel in the plurality of seismic channels; the first determining subunit is used for determining a first seismic wave period as the top boundary position period of each seismic channel if the absolute value of the difference between the amplitude ratio of the first seismic wave period starting from the top of the target oil reservoir in each seismic channel and 1 is smaller than a first threshold value for each seismic channel; the second determining subunit is used for determining the k +1 th seismic wave period as the top boundary position period of each seismic channel if the amplitude ratio of the k +1 th seismic wave period starting from the top of the target oil reservoir in each seismic channel is smaller than the amplitude ratio of the k +1 th seismic wave period and the absolute value of the difference between the amplitude ratio of the k +1 th seismic wave period and 1 is smaller than a first threshold, wherein k is a positive integer; the third determining subunit is used for determining the top boundary position and the bottom boundary position of the invasion position of diabase in the target oil reservoir according to the top boundary position periods of the plurality of seismic channels; and the third determining unit is used for determining the invasion position of the diabase in the target oil reservoir according to the top boundary position of the invasion position of the diabase in the target oil reservoir and the bottom boundary position of the invasion position of the diabase in the target oil reservoir.
6. An apparatus for determining the location of diabase invasion in a reservoir, the apparatus comprising:
a processor and a memory for storing processor-executable instructions;
wherein the processor is configured to perform the method of any one of claims 1-4.
7. A computer-readable storage medium, characterized in that the storage medium has stored therein a computer program which, when being executed by a processor, carries out the method of any one of claims 1-4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811231734.5A CN111077570B (en) | 2018-10-22 | 2018-10-22 | Method, device and storage medium for determining invasion position of diabase in oil reservoir |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811231734.5A CN111077570B (en) | 2018-10-22 | 2018-10-22 | Method, device and storage medium for determining invasion position of diabase in oil reservoir |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111077570A CN111077570A (en) | 2020-04-28 |
CN111077570B true CN111077570B (en) | 2021-11-30 |
Family
ID=70309898
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811231734.5A Active CN111077570B (en) | 2018-10-22 | 2018-10-22 | Method, device and storage medium for determining invasion position of diabase in oil reservoir |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111077570B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11977198B2 (en) | 2020-10-06 | 2024-05-07 | Saudi Arabian Oil Company | Isofrequency volumes ratio workflow to detect gas reservoirs in 3D domain |
US11333780B2 (en) | 2020-10-09 | 2022-05-17 | Saudi Arabian Oil Company | Method and system for processing a three-dimensional (3D) seismic dataset |
US11592589B2 (en) | 2021-01-14 | 2023-02-28 | Saudi Arabian Oil Company | Seismic attribute map for gas detection |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520518A (en) * | 2008-02-25 | 2009-09-02 | 中国石油集团东方地球物理勘探有限责任公司 | Method for recognizing the lithology of petrosilex by using the combined characteristics of gravity-magnetic-electronic anomaly |
CN107728203A (en) * | 2016-08-11 | 2018-02-23 | 中国石油天然气股份有限公司 | Method and device for determining development degree of carbonate reservoir |
CN110469315A (en) * | 2019-08-21 | 2019-11-19 | 中国石油天然气股份有限公司大港油田分公司 | Well spacing determines that method and note adopt array pitch and determine method |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160291183A1 (en) * | 2015-04-01 | 2016-10-06 | Chevron U.S.A. Inc. | System and method for seismic imaging of complex subsurface volumes |
CN105388511A (en) * | 2015-10-16 | 2016-03-09 | 辽宁工程技术大学 | Speed anisotropic microseismic monitoring positioning method, microseismic monitoring positioning terminal and microseismic monitoring positioning system |
-
2018
- 2018-10-22 CN CN201811231734.5A patent/CN111077570B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101520518A (en) * | 2008-02-25 | 2009-09-02 | 中国石油集团东方地球物理勘探有限责任公司 | Method for recognizing the lithology of petrosilex by using the combined characteristics of gravity-magnetic-electronic anomaly |
CN107728203A (en) * | 2016-08-11 | 2018-02-23 | 中国石油天然气股份有限公司 | Method and device for determining development degree of carbonate reservoir |
CN110469315A (en) * | 2019-08-21 | 2019-11-19 | 中国石油天然气股份有限公司大港油田分公司 | Well spacing determines that method and note adopt array pitch and determine method |
Non-Patent Citations (1)
Title |
---|
HZ地区火成岩地震识别、储层预测技术及其应用效果;杨德洲,等;《第三届隐蔽油气藏国际学术研讨会论文集》;20040430;347-357 * |
Also Published As
Publication number | Publication date |
---|---|
CN111077570A (en) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112922589B (en) | Pinch-out line determining method, pinch-out line determining device, terminal and storage medium | |
CN111077570B (en) | Method, device and storage medium for determining invasion position of diabase in oil reservoir | |
CN110174699B (en) | Method, device and storage medium for determining fault closure | |
CN109696706B (en) | Salt structure detection method, device and storage medium | |
CN114970919A (en) | Drilling rate determining method and device, computer equipment and storage medium | |
CN110673214B (en) | Method and device for predicting depths of inlet target point and end point of horizontal well | |
CN111255434A (en) | Well testing method and device for gas well and computer storage medium | |
CN109281664B (en) | Method and device for predicting oil production of horizontal well | |
CN108733902B (en) | Method and device for determining permeability of oil reservoir and storage medium | |
CN111736215B (en) | Fault fault distance determining method and device | |
CN109343113B (en) | Method and device for predicting oil and gas reservoir position | |
CN115542389A (en) | Method and device for predicting formation temperature and computer storage medium | |
CN113931621A (en) | Method and device for determining gas well accumulated liquid information and storage medium | |
CN112698391B (en) | Sand body describing method, device, equipment and storage medium | |
CN113586043B (en) | Method and device for determining saturation parameters of bound water and computer equipment | |
CN112096373B (en) | Method and device for selecting shale gas area and computer storage medium | |
CN112987105B (en) | Method, device, terminal and storage medium for quantitatively predicting underground rock salt layer distribution | |
CN115508894B (en) | Method, device, equipment and storage medium for determining ground grid parameters | |
CN115263244B (en) | Method and device for controlling drainage of air-water layer and computer storage medium | |
CN111088969B (en) | Method and device for determining separate injection scheme of water injection well and storage medium | |
CN110058306B (en) | Method, apparatus and computer storage medium for determining three-dimensional velocity volume | |
CN110630254B (en) | Oil-water layer distinguishing method and device and computer readable storage medium | |
CN114861961A (en) | Method and device for determining development scale parameters of silt falling layer | |
CN115163056A (en) | Method and device for determining dynamic reserves, computer equipment and storage medium | |
CN115708001A (en) | Method, device and equipment for predicting oil and gas reservoir and storage medium |
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 |