CN113109867A - Seismic attribute extraction method and device for wedge-shaped stratum braided water channel - Google Patents
Seismic attribute extraction method and device for wedge-shaped stratum braided water channel Download PDFInfo
- Publication number
- CN113109867A CN113109867A CN202010030153.6A CN202010030153A CN113109867A CN 113109867 A CN113109867 A CN 113109867A CN 202010030153 A CN202010030153 A CN 202010030153A CN 113109867 A CN113109867 A CN 113109867A
- Authority
- CN
- China
- Prior art keywords
- seismic
- data
- wedge
- water channel
- tracking
- 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.)
- Granted
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 23
- 238000004590 computer program Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims 2
- 239000004576 sand Substances 0.000 abstract description 21
- 239000000284 extract Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 16
- 238000011161 development Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 230000008021 deposition Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 201000003373 familial cold autoinflammatory syndrome 3 Diseases 0.000 description 2
- 238000005755 formation reaction Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000007480 spreading Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 230000008033 biological extinction Effects 0.000 description 1
- 230000032798 delamination Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 230000003828 downregulation Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000003827 upregulation Effects 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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/308—Time lapse or 4D effects, e.g. production related effects to the formation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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)
- Fluid Mechanics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method and a device for extracting seismic attributes of a wedge-shaped stratum braided water channel, wherein the method comprises the following steps: acquiring seismic data of a wedge-shaped stratum braided water channel; generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data; generating horizon tracking data according to the seismic data of the wedge-shaped stratigraphic braided water channel, a preset tracking range and a preset tracking mode; and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data. The embodiment of the invention effectively extracts the seismic attribute of the wedge-shaped stratum braided water channel and improves the extraction precision of the seismic attribute, thereby finely depicting the water channel sand body.
Description
Technical Field
The invention relates to the technical field of geophysical exploration seismic data processing, in particular to a method and a device for extracting seismic attributes of a wedge-shaped stratum braided water channel.
Background
In the stratum of the lithologic oil and gas reservoir of the braided river or the braided river delta, the braided water channel has the deposition characteristics of quick transverse lithofacies change and frequent longitudinal spatial and temporal evolution. Sandstone bodies on seismic sections are generally characterized by strong amplitude, because the sand bodies generally sediment under strong hydrodynamic conditions and the mudstone deposition environment is quieter.
In order to finely depict the plane spread characteristics of the sand body, plane amplitude attributes are generally extracted, and the attribute extraction has extremely high requirements on stratum interpretation, because the precision of the stratum interpretation determines the precision of the attribute extraction along the stratum and further determines the precision of sand body carving.
The deposition thickness of the wedge-shaped stratum usually shows rapid and uneven change, the vertical evolution of the braided water channel is rapid, the thickness of the sand body is different from several meters to dozens of meters, the attribute extracted by the difference of several milliseconds on seismic data can show obvious difference, the existing seismic attribute extraction method can not achieve the precision, and the water channel sand body is difficult to be finely carved.
Disclosure of Invention
The embodiment of the invention provides a seismic attribute extraction method of a wedge-shaped stratum braided water channel, which is used for extracting seismic attributes of the wedge-shaped stratum braided water channel and improving the extraction precision of the seismic attributes so as to finely depict water channel sand bodies, and comprises the following steps:
acquiring seismic data of a wedge-shaped stratum braided water channel;
generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data;
generating horizon tracking data according to the seismic data of the wedge-shaped stratigraphic braided water channel, a preset tracking range and a preset tracking mode;
and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data.
The embodiment of the invention provides a seismic attribute extraction device of a wedge-shaped stratum braided water channel, which is used for extracting seismic attributes of the wedge-shaped stratum braided water channel and improving the extraction precision of the seismic attributes so as to finely depict water channel sand bodies, and comprises the following components:
the data acquisition module is used for acquiring seismic data of the wedge-shaped stratum braided water channel;
the interpretation data generation module is used for generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data;
the tracking data generation module is used for generating horizon tracking data according to the seismic data of the wedge-shaped stratum braided water channel, a preset tracking range and a preset tracking mode;
and the seismic attribute extraction module is used for extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the seismic attribute extraction method of the wedge-shaped stratigraphic braided water channel when executing the computer program.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the seismic attribute extraction method of the wedge-shaped stratigraphic braided water channel when executing the computer program.
The embodiment of the invention obtains the seismic data of the wedge-shaped stratum braided water channel; generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data; generating horizon tracking data according to the seismic data of the wedge-shaped stratigraphic braided water channel, a preset tracking range and a preset tracking mode; and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data. According to the method and the device, the characteristics of longitudinal sparsity of the horizon interpretation data, large interpretation blank window and large longitudinal interpretation density of the horizon tracking data are utilized, and the horizon interpretation data and the horizon tracking data are combined to extract the seismic attribute of the wedge-shaped stratum braided water channel, so that the seismic attribute of the wedge-shaped stratum braided water channel is effectively extracted, the extraction precision of the seismic attribute is improved, and water channel sand bodies are finely carved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a schematic diagram of a seismic attribute extraction method for a wedge-shaped stratigraphic braided water channel in an embodiment of the invention;
FIG. 2 is a structural diagram of a seismic attribute extraction device of a wedge-shaped stratigraphic braided water channel in an embodiment of the invention;
FIG. 3 is a schematic diagram of X1 well depth calibration in accordance with an embodiment of the present invention;
FIG. 4 is a schematic illustration of horizon-interpreting data in a seismic section in an embodiment of the invention;
FIG. 5 is a diagram illustrating a default tracking mode according to an embodiment of the present invention;
FIG. 6 is a diagram illustrating high density horizon tracking data in accordance with an embodiment of the present invention;
FIG. 7 is a diagram illustrating selection of horizons to be constraints in an embodiment of the present invention;
FIG. 8 is a diagram illustrating the combination of horizon-based interpretation data and horizon-based tracking data in accordance with an embodiment of the present invention;
FIGS. 9-10 are plots of the seismic attributes along the horizon for adjacent troughs and peaks in an embodiment of the present invention;
FIG. 11 is a schematic illustration of the location of a seismic section geodesic 414 along a layer attribute in an embodiment of the invention;
FIG. 12 is a plan view of the attributes along the strata during actual formation development in an embodiment of the present invention;
fig. 13 is a plan-lay layout of a braided waterway resulting from the scribing of a small layer B in an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the embodiments of the present invention are further described in detail below with reference to the accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
As mentioned above, in order to precisely characterize the plane spread of the sand body, the plane amplitude-type attribute is usually extracted, and the attribute extraction requires a very high level of stratigraphic interpretation, because the precision of stratigraphic interpretation determines the precision of extracting the attribute along the layer, and thus determines the precision of sand body carving. Common seismic attribute extraction methods are of four types: firstly, slicing at equal time, and quickly checking the attributes of slices at different time; secondly, slicing along the horizon to obtain the seismic attribute of the position of the interpretation horizon; thirdly, equal-proportion slicing, namely, interpolating a plurality of horizons in equal proportion according to the number of the set slices based on the thickness of the seismic strata of the two sets of interpretation horizons, thereby extracting seismic attributes; and fourthly, extracting seismic attributes by using different time or different interpretation horizon open time windows. The time slicing method is only suitable for a horizontal ideal stratum, seismic attributes of the bedding slices at positions with horizon interpretation are accurate, the stratum is slightly changed or structural positions with horizons are not interpreted, and the seismic attributes of the bedding are not prepared due to the time-through problem. Equal-scale slicing is only suitable for sedimentary formations of uniformly varying thickness. The seismic attributes extracted by setting the time window are comprehensively reflected on the seismic attributes of the sedimentary geologic body within a certain time period, the water channel can be accurately drawn for the stratum with stable sedimentary water channel development, and the sedimentary environment developing in the plaid river and the plaid river delta is not capable of accurately and continuously displaying the plane spread characteristics of the water channel and the process of the evolution of the plane spread characteristics along with the time because the water channel evolves very fast in different periods and the seismic attributes extracted by setting the time window can only be comprehensively reflected on a section of vertical stratum. The deposition thickness of the wedge-shaped stratum usually shows rapid and uneven change, the vertical evolution of the braided water channel is rapid, the thickness of the sand body is different from several meters to dozens of meters, the attribute extracted by the difference of several milliseconds on seismic data can show obvious difference, the existing seismic attribute extraction method can not achieve the precision, and the water channel sand body is difficult to be finely carved.
In order to extract seismic attributes of a wedge-shaped stratigraphic braided water channel and improve the extraction precision of the seismic attributes, so as to finely depict water channel sand bodies, an embodiment of the present invention provides a seismic attribute extraction method of a wedge-shaped stratigraphic braided water channel, as shown in fig. 1, the method may include:
102, generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data;
103, generating horizon tracking data according to the seismic data of the wedge-shaped stratigraphic braided water channel, a preset tracking range and a preset tracking mode;
and 104, extracting the seismic attribute of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data.
As shown in FIG. 1, the embodiment of the invention obtains seismic data of a wedge-shaped stratigraphic braided water channel; generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data; generating horizon tracking data according to the seismic data of the wedge-shaped stratigraphic braided water channel, a preset tracking range and a preset tracking mode; and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data. According to the method and the device, the characteristics of longitudinal sparsity of the horizon interpretation data, large interpretation blank window and large longitudinal interpretation density of the horizon tracking data are utilized, and the horizon interpretation data and the horizon tracking data are combined to extract the seismic attribute of the wedge-shaped stratum braided water channel, so that the seismic attribute of the wedge-shaped stratum braided water channel is effectively extracted, the extraction precision of the seismic attribute is improved, and water channel sand bodies are finely carved.
In specific implementation, seismic data of the wedge-shaped stratum braided water channel are obtained.
In specific implementation, the horizon interpretation data are generated according to the seismic data of the wedge-shaped stratigraphic braided water channel and the pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data.
In the embodiment, the related information and data of drilled wells, layered data, historical seismic data and existing horizon interpretation data are input, well seismic synthetic records are generated according to a conventional synthetic record manufacturing method, and well seismic relations are contrastingly matched with typical geological interface wave group characteristics and transverse well connections, so that synthetic records and seismic goodness of fit are observed. The calibration of the synthetic record can be used for checking the historical horizon interpretation result and determining the position of the target horizon on the seismic time section, so that the research work is developed in a targeted manner.
In an embodiment, if no horizon interpretation data exists, horizon interpretation data is generated according to the seismic data of the wedge stratigraphic braided water channel and a pre-generated well seismic synthetic record, and if the horizon interpretation data exists, no additional generation is needed.
And during specific implementation, generating horizon tracking data according to the seismic data of the wedge-shaped stratum braided water channel, a preset tracking range and a preset tracking mode.
In an embodiment, the preset tracking range includes: and presetting a vertical tracking range and a space sampling point tracking range.
In an embodiment, the preset tracking mode includes: peak and valley tracking modes, and peak, valley and zero phase tracking modes.
In the embodiment, based on the principle of seismic stratigraphy, a preset tracking range and a preset tracking mode of seismic data are designated, the interpretation system automatically matches peaks, troughs and zero phases with the maximum similarity of adjacent seismic props, tracking is carried out, and horizon tracking data are generated.
And during specific implementation, extracting the seismic attribute of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data.
In an embodiment, the extracting seismic attributes of the wedge-shaped stratigraphic braided water channel according to the horizon interpretation data and the horizon tracking data includes: constraining the horizon tracking data by using the horizon interpretation data to obtain the horizon tracking data of the isochronous interface; and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the horizon tracking data of the isochronous interface.
The inventor finds that the longitudinal density of the horizon interpretation data can not meet the requirement of braided water channel identification, the horizon tracking data usually has a local time string phenomenon, and if the horizon interpretation data and the horizon tracking data can be perfectly combined, the dual requirements of horizon fine interpretation, longitudinal upper horizon density and great increase can be met. Judging whether the obtained horizon tracking data under the constraint of the horizon interpretation data is an isochronous interface, if so, considering the horizon result to be reasonable, if still having local penetration, debugging related parameters, or because the constraint force of the horizon interpretation data on the horizon tracking data is insufficient, under the condition, the horizon interpretation data needs to be added, and the constraint on the horizon tracking data is strengthened until the horizon tracking data is the isochronous interface completely. When the horizon meets the condition that peaks, troughs or zero phases on the 3D seismic profile are connected and meet the geological deposition characteristics of the local area, the horizon is considered as an isochronous interface, and the result can be used for seismic attribute extraction. And further, extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the horizon tracking data of the isochronous interface. Based on the change characteristics of the horizon seismic attributes and the pinch-off point position of the stratum of the typical section, the effective range of the horizon seismic attributes is implemented, and the attributes obtained when the horizon seismic attributes and the pinch-off point position are consistent can be used for identifying and finely depicting water channel sand bodies in the wedge-shaped stratum of a plait river or a plait river delta.
Based on the same inventive concept, the embodiment of the invention also provides a device for extracting the seismic attribute of the wedge-shaped stratigraphic braided water channel, which is described in the following embodiment. Because the principles for solving the problems are similar to the seismic attribute extraction method of the wedge-shaped stratum braided water channel, the implementation of the device can refer to the implementation of the method, and repeated parts are not described again.
Fig. 2 is a structural diagram of a seismic property extraction device for a wedge-shaped stratigraphic braided water channel in an embodiment of the invention, and as shown in fig. 2, the device includes:
the data acquisition module 201 is used for acquiring seismic data of a wedge-shaped stratum braided water channel;
an interpretation data generation module 202, configured to generate horizon interpretation data according to the seismic data of the wedge-shaped stratigraphic braided water channel and a pre-generated well seismic synthetic record, where the well seismic synthetic record is pre-generated according to historical seismic data;
the tracking data generation module 203 is used for generating horizon tracking data according to the seismic data of the wedge-shaped stratigraphic braided water channel, a preset tracking range and a preset tracking mode;
and the seismic attribute extraction module 204 is configured to extract the seismic attributes of the wedge-shaped stratigraphic braided water channel according to the layer interpretation data and the layer tracking data.
In one embodiment, the preset tracking range includes: and presetting a vertical tracking range and a space sampling point tracking range.
In one embodiment, the preset tracking mode includes: peak and valley tracking modes, and peak, valley and zero phase tracking modes.
In one embodiment, the seismic attribute extraction module 204 is further configured to:
constraining the horizon tracking data by using the horizon interpretation data to obtain the horizon tracking data of the isochronous interface;
and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the horizon tracking data of the isochronous interface.
In summary, in the embodiments of the present invention, seismic data of a wedge-shaped stratigraphic braided water channel is obtained; generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data; generating horizon tracking data according to the seismic data of the wedge-shaped stratigraphic braided water channel, a preset tracking range and a preset tracking mode; and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data. According to the method and the device, the characteristics of longitudinal sparsity of the horizon interpretation data, large interpretation blank window and large longitudinal interpretation density of the horizon tracking data are utilized, and the horizon interpretation data and the horizon tracking data are combined to extract the seismic attribute of the wedge-shaped stratum braided water channel, so that the seismic attribute of the wedge-shaped stratum braided water channel is effectively extracted, the extraction precision of the seismic attribute is improved, and water channel sand bodies are finely carved.
A specific example is given below to illustrate a specific application of the method for extracting seismic attributes of a wedge-shaped stratigraphic braided water channel in the embodiment of the present invention. In this embodiment, the method first performs synthetic record calibration, verifies or calibrates the target layer wave group characteristics, loads the layer position interpretation data of the target layer position obtained by existing or manual interpretation, starts the intelligent automatic tracking function of the interpretation system to obtain the longitudinal high-density layer position tracking data, and corrects the layer position tracking data of the previous step under the constraint condition of the layer position interpretation data. And finally, extracting high-density and high-precision seismic attributes along the layer under the condition of meeting the isochronous conditions, and further identifying and releasing the water channel sand bodies. The method comprises the following specific steps:
1. making well seismic synthetic records
Firstly, the method needs to utilize the drilled information and data to calibrate the horizon in a research area, and determine the characteristics of different target stratum wave groups and the bright point reflection characteristics and positions of different sand layer groups. As fig. 3 is a deep calibration for the X1 well, indicating that K1tg delamination is a strong valley reflection characteristic, J3q, J2t2 and J2t are all peak displays.
2. Horizon interpretation data
If no horizon interpretation data exists, the characteristics of the main target horizon wave group need to be interpreted on the basis of the calibration of the synthetic record to obtain horizon interpretation data in the seismic profile as shown in FIG. 4, and the interpretation horizon range is the actual development range of the horizon.
3. Horizon tracking data
Starting an automatic tracking function of an interpretation system, according to a tracking mode set in fig. 5, on the premise of no constraint condition, designating an upper time range and a lower time range of a target layer, and tracking a peak, a trough and a zero phase according to a wave characteristic in the range to obtain high-density horizon tracking data as shown in fig. 6. It can be seen that, due to the automatic tracking result, the phenomenon of horizon up-regulation or down-regulation exists locally, and is not completely equal to the time. By considering the characteristics of the wedge-shaped stratum, the automatic tracking can expand the development range of the actual layer according to the extending directions of wave crests and wave troughs, and the condition that the seismic attribute range extracted due to the fact that the ranges of the upper layer and the lower layer of the wedge-shaped stratum are different in the subsequent attribute extraction process can be made up.
4. Combining horizon interpretation data with horizon tracking data
Selecting horizons to be used as constraints, as shown in FIG. 7, results in a combination of horizon-interpretation data and horizon-tracking data as shown in FIG. 8. Because the horizon scheme is explained manually while automatic tracking is carried out, the full-automatic high-density tracking horizon can be adjusted and corrected, and the isochronism is realized on the basis of local non-isochronism of horizon tracking data. If the local inequality exists, whether the horizon interpretation data are completely isochronously closed or not is checked, and the parameters are adjusted until the horizon data are completely isochronously closed, so that the horizon interpretation data can be used for extracting seismic attributes. In FIG. 8, the thickened horizons are horizon-oriented data, the thinner horizons are horizon-oriented data under the constraint of horizon-oriented data, and it can be seen that the trends of the two are substantially the same, and the two are equal in seismic profile.
5. Extracting along-layer seismic attributes
In general, the amplitude-type attribute can better reflect the spreading characteristics of water channel sand bodies, taking the maximum amplitude attribute as an example, fig. 9 and fig. 10 are respectively the along-layer seismic attribute maps of adjacent troughs and peaks, the layer position of a bold solid line in fig. 11 is the position of the along-layer attribute of fig. 9 on a seismic section main survey line 414, and the layer position of a bold dotted line is the position of the attribute of fig. 10 on the seismic section main survey line 414, and the two maps not only clearly depict the planar spreading condition of a braided water channel, but also can observe the evolution process of the water channels through the disappearance and appearance of partial water channels and the change of the amplitude intensity of the water channels.
6. Position of sharp extinction point of actual layer
The thin horizon in fig. 8 is the artificial + intelligent tracking horizon, the thick horizon is the artificial interpretation horizon, it can be seen that the artificial + intelligent tracking horizon actual tracking range is greater than the horizon actual development range, therefore, the horizon seismic attribute range extracted based on the artificial + intelligent tracking horizon in fig. 9 and 10 is also greater than the actual stratigraphic attribute range, the west side of the white line is the actual stratigraphic development range, the east side is the range in which the horizon is increased in the (4) th step, because the two horizons of the east side are completely overlapped, the east side of fig. 9 and 10 is also completely consistent with the horizon seismic attribute. Based on the attribute change characteristics of the two sides of the pinch-out line along the stratum attributes and in combination with the positions of the pinch-out points of the stratum of the seismic section, a stratum attribute plane map in the actual stratum development range can be easily obtained, and as shown in fig. 12, the seismic attribute plane range is the actual development range of the stratum, and the black frame range is the research area range.
7. Sand body picture for water channel
According to the seismic attribute plane map of the stratums obtained in fig. 12, all water channels can be easily carved, and different target stratums of water channel plane layouts can be obtained, for example, fig. 13 is a plaited water channel plane superposition layout obtained by B small-layer carving.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A seismic attribute extraction method of a wedge-shaped stratum braided water channel is characterized by comprising the following steps:
acquiring seismic data of a wedge-shaped stratum braided water channel;
generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data;
generating horizon tracking data according to the seismic data of the wedge-shaped stratigraphic braided water channel, a preset tracking range and a preset tracking mode;
and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data.
2. A method for extracting seismic attributes of a wedge stratigraphic braid flume according to claim 1, wherein the preset tracking range comprises: and presetting a vertical tracking range and a space sampling point tracking range.
3. A method for extracting seismic attributes of a wedge stratigraphic braid flume according to claim 1, wherein the preset tracking mode comprises: peak and valley tracking modes, and peak, valley and zero phase tracking modes.
4. A method for extracting seismic attributes of a wedge stratigraphic braided water channel according to claim 1, wherein the extracting seismic attributes of a wedge stratigraphic braided water channel from the horizon interpretation data and the horizon tracking data comprises:
constraining the horizon tracking data by using the horizon interpretation data to obtain the horizon tracking data of the isochronous interface;
and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the horizon tracking data of the isochronous interface.
5. A seismic attribute extraction device of wedge stratum braided water channel, characterized by comprising:
the data acquisition module is used for acquiring seismic data of the wedge-shaped stratum braided water channel;
the interpretation data generation module is used for generating horizon interpretation data according to the seismic data of the wedge-shaped stratum braided water channel and a pre-generated well seismic synthetic record, wherein the well seismic synthetic record is pre-generated according to historical seismic data;
the tracking data generation module is used for generating horizon tracking data according to the seismic data of the wedge-shaped stratum braided water channel, a preset tracking range and a preset tracking mode;
and the seismic attribute extraction module is used for extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the layer interpretation data and the layer tracking data.
6. A wedge stratigraphic braid water course seismic attribute extraction device of claim 5 wherein the preset tracking range comprises: and presetting a vertical tracking range and a space sampling point tracking range.
7. A wedge stratigraphic braid flume seismic attribute extraction device according to claim 5 wherein the preset tracking mode comprises: peak and valley tracking modes, and peak, valley and zero phase tracking modes.
8. A wedge stratigraphic braid water course seismic attribute extraction device according to claim 5, wherein the seismic attribute extraction module is further configured to:
constraining the horizon tracking data by using the horizon interpretation data to obtain the horizon tracking data of the isochronous interface;
and extracting the seismic attributes of the wedge-shaped stratum braided water channel according to the horizon tracking data of the isochronous interface.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 4 when executing the computer program.
10. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010030153.6A CN113109867B (en) | 2020-01-13 | 2020-01-13 | Seismic attribute extraction method and device for wedge-shaped stratum plait-shaped water channel |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010030153.6A CN113109867B (en) | 2020-01-13 | 2020-01-13 | Seismic attribute extraction method and device for wedge-shaped stratum plait-shaped water channel |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113109867A true CN113109867A (en) | 2021-07-13 |
| CN113109867B CN113109867B (en) | 2023-06-30 |
Family
ID=76709818
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010030153.6A Active CN113109867B (en) | 2020-01-13 | 2020-01-13 | Seismic attribute extraction method and device for wedge-shaped stratum plait-shaped water channel |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113109867B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010096759A1 (en) * | 2009-02-23 | 2010-08-26 | Schlumberger Canada Limited | Dynamic extraction of extrema-based geometric primitives in 3d seismic voxel volume |
| CN102520446A (en) * | 2012-01-06 | 2012-06-27 | 中国石油化工股份有限公司胜利油田分公司西部新区研究中心 | High-precision automatic identification system for three-dimensional sequence |
| CN102981182A (en) * | 2012-11-20 | 2013-03-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 2D seismic data all-horizon automatic tracking method based on unsupervised classification |
| CN103135136A (en) * | 2011-11-25 | 2013-06-05 | 中国石油化工股份有限公司 | Automatic fault interpretation device for three-dimensional seismic data body |
| CN109709603A (en) * | 2018-11-23 | 2019-05-03 | 中国石油天然气股份有限公司 | Seismic horizon identification and method for tracing, system |
-
2020
- 2020-01-13 CN CN202010030153.6A patent/CN113109867B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010096759A1 (en) * | 2009-02-23 | 2010-08-26 | Schlumberger Canada Limited | Dynamic extraction of extrema-based geometric primitives in 3d seismic voxel volume |
| CN103135136A (en) * | 2011-11-25 | 2013-06-05 | 中国石油化工股份有限公司 | Automatic fault interpretation device for three-dimensional seismic data body |
| CN102520446A (en) * | 2012-01-06 | 2012-06-27 | 中国石油化工股份有限公司胜利油田分公司西部新区研究中心 | High-precision automatic identification system for three-dimensional sequence |
| CN102981182A (en) * | 2012-11-20 | 2013-03-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 2D seismic data all-horizon automatic tracking method based on unsupervised classification |
| CN109709603A (en) * | 2018-11-23 | 2019-05-03 | 中国石油天然气股份有限公司 | Seismic horizon identification and method for tracing, system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113109867B (en) | 2023-06-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Eberli et al. | Seismic Imaging of Carbonate Reservoirs and Systems: AAPG Memoir 81 | |
| CN109870719B (en) | Well position layout method, device and system for carbonate rock tight thin reservoir | |
| EP3368921B1 (en) | Interactive salt model modification | |
| CN110632658B (en) | Lateral sealing analysis method and device for fault | |
| US11885922B2 (en) | Adaptive tracking of geological objects | |
| Schneider et al. | Interpretation of fractured zones using seismic attributes—Case study from Teapot Dome, Wyoming, USA | |
| Imam et al. | Mapping the geological structures in the Ras El Ush field (Gulf of Suez, Egypt), based on seismic interpretation and 3D modeling techniques | |
| CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
| CN112180468B (en) | Gravel rock lithofacies prediction method under control of isochronous stratum model | |
| Spina et al. | Characterization of the Devonian Kharyaga carbonate platform (Russia): Integrated and multiscale approach | |
| CN114462113A (en) | Concealed river depicting method and device based on amplitude ratio, electronic equipment and medium | |
| CN113109867A (en) | Seismic attribute extraction method and device for wedge-shaped stratum braided water channel | |
| Gorbunov et al. | SAI-advanced solution to geologic objectives | |
| Zheng et al. | Geological interpretation of well test analysis: a case study from a fluvial reservoir in the Gulf of Thailand | |
| CN113885096B (en) | High-frequency layer sequence division and small layer comparison method and device, electronic equipment and medium | |
| CN112433248B (en) | Method for detecting hidden reservoir stratum in carbonate rock deposition environment | |
| CN113945993B (en) | Beach reservoir prediction method and device | |
| Asheibi | Seismic-based paleoenvironmental analysis of the Paleocene carbonate shelf in Ajdabiya Trough, north-central of Libya | |
| CN114859414B (en) | Method and device for automatically extracting stratum dip angle information from seismic data | |
| Al Musawi et al. | Study of different geostatistical methods to model formation porosity (Cast study of Zubair formation in Luhais oil field) | |
| Bouanga et al. | Implications of HorizonCubes in shallow hazards interpretation | |
| CN111596348A (en) | Stratum velocity correction method and device based on acoustic time difference logging data | |
| Meek | Applications of 3D seismic attribute analysis workflows: a case study from Ness County, Kansas, USA | |
| CN109752760B (en) | Method and device for determining reliability of seismic attribute | |
| Dirstein et al. | Insights from the automated extraction of surfaces from the Bunda 3D seismic survey |
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 |