Method and device for obtaining frost heaving amount in permafrost region, electronic equipment and storage medium

Abstract

This invention relates to the field of natural gas pipeline safety treatment technology, and provides a method, apparatus, electronic device, and storage medium for obtaining frost heave in permafrost areas. The method includes: determining a surface temperature function; determining a constant temperature at a first depth below the surface; constructing a finite element model of the permafrost region calculation area based on the surface temperature function, the first depth, and the constant temperature; applying boundary conditions to the finite element model to obtain the time it takes for a freezing front to reach the lower edge of each grid layer in an M-layer grid from the surface down within the time range of the surface temperature function, as well as the temperature gradient at the location of the freezing front; and determining the frost heave based on the time and the temperature gradient. This method enables rapid and accurate determination of frost heave, providing fundamental data for the safety assessment of buried pipelines.

Description

Technical Field

[0001] This invention relates to the field of natural gas transportation pipeline safety treatment technology, and in particular to a method, apparatus, electronic device and storage medium for obtaining frost heave in permafrost areas. Background Technology

[0002] Permafrost refers to various rocks and soils containing ice and located below zero degrees Celsius. my country is the third largest country in the world with permafrost, covering approximately 21.5% of its land area. Permafrost is mainly distributed in the mid- and low-latitude Qinghai-Tibet Plateau, with other distributions in the Pamir Mountains, western high mountains, the Greater and Lesser Khingan Mountains in the northeast, and the east. Oil and gas pipelines inevitably traverse permafrost zones. For example, the Mohe-Daqing section of the China-Russia crude oil pipeline, one of my country's four major energy strategic channels, is 953 kilometers long, with its northern 512 kilometers traversing permafrost regions. Pipelines are subjected to the pressure of the medium inside the pipe, the weight of the medium, soil pressure, and stress generated by ground subsidence and deformation due to temperature changes, placing them in a highly complex stress state. Differential frost heave is the main cause of pipeline failure in permafrost regions. When moisture in permafrost freezes in winter, the accumulation of ice lenses causes the ground surface to bulge; this phenomenon is called frost heave. If the frost heave is too large, it can cause excessive deformation of the pipeline, leading to failure.

[0003] When designing or conducting safety assessments of pipelines in permafrost regions, it is necessary to obtain key parameters and characteristics of the permafrost area, such as the changes in frost heave depth and frost heave over time. Existing methods struggle to obtain these parameters in permafrost regions. Summary of the Invention

[0004] To address the problems existing in the prior art, the present invention provides a method, apparatus, electronic device, and storage medium for obtaining frost heave in permafrost regions.

[0005] In a first aspect, the present invention provides a method for obtaining frost heave in permafrost regions, comprising:

[0006] Determine the surface temperature function;

[0007] Determine the constant temperature at the first depth below the Earth's surface;

[0008] A finite element model of the calculation area of ​​the permafrost region is constructed based on the surface temperature function, the first depth, and the constant temperature; wherein, the upper boundary of the finite element model is the surface, and the lower boundary is the first depth; the principle for determining the left and right boundaries is that the temperature gradient level change at the left and right boundaries is less than K℃ / m, where K is a preset value; within the calculation area of ​​the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary;

[0009] Boundary conditions are applied to the finite element model to obtain the time it takes for the freezing front to reach the bottom of each grid layer in the M-layer grid from the surface to the bottom within the time range of the surface temperature function, as well as the temperature gradient at the location of the freezing front.

[0010] The amount of frost heave is determined based on the time and the temperature gradient.

[0011] In one embodiment, the surface temperature function includes:

[0012]

[0013] Where T(x,0,t) is the ambient temperature, T m Let A be the annual average atmospheric temperature, A be the amplitude of atmospheric temperature change, and t be time.

[0014] In one embodiment, the boundary conditions include: adiabatic left and right boundaries, a constant temperature lower boundary, and an upper boundary applied based on the surface temperature function.

[0015] In one embodiment, determining the frost heave amount based on the time and the temperature gradient includes:

