CN109143334A - A kind of layer of position interpolation method, apparatus and system - Google Patents
A kind of layer of position interpolation method, apparatus and system Download PDFInfo
- Publication number
- CN109143334A CN109143334A CN201810796224.6A CN201810796224A CN109143334A CN 109143334 A CN109143334 A CN 109143334A CN 201810796224 A CN201810796224 A CN 201810796224A CN 109143334 A CN109143334 A CN 109143334A
- Authority
- CN
- China
- Prior art keywords
- layer
- interpolation
- seed point
- site
- grid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/643—Horizon tracking
Abstract
This specification embodiment discloses a kind of layer of position interpolation method, apparatus and system, and the method includes carrying out gridding processing to layer position data according to default layer position sizing grid；Obtain the layer position grid vertex for meeting preset threshold condition with interpolation layer site distance；The preparation seed point in interpolation layer site is searched in the adjacent layer position grid centered on the vertex, obtains the first preparation seed point set；The line and tomography polygon for judging the preparation seed point and interpolation layer site that the first preparation seed point is concentrated with the presence or absence of intersection point, if it is judged that be it is no, by the corresponding seed point prepared seed point and be determined as interpolation layer site；The value that interpolation calculation obtains interpolation layer site is carried out according to the seed point.Using each embodiment of this specification, the efficiency and accuracy of layer position interpolation can be improved.
Description
Technical field
The present invention relates to seismic exploration technique fields, particularly, are related to a kind of layer of position interpolation method, apparatus and system.
Background technique
Layer position interpolation is one of common method in seismic interpretation, is carried out by layer position interpolation along layer attributes extraction, layer position
Operation, at figure etc., the key of layer position interpolation is interpolation precision and efficiency.At present mainstream interpretation software (Landmark,
GeoFrame etc.) layer position interpolation algorithm tomography blocking is not accounted for when interpolation, after the completion of interpolation, then polygon with tomography
Shape shearing.Layer position is caused not extend on tomography according to actual trend, but the layer position of tomography two sides is connected directly, and is caused
Near Faults layer position interpolation inaccuracy.And the oil gas of many geological conditions is mainly gathered near Faults, therefore not using this
Accurate layer position, which carries out attributes extraction, inverting, oil and gas detection etc., will certainly cause result inaccurate, be not able to satisfy oil field prospecting and open
Hair demand.
Summary of the invention
This specification embodiment is designed to provide a kind of layer of position interpolation method, apparatus and system, and layer position can be improved
The accuracy and efficiency of interpolation.
This specification provides a kind of layer of position interpolation method, apparatus and system includes that under type such as is realized:
A kind of layer of position interpolation method, comprising:
Gridding processing is carried out to layer position data according to default layer position sizing grid；
Obtain the layer position grid vertex for meeting preset threshold condition with interpolation layer site distance；
The preparation seed point that interpolation layer site is searched in the adjacent layer position grid centered on the vertex, is obtained
Obtain the first preparation seed point set；
Judge that the line and tomography of the preparation seed point and interpolation layer site that the first preparation seed point is concentrated are polygon
Shape whether there is intersection point, if it is judged that be it is no, corresponding seed point of preparing is determined as the kind in interpolation layer site
Sub- point；
The value that interpolation calculation obtains interpolation layer site is carried out according to the seed point.
In another embodiment of the method that this specification provides, judgement the first preparation seed point is concentrated
Preparation seed point and interpolation layer site line and tomography polygon whether there is intersection point, comprising:
Gridding processing is carried out to layer data according to layer bit interval size, obtains what the first preparation seed point was concentrated
The tomography grid that the line in preparation seed point and interpolation layer site passes through；
Tomography number in the tomography grid is merged and obtains tomography number union；
The tomography number for judging that the line in the prepared seed point and interpolation layer site is numbered and concentrated with tomography is corresponding
Tomography polygon whether there is intersection point.
In another embodiment of the method that this specification provides, institute is described to obtain the first preparation seed point set
In preparation seed point and interpolation layer site the tomography grid that passes through of line, comprising:
The layer site that record seismic survey lines intersect with tomography polygon；
It removes the layer site that seismic survey lines described in the first preparation seed point collection intersect with tomography polygon, obtains the
Two preparation seed point sets；
Obtain the tomography that the line of preparation seed point and interpolation layer site that the second preparation seed point is concentrated passes through
Grid.
In another embodiment of the method that this specification provides, the basis presets layer position sizing grid to layer position
Data carry out gridding processing, comprising:
Gridding processing is carried out to layer position data according to the layer position interpretation interval size of half.
In another embodiment of the method that this specification provides, the preparation in interpolation layer site described in described search
Before seed point, further includes: delete the layer position data in closed fault polygon.
In another embodiment of the method that this specification provides, the preparation in interpolation layer site described in described search
Before seed point, further includes:
Interpolation calculation is carried out to the layer position grid vertex of no value, obtains the value of layer position grid vertex.
This specification embodiment also provides a kind of layer of position interpolation method, comprising:
Gridding processing is carried out to layer position data according to default layer position sizing grid；
The rectangular extent that a layer position sizing grid is determined centered on the grid vertex of layer position, obtains in the rectangular extent
Interpolation layer site obtain value cell to be inserted；
The pre- of value cell to be inserted is searched in the adjacent layer position grid centered on the grid vertex of the layer position to make preparations for sowing
It is sub-, obtain third preparation seed point set；
The value in interpolation layer site in value cell to be inserted is calculated according to the third preparation seed point set, comprising:
Judge the interpolation layer site in the preparation seed point and value cell to be inserted of the third preparation seed point concentration
Line and tomography polygon whether there is intersection point, if it is judged that be it is no, by corresponding preparation seed point be determined as correspondence to
The seed point in interpolated layer site carries out the value that interpolation calculation obtains corresponding interpolation layer site according to the seed point.
It is described to calculate interpolation layer position in value cell to be inserted in another embodiment of the method that this specification provides
The value of point, comprising:
Using multithreads computing, the value in interpolation layer site in each value cell to be inserted is calculated separately.
This specification embodiment also provides a kind of layer of position interpolating apparatus, comprising:
First layer position gridding module, for carrying out gridding processing to layer position data according to default layer position sizing grid；
First search module, for obtaining the layer position grid top for meeting preset threshold condition with interpolation layer site distance
Point, searches for the preparation seed point in interpolation layer site in the adjacent layer position grid centered on the vertex, obtains the
One preparation seed point set；
Seed point determining module, for judging preparation seed point and the interpolation layer position of the first preparation seed point concentration
Point line and tomography polygon whether there is intersection point, if it is judged that be it is no, by it is corresponding prepare seed point be determined as institute
State the seed point in interpolation layer site
Interpolation calculation module, for carrying out the value that interpolation calculation obtains interpolation layer site according to the seed point.
This specification embodiment also provides a kind of layer of position interpolating apparatus, comprising:
Second layer position gridding module, for carrying out gridding processing to layer position data according to default layer position sizing grid；
Interpolating unit division module, for determining the rectangle model of a layer position sizing grid centered on the grid vertex of layer position
It encloses, the interpolation layer site obtained in the rectangular extent obtains value cell to be inserted；
Second search module, for searching for the value cell to be inserted in the adjacent layer position grid centered on the vertex
Preparation seed point, obtain third preparation seed point set；
Concurrent operation module calculates separately interpolation layer position in each value cell to be inserted for utilizing multithreads computing
The value of point, and seed point determination unit and interpolation calculation unit including multiple concurrent operations, wherein
The seed point determination unit, for judging the preparation seed point and interpolation of the third preparation seed point concentration
The line in the interpolation layer site in unit and tomography polygon whether there is intersection point, if it is judged that be it is no, will be corresponding
Preparation seed point is determined as the seed point in corresponding interpolation layer site；
The interpolation calculation unit obtains corresponding interpolation layer site for carrying out interpolation calculation according to the seed point
Value.
The layer position interpolating apparatus of this specification embodiment, including processor and depositing for storage processor executable instruction
Reservoir, when described instruction is executed by the processor realize the following steps are included:
Gridding processing is carried out to layer position data according to default layer position sizing grid；
Obtain the layer position grid vertex for meeting preset threshold condition with interpolation layer site distance；
The preparation seed point that interpolation layer site is searched in the adjacent layer position grid centered on the vertex, is obtained
Obtain the first preparation seed point set；
Judge that the line and tomography of the preparation seed point and interpolation layer site that the first preparation seed point is concentrated are polygon
Shape whether there is intersection point, if it is judged that be it is no, corresponding seed point of preparing is determined as the kind in interpolation layer site
Sub- point；
The value that interpolation calculation obtains interpolation layer site is carried out according to the seed point.
This specification embodiment also provides a kind of layer of position interplotation system, including at least one processor and storage computer
The memory of executable instruction, the processor realize method described in any one above-mentioned embodiment when executing described instruction
Step.
A kind of layer of position interpolation method, the apparatus and system that this specification one or more embodiment provides, can pass through head
Gridding processing first is carried out to layer position data, is determined by the splicing of four grid nearest apart from interpolation layer site distance
Preparation seed point, substantially increases the efficiency for searching for preparation seed point.Then, faults control is further increased to determine seed
Point, to substantially increase the accuracy of interpolation calculation result.Using the scheme of this specification embodiment, layer position can be improved and insert
The accuracy and efficiency of value.
Detailed description of the invention
In order to illustrate more clearly of this specification embodiment or technical solution in the prior art, below will to embodiment or
Attached drawing needed to be used in the description of the prior art is briefly described, it should be apparent that, the accompanying drawings in the following description is only
The some embodiments recorded in this specification, for those of ordinary skill in the art, in not making the creative labor property
Under the premise of, it is also possible to obtain other drawings based on these drawings.In the accompanying drawings:
Fig. 1 is the flow diagram for a kind of layer of position interpolation method embodiment that this specification provides；
Fig. 2 is the flow diagram for another layer position interpolation method embodiment that this specification provides；
Line near Faults interpolation result contrast schematic diagram has been explained in the example that Fig. 3 provides for this specification；
Line near Faults interpolation result contrast schematic diagram is not explained in the example that Fig. 4 provides for this specification；
Fig. 5 is the modular structure schematic diagram for a kind of layer of position interpolating apparatus embodiment that this specification provides；
Fig. 6 is the modular structure schematic diagram for another layer position interpolating apparatus embodiment that this specification provides.
Specific embodiment
In order to make those skilled in the art more fully understand the technical solution in this specification, below in conjunction with this explanation
Attached drawing in book one or more embodiment carries out the technical solution in this specification one or more embodiment clear, complete
Site preparation description, it is clear that described embodiment is only specification a part of the embodiment, instead of all the embodiments.Based on saying
Bright book one or more embodiment, it is obtained by those of ordinary skill in the art without making creative efforts all
The range of this specification example scheme protection all should belong in other embodiments.
Layer position interpolation algorithm does not account for tomography blocking when interpolation at present, but uses tomography again after the completion of interpolation
Polygon shearing.Layer position is caused not extend on tomography according to actual trend, the layer position of tomography two sides is connected directly, and is caused
Near Faults layer position interpolation inaccuracy.In view of the above technical problems, this specification embodiment provides a kind of layer of position interpolation method,
Gridding processing is carried out to layer position data first, may include multiple layers of position data in each layer of position grid, then utilizes distance
The splicing of four nearest grids of interpolation layer site distance substantially increases to determine prepared seed point and searches for preparation seed point
Efficiency.Later, faults control is further increased to determine seed point, to make full use of useful number in interpolation calculation
According to effectivelying prevent influence of another side data of tomography to calculated result, the interpolation processing of the realization of efficiently and accurately to target work area.
Using the scheme of this specification embodiment, the accuracy and efficiency of layer position interpolation can be improved.
Fig. 1 is a kind of layer of position interpolation method embodiment flow diagram that this specification provides.Although this specification
It provides as the following examples or method operating procedure shown in the drawings or apparatus structure, but based on conventional or without creativeness
Labour may include more in the method or device or part merge after less operating procedure or modular unit.?
In the step of there is no necessary causalities in logicality or structure, the modular structure of the execution sequences of these steps or device is not
It is limited to this specification embodiment or execution shown in the drawings sequence or modular structure.The method or modular structure in reality
In device, server or end product in application, can according to embodiment or method shown in the drawings or modular structure into
Row sequence execute or it is parallel execute (such as parallel processor or multiple threads environment, even include distributed treatment,
The implementation environment of server cluster).
Specific one embodiment is as shown in Figure 1, in the one embodiment for the layer position interpolation method that this specification provides, institute
The method of stating may include:
S102: gridding processing is carried out to layer position data according to default layer position sizing grid.
Gridding processing can be carried out to layer position data, correspondingly, may include multiple layers of digit in each layer of position grid
According to.Such as, the layer position data in certain target work area can be projected in rectangular coordinate system, so as to pass through in rectangular coordinate system
The value in discrete layer site reflects the spatial distributions of target work area layer position data.Wherein, in some embodiments of this specification, institute
It states the shortest distance between layer site and is properly termed as a layer bit interval.Layer position data in rectangular coordinate system are carried out at gridding
Reason, correspondingly, may include multiple layers of site in each layer of position grid.
In some embodiments, gridding processing can be carried out to layer position data according to default layer position sizing grid, it is described
Default layer position sizing grid can be obtained according to real data and empirical analysis.Preferably, in one embodiment of this specification,
Gridding processing can be carried out to layer position data according to the layer position interpretation interval size of half.Layer position interpretation interval can
To include the distance between target work area survey line.For example, can be using the half slightly larger than layer position interpretation interval as layer position grid
Size so that the grid cell that adjacent four layer positions grid by vertex of single layer position grid vertex forms has
It imitates coating position and explains the actual measurement layer position data obtained, improve the accuracy of succeeding layer position interpolation.
S104: the layer position grid vertex for meeting preset threshold condition with interpolation layer site distance is obtained, with the top
The preparation seed point in interpolation layer site is searched in adjacent layer position grid centered on point, obtains the first preparation seed point
Collection.
Interpolation layer site may include the layer site of no value, and the prepared seed point may include having the layer position of value
Point.The available layer position grid vertex for meeting preset threshold condition with interpolation layer site distance, can such as calculate interpolation
Layer site obtains at a distance from each grid vertex apart from the corresponding grid vertex of minimum value.It is then possible to this layer of position grid top
Adjacent four layer position grids centered on point are as a grid cell, and search has the layer of value in this adjacent four layer positions grid
Preparation seed point of the site as the interpolation layer site obtains the first preparation seed point set in corresponding interpolation layer site.Into
One step, if searched in adjacent four layer positions grid preparation seed points it is less, influence final interpolation as a result, can
Further to be searched in adjacent layer position grid, to obtain the preparation seed point that can satisfy interpolation needs.
In another embodiment of this specification, can layer position grid vertex to no value carry out interpolation pretreatment and obtain layer
The value of position grid vertex.Correspondingly, the prepared seed point may include layer site and the layer position grid vertex for having value.It utilizes
The scheme of the present embodiment, it is ensured that the preparation seed point in each layer of position grid is enough, reduce search again for it is neighbouring other
Preparation seed point probability in the grid of layer position, further increases efficiency.
S106: judge the line and tomography of the preparation seed point and interpolation layer site that the first preparation seed point is concentrated
Polygon whether there is intersection point, if it is judged that be it is no, by it is corresponding prepare seed point be determined as interpolation layer site
Seed point.
The seed point in corresponding interpolation layer site can be filtered out from preparation seed point according to following judgment steps:
Some the preparation seed point for obtaining interpolation layer site, judges the line of the preparation seed point and interpolation layer site
It whether there is intersection point with each tomography.Such as fault recognizing first can be carried out to target work area, draw tomography polygon, tomography is polygon
Graphic data projects under the rectangular coordinate system of layer position data.Judge the line and each tomography in preparation seed point and interpolation layer site
Polygon whether there is intersection point.
If it is judged that be no, i.e., no intersection point, then it can be using the preparation seed point as the kind in interpolation layer site
Sub- point；If there is intersection point, which is cast out.
The whole preparation seed point for treating interpolated layer site repeats above-mentioned judgment step, so as to filter out interpolation layer
The seed point in site.
In another embodiment of this specification, the line and tomography for judging preparation seed point and interpolation layer site
Polygon whether there is intersection point, can also include:
S1062: gridding processing is carried out to layer data according to layer bit interval size, obtains the first preparation seed point
The tomography grid that the preparation seed point of concentration and the line in interpolation layer site pass through；
S1064: the tomography number in the tomography grid is merged and obtains tomography number union；
S1066: judge that the line of the prepared seed point and interpolation layer site and tomography are numbered and the tomography concentrated is compiled
Number corresponding tomography polygon whether there is intersection point.
Tomography polygon can be numbered in advance, using layer bit interval size as tomography sizing grid to layer position right angle
Layer data under coordinate system carries out gridding processing, obtains tomography grid.Then, preparation seed point and interpolation layer position are obtained
The tomography grid that the line of point passes through, and obtain the number of the tomography in corresponding tomography grid.It will be each disconnected in corresponding tomography grid
Layer number merges, and calculates the union of the tomography number after merging, obtains tomography number union.
Then, it is numbered using tomography and each tomography number is concentrated to filter out corresponding tomography polygon, then with preparing seed
Point and the line in interpolation layer site number with tomography and each tomography are concentrated to number corresponding tomography polygon intersection, and judgement is
It is no that there are intersection points.If the preparation seed point to be all determined as to the seed point in interpolation layer site without intersection point；Otherwise, cast out
The preparation seed point.Treat the corresponding all preparations seed points in interpolated layer site and repeat above-mentioned judgment step, thus filter out to
The seed point set in interpolated layer site.
Using the union for the tomography grid interrupting layer number that the line of preparation seed point and interpolation layer site passes through, come true
Determine the tomography polygon for asking friendship, it can be while reducing the tomography polygon quantity for asking friendship, it can also be ensured that use
In ask friendship tomography polygon will not repeat and will not miss.Meanwhile the size of tomography grid is identical with layer bit interval, it can be with
It further increases and merges the simplicity that tomography grid obtains tomography number union.It, can be further using the scheme of the present embodiment
Quickly determine that may influence subsequent interpolation calculation obtains tomography polygon, to improve the efficiency of screening seed point and accurate
Property.
In another embodiment of this specification, the preparation seed point of the first preparation seed point concentration and to be inserted is obtained
The tomography grid that the line in value layer site passes through may include: the layer site for recording seismic survey lines and intersecting with tomography polygon, so
Afterwards, the layer site that seismic survey lines described in the first preparation seed point collection intersect with tomography polygon is removed, it is pre- to obtain second
Make preparations for sowing sub- point set, and the line that obtains preparation seed point and interpolation layer site that the second preparation seed point is concentrated passes through
Tomography grid.
To layer data carry out gridding processing when, can recorde seismic survey lines (such as each main profile (Inline) and
Cross-track (Crossline)) the layer site intersected with tomography polygon.It is then possible to first remove the first preparation seed point set
In include the layer site intersected with tomography of the seismic survey lines, obtain the second preparation seed point set.Then, according to the second preparation
Seed point set determines the seed point in corresponding interpolation layer site.By remembering while carrying out gridding processing to layer data
The layer site that lower seismic survey lines intersect with tomography is recorded, and part layer site is filtered.It is equivalent to filter out in advance and partially be discontented with
The preparation seed point required enough, so as to further increase the efficiency that subsequent seed point determines.
In another embodiment of this specification, tomography polygon data can also projected to layer position rectangular coordinate system
After lower, the layer position data of closed fault polygonal internal are deleted.For example, the method for graphics, quick place to go closing can be passed through
The layer position data of tomography polygonal internal.So as to remove the layer site data for being unsatisfactory for seed point requirement in advance, after reduction
The continuous operand for determining seed point, improves operation efficiency.
S108: the value that interpolation calculation obtains interpolation layer site is carried out according to the seed point.
It can use the seed point interpolation calculation filtered out and obtain the layer position data in the interpolation layer site.For example, can be with
Interpolation calculation is carried out by such as inverse distance-weighting, least square scheduling algorithm, obtains the value in corresponding interpolation layer site.
Repeat the above steps the value in all layers of site in available entire target work area, to generate entire target work area
Layer position.
Using this specification scheme provided by the above embodiment, by carrying out gridding processing to layer position data first, so
Afterwards, it obtains and apart from interpolation layer site apart from nearest layer position grid vertex.It utilizes using this layer of position grid vertex as vertex
The splicing of four adjacent mesh determines prepared seed point, can substantially increase the efficiency for searching for preparation seed point.Then, into
One step increases faults control to determine seed point, to make full use of useful data in interpolation calculation, effectively prevent tomography
Influence of another side data to calculated result, the value in the determination interpolation layer site of efficiently and accurately.
Fig. 2 indicates another layer position interpolation method embodiment flow diagram that this specification provides.As shown in Fig. 2, this
Specification embodiment also provides a kind of layer of position interpolation method, the method may include:
S202: gridding processing is carried out to layer position data according to default layer position sizing grid；
S204: the rectangular extent of a layer position sizing grid is determined centered on the grid vertex of layer position, obtains the rectangle
Interpolation layer site in range obtains value cell to be inserted；
S206: the pre- of the value cell to be inserted is searched in the adjacent layer position grid centered on the grid vertex of the layer position
Standby seed point, obtains third preparation seed point set；
S208: the value in interpolation layer site in value cell to be inserted is calculated according to the third preparation seed point set, comprising:
Judge the interpolation layer site in the preparation seed point and value cell to be inserted of the third preparation seed point concentration
Line and tomography polygon whether there is intersection point, if it is judged that be it is no, by corresponding preparation seed point be determined as correspondence to
The seed point in interpolated layer site carries out the value that interpolation calculation obtains corresponding interpolation layer site according to the seed point.
In the present embodiment, to the method that layer position data carry out gridding processing, side provided by the above embodiment can be referred to
Case is implemented, and does not do tired state here.It is then possible to determine a collection of interpolation using four adjacent layer position grids for possessing common vertex
Layer site, the range in this batch of interpolation layer site are centered on this common vertex, the rectangle model of layer position sizing grid
It encloses.It can be using the interpolation layer site in the rectangular extent as a value cell to be inserted.
The preparation seed point of the value cell to be inserted is searched in the adjacent layer position grid centered on the vertex, is obtained
The third preparation seed point set of the value cell to be inserted.Accordingly, the method for searching for preparation seed point can also refer to above-mentioned implementation
The scheme that example provides is implemented, and does not do tired state here.
The preparation seed point that all interpolation layers site in value cell to be inserted can be concentrated with third preparation seed point
As the preparation seed point of itself, so as to once determine pre- the making preparations for sowing in all interpolation layers site in this rectangular extent
Sub- point.To reduce prepared seed point volumes of searches, improve the efficiency of prepared seed point search.
Then, the value in interpolation layer site in value cell to be inserted is calculated according to the third preparation seed point set.Such as
Certain value cell ∑ A to be insertedi, ∑ AiComprising m interpolation layer site, correspond to ∑ AiThird preparation seed point set be ∑ Bi.It is right
J-th of interpolation layer site A in m interpolation layer sitej, judge j-th of interpolation layer site and ∑ AiCorresponding ∑ Bi
In each prepared seed point line whether with tomography polygon there are intersection point, and then determine the kind in j-th of interpolation layer site
It is sub-, obtain the seed point set in j-th of interpolation layer site.It is to be inserted according to j-th of acquisition of seed point set progress interpolation calculation
It is worth the value in layer site.According to the above method, value cell ∑ A to be inserted is determinediThe value in middle m interpolation layer site.
In one embodiment of this specification, multithreads computing can also be utilized, each value cell to be inserted is calculated separately
The value in interior interpolation layer site.
Specifically, the maximum number N of a parallel computation can be calculated according to current computer CPU number.
It is then possible to obtain N number of value cell to be inserted, N number of value cell to be inserted and corresponding prepared seed point set are divided respectively
It is fitted in the calculating equipment of N number of parallel computation, each interpolation layer site in corresponding interpolating unit is determined in each calculating equipment
Seed point, and obtain according to seed point interpolation calculation the value in corresponding interpolation layer site.
Multi-threaded parallel operation is carried out to N number of value cell to be inserted and calculates separately out N number of interpolation list according to the method described above
The value in each interpolation layer site in member.
Wherein it is determined that the seed point in interpolation layer site and according to seed point interpolation calculation obtain corresponding interpolation layer position
The scheme of the value of point can be implemented with reference to scheme provided by the above embodiment, not do tired state here.
It is searched for by that can be greatly reduced using value cell to be inserted the time of preparation seed point, meanwhile, by interpolation list
Member can also further increase the efficiency of data processing as a unit to carry out concurrent operation.And it can also make full use of
Computer resource, all cpu resources for the current computer that avoids perfecting monopoly, influences the reaction speed of other application.It can be with
Reduction needs the requirement of the big work area, complicated fault polygon situation of larger space to Cyberspace and arithmetic speed, further
Improve the efficiency of overall data process.
In order to enable further verifying the practicability and feasibility of this method, this specification additionally provided using above-mentioned side
One specific example of case, by taking certain work area as an example, it is known that the work area situation is as follows:
(1) gridding processing is carried out to the layer position data in target work area.
Layer position data can be projected in three-dimensional cartesian coordinate system, reflect the three of target work area by discrete layer site
Layer position data distribution is tieed up, the distance between layer site is known as layer bit interval.It is then possible to be slightly larger than layer position interpretation interval half
As the size of layer position grid, gridding processing is carried out to layer position data, layer position grid Hg is obtained, is wrapped in each layer of position grid Hg
Containing multiple layers of site.
(2) gridding processing is carried out to tomography.
Before carrying out gridding processing to tomography, all tomography polygons are numbered, and pass through coordinate transform,
It projects under the rectangular coordinate system of layer position data.
Then, the layer position data of closed fault polygonal internal are quickly removed by the methods of graphics, it is subsequent to reduce
Operand, improve operation efficiency.
Tomography is subjected to gridding processing, the size of tomography grid is identical with layer bit interval, facilitates subsequent tomography grid
Merge, and records the number of the tomography in each tomography grid.
All tomography polygons are recorded simultaneously intersects corresponding layer site ∑ C with Inline, Crossline survey line, it can
With quickly filtered out in subsequent determining seed point a batch preparation seed point, with optimize it is subsequent ask tomography polygon and preparation seed
The efficiency whether point and interpolation layer site line intersect.
(3) preparation seed point is searched for.
Then interpolation layer site A is obtained, searching for around it has the layer site of value as preparation seed point.It can be with when search
It obtains with interpolation layer site apart from nearest layer position grid vertex D, in search by using D as four adjacent layers on vertex
Preparation seed point B is searched in the grid of position, obtains preparation seed point set ∑ B1。
Obtain preparation seed point set ∑ B1In be overlapped with ∑ C layer site (or calculate ∑ B1With the intersection of ∑ C) ∑ B1C,
Then, ∑ B is removed1Included in ∑ B1Layer site in C obtains preparation seed point set ∑ B2.Break to quickly remove to be located at
Preparation seed point on layer polygon, improves the efficiency of subsequent determining seed point.
(4) tomography, which constrains, determines seed point.
Required seed point is finally determined by tomography constraint, obtains current interpolation layer site A and ∑ B2In preparation
The line segment AB of seed point B composition, judges whether line segment AB and each tomography polygon Fn have intersection point.Specifically, available AB institute
Then the grid intersection ∑ fm (m=A, B) in tomography grid Fg passed through calculates each tomography grid in grid intersection ∑ fm and is wrapped
The union of the tomography number contained, obtains tomography number union.To obtain the tomography for being possible to intersection according to tomography number union
Polygon.
It is numbered again with AB and tomography and each tomography is concentrated to number corresponding tomography polygon intersection, if all not handing over
Point, then B is seed point required for A, and otherwise B is not seed point.
To ∑ B2In all prepared seed points repeat above-mentioned judgment step, obtain all seed point T of Ab。
(5) interpolation calculation.
According to the resulting all seed point T of step (4)b, obtained using inverse distance-weighting, least square scheduling algorithm interpolation calculation
To the value of A.
(6) step (3)-(5) are repeated, all layers of place value in entire work area are calculated.
As shown in Figure 3 and 4, Fig. 3 and Fig. 4 indicates that the scheme that this specification embodiment provides is inserted with prior art layer position
It is worth the contrast schematic diagram of result.Wherein, line near Faults interpolation result has been explained in Fig. 3 expression, and the dotted line in Fig. 3 is indicated using existing
The layer bit line for having technical method interpolation to obtain, solid line indicate the layer bit line obtained using this specification scheme interpolation.Fig. 4 is indicated not
Explain line near Faults interpolation result, the line of the arrow meaning in Fig. 4 is the layer bit line obtained using this specification scheme interpolation,
Another line is the layer bit line obtained using art methods interpolation.By comparison, it was found that being mentioned using this specification embodiment
The scheme of confession carries out layer position interpolation, and in tomography attachment layer position, the result of interpolation extends according to the trend of layer, do not occur going here and there layer and
The phenomenon that shake, meets geologic rule.
This specification additionally provides another specific example using above scheme, can be further improved entire work area layer
The efficiency of position interpolation.
Firstly, accordingly and tomography carries out gridding processing respectively, the same above-mentioned steps of concrete scheme (1) and (2) to layer digit.
Then, a collection of interpolation layer site is determined using four adjacent layer position grids for possessing common vertex D, this batch to be inserted
The range in value layer site is centered on this common vertex D, the rectangular extent of layer position sizing grid.By the rectangle model
Interpolation layer site in enclosing is as a value cell ∑ A to be inserted.
Pre- the making preparations for sowing of the value cell ∑ A to be inserted is searched in adjacent four layer positions grid centered on the vertex D
Sub- point B obtains preparation seed point set ∑ B1.And filtering ∑ B1Layer site in C obtains preparation seed point set ∑ B2。
Later, the maximum number N of a parallel computation can be calculated according to current computer CPU number.Obtain N number of interpolation
Unit ∑ A, by N number of value cell ∑ A to be inserted and corresponding prepared seed point set ∑ B2It is separately dispensed into the calculating of N number of parallel computation
In equipment.The seed point in each interpolation layer site in ∑ A is determined in each calculating equipment, and according to seed point interpolation calculation
Obtain the value in corresponding interpolation layer site.
Wherein it is determined that the seed point in interpolation layer site and according to seed point interpolation calculation obtain corresponding interpolation layer position
The scheme of the value of point can be implemented with reference to step (4) and (5).
Actual operation shows to implement the scheme provided using this specification, can be with further by the way of concurrent operation
Significantly promote operation efficiency.For the less complicated or less big range work area (within thousand square kilometres) of tomography, substantially
It is that moment operation is completed；For big work area (up to ten thousand square kilometres), and the situation that tomography is more and complicated, it can also complete faster
(test had the work area of thousands of tomographies to do layer position interpolation at more than 14000 square kilometres, it is only necessary to less than 4 minutes for operation
Operation is completed, and interpolation is also fine).
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Specifically it is referred to
The description of aforementioned relevant treatment related embodiment, does not do repeat one by one herein.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims
It is interior.In some cases, the movement recorded in detail in the claims or step can be come according to the sequence being different from embodiment
It executes and desired result still may be implemented.In addition, process depicted in the drawing not necessarily require show it is specific suitable
Sequence or consecutive order are just able to achieve desired result.In some embodiments, multitasking and parallel processing be also can
With or may be advantageous.
A kind of layer of position interpolation method that this specification one or more embodiment provides, can be by first to layer position data
Gridding processing is carried out, prepared seed point is determined by the splicing of four grid nearest apart from interpolation layer site distance,
Substantially increase the efficiency for searching for preparation seed point.Then, faults control is further increased to determine seed point, to mention significantly
The high accuracy of interpolation calculation result.Using the scheme of this specification embodiment, can be improved the accuracy of layer position interpolation with
And efficiency.
Based on layer position interpolation method described above, this specification one or more embodiment also provides a kind of layer of position interpolation
Device.The device may include the system for having used this specification embodiment the method, software (application), module, group
Part, server etc. simultaneously combine the necessary device for implementing hardware.Based on same innovation thinking, this specification embodiment provide one
Device in a or multiple embodiments is as described in the following examples.The implementation and method phase solved the problems, such as due to device
Seemingly, therefore the implementation of the specific device of this specification embodiment may refer to the implementation of preceding method, and overlaps will not be repeated.
Used below, the combination of the software and/or hardware of predetermined function may be implemented in term " unit " or " module ".Although with
Device described in lower embodiment is preferably realized with software, but the combined realization of hardware or software and hardware
It may and be contemplated.Specifically, Fig. 5 indicates the modular structure for a kind of layer of position interpolating apparatus embodiment that this specification provides
Schematic diagram, such as Fig. 5, the apparatus may include:
First layer position gridding module 302 can be used for carrying out grid to layer position data according to default layer position sizing grid
Change processing；
First search module 304 can be used for obtaining the layer position for meeting preset threshold condition with interpolation layer site distance
Grid vertex searches for the preparation seed point in interpolation layer site in the adjacent layer position grid centered on the vertex,
Obtain the first preparation seed point set；
Seed point determining module 306, can be used for judging preparation seed point that the first preparation seed point is concentrated and to
The line in interpolated layer site and tomography polygon whether there is intersection point, if it is judged that be it is no, seed point will be prepared accordingly
It is determined as the seed point in interpolation layer site；
Interpolation calculation module 308 can be used for carrying out interpolation calculation according to the seed point obtaining interpolation layer position
The value of point.
Using the offer scheme of this specification above-described embodiment, the accuracy and efficiency of layer position interpolation can be improved.
Fig. 6 indicates the modular structure schematic diagram for another layer position interpolating apparatus embodiment that this specification provides, such as Fig. 6 institute
Show, the apparatus may include:
Second layer position gridding module 402 can be used for carrying out grid to layer position data according to default layer position sizing grid
Change processing；
Interpolating unit division module 404 can be used for determining a layer position sizing grid centered on the grid vertex of layer position
Rectangular extent, obtain interpolation layer site in the rectangular extent and obtain value cell to be inserted；
Second search module 406, can be used for searching in the adjacent layer position grid centered on the vertex it is described to
The preparation seed point of interpolating unit obtains third preparation seed point set；
Concurrent operation module 408 can be used for calculating separately to be inserted in each value cell to be inserted using multithreads computing
It is worth the value in layer site.
Wherein, the concurrent operation module 408 may include multiple parallel seed point determination unit and interpolation calculation
Unit, wherein
The seed point determination unit, can be used for judging preparation seed point that the third preparation seed point is concentrated and to
The line in the interpolation layer site in interpolating unit and tomography polygon whether there is intersection point, if it is judged that be it is no, by phase
The preparation seed point answered is determined as the seed point in corresponding interpolation layer site；
The interpolation calculation unit can be used for carrying out interpolation calculation according to the seed point and obtain corresponding interpolation layer site
Value.
Using the offer scheme of this specification above-described embodiment, the accuracy and effect of layer position interpolation can be further improved
Rate.
It should be noted that device described above can also include other embodiment party according to the description of embodiment of the method
Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
A kind of layer of position interpolating apparatus that this specification one or more embodiment provides, can be by first to layer position data
Gridding processing is carried out, prepared seed point is determined by the splicing of four grid nearest apart from interpolation layer site distance,
Substantially increase the efficiency for searching for preparation seed point.Then, faults control is further increased to determine seed point, to mention significantly
The high accuracy of interpolation calculation result.Using the scheme of this specification embodiment, can be improved the accuracy of layer position interpolation with
And efficiency.
Method or apparatus described in above-described embodiment that this specification provides can realize that business is patrolled by computer program
It collects and records on a storage medium, the storage medium can be read and be executed with computer, realize this specification embodiment institute
The effect of description scheme.Therefore, this specification also provides a kind of layer of position interpolating apparatus, including processor and storage processor can be held
Row instruction memory, when described instruction is executed by the processor realization the following steps are included:
Gridding processing is carried out to layer position data according to default layer position sizing grid；
Obtain the layer position grid vertex for meeting preset threshold condition with interpolation layer site distance；
The preparation seed point that interpolation layer site is searched in the adjacent layer position grid centered on the vertex, is obtained
Obtain the first preparation seed point set；
Judge that the line and tomography of the preparation seed point and interpolation layer site that the first preparation seed point is concentrated are polygon
Shape whether there is intersection point, if it is judged that be it is no, corresponding seed point of preparing is determined as the kind in interpolation layer site
Sub- point；
The value that interpolation calculation obtains interpolation layer site is carried out according to the seed point.
The storage medium may include the physical unit for storing information, usually by after information digitalization again with benefit
The media of the modes such as electricity consumption, magnetic or optics are stored.It may include: that letter is stored in the way of electric energy that the storage medium, which has,
The device of breath such as, various memory, such as RAM, ROM；The device of information is stored in the way of magnetic energy such as, hard disk, floppy disk, magnetic
Band, core memory, magnetic bubble memory, USB flash disk；Using optical mode storage information device such as, CD or DVD.Certainly, there are also it
Readable storage medium storing program for executing of his mode, such as quantum memory, graphene memory etc..
It should be noted that device described above can also include other embodiment party according to the description of embodiment of the method
Formula.Concrete implementation mode is referred to the description of related method embodiment, does not repeat one by one herein.
A kind of layer of position interpolating apparatus described in above-described embodiment, can be by first carrying out at gridding layer position data
Reason determines prepared seed point by the splicing of four grid nearest apart from interpolation layer site distance, substantially increases and search
The efficiency of rope preparation seed point.Then, faults control is further increased to determine seed point, to substantially increase interpolation calculation
As a result accuracy.Using the scheme of this specification embodiment, the accuracy and efficiency of layer position interpolation can be improved.
This specification also provides a kind of layer of position interplotation system, and the system can be individual layer position interplotation system, can also
To apply in a plurality of types of seismic data analysis systems.The system can be individual computer, also may include
One or more the methods of this specification or the server cluster of one or more embodiment devices, system (packet are used
Include distributed system), software (application), practical operation device, logic gates device, quantum computer etc. and combine necessary
Implement the terminal installation of hardware.Layer position interplotation system may include that at least one processor and storage computer are executable
The memory of instruction, the processor realize method described in above-mentioned any one or multiple embodiments when executing described instruction
The step of.
It should be noted that system described above can also include others according to the description of method or Installation practice
Embodiment, concrete implementation mode are referred to the description of related method embodiment, do not repeat one by one herein.
A kind of layer of position interplotation system described in above-described embodiment, can be by first carrying out at gridding layer position data
Reason determines prepared seed point by the splicing of four grid nearest apart from interpolation layer site distance, substantially increases and search
The efficiency of rope preparation seed point.Then, faults control is further increased to determine seed point, to substantially increase interpolation calculation
As a result accuracy.Using the scheme of this specification embodiment, the accuracy and efficiency of layer position interpolation can be improved.
It should be noted that this specification device or system described above according to the description of related method embodiment also
It may include other embodiments, concrete implementation mode is referred to the description of embodiment of the method, does not go to live in the household of one's in-laws on getting married one by one herein
It states.All the embodiments in this specification are described in a progressive manner, and same and similar part is mutual between each embodiment
Mutually referring to each embodiment focuses on the differences from other embodiments.Especially for hardware+program
For class, storage medium+program embodiment, since it is substantially similar to the method embodiment, so be described relatively simple, it is related
Place illustrates referring to the part of embodiment of the method.
It is above-mentioned that this specification specific embodiment is described.Other embodiments are in the scope of the appended claims
It is interior.In some cases, the movement recorded in detail in the claims or step can be come according to the sequence being different from embodiment
It executes and desired result still may be implemented.In addition, process depicted in the drawing not necessarily require show it is specific suitable
Sequence or consecutive order are just able to achieve desired result.In some embodiments, multitasking and parallel processing be also can
With or may be advantageous.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity,
Or it is realized by the product with certain function.It is a kind of typically to realize that equipment is computer.Specifically, computer for example may be used
Think personal computer, laptop computer, vehicle-mounted human-computer interaction device, cellular phone, camera phone, smart phone, individual
Digital assistants, media player, navigation equipment, electronic mail equipment, game console, tablet computer, wearable device or
The combination of any equipment in these equipment of person.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module can be realized in the same or multiple software and or hardware when specification one or more, it can also be with
The module for realizing same function is realized by the combination of multiple submodule or subelement etc..Installation practice described above is only
It is only illustrative, for example, in addition the division of the unit, only a kind of logical function partition can have in actual implementation
Division mode, such as multiple units or components can be combined or can be integrated into another system or some features can be with
Ignore, or does not execute.Another point, shown or discussed mutual coupling, direct-coupling or communication connection can be logical
Some interfaces are crossed, the indirect coupling or communication connection of device or unit can be electrical property, mechanical or other forms.
It is also known in the art that other than realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come so that controller with logic gate, switch, specific integrated circuit, programmable
Logic controller realizes identical function with the form for being embedded in microcontroller etc..Therefore this controller is considered one kind
Hardware component, and the structure that the device for realizing various functions that its inside includes can also be considered as in hardware component.Or
Person even, can will be considered as realizing the device of various functions either the software module of implementation method can be hardware again
Structure in component.
The present invention be referring to according to the method for the embodiment of the present invention, the process of equipment (system) and computer program product
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap
Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including described want
There is also other identical elements in the process, method or equipment of element.
It will be understood by those skilled in the art that this specification one or more embodiment can provide as method, system or calculating
Machine program product.Therefore, this specification one or more embodiment can be used complete hardware embodiment, complete software embodiment or
The form of embodiment combining software and hardware aspects.Moreover, this specification one or more embodiment can be used at one or
It is multiple wherein include computer usable program code computer-usable storage medium (including but not limited to magnetic disk storage,
CD-ROM, optical memory etc.) on the form of computer program product implemented.
This specification one or more embodiment can computer executable instructions it is general on
It hereinafter describes, such as program module.Generally, program module includes executing particular task or realization particular abstract data type
Routine, programs, objects, component, data structure etc..This this specification one can also be practiced in a distributed computing environment
Or multiple embodiments, in these distributed computing environments, by being held by the connected remote processing devices of communication network
Row task.In a distributed computing environment, program module can be located at the local and remote computer including storage equipment
In storage medium.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ",
The description of " specific example " or " some examples " etc. means specific features described in conjunction with this embodiment or example, structure, material
Or feature is contained at least one embodiment or example of this specification.In the present specification, to the signal of above-mentioned term
Property statement must not necessarily be directed to identical embodiment or example.Moreover, specific features, structure, material or the spy of description
Point may be combined in any suitable manner in any one or more of the embodiments or examples.In addition, without conflicting with each other,
Those skilled in the art can be by different embodiments or examples described in this specification and different embodiments or examples
Feature is combined.
The foregoing is merely the embodiments of this specification, are not limited to this specification.For art technology
For personnel, this specification can have various modifications and variations.It is all made any within the spirit and principle of this specification
Modification, equivalent replacement, improvement etc., should be included within the scope of the claims of this specification.
Claims (12)
1. a kind of layer of position interpolation method characterized by comprising
Gridding processing is carried out to layer position data according to default layer position sizing grid；
Obtain the layer position grid vertex for meeting preset threshold condition with interpolation layer site distance；
It searches for the preparation seed point in interpolation layer site in the adjacent layer position grid centered on the vertex, obtains the
One preparation seed point set；
Judge the line for the preparation seed point and interpolation layer site that the first preparation seed point is concentrated is with tomography polygon
It is no that there are intersection points, if it is judged that be it is no, corresponding preparation seed point is determined as the seed point in interpolation layer site；
The value that interpolation calculation obtains interpolation layer site is carried out according to the seed point.
2. layer position according to claim 1 interpolation method, which is characterized in that judgement the first preparation seed point set
In preparation seed point and interpolation layer site line and tomography polygon whether there is intersection point, comprising:
Gridding processing is carried out to layer data according to layer bit interval size, obtains the preparation that the first preparation seed point is concentrated
The tomography grid that the line in seed point and interpolation layer site passes through；
Tomography number in the tomography grid is merged and obtains tomography number union；
Judge that the tomography number that the line in the prepared seed point and interpolation layer site is numbered and concentrated with tomography is corresponding disconnected
Layer polygon whether there is intersection point.
3. layer position according to claim 2 interpolation method, which is characterized in that described to obtain the first preparation seed point set
In preparation seed point and interpolation layer site the tomography grid that passes through of line, comprising:
The layer site that record seismic survey lines intersect with tomography polygon；
The layer site that seismic survey lines described in the first preparation seed point collection intersect with tomography polygon is removed, it is pre- to obtain second
Make preparations for sowing sub- point set；
Obtain the tomography grid that the line of preparation seed point and interpolation layer site that the second preparation seed point is concentrated passes through.
4. layer position according to claim 1 interpolation method, which is characterized in that the basis presets layer position sizing grid to layer
Position data carry out gridding processing, comprising:
Gridding processing is carried out to layer position data according to the layer position interpretation interval size of half.
5. layer position according to claim 1 interpolation method, which is characterized in that interpolation layer site is pre- described in described search
Before standby seed point, further includes: delete the layer position data in closed fault polygon.
6. layer position according to claim 1 interpolation method, which is characterized in that interpolation layer site is pre- described in described search
Before standby seed point, further includes:
Interpolation calculation is carried out to the layer position grid vertex of no value, obtains the value of layer position grid vertex.
7. a kind of layer of position interpolation method characterized by comprising
Gridding processing is carried out to layer position data according to default layer position sizing grid；
Centered on the grid vertex of layer position determine a layer position sizing grid rectangular extent, obtain in the rectangular extent to
Interpolated layer site obtains value cell to be inserted；
The preparation seed point of the value cell to be inserted is searched in the adjacent layer position grid centered on the grid vertex of the layer position,
Obtain third preparation seed point set；
The value in interpolation layer site in value cell to be inserted is calculated according to the third preparation seed point set, comprising:
Judge the line in the interpolation layer site in the preparation seed point and value cell to be inserted of the third preparation seed point concentration
With tomography polygon whether there is intersection point, if it is judged that be it is no, by it is corresponding prepare seed point be determined as correspond to interpolation
The seed point in layer site carries out the value that interpolation calculation obtains corresponding interpolation layer site according to the seed point.
8. layer position according to claim 7 interpolation method, which is characterized in that described to calculate interpolation layer in value cell to be inserted
The value in site, comprising:
Using multithreads computing, the value in interpolation layer site in each value cell to be inserted is calculated separately.
9. a kind of layer of position interpolating apparatus characterized by comprising
First layer position gridding module, for carrying out gridding processing to layer position data according to default layer position sizing grid；
First search module, for obtaining the layer position grid vertex for meeting preset threshold condition with interpolation layer site distance,
The preparation seed point in interpolation layer site is searched in adjacent layer position grid centered on the vertex, obtains the first preparation
Seed point set；
Seed point determining module, for judging preparation seed point that the first preparation seed point is concentrated and interpolation layer site
Line and tomography polygon whether there is intersection point, if it is judged that be it is no, by it is corresponding prepare seed point be determined as described in
The seed point in interpolated layer site
Interpolation calculation module, for carrying out the value that interpolation calculation obtains interpolation layer site according to the seed point.
10. a kind of layer of position interpolating apparatus characterized by comprising
Second layer position gridding module, for carrying out gridding processing to layer position data according to default layer position sizing grid；
Interpolating unit division module, for determining the rectangular extent of a layer position sizing grid centered on the grid vertex of layer position,
The interpolation layer site obtained in the rectangular extent obtains value cell to be inserted；
Second search module, for searching for the pre- of the value cell to be inserted in the adjacent layer position grid centered on the vertex
Standby seed point, obtains third preparation seed point set；
Concurrent operation module calculates separately interpolation layer site in each value cell to be inserted for utilizing multithreads computing
Value, and seed point determination unit and interpolation calculation unit including multiple concurrent operations, wherein
The seed point determination unit, for judging preparation seed point and value cell to be inserted that the third preparation seed point is concentrated
In interpolation layer site line and tomography polygon whether there is intersection point, if it is judged that be it is no, will prepare accordingly
Seed point is determined as the seed point in corresponding interpolation layer site；
The interpolation calculation unit, for carrying out the value that interpolation calculation obtains corresponding interpolation layer site according to the seed point.
11. a kind of layer of position interpolating apparatus, which is characterized in that including processor and for the storage of storage processor executable instruction
Device, when described instruction is executed by the processor realize the following steps are included:
Gridding processing is carried out to layer position data according to default layer position sizing grid；
Obtain the layer position grid vertex for meeting preset threshold condition with interpolation layer site distance；
It searches for the preparation seed point in interpolation layer site in the adjacent layer position grid centered on the vertex, obtains the
One preparation seed point set；
Judge the line for the preparation seed point and interpolation layer site that the first preparation seed point is concentrated is with tomography polygon
It is no that there are intersection points, if it is judged that be it is no, corresponding preparation seed point is determined as the seed point in interpolation layer site；
The value that interpolation calculation obtains interpolation layer site is carried out according to the seed point.
12. a kind of layer of position interplotation system, which is characterized in that including at least one processor and storage computer executable instructions
Memory, the processor realizes any one of claim 1-8 the method when executing described instruction the step of.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810796224.6A CN109143334B (en) | 2018-07-19 | 2018-07-19 | Level interpolation method, device and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810796224.6A CN109143334B (en) | 2018-07-19 | 2018-07-19 | Level interpolation method, device and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109143334A true CN109143334A (en) | 2019-01-04 |
| CN109143334B CN109143334B (en) | 2020-08-07 |
Family
ID=64800973
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810796224.6A Active CN109143334B (en) | 2018-07-19 | 2018-07-19 | Level interpolation method, device and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109143334B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222365A (en) * | 2011-07-29 | 2011-10-19 | 电子科技大学 | Method for reconstructing curved surface of complex space |
| CN102385067A (en) * | 2011-08-03 | 2012-03-21 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Drawing method for isoline containing reverse fault |
| CN102819039A (en) * | 2012-08-22 | 2012-12-12 | 电子科技大学 | Interval velocity model building method under complicated geological conditions |
-
2018
- 2018-07-19 CN CN201810796224.6A patent/CN109143334B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102222365A (en) * | 2011-07-29 | 2011-10-19 | 电子科技大学 | Method for reconstructing curved surface of complex space |
| CN102385067A (en) * | 2011-08-03 | 2012-03-21 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Drawing method for isoline containing reverse fault |
| CN102819039A (en) * | 2012-08-22 | 2012-12-12 | 电子科技大学 | Interval velocity model building method under complicated geological conditions |
Non-Patent Citations (4)
| Title |
|---|
| SUYUN LI ET AL.: "Geologic Surface Reconstruction Based on Fault Constraints", 《2007 IEEE》 * |
| 孙涛等: "GEOEAST高效工业制图技术", 《石油工业计算机应用》 * |
| 聂小燕等: "基于边界约束的空间复杂曲面重构算法", 《计算机应用研究》 * |
| 蒋晓路等: "基于曲面约束的地质层面重构方法研究", 《计算机应用研究》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109143334B (en) | 2020-08-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109255829B (en) | Method for generating hierarchical acceleration structure and performing intersection test and ray tracing system | |
| US11562526B2 (en) | 3-D graphics rendering with implicit geometry | |
| CN104699601B (en) | For the method and system for the injection for performing state machine driving | |
| Waser et al. | World lines | |
| CN107976713A (en) | A kind of lower method and device for removing sedimentation setting of higher-dimension seismic data input | |
| CN103177417B (en) | Based on the mathematical morphology LiDAR point cloud fast filtering method of GPGPU | |
| CN108286976A (en) | The fusion method and device and hybrid navigation system of a kind of point cloud data | |
| CN109100795A (en) | A kind of inspection point distribution method, the apparatus and system of face element | |
| KR102325047B1 (en) | Grahp data processing methdo and apparatus thereof | |
| Audet et al. | Robust and efficient polygon overlay on parallel stream processors | |
| CN110378175A (en) | The recognition methods of road edge and device | |
| CN106709503A (en) | Large spatial data clustering algorithm K-DBSCAN based on density | |
| CN109033127A (en) | A kind of synchrodata method of calibration, device and equipment | |
| CN108267783B (en) | A kind of method, apparatus and system of determining buried-hill trap | |
| CN109725351A (en) | A kind of the determination method, apparatus and system of 3D elastic wave mixed absorbing boundary | |
| CN112395745A (en) | Underground river reservoir body geological model establishing method and processing equipment | |
| CN108897036A (en) | A kind of seismic data processing technique and device | |
| CN111413741B (en) | Sandstone-type uranium ore resource amount calculation method and device | |
| CN108107471A (en) | The acquisition methods and device of a kind of point of orientation first arrival data volume | |
| CN109783381A (en) | A kind of test data generating method, apparatus and system | |
| CN109388306A (en) | Information display method and device | |
| CN109189343A (en) | A kind of metadata rule method, apparatus, equipment and computer readable storage medium | |
| CN106484532B (en) | GPGPU parallel calculating method towards SPH fluid simulation | |
| CN108828669A (en) | A kind of two-dimensional intersection survey line static corrections processing method, apparatus and system | |
| CN109143334A (en) | A kind of layer of position interpolation method, apparatus and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |