CN110610539A - Stratum curved surface construction method, device, equipment and storage medium - Google Patents
Stratum curved surface construction method, device, equipment and storage medium Download PDFInfo
- Publication number
- CN110610539A CN110610539A CN201910745004.5A CN201910745004A CN110610539A CN 110610539 A CN110610539 A CN 110610539A CN 201910745004 A CN201910745004 A CN 201910745004A CN 110610539 A CN110610539 A CN 110610539A
- Authority
- CN
- China
- Prior art keywords
- drilling data
- curved surface
- data
- point
- stratum
- 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.)
- Pending
Links
- 238000010276 construction Methods 0.000 title claims abstract description 57
- 238000005553 drilling Methods 0.000 claims abstract description 157
- 238000000034 method Methods 0.000 claims abstract description 35
- 230000001788 irregular Effects 0.000 claims abstract description 33
- 238000007781 pre-processing Methods 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 48
- 238000006073 displacement reaction Methods 0.000 claims description 46
- 230000006870 function Effects 0.000 claims description 42
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 230000010354 integration Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 238000003786 synthesis reaction Methods 0.000 claims description 5
- 230000001131 transforming effect Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims 2
- 230000000694 effects Effects 0.000 abstract description 5
- 230000003993 interaction Effects 0.000 abstract description 4
- 238000013461 design Methods 0.000 abstract description 3
- 238000012800 visualization Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000002452 interceptive effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a method for constructing a stratum curved surface, which comprises the following steps: acquiring drilling data of a stratum curved surface to be constructed, and constructing an irregular triangular net based on the acquired drilling data; the number of the drilling data is multiple; comprehensively preprocessing a plurality of drilling data based on an irregular triangulation network to obtain a preprocessing result; the preprocessing result comprises drilling data divided into n levels; performing rarefaction treatment on the drilling data in the pretreatment result to obtain the drilling data after rarefaction treatment; wherein, the drilling data after rarefaction processing participates in the construction of the stratum curved surface; the drilling data after rarefaction treatment is fitted by using the curved surface spline function, and a stratum curved surface is constructed, so that the interaction function of the system is enhanced, the humanized design of the system is reflected, and the three-dimensional visualization effect of the geological interface is more reasonable and accurate.
Description
Technical Field
The disclosure relates to the field of geological engineering, and in particular to a method, a device, equipment and a storage medium for constructing a stratum curved surface.
Background
Due to the complex geological structure, large amount of geological information and uneven data distribution of large-scale engineering, the method is particularly suitable for the engineering fields of water conservancy, hydropower, mines and the like. Two-dimensional and static geological engineering data are poor in intuition, cannot fully reveal spatial change rules, cannot meet the spatial analysis requirements of engineering technicians more and more, and bring great difficulty to engineering construction, mining resources and other aspects. With the development of scientific calculation visualization technology and geological information computer simulation technology, a new way is provided for establishing a very vivid three-dimensional graphic effect and constructing a three-dimensional geological model based on computable information. Three-dimensional geological modeling and analysis has become one of the leading edges and hot spots of research in many fields such as current engineering geology, mathematical geology and computer science.
At present, the general surface modeling method is suitable for geological interface interpolation processing based on a drilling data set. However, the obtained curved surface cannot be matched with other data in local areas, and the actual effect is limited. A visual curved surface modeling method developed in recent years and a man-machine interaction technology support fusion of multi-source geological data interpolation and simplification of a geological curved surface modeling process. On the premise of ensuring the smoothness of the curved surface, the method supports the local adjustment of the spatial distribution form of the curved surface so as to fit various geological observation data, such as elevation, attitude, fault and the like. The method can fully utilize various geological information collected, but the accuracy of the curved surface is low.
Disclosure of Invention
In view of this, the present disclosure provides a method for constructing a formation curved surface, including:
acquiring drilling data of a stratum curved surface to be constructed, and constructing an irregular triangular net based on the acquired drilling data; the number of the drilling data is multiple;
comprehensively preprocessing the plurality of drilling data based on the irregular triangulation network to obtain a preprocessing result; wherein the pre-processing result comprises drilling data divided into n levels;
performing rarefaction treatment on the drilling data in the preprocessing result to obtain the drilling data after rarefaction treatment; wherein the drilling data after the rarefaction treatment participates in the construction of the stratum curved surface;
and fitting the drilling data after the thinning treatment by using a curved surface spline function to construct and obtain the stratum curved surface.
In one possible implementation, the constructing an irregular triangulation network based on the acquired drilling data includes:
taking the obtained coordinates of the drilling data as discrete points, selecting any one point from the discrete points as a first point, searching the discrete point closest to the first point from the discrete points as a second point, and taking a connecting line between the first point and the second point as an initial baseline;
selecting a third point from the plurality of discrete points by applying the Delaunay rule on the right side of the initial baseline, and sequentially connecting the first point, the second point and the third point to generate a corresponding triangle;
selecting the other two sides except the initial baseline from the generated triangles as the reselected initial baselines respectively to generate corresponding triangles until the right side of the reselected initial baseline has no discrete point;
wherein the generated triangles are constructed into the irregular triangulation network.
In one possible implementation, the comprehensive preprocessing of the plurality of borehole data based on the irregular triangulation network includes:
acquiring three side lengths of each triangle in the irregular triangulation network;
when the three side lengths of the currently acquired triangle are all smaller than a preset threshold value, integrating three vertexes of the currently acquired triangle into a point, and taking the point obtained after integration as integrated drilling data;
the preset threshold represents the resolution of the stratum curved surface to be constructed under each level;
the grade of the stratum curved surface is obtained by segmenting through a Google tile map algorithm, the grade of the stratum curved surface is n grade, and the value range of n is as follows: n is more than 0 and less than or equal to 22.
In one possible implementation, the drill hole data in the preprocessing result is subjected to a thinning process, which includes:
converting the longitude and latitude coordinates of the drilling data divided into n levels into Google pixel indexes by using an application programming interface of the Google map;
calculating upper left corner coordinates (x1, y1) and lower right corner coordinates (x2, y2) of the mth level lower borehole data using the Google pixel index;
generating a sparse grid of p rows and q columns of borehole data at level n from the upper left corner coordinates (x1, y1) and the lower right corner coordinates (x2, y 2);
traversing the drilling data in each thinning grid, and when a plurality of drilling data exist in the currently acquired thinning grid, keeping the drilling data closest to the center of the currently acquired thinning grid to obtain the drilling data after thinning processing;
wherein, the value range of m is as follows: m is more than or equal to 0 and less than or equal to n-4.
In a possible implementation manner, fitting the drilling data after the rarefaction process by using a surface spline function to construct and obtain the formation surface includes:
reading a displacement formula deduced according to the plane elasticity theory:
wherein A isiRepresenting the vertical displacement set on the i point of the stratum curved surface; zA(r) is characterized inSetting vertical displacement A on i point of the stratum curved surfaceiAnd limiting the displacement to 0 except the radius R from the point iiCausing a displacement at r from point i; i is any drilling data participating in the construction of the drilling data of the stratum curved surface;
wherein r isi 2＝(x-xi)2+(y-yi)2R is ri, x, y are independent variables;
and overlapping the drilling data after the rarefaction treatment by using the displacement formula to obtain the curved surface spline function:
wherein N is the number of the drilling data after the thinning treatment;
reading the established measured elevation equation:
solving the actually measured elevation equation based on the drilling data after rarefaction treatment to obtain the coefficient in the curved surface spline function;
and fitting the drilling data after the thinning treatment by using the curve spline function after the coefficient is determined, and constructing to obtain the stratum curved surface.
In one possible implementation, reading the established measured elevation equation, and solving the measured elevation equation based on the drilling data after rarefaction processing includes:
transforming the measured elevation equation into a matrix form:
AX＝B；
and obtaining a linear algebraic equation set with the coefficient matrix being a symmetric positive definite matrix, and solving the linear algebraic equation set by using the square root of the solved symmetric positive definite equation set to obtain the coefficient in the curved surface spline function.
According to another aspect of the present disclosure, there is also provided a formation curved surface construction apparatus, including a triangulation network construction module, a data integration module, a data rarefying module and a surface construction module;
the triangular net constructing module is configured to acquire drilling data of a stratum curved surface to be constructed and construct an irregular triangular net based on the acquired drilling data; the number of the drilling data is multiple;
the data synthesis module is configured to comprehensively preprocess the plurality of drilling data based on the irregular triangulation network to obtain a preprocessing result; wherein the pre-processing result comprises drilling data divided into n levels;
the data rarefying module is configured to rarefy the drilling data in the preprocessing result to obtain the drilling data after rarefying; wherein the drilling data after the rarefaction treatment participates in the construction of the stratum curved surface;
and the earth surface construction module is configured to fit the drilling data after the rarefaction treatment by using a surface spline function to construct and obtain the stratum curved surface.
In one possible implementation manner, the earth surface construction module comprises a displacement formula reading unit, a drilling data superposition unit, an equation reading solving unit and a data fitting unit;
the displacement formula reading unit is configured to read a displacement formula derived from a plane elasticity theory:
wherein A isiRepresenting the vertical displacement set on the i point of the stratum curved surface; zA(r) characterizing setting a vertical displacement A on the i point of the formation surfaceiAnd limiting the displacement to 0 except the radius R from the point iiCausing a displacement at r from point i; i is any drilling data participating in the construction of the drilling data of the stratum curved surface;
wherein r isi 2＝(x-xi)2+(y-yi)2R is ri, x, y areAn independent variable;
the drilling data superposition unit is configured to superpose the drilling data after rarefaction processing by using the displacement formula to obtain the curved surface spline function:
wherein N is the number of the drilling data after the thinning treatment;
the equation reading solving unit is configured to read the established measured elevation equation:
wherein h isiIn order to measure the elevation,
the equation reading solving unit is further configured to solve the actually measured elevation equation based on the drilling data after rarefaction processing to obtain coefficients in the curved surface spline function;
and the data fitting unit is configured to fit the drilling data after the rarefaction treatment by using the curve spline function after the coefficient is determined, so as to construct and obtain the stratum curved surface.
According to another aspect of the present disclosure, there is provided a formation curved surface constructing apparatus including:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to execute the executable instructions to implement any of the methods described above.
According to another aspect of the present disclosure, there is provided a non-transitory computer readable storage medium having computer program instructions stored thereon, wherein the computer program instructions, when executed by a processor, implement the method of any of the preceding.
According to the stratum curved surface construction method, initial drilling data are read, an irregular triangular net is constructed based on the obtained drilling data, comprehensive pretreatment is conducted on a plurality of drilling data based on the irregular triangular net to obtain a pretreatment result, thinning treatment is conducted on the drilling data in the pretreatment result to obtain the drilling data after thinning treatment, the drilling data after thinning treatment are fitted through a curved surface spline function, a stratum curved surface is constructed, dynamic rapid interpolation can be achieved on massive drilling data, the stratum curved surface construction method is suitable for geological curved surface reconstruction under different scales, the interaction function of the system is enhanced, the humanized design of the system is reflected, and the three-dimensional visualization effect of a geological interface is more reasonable and accurate. The interactive interpolation control provides an effective solution for the contradiction between the scattered geological point data and the high requirement of geological analysis.
Other features and aspects of the present disclosure will become apparent from the following detailed description of exemplary embodiments, which proceeds with reference to the accompanying drawings.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate exemplary embodiments, features, and aspects of the disclosure and, together with the description, serve to explain the principles of the disclosure.
FIG. 1 is a flow chart of a method for constructing a formation surface according to an embodiment of the present disclosure;
fig. 2 is a schematic diagram of initial borehole data acquired in the formation curved surface construction method according to the embodiment of the present application;
fig. 3 is a schematic diagram of an irregular triangulation network constructed based on borehole data in the formation curved surface construction method according to the embodiment of the present application;
FIG. 4 is a schematic diagram of borehole data after synthesis in a formation surface construction method according to an embodiment of the present application;
fig. 5 is a schematic diagram of a thinning grid at the nth stage in the formation curved surface construction method according to the embodiment of the present application;
fig. 6 is a schematic diagram of borehole data after the nth-level grid is diluted in the formation curved surface construction method according to the embodiment of the present application;
FIG. 7 is a schematic diagram of borehole data after thinning in a formation surface construction method according to an embodiment of the present application;
FIG. 8 is a block diagram of an apparatus for formation surface construction according to an embodiment of the present application;
fig. 9 is a block diagram of a formation surface construction apparatus according to an embodiment of the present application.
Detailed Description
Various exemplary embodiments, features and aspects of the present disclosure will be described in detail below with reference to the accompanying drawings. In the drawings, like reference numbers can indicate functionally identical or similar elements. While the various aspects of the embodiments are presented in drawings, the drawings are not necessarily drawn to scale unless specifically indicated.
The word "exemplary" is used exclusively herein to mean "serving as an example, embodiment, or illustration. Any embodiment described herein as "exemplary" is not necessarily to be construed as preferred or advantageous over other embodiments.
Furthermore, in the following detailed description, numerous specific details are set forth in order to provide a better understanding of the present disclosure. It will be understood by those skilled in the art that the present disclosure may be practiced without some of these specific details. In some instances, methods, means, elements and circuits that are well known to those skilled in the art have not been described in detail so as not to obscure the present disclosure.
In order to make the technical solution of the present disclosure clearer, the present disclosure is further described in detail below with reference to the accompanying drawings and specific embodiments.
Referring to fig. 1, as a specific embodiment of the formation curved surface constructing method of the present disclosure, the method includes the following steps:
and S100, acquiring drilling data of a stratum curved surface to be constructed, and constructing an irregular triangular net based on the acquired drilling data. And S200, comprehensively preprocessing a plurality of drilling data based on the irregular triangulation network to obtain a preprocessing result. And step S300, performing rarefaction treatment on the drilling data in the preprocessing result to obtain the drilling data after rarefaction treatment. And S400, fitting the drilling data after the rarefaction treatment by using a curved surface spline function to construct a stratum curved surface.
The curved surface spline function disclosed by the invention can realize dynamic rapid interpolation on mass drilling data, is suitable for geological curved surface reconstruction under different scales, enhances the interaction function of the system, embodies the humanized design of the system and enables the constructed curved surface to be more reasonable and accurate. And the interactive interpolation control provides an effective solution for the contradiction between scattered geological point data and high geological analysis requirements.
Specifically, referring to fig. 2, the obtained original borehole data distribution, that is, the borehole data of the formation curved surface to be constructed, is shown, where the number of the obtained original borehole data is multiple.
In a possible implementation manner, in step S100, when drilling data of a curved surface of a formation to be constructed is acquired, and an irregular triangulation network is constructed based on the acquired drilling data, the following manner may be performed.
Firstly, the coordinates of the acquired drilling data are used as discrete points, any one point is selected from the discrete points to be used as a first point, the discrete point closest to the first point is searched from the discrete points to be used as a second point, and a connecting line between the first point and the second point is used as an initial baseline.
Such as: referring to fig. 2, in the discrete points of the borehole data, the j point is selected as a first point, a point closest to the j point is searched in the discrete points as a second point, the i point is the second point, and a line connecting the i point and the j point, and the i point and the j point is used as an initial baseline.
Then, a third point is selected from the plurality of discrete points by applying the Delaunay rule on the right side of the initial baseline, and the first point, the second point and the third point are connected in sequence to generate a corresponding triangle. And then, selecting the other two sides except the initial baseline from the generated triangles as the reselected initial baselines respectively, and generating the corresponding triangles until no discrete point exists on the right side of the reselected initial baselines.
For example, referring to fig. 2, in the discrete points on the right side of the initial baseline formed by connecting the i point and the j point, the Delaunay rule is used to find the point closest to the initial baseline formed by connecting the i point and the j point, and the point is used as the third point, then the n point is the third point, the n point is connected with the i point, the n point is connected with the j point, a triangle is generated, the edge formed by connecting the n point and the i point of the triangle and the edge formed by connecting the n point and the j point are respectively used as the reselected initial baseline, and a plurality of triangles are generated until there is no discrete point on the right side of the reselected initial baseline. That is, all baselines are processed through the above steps.
Thus, the generated triangles can be constructed into an irregular triangulation network, as shown in fig. 3.
Here, it should be noted that other triangulation network construction algorithms may also be used when constructing the irregular triangulation network, and are not specifically limited herein.
After the corresponding irregular triangulation network is constructed in the above manner, the plurality of drilling data may be comprehensively preprocessed based on the irregular triangulation network in step S200 to obtain a preprocessing result. Wherein the pre-processing result comprises the borehole data divided into n levels.
In a possible implementation manner, when comprehensively processing a plurality of drilling data based on the constructed irregular triangulation network, the method may include the following steps:
and acquiring three side lengths of each triangle in the irregular triangulation network. And when the three side lengths of the currently acquired triangle are all smaller than a preset threshold value, integrating three vertexes of the currently acquired triangle into one point, and taking the point obtained after integration as the drilling data after integration. It should be noted that the preset threshold value represents the resolution of the formation surface to be constructed at each level. The grade of the stratum curved surface is obtained by segmenting through a Google tile map algorithm, the grade of the stratum curved surface is n grade, and the value range of n is as follows: n is more than 0 and less than or equal to 22.
That is, in order to facilitate hierarchical display of the generated stratigraphic curved surface under different scales, the google tile map algorithm is utilized to divide the generated stratigraphic curved surface model into n levels (for example, the value of n can be 22, that is, the division is 0-22 levels), the resolution of each level is used as a threshold of the side length of the irregular triangulation network, if the three side lengths of the triangle in the triangulation network are all smaller than the threshold, three vertexes of the triangle are integrated into one point, the point is the vertical center of the triangle, the integrated drilling data is as shown in fig. 4, and the point l, the point m and the point n are integrated into one point q. Therefore, through comprehensive pretreatment on the plurality of drilling data, the drilling data of each level of the stratum curved surface model are obtained after primary screening of the drilling data is realized.
However, after the comprehensive preprocessing of the drilling data through the above steps, there still exists excessive drilling data. Referring to FIG. 4, e.g., borehole data e, f, g, etc. At this time, further thinning of the drilling data is needed, so that the drilling data of the stratum curved surface model participating in construction at each level are determined.
And when the comprehensive drilling data is further subjected to thinning treatment, a GridTile algorithm can be adopted to realize the thinning treatment. That is, in step S300, the drill hole data in the preprocessing result is subjected to thinning processing using GridTile algorithm, and the drill hole data after thinning processing is obtained.
The method specifically comprises the following steps of adopting a GridTile algorithm to carry out rarefaction treatment on the preprocessed drilling data (namely, the drilling data in the preprocessing result):
converting longitude and latitude coordinates of the drilling data divided into n levels into Google pixel indexes by using an application programming interface of the Google map;
calculating the coordinates of the upper left corner (x1, y1) and the coordinates of the lower right corner (x2, y2) of the drill hole data under the mth level by using Google pixel indexes;
generating a sparse grid of p rows and q columns of the drilling data under the nth level according to the coordinates of the upper left corner (x1, y1) and the coordinates of the lower right corner (x2, y 2);
traversing the drilling data in each rarefying grid, and when a plurality of drilling data exist in the currently acquired rarefying grid, reserving the drilling data closest to the center of the currently acquired rarefying grid to obtain the drilling data after rarefying treatment;
it should be noted that, the value range of m is: m is more than or equal to 0 and less than or equal to n-4.
For example, n may take the value of 22. Therefore, in a possible implementation manner, when the GridTile algorithm is used for rarefiing the drill holes and determining the drill hole data participating in the construction of the formation curved surface model at each of the 0-22 levels, the method specifically includes the following steps:
and converting the longitude and latitude coordinates of the 22-level drilling data into Google pixel indexes by utilizing a Google map application program interface.
The coordinates (x1, y1) of the upper left corner and the coordinates (x2, y2) of the envelope rectangle of the drill data at the nth' stage are calculated. Traversing all borehole data Google pixel indexes, x1 being the minimum lateral index value of borehole data, x2 being the maximum lateral pixel index value of borehole data, y1 being the minimum longitudinal index value of borehole data, and y2 being the maximum longitudinal index value of borehole data.
And generating an nth' lower drilling thinning grid. If the number of rows and columns of the thinning grid is p and q, the values of p and q are:
referring to fig. 5, after the envelope rectangle of the nth' stage drilling data is divided into grids of p rows and q columns, the drilling data in each grid is traversed, if a plurality of drilling data exist in the grids, only the drilling data closest to the center of the grids are retained, as shown in fig. 6, three points c, e, and g are deleted, and finally the thinning result is shown in fig. 7.
And according to the steps, sequentially thinning the drilling data of each level below the 0-18 levels. In order to ensure modeling accuracy under a large scale, the 19-22-level drilling data can be not subjected to thinning treatment.
Further, as a specific embodiment of the formation curved surface construction method of the present disclosure, after the obtained drilling data is processed in any of the above manners, step 400 may be executed, and the formation curved surface is constructed by fitting the drilling data after the thinning processing by using a curved surface spline function.
In a possible implementation manner, when fitting the drilling data after the rarefaction process by using a curved surface spline function, the following implementation manner can be used.
First, the displacement formula derived from the plane elasticity theory is read:
wherein A isiRepresenting vertical displacement set on a point i of the stratum curved surface; zA(r) characterizing the setting of vertical displacement A on the i point of the formation surfaceiAnd limiting the displacement to 0 except the radius R from the point iiCausing a displacement at r from point i; and i is any drilling data in the drilling data participating in the construction of the formation curved surface.
In addition, r isi 2＝(x-xi)2+(y-yi)2R is ri, x and y are independent variables.
Secondly, by utilizing a displacement formula, the drilling data after rarefaction treatment is superposed to obtain a curved surface spline function:
n is the number of drill data after the thinning process.
And then, reading the established actual measurement elevation equation, and solving the actual measurement elevation equation based on the drilling data after rarefaction processing so as to determine undetermined coefficients in the curved surface spline function.
Wherein the established actual measurement elevation equation is as follows:
in addition, h isiIn order to measure the elevation,
in a possible implementation manner, the actual measurement elevation equation is solved to determine the undetermined coefficient in the surface spline function, and the determination can be realized in a matrix manner.
That is, the elevation equation is transformed into a matrix:
AX＝B；
the method comprises the steps of obtaining a linear algebraic equation set with a coefficient matrix being a symmetrical positive definite matrix, solving the linear algebraic equation set by utilizing a square root of the solved symmetrical positive definite equation set to obtain a coefficient in a curved surface spline function, and fitting the drill hole data after thinning treatment by utilizing the curved surface spline function after determining the coefficient to construct a stratum curved surface.
Here, when the formation curved surface is obtained by fitting the curved surface spline function, the curved surface spline function can be regarded as deformation when an infinite flat plate is bent, and given a set of curved surface coordinates (Xi, Yi, Zi) (i ═ 1, 2, … …, n), the function of two independent variables can be approximated and estimated, so that a curved surface which better conforms to reality can be obtained. Giving a vertical displacement A at point i in the planeiWhen the displacement is 0 except the radius R from the point i, the elastic theory of the plane can be used to deduce the displacement from AiCause the distance i point to be r out (r)<R) displacement ZA(r) is
In the formula, ri 2＝(x-xi)2+(y-yi)2The radius is affected for the flat area. r ═ ri(ii) a x and y are independent variables.
The single-value surface fitting function to the discrete point data by superposing the N measured data by using the formula (1)
In the formula, AiIs the undetermined coefficient. So that x is equal to xi，y＝yiWhere Z (x, y) is equal to the measured elevation hiThe following equation is established
In the formula (I), the compound is shown in the specification,n equations are shared by N actual measurement points (3), and N undetermined coefficients A can be obtained by solving the equationsjThese values are substituted into the formula (2), and it is obvious that the formula (2) passes through the N measured surface elevation points and is smooth and continuous everywhere.
General formula hiWritten in matrix form as: AX is B. In addition
Therein is shown
X＝(A1，A2，…，An)T
B＝(h1，h2，…，hn)T
The system of equations is a linear algebraic system of equations with a coefficient matrix being a symmetric positive definite matrix, which can be solved by solving the square root of the symmetric positive definite system of equations.
It should be noted that, although the formation surface construction method is described by taking the above steps as an example, those skilled in the art will understand that the disclosure should not be limited thereto. In fact, a user can flexibly set the formation curved surface construction method according to personal preference and/or practical application scenes as long as the required functional effect is achieved.
Therefore, the constructed curved surface is more practical by improving the spline function of the curved surface, and the accuracy of the constructed curved surface of the stratum can be improved according to the construction method of the curved surface of the stratum disclosed by the embodiment of the disclosure.
In accordance with another aspect of the present disclosure, a formation surface construction apparatus 100 is also provided. Since the working principle of the formation curved surface constructing apparatus 100 according to the embodiment of the present disclosure is the same as or similar to that of the formation curved surface constructing method according to the embodiment of the present disclosure, repeated descriptions are omitted. Referring to fig. 8, the formation surface construction apparatus 100 of the embodiment of the present disclosure includes a triangulation network construction module 110, a data synthesis module 120, a data thinning module 130, and a surface construction module 140.
The triangulation network construction module 110 is configured to acquire drilling data of a formation curved surface to be constructed, and construct an irregular triangulation network based on the acquired drilling data; wherein, the number of the drilling data is a plurality.
The data synthesis module 120 is configured to perform comprehensive preprocessing on the plurality of drilling data based on the irregular triangulation network to obtain a preprocessing result; wherein the pre-processing result comprises the borehole data divided into n levels.
The data thinning module 130 is configured to perform thinning processing on the drilling data in the preprocessing result to obtain the drilling data after thinning processing; and the drilling data after the thinning treatment participates in the construction of the stratum curved surface.
And the surface construction module 140 is configured to fit the drilling data after the rarefaction treatment by using a surface spline function to construct a stratum surface.
Further, the earth surface construction module 140 in the formation construction curved surface device 100 of the present disclosure includes a displacement formula reading unit, a drilling data superposition unit, an equation reading solving unit, and a data fitting unit.
A displacement formula reading unit configured to read a displacement formula derived from the plane elasticity theory:
wherein A isiRepresenting the vertical displacement set on the i point of the stratum curved surface; zA(r) characterizing the setting of vertical displacement A on the i point of the formation surfaceiAnd limiting the displacement to 0 except the radius R from the point iiCausing a displacement at r from point i; and i is any drilling data in the drilling data participating in the construction of the formation curved surface.
In addition, r isi 2＝(x-xi)2+(y-yi)2R is ri, x and y are independent variables.
The drilling data superposition unit is configured to superpose the drilling data after rarefaction treatment by using a displacement formula to obtain a curved surface spline function:
n is the number of drill data after the thinning process.
An equation reading solving unit configured to read the established measured elevation equation:
in addition, h isiIn order to measure the elevation,
the equation reading solving unit is also configured to solve the actually measured elevation equation based on the drilling data after rarefaction processing to obtain a coefficient in a curved surface spline function;
and the data fitting unit is configured to fit the drilling data after the rarefaction treatment by using the curve spline function after the coefficient is determined, so as to construct a stratum curved surface.
Still further, based on the same inventive concept, in accordance with another aspect of the present disclosure, there is also provided a formation curved surface constructing apparatus 200. Referring to fig. 9, the stratigraphic surface construction apparatus 200 according to the embodiment of the present disclosure includes a processor 210 and a memory 220 for storing instructions executable by the processor 210. Wherein the processor 210 is configured to execute the executable instructions to implement any of the formation surface construction methods described above.
Here, it should be noted that the number of the processors 210 may be one or more. Meanwhile, in the stratigraphic surface construction equipment 200 of the embodiment of the present disclosure, an input device 230 and an output device 240 may be further included. The processor 210, the memory 220, the input device 230, and the output device 240 may be connected via a bus, or may be connected via other methods, which is not limited in detail herein.
The memory 220, which is a computer-readable storage medium, may be used to store software programs, computer-executable programs, and various modules, such as: the method for constructing the stratum curved surface of the embodiment of the disclosure corresponds to a program or a module. The processor 210 executes various functional applications and data processing of the formation surface construction apparatus 200 by executing software programs or modules stored in the memory 220.
The input device 230 may be used to receive an input number or signal. Wherein the signal may be a key signal generated in connection with user settings and function control of the device/terminal/server. The output device 240 may include a display device such as a display screen.
According to another aspect of the present disclosure, there is also provided a non-transitory computer readable storage medium having stored thereon computer program instructions, which when executed by the processor 210, implement the formation surface construction method as described in any of the preceding.
Having described embodiments of the present disclosure, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the disclosed embodiments. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The terminology used herein is chosen in order to best explain the principles of the embodiments, the practical application, or improvements made to the technology in the marketplace, or to enable others of ordinary skill in the art to understand the embodiments disclosed herein.
Claims (10)
1. A method for constructing a formation curved surface, comprising:
acquiring drilling data of a stratum curved surface to be constructed, and constructing an irregular triangular net based on the acquired drilling data; the number of the drilling data is multiple;
comprehensively preprocessing the plurality of drilling data based on the irregular triangulation network to obtain a preprocessing result; wherein the pre-processing result comprises drilling data divided into n levels;
performing rarefaction treatment on the drilling data in the preprocessing result to obtain the drilling data after rarefaction treatment; wherein the drilling data after the rarefaction treatment participates in the construction of the stratum curved surface;
and fitting the drilling data after the thinning treatment by using a curved surface spline function to construct and obtain the stratum curved surface.
2. The method of claim 1, wherein constructing an irregular triangulation network based on the acquired borehole data comprises:
taking the obtained coordinates of the drilling data as discrete points, selecting any one point from the discrete points as a first point, searching the discrete point closest to the first point from the discrete points as a second point, and taking a connecting line between the first point and the second point as an initial baseline;
selecting a third point from the plurality of discrete points by applying the Delaunay rule on the right side of the initial baseline, and sequentially connecting the first point, the second point and the third point to generate a corresponding triangle;
selecting the other two sides except the initial baseline from the generated triangles as the reselected initial baselines respectively to generate corresponding triangles until the right side of the reselected initial baseline has no discrete point;
wherein the generated triangles are constructed into the irregular triangulation network.
3. The method of claim 1, wherein comprehensively preprocessing a plurality of the borehole data based on the irregular triangulation network comprises:
acquiring three side lengths of each triangle in the irregular triangulation network;
when the three side lengths of the currently acquired triangle are all smaller than a preset threshold value, integrating three vertexes of the currently acquired triangle into a point, and taking the point obtained after integration as integrated drilling data;
the preset threshold represents the resolution of the stratum curved surface to be constructed under each level;
the grade of the stratum curved surface is obtained by segmenting through a Google tile map algorithm, the grade of the stratum curved surface is n grade, and the value range of n is as follows: n is more than 0 and less than or equal to 22.
4. The method of claim 3, wherein the rarefying of the borehole data in the pre-processing results comprises:
converting the longitude and latitude coordinates of the drilling data divided into n levels into Google pixel indexes by using an application programming interface of the Google map;
calculating upper left corner coordinates (x1, y1) and lower right corner coordinates (x2, y2) of the mth level lower borehole data using the Google pixel index;
generating a sparse grid of p rows and q columns of borehole data at level n from the upper left corner coordinates (x1, y1) and the lower right corner coordinates (x2, y 2);
traversing the drilling data in each thinning grid, and when a plurality of drilling data exist in the currently acquired thinning grid, keeping the drilling data closest to the center of the currently acquired thinning grid to obtain the drilling data after thinning processing;
wherein, the value range of m is as follows: m is more than or equal to 0 and less than or equal to n-4.
5. The method according to any one of claims 1 to 4, wherein fitting the rarefaction-processed borehole data with a surface spline function to construct the formation surface comprises:
reading a displacement formula deduced according to the plane elasticity theory:
wherein A isiRepresenting the vertical displacement set on the i point of the stratum curved surface; zA(r) characterizing setting a vertical displacement A on the i point of the formation surfaceiAnd limiting the displacement to 0 except the radius R from the point iiCausing a displacement at r from point i; i is any drilling data participating in the construction of the drilling data of the stratum curved surface;
wherein the content of the first and second substances,r is ri, x, y are independent variables;
and overlapping the drilling data after the rarefaction treatment by using the displacement formula to obtain the curved surface spline function:
wherein N is the number of the drilling data after the thinning treatment;
reading the established measured elevation equation:
wherein h isiIn order to measure the elevation,
solving the actually measured elevation equation based on the drilling data after rarefaction treatment to obtain the coefficient in the curved surface spline function;
and fitting the drilling data after the thinning treatment by using the curve spline function after the coefficient is determined, and constructing to obtain the stratum curved surface.
6. The method of claim 5, wherein reading the established measured elevation equations and solving the measured elevation equations based on the rarefied borehole data comprises:
transforming the measured elevation equation into a matrix form:
AX＝B；
and obtaining a linear algebraic equation set with the coefficient matrix being a symmetric positive definite matrix, and solving the linear algebraic equation set by using the square root of the solved symmetric positive definite equation set to obtain the coefficient in the curved surface spline function.
7. The stratum curved surface construction device is characterized by comprising a triangular net construction module, a data integration module, a data rarefying module and a stratum construction module;
the triangular net constructing module is configured to acquire drilling data of a stratum curved surface to be constructed and construct an irregular triangular net based on the acquired drilling data; the number of the drilling data is multiple;
the data synthesis module is configured to comprehensively preprocess the plurality of drilling data based on the irregular triangulation network to obtain a preprocessing result; wherein the pre-processing result comprises drilling data divided into n levels;
the data rarefying module is configured to rarefy the drilling data in the preprocessing result to obtain the drilling data after rarefying; wherein the drilling data after the rarefaction treatment participates in the construction of the stratum curved surface;
and the earth surface construction module is configured to fit the drilling data after the rarefaction treatment by using a surface spline function to construct and obtain the stratum curved surface.
8. The stratigraphic surface construction device according to claim 7, wherein the surface construction module comprises a displacement formula reading unit, a borehole data superposition unit, an equation reading solution unit and a data fitting unit;
the displacement formula reading unit is configured to read a displacement formula derived from a plane elasticity theory:
wherein A isiRepresenting the vertical displacement set on the i point of the stratum curved surface; zA(r) characterizing setting a vertical displacement A on the i point of the formation surfaceiAnd limiting the displacement to 0 except the radius R from the point iiCausing a displacement at r from point i; i is any drilling data participating in the construction of the drilling data of the stratum curved surface;
wherein the content of the first and second substances,r is ri, x, y are independent variables;
the drilling data superposition unit is configured to superpose the drilling data after rarefaction processing by using the displacement formula to obtain the curved surface spline function:
wherein N is the number of the drilling data after the thinning treatment;
the equation reading solving unit is configured to read the established measured elevation equation:
wherein h isiIn order to measure the elevation,
the equation reading solving unit is further configured to solve the actually measured elevation equation based on the drilling data after rarefaction processing to obtain coefficients in the curved surface spline function;
and the data fitting unit is configured to fit the drilling data after the rarefaction treatment by using the curve spline function after the coefficient is determined, so as to construct and obtain the stratum curved surface.
9. A formation surface construction apparatus comprising:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to carry out the method of any one of claims 1 to 6 when executing the executable instructions.
10. A non-transitory computer readable storage medium having computer program instructions stored thereon, wherein the computer program instructions, when executed by a processor, implement the method of any of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910745004.5A CN110610539A (en) | 2019-08-13 | 2019-08-13 | Stratum curved surface construction method, device, equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910745004.5A CN110610539A (en) | 2019-08-13 | 2019-08-13 | Stratum curved surface construction method, device, equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110610539A true CN110610539A (en) | 2019-12-24 |
Family
ID=68889815
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910745004.5A Pending CN110610539A (en) | 2019-08-13 | 2019-08-13 | Stratum curved surface construction method, device, equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110610539A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112037329A (en) * | 2020-09-03 | 2020-12-04 | 华东勘测设计院(福建)有限公司 | Irregular curved surface construction method |
| CN112446956A (en) * | 2021-02-01 | 2021-03-05 | 广汉前沿岩土科技有限公司 | Geological three-dimensional model construction method based on drilling data |
| CN113239136A (en) * | 2021-05-17 | 2021-08-10 | 北京车和家信息技术有限公司 | Data processing method, device, equipment and medium |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102254349A (en) * | 2011-06-30 | 2011-11-23 | 华东师范大学 | Method for constructing three-dimensional physical model of sedimentary stratum system by using drilling data |
| CN102346923A (en) * | 2010-07-30 | 2012-02-08 | 中国科学院遥感应用研究所 | Data grading organization method based on longitude and latitude grid |
| CN104572755A (en) * | 2013-10-24 | 2015-04-29 | 高德软件有限公司 | Method for creating data index, data searching method and related device |
| CN104821013A (en) * | 2015-05-11 | 2015-08-05 | 武汉大学 | Method and system for specific surface area extraction based on geodetic coordinate system digital elevation model |
| CN105184864A (en) * | 2015-08-14 | 2015-12-23 | 中国能源建设集团安徽省电力设计院有限公司 | Site stratum three-dimensional geological structure model generation method for natural foundation replacement quantities calculation |
| CN109002451A (en) * | 2017-06-07 | 2018-12-14 | 杭州海康威视系统技术有限公司 | Map datum vacuates method and device |
-
2019
- 2019-08-13 CN CN201910745004.5A patent/CN110610539A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102346923A (en) * | 2010-07-30 | 2012-02-08 | 中国科学院遥感应用研究所 | Data grading organization method based on longitude and latitude grid |
| CN102254349A (en) * | 2011-06-30 | 2011-11-23 | 华东师范大学 | Method for constructing three-dimensional physical model of sedimentary stratum system by using drilling data |
| CN104572755A (en) * | 2013-10-24 | 2015-04-29 | 高德软件有限公司 | Method for creating data index, data searching method and related device |
| CN104821013A (en) * | 2015-05-11 | 2015-08-05 | 武汉大学 | Method and system for specific surface area extraction based on geodetic coordinate system digital elevation model |
| CN105184864A (en) * | 2015-08-14 | 2015-12-23 | 中国能源建设集团安徽省电力设计院有限公司 | Site stratum three-dimensional geological structure model generation method for natural foundation replacement quantities calculation |
| CN109002451A (en) * | 2017-06-07 | 2018-12-14 | 杭州海康威视系统技术有限公司 | Map datum vacuates method and device |
Non-Patent Citations (9)
| Title |
|---|
| U011393661: "Google地图经度纬度到像素索引值的转换公式", 《BLOG.CSDN.NET/U011393661/ARTICLE/DETAILS/14055439》, 2 November 2013 (2013-11-02), pages 1 - 2 * |
| 唐泽圣 等: "《三维数据场可视化》", 31 May 2000, 清华大学出版社, pages: 196 - 200 * |
| 宋仁波等: "基于限定Delaunay三角剖分和OpenGL的三维地层可视化系统设计与实现", 《工程勘察》 * |
| 宋仁波等: "基于限定Delaunay三角剖分和OpenGL的三维地层可视化系统设计与实现", 《工程勘察》, no. 09, 1 September 2013 (2013-09-01) * |
| 方子岩 等: "《摄影测量学》", 31 October 2013, 长江出版社, pages: 364 - 368 * |
| 赵鹏大 等: "《数学地质和地质信息》", 31 July 1999, 地质出版社, pages: 69 - 77 * |
| 郭艳军等: "基于钻孔数据和交叉折剖面约束的三维地层建模方法研究", 《地理与地理信息科学》 * |
| 郭艳军等: "基于钻孔数据和交叉折剖面约束的三维地层建模方法研究", 《地理与地理信息科学》, no. 02, 15 March 2009 (2009-03-15), pages 23 - 27 * |
| 陆守一 编: "《地理信息系统》", 北京：高等教育出版社, pages: 166 - 170 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112037329A (en) * | 2020-09-03 | 2020-12-04 | 华东勘测设计院(福建)有限公司 | Irregular curved surface construction method |
| CN112037329B (en) * | 2020-09-03 | 2023-10-13 | 华东勘测设计院(福建)有限公司 | Irregular curved surface construction method |
| CN112446956A (en) * | 2021-02-01 | 2021-03-05 | 广汉前沿岩土科技有限公司 | Geological three-dimensional model construction method based on drilling data |
| CN113239136A (en) * | 2021-05-17 | 2021-08-10 | 北京车和家信息技术有限公司 | Data processing method, device, equipment and medium |
| CN113239136B (en) * | 2021-05-17 | 2023-12-19 | 北京车和家信息技术有限公司 | Data processing method, device, equipment and medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huang et al. | Explorations of the implementation of a parallel IDW interpolation algorithm in a Linux cluster-based parallel GIS | |
| WO2022007398A1 (en) | Unstructured grid flood simulation system based on gpu acceleration thechnology | |
| RU2603978C1 (en) | Creation of formation power grid to determine estimates of margins for formation | |
| EP2869096A2 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| CN110610539A (en) | Stratum curved surface construction method, device, equipment and storage medium | |
| CN104821013A (en) | Method and system for specific surface area extraction based on geodetic coordinate system digital elevation model | |
| CN105631935A (en) | BIM-based three-dimensional geological automatic modeling method | |
| US9239398B2 (en) | Method of generating and combining multiple horizons to determine a seismic horizon and its uncertainty | |
| Fadakar Alghalandis et al. | The RANSAC method for generating fracture networks from micro-seismic event data | |
| Liu et al. | Change detection of mobile LIDAR data using cloud computing | |
| US20220244424A1 (en) | Geological Grid Analysis | |
| CN116774292B (en) | Seismic wave travel time determining method, system, electronic equipment and storage medium | |
| Mojica et al. | Regularization parameter selection in the 3D gravity inversion of the basement relief using GCV: A parallel approach | |
| Berre et al. | Call for participation: Verification benchmarks for single-phase flow in three-dimensional fractured porous media | |
| CN107507179B (en) | Rock-soil mass quantitative analysis method based on GOCAD | |
| CN105931297A (en) | Data processing method applied to three-dimensional geological surface model | |
| Brunsvik et al. | Three-dimensional paganica fault morphology obtained from hypocenter clustering (L'Aquila 2009 seismic sequence, Central Italy) | |
| CN117409157A (en) | Online sectioning and graph cutting method for three-dimensional geological model | |
| Liao et al. | Data assimilation for strongly nonlinear problems by transformed ensemble Kalman filter | |
| CN115375867B (en) | Method, system, equipment and medium for calculating geothermal resource quantity by using grid model | |
| Lieberwirth et al. | Applying low budget equipment and open source software for high resolution documentation of archaeological stratigraphy and features | |
| CN110610544A (en) | Drilling data thinning method and device, equipment and storage medium | |
| CN105869209A (en) | Deformed triangular data processing method in three-dimensional geological surface model | |
| Senyurdusev et al. | Exploring the opportunities of open source data use in creation 3d procedural city models | |
| RU2812143C1 (en) | Method and device for measuring characteristics of rock column to create model of pore system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191224 |