[0016] The amount of frost heave is determined using the following formula based on the time and the temperature gradient;

[0017]

[0018] Among them, h f P represents the frost heave. e The external load is represented by a and b, which are soil-related parameters, gradT(t) is the temperature gradient, and V is the temperature gradient. t This represents the change over time.

[0019] In a second aspect, the present invention provides a device for obtaining frost heave in permafrost areas, comprising:

[0020] The fitting module is used to determine the surface temperature function;

[0021] The acquisition module is used to determine the constant temperature at the first depth below the Earth's surface.

[0022] The construction module is used to construct a finite element model of the permafrost region calculation area based on the surface temperature function, the first depth, and the constant temperature; wherein, the upper boundary of the finite element model is the surface, and the lower boundary is the first depth; the principle for determining the left and right boundaries is that the temperature gradient level change at the left and right boundaries is less than K℃ / m, where K is a preset value; within the calculation area of ​​the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary;

[0023] The calculation module is used to apply boundary conditions to the finite element model and obtain the time from the surface temperature function to the bottom of each grid in the M-layer grid from the surface to the bottom, as well as the temperature gradient at the location of the freezing front.

[0024] A determination module is used to determine the amount of frost heave based on the time and the temperature gradient.

[0025] Thirdly, the present invention provides an electronic device, including a memory and a memory storing a computer program, wherein the processor executes the program to implement the steps of the method for obtaining frost heave in permafrost regions as described in the first aspect.

[0026] Fourthly, the present invention provides a processor-readable storage medium storing a computer program for causing the processor to perform the steps of the method for obtaining frost heave in permafrost regions as described in the first aspect.

[0027] The present invention provides a method, apparatus, electronic device, and storage medium for obtaining frost heave in permafrost regions. By constructing a finite element model of the calculation area of ​​the permafrost region based on the surface temperature function, a first depth, and a constant temperature, and applying boundary conditions to the finite element model, the temperature field changes of the calculation area within the calculation time range of the surface temperature function are simulated, and the changes of frost heave depth and frost heave amount over time are obtained. Based on these changes, the maximum frost heave amount can be quickly and accurately determined, providing basic data for the safety assessment and design of buried pipelines. The method is low-cost, easy to operate, quick to calculate, and accurate in providing results. Attached Figure Description

[0028] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0029] Figure 1 This is a flowchart illustrating the method for obtaining frost heave in permafrost regions provided by the present invention.

[0030] Figure 2 This is a schematic diagram of the temperature field changes on the 12th day provided by the present invention;

[0031] Figure 3 This is a schematic diagram of the temperature field changes on day 152 provided by the present invention;

[0032] Figure 4 This is a diagram showing the location of the freezing front and the temperature gradient at the freezing front in the temperature field provided by the present invention.

[0033] Figure 5 This is a schematic diagram of the device for obtaining frost heave in permafrost regions provided by the present invention;

[0034] Figure 6 This is a schematic diagram of the structure of the electronic device provided by the present invention. Detailed Implementation

[0035] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0036] The following is combined Figures 1-6 The present invention describes a method, apparatus, electronic device, and storage medium for obtaining frost heave in permafrost regions.

[0037] Figure 1 A flowchart illustrating a method for obtaining frost heave in permafrost regions according to the present invention is shown below. Figure 1 The method includes:

[0038] 11. Determine the surface temperature function;

[0039] 12. Determine the constant temperature at the first depth below the Earth's surface;

[0040] 13. Construct a finite element model of the calculation area of ​​the permafrost region based on the surface temperature function, the first depth, and the constant temperature; wherein, the upper boundary of the finite element model is the surface, and the lower boundary is the first depth; the principle for determining the left and right boundaries is that the temperature gradient at the left and right boundaries is less than K℃ / m, where K is a preset value; within the calculation area of ​​the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary;

[0041] 14. Apply boundary conditions to the finite element model to obtain the time it takes for the freezing front to reach the bottom of each grid layer in the M-layer grid from the surface to the bottom within the time range of the surface temperature function, as well as the temperature gradient at the location of the freezing front.

[0042] 15. Determine the amount of frost heave based on the time and temperature gradient.

[0043] Regarding steps 11 to 14, it should be noted that in this invention, temperature monitoring is used to record the changes in surface temperature over time in different regions (including permafrost and non-permafrost regions). The changes in surface temperature are then fitted to obtain a fitting function, which serves as the surface temperature function.

[0044] A further method of the present invention, wherein the surface temperature function includes:

[0045]

[0046] Where T(x,0,t) is the ambient temperature, T m Let A be the annual average atmospheric temperature, A be the amplitude of atmospheric temperature change, and t be time.

[0047] In this invention, it should be noted that by ignoring horizontal temperature variations at the same depth, the Earth can be approximated as a semi-infinite object, with temperature changes occurring only in the depth direction. Under the influence of the aforementioned boundary condition of periodic atmospheric temperature variations, the Earth's temperature field also exhibits the characteristic of temperature wave attenuation. That is, the temperature variation at any depth within the Earth over time follows a sinusoidal function with the same period as the temperature variation at the surface, but the amplitude decreases with increasing depth. The deeper the depth, the greater the amplitude attenuation; therefore, when the depth is sufficiently deep, the amplitude of temperature fluctuations attenuates to a negligible level. When the temperature value at a certain location below a certain distance from the pipeline's burial depth differs from the natural Earth temperature field value at the same depth by less than a (e.g., a = 0.5℃), the Earth's temperature field at this depth is no longer affected by the pipeline temperature. Let this depth be H, and the constant temperature be T0.

[0048] For example, tests have shown that a certain location has a constant temperature of 5°C at a depth of 10m below the surface.

[0049] In this invention, a finite element model of the calculation area of ​​the permafrost region is constructed based on the surface temperature function, the first depth, and the constant temperature. The finite element model is a model established using the finite element analysis method, which is a combination of elements connected only at the nodes, transmitting forces only through the nodes, and being constrained only at the nodes.

[0050] The upper boundary of the finite element model is the Earth's surface, and the lower boundary is the depth H of the determined isothermal point below the Earth's surface. The principles for determining the left and right boundaries are as follows: the horizontal change in temperature gradient at the left and right boundaries is less than K℃ / m, where K is a preset value, for example, K is 0.5℃ / m. The natural temperature field of the earth at this point is no longer affected by the pipe temperature. This allows us to determine the horizontal influence range of the temperature field around the pipe that transports the heat medium.

[0051] The temperature field of the finite element model is discretized into M layers of mesh from the upper boundary to the lower boundary. The collection of temperatures at each point within a material system is called the temperature field. It is a function of time and spatial coordinates, reflecting the spatial and temporal distribution of temperature.

[0052] For example, starting in October each year, simulate 152 days of temperature changes. Figure 2 This is a schematic diagram of the temperature field changes on day 12. Figure 3 This is a schematic diagram of the temperature field changes on day 152. From... Figure 2 and Figure 3 As can be seen, the temperature at different times and depths is constantly changing, and the position of the 0-degree isotherm in the temperature field gradually shifts downward over time.

[0053] In this invention, boundary conditions are applied to the finite element model. By solving the model, the time it takes for the freezing front to reach the lower boundary of each grid in the M-layer grid from the upper boundary to the lower boundary is obtained within the time range of the surface temperature function, as well as the temperature gradient at the location of the freezing front. At the same time, the depth of frozen soil is obtained.

[0054] A freezing front is a movable contact interface between frozen and unfrozen soil, specifically the 0°C interface. During the freezing process, the freezing front moves. Based on finite element analysis results, for each of the M-layer meshes, the time it takes for the freezing front to reach the lower edge of each mesh layer is determined, and the temperature gradient at that time is calculated. When meshing, M should be sufficiently large to minimize computational errors.

