CN110244021A - A kind of strata division method based on anisotropy interpolation - Google Patents
A kind of strata division method based on anisotropy interpolation Download PDFInfo
- Publication number
- CN110244021A CN110244021A CN201910516076.2A CN201910516076A CN110244021A CN 110244021 A CN110244021 A CN 110244021A CN 201910516076 A CN201910516076 A CN 201910516076A CN 110244021 A CN110244021 A CN 110244021A
- Authority
- CN
- China
- Prior art keywords
- point
- attribute value
- interpolation
- small lattice
- centric
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention belongs to strata division technical fields, there is provided a kind of strata division method based on anisotropy interpolation: exploratory bore-hole records formation lithology title and each formation lithology depth bounds, cube net is constructed according to exploratory bore-hole, drilling is separated into multiple data points along depth direction.A formation lithology title is selected, whether is fallen in the formation lithology depth bounds according to discrete data point in each drilling, differentiation assignment is carried out to the discrete data point in each drilling, as sample point.According to sample point attribute value, anisotropy interpolation further is carried out according to specified trend in cube net, generates attribute contour surface.By aforesaid operations, the corresponding attribute contour surface of each lithology title is sequentially generated, strata division is divided according to each lithology contour surface.This method repeats the case where more, alternating layers development etc. causes manual zoning difficult or can not be layered for place lithology is encountered during practical prospecting.The layering result of complex geological condition can quickly and accurately be obtained.
Description
Technical field
The invention belongs to strata division technical fields, and in particular to a kind of strata division side based on anisotropy interpolation
Method.
Background technique
Soil layer is practical to reconnoitre work because of origin cause of formation multiplicity, complicated, the typical feature that generally existing layer position repeats, lenticular body is developed
In work, when encountering this kind of geological conditions, great difficulty is brought to manual zoning, side is examined in traditional single hole layering, two-dimentional school
Formula does not ensure that space reasonability, the correctness of layering result.The formation geology epoch are relatively clear, and deposition rule is preferable
Rock stratum can complete stratified operation by manual single hole operation mode, when meeting construction influence or complex geologic conditions, artificial point
There is also bigger difficulty for layer, or even have the case where can not being layered.Such as the generally existing red beds in Southwest China, South China
Layer, lithology are mainly mud stone, sandstone alternating layers, there are consistent occurrence in two kinds of rock stratum and interval occurs, spatially each exploratory hole
How lithologic interface is layered, line is extremely difficult, carries out stratigraphic horizon division by hand by geological personnel, generally requires to expend a large amount of
Time and efforts, or even layer position can not be marked off at all.The lamination problem encountered in real work significantly reduces prospecting
Working efficiency causes difficulty to subsequent three-dimensional geological modeling, ground design, influences entire construction period.
Summary of the invention
The purpose of the present invention is overcome layering error rate height, low efficiency, special case in traditional batch technology that can not be layered
The problem of.
For this purpose, S1: carrying out ground in Project Areas the present invention provides a kind of strata division method based on anisotropy interpolation
Matter exploration, acquires the lithologic interface of each exploratory bore-hole, and the information data of the lithologic interface includes formation lithology title and stratum
Lithology depth bounds；
S2: importing the information data of the lithologic interface of each drilling, founds big cube of net in the building of entire exploration space, should
Big cube of net includes multiple small lattices, with reference to the big wisp drilling of small lattice along depth direction it is discrete obtain it is multiple discrete
Data point；
S3: one formation lithology title of selection, it is discrete in the formation lithology depth bounds by falling in each drilling
The attribute value of data point is assigned a value of 1, and the attribute value of remaining data points is assigned a value of 0, finally by the discrete data point of these initialization
As sample point；
S4: each small cube is obtained using anisotropy Kriging interpolation method according to the attribute value of the sample point
The attribute value of grid position central point traverses the attribute value of all small lattice place-centric points, attribute value is greater than pre-
If the attribute value of the small lattice place-centric point of median is changed to 1, less than in the small lattice position of this median
The attribute value of heart point is changed to 0, and the attribute value of the sample point is assigned to the position of the nearest small lattice of distance sample point
Central point；
S5: in big cube of net, all properties value is connected generation etc. for the place-centric point of 1 small lattice
Value face；
S6: repeating step S3 to S5, successively creates the contour surface of each formation lithology in the drilling；
S7: according to the contour surface of each formation lithology, the strata division of the exploration area is completed.
Preferably, the default median is 0.5.
Preferably, the step S4 is specifically included: according to the sample point attribute value, being inserted using anisotropy Kriging
Value method, obtains the attribute value of each small lattice place-centric point, and the calculating of each attribute value includes and sample data
Correlation degree, the correlation degree include weight relationship horizontally and vertically, the distribution of surrounding sample position；Traversal institute
Attribute value is greater than the small lattice place-centric point of default median by the attribute value for the small lattice place-centric point having
Attribute value be changed to 1, the attribute value less than the small lattice place-centric point of default median is changed to 0, and the sample
The attribute value of point is assigned to the place-centric point of the nearest small lattice of distance sample point；If each small lattice place-centric point
Attribute value be 1, then it is assumed that the small lattice place-centric point belongs to the corresponding stratum of the lithologic interface, if each small cube of net
The attribute value of lattice place-centric point is 0, then it is assumed that the small lattice place-centric point is not belonging to the lithologic interface correspondingly
Layer.
Preferably, the step S4 further include: obtained respectively by plane of reference interpolation calculation of the occurrence face at formation lithology interface
The attribute value of the small lattice place-centric point；When lithologic interface does not have occurrence face, using earth's surface trend surface as the plane of reference
Interpolation calculation obtains the attribute value of each small lattice place-centric point.
Preferably, the formation lithology title such as plain fill, mud, silty clay, silt and flour sand, sandstone, mud stone etc..
Preferably, a common contour surface is constituted if the contour surface in two neighboring drilling is connected.
Preferably, small lattice place-centric point is linked to be generation contour surface, the space in the contour surface forms one
A equivalence layer.
Preferably, after the step S5 further include: using smooth discrete interpolation method that contour surface progress is excellent
Change handles to obtain smooth contour surface.
Preferably, the strata division method further include: by lithologic interface with a thickness of zero bore position be added to etc.
In value face, to determine the pinching position on stratum.
Beneficial effects of the present invention: this strata division method based on anisotropy interpolation provided by the invention, in work
Journey area carries out geological prospecting, and the lithologic interface and catalog data, lithologic interface acquired in each exploratory spot drilling contains stratum
Lithology Name & Location information, the formation lithology disclosed exploratory spot are sorted out, and all stratum that engineering ground discloses are obtained
Lithology." cube net " is constructed within the scope of exploration area, " cube net " is a kind of space lattice set, by certain amount, uniformly big
The small lattice filling of small hexahedron is constituted.A certain lithology title is chosen, assigns specific geological property value, in each drilling
The lithologic interface position will be assigned this attribute value, and the corresponding attribute value in lithologic interface position in each drilling is passed to most
Close small lattice becomes to attribute value according to specified in small lattice according to anisotropy Kriging interpolation method
Gesture carries out each to differentiation interpolation, and after the completion of interpolation, each small lattice is attached to geological property, attached according to small lattice
The geological property value of band generates contour surface.Continue to choose next lithology title, assign specific geological property value, repeats above-mentioned insert
Value Operations, successively creation obtains the corresponding contour surface of each lithology title, divides strata division code name according to each lithology contour surface.It should
Method for encountered during practical prospecting place lithology repeat more, alternating layers development (red beds) etc. cause manual zoning difficult or
The case where can not being layered.The layering of complex geological condition can quickly and accurately be obtained as a result, serve future three-dimensional geological modeling
And analysis work.
The present invention is described in further details below with reference to attached drawing.
Detailed description of the invention
Fig. 1 is the flow diagram of the strata division method the present invention is based on anisotropy interpolation；
Fig. 2 is that the present invention is based on the occurrence face figures in the reference trend of the strata division method of anisotropy interpolation；
Fig. 3 is the reference isogram of the strata division method the present invention is based on anisotropy interpolation；
Fig. 4 is the interpolation contour surface without reference to trend surface of the strata division method the present invention is based on anisotropy interpolation
Figure；
Fig. 5 is the interpolation contour surface of the strata division method the present invention is based on anisotropy interpolation having with reference to trend surface
Figure；
Fig. 6 is the formation lithology interface distributions ordering chart of the strata division method the present invention is based on anisotropy interpolation；
Fig. 7 is the single interface debonding reference model of the strata division method the present invention is based on anisotropy interpolation；
Fig. 8 is a variety of interface debonding reference models of the strata division method the present invention is based on anisotropy interpolation.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts all other
Embodiment shall fall within the protection scope of the present invention.
In the description of the present invention, it is to be understood that, term " center ", "upper", "lower", "front", "rear", " left side ",
The orientation or positional relationship of the instructions such as " right side ", "vertical", "horizontal", "top", "bottom", "inner", "outside" is based on the figure
Orientation or positional relationship is merely for convenience of description of the present invention and simplification of the description, rather than the device of indication or suggestion meaning or
Element must have a particular orientation, be constructed and operated in a specific orientation, therefore be not considered as limiting the invention.
Term " first ", " second " be used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance or
Implicitly indicate the quantity of indicated technical characteristic." first " is defined as a result, the feature of " second " can be expressed or imply
Ground includes one or more of the features；In the description of the present invention, unless otherwise indicated, the meaning of " plurality " is two or
It is more than two.
The embodiment of the invention provides a kind of strata division method based on anisotropy interpolation, S1: carried out in Project Areas
Geological prospecting, acquires the lithologic interface of each exploratory bore-hole, and the information data of the lithologic interface includes formation lithology title and ground
Layer lithology depth bounds.Then the first drilling in exploration area acquires the geologic data in drilling, has in each drilling very much
Layer, one lithologic interface of each layer of correspondence, each lithologic interface are made of the identical discrete point grid of many geological properties, often
A lithologic interface has a formation lithology title and formation lithology depth bounds.Each lithologic interface can correspond to an occurrence
Face.
S2: importing the information data of the lithologic interface of each drilling, constructs cube net in entire exploration space, this is big vertical
Square net includes multiple small lattices, obtains multiple discrete datas along depth direction is discrete with reference to the drilling of small lattice big wisp
Point.It follows that the corresponding small lattice of each discrete data point.The lithologic interface of all drillings all corresponds to a depth
Position, but the depth that can have the identical lithologic interface of different drillings is different, thickness is also different, it is therefore desirable to which classification is named
Word realizes the initialization of bed boundary data, convenient for establishing equivalent layer below.As shown in Table 1 and Table 2, table 1 is some drilling
Interior stratigraphic record data, table 2 are the tables for the data initialization for after importing for all borehole formation data and naming formation.
The 1 initial layered sheet of Project Areas borehole formation of table
| Layer bottom depth | Stratum code name | Stratigraphic Time | Lithology title |
| 1.5 | 1 | Q_0_ml | Plain fill |
| 8 | 3 | Q_h_al+1 | Silty clay |
| 11.2 | 5 | Q_p_al | Flour sand |
| 16.8 | 3 | Q_h_al+1 | Silty clay |
| 25 | 5 | Q_p_al | Flour sand |
The initial lithology layered sheet of 2 exploratory spot of table
| Affiliated engineering | Stratum generation | Stratigraphic Time | Ground class | Ground title | Boundary | System |
| Lithology modeling | 1 | Q_0_ml | Soil layer | Plain fill | Cenozoic | 4th system |
| Lithology modeling | 2 | Q_h_1 | Soil layer | Mud | Cenozoic | 4th system |
| Lithology modeling | 3 | Q_h_al+1 | Soil layer | Silty clay | Cenozoic | 4th system |
| Lithology modeling | 4 | Q_p_al | Soil layer | Silt | Cenozoic | 4th system |
| Lithology modeling | 5 | Q_p_al | Soil layer | Flour sand | Cenozoic | 4th system |
S3: one formation lithology title of selection, it is discrete in the formation lithology depth bounds by falling in each drilling
The attribute value of data point is assigned a value of 1, and the attribute value of remaining data points is assigned a value of 0, finally by the discrete data point of these initialization
As sample point.A formation lithology title is selected, then the corresponding formation lithology depth bounds of the formation lithology title will
The attribute value of the discrete data point fallen in the formation lithology depth bounds in others drilling is assigned a value of 1, remaining data points
Attribute value be assigned a value of 0, the identical all discrete data points of formation lithology depth bounds of the formation lithology title can be obtained.
S4: each small cube is obtained using anisotropy Kriging interpolation method according to the attribute value of the sample point
The attribute value of grid position central point traverses the attribute value of all small lattice place-centric points, attribute value is greater than pre-
If the attribute value of the small lattice place-centric point of median is changed to 1, less than in the small lattice position of this median
The attribute value of heart point is changed to 0, and the attribute value of the sample point is assigned to the position of the nearest small lattice of distance sample point
Central point.It follows that also assigning the point around discrete point to attribute value 1 or 0 according to Kriging (Ke Lijin) interpolation method, belong to
Property value be 1 to indicate this discrete point and reference sample point is the same formation lithology title, attribute value is 0 not to be then.
S5: in big cube of net, all properties value is connected generation etc. for the place-centric point of 1 small lattice
Value face.It is also used when being attached between the discrete point (i.e. particle) and the discrete point of surrounding same names in big cube of net different
Property Kriging interpolation method realize smooth transition.
S6: repeating step S3 to S5, successively creates the contour surface of each formation lithology in the drilling；
S7: according to the contour surface of each formation lithology, the strata division of the exploration area is completed.
Specifically, it as shown in fig. 7, first passing through the formation lithology title and formation lithology depth bounds of drilling, and utilizes
Kriging (Ke Lijin) interpolation algorithm obtains the identical discrete points data of a formation lithology title, that is, selects a formation lithology
The attribute value of the discrete data point fallen in the formation lithology depth bounds is assigned a value of 1 in each drilling by title,
The attribute value at remainder strong point is assigned a value of 0, finally using the discrete data point of these initialization as sample point.The sample point includes boring
In hole and outside drilling, sample point connection then can be formed into the corresponding contour surface of formation lithology title.So circulation, just
Each formation lithology title can be connected into corresponding contour surface, can so complete the strata division to exploration area.
The corresponding attribute of all small lattices in spatial dimension is calculated according to lithologic interface, and lithology circle
The attribute of the location point in face is mapped in small lattice, generates multiple equivalent spaces to get stratigraphic model approximate shape is arrived.
As shown in figure 8, realizing all stratum in drilling according to each layer of incidence relation and the overlying relation on stratum
Layering on the whole obtains layering reference result.All lithologic interfaces all form respectively corresponding contour surface, so
The geology of the exploration area is presented in a manner of strata division.
The present invention is analyzed and has been studied to the distribution character on stratum or mineral reserve under natural conditions, it is found that it is usually in
Band-like, block distribution, or there are any trend surfaces, such as underground reservoir, and a main production can be determined according to the trend rule
Shape direction.If simulated object is distributed along the occurrence, then it is assumed that association of this object between its attribute of this direction is obvious, other
Direction is then opposite；According to occurrence face or earth's surface trend surface, the formation data of different directions is different to unknown point weighing factor to be inserted
Value method, i.e. anisotropy interpolation.In a specific implement scene, as shown in Figure 1, geological prospecting point data is first imported,
The formation lithology interface data that drilling discloses is obtained, the borehole formation lithologic interface data of importing are then automatically processed.Selection one
When a formation lithology title, the category of the discrete data point in the corresponding formation lithology depth bounds of the formation lithology title will be fallen in
Property value be assigned a value of 1, the attribute value of remaining data points is assigned a value of 0, using these initialization discrete data point as sample point, finally
Place-centric point all properties value for 1 small lattice connects generation contour surface.Wherein, such as Fig. 2 to Fig. 5 institute
Show, when to small lattice attribute value assignment around sample point, to be inserted by the plane of reference of the occurrence face at formation lithology interface
The attribute value that each small lattice place-centric point is calculated in value is the most accurate.
Preferred scheme, the default median are 0.5.The value is empirical data.
Preferred scheme, in the step S3 further include: if not searching out the lithology with adjacent holes in a drilling
The thickness of the one corresponding lithologic interface to drill is then set to by the Name & Location attribute at interface lithologic interface all the same
Zero.It follows that the number of plies of the lithologic interface in each drilling is identical under normal circumstances, only thickness is had differences, but special
In the case of, such as there are when pinching situation, then the lithologic interface thickness of this layer is defined as zero, i.e., this with a thickness of zero lithology
The top layer at interface and stratum are overlapped will should be with a thickness of zero formation lithology interface and same position attribute and thickness position zero
It when formation lithology interface connects, is attached also with Kriging regression algorithm, is only somebody's turn to do the formation lithology with a thickness of zero
The top and bottom of interface sample are coplanar.
Preferred scheme, the formation lithology title such as plain fill, mud, silty clay, silt and flour sand, sandstone, mud
Rock etc..As shown in fig. 6, it is known from figures that, formation lithology title in case drilling there are plain fill, silty clay,
Three kinds of flour sand etc., there may be others in other cases, be not listed one by one herein.
Preferred scheme, the step S4 are specifically included: according to the sample point attribute value, using anisotropy
Kriging interpolation method, obtains the attribute value of each small lattice place-centric point, and the calculating of each attribute value includes and sample
The correlation degree of notebook data, the correlation degree include weight relationship horizontally and vertically, surrounding sample position point
Cloth；Attribute value is greater than the small lattice of default median by the attribute value for traversing all small lattice place-centric points
The attribute value of place-centric point is changed to 1, and the attribute value less than the small lattice place-centric point of default median is changed to 0, and
The attribute value of the sample point is assigned to the place-centric point of the nearest small lattice of distance sample point；If each small lattice
The attribute value of place-centric point is 1, then it is assumed that the small lattice place-centric point belongs to the corresponding stratum of the lithologic interface, if
The attribute value of each small lattice place-centric point is 0, then it is assumed that the small lattice place-centric point is not belonging to lithology circle
The corresponding stratum in face.
It follows that longitudinal specific setting method of relating value is: being first set to net according to the minimum thickness in all stratum
The several times of lattice (i.e. small lattice) longitudinal direction size, longitudinal size of grid can serve as longitudinal relating value, with described big
The average value of the position attribution of all discrete points in cube net is longitudinal relating value, by each small lattice place-centric point
The attribute value of the identical environment discrete point of attribute value is compared with longitudinal Attribute Association value of the small lattice, if more than or
Equal to the longitudinal direction Attribute Association value, then it is assumed that the place-centric point of the small lattice belongs to rock corresponding to the small lattice
Property interface, i.e., attribute value be 1, if be less than the longitudinal direction relating value, then it is assumed that the place-centric point of the small lattice is not belonging to this
Corresponding lithologic interface described in small lattice, i.e. the property value of the small lattice place-centric point are 0, then according to each to different
Property Kriging interpolation method is by the place-centric point of each small lattice and belongs to the corresponding lithologic interface of the small lattice
Environment discrete point connection generates equivalent layer.It follows that all small lattices of traversal, each small lattice position, root
According to anisotropy Kriging interpolation method, position attribution is obtained, the calculating of each position attribution contains the association of sample data
Degree, the level of setting and longitudinal weight relationship, the position distribution etc. of surrounding sample.The attribute of each mesh point is each to different
Property Kriging interpolation algorithm obtain, add original sample data, each position has category in entire cube of net space
Property value.Each attribute value is compared with median, greater than the point for being considered as belonging to stratum of median, attribute modification 1；In being less than
Between the point for being taken as being not belonging to stratum that is worth, attribute modification 0.In big cube of net, created according to small lattice equivalent
For face algorithm by strata division, each contour surface finally obtained seeks to the strata division model of creation.
Preferred scheme, the step S4 further include: obtained by plane of reference interpolation calculation of the occurrence face at formation lithology interface
To the attribute value of each small lattice place-centric point；It is ginseng with earth's surface trend surface when lithologic interface does not have occurrence face
The face interpolation calculation of examining obtains the attribute value of each small lattice place-centric point.Come in terms of interpolation using occurrence face as the plane of reference
Calculate the attribute value of other ground particles (i.e. discrete point) on a discrete point periphery.
Preferred scheme constitutes a common contour surface if the contour surface in two neighboring drilling is connected.If not yet
There is intersection, as the distance of extension is elongated, top surface equivalence layer and bottom surface equivalence layer can also intersect.
Preferred scheme constructs the correlation model of Kriging interpolation, and the correlation model is along earth's surface trend surface direction
Correlation degree is better than along the earth's surface Trend simulating method to correlation degree.The key of anisotropy interpolation is that Kriging to be constructed is inserted
The correlation model (variation function) of value, correlation model is strong along trend surface directional correlation degree, and correlation distance is big, and Trend simulating method to
Correlation degree is weak, and correlation distance is small, so that the interpolation range of bed boundary position attribution, the just region near trend surface.
Preferred scheme, the strata division method further include: by being added with a thickness of zero bore position for lithologic interface
Into equivalent layer.Bore position with a thickness of 0, that is, pinching drilling is also added in sample, so that in model adjust automatically stratum
Pinching position.Pinching is a kind of geological phenomenon of geologic horizon distribution, is present in sedimentation, it is the art of geological prospecting
Language.It is to say that a certain clay distribution is extremely uneven, is surveyed in range at us, there is deficient phenomena.It is visible on log sheet to be somebody's turn to do
Soil layer fades away, and the sharp figure sharply terminated is presented in thickness.
Beneficial effects of the present invention: this strata division method based on anisotropy interpolation provided by the invention, in work
Journey area carries out geological prospecting, and the lithologic interface and catalog data, lithologic interface acquired in each exploratory spot drilling contains stratum
Lithology Name & Location information, the formation lithology disclosed exploratory spot are sorted out, and all stratum that engineering ground discloses are obtained
Lithology." cube net " is constructed within the scope of exploration area, " cube net " is a kind of space lattice set, by certain amount, uniformly big
The small lattice filling of small hexahedron is constituted.A certain lithology title is chosen, assigns specific geological property value, in each drilling
The lithologic interface position will be assigned this attribute value, and the corresponding attribute value in lithologic interface position in each drilling is passed to most
Close small lattice becomes to attribute value according to specified in small lattice according to anisotropy Kriging interpolation method
Gesture carries out each to differentiation interpolation, and after the completion of interpolation, each small lattice is attached to geological property, attached according to small lattice
The geological property value of band generates contour surface.Continue to choose next lithology title, assign specific geological property value, repeats above-mentioned insert
Value Operations, successively creation obtains the corresponding contour surface of each lithology title, divides strata division code name according to each lithology contour surface.It should
Method for encountered during practical prospecting place lithology repeat more, alternating layers development (red beds) etc. cause manual zoning difficult or
The case where can not being layered.The layering of complex geological condition can quickly and accurately be obtained as a result, serve future three-dimensional geological modeling
And analysis work.
The foregoing examples are only illustrative of the present invention, does not constitute the limitation to protection scope of the present invention, all
It is within being all belonged to the scope of protection of the present invention with the same or similar design of the present invention.
Claims (8)
1. a kind of strata division method based on anisotropy interpolation, it is characterised in that:
S1: geological prospecting is carried out in Project Areas, acquires the lithologic interface of each exploratory bore-hole, the information packet of the lithologic interface
Title containing formation lithology and formation lithology depth bounds；
S2: importing the information data of the lithologic interface of each drilling, founds big cube of net in the building of entire exploration space, this is big vertical
Square net includes multiple small lattices, obtains multiple discrete datas along depth direction is discrete with reference to the drilling of small lattice big wisp
Point；
S3: one formation lithology title of selection will fall in the discrete data in the formation lithology depth bounds in each drilling
Point attribute value be assigned a value of 1, the attribute value of remaining data points is assigned a value of 0, finally using these initialization discrete data point as
Sample point；
S4: each small lattice is obtained using anisotropy Kriging interpolation method according to the attribute value of the sample point
The attribute value of place-centric point traverses the attribute value of all small lattice place-centric points, and attribute value is greater than in default
Between the attribute value of small lattice place-centric point that is worth be changed to 1, less than the small lattice place-centric point of this median
Attribute value be changed to 0, and the attribute value of the sample point is assigned to the place-centric of the nearest small lattice of distance sample point
Point；
S5: in big cube of net, all properties value is connected generation equivalence for the place-centric point of 1 small lattice
Face；
S6: repeating step S3 to S5, successively creates the contour surface of each formation lithology in the drilling；
S7: according to the contour surface of each formation lithology, the strata division of the exploration area is completed.
2. the strata division method according to claim 1 based on anisotropy interpolation, it is characterised in that: it is described it is default in
Between value be 0.5.
3. the strata division method according to claim 1 based on anisotropy interpolation, which is characterized in that the step S4
It specifically includes: each small lattice is obtained using anisotropy Kriging interpolation method according to the sample point attribute value
The calculating of the attribute value of place-centric point, each attribute value includes the correlation degree with sample data, and the correlation degree includes
Weight relationship horizontally and vertically, the distribution of surrounding sample position；Traverse all small lattice place-centric points
Attribute value, the attribute value that attribute value is greater than the small lattice place-centric point of default median is changed to 1, is less than in default
Between the attribute value of small lattice place-centric point that is worth be changed to 0, and the attribute value of the sample point is assigned to distance sample point
The place-centric point of nearest small lattice；If the attribute value of each small lattice place-centric point is 1, then it is assumed that this is small vertical
Square grid position central point belongs to the corresponding stratum of the lithologic interface, if the attribute value of each small lattice place-centric point is 0,
Then think that the small lattice place-centric point is not belonging to the corresponding stratum of the lithologic interface.
4. the strata division method according to claim 3 based on anisotropy interpolation, which is characterized in that the step S4
Further include: each small lattice place-centric point is obtained by plane of reference interpolation calculation of the occurrence face at formation lithology interface
Attribute value；When lithologic interface does not have occurrence face, each small cube of net is obtained by plane of reference interpolation calculation of earth's surface trend surface
The attribute value of lattice place-centric point.
5. the strata division method according to claim 1 based on anisotropy interpolation, it is characterised in that: the formation rock
Property title such as plain fill, mud, silty clay, silt and flour sand, sandstone, mud stone etc..
6. the strata division method according to claim 1 based on anisotropy interpolation, it is characterised in that: by small cube of net
Lattice place-centric point is linked to be generation contour surface, and the space in the contour surface forms an equivalent layer.
7. the strata division method according to claim 7 based on anisotropy interpolation, which is characterized in that the step S5
Later further include: the contour surface is optimized into processing using smooth discrete interpolation method and obtains smooth contour surface.
8. the strata division method according to claim 1 based on anisotropy interpolation, which is characterized in that the stratum point
Layer method further include: by being added with a thickness of zero bore position into contour surface for lithologic interface, to determine the pinching position on stratum
It sets.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910516076.2A CN110244021B (en) | 2019-06-14 | 2019-06-14 | Stratum layering method based on anisotropic interpolation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910516076.2A CN110244021B (en) | 2019-06-14 | 2019-06-14 | Stratum layering method based on anisotropic interpolation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110244021A true CN110244021A (en) | 2019-09-17 |
| CN110244021B CN110244021B (en) | 2021-08-20 |
Family
ID=67887151
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910516076.2A Active CN110244021B (en) | 2019-06-14 | 2019-06-14 | Stratum layering method based on anisotropic interpolation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110244021B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110837668A (en) * | 2019-10-21 | 2020-02-25 | 中铁第四勘察设计院集团有限公司 | BIM software-based method and device for calculating earth volume of different soil qualities |
| CN112419500A (en) * | 2020-12-09 | 2021-02-26 | 上海申元岩土工程有限公司 | Three-dimensional geological model modeling method |
| CN112686779A (en) * | 2020-12-16 | 2021-04-20 | 天津泰勘工程技术咨询有限公司 | Method for dividing sequence of standard soil layers of foundation soil |
| CN117152369A (en) * | 2023-09-08 | 2023-12-01 | 中南大学 | Construction method, device, equipment and computer readable storage medium of frozen soil model |
| CN118015219A (en) * | 2024-04-09 | 2024-05-10 | 中国地质大学（武汉） | Geological model generation method, device and equipment based on qualitative kriging interpolation |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102147933A (en) * | 2010-02-05 | 2011-08-10 | 王远功 | Method for displaying stratum structure in three-dimensional scene |
| CN106558100A (en) * | 2016-10-25 | 2017-04-05 | 中国电建集团华东勘测设计研究院有限公司 | A kind of stratigraphic model method for automatic modeling based on borehole data |
| CN106837324A (en) * | 2017-04-06 | 2017-06-13 | 山东科技大学 | The fast determination method and borehole data acquisition system of rock lithology title and intensity |
| CN107633556A (en) * | 2017-09-22 | 2018-01-26 | 辽宁工程技术大学 | It is a kind of quantitatively to obtain the probabilistic method of three dimensional ore deposit geological model |
| CN107944086A (en) * | 2017-10-27 | 2018-04-20 | 江苏普旭软件信息技术有限公司 | A kind of fast modeling method based on borehole data |
| WO2018148492A1 (en) * | 2017-02-09 | 2018-08-16 | Schlumberger Technology Corporation | Geophysical deep learning |
| CN109147045A (en) * | 2018-08-08 | 2019-01-04 | 中国电建集团华东勘测设计研究院有限公司 | City section geology three-dimensional method for automatic modeling based on topological subregion |
-
2019
- 2019-06-14 CN CN201910516076.2A patent/CN110244021B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102147933A (en) * | 2010-02-05 | 2011-08-10 | 王远功 | Method for displaying stratum structure in three-dimensional scene |
| CN106558100A (en) * | 2016-10-25 | 2017-04-05 | 中国电建集团华东勘测设计研究院有限公司 | A kind of stratigraphic model method for automatic modeling based on borehole data |
| WO2018148492A1 (en) * | 2017-02-09 | 2018-08-16 | Schlumberger Technology Corporation | Geophysical deep learning |
| CN106837324A (en) * | 2017-04-06 | 2017-06-13 | 山东科技大学 | The fast determination method and borehole data acquisition system of rock lithology title and intensity |
| CN107633556A (en) * | 2017-09-22 | 2018-01-26 | 辽宁工程技术大学 | It is a kind of quantitatively to obtain the probabilistic method of three dimensional ore deposit geological model |
| CN107944086A (en) * | 2017-10-27 | 2018-04-20 | 江苏普旭软件信息技术有限公司 | A kind of fast modeling method based on borehole data |
| CN109147045A (en) * | 2018-08-08 | 2019-01-04 | 中国电建集团华东勘测设计研究院有限公司 | City section geology three-dimensional method for automatic modeling based on topological subregion |
Non-Patent Citations (3)
| Title |
|---|
| CHARLES H. GREENE等: "Three-dimensional acoustic visualization", 《DEEP-SEA RESEARCH II》 * |
| 唐彤彤等: "基于PSO － GA 的Kriging 插值法建立透地通信", 《煤炭学报》 * |
| 陈鹏等: "分层地层中的电磁波传播数值模拟", 《计算机工程》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110837668A (en) * | 2019-10-21 | 2020-02-25 | 中铁第四勘察设计院集团有限公司 | BIM software-based method and device for calculating earth volume of different soil qualities |
| CN112419500A (en) * | 2020-12-09 | 2021-02-26 | 上海申元岩土工程有限公司 | Three-dimensional geological model modeling method |
| CN112419500B (en) * | 2020-12-09 | 2023-03-14 | 上海申元岩土工程有限公司 | Three-dimensional geological model modeling method |
| CN112686779A (en) * | 2020-12-16 | 2021-04-20 | 天津泰勘工程技术咨询有限公司 | Method for dividing sequence of standard soil layers of foundation soil |
| CN112686779B (en) * | 2020-12-16 | 2022-08-02 | 天津泰勘工程技术咨询有限公司 | Method for dividing sequence of standard soil layers of foundation soil |
| CN117152369A (en) * | 2023-09-08 | 2023-12-01 | 中南大学 | Construction method, device, equipment and computer readable storage medium of frozen soil model |
| CN118015219A (en) * | 2024-04-09 | 2024-05-10 | 中国地质大学（武汉） | Geological model generation method, device and equipment based on qualitative kriging interpolation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110244021B (en) | 2021-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110244021A (en) | A kind of strata division method based on anisotropy interpolation | |
| Ringrose et al. | Reservoir model design | |
| CN106934858B (en) | Three-dimensional geological modeling method and system for scale region of mining area | |
| CN106023313B (en) | A kind of three-dimensional geological number exploration report generation method | |
| CN102222365B (en) | Method for reconstructing curved surface of complex space | |
| CN104375179B (en) | Method for looking for oil and gas reservoir based on TRAP-3D software | |
| EP3631528A1 (en) | Constructing structural models of the subsurface | |
| CN109003330A (en) | A kind of three dimensional contour line method based on basement rock boundary constraint | |
| CN106780730A (en) | The construction method and device of three-dimensional geological model | |
| CN106097445A (en) | A kind of method for drafting of novel three-dimensional stratum curved surface | |
| EP3374596B1 (en) | Fracture network triangle mesh adjustment | |
| CN104834008B (en) | Method for calculating filling degree of fractured-vuggy carbonate reservoir | |
| CN109979011A (en) | Plains region three-dimensional geological model building method based on multi-source heterogeneous data | |
| Al Rassas et al. | Application of 3d reservoir geological model on es1 formation, block nv32, shenvsi oilfield, China | |
| CN110322982A (en) | A kind of high-level waste geology treatment gun drilling method for arranging | |
| CN111859687A (en) | Mixed geological modeling method and system for depicting geological structure of uranium-bearing sand layer | |
| CN106294900A (en) | A kind of structure and methods for using them of the digital rock mass for drilling well | |
| US20180276886A1 (en) | Target object simulation using undulating surfaces | |
| CN110443890A (en) | Situ Leaching mineral deposit Stratum Modeling | |
| Majeed et al. | Constructing 3D Geological Model for Tertiary Reservoir in Khabaz Oil Field by using Petrel software. | |
| Hamdulla et al. | Building geological model for tertiary reservoir of exploration Ismail oil field, North Iraq | |
| Abdul Majeed et al. | Constructing 3D Geological Model for Tertiary Reservoir in Khabaz Oil Field by using Petrel software. | |
| CA3001127C (en) | Target object simulation using orbit propagation | |
| Fisher et al. | Rational splines and multidimensional geologic modeling | |
| Hiltunen et al. | Ground proving three seismic refraction tomography programs |
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 |