CN110244021B - Stratum layering method based on anisotropic interpolation - Google Patents
Stratum layering method based on anisotropic interpolation Download PDFInfo
- Publication number
- CN110244021B CN110244021B CN201910516076.2A CN201910516076A CN110244021B CN 110244021 B CN110244021 B CN 110244021B CN 201910516076 A CN201910516076 A CN 201910516076A CN 110244021 B CN110244021 B CN 110244021B
- Authority
- CN
- China
- Prior art keywords
- stratum
- lithologic
- attribute value
- small cubic
- isosurface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
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 the technical field of stratum layering, and provides a stratum layering method based on anisotropic interpolation, which comprises the following steps: the exploration drilling hole records the stratum lithologic name and the lithologic depth range of each stratum, a square net is established according to the exploration drilling hole structure, and the drilling hole is dispersed into a plurality of data points along the depth direction. And selecting a stratum lithology name, and performing differentiated assignment on the discrete data points in each drill hole according to whether the discrete data points in each drill hole fall within the stratum lithology depth range or not to serve as sample points. And according to the sample point attribute value, further carrying out anisotropic interpolation in the cubic network according to a specified trend to generate an attribute isosurface. And sequentially generating attribute isosurface corresponding to each lithologic name according to the operation, and dividing stratum layers according to each lithologic isosurface. The method aims at the situation that manual layering is difficult or cannot be performed due to the fact that in the actual exploration process, the number of repeated lithologies of a field is large, and the interbedded development and the like are caused. The layering result of the complex geological condition can be quickly and accurately obtained.
Description
Technical Field
The invention belongs to the technical field of stratum layering, and particularly relates to a stratum layering method based on anisotropic interpolation.
Background
Due to the fact that the soil layers are various and complex in cause and have the typical characteristics of repeated layer positions and lens body development, great difficulty is brought to manual layering in actual reconnaissance work when geological conditions of the type are met, and the space rationality and the accuracy of layering results cannot be guaranteed by a traditional single-hole layering and two-dimensional checking mode. The stratum with relatively clear geological times and better deposition rules can finish the layering work in a manual single-hole operation mode, and when the structural influence or the geological condition is complex, the manual layering has great difficulty and even the situation that the layering cannot be realized exists. For example, the lithology of red stratum generally existing in southwest and south China is mainly mudstone and sandstone interbedding, the two stratums have consistent occurrence and intervals, how to layer and connect lithologic interfaces of exploration holes in space is very difficult, stratum horizon division is performed manually by geologists, a large amount of time and energy are consumed, and even the horizon cannot be divided at all. The layering problem in actual work greatly reduces the investigation work efficiency, causes difficulty in subsequent three-dimensional geological modeling and geotechnical design, and influences the whole engineering period.
Disclosure of Invention
The invention aims to solve the problems of high layering error rate, low efficiency and incapability of layering in special situations in the traditional layering technology.
Therefore, the invention provides a stratigraphic layering method based on anisotropic interpolation, S1: carrying out geological exploration in an engineering area, and collecting lithologic interfaces of exploration drill holes, wherein information data of the lithologic interfaces comprise stratum lithologic names and stratum lithologic depth ranges;
s2: importing information data of lithologic interfaces of the drill holes, establishing a large cubic network in the whole exploration space structure, wherein the large cubic network comprises a plurality of small cubic grids, and dispersing the drill holes along the depth direction by referring to the size of the small cubic grids to obtain a plurality of discrete data points;
s3, selecting a formation lithology name, assigning the attribute value of the discrete data point falling in the formation lithology depth range to be 1 and the attribute values of the other data points to be 0 in each drill hole, and finally taking the initialized discrete data points as sample points;
s4: according to the attribute values of the sample points, an anisotropic Kriging interpolation method is adopted to obtain the attribute value of the position center point of each small cubic grid, the attribute values of the position center points of all the small cubic grids are traversed, the attribute value of the position center point of the small cubic grid with the attribute value larger than a preset intermediate value is changed into 1, the attribute value of the position center point of the small cubic grid with the attribute value smaller than the intermediate value is changed into 0, and the attribute value of the sample point is assigned to the position center point of the small cubic grid closest to the sample point;
s5: in the large cubic network, connecting the position center points of all the small cubic grids with the attribute value of 1 to generate an isosurface;
s6: repeating the steps S3 to S5, and sequentially creating isosurface of lithology of each stratum in the drill hole;
s7: and completing stratum layering of the exploration area according to the lithologic isosurface of each stratum.
Preferably, the preset intermediate value is 0.5.
Preferably, the step S4 specifically includes: obtaining an attribute value of the central point of each small cubic grid position by adopting an anisotropic Kriging interpolation method according to the sample point attribute value, wherein the calculation of each attribute value comprises the correlation degree with sample data, the correlation degree comprises the weight relation between the horizontal direction and the vertical direction, and the position distribution of surrounding samples; traversing the attribute values of the position center points of all the small cubic grids, changing the attribute value of the position center point of the small cubic grid with the attribute value larger than a preset intermediate value into 1, changing the attribute value of the position center point of the small cubic grid with the attribute value smaller than the preset intermediate value into 0, and assigning the attribute value of the sample point to the position center point of the small cubic grid closest to the sample point; and if the attribute value of the position center point of each small cubic grid is 1, the position center point of each small cubic grid is considered to belong to the stratum corresponding to the lithologic interface, and if the attribute value of the position center point of each small cubic grid is 0, the position center point of each small cubic grid is considered not to belong to the stratum corresponding to the lithologic interface.
Preferably, the step S4 further includes: performing interpolation calculation by taking a state surface of a stratum lithologic interface as a reference surface to obtain an attribute value of the position center point of each small cubic grid; and when the lithologic interface has no occurrence face, performing interpolation calculation by taking the earth surface trend face as a reference face to obtain the attribute value of the position center point of each small cubic grid.
Preferably, the formation lithology designations are plain fill, silt, silty clay, silt and silt, sandstone, mudstone, and the like.
Preferably, if the isosurfaces in two adjacent boreholes are connected, a common isosurface is formed.
Preferably, the small cubic grid location center points are connected to generate iso-surfaces, the space within the iso-surfaces forming an iso-layer.
Preferably, the step S5 is followed by: and optimizing the isosurface by using a smooth discrete interpolation method to obtain the smooth isosurface.
Preferably, the stratigraphic layering method further comprises: and adding the drilling position with the thickness of the lithologic interface being zero into the isosurface to determine the pinch-out position of the stratum.
The invention has the beneficial effects that: the stratum layering method based on anisotropic interpolation provided by the invention is characterized in that geological exploration is carried out in an engineering area, a lithology interface in a drill hole of each exploration point is collected and data is recorded, the lithology interface comprises a stratum lithology name and position information, and the lithology of the stratum disclosed by the exploration points is classified to obtain all the lithology of the stratum disclosed by an engineering field. A cubic net is constructed in an exploration area range, wherein the cubic net is a space grid set and is formed by filling a certain number of hexahedral small cubic grids with uniform sizes. Selecting a certain lithologic name, giving a specific geological attribute value, distributing the attribute value to the lithologic interface position in each drill hole, transmitting the attribute value corresponding to the lithologic interface position in each drill hole to the nearest small cubic grid, carrying out anisotropic differential interpolation on the attribute value in the small cubic grid according to an anisotropic Kriging interpolation method according to an appointed trend, attaching geological attributes to each small cubic grid after the interpolation is finished, and generating an isosurface according to the geological attribute values attached to the small cubic grid. And continuously selecting the next lithologic name, giving a specific geological attribute value, repeating the interpolation operation, sequentially creating and obtaining an isosurface corresponding to each lithologic name, and dividing the stratigraphic layering code according to each lithologic isosurface. The method aims at the situation that manual layering is difficult or impossible due to more repeated lithologies of the field, interbedded development (red layer) and the like in the actual exploration process. The layering result of the complex geological condition can be quickly and accurately obtained, and the follow-up three-dimensional geological modeling and analysis work can be served.
The present invention will be described in further detail below with reference to the accompanying drawings.
Drawings
FIG. 1 is a schematic flow chart of the method for stratigraphic layering based on anisotropic interpolation of the present invention;
FIG. 2 is a plot of a production face in a reference trend for the anisotropic interpolation based stratigraphic layering method of the present invention;
FIG. 3 is a reference trend surface plot for the anisotropic interpolation based stratigraphic layering method of the present invention;
FIG. 4 is an interpolated isosurface map without a reference trend surface for the anisotropic interpolation based stratigraphic layering method of the present invention;
FIG. 5 is an interpolated isosurface map with a reference trend surface for the anisotropic interpolation based stratigraphic layering method of the present invention;
FIG. 6 is a distribution ranking diagram of a stratigraphic lithology interface of the stratigraphic layering method based on anisotropic interpolation of the present invention;
FIG. 7 is a single interface layered reference model of the present invention based on the anisotropic interpolation stratum layering method;
FIG. 8 is a plurality of interface layering reference models of the formation layering method based on anisotropic interpolation.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the description of the present invention, it is to be understood that the terms "center", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", and the like indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplicity of description, and do not indicate or imply that the referenced devices or elements must have a particular orientation, be constructed and operated in a particular orientation, and thus, are not to be construed as limiting the present invention.
The terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature; in the description of the present invention, "a plurality" means two or more unless otherwise specified.
The embodiment of the invention provides a stratigraphic layering method based on anisotropic interpolation, which comprises the following steps of S1: geological exploration is carried out in an engineering area, lithology interfaces of all exploration drill holes are collected, and information data of the lithology interfaces comprise stratum lithology names and stratum lithology depth ranges. Drilling holes in an exploration area, then collecting geological data in the holes, wherein each hole is provided with a plurality of layers, each layer corresponds to a lithology interface, each lithology interface is composed of a plurality of discrete point grids with the same geological property, and each lithology interface is provided with a stratum lithology name and a stratum lithology depth range. Each lithology interface corresponds to a producing face.
S2: and importing information data of lithologic interfaces of the drill holes, establishing a square net in the whole exploration space structure, wherein the large cubic net comprises a plurality of small cubic grids, and discretizing the drill holes along the depth direction by referring to the size of the small cubic grids to obtain a plurality of discrete data points. It follows that each discrete data point corresponds to a small cubic grid. The lithologic interfaces of all the drill holes correspond to a depth position, but the same lithologic interfaces of different drill holes have different depths and different thicknesses, so that the names need to be classified and taken, the initialization of stratum interface data is realized, and the establishment of an equivalent layer is facilitated later. As shown in tables 1 and 2, table 1 is data of a formation log in a certain borehole, and table 2 is a table initialized with respect to data formed by naming all borehole formation data after being imported.
TABLE 1 initial layering table for a certain borehole stratum in an engineering area
Depth of bottom of layer | Stratum code number | Stratum era | Lithology name |
1.5 | 1 | Q_0_ml | Plain filling |
8 | 3 | Q_h_al+1 | Powdery clay |
11.2 | 5 | Q_p_al | Silt |
16.8 | 3 | Q_h_al+1 | Powdery clay |
25 | 5 | Q_p_al | Silt |
TABLE 2 initial lithology stratigraphic table of exploration points
Belonging to the engineering | Stratum generation | Stratum era | Rock and soil | Name of rock and soil | Boundary of China | Is a |
Lithology modeling | ||||||
1 | Q_0_ml | Soil layer | Plain filling | New world of life | Fourth series | |
Lithology modeling | 2 | Q_h_1 | Soil layer | Sludge | New world of life | Fourth |
Lithology modeling | ||||||
3 | Q_h_al+1 | Soil layer | Powdery clay | New world of life | Fourth series | |
Lithology modeling | 4 | Q_p_al | Soil layer | Silt | New world of life | Fourth |
Lithology modeling | ||||||
5 | Q_p_al | Soil layer | Silt | New world of life | Fourth series |
S3: and selecting a formation lithology name, assigning the attribute value of the discrete data point falling in the formation lithology depth range to be 1, assigning the attribute values of the other data points to be 0 in each drill hole, and finally taking the initialized discrete data points as sample points. And selecting a stratum lithology name, wherein the stratum lithology name corresponds to a stratum lithology depth range, then assigning the attribute values of the discrete data points falling in the stratum lithology depth range in other boreholes to be 1, and assigning the attribute values of the other data points to be 0, so that all the discrete data points with the same stratum lithology depth range of the stratum lithology name can be obtained.
S4: and according to the attribute values of the sample points, obtaining the attribute value of the position central point of each small cubic grid by adopting an anisotropic Kriging interpolation method, traversing the attribute values of the position central points of all the small cubic grids, changing the attribute value of the position central point of the small cubic grid with the attribute value larger than a preset intermediate value into 1, changing the attribute value of the position central point of the small cubic grid with the attribute value smaller than the intermediate value into 0, and assigning the attribute value of the sample point to the position central point of the small cubic grid closest to the sample point. Therefore, points around the discrete point are also given an attribute value of 1 or 0 according to the Kriging (Kriging) interpolation method, the attribute value of 1 indicates that the discrete point and the reference sample point are the same formation lithology name, and the attribute value of 0 does not indicate that the discrete point and the reference sample point are the same formation lithology name.
S5: in the large cubic network, the position center points of all the small cubic grids with the attribute value of 1 are connected to generate an isosurface. When the discrete points (namely geological points) in the large cubic network are connected with the surrounding discrete points with the same name, the smooth transition is realized by adopting a different Kriging interpolation method.
S6: repeating the steps S3 to S5, and sequentially creating isosurface of lithology of each stratum in the drill hole;
s7: and completing stratum layering of the exploration area according to the lithologic isosurface of each stratum.
Specifically, as shown in fig. 7, discrete data points with the same lithology name of a formation are obtained by using the lithology name and the lithology depth range of a borehole and by using a Kriging (Kriging) interpolation algorithm, that is, a lithology name of a formation is selected, in each borehole, the attribute values of the discrete data points falling within the lithology depth range of the formation are assigned to 1, the attribute values of the remaining data points are assigned to 0, and finally the initialized discrete data points are used as sample points. The sample points include the inside and outside of the borehole, and the isosurface corresponding to the formation lithology name can be formed by connecting the sample points. And circulating in this way, connecting the lithologic names of each stratum into respective corresponding isosurface, and thus completing stratum layering of the exploration area.
Calculating the attributes corresponding to all the small cubic grids in the space range according to the lithological interface, mapping the attributes of the position points of the lithological interface into the small cubic grids, and generating a plurality of equivalent spaces to obtain the approximate shape of the stratum model.
As shown in fig. 8, according to the association relationship of each layer and the upper and lower position relationship of the stratum, the overall layering of all the strata in the borehole is realized, and the layering reference result is obtained. All lithologic interfaces form corresponding isosurface respectively, so that the geology of the exploration area can be presented in a stratigraphic layering mode.
The invention analyzes and researches the distribution characteristics of the stratum or the mineral deposit under natural conditions, finds that the stratum or the mineral deposit is generally distributed in a strip shape or a block shape, or has any trend surface such as an underground aquifer, and can determine a main attitude direction according to the trend rule. If the simulation object is distributed along the occurrence, the correlation between the attributes of the object in the direction is obvious, and the other directions are opposite; and according to the attitude surface or the earth surface trend surface, the stratum data in different directions have different influence weights on the unknown points, namely anisotropic interpolation. In a specific implementation scenario, as shown in fig. 1, the geological exploration point data is imported, the formation lithology interface data revealed by the borehole is obtained, and then the imported borehole formation lithology interface data is automatically processed. When a stratum lithology name is selected, the attribute values of discrete data points falling in a stratum lithology depth range corresponding to the stratum lithology name are assigned to be 1, the attribute values of the other data points are assigned to be 0, the initialized discrete data points are used as sample points, and finally the position center points of all small cubic grids with the attribute values of 1 are connected to generate an isosurface. As shown in fig. 2 to 5, when the attribute values of the microcube grids around the sample points are assigned, the attribute values of the central points of the positions of the microcube grids are obtained by interpolation calculation with the occurrence surface of the stratigraphic lithology interface as a reference surface.
Preferably, the preset intermediate value is 0.5. This value is empirical data.
Preferably, the step S3 further includes: and if the lithology interface with the same name and position attributes as those of the lithology interface of the adjacent hole is not found in the drill hole, setting the thickness of the corresponding lithology interface of the drill hole as zero. It can be seen from this that, in general, the number of layers of the lithologic interface in each borehole is the same, and only the thickness is different, but in special cases, such as pinch-out, the thickness of the lithologic interface of this layer is defined as zero, that is, when the top layer of the lithologic interface with the thickness of zero and the stratum are overlapped and the lithologic interface of the stratum with the thickness of zero is connected with the lithologic interface of the stratum with the same position attribute and the thickness of zero, the connection is performed by the kriging interpolation algorithm in the same way, but the top surface and the bottom surface of the lithologic interface sample of the stratum with the thickness of zero are coplanar.
Preferably, the lithologic names of the stratum are plain fill, silt, silty clay, silt and silty sand, sandstone, mudstone and the like. As shown in fig. 6, it can be known that the lithologic names of the strata in the borehole in this case include three types, namely plain fill, silty clay, silt, etc., and in other cases, there may be others, which are not listed here.
Preferably, the step S4 specifically includes: obtaining an attribute value of the central point of each small cubic grid position by adopting an anisotropic Kriging interpolation method according to the sample point attribute value, wherein the calculation of each attribute value comprises the correlation degree with sample data, the correlation degree comprises the weight relation between the horizontal direction and the vertical direction, and the position distribution of surrounding samples; traversing the attribute values of the position center points of all the small cubic grids, changing the attribute value of the position center point of the small cubic grid with the attribute value larger than a preset intermediate value into 1, changing the attribute value of the position center point of the small cubic grid with the attribute value smaller than the preset intermediate value into 0, and assigning the attribute value of the sample point to the position center point of the small cubic grid closest to the sample point; and if the attribute value of the position center point of each small cubic grid is 1, the position center point of each small cubic grid is considered to belong to the stratum corresponding to the lithologic interface, and if the attribute value of the position center point of each small cubic grid is 0, the position center point of each small cubic grid is considered not to belong to the stratum corresponding to the lithologic interface.
Therefore, the specific setting method of the longitudinal correlation value is as follows: firstly, according to the minimum thickness of all stratums, determining the position center point of each small cubic grid and the lithologic interface corresponding to the small cubic grid as 0, then, according to the anisotropic Kriging interpolation method, determining the position center point of each small cubic grid and the lithologic interface corresponding to the small cubic grid as the lithologic interface corresponding to the small cubic grid, wherein the longitudinal size of the grid can be used as a longitudinal correlation value, the average value of the position attributes of all discrete points in the large cubic grid is used as the longitudinal correlation value, the attribute values of the surrounding discrete points with the same attribute value of the position center point of each small cubic grid are compared with the longitudinal attribute correlation value of the small cubic grid, if the attribute values are larger than or equal to the longitudinal attribute correlation value, the position center point of each small cubic grid is considered to belong to the lithologic interface corresponding to the small cubic grid, namely, the property value of the position center point of the small cubic grid is considered to be 1, and if the attribute values of the position center point of each small cubic grid and the lithologic interface corresponding to the small cubic grid are considered to be smaller than the longitudinal correlation value The connection of the surrounding discrete points generates an iso-layer. Therefore, the position attributes are obtained by traversing all the small cubic grids and each small cubic grid position according to the anisotropic Kriging interpolation method, and the calculation of each position attribute comprises the association degree of sample data, the set horizontal and longitudinal weight relationship, the position distribution of surrounding samples and the like. The attribute of each grid point is obtained by an anisotropic Kriging interpolation algorithm, and each position in the whole cubic network space has an attribute value by adding original sample data. Each attribute value is compared with the middle value, and if the attribute value is larger than the middle value, the point is considered to belong to the stratum, and the attribute is modified to be 1; points less than the intermediate value are considered to be points not belonging to the formation, and the attribute is modified to 0. In a large cubic network, an isosurface algorithm is created according to a small cubic grid to layer the stratum, and each isosurface obtained finally is the stratum layering model to be created.
Preferably, the step S4 further includes: performing interpolation calculation by taking a state surface of a stratum lithologic interface as a reference surface to obtain an attribute value of the position center point of each small cubic grid; and when the lithologic interface has no occurrence face, performing interpolation calculation by taking the earth surface trend face as a reference face to obtain the attribute value of the position center point of each small cubic grid. And (3) taking the attitude surface as a reference surface to interpolate attribute values of other geological points (namely discrete points) around one discrete point.
Preferably, the contour surfaces in two adjacent boreholes form a common contour surface if they are connected. Without an intersection, the top and bottom surface contour layers will also intersect as the distance of extension increases.
Preferably, a correlation model of Kriging interpolation is constructed, and the correlation degree of the correlation model along the direction of the surface trend surface is stronger than the correlation degree along the direction of the surface trend surface. The key of anisotropic interpolation is to construct a correlation model (variogram) of Kriging interpolation, wherein the correlation model has strong correlation degree and large correlation distance along the direction of a trend surface, the normal correlation degree of the trend surface is weak, and the correlation distance is small, so that the interpolation range of the stratum interface position attribute is in an area near the trend surface.
In a preferred embodiment, the stratigraphic layering method further comprises: the drill hole location with zero thickness of the lithologic interface is added to the iso-layer. The drill hole position with the thickness of 0, namely the pinch-out drill hole, is also added into the sample, so that the model automatically adjusts the pinch-out position in the stratum. Pinch-out is a geological phenomenon of the distribution of geologic horizons, present in sedimentary geology, which is the term for geological exploration. The soil layer is extremely unevenly distributed, and a defect phenomenon exists in the survey range of people. The soil layer was seen to fade away on the geological profile and the thickness appeared as a sharp pattern ending sharply.
The invention has the beneficial effects that: the stratum layering method based on anisotropic interpolation provided by the invention is characterized in that geological exploration is carried out in an engineering area, a lithology interface in a drill hole of each exploration point is collected and data is recorded, the lithology interface comprises a stratum lithology name and position information, and the lithology of the stratum disclosed by the exploration points is classified to obtain all the lithology of the stratum disclosed by an engineering field. A cubic net is constructed in an exploration area range, wherein the cubic net is a space grid set and is formed by filling a certain number of hexahedral small cubic grids with uniform sizes. Selecting a certain lithologic name, giving a specific geological attribute value, distributing the attribute value to the lithologic interface position in each drill hole, transmitting the attribute value corresponding to the lithologic interface position in each drill hole to the nearest small cubic grid, carrying out anisotropic differential interpolation on the attribute value in the small cubic grid according to an anisotropic Kriging interpolation method according to an appointed trend, attaching geological attributes to each small cubic grid after the interpolation is finished, and generating an isosurface according to the geological attribute values attached to the small cubic grid. And continuously selecting the next lithologic name, giving a specific geological attribute value, repeating the interpolation operation, sequentially creating and obtaining an isosurface corresponding to each lithologic name, and dividing the stratigraphic layering code according to each lithologic isosurface. The method aims at the situation that manual layering is difficult or impossible due to more repeated lithologies of the field, interbedded development (red layer) and the like in the actual exploration process. The layering result of the complex geological condition can be quickly and accurately obtained, and the follow-up three-dimensional geological modeling and analysis work can be served.
The above examples are merely illustrative of the present invention and should not be construed as limiting the scope of the invention, which is intended to be covered by the claims and any design similar or equivalent to the scope of the invention.
Claims (7)
1. A stratum layering method based on anisotropic interpolation is characterized in that:
s1: performing geological exploration in an engineering area, and acquiring information data of lithologic interfaces of exploration drilling holes, wherein the information data of the lithologic interfaces comprise stratum lithologic names and stratum lithologic depth ranges;
s2: importing information data of lithologic interfaces of the drill holes, constructing a large cubic network in the whole exploration space, wherein the large cubic network comprises a plurality of small cubic grids, and dispersing the drill holes along the depth direction by referring to the size of the small cubic grids to obtain a plurality of discrete data points;
s3, selecting a formation lithology name, assigning the attribute value of the discrete data point falling in the formation lithology depth range to be 1 and the attribute values of the other data points to be 0 in each drill hole, and finally taking the initialized discrete data points as sample points;
s4: obtaining an attribute value of the central point of each small cubic grid position by adopting an anisotropic Kriging interpolation method according to the attribute value of the sample point, wherein the calculation of each attribute value comprises the correlation degree with sample data, the correlation degree comprises the weight relation between the horizontal direction and the vertical direction, and the position distribution of surrounding samples; traversing the attribute values of the position center points of all the small cubic grids, changing the attribute value of the position center point of the small cubic grid with the attribute value larger than a preset intermediate value into 1, changing the attribute value of the position center point of the small cubic grid with the attribute value smaller than the intermediate value into 0, and assigning the attribute value of the sample point to the position center point of the small cubic grid closest to the sample point; if the attribute value of the position center point of each small cubic grid is 1, the position center point of each small cubic grid is considered to belong to the stratum corresponding to the lithologic interface, and if the attribute value of the position center point of each small cubic grid is 0, the position center point of each small cubic grid is considered not to belong to the stratum corresponding to the lithologic interface;
s5: in the large cubic network, connecting the position center points of all the small cubic grids with the attribute value of 1 to generate an isosurface;
s6: repeating the steps S3 to S5, and sequentially creating isosurface of lithology of each stratum in the drill hole;
s7: and completing stratum layering of the exploration area according to the lithologic isosurface of each stratum.
2. The method of claim 1, wherein the method comprises: the preset intermediate value is 0.5.
3. The method for stratigraphic layering based on anisotropic interpolation of claim 1, wherein said step S4 further comprises: performing interpolation calculation by taking a state surface of a stratum lithologic interface as a reference surface to obtain an attribute value of the position center point of each small cubic grid; and when the lithologic interface has no occurrence face, performing interpolation calculation by taking the earth surface trend face as a reference face to obtain the attribute value of the position center point of each small cubic grid.
4. The method of claim 1, wherein the method comprises: the stratum lithology names comprise plain filling soil, silt, silty clay, silty soil, silty sand, sandstone and mudstone.
5. The method of claim 1, wherein the method comprises: and connecting the position center points of the small cubic grids to generate an isosurface, wherein the space in the isosurface forms an isolayer.
6. The method for stratigraphic layering based on anisotropic interpolation of claim 1, wherein said step S5 is followed by further comprising: and optimizing the isosurface by using a smooth discrete interpolation method to obtain the smooth isosurface.
7. The method of claim 1, further comprising: and adding the drilling position with the thickness of the lithologic interface being zero into the isosurface to determine the pinch-out position of the stratum.
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 CN110244021A (en) | 2019-09-17 |
CN110244021B true 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) |
Families Citing this family (4)
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 |
CN112419500B (en) * | 2020-12-09 | 2023-03-14 | 上海申元岩土工程有限公司 | Three-dimensional geological model modeling method |
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 |
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 |
---|
Three-dimensional acoustic visualization;Charles H. Greene等;《Deep-Sea Research II》;19981231;全文 * |
分层地层中的电磁波传播数值模拟;陈鹏等;《计算机工程》;20100831;全文 * |
基于PSO - GA 的Kriging 插值法建立透地通信;唐彤彤等;《煤炭学报》;20161130;第36卷(第15期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110244021A (en) | 2019-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110244021B (en) | Stratum layering method based on anisotropic interpolation | |
USRE49507E1 (en) | Faulted geological structures having unconformities | |
EP3217192B1 (en) | Method for looking for oil and gas reservoir on the basis of trap-3d software | |
US10228478B2 (en) | Method of searching for oil-gas reservoir based on trap-3D software | |
US11209560B2 (en) | Assignment of systems tracts | |
Radwan | Three-dimensional gas property geological modeling and simulation | |
US10571601B2 (en) | Geologic fault representation for flow simulation | |
US20150113379A1 (en) | Representation of unstructured grids | |
WO2016069706A1 (en) | Generation of structural elements for subsurface formation using stratigraphic implicit function | |
CN105651676A (en) | Reservoir heterogeneity characterization method under regular development well pattern of horizontal wells | |
CN103077558A (en) | Method for establishing large-sized karst cave reservoir body distribution model of fractured-vuggy carbonate rock reservoir | |
CN110056346B (en) | Oil reservoir three-dimensional original water saturation simulation method based on trend change function | |
CN109541685B (en) | River channel sand body identification method | |
US20220163692A1 (en) | Modeling and simulating faults in subterranean formations | |
CN111768503B (en) | Sea sand resource amount estimation method based on three-dimensional geological model | |
Al-Rubaye et al. | A high resolution 3D geomodel for giant carbonate reservoir-A field case study from an Iraqi oil field | |
EP3785233B1 (en) | Fault radiation based grid compartmentalization | |
Cox et al. | Integrated modeling for optimum management of a giant gas condensate reservoir, Jurassic eolian Nugget Sandstone, Anschutz Ranch East field, Utah overthrust (USA) | |
Metwalli et al. | Reservoir petrophysical modeling and risk analysis in reserve estimation; a case study from Qasr Field, North Western Desert, Egypt | |
Hamdulla et al. | Building geological model for tertiary reservoir of exploration Ismail oil field, North Iraq | |
Podnebesnykh et al. | New approach to the evaluation of the structure of initial reserves in Ozhginskoe gas-oil field | |
CN106324668B (en) | A kind of thin seismic reservoir the Forward Modelings based on double change Review of geologic model building techniques | |
Portilla-Lazo et al. | Facies Simulation Model Construction Implementing Lithological Logging For The Sacha Field Formations | |
Shodmonovna | Structural 3D Modeling and Creation of a ThreeDimensional Grid | |
Ma et al. | Constructing 3D Model Framework and Change of Support in Data Mapping |
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 |