[0055] It should be noted that the freezing front will first appear at the upper boundary of the computational domain, then move towards the lower boundary over time. After reaching a certain depth, it will stop moving downwards and instead move upwards; this depth is the frost depth. Therefore, it is only necessary to find the grid cells containing the freezing front and calculate the temperature gradient. Figure 4 The figure shows the location of the freezing front and the temperature gradient at the freezing front in the temperature field. As can be seen from the figure, the location of the freezing front gradually shifts downwards over time, with the frozen soil depth reaching approximately 4.25m.

[0056] In this invention, determining the frost heave amount based on time and temperature gradients includes:

[0057] The frost heave amount is determined using the following formula based on the time and temperature gradient;

[0058]

[0059] Among them, h f P represents the frost heave. e The external load is represented by a and b, which are soil-related parameters, gradT(t) is the temperature gradient, and V is the temperature gradient. t This represents the change over time.

[0060] For example, starting in October each year, the temperature changes over 152 days were simulated. Based on the temperature field calculations of the finite element model, Table 1 lists the locations of the freezing fronts at different times and calculates the surface frost heave. According to Table 1, starting from day 0 of the surface temperature function, the maximum frost heave is reached on day 152.

[0061] Table 1 shows the location of the freezing front and the corresponding surface frost heave at different times.

[0062]

[0063]

[0064] In this invention, after obtaining the temperature field and frost heave results, the safety of oil and gas pipelines under frozen soil conditions can be evaluated by simulating the stress on the pipeline throughout the frost heave process, thus providing support for pipeline design and operation.

[0065] The method for obtaining frost heave in permafrost regions provided by this invention constructs a finite element model of the calculation area of ​​the permafrost region based on the surface temperature function, a first depth, and a constant temperature. Boundary conditions are applied to the finite element model to simulate the temperature field changes in the calculation area within the calculation time range of the surface temperature function, and the changes in frost heave depth and frost heave amount over time are obtained. Based on these changes, the maximum frost heave amount can be quickly and accurately determined, providing basic data for the safety assessment and design of buried pipelines. This method is low in cost, easy to operate, quick to calculate, and accurate in results.

[0066] The apparatus for obtaining frost heave in permafrost regions provided by the present invention is described below. The apparatus for obtaining frost heave in permafrost regions described below can be referred to in correspondence with the method for obtaining frost heave in permafrost regions described above.

[0067] Figure 5 A schematic flowchart of a device for obtaining frost heave in permafrost regions provided by the present invention is shown below. Figure 5 The device includes a fitting module 51, an acquisition module 52, a construction module 53, a calculation module 54, and a determination module 55, wherein:

[0068] The fitting module is used to determine the surface temperature function;

[0069] The acquisition module is used to determine the constant temperature at the first depth below the Earth's surface.

[0070] The construction module is used to construct a finite element model of the calculation area of ​​the permafrost region based on the surface temperature function, the first depth, and the constant temperature. The upper boundary of the finite element model is the surface, and the lower boundary is the first depth. The principle for determining the left and right boundaries is that the temperature gradient at the left and right boundaries is less than K℃ / m, where K is a preset value. Within the calculation area of ​​the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary.

[0071] The calculation module is used to apply boundary conditions to the finite element model and obtain the time from the surface temperature function to the bottom of each grid in the M-layer grid from the surface to the bottom, as well as the temperature gradient at the location of the freezing front.

[0072] The determination module is used to determine the amount of frost heave based on time and temperature gradients.

[0073] In a further embodiment of the above-mentioned device, the surface temperature function includes:

[0074]

[0075] Where T(x,0,t) is the ambient temperature, T m Let A be the annual average atmospheric temperature, A be the amplitude of atmospheric temperature change, and t be time.

[0076] In a further embodiment of the above-described device, the boundary conditions include: the left and right boundaries are adiabatic, the lower boundary is at a constant temperature, and the upper boundary is applied based on the surface temperature function.

[0077] In a further embodiment of the above apparatus, the determining module is specifically used for:

[0078] The amount of frost heave is determined using the following formula based on the time and the temperature gradient;

[0079]

[0080] Among them, h f P represents the frost heave. e The external load is represented by a and b, which are soil-related parameters, gradT(t) is the temperature gradient, and V is the temperature gradient. t This represents the change over time.

[0081] Since the device described in this embodiment of the invention is based on the same principle as the method described in the above embodiments, more detailed explanations will not be repeated here.

[0082] It should be noted that, in the embodiments of the present invention, the relevant functional modules can be implemented by a hardware processor.

[0083] The device for obtaining frost heave in permafrost regions provided by this invention constructs a finite element model of the calculation area of ​​the permafrost region based on the surface temperature function, a first depth, and a constant temperature. Boundary conditions are applied to the finite element model to simulate the temperature field changes in the calculation area within the calculation time range of the surface temperature function, and the changes in frost heave depth and frost heave amount over time are obtained. Based on these changes, the maximum frost heave amount can be quickly and accurately determined, providing basic data for the safety assessment and design of buried pipelines. It is low in cost, easy to operate, quick to calculate, and accurate in results.

[0084] Figure 6 An example is a schematic diagram of the physical structure of an electronic device, such as... Figure 6As shown, the electronic device may include: a processor 61, a communication interface 62, a memory 63, and a communication bus 64, wherein the processor 61, the communication interface 62, and the memory 63 communicate with each other through the communication bus 64. The processor 61 can call the computer program in the memory 63 to execute the steps of the method for obtaining the frost heave amount in the permafrost region, such as: determining the surface temperature function; determining the constant temperature at the first depth below the surface; constructing a finite element model of the calculation area of ​​the permafrost region based on the surface temperature function, the first depth, and the constant temperature; wherein the upper boundary of the finite element model is the surface, and the lower boundary is the first depth; the principle for determining the left and right boundaries is that the temperature gradient level change at the left and right boundaries is less than K℃ / m, where K is a preset value; within the calculation area of ​​the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary; boundary conditions are applied to the finite element model to obtain the time from the surface to the lower boundary of each layer of mesh in the M layers of mesh within the time range of the surface temperature function, as well as the temperature gradient at the location of the freezing front; and the frost heave amount is determined based on the time and temperature gradient.

[0085] Furthermore, the logical instructions in the aforementioned memory 63 can be implemented as software functional units and, when sold or used as independent products, can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0086] On the other hand, the present invention also provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions, and when the program instructions are executed by the computer, the computer can perform the steps of a method for obtaining frost heave in permafrost regions, such as: determining a surface temperature function; determining a constant temperature at a first depth below the surface; constructing a finite element model of the calculation area of ​​the permafrost region based on the surface temperature function, the first depth, and the constant temperature; wherein the upper boundary of the finite element model is the surface, and the lower boundary is the first depth; the principle for determining the left and right boundaries is that the temperature gradient level change at the left and right boundaries is less than K℃ / m, where K is a preset value; within the calculation area of ​​the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary; boundary conditions are applied to the finite element model to obtain the time from the surface to the lower boundary of each layer of mesh in the M layers within the time range of the surface temperature function, as well as the temperature gradient at the location of the freezing front; and the frost heave is determined based on the time and temperature gradient.

[0087] On the other hand, embodiments of the present invention also provide a processor-readable storage medium storing a computer program for causing the processor to execute steps of a method for obtaining frost heave in permafrost regions, such as: determining a surface temperature function; determining a constant temperature at a first depth below the surface; constructing a finite element model of a calculation region for the permafrost region based on the surface temperature function, the first depth, and the constant temperature; wherein the upper boundary of the finite element model is the surface, and the lower boundary is the first depth; the principle for determining the left and right boundaries is that the temperature gradient level change at the left and right boundaries is less than K℃ / m, where K is a preset value; within the calculation region of the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary; boundary conditions are applied to the finite element model to obtain the time from the surface to the lower boundary of each layer of mesh in the M layers within the time range of the surface temperature function, as well as the temperature gradient at the location of the freezing front; and the frost heave is determined based on the time and temperature gradient.

