CN110873892B - Horizon extraction method and device - Google Patents
Horizon extraction method and device Download PDFInfo
- Publication number
- CN110873892B CN110873892B CN201811013943.2A CN201811013943A CN110873892B CN 110873892 B CN110873892 B CN 110873892B CN 201811013943 A CN201811013943 A CN 201811013943A CN 110873892 B CN110873892 B CN 110873892B
- Authority
- CN
- China
- Prior art keywords
- velocity
- stratum
- longitudinal wave
- speed
- coordinate point
- 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
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. analysis, for interpretation, for correction
Abstract
The invention provides a layer extraction method and a layer extraction device, wherein the method comprises the following steps: acquiring the longitudinal wave speed of a well depth sampling point in the well depth direction; obtaining a stratum interface according to the longitudinal wave speed of the current sampling point or the current stratum thickness; and determining an equal-velocity coordinate point set in a preset region of a stratum interface and taking the equal-velocity coordinate point set as an extracted horizon, wherein the equal-velocity coordinate point set comprises at least one equal-velocity coordinate point, and the difference between the longitudinal wave velocity and the horizon velocity of the equal-velocity coordinate point is within a preset velocity difference threshold interval range. Therefore, the situation that the horizon picking is inaccurate and multiple solutions occur when the wave crest or the wave trough is directly picked from the seismic data is avoided, and more stable depth domain horizon data can be obtained.
Description
Technical Field
The invention relates to the technical field of seismic data interpretation, in particular to a level extraction method and a level extraction device.
Background
Horizon extraction is a key seismic data interpretation technology, and horizon data for structural interpretation at present is mainly obtained from a time domain imaging section. The application of horizon data is not limited to structural interpretation, and is widely applied to seismic inversion. Due to the multi-solution problem of seismic inversion, various known geological data, well logging data and the like are required to be used as constraint conditions. The depth domain horizon data has an important role in data-driven depth domain seismic inversion, and currently, a method for obtaining a depth domain horizon directly picks peaks or troughs from seismic data, but the mode can cause the situations of inaccurate horizon pickup and multiple solutions, so that a new horizon extraction method is needed to be provided to improve the accuracy of horizon extraction.
Disclosure of Invention
The invention provides a horizon extraction method and a horizon extraction device, which are used for solving the technical problem of low accuracy of horizon extraction in the prior art.
The invention provides a layer extraction method on one hand, which comprises the following steps:
acquiring the longitudinal wave speed of a well depth sampling point in the well depth direction;
obtaining stratum interfaces according to the longitudinal wave speed or the current stratum thickness of the current sampling point, wherein the well depth sampling point comprises the current sampling point, the current stratum thickness is the distance from the determined stratum interface closest to the current sampling point in all the determined stratum interfaces;
determining the obtained horizon speed of the stratum interface;
and determining an equal-velocity coordinate point set in a preset region of the stratum interface, and taking the equal-velocity coordinate point set as the extracted layer, wherein the equal-velocity coordinate point set comprises at least one equal-velocity coordinate point, and the difference between the longitudinal wave velocity of the equal-velocity coordinate point and the layer velocity is within a preset velocity difference threshold interval range.
Further, the step of obtaining a formation interface according to the longitudinal wave velocity or the current formation thickness of the current sampling point includes:
acquiring the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point in the well depth direction, and determining the position of the current sampling point as the stratum interface when the difference between the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point and the longitudinal wave velocity of the current sampling point is greater than a first preset threshold value;
or when the thickness of the current stratum is larger than a second preset threshold value, determining the position of the current sampling point as the stratum interface.
Further, the step of determining a set of constant velocity coordinate points in a preset region of the stratum interface and using the set of constant velocity coordinate points as the extracted horizon includes:
obtaining longitudinal wave velocity of a coordinate point to be selected in a preset area of the stratum interface, wherein the coordinate point to be selected is a point obtained by taking the obtained stratum interface or the coordinate point with the equal velocity as a datum point and taking a preset step length as a step length,
and when the difference between the longitudinal wave velocity of the coordinate point to be selected and the horizon velocity is within the range of the preset velocity difference threshold value interval, adding the coordinate point to be selected in the constant velocity coordinate point set, and taking the constant velocity coordinate point set as the extracted horizon.
Further, the step of determining the obtained horizon speed of the formation interface comprises:
acquiring the obtained longitudinal wave speeds of all well depth sampling points in a layer where the stratum interface is located;
and taking the average value of the longitudinal wave speeds of all the well depth sampling points as the obtained horizon speed of the stratum interface.
Further, the preset speed difference threshold interval range at least comprises 0.
Another aspect of the present invention provides a level extraction apparatus, comprising:
the longitudinal wave velocity acquisition module is used for acquiring the longitudinal wave velocity of a well depth sampling point in the well depth direction;
the stratum interface acquisition module is used for acquiring stratum interfaces according to the longitudinal wave speed or the current stratum thickness of the current sampling point, wherein the well depth sampling point comprises the current sampling point, and the current stratum thickness is the distance from the determined stratum interface closest to the current sampling point in all the determined stratum interfaces;
the horizon speed acquisition module is used for determining the obtained horizon speed of the stratum interface;
the acquisition module of the constant velocity coordinate point set is used for determining the constant velocity coordinate point set in a preset area of the stratum interface and taking the constant velocity coordinate point set as the extracted layer, the constant velocity coordinate point set comprises at least one constant velocity coordinate point, and the difference between the longitudinal wave velocity of the constant velocity coordinate point and the layer velocity is positioned in the range of a preset velocity difference threshold value interval.
Further, the formation interface acquiring module is specifically configured to:
acquiring the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point in the well depth direction, and determining the position of the current sampling point as the stratum interface when the difference between the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point and the longitudinal wave velocity of the current sampling point is greater than a first preset threshold value;
or when the thickness of the current stratum is larger than a second preset threshold value, determining the position of the current sampling point as the stratum interface.
Further, the constant velocity coordinate point set obtaining module is specifically configured to:
obtaining longitudinal wave velocity of a coordinate point to be selected in a preset area of the stratum interface, wherein the coordinate point to be selected is a point obtained by taking the obtained stratum interface or the coordinate point with the equal velocity as a datum point and taking a preset step length as a step length,
and when the difference between the longitudinal wave velocity of the coordinate point to be selected and the horizon velocity is within the range of the preset velocity difference threshold value interval, adding the coordinate point to be selected in the constant velocity coordinate point set, and taking the constant velocity coordinate point set as the extracted horizon.
Further, the horizon speed obtaining module is specifically configured to:
acquiring the obtained longitudinal wave speeds of all well depth sampling points in a layer where the stratum interface is located;
and taking the average value of the longitudinal wave speeds of all the well depth sampling points as the obtained horizon speed of the stratum interface.
Further, the preset speed difference threshold interval range at least comprises 0.
According to the horizon extraction method and device, the longitudinal wave velocity of a well depth sampling point in the well depth direction is obtained through the relation between the acoustic wave time difference and the longitudinal wave velocity; obtaining a stratum interface according to the longitudinal wave speed of the current sampling point or the current stratum thickness; then determining the obtained horizon speed of the stratum interface; and finally, determining an equal-speed coordinate point set in a preset region of the stratum interface, and taking the equal-speed coordinate point set as the extracted layer, wherein the equal-speed coordinate point set comprises at least one equal-speed coordinate point, and the difference between the longitudinal wave speed of the equal-speed coordinate point and the layer speed is within a preset speed difference threshold interval range. The method for extracting the horizon through the constant velocity value avoids the situations of inaccurate horizon pickup and multiple solutions when peaks or valleys are directly picked up from seismic data, and can acquire more stable depth domain horizon data.
Drawings
The invention will be described in more detail hereinafter on the basis of embodiments and with reference to the accompanying drawings. Wherein:
FIG. 1 is a schematic flow chart of a horizon extracting method according to an embodiment of the present invention;
fig. 2 is a schematic flowchart of a data access method according to a second embodiment of the present invention;
fig. 3 is a schematic diagram illustrating a blocking processing result of logging data in the data access method according to the second embodiment of the present invention;
FIG. 4 is a diagram illustrating an effect of matching a level line with a three-dimensional data volume in the data access method according to the second embodiment of the present invention;
fig. 5 is a schematic structural diagram of a data access device according to a third embodiment of the present invention.
In the drawings, like parts are provided with like reference numerals. The figures are not drawn to scale.
Detailed Description
The invention will be further explained with reference to the drawings.
The following detailed description will be provided for the embodiments of the present invention with reference to the accompanying drawings and examples, so that the practitioner of the present invention can fully understand how to apply the technical means to solve the technical problems, achieve the technical effects, and implement the present invention according to the implementation procedures. It should be noted that, as long as there is no conflict, the embodiments and the features of the embodiments of the present invention may be combined with each other, and the technical solutions formed are within the scope of the present invention.
Example one
FIG. 1 is a schematic flow chart of a horizon extracting method according to an embodiment of the present invention; as shown in fig. 1, the present embodiment provides a method for extracting a layer, which includes the following steps:
A plurality of well depth sampling points are preset in the well depth direction, and well logging data of each well depth sampling point are obtained. The acoustic wave time difference data in the logging data can be obtained through measurement, then the acoustic wave time difference data is converted into the longitudinal wave velocity by utilizing the relation between the acoustic wave time difference and the velocity, and the longitudinal wave velocity of the well depth sampling point can be obtained through the method.
102, obtaining stratum interfaces according to the longitudinal wave speed or the current stratum thickness of the current sampling point, wherein the well depth sampling point comprises the current sampling point, the current stratum thickness is the distance from the determined stratum interface closest to the current sampling point from all the determined stratum interfaces.
Specifically, each point in the well depth sampling points is sequentially regarded as the current sampling point to be analyzed. Each current sampling point is a position point in the well depth direction and corresponds to a longitudinal wave speed. Determining a stratum interface according to the longitudinal wave speed of the current sampling point, for example, acquiring the longitudinal wave speed of a well depth sampling point adjacent to the current sampling point, and when the current sampling point is a first point and a last point in the well depth direction, not setting the interface. When the current sampling points are the first point and the last point in the well depth sampling point in the well depth direction, two points adjacent to the current sampling point are provided, and if the longitudinal wave speed of at least one well depth sampling point adjacent to the current sampling point and the longitudinal wave speed of the current sampling point are greater than a first preset threshold value, the positions of the two well depth sampling points are the stratum interface.
And 103, determining the obtained horizon speed of the stratum interface.
The horizon velocity of the formation interfaces may be determined by:
acquiring the obtained longitudinal wave speeds of all well depth sampling points in a layer where the stratum interface is located;
and taking the average value of the longitudinal wave speeds of all the well depth sampling points as the obtained horizon speed of the stratum interface.
The formation in which the formation interface is located is first determined, for example, the layer in which the formation interface is located may be a portion between the formation interface (which may be considered as the current formation interface) and an adjacent formation interface, which is a formation interface located on the upper side of the current formation interface (which is closer to the wellhead location).
And averaging the longitudinal wave speeds of all the well depth sampling points in the layer where the stratum interface is located, and taking the average value as the layer speed of the stratum interface.
And 104, determining an equal-velocity coordinate point set in a preset region of the stratum interface, wherein the equal-velocity coordinate point set is used as the extracted layer, the equal-velocity coordinate point set comprises at least one equal-velocity coordinate point, and the difference between the longitudinal wave velocity of the equal-velocity coordinate point and the layer velocity is within a preset velocity difference threshold interval range.
The preset area can be selected according to the actual situation, and is not limited herein, and the preset area is a three-dimensional space area. The stratum interface is the coordinates of a certain well depth sampling point in the well depth direction, so the stratum interface can be taken as a datum point, and an equal-speed coordinate point is determined in a preset area of the datum point. For example, the longitudinal wave velocity of a certain position point in the preset area of the reference point is obtained, and if the difference between the longitudinal wave velocity and the horizon velocity is within the preset velocity difference threshold interval, the position point is the coordinate point of the constant velocity. The preset speed difference threshold interval may be set according to actual conditions, and is not limited herein, and preferably, the preset speed difference threshold interval is an interval with 0 as a symmetric point.
In the horizon extraction method provided by this embodiment, the longitudinal wave velocity of a well depth sampling point in the well depth direction is obtained through the relationship between the acoustic wave time difference and the longitudinal wave velocity; obtaining a stratum interface according to the longitudinal wave speed of the current sampling point or the current stratum thickness; then determining the obtained horizon speed of the stratum interface; and finally, determining an equal-speed coordinate point set in a preset region of the stratum interface, and taking the equal-speed coordinate point set as the extracted layer, wherein the equal-speed coordinate point set comprises at least one equal-speed coordinate point, and the difference between the longitudinal wave speed of the equal-speed coordinate point and the layer speed is within a preset speed difference threshold interval range. The method for extracting the horizon through the constant velocity value avoids the situations of inaccurate horizon pickup and multiple solutions when peaks or valleys are directly picked up from seismic data, and can acquire more stable depth domain horizon data.
Example two
This embodiment is a supplementary explanation based on the first embodiment.
FIG. 2 is a schematic flowchart of a horizon extracting method according to a second embodiment of the present invention; as shown in fig. 2, the present embodiment provides a method for extracting a layer, which includes the following steps:
202, acquiring the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point in the well depth direction, and determining the position of the current sampling point as the stratum interface when the difference between the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point and the longitudinal wave velocity of the current sampling point is greater than a first preset threshold value; or when the thickness of the current stratum is larger than a second preset threshold value, determining the position of the current sampling point as the stratum interface.
The first preset threshold and the second preset threshold may be set according to actual conditions, and are not limited herein. And carrying out logging data layering according to the first preset threshold and the second preset threshold. If Δ v > dv0Then Block 1, the current sample point is located at the formation interface, where Δ v | v ═ v2-v1|，v1,v2Corresponding longitudinal wave velocity v for two well depth sampling points adjacent to logging information1,v2One of them is the current sample point, dv0Is a first predetermined threshold, also called the minimum speed difference. If h > h0Then Block is 1, where h is the current formation thickness, h0The sampling point is a second preset threshold value, also called the maximum stratum thickness, and the position of the current sampling point is a stratum interface.
And step 203, determining the obtained horizon speed of the stratum interface.
Further, step 203 includes steps 2031 to 2032:
Specifically, the longitudinal wave velocities of all well depth sampling points in the layer where all the stratum interfaces are located are averaged, and the average value is the layer velocity of the stratum interfaces. And determining the horizon speed of the stratum interface as the standard of horizon picking.
In order to prevent the layer crossing phenomenon between layers, the step adopts a base point tracking technology to pick up the layers: using stratum boundary surface as base point (x)0,y0,z0) And performing constant velocity tracking in a preset area, namely: u. of1∈(x0±dx,y0±dy,z0D z) where u1Representing the point on the horizon obtained by the base point tracking, wherein the point is a coordinate point to be selected. By using the technology, the search space is reduced, and the calculation cost is reduced.
Meanwhile, a compatible speed value is introduced in the step to ensure the stability of the picked constant-speed (longitudinal wave speed value) horizon: v (u)1)∈va+ -. DELTA.v, wherein vaRepresenting the horizon speed corresponding to the horizon in which the current stratum interface is positioned, Deltav represents an allowable speed difference value, v (u)1) Represents u1Velocity of longitudinal waves at a point. Resulting constant velocity volume based depth domain horizon data is formed. Finally obtaining the constant velocity value vi(i ═ 1,2 … … n) is a standard three-dimensional space coordinate file, where n is the number of horizons.
In the horizon extraction method provided by this embodiment, the longitudinal wave velocity of a well depth sampling point in the well depth direction is obtained through the relationship between the acoustic wave time difference and the longitudinal wave velocity; acquiring the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point in the well depth direction, and determining the position of the current sampling point as the stratum interface when the difference between the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point and the longitudinal wave velocity of the current sampling point is greater than a first preset threshold value; or when the thickness of the current stratum is larger than a second preset threshold value, determining the position of the current sampling point as the stratum interface; determining the obtained horizon speed of the stratum interface; and in a preset area of the stratum interface, acquiring the longitudinal wave speed of a coordinate point to be selected, and when the difference between the longitudinal wave speed of the coordinate point to be selected and the horizon speed is within the range of a preset speed difference threshold interval, adding the coordinate point to be selected in the constant speed coordinate point set, and taking the constant speed coordinate point set as the extracted horizon. The method for extracting the horizon through the constant velocity value avoids the situations of inaccurate horizon pickup and multiple solutions when peaks or valleys are directly picked up from seismic data, and can acquire more stable depth domain horizon data.
The method is verified by the actual data of a certain area. Fig. 3 is a block diagram of the logging data obtained by using the minimum speed difference (i.e., the first preset threshold) and the maximum thickness (i.e., the second preset threshold), where the blue line (the line with large fluctuation in the diagram) represents the actual logging data, the red line (the straight line penetrating through the blue line in the diagram) is a curve obtained after fast transformation, and the start and end positions of each block are the horizon positions at the coordinates of the logging positions. Fig. 4 is a diagram showing the reflection of horizon data obtained by performing horizon tracking using a depth domain three-dimensional velocity volume on a vertical and horizontal measurement line. It can be seen from FIG. 4 that the horizon data fully conforms to the structural morphology of the original seismic volume. In the finally obtained depth domain horizon data map, the same horizon data are consistent in color. The method can accurately pick up the horizon data of the depth domain by utilizing the logging data and the velocity volume data of the three-dimensional depth domain, has obvious effect and has stability and practicability.
EXAMPLE III
The present embodiment is used for executing the method of the first embodiment.
Fig. 5 is a schematic structural diagram of a data access device according to a third embodiment of the present invention; as shown in fig. 5, the present embodiment provides a data access apparatus 300, which includes a longitudinal wave velocity obtaining module 301, a formation interface obtaining module 302, a horizon velocity obtaining module 303, and a constant velocity coordinate point set obtaining module 304.
The longitudinal wave velocity obtaining module 301 is configured to obtain a longitudinal wave velocity of a well depth sampling point in the well depth direction;
the stratum interface acquisition module 302 is configured to acquire a stratum interface according to the longitudinal wave speed or the current stratum thickness of the current sampling point, where the well depth sampling point includes the current sampling point, and the current stratum thickness is a distance from the determined stratum interface closest to the current sampling point, among all the determined stratum interfaces, to the current sampling point;
a horizon speed obtaining module 303, configured to determine an obtained horizon speed of the formation interface;
an equal-velocity coordinate point set obtaining module 304, configured to determine an equal-velocity coordinate point set in a preset region of the formation interface, and use the equal-velocity coordinate point set as the extracted layer, where the equal-velocity coordinate point set includes at least one equal-velocity coordinate point, and a difference between a longitudinal wave velocity of the equal-velocity coordinate point and the layer velocity is within a preset velocity difference threshold interval range.
Further, the formation interface obtaining module 302 is specifically configured to:
acquiring the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point in the well depth direction, and determining the position of the current sampling point as the stratum interface when the difference between the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point and the longitudinal wave velocity of the current sampling point is greater than a first preset threshold value;
or when the thickness of the current stratum is larger than a second preset threshold value, determining the position of the current sampling point as the stratum interface.
Further, the constant velocity coordinate point set obtaining module 304 is specifically configured to:
obtaining longitudinal wave velocity of a coordinate point to be selected in a preset area of the stratum interface, wherein the coordinate point to be selected is a point obtained by taking the obtained stratum interface or the coordinate point with the equal velocity as a datum point and taking a preset step length as a step length,
and when the difference between the longitudinal wave velocity of the coordinate point to be selected and the horizon velocity is within the range of the preset velocity difference threshold value interval, adding the coordinate point to be selected in the constant velocity coordinate point set, and taking the constant velocity coordinate point set as the extracted horizon.
Further, the horizon speed obtaining module 303 is specifically configured to:
acquiring the obtained longitudinal wave speeds of all well depth sampling points in a layer where the stratum interface is located;
and taking the average value of the longitudinal wave speeds of all the well depth sampling points as the obtained horizon speed of the stratum interface.
Further, the preset speed difference threshold interval range at least comprises 0.
The data access device in this embodiment can implement each process in the method embodiments of fig. 1 to fig. 2, and is not described here again to avoid repetition.
The data access device 300 provided in this embodiment obtains the longitudinal wave velocity of the well depth sampling point in the well depth direction according to the relationship between the acoustic time difference and the longitudinal wave velocity; obtaining a stratum interface according to the longitudinal wave speed of the current sampling point or the current stratum thickness; then determining the obtained horizon speed of the stratum interface; and finally, determining an equal-speed coordinate point set in a preset region of the stratum interface, and taking the equal-speed coordinate point set as the extracted layer, wherein the equal-speed coordinate point set comprises at least one equal-speed coordinate point, and the difference between the longitudinal wave speed of the equal-speed coordinate point and the layer speed is within a preset speed difference threshold interval range. The method for extracting the horizon through the constant velocity value avoids the situations of inaccurate horizon pickup and multiple solutions when peaks or valleys are directly picked up from seismic data, and can acquire more stable depth domain horizon data.
It is to be understood that the disclosed embodiments of the invention are not limited to the particular structures, process steps, or materials disclosed herein but are extended to equivalents thereof as would be understood by those ordinarily skilled in the relevant arts. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting.
Reference in the specification to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the invention. Thus, the appearances of the phrase "one embodiment" or "an embodiment" in various places throughout this specification are not necessarily all referring to the same embodiment.
Furthermore, the described features, structures, or characteristics may be combined in any other suitable manner in one or more embodiments. In the previous descriptions, some specific details are provided, such as materials, etc., to provide a thorough understanding of embodiments of the invention. One skilled in the relevant art will recognize, however, that the invention may be practiced without one or more of the specific details, or with other methods, components, materials, and so forth. In other instances, well-known structures, materials, or operations are not shown or described in detail to avoid obscuring aspects of the invention.
While the above examples are illustrative of the principles of the present invention in one or more applications, it will be apparent to those of ordinary skill in the art that various changes in form, usage and details of implementation can be made without departing from the principles and concepts of the invention. Accordingly, the invention is defined by the appended claims.
Claims (6)
1. A method of extracting a layer, comprising:
acquiring the longitudinal wave speed of a well depth sampling point in the well depth direction;
obtaining stratum interfaces according to the longitudinal wave speed or the current stratum thickness of the current sampling point, wherein the well depth sampling point comprises the current sampling point, the current stratum thickness is the distance from the determined stratum interface closest to the current sampling point in all the determined stratum interfaces;
determining the obtained horizon speed of the stratum interface;
determining an equal-velocity coordinate point set in a preset region of the stratum interface, and taking the equal-velocity coordinate point set as the extracted layer, wherein the equal-velocity coordinate point set comprises at least one equal-velocity coordinate point, and the difference between the longitudinal wave velocity of the equal-velocity coordinate point and the layer velocity is within a preset velocity difference threshold interval range;
the step of obtaining the stratum interface according to the longitudinal wave speed or the current stratum thickness of the current sampling point comprises the following steps: acquiring the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point in the well depth direction, and determining the position of the current sampling point as the stratum interface when the difference between the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point and the longitudinal wave velocity of the current sampling point is greater than a first preset threshold value; or when the thickness of the current stratum is larger than a second preset threshold value, determining the position of the current sampling point as the stratum interface;
the step of determining a constant velocity coordinate point set in a preset region of the stratum interface and taking the constant velocity coordinate point set as the extracted horizon includes: in the preset area of the stratum interface, longitudinal wave speed of a coordinate point to be selected is obtained, wherein the coordinate point to be selected is obtained by taking the obtained stratum interface or the constant speed coordinate point as a datum point and taking a preset step length as a step length, when the difference between the longitudinal wave speed of the coordinate point to be selected and the horizon speed is within the range of a preset speed difference threshold value interval, the coordinate point to be selected is added in the constant speed coordinate point set, and the constant speed coordinate point set is taken as the extracted horizon.
2. The method of horizon extraction of claim 1 wherein the step of determining the obtained horizon velocity of the stratigraphic interface comprises:
acquiring the obtained longitudinal wave speeds of all well depth sampling points in a layer where the stratum interface is located;
and taking the average value of the longitudinal wave speeds of all the well depth sampling points as the obtained horizon speed of the stratum interface.
3. The horizon extraction method according to claim 1 or 2, characterized in that the preset speed difference threshold interval range comprises at least 0.
4. A level extraction device, comprising:
the longitudinal wave velocity acquisition module is used for acquiring the longitudinal wave velocity of a well depth sampling point in the well depth direction;
the stratum interface acquisition module is used for acquiring stratum interfaces according to the longitudinal wave speed or the current stratum thickness of the current sampling point, wherein the well depth sampling point comprises the current sampling point, and the current stratum thickness is the distance from the determined stratum interface closest to the current sampling point in all the determined stratum interfaces;
the horizon speed acquisition module is used for determining the obtained horizon speed of the stratum interface;
the constant velocity coordinate point set acquisition module is used for determining a constant velocity coordinate point set in a preset region of the stratum interface and taking the constant velocity coordinate point set as the extracted layer, wherein the constant velocity coordinate point set comprises at least one constant velocity coordinate point, and the difference between the longitudinal wave velocity of the constant velocity coordinate point and the layer velocity is within a preset velocity difference threshold interval range;
the formation interface acquisition module is specifically configured to: acquiring the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point in the well depth direction, and determining the position of the current sampling point as the stratum interface when the difference between the longitudinal wave velocity of the well depth sampling point adjacent to the current sampling point and the longitudinal wave velocity of the current sampling point is greater than a first preset threshold value; or when the thickness of the current stratum is larger than a second preset threshold value, determining the position of the current sampling point as the stratum interface;
the constant velocity coordinate point set acquisition module is specifically configured to: in the preset area of the stratum interface, longitudinal wave speed of a coordinate point to be selected is obtained, wherein the coordinate point to be selected is obtained by taking the obtained stratum interface or the constant speed coordinate point as a datum point and taking a preset step length as a step length, when the difference between the longitudinal wave speed of the coordinate point to be selected and the horizon speed is within the range of a preset speed difference threshold value interval, the coordinate point to be selected is added in the constant speed coordinate point set, and the constant speed coordinate point set is taken as the extracted horizon.
5. The horizon extraction apparatus of claim 4, wherein the horizon speed acquisition module is specifically configured to:
acquiring the obtained longitudinal wave speeds of all well depth sampling points in a layer where the stratum interface is located;
and taking the average value of the longitudinal wave speeds of all the well depth sampling points as the obtained horizon speed of the stratum interface.
6. The horizon extraction apparatus according to claim 4 or 5, characterized in that the preset speed difference threshold interval range comprises at least 0.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811013943.2A CN110873892B (en) | 2018-08-31 | 2018-08-31 | Horizon extraction method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811013943.2A CN110873892B (en) | 2018-08-31 | 2018-08-31 | Horizon extraction method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110873892A CN110873892A (en) | 2020-03-10 |
| CN110873892B true CN110873892B (en) | 2021-10-22 |
Family
ID=69715347
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811013943.2A Active CN110873892B (en) | 2018-08-31 | 2018-08-31 | Horizon extraction method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110873892B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113687423B (en) * | 2021-08-04 | 2024-03-01 | 潜能恒信能源技术股份有限公司 | Seismic data interpretation method and system based on interface extraction |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5648937A (en) * | 1995-01-18 | 1997-07-15 | Atlantic Richfield Company | Method and apparatus for correlating geological structure horizons from velocity data to well observations |
| CN101506686A (en) * | 2006-06-21 | 2009-08-12 | 特拉斯帕克地球科学公司 | Interpretation of geologic depositional systems |
| CN103149588A (en) * | 2013-02-20 | 2013-06-12 | 中国石油天然气股份有限公司 | Method and system for calculating VTI (Velocity Time Integral) anisotropic parameter by utilizing well seismic calibration |
| CN105301638A (en) * | 2015-10-12 | 2016-02-03 | 中国石油天然气集团公司 | Method and device for extracting bottom interface of weathered layer |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7917295B2 (en) * | 2008-04-30 | 2011-03-29 | Westerngeco L.L.C. | Modeling and filtering coherent noise in seismic surveying |
-
2018
- 2018-08-31 CN CN201811013943.2A patent/CN110873892B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5648937A (en) * | 1995-01-18 | 1997-07-15 | Atlantic Richfield Company | Method and apparatus for correlating geological structure horizons from velocity data to well observations |
| CN101506686A (en) * | 2006-06-21 | 2009-08-12 | 特拉斯帕克地球科学公司 | Interpretation of geologic depositional systems |
| CN103149588A (en) * | 2013-02-20 | 2013-06-12 | 中国石油天然气股份有限公司 | Method and system for calculating VTI (Velocity Time Integral) anisotropic parameter by utilizing well seismic calibration |
| CN105301638A (en) * | 2015-10-12 | 2016-02-03 | 中国石油天然气集团公司 | Method and device for extracting bottom interface of weathered layer |
Non-Patent Citations (4)
| Title |
|---|
| 3D acoustic isotropic full waveform inversion of land data application;Guanghui Hu et al.;《SEG 2017 Workshop: Full-waveform Inversion and Beyond》;20171231;第1-4页 * |
| 井间地震资料层位标定方法研究;杨国权等;《石油物探》;20050531;第44卷(第3期);第217页 * |
| 构造约束全波形反演及其海上资料应用;胡光辉等;《石油物探》;20180731;第57卷(第4期);第592-596页 * |
| 特征点法用于合成声测井解释;边西燕等;《石油物探》;19970630;第36卷(第2期);第109-115页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110873892A (en) | 2020-03-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110208859B (en) | Intelligent calculation method for quantitative parameters of oil-based mud well fractures based on ultrasonic imaging | |
| CN110308444B (en) | Road horizon intelligent identification and interference source elimination method | |
| CN105093307B (en) | The oblique true formation thickness acquiring method of Lower Paleozoic strata | |
| CN104808248A (en) | Attitude of stratum detecting method based on electric imaging logging full hole image | |
| CN111551927B (en) | Underground pipeline diameter measuring method based on three-dimensional ground penetrating radar | |
| CN108375785B (en) | Method and device for correcting position of crack belt | |
| CN104375172A (en) | Volcanic underlayer structural configuration correct imaging method | |
| CN115373024A (en) | Method and device for inverting passive land edge crustal structure based on stratum recording settlement | |
| CN110873892B (en) | Horizon extraction method and device | |
| Di et al. | 3D structural-orientation vector guided autotracking for weak seismic reflections: A new tool for shale reservoir visualization and interpretation | |
| CN105093318B (en) | A kind of adaptive wave equation wave field extrapolation static correcting method | |
| CN110376661A (en) | The control methods of single river channel sand and device between a kind of missile silo | |
| CN104297753B (en) | A kind of pathfinder image inverting wind direction of ocean surface method based on adaptive shrinkage operator | |
| CN106443674B (en) | A kind of Ground Penetrating Radar Velocity Estimation based on diffraction and imaging and minimum entropy technique | |
| CN102338885A (en) | First-arrival-time automatic picking method of tricomponent VSP (Vertical Seismic Profile) data | |
| CN112684498A (en) | Reservoir fracture prediction method and system based on wide-azimuth seismic data | |
| CN110261905A (en) | Complex value based on pitch angle control is concerned with microfault recognition methods | |
| CN109490962A (en) | A kind of removing method of shallow seismic profile Non-zero Offset | |
| CN105866834A (en) | Single branch refraction data interpretation method based on automatic identification refraction wave starting point and system thereof | |
| CN109116416A (en) | Inclination ocean bottom seismometer secondary localization method based on three dimensional constitution | |
| CN109190318B (en) | Inclined shaft trajectory calculation method based on time measurement | |
| CN109143364B (en) | First arrival time picking method | |
| CN113552623A (en) | Method and device for determining boundary of sediment sand body of meandering river | |
| CN111474580A (en) | Azimuth angle gather extraction method and system based on offset vector piece | |
| CN107356970B (en) | High-precision well seismic data matching method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |