CN117152369A - Construction method, device, equipment and computer readable storage medium of frozen soil model - Google Patents
Construction method, device, equipment and computer readable storage medium of frozen soil model Download PDFInfo
- Publication number
- CN117152369A CN117152369A CN202311156890.0A CN202311156890A CN117152369A CN 117152369 A CN117152369 A CN 117152369A CN 202311156890 A CN202311156890 A CN 202311156890A CN 117152369 A CN117152369 A CN 117152369A
- Authority
- CN
- China
- Prior art keywords
- interpolation
- grid
- frozen soil
- data
- point cloud
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 71
- 238000010276 construction Methods 0.000 title claims description 14
- 238000000034 method Methods 0.000 claims abstract description 42
- 238000005553 drilling Methods 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 20
- 230000006870 function Effects 0.000 claims description 23
- 238000010586 diagram Methods 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 13
- 238000012216 screening Methods 0.000 claims description 8
- 238000013139 quantization Methods 0.000 abstract description 6
- 238000004891 communication Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- YZMCKZRAOLZXAZ-UHFFFAOYSA-N sulfisomidine Chemical compound CC1=NC(C)=CC(NS(=O)(=O)C=2C=CC(N)=CC=2)=N1 YZMCKZRAOLZXAZ-UHFFFAOYSA-N 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 101100129500 Caenorhabditis elegans max-2 gene Proteins 0.000 description 1
- 208000025174 PANDAS Diseases 0.000 description 1
- 208000021155 Paediatric autoimmune neuropsychiatric disorders associated with streptococcal infection Diseases 0.000 description 1
- 240000000220 Panda oleosa Species 0.000 description 1
- 235000016496 Panda oleosa Nutrition 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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
Embodiments of the present application provide methods, apparatuses, devices, and computer-readable storage media for constructing a frozen soil model. The method comprises the steps of obtaining ground temperature data; setting space three-dimensional interpolation grid parameters based on the ground temperature data to obtain grid information to be interpolated; performing interpolation processing on the horizontal grids corresponding to the grid information to obtain frozen earth point cloud data; performing tetrahedral subdivision on the frozen earth point cloud data; and establishing a three-dimensional frozen soil model according to the subdivision result of the frozen soil point cloud data. In this way, the frozen soil distribution in the reservoir area can be represented from the three-dimensional angle under the condition of fully utilizing the drilling data, so that the processing flow is simplified, and meanwhile, the accuracy of data quantization is greatly improved.
Description
Technical Field
Embodiments of the present application relate to the field of data processing, and in particular, to a method, an apparatus, a device, and a computer readable storage medium for constructing a frozen soil model.
Background
Aiming at the quantification work of frozen soil in alpine regions, the research scope of the traditional method is often single, namely the frozen soil in the plane is quantified from the angle of a certain horizontal plane or a vertical section through drilling data, so that the follow-up work is carried out.
The quantization of the horizontal plane usually adopts methods such as Kerling interpolation, taylor polygon interpolation and the like; for vertical profiles, the ice layer quantization is mostly performed by means of a substantially linear interpolation. The existing survey data cannot be fully utilized by adopting the processing mode, and the quantization range is limited.
For example:
1. if frozen soil quantification of the vertical section is carried out, the arrangement form of drilling holes is greatly depended, and collinear drilling points with a certain number are often required to be searched;
2. the quantification of frozen soil for a horizontal plane often reflects only a certain height of the frozen soil distribution range.
However, stability analysis and the like for the entire engineering often cannot be performed only by a certain plane. In the two-dimensional method, if the designated section has insufficient drilling data, the ice layer distribution condition of the section cannot be directly obtained, and the accuracy of the analysis result is directly affected. The frozen soil model constructed by the existing mode cannot be matched with three-dimensional analysis, the model format of the frozen soil model also has problems, and the follow-up analysis work can be greatly influenced by the lack of three-dimensional information of frozen soil distribution.
In summary, how to fully utilize drilling data to characterize the distribution of frozen soil in a reservoir from a three-dimensional rather than planar perspective is a current need to be addressed.
Disclosure of Invention
According to the embodiment of the application, a construction scheme of a frozen soil model is provided, frozen soil is represented by a three-dimensional form, data are filled by utilizing kringing, a tetrahedron algorithm is optimized through a preset grid and side length screening, the common boundary problem of the method is solved, a three-dimensional frozen soil grid which can represent three-dimensional characteristics and can be used for subsequent numerical calculation is obtained, frozen soil distribution in a storehouse area can be represented from a three-dimensional rather than planar angle under the condition of fully utilizing drilling data, and the accuracy of data quantization is greatly improved while the processing flow is simplified.
In a first aspect of the application, a method for constructing a frozen soil model is provided. The method comprises the following steps:
acquiring ground temperature data;
setting space three-dimensional interpolation grid parameters based on the ground temperature data to obtain grid information to be interpolated;
performing interpolation processing on the horizontal grids corresponding to the grid information to obtain frozen earth point cloud data;
performing tetrahedral subdivision on the frozen earth point cloud data; and establishing a three-dimensional frozen soil model according to the subdivision result of the frozen soil point cloud data.
Further, the ground temperature data includes borehole horizontal position, borehole depth, and borehole temperature information.
Further, setting a spatial three-dimensional interpolation grid parameter based on the ground temperature data, and obtaining grid information to be interpolated includes:
setting an xy plane interpolation range and an interpolation interval in the z direction based on the ground temperature data;
setting an interpolation distance based on the interpolation interval in the z direction;
according to the interpolation distance, equally dividing the interpolation interval to obtain heights corresponding to one or more planes needing interpolation;
respectively setting the density of the plane interpolation grid in the x and y directions based on the xy plane interpolation range;
and obtaining grid information to be interpolated based on the effective height and the density of the plane interpolation grid in the x and y directions.
Further, the interpolation section in the z direction is set by:
the number of effective drilling points is determined by the following method:
N(z)>N eff ＝N×a
wherein N is eff A threshold value for effective drill points;
n is the total number of drilling points;
a is a condition percentage parameter;
and determining an interpolation interval in the z direction based on the number of the effective drilling points and the maximum value and the minimum value of the heights in each drilling.
Further, performing interpolation processing on the horizontal grid corresponding to the grid information to obtain frozen earth point cloud data includes:
based on the grid information to be interpolated, obtaining temperature values and position information of effective data points of corresponding planes of different heights;
sequentially acquiring data scatter diagrams of known point distances and half variances of the temperature values of the known point distances in different planes based on the temperature values and the position information;
selecting a corresponding variation function based on the data scatter diagram; the variation function comprises a spherical model, a Gaussian model and an exponential model;
interpolating a horizontal grid in a current plane based on the variation function to obtain temperature fitting values of all points in the horizontal grid;
obtaining frozen soil point cloud data based on the temperature fitting values of all the points;
and determining the parameters of the variation function by adopting a least square method.
Further, performing tetrahedral subdivision on the frozen earth point cloud data; according to the subdivision result of the frozen soil point cloud data, the establishing of the three-dimensional frozen soil model comprises the following steps:
defining a tetrahedron of which the outer sphere contains all frozen earth point cloud data;
adding frozen soil point cloud data into the tetrahedrons one by one, finding the tetrahedrons where new adding points are located, and deleting the tetrahedrons where the new adding points are located;
calculating the externally connected balls of all tetrahedrons adjacent to the deleted tetrahedrons, and if the externally connected balls comprise newly added points, connecting the newly added points with the points on the externally connected balls to obtain three new tetrahedrons;
repeating the steps until all the frozen soil point cloud data are added into tetrahedral subdivision to obtain a frozen soil grid set;
with tetrahedron as screening condition, calculating the side length set { L } of any tetrahedron i |l i1 ,l i2 ,l i3 ,l i4 ,l i5 ,l i6 }：
Wherein j and k are sides l respectively im Point numbers at both ends;
m=1, 2,3 … … 6; screening all tetrahedrons based on the side length set to obtain tetrahedrons conforming to grid topology conditions;
and establishing a three-dimensional frozen soil model based on the tetrahedron conforming to the grid topological condition.
Further, the screening all tetrahedrons based on the side length set to obtain tetrahedrons meeting the grid topology condition includes:
deleting tetrahedrons meeting the following conditions to obtain tetrahedrons meeting the grid topology conditions:
wherein a, b and c are interpolation densities in three directions of x, y and z.
In a second aspect of the application, a construction device for a frozen soil model is provided. The device comprises:
the acquisition module is used for acquiring the ground temperature data;
the setting module is used for setting space three-dimensional interpolation grid parameters based on the ground temperature data to obtain grid information to be interpolated;
the interpolation module is used for carrying out interpolation processing on the horizontal grids corresponding to the grid information to obtain frozen soil point cloud data;
the construction module is used for carrying out tetrahedral subdivision on the frozen earth point cloud data; and establishing a three-dimensional frozen soil model according to the subdivision result of the frozen soil point cloud data.
In a third aspect of the application, an electronic device is provided. The electronic device includes: a memory and a processor, the memory having stored thereon a computer program, the processor implementing the method as described above when executing the program.
In a fourth aspect of the application, there is provided a computer readable storage medium having stored thereon a computer program which when executed by a processor implements a method as according to the first aspect of the application.
According to the construction method of the frozen soil model, provided by the embodiment of the application, the ground temperature data are obtained; setting space three-dimensional interpolation grid parameters based on the ground temperature data to obtain grid information to be interpolated; performing interpolation processing on the horizontal grids corresponding to the grid information to obtain frozen earth point cloud data; performing tetrahedral subdivision on the frozen earth point cloud data; according to the subdivision result of the frozen soil point cloud data, a three-dimensional frozen soil model is established, frozen soil distribution in a reservoir area can be represented from a three-dimensional rather than planar angle under the condition of fully utilizing drilling data, the processing flow is simplified, and meanwhile, the accuracy of data quantization is greatly improved.
It should be understood that the description in this summary is not intended to limit the critical or essential features of the embodiments of the application, nor is it intended to limit the scope of the application. Other features of the present application will become apparent from the description that follows.
Drawings
The above and other features, advantages and aspects of embodiments of the present application will become more apparent by reference to the following detailed description when taken in conjunction with the accompanying drawings. In the drawings, wherein like or similar reference numerals denote like or similar elements, in which:
FIG. 1 is a flow chart of a method of constructing a frozen soil model according to an embodiment of the application;
FIG. 2 is a flowchart of interpolation interval calculation according to an embodiment of the present application;
FIG. 3 is a schematic diagram of borehole data points in accordance with an embodiment of the present application;
FIG. 4 is a schematic diagram of an active data point set according to an embodiment of the present application;
FIG. 5 is a schematic illustration of an interpolation grid according to an embodiment of the present application;
FIG. 6 is a grid contrast schematic in accordance with an embodiment of the present application;
FIG. 7 is a schematic diagram of a three-dimensional ice layer model according to an embodiment of the application;
FIG. 8 is a schematic view of a construction apparatus of a frozen soil model according to an embodiment of the application;
fig. 9 is a schematic diagram of a structure of a terminal device or server suitable for implementing an embodiment of the present application.
Detailed Description
For the purposes of making the objects, technical solutions and advantages of the embodiments of the present disclosure more apparent, the technical solutions of the embodiments of the present disclosure will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present disclosure, and it is apparent that the described embodiments are some embodiments of the present disclosure, but not all embodiments. All other embodiments, which can be made by one of ordinary skill in the art based on the embodiments in this disclosure without inventive faculty, are intended to be within the scope of this disclosure.
In addition, the term "and/or" herein is merely an association relationship describing an association object, and means that three relationships may exist, for example, a and/or B may mean: a exists alone, A and B exist together, and B exists alone. In addition, the character "/" herein generally indicates that the front and rear associated objects are an "or" relationship.
Fig. 1 shows a flowchart of a method of constructing a frozen soil model according to an embodiment of the disclosure. The method comprises the following steps:
s110, acquiring ground temperature data.
In this embodiment, the execution body of the construction method for the frozen soil model may acquire the ground temperature data in a wired manner or a wireless connection manner.
Further, the execution body may acquire the ground temperature data transmitted from the electronic device connected to the execution body in communication, or may store the ground temperature data locally in advance.
The ground temperature data comprise horizontal position, depth, temperature and other information of the drilling hole.
In some embodiments, to increase the subsequent calculation speed, the directly acquired ground temperature data may be cleaned to remove redundant interference data. For example, the geothermal data is extracted to an excel file, and the data is processed by pandas to delete invalid borehole data.
In some embodiments, a drilling information database can be constructed based on the obtained ground temperature data, and the database can be directly called for data processing in the subsequent steps, so that the data processing speed is improved, and the operations of adding, deleting, checking, changing and the like can be conveniently performed on the data by a manager.
Specifically, the borehole information database may be represented as:
{K|k 1 ,k 2 ,k 3 ,…,k i ,…,k N }
wherein N is the total number of drilling points;
k i a data set for any one of the drill points;
for each drilling point data set k i Consisting of a series of points containing depth information and temperature information, can be expressed as:
{k i |p i1 ,p i2 ,p i3 ,…,p ij ,…,p in }
wherein p is ij For the drilling point (k i ) Information of the jth data point;
further, define X coordinate information X ij ；
Y-coordinate information Y ij ；
Height coordinate information z ij ；
Temperature information t ij ；
That is, { p ij |x ij ,y ij ,z ij ,t ij }}；
For each borehole data set k i According to where z ij And t ij In the interval { z } min―i ,z max―i Linear interpolation is carried out in the interval to obtain the functional relation between the height of each drilling point in the interval and the temperature:
t＝f i (z)
wherein t is a temperature value corresponding to z height;
that is, a borehole interpolation function set is obtained:
{F|f 1 (z),f 2 (z),…,f i (z),…,f N (z)}
further, the method comprises the steps of,
and S120, setting a spatial three-dimensional interpolation grid parameter based on the ground temperature data to obtain grid information to be interpolated.
In some embodiments, parameters of the spatial three-dimensional interpolation grid are set based on the borehole information database constructed in step S110. That is, an interpolation section in the vertical z direction and an interpolation region in the xy plane are defined, and on the basis of the interpolation section, interpolation densities a, b, and c of the three-dimensional interpolation grid in the x, y, and z directions are set, respectively, thereby obtaining grid information to be interpolated.
Specifically, setting an interpolation range O of an xy plane according to an actual application scene; in order to ensure the reliability of the interpolation result in each plane, the number N (z) of effective drilling points in the plane needs to be set:
N(Z)>N eff ＝N×a
wherein N is eff A threshold value for effective drill points;
n is the total number of drilling points;
a is a condition percentage parameter, which can be set according to the actual application scene, and is preferably 0.8;
as shown in fig. 2, fig. 2 shows a specific calculation manner of Domian, and the set vertical interpolation interval may be expressed as:
Domian＝{Z min ,Z max }
as shown in fig. 3, fig. 3 is a schematic diagram of borehole data points. For any height z ε Domian, N (z) needs to be satisfied>N eff For all drilling points k i E, K, obtaining the minimum value and the maximum value of the height in each drilling hole:
minimum value:
{z min |z min―1 ,z min―2 ,…,z min―i ,…,z min―N }
maximum value:
{z max |z max―1 ,z max―2 ,…,z max―i ,…,z max―N }
further, after determining an interpolation interval Domian, setting an interpolation distance d (grid vertical density), and equally dividing the interpolation interval according to the interpolation distance d to obtain heights corresponding to m planes needing interpolation:
{H|h 0 ,h 1 ,…,h k ,…,h m }
further, referring to fig. 4, fig. 4 shows a height h k Is provided for the flat active data point set. Based on each interpolation height, obtaining effective data points of the interpolation plane:
{P eff (h k )|P eff―1 (h k ),P eff―2 (h k ),…,P eff―i (h k ),…,P eff―N (h k )}
wherein:
P eff―i (h k )＝f i (h k )
as shown in fig. 5, densities a and b of the planar interpolation grid in the x and y directions are set in the interpolation range O, and the interpolation grid is obtained:
{G|G 0 ,G 1 ,…,G i ,…,G m }
wherein G is i For the horizontal grid point set corresponding to the interpolation plane i:
wherein,the total number of interpolation points in the plane;
p ij is the interpolation point with the sequence number j in the plane.
Further, referring to fig. 6, fig. 6 is a grid contrast schematic diagram showing a grid constructed by the method of the present disclosure and a grid constructed by the prior art, respectively.
And S130, carrying out interpolation processing on the horizontal grids corresponding to the grid information to obtain frozen soil point cloud data.
In some embodiments, temperature values and location information for valid data points for different heights corresponding to planes may be obtained from a borehole information database.
Sequentially obtaining data scatter diagrams of known point distances and half variances of the temperature values of the known point distances in different planes based on the temperature values and the position information, and simultaneously selecting a proper variation function; wherein, the least square method can be adopted to determine the parameters of the variation function;
further, interpolation is carried out on the horizontal grids in the plane through the variation function, temperature field point cloud data are obtained, and points meeting the conditions are screened out according to the temperature field point cloud data to serve as frozen soil data points.
Specifically, for any height h k Interpolation planes corresponding to the interpolation planes, the interpolation includes:
a，drilling data P based on current altitude eff (h k ) And calculating distance and half variance for N data in the data in pairs to obtain a distance matrix and half variance:
wherein, the distance matrix is:
{d|[d 11 ,d 12 ,…,d 1N ],[d 21 ,d 22 ,…,d 2N ],…,[d N1 ,d N2 ,…,d NN ]}
the half variance is:
{r|[r 11 ,r 12 ,…,r 1N ],[r 21 ,r 22 ,…,r 2N ],…,[r N1 ,r N2 ,…,r NN ]}
b, drawing r and d scatter diagrams according to the values and the corresponding relations in the d and r matrixes;
c, selecting an optimal variation function type according to the scatter diagram, and calculating variation function parameters through a least square method to obtain a function r=gamma (d);
the variation function may be selected as follows:
determining the relative position relation among the data points according to the scatter diagram;
and selecting an optimal variation function based on the position relation, wherein the variation function comprises a spherical model, a Gaussian model, an exponential model and the like. For example, when the current position relationship satisfies a preset condition, selecting a spherical model, otherwise, selecting an index model and the like, which are not particularly limited herein; the preset conditions can be set according to actual application scenes;
d, for at height h k To be interpolated G in i Unknown point p of (2) ij Calculate it to all known points P eff (h k ) Distance { L|L } 1 ,L 2 ,…,L N -a }; at the same time, the half variance estimation value { R|R } of the point is calculated by the following equation set 1 ,R 2 ,…,R N }：
Further, solving the equation set to obtain an optimal coefficient lambda i ；
e, based on the optimal coefficient lambda i And carrying out weighted summation on the attribute values of the known points to obtain an estimated value of the temperature of the unknown points:
f, repeating the steps a-e, and traversing all points to be interpolated of the current height. That is, walk { H|h 0 ,h 1 ,…,h k ,h m And (3) carrying out interpolation processing on all the heights in the grid G to obtain temperature fitting values of all the points in the interpolation grid G, and taking the points with the temperature fitting values smaller than 1 as east earth stores to obtain a frozen earth point set (frozen earth point cloud data).
S140, performing tetrahedral subdivision on the frozen earth point cloud data; and establishing a three-dimensional frozen soil model according to the subdivision result of the frozen soil point cloud data.
a, defining a tetrahedron of which the outer sphere contains all frozen soil point cloud data;
b, adding frozen earth point cloud data into the tetrahedrons one by one, finding the tetrahedrons where the new adding points are located, and deleting the tetrahedrons where the new adding points are located;
c, calculating the outer balls of all tetrahedrons adjacent to the deleted tetrahedrons, and if the outer balls comprise newly added points, connecting the newly added points with the points on the outer balls to obtain three new tetrahedrons; if the outer ball does not include the newly added point, no treatment is carried out;
d, repeating the steps b and c until all the frozen soil point cloud data are added into the tetrahedral subdivision to obtain a frozen soil grid set:
{M|m 1 ,m 2 ,…,m i ,…,m u }
e, taking tetrahedron as screening condition, calculating the side length set { L } of any tetrahedron i |l i1 ,l i2 ,l i3 ,l i4 ,l i5 ,l i6 }：
Wherein j and k are sides l respectively im Point numbers at both ends;
m＝1，2，3……6；
further, based on the side length set, deleting the tetrahedrons meeting the following conditions to obtain the tetrahedrons meeting the grid topology conditions:
wherein a, b and c are interpolation densities in the directions of x, y and z;
f, based on the tetrahedron conforming to the grid topology condition, a three-dimensional frozen soil model can be established by means of modeling software and the like, and referring to fig. 7, fig. 7 is a schematic diagram of the three-dimensional frozen soil model according to the embodiment of the application.
Different from the existing tetrahedron segmentation technology, the tetrahedron segmentation method solves the problem of common boundary information processing in tetrahedron segmentation, and greatly improves the construction precision of the model.
According to the embodiment of the disclosure, the following technical effects are achieved:
the three-dimensional frozen soil model is automatically created, and the whole working flow can be completed only by inputting parameters, so that the generation steps of the model are greatly simplified;
according to the frozen soil model obtained by the method, the frozen soil distribution condition of any horizontal plane or vertical section can be obtained without depending on the distribution condition of drilling data, and the problem of insufficient plane observation data is solved.
It should be noted that, for simplicity of description, the foregoing method embodiments are all described as a series of acts, but it should be understood by those skilled in the art that the present application is not limited by the order of acts described, as some steps may be performed in other orders or concurrently in accordance with the present application. Further, those skilled in the art will also appreciate that the embodiments described in the specification are alternative embodiments, and that the acts and modules referred to are not necessarily required for the present application.
The above description of the method embodiments further describes the solution of the present application by means of device embodiments.
Fig. 8 shows a construction apparatus 800 of a frozen soil model according to an embodiment of the application, including, as shown in fig. 8:
an acquisition module 810 for acquiring ground temperature data;
a setting module 820, configured to set a spatial three-dimensional interpolation grid parameter based on the ground temperature data, to obtain grid information to be interpolated;
the interpolation module 830 is configured to perform interpolation processing on the horizontal grid corresponding to the grid information to obtain frozen earth point cloud data;
a construction module 840, configured to perform tetrahedral subdivision on the frozen earth point cloud data; according to the subdivision result of the frozen soil point cloud data, a three-dimensional frozen soil model is established
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the described modules may refer to corresponding procedures in the foregoing method embodiments, which are not described herein again.
Fig. 9 shows a schematic diagram of a structure of a terminal device or server suitable for implementing an embodiment of the application.
As shown in fig. 9, the terminal device or the server includes a Central Processing Unit (CPU) 901, which can execute various appropriate actions and processes in accordance with a program stored in a Read Only Memory (ROM) 902 or a program loaded from a storage section 908 into a Random Access Memory (RAM) 903. In the RAM903, various programs and data required for the operation of the terminal device or the server are also stored. The CPU901, ROM 902, and RAM903 are connected to each other through a bus 904. An input/output (I/O) interface 905 is also connected to the bus 904.
The following components are connected to the I/O interface 905: an input portion 806 including a keyboard, mouse, etc.; an output portion 807 including a display such as a Cathode Ray Tube (CRT), a Liquid Crystal Display (LCD), and a speaker; a storage portion 908 including a hard disk or the like; and a communication section 909 including a network interface card such as a LAN card, a modem, or the like. The communication section 909 performs communication processing via a network such as the internet. The drive 910 is also connected to the I/O interface 905 as needed. A removable medium 911 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is installed as needed on the drive 910 so that a computer program read out therefrom is installed into the storage section 908 as needed.
In particular, the above method flow steps may be implemented as a computer software program according to an embodiment of the application. For example, embodiments of the application include a computer program product comprising a computer program embodied on a machine-readable medium, the computer program comprising program code for performing the method shown in the flowcharts. In such an embodiment, the computer program may be downloaded and installed from the network via the communication portion 909 and/or installed from the removable medium 911. The above-described functions defined in the system of the present application are performed when the computer program is executed by a Central Processing Unit (CPU) 901.
The computer readable medium shown in the present application may be a computer readable signal medium or a computer readable storage medium, or any combination of the two. The computer readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any of the foregoing. More specific examples of the computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. In the present application, however, the computer-readable signal medium may include a data signal propagated in baseband or as part of a carrier wave, with the computer-readable program code embodied therein. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: wireless, wire, fiber optic cable, RF, etc., or any suitable combination of the foregoing.
The flowcharts and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The units or modules involved in the embodiments of the present application may be implemented in software or in hardware. The described units or modules may also be provided in a processor. Wherein the names of the units or modules do not in some cases constitute a limitation of the units or modules themselves.
As another aspect, the present application also provides a computer-readable storage medium that may be contained in the electronic device described in the above embodiment; or may be present alone without being incorporated into the electronic device. The computer-readable storage medium stores one or more programs that when executed by one or more processors perform the methods described herein.
The above description is only illustrative of the preferred embodiments of the present application and of the principles of the technology employed. It will be appreciated by persons skilled in the art that the scope of the application is not limited to the specific combinations of the features described above, but also covers other embodiments which may be formed by any combination of the features described above or their equivalents without departing from the spirit of the application. Such as the above-mentioned features and the technical features having similar functions (but not limited to) applied for in the present application are replaced with each other.
Claims (10)
1. The construction method of the frozen soil model is characterized by comprising the following steps of:
acquiring ground temperature data;
setting space three-dimensional interpolation grid parameters based on the ground temperature data to obtain grid information to be interpolated;
performing interpolation processing on the horizontal grids corresponding to the grid information to obtain frozen earth point cloud data;
performing tetrahedral subdivision on the frozen earth point cloud data; and establishing a three-dimensional frozen soil model according to the subdivision result of the frozen soil point cloud data.
2. The method of claim 1, wherein the geothermal data comprises borehole level, borehole depth, and borehole temperature information.
3. The method of claim 2, wherein setting spatial three-dimensional interpolation grid parameters based on the ground temperature data, and obtaining grid information to be interpolated comprises:
setting an xy plane interpolation range and an interpolation interval in the z direction based on the ground temperature data;
setting an interpolation distance based on the interpolation interval in the z direction;
according to the interpolation distance, equally dividing the interpolation interval to obtain heights corresponding to one or more planes needing interpolation;
respectively setting the density of the plane interpolation grid in the x and y directions based on the xy plane interpolation range;
and obtaining grid information to be interpolated based on the height and the density of the plane interpolation grid in the x and y directions.
4. A method according to claim 3, characterized in that the interpolation interval in the z-direction is set by:
the number of effective drilling points is determined by the following method:
N(z)>N eff ＝N×a
wherein N is eff A threshold value for effective drill points;
n is the total number of drilling points;
a is a condition percentage parameter;
and determining an interpolation interval in the z direction based on the number of the effective drilling points and the maximum value and the minimum value of the heights in each drilling.
5. The method of claim 4, wherein the interpolating the horizontal grid corresponding to the grid information to obtain the frozen earth point cloud data comprises:
based on the grid information to be interpolated, obtaining temperature values and position information of effective data points of corresponding planes of different heights;
sequentially acquiring data scatter diagrams of known point distances and half variances of the temperature values of the known point distances in different planes based on the temperature values and the position information;
selecting a corresponding variation function based on the data scatter diagram; the variation function comprises a spherical model, a Gaussian model and an exponential model;
interpolating a horizontal grid in a current plane based on the variation function to obtain temperature fitting values of all points in the horizontal grid;
obtaining frozen soil point cloud data based on the temperature fitting values of all the points;
and determining the parameters of the variation function by adopting a least square method.
6. The method of claim 5, wherein the tetrahedrally dissecting the frozen earth point cloud data; according to the subdivision result of the frozen soil point cloud data, the establishing of the three-dimensional frozen soil model comprises the following steps:
defining a tetrahedron of which the outer sphere contains all frozen earth point cloud data;
adding frozen soil point cloud data into the tetrahedrons one by one, finding the tetrahedrons where new adding points are located, and deleting the tetrahedrons where the new adding points are located;
calculating the externally connected balls of all tetrahedrons adjacent to the deleted tetrahedrons, and if the externally connected balls comprise newly added points, connecting the newly added points with the points on the externally connected balls to obtain three new tetrahedrons;
repeating the steps until all the frozen soil point cloud data are added into tetrahedral subdivision to obtain a frozen soil grid set;
with tetrahedron as screening condition, calculating the side length set { L } of any tetrahedron i |l i1 ,l i2 ,l i3 ,l i4 ,l i5 ,l i6 }：
Wherein j and k are sides l respectively im Point numbers at both ends;
m＝1，2，3……6；
screening all tetrahedrons based on the side length set to obtain tetrahedrons conforming to grid topology conditions;
and establishing a three-dimensional frozen soil model based on the tetrahedron conforming to the grid topological condition.
7. The method of claim 6, wherein screening all tetrahedrons based on the set of side lengths to obtain tetrahedrons that meet a mesh topology comprises:
deleting tetrahedrons meeting the following conditions to obtain tetrahedrons meeting the grid topology conditions:
wherein a, b and c are interpolation densities in three directions of x, y and z.
8. The construction device of the frozen soil model is characterized by comprising the following components:
the acquisition module is used for acquiring the ground temperature data;
the setting module is used for setting space three-dimensional interpolation grid parameters based on the ground temperature data to obtain grid information to be interpolated;
the interpolation module is used for carrying out interpolation processing on the horizontal grids corresponding to the grid information to obtain frozen soil point cloud data;
the construction module is used for carrying out tetrahedral subdivision on the frozen earth point cloud data; and establishing a three-dimensional frozen soil model according to the subdivision result of the frozen soil point cloud data.
9. An electronic device comprising a memory and a processor, the memory having stored thereon a computer program, characterized in that the processor, when executing the computer program, implements the method according to any of claims 1-7.
10. A computer readable storage medium, on which a computer program is stored, which computer program, when being executed by a processor, implements the method according to any of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311156890.0A CN117152369A (en) | 2023-09-08 | 2023-09-08 | Construction method, device, equipment and computer readable storage medium of frozen soil model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311156890.0A CN117152369A (en) | 2023-09-08 | 2023-09-08 | Construction method, device, equipment and computer readable storage medium of frozen soil model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117152369A true CN117152369A (en) | 2023-12-01 |
Family
ID=88909778
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311156890.0A Pending CN117152369A (en) | 2023-09-08 | 2023-09-08 | Construction method, device, equipment and computer readable storage medium of frozen soil model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117152369A (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110097637A (en) * | 2019-05-09 | 2019-08-06 | 正元地理信息集团股份有限公司 | A kind of three-dimensional geological attribute model temporal-spatial interpolating method and system |
| CN110244021A (en) * | 2019-06-14 | 2019-09-17 | 深圳市秉睦科技有限公司 | A kind of strata division method based on anisotropy interpolation |
| CN111127630A (en) * | 2019-12-17 | 2020-05-08 | 核工业北京地质研究院 | Random modeling method for hyperspectral alteration information of drill core in three-dimensional space |
| CN112669463A (en) * | 2020-12-25 | 2021-04-16 | 河南信大融通信息科技有限公司 | Method for reconstructing curved surface of three-dimensional point cloud, computer device and computer-readable storage medium |
| CN113487733A (en) * | 2021-06-02 | 2021-10-08 | 长江勘测规划设计研究有限责任公司 | Method for quickly establishing three-dimensional model of bottom surface of covering layer |
| CN115830252A (en) * | 2022-11-09 | 2023-03-21 | 河海大学 | B/S framework site underground water organic pollution simulation system and online simulation method |
| CN116645484A (en) * | 2023-07-26 | 2023-08-25 | 航天宏图信息技术股份有限公司 | Geological curved surface model construction method and device, electronic equipment and storage medium |
-
2023
- 2023-09-08 CN CN202311156890.0A patent/CN117152369A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110097637A (en) * | 2019-05-09 | 2019-08-06 | 正元地理信息集团股份有限公司 | A kind of three-dimensional geological attribute model temporal-spatial interpolating method and system |
| CN110244021A (en) * | 2019-06-14 | 2019-09-17 | 深圳市秉睦科技有限公司 | A kind of strata division method based on anisotropy interpolation |
| CN111127630A (en) * | 2019-12-17 | 2020-05-08 | 核工业北京地质研究院 | Random modeling method for hyperspectral alteration information of drill core in three-dimensional space |
| CN112669463A (en) * | 2020-12-25 | 2021-04-16 | 河南信大融通信息科技有限公司 | Method for reconstructing curved surface of three-dimensional point cloud, computer device and computer-readable storage medium |
| CN113487733A (en) * | 2021-06-02 | 2021-10-08 | 长江勘测规划设计研究有限责任公司 | Method for quickly establishing three-dimensional model of bottom surface of covering layer |
| CN115830252A (en) * | 2022-11-09 | 2023-03-21 | 河海大学 | B/S framework site underground water organic pollution simulation system and online simulation method |
| CN116645484A (en) * | 2023-07-26 | 2023-08-25 | 航天宏图信息技术股份有限公司 | Geological curved surface model construction method and device, electronic equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9959670B2 (en) | Method for rendering terrain | |
| Souza et al. | Occupancy-elevation grid: an alternative approach for robotic mapping and navigation | |
| CN112652036A (en) | Road data processing method, device, equipment and storage medium | |
| CN110827405A (en) | Digital remote sensing geological mapping method and system | |
| CN114332291A (en) | Oblique photography model building outer contour rule extraction method | |
| CN114693855B (en) | Point cloud data processing method and device | |
| CN116051635A (en) | Pose determining method and device | |
| CN117152369A (en) | Construction method, device, equipment and computer readable storage medium of frozen soil model | |
| CN111080080B (en) | Village geological disaster risk prediction method and system | |
| CN111553986B (en) | Construction method and construction device of triangulation network and generation method of digital surface model | |
| CN116246069B (en) | Method and device for self-adaptive terrain point cloud filtering, intelligent terminal and storage medium | |
| CN110363863B (en) | Input data generation method and system of neural network | |
| CN111598941A (en) | Method, device and equipment for measuring gradient of tower and storage medium | |
| CN116307325A (en) | Line planning method and device for power distribution network, electronic equipment and storage medium | |
| US10101497B2 (en) | Method for determining a karstic region | |
| Gökgöz et al. | Comparison of two methods for multiresolution terrain modelling in GIS | |
| Ridzuan et al. | Voxelization techniques: data segmentation and data modelling for 3d building models | |
| CN114998547B (en) | Grid random sampling method, system and storage medium for power line and pole tower points | |
| CN116957414B (en) | Village planning analysis method and device based on artificial intelligence | |
| CN117216936A (en) | Land and stone visualization deployment method based on unmanned aerial vehicle multi-source data twinning | |
| CN117075196A (en) | Crack morphology analysis method, device, equipment and medium | |
| CN116030131A (en) | Pose determining method and device | |
| CN117788741A (en) | Building model reconstruction method, device, equipment and storage medium based on point cloud | |
| CN118115670A (en) | Unmanned aerial vehicle viewpoint generation method, device and equipment based on incremental point cloud data | |
| KR101381811B1 (en) | Method and apparatus for simplifying a height map and storage medium storing the same |
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 |