[0088] The processor-readable storage medium can be any available medium or data storage device that the processor can access, including but not limited to magnetic memory (e.g., floppy disk, hard disk, magnetic tape, magneto-optical disk (MO)), optical memory (e.g., CD, DVD, BD, HVD), and semiconductor memory (e.g., ROM, EPROM, EEPROM, non-volatile memory (NAND FLASH), solid-state drive (SSD)).

[0089] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs. Those skilled in the art can understand and implement this without any creative effort.

[0090] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus necessary general-purpose hardware platforms, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments.

[0091] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for obtaining frost heave in permafrost regions, characterized in that, include: Determine the surface temperature function; Determine the constant temperature at the first depth below the Earth's surface; A finite element model of the calculation area of ​​the permafrost region is constructed based on the surface temperature function, the first depth, and the constant temperature; wherein, the upper boundary of the finite element model is the surface, and the lower boundary is the first depth; the principle for determining the left and right boundaries is that the temperature gradient level change at the left and right boundaries is less than K℃ / m, where K is a preset value; within the calculation area of ​​the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary; Boundary conditions are applied to the finite element model to obtain the time it takes for the freezing front to reach the bottom of each grid layer in the M-layer grid from the surface to the bottom within the time range of the surface temperature function, as well as the temperature gradient at the location of the freezing front. The amount of frost heave is determined based on the time and the temperature gradient. The step of determining the amount of frost heave based on the time and the temperature gradient includes: The amount of frost heave is determined using the following formula based on the time and the temperature gradient; ； in, This is the amount of frost heave. The load is the external load, and a and b are parameters related to the soil properties. For temperature gradient, This represents the change over time.

2. The method for obtaining frost heave in permafrost areas according to claim 1, characterized in that, The surface temperature function includes: ； in, Atmospheric ambient temperature The annual average temperature of the atmosphere. This represents the amplitude of atmospheric temperature change. For time.

3. The method for obtaining frost heave in permafrost areas according to claim 1, characterized in that, The boundary conditions include: the left and right boundaries are adiabatic, the lower boundary is at a constant temperature, and the upper boundary is applied based on the surface temperature function.

4. A device for obtaining frost heave in permafrost regions, characterized in that, include: The fitting module is used to determine the surface temperature function; The acquisition module is used to determine the constant temperature at the first depth below the Earth's surface. The construction module is used to construct a finite element model of the permafrost region calculation area based on the surface temperature function, the first depth, and the constant temperature; wherein, the upper boundary of the finite element model is the surface, and the lower boundary is the first depth; the principle for determining the left and right boundaries is that the temperature gradient level change at the left and right boundaries is less than K℃ / m, where K is a preset value; within the calculation area of ​​the finite element model, it is discretized into M layers of mesh from the upper boundary to the lower boundary; The calculation module is used to apply boundary conditions to the finite element model and obtain the time from the surface temperature function to the bottom of each grid in the M-layer grid from the surface to the bottom, as well as the temperature gradient at the location of the freezing front. The determining module is used to determine the amount of frost heave based on the time and the temperature gradient; The determining module is specifically used for: The amount of frost heave is determined using the following formula based on the time and the temperature gradient; ； in, This is the amount of frost heave. The load is the external load, and a and b are parameters related to the soil properties. For temperature gradient, This represents the change over time.

5. The device for obtaining frost heave in permafrost areas according to claim 1, characterized in that, The surface temperature function includes: ； in, Atmospheric ambient temperature The annual average temperature of the atmosphere. This represents the amplitude of atmospheric temperature change. For time.

6. The device for obtaining frost heave in permafrost areas according to claim 1, characterized in that, The boundary conditions include: the left and right boundaries are adiabatic, the lower boundary is at a constant temperature, and the upper boundary is applied based on the surface temperature function.

7. An electronic device comprising a processor and a memory storing a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method for obtaining frost heave in permafrost regions as described in any one of claims 1 to 3.

8. A processor-readable storage medium, characterized in that, The processor-readable storage medium stores a computer program for causing the processor to perform the steps of the method for obtaining frost heave in permafrost regions according to any one of claims 1 to 3.

Images