CN104166163B - Tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube - Google Patents
Tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube Download PDFInfo
- Publication number
- CN104166163B CN104166163B CN201410425325.4A CN201410425325A CN104166163B CN 104166163 B CN104166163 B CN 104166163B CN 201410425325 A CN201410425325 A CN 201410425325A CN 104166163 B CN104166163 B CN 104166163B
- Authority
- CN
- China
- Prior art keywords
- point
- tomography
- data
- curved surface
- fault
- 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
Abstract
The invention discloses a kind of tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube.It comprises the following steps: the conversion of Formica fusca volume data, data binaryzation, data de-noising process, set up initial seed point queue, set up connected component, determine tomography point, ask for tomography point normal vector, fault plane divide, fault plane secondary divide and matching structure face.The invention has the beneficial effects as follows: the present invention uses method based on space lattice distance, cross section can well be separated by the fault plane obtained, and has more integrity, can preferably match with the tomography in initial data；Use big data processing algorithm simultaneously, make the present invention realize more convenient efficiently.
Description
Technical field
The invention belongs to tomography curved surface extraction method technical field, particularly relate to a kind of based on three-dimensional big data quantity earthquake number
Tomography curved surface extraction method according to body.
Background technology
Seismic data interpretation is a very important part in geological prospecting, the explanation of its interrupting layer one of core especially.
Earth formation or rock mass can rupture when the External Force Acting sufficiently large by intensity, occur substantially then along the plane of fracture
Relative displacement, the most just define rift structure.Earthquake Faulting is the most common the most particularly significant in the earth's crust, and it generally divides
For joint and tomography two class, wherein earthquake fault can affect formation and the distribution of oil and natural gas etc., therefore explains and understands fully
The distribution of tomography all has particularly important meaning to exploration and the exploitation of the resource such as oil, natural gas.The eighties in last century it
Before, geologist is mostly to explain tomography by the way of artificial, utilizes two dimensional cross-section and the section of three-dimensional data
Information manually explains tomography.Along being perpendicular to any line direction of fault strike during the trailing of faults, or along main profile side
It is tracked to by-line, carries out contrasting, extending by tomography then along stratum vertical direction or horizontal direction, thus extend
In three dimensions.This method have that workload is big, cycle length, subjectivity strong and result can not the shortcoming such as checking repeatedly,
Its processing procedure is numerous and diverse but also easily causes bigger error.Effectively using of the method needs explanation personnel to grasp abundant ground
Matter knowledge, and interpretation process carries out manual intervention frequently, therefore this process is too dependent on explanation personnel relevant geology section
The relevant knowledge learned and rich experience.Along with the innovation of computer hardware technique and developing rapidly of image processing techniques,
Research worker attempts to apply to image processing techniques during earthquake fault explains, three dimensions detected from two dimension slicing by road
The attributive analysis of body, from manual interpretation to automatic tracing, analyzes 3D data volume from two dimensional cross-section and explains, fault recognizing skill
Art is developed on an unprecedented scale, but still also exist many difficult problems need solve, this allow for tomography identification explain remain ground
One weight difficult point of seismic exploration research field.And developing rapidly along with Computerized three-dimensional image processing techniques, engender
High-resolution coherent analysis technology so that the precision of the three-dimension disclocation identification under raising mass seismic data and efficiency become can
Can, the most therefore become a highly important tackling key problem field.Along with constantly complicating and exploration engineering of geological prospecting
Making constant progress, the identification of tomographic systems has been carried out deep with explanation by vast geology investigation and prospecting person and scientific research personnel
Research, it is proposed that more and more convenient, more practical, more careful describing and the method explaining tomography.For three
Dimension fault recognizing, in nineteen ninety-five, from the Bahorich M. and Farmer S of Amoco oil company at the 65th SEG
Coherent body technique is formally proposed in annual meeting.1999, Gersztenkorn[3]Propose Deng on the basis of existing coherent body technique
The also development for coherent body technique later of a kind of variation coherent body method based on covariance matrix provides technical support and reason
Opinion foundation.2002, Cohen etc. proposed the high-level data statistical approach of utilization and the discontinuity etc. by layer position
More accurate and effective extracting method, also the fast development for the interpretation technique of three-dimension disclocation identification later is laid a good foundation.
Randen etc. proposed the tomography automatically extracting out in seismic volume by the way of " human oasis exploited " carries out detection suppression in 2002.
First fault attributes in 3-d seismic data set is strengthened by the method, including variance attribute, dip and azimuth attribute etc.,
Then combining " human oasis exploited " for attribute character and it is carried out noise suppressed, finally combine tomography moves towards information interactive
Ground extracts tomography.The method reasonable can play compacting noise and the purpose of non-faulting response.Gibson in 2003 etc. carry
Having gone out a kind of HCF tomography automatic identifying method, the method uses coherent body to weigh the discontinuity of geological data, by advance
If seed points and threshold value from coherent body, obtain some tomography fettucelles, then by the preferential (highest of maximum confidence
Confidence first, HCF) merger strategy obtains final tomography curved surface from tomography fettucelle.2005, Dorn and
James, Tingdahl, Pierre Jacquemin decile you can well imagine out by signal processing technology, artificial neural network technology and
The method of double Hough transformations (double hough transform) realizes the automatic of tomography and semi-automatic identification.And then 2006
Year, Admasu etc. proposes that tomography is highlighted and the united mode of Active contour models extractive technique achieves the semi-automatic of tomography and chases after
Track.In the same year, Won-ki Jeong etc. proposes and uses interactive operation to carry out based on GPU (Graphics Processing Unit)
The method of fault recognizing, the method, in the process field of mass seismic data, has the biggest reference.2008,
Benjamin J etc. also been proposed interactive tomography curved surface computational methods based on level set (level Sets).The method is by level set
Computational methods, clustering technique and three-dimensional visualization technique triplicity to together with, extract tomography curved surface.Above method is followed the trail of
Cross section can not well be separated by the fault plane arrived, and can not the preferably tomography in rain initial data identical, tomography it
Between there may be focal adhesion phenomenon.When carrying out big data and processing, due to the restriction of calculator memory resource, can be to tomography
The extraction of curved surface produces impact.
Summary of the invention
In order to solve problem above, it is automatic that the present invention proposes a kind of tomography curved surface based on three-dimensional big data quantity seismic data cube
Extracting method.
The technical scheme is that a kind of tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube,
It is characterized in that, comprise the following steps:
S1. utilize ant group algorithm that seismic amplitude volume data is converted into Formica fusca volume data；
S2. the Formica fusca volume data in step S1 is carried out data binaryzation, specifically includes following steps:
S21., in the range of the maximum and minima of Formica fusca volume data, an intermediate value is chosen as attribute thresholds:
S22. some spatial point corresponding more than the data of attribute thresholds in Formica fusca volume data being set as on the tomography of space, and mark
It is designated as 1；
S23. spatial point corresponding less than the data of attribute thresholds in Formica fusca volume data is set as background dot or noise spot, and marks
It is designated as 0, obtains 0, the two-value data body of 1；
S3. utilize and open operational approach, the two-value data body obtained in step S2 is carried out denoising；
S4. the point on a tomography chosen in Formica fusca volume data in space lattice is as seed points, and this seed points is stored in
In one array, obtain an initialized seed points queue；
S5. the neighbor point near nodes for research point, and neighbor point is added the seed points queue in step S4, obtain one with
The connected component that proximity relations is corresponding, specifically includes following steps:
S51. setting search step-length is 1, according to space lattice distance and proximity relations, the neighbor point of nodes for research point；
S52. the neighbor point in step S51 is added in the seed points queue in step S4 as seed points, generate new kind
Son point queue, obtains a connected component corresponding with proximity relations, the point being on original fault plane；
S6. repeat step S4 and S5, the point on tomographies all in space is all included into corresponding connected component, determines each
Point on precursor fault face；
S7. according to the point on precursor fault face in step S6, ask on precursor fault face normal vector a little, specifically include
Following steps:
S71. using the point on precursor fault face as original point, some fields of each original point on precursor fault face are asked for
Point；
S72. the field point in step S71 is fitted to a characteristic face, ask for the normal vector of characteristic face；
S73. using the normal vector of characteristic face in step S72 as the normal vector of original point, ask on precursor fault face a little
Normal vector；
S8. poor according to the inclination angle between point each on precursor fault face, with the co-hade threshold value that sets as standard, by each
Individual precursor fault face is divided into two mutually disjoint fault planes；
S9. the fault plane obtained in step S8 is carried out secondary division process；
S10. the fault plane after dividing secondary in step S9 is fitted structure face and processes, it is achieved tomography curved surface extracts.
Further, above-mentioned steps S3 utilizes and opens operational approach, two-value data body is carried out denoising and specifically includes following step
Rapid:
S31. setting collection and be combined into A, structural element is B, utilizes structural element B that set A is carried out corrosion treatmentCorrosion Science, specifically wraps
Include following steps:
S311. the pixel in set A is contrasted by the initial point of structural element B one by one；
If S312. all pixels of structural element B are all contained in gathering in the range of A, then by the respective pixel of set A
Point retains；
If S313. all pixels of structural element B are not comprised in gathering in the range of A, then by the corresponding picture of set A
Element shop is given up；
S32. utilize structural element B that the corrosion treatmentCorrosion Science result in step S21 is carried out expansion process, specifically include following steps:
S321. about initial point, structural element B is made reflection to process, obtain structural element
S322. by the structural element in step S221Initial point with set A pixel contrast one by one；
If S323. structural elementIn pixel do not have any one point set A in the range of, then will set A in
Respective pixel point retains；
If S324. structural elementIn pixel in any one point set A in the range of, then will set A in
Respective pixel point is given up.
Further, fault plane is carried out secondary division and processes and specifically include following steps by above-mentioned steps S9:
S91. choose some initial points on fault plane, constitute precursor fault curved surface；
S92. judge whether step S91 also has unallocated complete tomography point in precursor fault curved surface；
If S93. not having unallocated complete tomography point in precursor fault curved surface, then completing secondary division and process, algorithm terminates；
If S94. precursor fault curved surface has unallocated complete tomography point, then search for point to be divided on each existing tomography
Some neighbor points；
S95. neighbor point matching in step S94 is constituted micro-plane, and calculate the offset distance of point to be divided and the micro-plane of matching；
S96. judge whether point to be divided is in the offset distance threshold range of micro-plane；
If point the most to be divided is in the offset distance threshold range of micro-plane, then add corresponding tomography song by be divided
In face, repeat step S92；
If point the most to be divided is not in the offset distance threshold range of micro-plane, then by tomography point structure obtain one new
Tomography curved surface, repeats step S92.
The invention has the beneficial effects as follows: the present invention uses method based on space lattice distance, the fault plane obtained can be well
Cross section is separated, and has more integrity, can preferably match with the tomography in initial data；Use big simultaneously
Data processing algorithm, make the present invention realize more convenient efficiently.
Accompanying drawing explanation
Fig. 1 is the tomography curved surface extraction method schematic flow sheet based on three-dimensional big data quantity seismic data cube of the present invention.
Fig. 2 be the present invention tomography between focal adhesion schematic diagram.
Fig. 3 is that set A is opened operation chart by the structural element B of the present invention.
Fig. 4 is that the operation of opening of the present invention processes selecting structure element schematic diagram.
Fig. 5 is the two-value data body schematic diagram after opening operation and processing of the present invention.
Fig. 6 is the fault branch schematic diagram of the present invention.
Fig. 7 is that the fault branch of the present invention generates singular triangular schematic diagram.
Fig. 8 is that the tomography secondary division of the present invention processes schematic diagram.
Fig. 9 is the tomography of the present invention curved surface schematic diagram of tomography without bifurcated after secondary division processes.
Figure 10 is the tomography of present invention tomography curved surface scatterplot overall situation schematic diagram after secondary division processes.
Figure 11 is the tomography of present invention tomography curved surface scatterplot partial schematic diagram after secondary division processes.
Figure 12 is the big data quantity piecemeal Processing Algorithm schematic flow sheet of the present invention.
Figure 13 is the data equal-specification piecemeal schematic diagram of the present invention.
Figure 14 is that the index of the present invention maps schematic diagram.
Detailed description of the invention
In order to make the purpose of the present invention, technical scheme and advantage clearer, below in conjunction with drawings and Examples, to this
Invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, not
For limiting the present invention.
As it is shown in figure 1, be the tomography curved surface extraction method flow process based on three-dimensional big data quantity seismic data cube of the present invention
Schematic diagram.A kind of tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube, it is characterised in that include
Following steps:
S1. utilize ant group algorithm that seismic amplitude volume data is converted into Formica fusca volume data.
S2. the Formica fusca volume data in step S1 is carried out data binaryzation, specifically includes following steps:
S21., in the range of the maximum and minima of Formica fusca volume data, an intermediate value is chosen as attribute thresholds.
According to the property distribution situation of " Formica fusca body " data, in the range of the maximum and minima of " Formica fusca body " data,
Needing to choose an intermediate value as the case may be as arranging an attribute thresholds, the attribute thresholds of setting more levels off to maximum
Value, then finally follow the trail of the tomography point obtained the fewest.
S22. some spatial point corresponding more than the data of attribute thresholds in Formica fusca volume data being set as on the tomography of space, and mark
It is designated as 1.
S23. spatial point corresponding less than the data of attribute thresholds in Formica fusca volume data is set as background dot or noise spot, and marks
It is designated as 0, obtains 0, the two-value data body of 1.
S3. utilize and open operational approach, the two-value data body obtained is carried out denoising, specifically include following step in step S2
Rapid:
S31. setting collection and be combined into A, structural element is B, utilizes structural element B that set A is carried out corrosion treatmentCorrosion Science, specifically wraps
Include following steps:
S311. the pixel in set A is contrasted by the initial point of structural element B one by one；
If S312. all pixels of structural element B are all contained in gathering in the range of A, then by the respective pixel of set A
Point retains；
If S313. all pixels of structural element B are not comprised in gathering in the range of A, then by the corresponding picture of set A
Element shop is given up；
S32. utilize structural element B that the corrosion treatmentCorrosion Science result in step S21 is carried out expansion process, specifically include following steps:
S321. about initial point, structural element B is made reflection to process, obtain structural element
S322. by the structural element in step S221Initial point with set A pixel contrast one by one；
If S323. structural elementIn pixel do not have any one point set A in the range of, then will set A in
Respective pixel point retains；
If S324. structural elementIn pixel in any one point set A in the range of, then will set A in
Respective pixel point is given up.
As in figure 2 it is shown, be focal adhesion schematic diagram between the tomography of the present invention.Present invention employing is opened operation process and is eliminated disconnected
Focal adhesion phenomenon between Ceng.As it is shown on figure 3, set A is opened operation chart by the structural element B for the present invention.
Image is held operation typically to disconnect narrower interruption, eliminate tapering burr, and image outline is smoothed.As
Shown in Fig. 4, the operation of opening for the present invention processes selecting structure element schematic diagram.As it is shown in figure 5, for the present invention through opening behaviour
Two-value data body schematic diagram after dealing with.The two-value data body before and after operation processes is opened it can be seen that utilize and open by contrast
Operation can actually remove noise in the case of keeping image original configuration feature, eliminates the local between tomography and glues
Even phenomenon.
S4. the point on a tomography chosen in Formica fusca volume data in space lattice is as seed points, and this seed points is stored in
In one array, obtain an initialized seed points queue.
S5. the neighbor point near nodes for research point, and neighbor point is added the seed points queue in step S4, obtain one with
The connected component that proximity relations is corresponding, specifically includes following steps:
S51. setting search step-length is 1, according to space lattice distance and proximity relations, the neighbor point of nodes for research point.
Determine that step-size in search, step-size in search are defaulted as 1, by suitable regulation step-size in search, the fault plane of tracking can be made more
Add complete.Search is till all of seed points all cannot find the seed points made new advances in the range of step-size in search again.
S52. the neighbor point in step S51 is added in the seed points queue in step S4 as seed points, generate new kind
Son point queue, obtains a connected component corresponding with proximity relations, the point being on original fault plane.
S6. repeat step S4 and S5, the point on tomographies all in space is all included into corresponding connected component, determines each
Point on precursor fault face.
S7. according to the point on precursor fault face in step S6, ask on precursor fault face normal vector a little, specifically include
Following steps:
S71. using the point on precursor fault face as original point, some fields of each original point on precursor fault face are asked for
Point.
S72. the field point in step S71 is fitted to a characteristic face, ask for the normal vector of characteristic face.
S73. using the normal vector of characteristic face in step S72 as the normal vector of original point, ask on precursor fault face a little
Normal vector.
S8. poor according to the inclination angle between point each on precursor fault face, with the co-hade threshold value that sets as standard, by each
Individual precursor fault face is divided into two mutually disjoint fault planes.
Here tilt threshold≤45 degree.
S9. the fault plane obtained in step S8 is carried out secondary division process, specifically includes following steps:
S91. choose some initial points on fault plane, constitute precursor fault curved surface.
S92. judge whether step S91 also has unallocated complete tomography point in precursor fault curved surface.
If S93. not having unallocated complete tomography point in precursor fault curved surface, then completing secondary division and process, algorithm terminates.
If S94. precursor fault curved surface has unallocated complete tomography point, then search for point to be divided on each existing tomography
Some neighbor points.
S95. neighbor point matching in step S94 is constituted micro-plane, and calculate the offset distance of point to be divided and the micro-plane of matching.
S96. judge whether point to be divided is in the offset distance threshold range of micro-plane.
If point the most to be divided is in the offset distance threshold range of micro-plane, then add corresponding tomography song by be divided
In face, repeat step S92.
If point the most to be divided is not in the offset distance threshold range of micro-plane, then by tomography point structure obtain one new
Tomography curved surface, repeats step S92.
As shown in Figure 6, for the fault branch schematic diagram of the present invention.As it is shown in fig. 7, be that the fault branch of the present invention generates strange
Different triangle schematic diagram.As shown in Figure 8, the tomography secondary division for the present invention processes schematic diagram.As it is shown in figure 9, be this
The tomography of the invention curved surface schematic diagram of tomography without bifurcated after secondary division processes.As shown in Figure 10, for the tomography warp of the present invention
Tomography curved surface scatterplot overall situation schematic diagram after secondary division process.As shown in figure 11, for the tomography of the present invention at secondary division
Tomography curved surface scatterplot partial schematic diagram after reason.Tomography is carried out using part plan matching when secondary divides by the present invention
Mode, splits into some tomography fettucelles by the tomography curved surface of those complexity so that each fettucelle is approximately a space
Plane, solves the bifurcation problem of phantom with this.
S10. the fault plane after dividing secondary in step S9 is fitted structure face and processes, it is achieved tomography curved surface extracts.
The present invention uses big data quantity piecemeal Processing Algorithm.As shown in figure 12, for the big data quantity piecemeal Processing Algorithm of the present invention
Schematic flow sheet.Whole piecemeal Processing Algorithm can be divided into two big steps, the first first step to be the piecemeals to big data quantity
Process, then second step is exactly that the data after piecemeal are carried out scheduling memory, specifically includes following steps:
Work area is divided into the little rectangular block of equal size by step 1..
The present invention uses equal-specification dividing mode, and big data quantity carries out the division of same size size, it is ensured that every after division
One little tile data size is identical.As shown in figure 13, for the data equal-specification piecemeal schematic diagram of the present invention.
Step 2. sets up mark file with rectangular block for ultimate unit.
Step 3. sets up rectangular block index information.
Step 4. opens up buffer area at internal memory.
Step 5. judges whether required identification information reads in caching.
If the required identification information of step 6. has read in caching, then directly reading cache information, operation terminates.
If the required identification information of step 7. does not reads in caching, then judge that caching is the fullest；
If step 8. caching less than, then directly identification information block is called in caching, reads cache information, operation terminates；
If step 9. caching is the fullest, then find the rectangular block information being not used by the most at most；
Step 10. updates file, and information needed is called in caching, replaces old rectangular block, and operation terminates.
The present invention uses the most untapped replacement algorithm in the paged memory management mechanism of computer operating system,
The data block being not used at most in every time carrying out all replacing nearest a period of time during data block replacement.Algorithm running
In, when needing the label information accessing certain mesh point, initially to the labelling finding this mesh point place in internal memory cache region
Data block, if found, directly accesses, simultaneously need to the nearest use time of accessed data block is set to 0, otherwise just needs
Corresponding mark data block is called in internal memory.When carrying out data block scheduling memory, first determine whether that internal memory cache region is the most
Full, if less than, then directly desired data block is called in internal memory, the nearest use time of this data block is set to 0 simultaneously,
And the time that uses recently of other data with existing blocks in internal memory cache region is increased 1 certainly.If buffer area is the fullest, then choose buffer area
In use recently the maximum data block of time value, the data block that it is not the most accessed in representing nearest a period of time,
After the information updating that will comprise in this data block to file, then displaced by new data block, simultaneously by new data
The nearest use time of block is set to 0.
Those of ordinary skill in the art is it will be appreciated that embodiment described here is to aid in the reader understanding present invention's
Principle, it should be understood that protection scope of the present invention is not limited to such special statement and embodiment.This area common
It is various specifically that technical staff can make various other without departing from essence of the present invention according to these technology disclosed by the invention enlightenment
Deformation and combination, these deformation and combination are the most within the scope of the present invention.
Claims (3)
1. a tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube, it is characterised in that comprise the following steps:
S1. utilize ant group algorithm that seismic amplitude volume data is converted into Formica fusca volume data；
S2. the Formica fusca volume data in step S1 is carried out data binaryzation, specifically includes following steps:
S21., in the range of the maximum and minima of Formica fusca volume data, an intermediate value is chosen as attribute thresholds；
S22. some spatial point corresponding more than the data of attribute thresholds in Formica fusca volume data being set as on the tomography of space, and it is labeled as 1；
S23. spatial point corresponding less than the data of attribute thresholds in Formica fusca volume data is set as background dot or noise spot, and is labeled as 0, obtain 0, the two-value data body of 1；
S3. utilize and open operational approach, the two-value data body obtained in step S2 is carried out denoising；
S4. the point on a tomography chosen in Formica fusca volume data in space lattice is as seed points, and this seed points is stored in an array, obtains an initialized seed points queue；
S5. the neighbor point near nodes for research point, and neighbor point adds the seed points queue in step S4, obtains a connected component corresponding with proximity relations, specifically includes following steps:
S51. setting search step-length is 1, according to space lattice distance and proximity relations, the neighbor point of nodes for research point；
S52. the neighbor point in step S51 is added in the seed points queue in step S4 as seed points, generate new seed points queue, obtain a connected component corresponding with proximity relations, the point being on original fault plane；
S6. repeat step S4 and S5, the point on tomographies all in space is all included into corresponding connected component, determines the point on each precursor fault face；
S7. according to the point on precursor fault face in step S6, ask on precursor fault face normal vector a little, specifically include following steps:
S71. using the point on precursor fault face as original point, some fields point of each original point on precursor fault face is asked for；
S72. the field point in step S71 is fitted to a characteristic face, ask for the normal vector of characteristic face；
S73. using the normal vector of characteristic face in step S72 as the normal vector of original point, ask on precursor fault face normal vector a little；
S8. poor according to the inclination angle between point each on precursor fault face, with the co-hade threshold value that sets as standard, each precursor fault face is divided into two mutually disjoint fault planes；
S9. the fault plane obtained in step S8 is carried out secondary division process；
S10. the fault plane after dividing secondary in step S9 is fitted structure face and processes, it is achieved tomography curved surface extracts.
2. as claimed in claim 1 tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube, it is characterised in that: described step S3 utilizes and opens operational approach, two-value data body carries out denoising and specifically includes following steps:
S31. setting collection and be combined into A, structural element is B, utilizes structural element B that set A is carried out corrosion treatmentCorrosion Science, specifically includes following steps:
S311. the pixel in set A is contrasted by the initial point of structural element B one by one；
If S312. all pixels of structural element B are all contained in gathering in the range of A, then the respective pixel point of set A is retained；
If S313. all pixels of structural element B are not comprised in gathering in the range of A, then the respective pixel point of set A is given up；
S32. utilize structural element B that the corrosion treatmentCorrosion Science result in step S31 is carried out expansion process, specifically include following steps:
S321. about initial point, structural element B is made reflection to process, obtain structural element
S322. by the structural element in step S321Initial point with set A pixel contrast one by one；
If S323. structural elementIn pixel do not have any one point set A in the range of, then will set A in respective pixel point retain；
If S324. structural elementIn pixel in any one point set A in the range of, then by set A in respective pixel point give up.
3. as claimed in claim 1 tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube, it is characterised in that: fault plane is carried out secondary division and processes and specifically include following steps by described step S9:
S91. choose some initial points on fault plane, constitute precursor fault curved surface；
S92. judge whether step S91 also has unallocated complete tomography point in precursor fault curved surface；
If S93. not having unallocated complete tomography point in precursor fault curved surface, then completing secondary division and process, algorithm terminates；
If S94. precursor fault curved surface has unallocated complete tomography point, then search for the point to be divided some neighbor points on each existing tomography；
S95. neighbor point matching in step S94 is constituted micro-plane, and calculate the offset distance of point to be divided and the micro-plane of matching；
S96. judge whether point to be divided is in the offset distance threshold range of micro-plane；
If point the most to be divided is in the offset distance threshold range of micro-plane, then adds in corresponding tomography curved surface by be divided, repeat step S92；
If point the most to be divided is not in the offset distance threshold range of micro-plane, is then obtained a new tomography curved surface by tomography point structure, repeat step S92.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410425325.4A CN104166163B (en) | 2014-08-27 | 2014-08-27 | Tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410425325.4A CN104166163B (en) | 2014-08-27 | 2014-08-27 | Tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104166163A CN104166163A (en) | 2014-11-26 |
| CN104166163B true CN104166163B (en) | 2016-08-24 |
Family
ID=51910062
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410425325.4A Active CN104166163B (en) | 2014-08-27 | 2014-08-27 | Tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104166163B (en) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104570096B (en) * | 2015-01-22 | 2017-09-15 | 中国海洋石油总公司 | A kind of granite buried hill top surface fault interpretation method |
| CN105182415B (en) * | 2015-09-02 | 2018-02-13 | 中国石油化工股份有限公司 | A kind of fault recognizing method based on scaling down processing |
| CN105425296B (en) * | 2015-10-28 | 2018-04-06 | 中国石油天然气集团公司 | Geologic body recognition methods and device |
| CN106199727B (en) * | 2016-08-30 | 2018-07-10 | 中国华能集团清洁能源技术研究院有限公司 | A kind of method using seismic interval velocity identification shale gas dessert |
| CN107368839B (en) * | 2017-06-22 | 2020-05-12 | 南京师范大学 | Automatic fault plane extraction method based on DEM |
| CN108037527A (en) * | 2017-12-07 | 2018-05-15 | 中国石油大学(华东) | A kind of complicated fault identification and detection method based on magic square operator |
| CN109086773B (en) * | 2018-08-29 | 2022-03-04 | 电子科技大学 | Fault plane identification method based on full convolution neural network |
| CN110174700B (en) * | 2019-05-16 | 2021-01-15 | 中海石油(中国)有限公司 | Seismic attribute boundary line enhancement method for simulating root growth |
| CN112130200B (en) * | 2020-09-23 | 2021-07-20 | 电子科技大学 | Fault identification method based on grad-CAM attention guidance |
| CN112505758B (en) * | 2020-11-17 | 2021-10-26 | 中国石油集团工程咨询有限责任公司 | Method for processing seismic fault image of complex geological structure based on fault bifurcation structure model |
| CN113267815B (en) * | 2021-07-07 | 2022-05-10 | 中海油田服务股份有限公司 | Method and device for filtering repeated broken edge data |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222365A (en) * | 2011-07-29 | 2011-10-19 | 电子科技大学 | Method for reconstructing curved surface of complex space |
| CN102222366A (en) * | 2011-07-29 | 2011-10-19 | 电子科技大学 | Method for fitting complex space curved surfaces |
| CN102867330A (en) * | 2012-08-29 | 2013-01-09 | 电子科技大学 | Region-division-based spatial complex horizon reconstruction method |
| CN103412331A (en) * | 2013-08-30 | 2013-11-27 | 电子科技大学 | Automatic extraction method for three-dimensional earthquake fault |
| CN103489222A (en) * | 2013-09-06 | 2014-01-01 | 电子科技大学 | Target body surface reconstruction method in three-dimensional image |
-
2014
- 2014-08-27 CN CN201410425325.4A patent/CN104166163B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222365A (en) * | 2011-07-29 | 2011-10-19 | 电子科技大学 | Method for reconstructing curved surface of complex space |
| CN102222366A (en) * | 2011-07-29 | 2011-10-19 | 电子科技大学 | Method for fitting complex space curved surfaces |
| CN102867330A (en) * | 2012-08-29 | 2013-01-09 | 电子科技大学 | Region-division-based spatial complex horizon reconstruction method |
| CN103412331A (en) * | 2013-08-30 | 2013-11-27 | 电子科技大学 | Automatic extraction method for three-dimensional earthquake fault |
| CN103489222A (en) * | 2013-09-06 | 2014-01-01 | 电子科技大学 | Target body surface reconstruction method in three-dimensional image |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104166163A (en) | 2014-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104166163B (en) | Tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube | |
| CN103017739B (en) | Manufacturing method of true digital ortho map (TDOM) based on light detection and ranging (LiDAR) point cloud and aerial image | |
| CN102520401B (en) | Building area extraction method based on LiDAR data | |
| CN104049245B (en) | Urban building change detection method based on LiDAR point cloud spatial difference analysis | |
| CN102521884B (en) | Three-dimensional roof reconstruction method based on LiDAR data and ortho images | |
| CN110263111B (en) | Priori knowledge-based space-time monitoring method for land utilization/coverage information | |
| CN103268358B (en) | Multi-source control point image database builds and update method | |
| CN104899562B (en) | Radar remote sensing image culture's recognizer based on Texture Segmentation fusion | |
| CN109709603A (en) | Seismic horizon identification and method for tracing, system | |
| CN105550691B (en) | Adaptive mountain valley ridge line extracting method and system based on scale space | |
| CN104036550B (en) | Laser radar point-cloud interpretation and reconstruction method for building elevations on basis of shape semantics | |
| CN102609982B (en) | Topology discovery method of space geological data based on unstructured mode | |
| CN104376595A (en) | Method for generating three-dimensional roads on basis of airborne LiDAR (light detection and ranging) and GIS (geographic information system) cooperation | |
| CN109872311A (en) | A kind of Rock Mass Integrality sentences knowledge method | |
| CN106709987B (en) | Dynamic construction method of three-dimensional geological profile model | |
| CN113762090A (en) | Disaster monitoring and early warning method for extra-high voltage intensive power transmission channel | |
| CN111754623A (en) | Discrete fracture modeling method based on multi-point geostatistics | |
| CN103035006A (en) | High-resolution aerial image partition method based on LEGION and under assisting of LiDAR | |
| CN103797382A (en) | Method and system of subsurface horizon assignment | |
| CN105719344A (en) | Generation method three-dimensional land utilization present situation data model | |
| CN112489208A (en) | Crack sheet extraction method based on ant algorithm and three-dimensional geological model construction method | |
| CN106500594A (en) | Fusion reflected intensity and the railroad track method for semi-automatically detecting of geometric properties | |
| CN104951752A (en) | Method for extracting houses from airborne laser point cloud data | |
| Ren et al. | Mapping High-Resolution Global Impervious Surface Area: Status and Trends | |
| WO2018196214A1 (en) | Statistics system and statistics method for geographical influence on vernacular architectural form |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201125 Address after: 610000 Sichuan city of Chengdu province Tianfu Tianfu Avenue South Huayang Street No. 846 Patentee after: CHENGDU AIWEI BEISI TECHNOLOGY Co.,Ltd. Address before: 611731 Chengdu province high tech Zone (West) West source Avenue, No. 2006 Patentee before: UNIVERSITY OF ELECTRONIC SCIENCE AND TECHNOLOGY |