CN102867330A - Region-division-based spatial complex horizon reconstruction method - Google Patents
Region-division-based spatial complex horizon reconstruction method Download PDFInfo
- Publication number
- CN102867330A CN102867330A CN2012103115671A CN201210311567A CN102867330A CN 102867330 A CN102867330 A CN 102867330A CN 2012103115671 A CN2012103115671 A CN 2012103115671A CN 201210311567 A CN201210311567 A CN 201210311567A CN 102867330 A CN102867330 A CN 102867330A
- Authority
- CN
- China
- Prior art keywords
- interpolation
- tomography
- point
- macrolattice
- subregion
- 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
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The invention discloses a region-division-based spatial complex horizon reconstruction method. The horizon reconstruction in three-dimensional geological modeling under the complex topography conditions such as multi-fault overlapping can be realized without the needing of a fault polygon; simultaneously, the dynamic control of specified grid spacing is realized by a hybrid grid interpolation method, so that a geologic structure can be expressed finely; the reconstructed geologic horizon surface and the fault are closed strictly; and the processing, such as blocking processing, in later geological modeling is facilitated.
Description
Technical field
The invention belongs to the Review of geologic model building techniques field, particularly a kind of spatial complex layer position reconstructing method of dividing based on the zone.
Background technology
Geologic modeling is by certain technological means limited geologic data to be processed on computers, reproduces geologic structure model.Essence is to change limited discrete space sample points certificate into continuous, visual geology face or geologic body, and wherein the reconstruct of aspect is crucial.The most important thing is the model structure of interpolation method and aspect in the method for aspect reconstruct.
Interpolation is to utilize the method for surface fitting, sets up a continuous function by discrete input sample point, obtains the functional value of optional position with the function of this reconstruction.And layer position reconstruct is exactly the method for utilizing the discrete reconstructed surface of limited three dimensions.
Spatial interpolation methods commonly used has following several at present:
(1) method is taken advantage of in inverse distance: it is a weighted mean method of interpolation that method is taken advantage of in inverse distance.It is maximum that it thinks that several nearest known points of interpolation point are treated the contribution of interpolation point, and its contribution is inversely proportional to distance.The degree parameter controlling weight coefficient how along with leave a net point distance and descend.The given higher weight of nearer data point, the given less weight of data point far away.As follows:
In the formula, z is estimated value, z
iI known point value, d
iBe the distance of i known point and interpolation point, p is weighting coefficient, and its value is higher, shows that the influence factor of distance is larger, and interpolation result just more has smooth effect.Have faster computing velocity apart from the inverse ratio weighting algorithm, but easily be subjected to the impact of data point cluster, often in the practical application situation that a kind of isolated point is higher than its ambient data point occurs, so need to limit according to the different situations inverse ratio weighted method of adjusting the distance.
(2) natural neighbor interpolation method: the natural neighbor interpolation method is the gridding new method that Surfer7.0 just has.Its ultimate principle is for one group of Tyson shape changeable more, when adding a new data point in data centralization, will revise these Thiessen polygon, and uses the weighted mean of adjoint point will determine the weight of interpolation.It is homogeneous and unconverted, and evenly spaced data are carried out interpolation of great use, but for the interpolation of inhomogeneous data defectiveness then, and inhomogeneous original layers bit data is very common.
(3) Ke Lijin (Kriging) algorithm: claim again space autocovariance optimal interpolation method, it is a kind of optimal interpolation method with the naming of South Africa Mineral Engineering teacher D.G.Krige, is a kind of geological statistics gridding method of great use.The method is based upon on variation function theory and the structure analysis basis, calculates weights by introducing take distance as the variation function of independent variable, can reflect the space structure characteristic of variable, random distribution nature that again can response variable.Say from the statistical significance, it is from correlation of variables and variability, in limited area the value of regionalized variable is carried out without a kind of method inclined to one side, optimal estimation; From the interpolation angle, it is that the data of space distribution are asked partially a kind of method of interpolation estimation of linear optimal, nothing.Because Kriging method has optimum efficiency in the interpolation of space layer position, therefore this programme, solves the layer position interpolation problem of complex-terrain take Kriging method as the basis.
Above interpolation method all is confined to solve the aspect match of two dimensional surface, can not solve three-dimensional aspect fitting problems.In the geologic structure of reality, the geological condition complex, the aspect match of common two dimensional surface can not solve the fitting problems of this type of geological stratum.So the three dimensions problem need to be converted into two-dimensional problems, the above interpolation method of recycling carries out the aspect match.
The structure of FEM layer model is exactly the organizational form of aspect data, and different organizational forms is set up for model different characteristics.Main deck structure has level line, regular grid (GRID), TIN (TIN) in the existing solution that proposes.Wherein regular grid and TIN use more widely FEM layer model structure:
(1) regular grid (GRID)
Regular grid is to come the presentation layer bit model with regularly arranged rectangular node.Data structure is simple, is easy to network forming, and memory data output is little, can also compress storage, and various analysis and calculations are very convenient effectively.But because mesh spacing is certain, exist the shortcoming that is difficult to determine suitable size of mesh opening expression true to nature for the topography and geomorphology of complexity.
(2) irregular grid (TIN)
The space aspects sampled point that gathers is connected into the whole zone of covering and mutual nonoverlapping many triangles by certain rule, just consisted of TIN.TIN can express structure and the thin section of landform relatively accurately, takes preferably aspect unique point, characteristic curve into account, and the expression complex-terrain is more accurate than rectangular node.But the storage of the data of TIN and complicated operation.The triangle form point of TIN and the topological relation more complicated on limit be difficult to TIN is carried out terrain analysis, and storage space are larger, generally are only applicable among a small circle large-scale high precision landform modeling.
Prior art related to the present invention comprises:
On the reconstruction of three-dimensional geological aspect, many implementations are arranged at present both at home and abroad.Existing scheme is according to whether needs tomography polygon is divided into two kinds: a kind of aspect reconfiguration scheme that is based on the constraint of tomography polygon, another kind is based on the aspect reconfiguration scheme that the zone is divided.
(1) the aspect reconfiguration scheme that retrains based on the tomography polygon
The tomography polygon is exactly the intersection of three-dimension disclocation face and layer plane.Tomography polygon on a plurality of interfaces of aspect reconstruct need of work analysis-by-synthesis, the form of three-dimensional space interbed position is inferred in polygonal constraint according to tomography.Many interpre(ta)tive system algorithms also are to finish on the basis of each cross-sectional image, but not directly reconstruct on three-dimensional raw data.Generally adopt the manual trace mode at present in the combination of tomography polygon, namely adopt manual breakpoint method of attachment, this mode not only error is larger, and impact is configured to the precision of figure, and efficient is extremely low.Wang Zhaohu, Liu Cai etc. propose a kind of tomography polygon Automatic Combined technology of explaining based on layer position, but this technology only can be in the situation that layer position be very continuous, the meticulousr automatic tracing of data goes out the tomography polygon, higher for data demand, inapplicable common situation.
(2) the aspect reconfiguration scheme of dividing based on the zone
Less based on the existing scheme of three-dimensional geological aspect reconstruct that the zone is divided, the at present domestic grid cutting algorithm that only has a kind of overlapping subregion of Cai Qiang and Yang Qin proposition
[6]This algorithm is introduced the concept of bridge limit and crossing loop wire, by the aspect reconstruct under the interlock subdivision algorithm realization overlapping region.This scheme Shortcomings part: one, this scheme is by asking intersection to divide overlapping subregion in three dimensions, in the three dimensions intersection find the solution more complicated and also speed is slower.Two, this scheme does not propose to solve multiple trap-up and layer position pinching, the method for the three-dimensional geological aspect reconstruct in particular cases such as the end, ground, earth's surface is gone out in layer position.
Existing scheme is divided into two classes according to the institutional framework of reconstruct aspect: based on the aspect reconstruct of specification grid (GRID) with based on the aspect reconstruct of triangular mesh (TIN).Regular grid is to come the presentation layer bit model with regularly arranged rectangular node.Because mesh spacing is certain, for the topography and geomorphology of complexity, be difficult to determine suitable size of mesh opening expression true to nature.The triangle form point of triangular mesh and the topological relation more complicated on limit be difficult to TIN is carried out terrain analysis, and storage space are larger, generally are only applicable among a small circle large-scale high precision landform modeling.
The outer scheme for the three-dimensional aspect reconstruct of geology and interpolation of Present Domestic has a lot, but all has following limitation:
(1) the three-dimensional aspect reconfiguration scheme of present most geology all needs to obtain in advance the changeable row of tomography.Huge by human-edited's tomography polygon workload, and accuracy is difficult to guarantee inefficiency.Also there is not good solution and automatically generate the tomography polygon.
(2) the three-dimensional aspect reconfiguration scheme of most of geology is not considered the situation of tomography, although some schemes have been considered tomography, does not consider for the situation of multiple trap-up.
(3) mesh spacing of the three-dimensional aspect reconfiguration scheme of geology of existing rule-based grid processing is determined, can not dynamically be suitable for the needs of geologic structure, is difficult to meticulous expression geologic structure.
Summary of the invention
In order to overcome the above-mentioned shortcoming of prior art, the invention provides a kind of spatial complex layer position reconstructing method of dividing based on the zone, by dividing overlapping subregion, under not needing the polygonal situation of tomography, can solve the aspect reconstruction in the three-dimensional geological modeling under the MODEL OVER COMPLEX TOPOGRAPHY such as multiple tomography is overlapping; By the interpolation method of hybrid grid, realized the dynamic control to the specification mesh spacing simultaneously, thereby can meticulously give expression to geologic structure, the geological stratification plane that reconstructs and tomography are strict sealings, are convenient to the processing in the later stage Geologic modeling, process as becoming piece.
The technical solution adopted for the present invention to solve the technical problems is: a kind of spatial complex layer position reconstructing method of dividing based on the zone comprises the steps:
Step 1, hybrid grid interpolation method: at the little grid of macrolattice nested inside, first macrolattice is carried out interpolation during interpolation, and then utilize the new interpolation point on existing macroreticular Seed Points and the macrolattice that little grid is carried out interpolation;
Step 2, definition layer bit field data structure;
Step 3, fault surface match: get the two-dimentional envelope of tomography kind subsegment decryption on grid plan as the interpolation envelope scope of cross section correct, then adopt the interpolation method of hybrid grid that regmatic network lattice point in the tomography envelope scope is carried out interpolation;
Step 5, initialization original segment data: the incidence relation of initialization segments and restrained boundary;
Step 6, subregion are divided;
Step 7, the match of tomography subregion layer plane:
(1) initialization upper lower burrs Seed Points;
(2) tomography upper lower burrs subregion face interpolation;
(3) the up and down intercepting of actual card of tomography;
The layer plane match of step 8, non-tomography subregion:
(1) initialization seed point: both with the data point of non-fault region as Seed Points, again with the borderline region data point of tomography subregion as Seed Points; And for the tomography subregion after the interpolation, with all on two dimensional surface not heavily data points of value also as Seed Points;
(2) interpolation of non-tomography subregion: for the Seed Points that has generated, adopt the hybrid grid interpolation method that each net point of the subregion scope of non-tomography is carried out interpolation;
The processing at the end of ground, earth's surface is gone out in step 9, layer position pinching, layer position: remove the illegal pinching part in current layer position and go out ground, the earth's surface part at the end.
Compared with prior art, good effect of the present invention is: solved the layer position surface reconstruction problem under the complex-terrain in the three-dimensional geological modeling, provide new thinking for isoplethes drawing, geology become piece etc.The inventive method has following advantage:
(1) only needs original layer position sampled data and tomography sampled data, do not need the tomography polygon can be to reconstructing aspect.The tomography polygon that need not to edit has manually saved numerous and diverse human-edited.
(2) support all kinds of tomographies, comprise trap-down, trap-up, can process the situation of multiple tomography simultaneously, be applicable to the topographic structure of various complexity, can process the layer position interpolation under the various complex situations, have good adaptability.
(3) the multiple application in suitable three-dimensional geological modeling and the two-dimensional layer iso-surface patch has good versatility.
(4) the layer position interpolation method that adopts hybrid network to format had both guaranteed the quick demonstration of data, had guaranteed again the precision of match.
(5) processing that ground, earth's surface end complex situations are gone out in supporting layer position pinching, layer position.
(6) propose to utilize the section envelope to carry out the method that subregion is divided, and provided concrete realization.
(7) by boundary constraint, realized the seamless match of aspect and tomography.Geological stratum and fault surface after the reconstruct strictly seal.
Description of drawings
Examples of the present invention will be described by way of reference to the accompanying drawings, wherein:
Fig. 1 is the process flow diagram of the inventive method;
Fig. 2 is the schematic diagram of deletion error data;
Fig. 3 is the schematic diagram of initialization original segment data.
Embodiment
First some basic geologic structures and scheme term are defined:
Layer position: refer to a certain ad-hoc location in succession of strata.
Tomography: after the stressed effect fracture of rock mass, the rift structure of remarkable displacement occurs in the both sides sillar along rupture surface.
Segment data: layer position or tomography on a section by one group of continuous data acquisition that discrete point forms.
Hanging wall: be positioned at the fault surface top.
Footwall: be positioned at the fault surface below.
Stratigraphic pitch-out: refer to sedimentary deposit towards the sedimentary basin edge, its thickness gradually attenuation until deposition not.
Gridding: the zone that discrete points data is carried out is in logic divided, with the logical nets of formation rule, be convenient to layer position interpolation.
Interpolation: utilize known point to calculate the process of unknown point.
Match: utilize layer position interpolation to finish the process that data afterwards form aspect.
A kind of spatial complex layer position reconstructing method of dividing based on the zone as shown in Figure 1, comprises the steps:
Step 1, hybrid grid interpolation method
In existing grid interpolation method commonly used, if the too sparse meeting of grid causes precision inadequate, grid is too close then can to affect interpolation speed.This programme proposes a kind of hybrid grid interpolation method of secondary Grid Nesting, and the one-level grid has guaranteed the speed of interpolation, and the secondary grid has guaranteed the precision of interpolation.
Hybrid grid is the secondary grid, i.e. the little grid of macrolattice nested inside.First the point on the macrolattice is carried out interpolation during interpolation, and then utilize the new interpolation point on existing macroreticular Seed Points and the macrolattice that little grid is carried out interpolation.Because the ratio of big or small grid is searched the scope of Seed Points and the quantity of Seed Points and can be inputted manually, makes interpolation carry out different parameter settings for different situations.So both solve seed that direct Kriging supports very few and cause the inadequate problem of interpolation precision, improved again interpolation speed.
(1) macrolattice interpolation
For the macrolattice point of interpolation, come this interpolation point is carried out interpolation calculation by the legal Seed Points around the search value interposition point.In the sparse layer position of some Seed Points, only relying on Seed Points to carry out interpolation can't all calculate all interpolation points complete.In this case, reservation can't be calculated the interpolation point of Z value, and after all interpolation points that can rely on Seed Points to carry out interpolation calculation were disposed, the point that interpolation is finished was used as new Seed Points and is carried out layer position compensation interpolation, until all interpolation points all are disposed.The concrete grammar of macrolattice interpolation is as follows:
1) scans one by one the macrolattice point, determine whether macrolattice point needs to carry out interpolation, if there has been data point in this macrolattice point, does not then need to carry out interpolation, otherwise need interpolation.
2) to needing the macrolattice point of interpolation, according to artificial inputted search scope, the macrolattice number of namely searching for is again searched for current interpolation point macroreticular Seed Points data on every side.
3) according to the sector of artificial input Seed Points is carried out the sector classification.
4) Seed Points to finding in each sector is arranged apart from current macrolattice point distance from small to large according to it.
5) according to each sector Seed Points quantity selected threshold of inputting, and select to join in the Seed Points formation of current macrolattice interpolation apart from the Seed Points data in the nearest threshold range of interpolation net point in this Sector Range.
6) if current macrolattice point does not find effective Seed Points data to carry out interpolation, current macrolattice interpolation not then, the next macrolattice point of diversion treatments, otherwise turn to the 7th) step.
7) with the Seed Points in the formation of current macrolattice interpolation Seed Points current macrolattice point is carried out the Kriging interpolation.
8) if all macrolattice points processed once, whether then check has the failure of macrolattice point interpolation, if the failure of macrolattice point interpolation is arranged, then with the macrolattice point of interpolation success also as Seed Points, return the 1st) step is carried out again interpolation, until all net point interpolation successes.
(2) little net point interpolation
After the macrolattice interpolation, all macrolattice points have been plugged corresponding data on the grid of layer position.Recycle existing macrolattice point data and carry out the subdivision interpolation of little grid.In order to improve interpolation speed, the interpolation of little net point is processed take macrolattice as unit.For certain macrolattice, search for the interior Seed Points of current macrolattice and the Seed Points data in eight macrolattices on every side as the primordial seed point, then utilize the Kriging interpolation that these primordial seed points are fitted to a curved surface.
Can directly carry out interpolation by the value of getting on this curved surface for all little net points in the macrolattice, avoid repeatedly Kriging interpolation.Namely by a piece of data, a Kriging interpolation fitting is realized the interpolation to little net point in the whole macrolattice, greatly improves the efficient of interpolation.For certain interpolation point I
xIf the legal seed point set that finds is S'={S
1, S
2..., S
k, and satisfy k 〉=2, then can utilize Ke Lijin to carry out interpolation.The point set coordinate that S' is corresponding is SP={p
1, p
2..., p
k, p wherein
i(1≤i≤k) is three-dimensional coordinate point (x
i, y
i, z
i), I
xCoordinate be (x
0, y
0,
z), z is unknown.Gold layer position interpolation is exactly to utilize the xy coordinate (x of known point set SP and interpolation point in the gram
0, y
0) calculate the process of interpolation point z value.
Step 2, definition layer bit field data structure
Layer bit field data is the set of a string continuous discrete point.Because section is the set of discrete point, so the attribute of section all is unified for all Seed Points on this section, and the logic of relations of tomography is strong.
The data structure of definition phase decryption in this programme:
Seg={ID,Direct,UpOrDown,Start,End,Points,LeftInfo,RightInfo}
ID: the call number of the layer position of section ownership.
Direct: the section direction: X-direction the section be X_Direct, Y direction be Y_Direct.
UpOrDown: the relation of section and tomography: being positioned at hanging wall is UP, and being positioned at footwall is Down.
Start: section starting position.
End: section end position.
The paragraph header of definition phase and section tail related information structure are: Info={ID, UpOrDown}
HeadInfo: the related information Info of the left paragraph header of section.
RearInfo: the related information Info of right section tail of section.
Points: the discrete point set of section.
The set of definition layer position is: S={S
1, S
2... S
i, S
iBe i layer position set.Layer position S
iData structure be
S
i={ ID, XSegments, YSegments}, ID are layer position S
iLayer item, XSegments is for belonging to a layer S
iAll segment datas set of X-direction, YSegments is all segment datas set of Y direction.
Step 3, fault surface match
Tomography is to break because of the stressed some strength that reaches in earth's crust rock stratum, and along the plane of fracture structure that obviously relatively moves is arranged.The face that fault rupture staggers is called fault surface.The sillar of fault surface both sides is called fault wall, is positioned at the dish that is called on the fault surface, is called lower wall under the fault surface.Unique nature according to tomography is divided into: the trap-up of the trap-down rising relative to upper dish that upper dish descends relatively.
Cross section correct needs to determine first the interpolation scope of cross section correct.Because this programme is based on grid, so get the two-dimentional envelope of tomography kind subsegment decryption on grid plan as the interpolation envelope scope of cross section correct.Then adopt the interpolation method of hybrid grid that regmatic network lattice point in the tomography envelope scope is carried out interpolation.
For some layer bit data, possible its raw data has part to violate actual tectonic structure, mainly is following situation: as shown in Figure 2, upper dish or lower wall at tomography, there is and can only has a layer plane, and in the explanation of seismic data, misdata may occur.Therefore need automatic cutting and manual editing's layer position raw data, the data of deletion error are in order to carry out a layer bit space interpolation, thereby simulate correct aspect.
Definition layer bit boundary and tomography are restrained boundary, and then all layers position section between two restrained boundaries has same alike result.If on a two dimensional cross-section, a plurality of layers of position section are arranged between two restrained boundaries, then need these layers position section is combined into a layer position section, concrete combined method is as follows:
(1) the layer position section on the current section is pressed a section starting position sort ascending.
(2) tomography on the current section of statistics.
(3) each tomography is processed.Upper dish is processed: find the layer position section Seg of the current hanging wall of next-door neighbour, then take this layer position section as initial segment, searching the next one can merge with the layer position section of its merging.Seg upgrades to point to and merges back segment, and recurrence merges next section, until certain restrained boundary.Lower wall is processed too.
(4) if there is not tomography, then only need all layer position sections are merged into one section by initial rest position.
Step 5, initialization original segment data
Because being the restrained boundary of tegillum bit boundary and tomography definition, layer position section retrains, so paragraph header and section tail are inevitable and some restrained boundaries are associated.Initialization original segment data are the incidence relation of initialization segments and restrained boundary, as shown in Figure 3.
Step 6, subregion are divided
The purpose that subregion is divided is to make whole layer position be divided into several little subregions, and these subregions all heavily are not worth, thereby makes the heavy tomography of more complicated and the interpolation of trap-up become not the heavily interpolation of value.Subregion divide to be to determine according to the envelope scope of section, and each tomography can corresponding two sub regions, namely go up regional the and lower wall subregion of plate.
The envelope scope of definition section Fi is Pi, and Pi is the two-dimensional projection envelope scope of section Fi on grid.Upper plate zone corresponding to definition tomography Fi is Ui, and the lower wall subregion is Di.Then Ui and Di are the set of the envelope scope of a series of tomographies of being associated of current tomography envelope scope and current tomography Fi.
The concrete grammar that subregion is divided is as follows:
(1) if tomography Fi envelope scope Pi not and other tomography envelope scope overlaids, then the upper lower burrs subregion of tomography Fi is the envelope scope Si of tomography Fi, i.e. Di={Pi}, Ui={Pi}.
(2) if having a plurality of tomographies the envelope scope of Fs...Ft} the envelope scope overlaid of Ps...Pt} and tomography Fi and and tomography Fi upper (descend) coil direct neighbor, then upper (descending) plate zone of tomography Fi is the set of the envelope scope of all relevant tomographies, be Ui (Di)={ Pi, Ps...Pt}.
(3) after the upper lower burrs subregion of all tomographies had been divided, zone, remaining non-tomography layer position was as an independent non-tomography subregion LH.
Then the whole zone of layer position H can be expressed as: H={U1, D1...Ui, Di, LH}.Ui is respectively and lower wall subregion regional with layer upper plate of the tomography Fi that position H is related with Di, and LH is non-fault region.
Step 7, the match of tomography subregion layer plane
(1) initialization upper lower burrs Seed Points
To a certain subregion interpolation, to determine first the Seed Points data of this sub regions interpolation.These Seed Points data heavily are not worth at equatorial projection.
If the data segment structure of i layer position is S
i={ ID, XSegments, YSegments}, ID are the layer item of layer position, and XSegments is all segment data set of X-direction, and YSegments is all segment data set of Y direction.Definition UpSegs is hanging wall kind subsegment, and DownSegs is footwall kind subsegment.Then the set form of UpSegs or DownSegs is
UpSegs={seg1,seg2...segi}
DownSegs={seg1,seg2...segj}
Its stage casing segi or segj are layer position S
iSection set in a certain section, and have following attribute: segi or segj the paragraph header association be current hanging wall, the section tail association of segi or segj be current footwall.If the set of hanging wall subregion Seed Points is UpPoints, the set of lower wall subregion Seed Points is DownPoints.Then the expression formula form of UpPoints and DownPoints is:
UpPoints={p1,p2...pi}
DownPoints={p1,p2...Pj}
Pi, pj are respectively the points in Ui, the Di section in the following formula, and the two is in the two-dimensional projection on the mesh coordinate respectively in the upper and lower plate zone (UpPoints, DownPoints) at tomography.
(2) tomography upper lower burrs subregion face interpolation
With the interpolation method of hybrid grid, the Seed Points take UpPoints and DownPoints as interpolation respectively carries out interpolation to each net point of the upper lower burrs subregion scope of tomography.
(3) the up and down intercepting of actual card of tomography
The layer bitplane data of the tomography upper lower burrs subregion that obtains after the definition interpolation is initial up and down card, and initial up and down card data are larger than the upper lower burrs layer bitplane data of reality in whole tomography upper lower burrs subregion scope, passed through fault surface.Need to come according to the constraint of tomography initial up and down card data are intercepted.For card and tomography are strictly sealed, we seal processing in fact descending card in the cross section.
If line does not pass through tomography between defining, then claim to be communicated with at these 2 at 2.We take one must be on the up and down card data of reality point as initial source point, go out actual up and down card by connective recursive search.In fact the recursive search method of card is as follows:
1) the initial upper disk surface grid of traversal, find one must be on the fault surface grid some p, then putting p must be on the upper disk surface of reality, therefore take a p as initial source point.
2) set up a data base STACK, be used for depositing source point.At first source point p is pressed among the data base STACK.
3) judge whether data station STACK is empty, if be empty, then recursive search finishes.Otherwise from data base STACK, take out 1 q.Traversal point q up and down four some pi(i=1 on mesh coordinate, 2,3,4).
4) if pi can be communicated with q, then pi is joined in fact in the card, simultaneously pi is pressed among the data base STACK.If pi can not be communicated with q, illustrate that then q is the frontier point of card in fact, then changes to q point z value the z value of the fault surface of corresponding net point.
Also be the same for the recursive search of card under the reality and the processing of upper disk surface.
The layer plane match of step 8, non-tomography subregion
(1) initialization seed point
Carry out layer position interpolation fitting iff the data point of selecting to be positioned at non-fault region as Seed Points, then can make between layer plane of the layer plane of non-tomography subregion and tomography subregion transition unsmooth.The present invention both with the data point of non-fault region as Seed Points, again with the borderline region data point of tomography subregion as Seed Points, thereby solved the rough problem of transition between the subregion.For the tomography subregion after the interpolation, with all on two dimensional surface not heavily data points of value also as Seed Points.
(2) interpolation of non-tomography subregion
For the Seed Points that has generated, with the interpolation method of hybrid grid, each net point of the subregion scope of non-tomography is carried out interpolation.Whole layer position surface fitting was finished after interpolation was finished.
The processing at the end of ground, earth's surface is gone out in step 9, layer position pinching, layer position
Layer position pinching is the same with the situation that the end, ground, earth's surface is gone out in layer position in essence, all is to go current aspect is retrained with pinching layer position and bottom surface, ground, earth's surface.So pinching that only need to the current layer position is illegal part is removed just passable with the part of going out the end, ground, earth's surface.
Defining current interpolated layer position is main stor(e)y position S, and the current layer position has the layer position of pinching to be H
i, ground surface is UH, is DH at the bottom of the ground.Definition H
i, UH, DH be the restraint layer position, then have between restraint layer position and the main stor(e)y position relativeness be which layer position upper, which is lower.The disposal route that the end, ground, earth's surface is gone out in then layer position pinching, layer position is as follows:
(1) if current restraint layer position is the earth's surface, then all is gone out the some intercepting on earth's surface (be on the net point Z value less than the Z value of corresponding net point on the ground surface).
(2) if current restraint layer position is the end, ground, the point of then all being gone out the ground end (be on the net point Z value greater than the Z value of corresponding net point on the bottom surface, ground) intercepts.
(3) if current restraint layer position is pinching layer position.If this pinching layer position is on the S of main stor(e)y position, then with the some intercepting of Z value on the net point less than the Z value of the corresponding net point in restraint layer position.If this pinching layer position is under the S of main stor(e)y position, then with the some intercepting of Z value on the net point greater than the Z value of the corresponding net point in restraint layer position.
Claims (4)
1. a spatial complex layer position reconstructing method of dividing based on the zone is characterized in that: comprise the steps:
Step 1, hybrid grid interpolation method: at the little grid of macrolattice nested inside, first macrolattice is carried out interpolation during interpolation, and then utilize the new interpolation point on existing macroreticular Seed Points and the macrolattice that little grid is carried out interpolation;
Step 2, definition layer bit field data structure;
Step 3, fault surface match: get the two-dimentional envelope of tomography kind subsegment decryption on grid plan as the interpolation envelope scope of cross section correct, then adopt the interpolation method of hybrid grid that regmatic network lattice point in the tomography envelope scope is carried out interpolation;
Step 4, layer position raw data editor: the data of deletion error; Definition layer bit boundary and tomography are restrained boundary, if on a two dimensional cross-section, a plurality of layers of position section are arranged between two restrained boundaries, then need these layers position section is combined into a layer position section;
Step 5, initialization original segment data: the incidence relation of initialization segments and restrained boundary;
Step 6, subregion are divided;
Step 7, the match of tomography subregion layer plane:
(1) initialization upper lower burrs Seed Points;
(2) tomography upper lower burrs subregion face interpolation;
(3) the up and down intercepting of actual card of tomography;
The layer plane match of step 8, non-tomography subregion:
(1) initialization seed point: both with the data point of non-fault region as Seed Points, again with the borderline region data point of tomography subregion as Seed Points; And for the tomography subregion after the interpolation, with all on two dimensional surface not heavily data points of value also as Seed Points;
(2) interpolation of non-tomography subregion: for the Seed Points that has generated, adopt the hybrid grid interpolation method that each net point of the subregion scope of non-tomography is carried out interpolation;
The processing at the end of ground, earth's surface is gone out in step 9, layer position pinching, layer position: remove the illegal pinching part in current layer position and go out ground, the earth's surface part at the end.
2. the spatial complex layer position reconstructing method of dividing based on the zone according to claim 1 is characterized in that: the described method that macrolattice is carried out interpolation is:
1) scans one by one the macrolattice point, determine whether macrolattice point needs to carry out interpolation, if there has been data point in macrolattice point, does not then need to carry out interpolation, otherwise need interpolation;
2) to needing the macrolattice point of interpolation, according to artificial inputted search scope, search for current interpolation point macroreticular Seed Points data on every side;
3) according to the sector of artificial input Seed Points is carried out the sector classification;
4) Seed Points to finding in each sector is arranged apart from current macrolattice point distance from small to large according to it;
5) according to each sector Seed Points quantity selected threshold of inputting, and select to join in the Seed Points formation of current macrolattice interpolation apart from the Seed Points data in the nearest threshold range of interpolation net point in this Sector Range;
6) if current macrolattice point does not find effective Seed Points data to carry out interpolation, current macrolattice interpolation not then, the next macrolattice point of diversion treatments, otherwise turn to the 7th) step;
7) with the Seed Points in the formation of current macrolattice interpolation Seed Points current macrolattice point is carried out the Kriging interpolation;
8) if all macrolattice points processed once, whether then check has the failure of macrolattice point interpolation, if the failure of macrolattice point interpolation is arranged, then with the macrolattice point of interpolation success also as Seed Points, return the 1st) step is carried out again interpolation, until all net point interpolation successes.
3. the spatial complex layer position reconstructing method of dividing based on the zone according to claim 1 is characterized in that: the described method that a plurality of layers of position section are combined into a layer position section is:
(1) the layer position section on the current section is pressed a section starting position sort ascending;
(2) tomography on the current section of statistics;
(3) each tomography is processed;
(4) if there is not tomography, then only need all layer position sections are merged into one section by initial rest position.
4. the spatial complex layer position reconstructing method of dividing based on the zone according to claim 1, it is characterized in that: the method that described subregion is divided is:
(1) if certain tomography envelope scope does not have and other tomography envelope scope overlaids, then the upper lower burrs subregion of this tomography is the envelope scope of tomography;
(2) if having the envelope scope of a plurality of tomographies and a certain tomography envelope scope overlaid and and upper dish or the lower wall direct neighbor of this tomography, then the upper dish of this tomography or lower wall subregion are the set of the envelope scope of all relevant tomographies;
(3) after the upper lower burrs subregion of all tomographies had been divided, zone, remaining non-tomography layer position was as an independent non-tomography subregion.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210311567.1A CN102867330B (en) | 2012-08-29 | 2012-08-29 | Region-division-based spatial complex horizon reconstruction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210311567.1A CN102867330B (en) | 2012-08-29 | 2012-08-29 | Region-division-based spatial complex horizon reconstruction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102867330A true CN102867330A (en) | 2013-01-09 |
| CN102867330B CN102867330B (en) | 2014-10-01 |
Family
ID=47446186
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210311567.1A Expired - Fee Related CN102867330B (en) | 2012-08-29 | 2012-08-29 | Region-division-based spatial complex horizon reconstruction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102867330B (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103489222A (en) * | 2013-09-06 | 2014-01-01 | 电子科技大学 | Target body surface reconstruction method in three-dimensional image |
| CN104166163A (en) * | 2014-08-27 | 2014-11-26 | 电子科技大学 | Method for automatically extracting fault curved surface based on three-dimensional large-data-volume seismic data cube |
| CN104182571A (en) * | 2014-08-12 | 2014-12-03 | 电子科技大学 | Kriging interpolation method based on Delaunay and GPU |
| CN104240301A (en) * | 2014-09-11 | 2014-12-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Geological surface reconstruction method and device |
| CN104613945A (en) * | 2015-02-11 | 2015-05-13 | 国家***第二海洋研究所 | Reconstruction method for terrain of shallow-sea large-sized complicated sand wave area |
| CN104656131A (en) * | 2013-11-22 | 2015-05-27 | 中国石油天然气集团公司 | Method for improving seismic horizon automatic tracking precision based on waveform recovery technique |
| CN105612560A (en) * | 2013-08-16 | 2016-05-25 | 界标制图有限公司 | Generating a custom reservoir from multiple compartments representing one or more geological structures |
| CN104199099B (en) * | 2014-09-19 | 2016-08-24 | 东北石油大学 | A kind of method utilizing multiwindow serial section efficiently to carry out seismic interpretation |
| CN112634450A (en) * | 2021-01-06 | 2021-04-09 | 安徽理工大学 | Method for integrally constructing local complex fault through corner grid |
| CN113379907A (en) * | 2020-03-10 | 2021-09-10 | 中国石油天然气股份有限公司 | Method and device for constructing fault block geological model |
| CN113674296A (en) * | 2021-09-03 | 2021-11-19 | 广东三维家信息科技有限公司 | Region cutting method and device, electronic equipment and storage medium |
| CN115438604A (en) * | 2022-11-08 | 2022-12-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Grid identification method based on prime number system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303414A (en) * | 2008-05-22 | 2008-11-12 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
| WO2010082969A1 (en) * | 2009-01-13 | 2010-07-22 | Exxonmobil Upstream Research Company | Methods and systems to volumetrically conceptualize hydrocarbon plays |
-
2012
- 2012-08-29 CN CN201210311567.1A patent/CN102867330B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101303414A (en) * | 2008-05-22 | 2008-11-12 | 北京航空航天大学 | Method for generating ground layer surface and geologic body based on level set |
| WO2010082969A1 (en) * | 2009-01-13 | 2010-07-22 | Exxonmobil Upstream Research Company | Methods and systems to volumetrically conceptualize hydrocarbon plays |
Non-Patent Citations (1)
| Title |
|---|
| 朱良峰等: "地质断层三维可视化模型的构建方法与实现技术", 《软件学报》, vol. 19, no. 8, 15 August 2008 (2008-08-15), pages 2004 - 2017 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105612560A (en) * | 2013-08-16 | 2016-05-25 | 界标制图有限公司 | Generating a custom reservoir from multiple compartments representing one or more geological structures |
| CN103489222A (en) * | 2013-09-06 | 2014-01-01 | 电子科技大学 | Target body surface reconstruction method in three-dimensional image |
| CN103489222B (en) * | 2013-09-06 | 2016-06-22 | 电子科技大学 | Target body surface reconstruction method in 3-D view |
| CN104656131A (en) * | 2013-11-22 | 2015-05-27 | 中国石油天然气集团公司 | Method for improving seismic horizon automatic tracking precision based on waveform recovery technique |
| CN104182571A (en) * | 2014-08-12 | 2014-12-03 | 电子科技大学 | Kriging interpolation method based on Delaunay and GPU |
| CN104166163B (en) * | 2014-08-27 | 2016-08-24 | 电子科技大学 | Tomography curved surface extraction method based on three-dimensional big data quantity seismic data cube |
| CN104166163A (en) * | 2014-08-27 | 2014-11-26 | 电子科技大学 | Method for automatically extracting fault curved surface based on three-dimensional large-data-volume seismic data cube |
| CN104240301A (en) * | 2014-09-11 | 2014-12-24 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Geological surface reconstruction method and device |
| CN104199099B (en) * | 2014-09-19 | 2016-08-24 | 东北石油大学 | A kind of method utilizing multiwindow serial section efficiently to carry out seismic interpretation |
| CN104613945A (en) * | 2015-02-11 | 2015-05-13 | 国家***第二海洋研究所 | Reconstruction method for terrain of shallow-sea large-sized complicated sand wave area |
| CN113379907A (en) * | 2020-03-10 | 2021-09-10 | 中国石油天然气股份有限公司 | Method and device for constructing fault block geological model |
| CN112634450A (en) * | 2021-01-06 | 2021-04-09 | 安徽理工大学 | Method for integrally constructing local complex fault through corner grid |
| CN113674296A (en) * | 2021-09-03 | 2021-11-19 | 广东三维家信息科技有限公司 | Region cutting method and device, electronic equipment and storage medium |
| CN115438604A (en) * | 2022-11-08 | 2022-12-06 | 中国空气动力研究与发展中心计算空气动力研究所 | Grid identification method based on prime number system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102867330B (en) | 2014-10-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102867330B (en) | Region-division-based spatial complex horizon reconstruction method | |
| CN102222365B (en) | Method for reconstructing curved surface of complex space | |
| Ming et al. | GSIS: A 3D geological multi-body modeling system from netty cross-sections with topology | |
| CN103454686B (en) | Small scale sedimentary facies based on strata slicing carries out the method and system of reservoir prediction | |
| CN101582173B (en) | Block model building method for complex geological structure | |
| CN101906965B (en) | Technology for identifying and analyzing stratums in VR (Virtual Reality) geological environment | |
| US10995592B2 (en) | Method and system for analyzing the uncertainty of subsurface model | |
| EP3329307B1 (en) | Assignment of systems tracts | |
| CN101303414B (en) | Method for generating ground layer surface and geologic body based on level set | |
| CN103901478A (en) | Method for determining deposition characteristics and distribution of reservoirs by combining logging and seismic information | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| CN106687827B (en) | Stratum modeling method for fault | |
| CN103700141B (en) | Multistage dimension geologic modeling analysis method of coal seam floor water gushing passageway | |
| CN103699751A (en) | Sand body reservoir architecture modeling method and system based on space vectors | |
| CN102609982B (en) | Topology discovery method of space geological data based on unstructured mode | |
| CN102867332B (en) | Based on the multistage subdivided meshes curved surface fitting method of complex boundary constraint | |
| CN106097445A (en) | A kind of method for drafting of novel three-dimensional stratum curved surface | |
| CN106355640A (en) | Ground-underground three-dimensional integrated space processing method and system | |
| Wu et al. | A 3D modeling approach to complex faults with multi-source data | |
| CN107481320A (en) | A kind of gridless routing of three-dimensional geological modeling | |
| Zhang et al. | 3D parametric modeling of complex geological structures for geotechnical engineering of dam foundation based on T‐splines | |
| CN102903149B (en) | Block forming method and device for geologic model | |
| CN108646294B (en) | Deep water turbidite sandstone reservoir high-precision inversion method under a kind of complicated fault system | |
| Ji et al. | An automated method to build 3D multi-scale geological models for engineering sedimentary layers with stratum lenses | |
| CN102750738A (en) | Three-dimensional geologic body self-generating technology based on plane drilling |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141001 Termination date: 20170829 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |