CN112305596A - Automatic horizon tracking method and device under fault control - Google Patents
Automatic horizon tracking method and device under fault control Download PDFInfo
- Publication number
- CN112305596A CN112305596A CN201910681480.5A CN201910681480A CN112305596A CN 112305596 A CN112305596 A CN 112305596A CN 201910681480 A CN201910681480 A CN 201910681480A CN 112305596 A CN112305596 A CN 112305596A
- Authority
- CN
- China
- Prior art keywords
- horizon
- tracking
- fault
- seed
- seed 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000000903 blocking effect Effects 0.000 claims abstract description 23
- 238000007781 pre-processing Methods 0.000 claims abstract description 21
- 238000010008 shearing Methods 0.000 claims abstract description 15
- 230000002159 abnormal effect Effects 0.000 claims abstract description 14
- 238000004590 computer program Methods 0.000 claims description 18
- 238000005520 cutting process Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 8
- 230000011218 segmentation Effects 0.000 claims description 5
- 238000000638 solvent extraction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 abstract description 12
- 238000010586 diagram Methods 0.000 description 22
- 238000012545 processing Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 4
- 230000000717 retained effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000013213 extrapolation Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect 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
-
- 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
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)
- Image Analysis (AREA)
Abstract
The invention provides a method and a device for automatically tracking a horizon under fault control, wherein the method for automatically tracking the horizon under fault control comprises the following steps: preprocessing the horizon seed point data to remove abnormal values; automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets; respectively carrying out automatic horizon tracking on each horizon seed point set to obtain horizon tracking results corresponding to the blocks; and shearing the horizon tracking result by using the fault blocks to obtain a horizon tracking result corresponding to each fault block. The invention can adapt to various complex construction conditions, greatly improve the efficiency and meet the requirement of high-efficiency and high-precision automatic horizon tracking under the complex construction conditions.
Description
Technical Field
The application relates to the technical field of computer data processing, in particular to a method and a device for automatically tracking a horizon under fault control.
Background
Horizon interpretation is the basic work of seismic data interpretation, and all subsequent work depends on horizon interpretation data, and is one of the most time-consuming works in structure interpretation, and especially in fault development areas, the workload of horizon interpretation is huge.
The current main methods for horizon interpretation include automatic, semi-automatic and manual methods. When the fracture is complicated, it is difficult to simultaneously ensure horizon interpretation accuracy and efficiency. The reasons for influencing the efficiency and accuracy of horizon interpretation are mainly reflected in the following two aspects:
firstly, geological conditions and the quality of seismic data are adopted, and when the quality of the seismic data of a work area is poor and faults are complex, the tracking result of the conventional automatic tracking algorithm is poor, so that the precision is reduced, and the interpretation precision cannot be met.
Secondly, when the number of fault positions is large, the geological structure is complex, and semi-automatic interpretation and manual interpretation are adopted, a large amount of manual intervention is needed, the interpretation efficiency is greatly reduced, and the efficiency cannot be met.
Disclosure of Invention
Aiming at the problems in the prior art, the application provides a method and a device for automatically tracking the horizon under fault control, which can adapt to various complex construction conditions, greatly improve the efficiency and meet the requirement of high-efficiency and high-precision automatic tracking of the horizon under the complex construction conditions.
In order to solve the technical problem, the application provides the following technical scheme:
in a first aspect, the present application provides a method for automatically tracking a horizon under fault control, including:
preprocessing the horizon seed point data to remove abnormal values;
automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets;
respectively carrying out automatic horizon tracking on each horizon seed point set to obtain horizon tracking results corresponding to the blocks;
and shearing the horizon tracking result by using the fault blocks to obtain a horizon tracking result corresponding to each fault block.
Further, the preprocessing the horizon seed point data to remove abnormal values includes:
dividing each horizon seed point into a plurality of sections according to the position of the invalid value of the horizon seed point;
intersecting each of the multiple sections with a fault respectively, and cutting each section into two sections from the intersection point;
and removing the section with less data in the two sections, or controlling whether the section with less data is reserved or not through parameters.
Further, the automatically blocking the preprocessed horizon seed points under the condition of fault block constraint to obtain a plurality of horizon seed point sets, including:
interpolating and extrapolating the preprocessed data of the horizon seed points to obtain continuous horizon line segments;
intersecting the horizon line segment with a fault to obtain a plurality of intersection points;
dividing the horizon line segment into a plurality of regions according to the intersection point position;
and dividing the preprocessed horizon seed points into a plurality of sets according to the region.
Further, the cutting the horizon tracking result by using the fault blocks to obtain a horizon tracking result corresponding to each fault block includes:
calculating the intersection point of the horizon and the fault in the horizon tracking result;
dividing the horizon in the horizon tracking result into a plurality of intervals by using intersection information;
and obtaining the position of the fault block by using the horizon seed point of the horizon in the horizon tracking result, and reserving the interval data of the position.
In a second aspect, the present application provides an automatic horizon tracking apparatus under fault control, including:
the preprocessing unit is used for preprocessing the layer seed point data and removing abnormal values;
the blocking unit is used for automatically blocking the preprocessed horizon seed points under the condition of fault block constraint to obtain a plurality of horizon seed point sets;
the horizon tracking unit is used for respectively carrying out horizon automatic tracking on each horizon seed point set to obtain a horizon tracking result corresponding to each block;
and the shearing unit is used for shearing the horizon tracking result by using the fault blocks to obtain the horizon tracking result corresponding to each fault block.
Further, the preprocessing unit includes:
the segmentation module is used for dividing each horizon seed point into a plurality of segments according to the position of the invalid value of the horizon seed point;
the cutting module is used for respectively intersecting each section of the multiple sections with a fault and cutting each section into two sections from the intersection point;
and the filtering module is used for removing the section with less data in the two sections or controlling whether to keep the section with less data or not through parameters.
Further, the blocking unit includes:
the interpolation module is used for interpolating and extrapolating the preprocessed data of the horizon seed points to obtain continuous horizon line segments;
the intersection module is used for intersecting the horizon line segment with the fault to obtain a plurality of intersection points;
the partitioning module is used for dividing the horizon line segment into a plurality of areas according to the intersection point position;
and the diversity combining module is used for dividing the preprocessed horizon seed points into a plurality of sets according to the region.
In a third aspect, the present application provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of the automatic horizon tracking under fault control when executing the program.
In a fourth aspect, the present application provides a computer-readable storage medium having stored thereon a computer program which, when being executed by a processor, carries out the steps of the automatic horizon tracking method under fault control as described.
According to the above technical solution, the present application provides a method and an apparatus for automatically tracking a horizon under fault control, where the method for automatically tracking a horizon under fault control includes: preprocessing the horizon seed point data to remove abnormal values; automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets; respectively carrying out automatic horizon tracking on each horizon seed point set to obtain horizon tracking results corresponding to the blocks; and shearing the horizon tracking result by using the fault blocks to obtain the horizon tracking result corresponding to each fault block, so that the method can adapt to various complex construction conditions, can greatly improve the efficiency, and meets the requirement of high-efficiency and high-precision horizon automatic tracking under the complex construction conditions.
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.
FIG. 1 is a flow chart of automatic horizon tracking under fault control in an embodiment of the present application;
FIG. 2 is a flowchart illustrating an example of preprocessing the horizon seed point data to remove abnormal values;
FIG. 3 is a flowchart illustrating automatic blocking of preprocessed horizon seed points under fault block constraint conditions to obtain a plurality of horizon seed point sets according to an embodiment of the present application;
FIG. 4 is a schematic structural diagram of an apparatus for automatically tracking horizons under fault control according to an embodiment of the present application;
FIG. 5 is a schematic structural diagram of the pretreatment unit 1 in the embodiment of the present application;
fig. 6 is a schematic structural diagram of a partitioning unit 2 in the embodiment of the present application;
fig. 7 is a flowchart of S104 in the embodiment of the present application;
FIG. 8 shows a diagram of an originally interpreted horizon seed point and fault;
FIG. 9 is a schematic diagram of a layer position seed point pre-processed with an extra segment beyond a fault removed;
FIG. 10 is a schematic illustration of interpolation and extrapolation using pre-processed horizon seed points such that horizons fill the entire profile;
FIG. 11 is a schematic diagram of the intersection of the fully inserted horizon with the fault to obtain the precise breakpoint position;
FIG. 12 is a schematic diagram of a plurality of zonal site sets formed by breaking and blocking the source layer sites using the intersection points with the faults;
FIG. 13 is a schematic diagram illustrating the full-profile horizon auto-tracking using the horizon point set in the first segment as seed points to obtain the tracking result of the full-profile;
FIG. 14 is a schematic diagram of the intersection of the tracked horizon result with a fault, dividing the result into multiple segments;
FIG. 15 is a diagram illustrating the full-profile horizon auto-tracking using the horizon point set in the second segment as seed points to obtain the tracking result of the full-profile;
FIG. 16 is a schematic diagram of intersection of the horizon result tracked by the horizon seed points in the second block with the fault, dividing the result into multiple segments, and determining the valid segment;
FIG. 17 is a diagram illustrating the final effect of two horizon tracking by looping all fault blocks for tracking;
fig. 18 is a block diagram of a cutting unit according to an embodiment of the present invention;
fig. 19 is a schematic structural diagram of an electronic device in an embodiment of the invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Aiming at the problems that the accuracy of an automatic tracking algorithm is low and the efficiency is low when manual interpretation is adopted in the existing horizon interpretation method under the condition that the quality of seismic data is poor, the horizon automatic tracking method under fault control is provided, horizon seed point data is preprocessed, blocks are divided according to fault blocks, automatic tracking is respectively carried out, and then tracking results are obtained through shearing.
Fig. 1 is a flowchart of an automatic horizon tracking method under fault control according to an embodiment of the present application, and as shown in fig. 1, the automatic horizon tracking method under fault control includes:
s101: preprocessing the horizon seed point data to remove abnormal values;
for the reason of artificial explanation, part of the seed point data is a continuously distributed horizon section which possibly passes through the fault, the continuously distributed horizon section can be intersected with the fault, the excessive part is removed or split into two sections according to the geological condition, and the horizon automatic tracking is established later.
S102: and automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets.
S103: respectively carrying out automatic horizon tracking on each horizon seed point set to obtain horizon tracking results corresponding to the blocks;
specifically, each seed point set is used as input, and a general horizon automatic tracking algorithm is used for tracking the whole profile to obtain a horizon tracking result covering the whole profile.
S104: shearing the horizon tracking result by using the fault blocks to obtain a horizon tracking result corresponding to each fault block;
because a plurality of horizon tracking results are obtained in the last step, seed points of the horizon tracking results come from different fault blocks, the tracking result is correct in the fault block, and the tracking result is not necessarily correct outside the fault block, the excessive data outside the fault block needs to be cut off, and a correct result is left; specifically, a layer is selected, the intersection point of the layer and a fault is obtained, the layer is divided into multiple sections by using intersection point information, the position of a fault block is obtained by using a seed point for tracking the layer, data of the layer is reserved, and the rest of the data are not reserved.
As can be seen from the flow shown in FIG. 1, in the present application, the horizon seed point data is preprocessed to remove abnormal values; then, automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets; then, automatically tracking the horizon of each horizon seed point set to obtain a horizon tracking result corresponding to each block; and finally, shearing the horizon tracking result by using the fault blocks to obtain the horizon tracking result corresponding to each fault block. The method can adapt to various complex construction conditions, greatly improve the efficiency and meet the requirement of high-efficiency and high-precision automatic horizon tracking under the complex construction conditions.
In seismic data interpretation, manual interpretation is mostly adopted under the condition of complex geological conditions or poor seismic data, so that the quality of interpreted data is not high easily, for example, an interpreted horizon exceeds a fault, and in subsequent tracking, the tracking quality near the fault is poor easily, so that the nonstandard interpreted horizon needs to be automatically processed to meet the requirement of subsequent automatic tracking. In one embodiment, as shown in fig. 2, the step of preprocessing the horizon seed point data and removing the outlier includes:
s201: and dividing each interpretation horizon seed point into a plurality of sections according to the position of the invalid value of the horizon seed point.
For human interpretation reasons, some positions have interpretation data, and some positions have no interpretation data, so that an interpretation horizon can be divided into multiple segments (segments) according to the positions of invalid values of the horizons.
S202: intersecting each section with a fault, and cutting the section into two sections from the intersection point.
S203: removing the section with less data in the two sections, or controlling whether to reserve the section with less data or not through parameters; for example, if a segment with less data is larger than a certain length (which may be set according to the situation), the segment is retained and may be regarded as a valid segment (segment).
In an embodiment, as shown in fig. 3, the step of automatically blocking the preprocessed horizon seed points under the condition of the fault block constraint to obtain a plurality of horizon seed point sets includes:
s301: interpolating and extrapolating the data of the preprocessed horizon seed points to obtain continuous horizon line segments;
because the horizon seed points are discrete, only partial positions on the whole section have effective values, the horizon seed points are inserted into the whole section to prepare for subsequent segmentation block processing; specifically, the preprocessed horizon data is used to merge all segments into a new segment for full-region interpolation.
S302: intersecting the horizon line segment with a fault to obtain a plurality of intersection points; specifically, a common line-to-line intersection algorithm may be used to find the intersection points, which, after finding the intersection points, may divide the formation into multiple sections or zones.
S303: dividing the horizon line segment into a plurality of regions according to the intersection point position; specifically, the obtained intersections are sorted according to the spatial positions of the cross sections, and a plurality of line segments are assigned to the respective regions (blocks).
S304: dividing the preprocessed horizon seed points into a plurality of sets according to the region;
specifically, by calculating the raw point and separate region range relationships, the raw horizon seed points may be divided into a plurality of horizon seed point sets.
After the division into a plurality of sets, the preprocessed horizon seed points can be distributed into each fault block (one fault block is arranged between two faults), and the automatic blocking function of the horizon seed points is completed.
S103 mainly uses the seed point data of the partitioned layer to perform automatic tracking, which in an embodiment may include the following two steps:
step 1: after the preprocessing of S101, the horizon seed points are segmented into a plurality of sets by faults, and then one horizon seed point set is used as input to automatically track the horizon of the whole work area.
Step 2: and circularly tracking all the horizon seed point sets to finally form a plurality of horizon tracking results.
Since the S103 obtains a plurality of horizon tracking results, the seed points of the horizon tracking results come from different fault blocks, the tracking result is correct in the fault block, and the tracking result is not necessarily correct outside the fault block, so that the excess data outside the fault block needs to be cut off, and a correct result is left.
In one embodiment, as shown in fig. 7, S104 includes the following steps:
s701: calculating the intersection point of the horizon and the fault in the horizon tracking result;
s702: dividing the horizon in the horizon tracking result into a plurality of intervals by using intersection information;
s703: and obtaining the position of the fault block by using the horizon seed point of the horizon in the horizon tracking result, reserving the interval data of the position, and reserving the rest parts.
In the flow shown in fig. 7, one horizon is selected at a time to obtain corresponding interval data, and S701 to S703 are repeated until all horizons are processed.
The invention can not only adapt to various complex construction conditions, but also greatly improve the efficiency and meet the requirement of high-efficiency and high-precision automatic horizon tracking under the complex construction conditions.
The present invention is described in detail with reference to the specific embodiment, as shown in fig. 8, it can be seen from fig. 8 that the explained horizon seed data has a position beyond a part of the fault, has no intersection, has a gap, and has a part of the horizon in a fault block. After the pretreatment of the horizon seed points, the redundant sections exceeding the fault part are removed as shown in FIG. 9. As shown in FIG. 10, the pre-processed horizon seed points are used for interpolation and extrapolation so that the horizon fills the entire profile. As shown in fig. 11, the intersection of the full-inserted horizon and the fault is used to obtain the accurate breakpoint position. As shown in fig. 12, the original layer points are subjected to segmentation processing using the intersection points with the fault, and the layer points are allocated to each segment, that is, a plurality of regional layer point sets are formed. As shown in fig. 13, the horizon point set in the first segment is used as a seed point to perform full-profile horizon automatic tracking, so as to obtain a tracking result of the full-profile. As shown in fig. 14, the tracked horizon result is intersected with the fault, the result is divided into multiple segments, which segment is the horizon tracking result in the fault is determined according to the position of the original input horizon seed point set, and the result is retained, and the rest is not retained. As shown in fig. 15, the horizon point set in the second segment is used as a seed point to perform full-profile horizon automatic tracking, so as to obtain a tracking result of the full-profile. As shown in fig. 16, the horizon result tracked by using the horizon seed points in the second fault block is intersected with the fault, the result is divided into a plurality of segments, and the effective segment is determined. As shown in fig. 17, all fault blocks are cycled to trace, and the final effect of two horizon tracing is achieved.
Taking an example that the automatic horizon tracking method under fault control is applied to a certain work area in the east of China, the fault of the work area belongs to a positive fault, the horizon seed points and interpreted fault data are directly utilized without any manual participation, and the automatic horizon tracking is formed efficiently and accurately by utilizing the technology, so that a satisfactory result is obtained.
Based on the same inventive concept, the embodiment of the present application further provides an automatic horizon tracking apparatus under fault control, which can be used to implement the method described in the foregoing embodiment, as described in the following embodiment. Because the principle of the automatic horizon tracking device under fault control for solving the problems is similar to the automatic horizon tracking method under fault control, the implementation of the automatic horizon tracking device under fault control can refer to the implementation of the automatic horizon tracking method under fault control, and repeated parts are not described again. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. While the system described in the embodiments below is preferably implemented in software, implementations in hardware, or a combination of software and hardware are also possible and contemplated.
As shown in fig. 4, the present application provides an embodiment of an automatic horizon tracking apparatus under fault control, and as shown in fig. 4, the automatic horizon tracking apparatus includes:
the preprocessing unit 1 is used for preprocessing the horizon seed point data and removing abnormal values;
the blocking unit 2 is used for automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets;
a horizon tracking unit 3, configured to perform horizon automatic tracking on each horizon seed point set, respectively, to obtain a horizon tracking result corresponding to each block;
and the shearing unit 4 is used for shearing the horizon tracking result by using the fault blocks to obtain the horizon tracking result corresponding to each fault block.
In one embodiment, as shown in fig. 5, the preprocessing unit 1 includes:
the segmentation module 11 is configured to divide each interpretation horizon seed point into multiple segments according to the position of the invalid value of the horizon seed point;
a cutting module 12, configured to intersect each of the segments with a fault, and cut the segment into two segments from the intersection;
and the filtering module 13 is used for removing the section with less data in the two sections, or controlling whether to keep the section with less data or not through parameters.
In an embodiment, as shown in fig. 6, the blocking unit 2 includes:
the interpolation module 21 is configured to interpolate and extrapolate the data of the preprocessed horizon seed points to obtain a continuous horizon line segment;
the intersection module 22 is configured to use the horizon line segment to intersect with a fault to obtain a plurality of intersection points;
a partitioning module 23, configured to divide the horizon line segment into a plurality of regions according to the intersection positions;
and the diversity combining module 24 is configured to divide the preprocessed horizon seed points into a plurality of sets according to the region.
In one embodiment, as shown in fig. 18, the cutting unit 4 includes:
an intersection point calculating module 1801, configured to calculate an intersection point between the horizon and the fault in the horizon tracking result;
an interval division module 1802, configured to divide a horizon in the horizon tracking result into multiple intervals by using intersection information;
a position determining module 1803, configured to obtain a position of the fault block by using the horizon seed point for tracking the horizon in the horizon tracking result, and retain interval data of the position.
The invention can not only adapt to various complex construction conditions, but also greatly improve the efficiency and meet the requirement of high-efficiency and high-precision automatic horizon tracking under the complex construction conditions.
An embodiment of the present application further provides a specific implementation manner of an electronic device capable of implementing all steps in the automatic horizon tracking method under fault control in the foregoing embodiment, and referring to fig. 19, the electronic device specifically includes the following contents:
a processor (processor)1201, a memory (memory)1202, a communication Interface 1203, and a bus 1204;
the processor 1201, the memory 1202 and the communication interface 1203 complete communication with each other through the bus 1204; the communication interface 1203 is configured to implement information transmission between related devices, such as a server-side device, a measurement device, and a client device.
The processor 1201 is configured to call a computer program in the memory 1202, and the processor executes the computer program to implement all the steps in the automatic horizon tracking method under fault control in the above embodiments, for example, to implement the following steps when the processor executes the computer program:
step 1: preprocessing the horizon seed point data to remove abnormal values;
for the reason of artificial explanation, part of the seed point data is a continuously distributed horizon section which possibly passes through the fault, the continuously distributed horizon section can be intersected with the fault, the excessive part is removed or split into two sections according to the geological condition, and the horizon automatic tracking is established later.
Step 2: and automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets.
And step 3: respectively carrying out automatic horizon tracking on each horizon seed point set to obtain horizon tracking results corresponding to the blocks;
specifically, each seed point set is used as input, and a general horizon automatic tracking algorithm is used for tracking the whole profile to obtain a horizon tracking result covering the whole profile.
And 4, step 4: and shearing the horizon tracking result by using the fault blocks to obtain a horizon tracking result corresponding to each fault block.
From the above description, it can be seen that, with the electronic device in the embodiment of the present application, not only various complex structural conditions can be adapted, but also the efficiency can be greatly improved, and the requirement of high-efficiency and high-precision automatic horizon tracking under the complex structural conditions can be met.
Embodiments of the present application also provide a computer-readable storage medium capable of implementing all steps in the automatic horizon tracking under fault control method in the foregoing embodiments, where the computer-readable storage medium stores thereon a computer program, and when the computer program is executed by a processor, the computer program implements all steps of the automatic horizon tracking under fault control method in the foregoing embodiments, for example, when the processor executes the computer program, the processor implements the following steps:
step 1: preprocessing the horizon seed point data to remove abnormal values;
for the reason of artificial explanation, part of the seed point data is a continuously distributed horizon section which possibly passes through the fault, the continuously distributed horizon section can be intersected with the fault, the excessive part is removed or split into two sections according to the geological condition, and the horizon automatic tracking is established later.
Step 2: and automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets.
And step 3: respectively carrying out automatic horizon tracking on each horizon seed point set to obtain horizon tracking results corresponding to the blocks;
specifically, each seed point set is used as input, and a general horizon automatic tracking algorithm is used for tracking the whole profile to obtain a horizon tracking result covering the whole profile.
And 4, step 4: and shearing the horizon tracking result by using the fault blocks to obtain a horizon tracking result corresponding to each fault block.
From the above description, it can be seen that the computer-readable storage medium in the embodiment of the present application can not only adapt to various complex structural conditions, but also greatly improve efficiency, and meet the requirement of high-efficiency and high-precision automatic horizon tracking under the complex structural conditions.
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 principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
Claims (10)
1. A method for automatically tracking horizons under fault control is characterized by comprising the following steps:
preprocessing the horizon seed point data to remove abnormal values;
automatically blocking the preprocessed layer seed points under the condition of broken block constraint to obtain a plurality of layer seed point sets;
respectively carrying out automatic horizon tracking on each horizon seed point set to obtain horizon tracking results corresponding to the blocks;
and shearing the horizon tracking result by using the fault blocks to obtain a horizon tracking result corresponding to each fault block.
2. The method for automatically tracking horizons under fault control according to claim 1, wherein the preprocessing of the horizon seed point data to remove abnormal values comprises:
dividing each horizon seed point into a plurality of sections according to the position of the invalid value of the horizon seed point;
intersecting each of the multiple sections with a fault respectively, and cutting each section into two sections from the intersection point;
and removing the section with less data in the two sections, or controlling whether the section with less data is reserved or not through parameters.
3. The method for automatically tracking horizons under fault control according to claim 1, wherein the automatically blocking the preprocessed horizon seed points under fault block constraint conditions to obtain a plurality of horizon seed point sets comprises:
interpolating and extrapolating the preprocessed data of the horizon seed points to obtain continuous horizon line segments;
intersecting the horizon line segment with a fault to obtain a plurality of intersection points;
dividing the horizon line segment into a plurality of regions according to the intersection point position;
and dividing the preprocessed horizon seed points into a plurality of sets according to the region.
4. The method for automatically tracking the horizon under the fault control according to claim 1, wherein the cutting the horizon tracking result by using the fault blocks to obtain the horizon tracking result corresponding to each fault block comprises:
calculating the intersection point of the horizon and the fault in the horizon tracking result;
dividing the horizon in the horizon tracking result into a plurality of intervals by using intersection information;
and obtaining the position of the fault block by using the horizon seed point of the horizon in the horizon tracking result, and reserving the interval data of the position.
5. An automatic horizon tracking device under fault control, comprising:
the preprocessing unit is used for preprocessing the layer seed point data and removing abnormal values;
the blocking unit is used for automatically blocking the preprocessed horizon seed points under the condition of fault block constraint to obtain a plurality of horizon seed point sets;
the horizon tracking unit is used for respectively carrying out horizon automatic tracking on each horizon seed point set to obtain a horizon tracking result corresponding to each block;
and the shearing unit is used for shearing the horizon tracking result by using the fault blocks to obtain the horizon tracking result corresponding to each fault block.
6. The apparatus for automatically tracking horizon under fault control according to claim 5, wherein the preprocessing unit comprises:
the segmentation module is used for dividing each horizon seed point into a plurality of segments according to the position of the invalid value of the horizon seed point;
the cutting module is used for respectively intersecting each section of the multiple sections with a fault and cutting each section into two sections from the intersection point;
and the filtering module is used for removing the section with less data in the two sections or controlling whether to keep the section with less data or not through parameters.
7. The automatic horizon tracking device under fault control according to claim 5, wherein the blocking unit comprises:
the interpolation module is used for interpolating and extrapolating the preprocessed data of the horizon seed points to obtain continuous horizon line segments;
the intersection module is used for intersecting the horizon line segment with the fault to obtain a plurality of intersection points;
the partitioning module is used for dividing the horizon line segment into a plurality of areas according to the intersection point position;
and the diversity combining module is used for dividing the preprocessed horizon seed points into a plurality of sets according to the region.
8. The automatic horizon tracking device under fault control according to claim 5, wherein the clipping unit includes:
the intersection point calculating module is used for calculating the intersection point of the horizon and the fault in the horizon tracking result;
the layer interval dividing module is used for dividing the layer position in the layer position tracking result into a plurality of layer intervals by using intersection point information;
and the position determining module is used for obtaining the position of the fault block by using the horizon seed point of the horizon in the horizon tracking result and reserving the interval data of the position.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor when executing the program implements the steps of the method for automatic horizon tracking under fault control of any of claims 1 to 4.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the automatic horizon tracking method under fault control according to any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910681480.5A CN112305596A (en) | 2019-07-26 | 2019-07-26 | Automatic horizon tracking method and device under fault control |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910681480.5A CN112305596A (en) | 2019-07-26 | 2019-07-26 | Automatic horizon tracking method and device under fault control |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112305596A true CN112305596A (en) | 2021-02-02 |
Family
ID=74329724
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910681480.5A Pending CN112305596A (en) | 2019-07-26 | 2019-07-26 | Automatic horizon tracking method and device under fault control |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112305596A (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030112704A1 (en) * | 2001-12-14 | 2003-06-19 | Goff Douglas Francis | Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume |
| CN1797039A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for automatic tracking 3D geological horizon |
| CN102819040A (en) * | 2012-08-22 | 2012-12-12 | 电子科技大学 | Three-dimensional seismic horizon automatic tracking method based on central dispersion and dip angle attribute |
| CN104181596A (en) * | 2014-08-27 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | Geologic horizon automatic tracking method and device |
| CN104375175A (en) * | 2013-08-15 | 2015-02-25 | 中国石油天然气集团公司 | Automatic horizon tracking method adopting dip angle propagation method |
| CN104597494A (en) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | Method and device for analyzing seismic formation body |
| CN104614761A (en) * | 2013-11-04 | 2015-05-13 | 中国石油天然气集团公司 | Method and device for double-layer flood fill seismic horizon space automatic tracking |
| US20160209530A1 (en) * | 2013-10-08 | 2016-07-21 | Landmark Graphics Corporation | Automated horizon auto-picking on multiple volumes |
| WO2018093531A1 (en) * | 2016-11-18 | 2018-05-24 | Landmark Graphics Corporation | Multi-z horizon auto-tracking |
-
2019
- 2019-07-26 CN CN201910681480.5A patent/CN112305596A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030112704A1 (en) * | 2001-12-14 | 2003-06-19 | Goff Douglas Francis | Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume |
| CN1797039A (en) * | 2004-12-29 | 2006-07-05 | 中国石油天然气集团公司 | Method for automatic tracking 3D geological horizon |
| CN102819040A (en) * | 2012-08-22 | 2012-12-12 | 电子科技大学 | Three-dimensional seismic horizon automatic tracking method based on central dispersion and dip angle attribute |
| CN104375175A (en) * | 2013-08-15 | 2015-02-25 | 中国石油天然气集团公司 | Automatic horizon tracking method adopting dip angle propagation method |
| US20160209530A1 (en) * | 2013-10-08 | 2016-07-21 | Landmark Graphics Corporation | Automated horizon auto-picking on multiple volumes |
| CN104597494A (en) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | Method and device for analyzing seismic formation body |
| CN104614761A (en) * | 2013-11-04 | 2015-05-13 | 中国石油天然气集团公司 | Method and device for double-layer flood fill seismic horizon space automatic tracking |
| CN104181596A (en) * | 2014-08-27 | 2014-12-03 | 中国石油集团东方地球物理勘探有限责任公司 | Geologic horizon automatic tracking method and device |
| WO2018093531A1 (en) * | 2016-11-18 | 2018-05-24 | Landmark Graphics Corporation | Multi-z horizon auto-tracking |
Non-Patent Citations (4)
| Title |
|---|
| DONALD A. HERRON等: "Horizon autopicking", PITFALLS, 31 May 2000 (2000-05-31) * |
| 殷文;郭加树;李援;朱剑兵;李长红;: "基于层控蚁群层位自动追踪的Wheeler转换方法及其应用", 中国石油大学学报(自然科学版), no. 01, 20 February 2017 (2017-02-20) * |
| 郑公营;曾婷婷;: "地震层位自动追踪技术研究", 物探化探计算技术, no. 06, 15 November 2013 (2013-11-15) * |
| 陈茂山;张雷;金瑞锋;王飞;: "适用于超大规模三维地震层位自动追踪的双尺度洪水填充方法", 石油地球物理勘探, no. 05, 15 October 2017 (2017-10-15) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108801273B (en) | Method and device for generating road reference line | |
| CN105657364B (en) | A kind of display methods of image processor, apparatus and system | |
| CN114170360B (en) | Merging and rendering method and system suitable for multiple terrain data and electronic equipment | |
| US10255716B1 (en) | Multi-resolution tiled 2.5D delaunay triangulation stitching | |
| US11327466B2 (en) | Command-value generation apparatus | |
| CN110688442A (en) | Electronic map point location aggregation method based on multilayer grid blocks | |
| CN106454155A (en) | Video shade trick processing method and device | |
| CN112305596A (en) | Automatic horizon tracking method and device under fault control | |
| CN110782774A (en) | Crowdsourcing road data distributed processing method and device | |
| US20170002643A1 (en) | Stratigraphic Correlation Method And Apparatus Based on Uncertainty | |
| CN107424583B (en) | Display data processing method and system for special-shaped image | |
| AU2013406198A1 (en) | Global grid building unfaulting sequence for complex fault-network topologies | |
| CN111637897B (en) | Map updating method, map updating device, storage medium, and processor | |
| CN114842159B (en) | Three-dimensional geologic model splicing method for complex structure of sedimentary stratum | |
| CN114202642A (en) | Method and device for filling oversized holes in three-dimensional triangulation network model | |
| CN111914372B (en) | Crack propagation depth calculation method and device and electronic equipment | |
| CN112560558B (en) | Lane line correction method and device and electronic equipment | |
| CN108599976B (en) | Network quality evaluation model construction method and device and computer equipment | |
| CN107832376B (en) | Method and device for selecting position reference point | |
| CN113076620A (en) | Water drive recoverable reserve determining method and device | |
| CN109544670B (en) | Cloud rendering method for improving lost graph | |
| CN117011479B (en) | Well modeling notch parameterization processing method and device | |
| CN105740319B (en) | A kind of method and device for correcting sampled data | |
| CN117744193B (en) | Multi-loop line breaking method, device, equipment and storage medium | |
| CN116880404B (en) | Production control method, device, equipment and medium based on constant model |
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 |