A temperature field calculation method, device, electronic device and storage medium thereof

Abstract

The present application relates to the technical field of dam safety control, and specifically to a temperature field calculation method, device, electronic device and storage medium thereof. The temperature field calculation method obtains the temperature boundary of one side of the predicted part; and establishes a heat conduction mathematical model of the predicted part according to the heat conduction equation; obtains the measured temperature of at least one characteristic point located at the predicted part; performs parameter inversion on the thermal conductivity of the predicted part according to the measured temperature of the characteristic point to obtain the thermal conductivity of the predicted part; and obtains the temperature field distribution of the predicted part according to the thermal conductivity of the predicted part and the heat conduction mathematical model of the predicted part. The temperature field calculation method obtains the temperature field boundary conditions and obtains the thermal conductivity of the predicted part through parameter inversion, so as to more accurately determine the temperature field distribution of the dam, and analyze the problems existing in the dam more accurately according to the temperature field distribution.

Description

Technical Field

[0001] The present application relates to the technical field of dam safety control, and in particular to a temperature field calculation method, device, electronic device and storage medium thereof. Background Art

[0002] In order to ensure the safety of the constructed dam, the bottom elevation of the high arch dam is often lower than the bottom elevation of the original riverbed. Therefore, when dealing with the foundation, the original riverbed in the dam site area needs to be excavated downward for dozens of meters. On the one hand, during the concrete pouring of the riverbed dam section, in order to prevent the dam body from cracking due to excessive axial stress on the upstream surface of the dam, it is necessary to insulate the upstream bottom, and the slag heap method can effectively ensure the temperature of the dam bottom; on the other hand, during the construction of the upstream dam, in order to meet the requirements of water retention, it is necessary to build a cofferdam upstream of the dam foundation. When the dam begins to store water, it is necessary to blast and dismantle the cofferdam. After the blasting, the aggregate of the cofferdam remains in the reservoir basin upstream of the dam, which increases the depth of the slag heap upstream of the dam.

[0003] At present, due to the lack of corresponding measured temperature data at the slag pile site, the bottom water temperature or ground temperature is usually used to calculate the temperature at the slag pile site during structural calculations. This creates a blind spot in dam structure analysis and makes it difficult to promptly discover potential problems in the dam. Summary of the invention

[0004] The purpose of the embodiments of the present application is to provide a temperature field calculation method, device, electronic device and storage medium thereof to improve the above-mentioned technical problems.

[0005] To achieve the above objectives, this application provides the following technical solutions:

[0006] In a first aspect, an embodiment of the present application provides a temperature field calculation method, the method comprising:

[0007] Obtain the temperature boundary on one side of the predicted part;

[0008] According to the temperature boundary on one side and the heat conduction equation Establishing a mathematical model of heat conduction of the predicted part; wherein a is the thermal conductivity of the predicted part;

[0009] Obtaining a measured temperature of at least one characteristic point; wherein the characteristic point is located at the predicted position;

[0010] Perform parameter inversion on the thermal conductivity of the predicted location according to the measured temperature of the characteristic point to obtain the thermal conductivity of the predicted location;

[0011] The temperature field distribution of the predicted part is obtained according to the thermal conductivity of the predicted part and the heat conduction mathematical model of the predicted part.

[0012] In the above method, a mathematical model of heat conduction of the predicted part is established; the measured temperature of the characteristic point located at the predicted part is obtained; the thermal conductivity of the predicted part is parameter inverted according to the measured temperature of the characteristic point to obtain the thermal conductivity of the predicted part. This temperature field calculation method can more accurately determine the temperature field distribution of the dam by obtaining accurate temperature field boundary conditions and obtaining the thermal conductivity of the predicted part through parameter inversion, and further, can more accurately analyze the problems existing in the dam.

[0013] Optionally, the predicted part is the slag pile part of the dam; the temperature boundary on one side of the predicted part is the upper water temperature boundary of the slag pile part; and the thermal conductivity of the predicted part is the thermal conductivity of the slag pile part.

[0014] In the above method, the upper water temperature boundary of the slag pile part can be obtained, and a heat conduction mathematical model of the slag pile part can be established according to the upper water temperature boundary of the slag pile part and the heat conduction equation, and the thermal conductivity coefficient of the slag pile part can be parameter inverted to obtain the temperature field distribution of the slag pile part.

[0015] Optionally, the predicted part is the dam foundation rock mass part; the temperature boundary on one side of the predicted part is the upper slag temperature boundary of the dam foundation rock mass part; and the thermal conductivity of the predicted part is the thermal conductivity of the dam foundation rock mass part.

[0016] In the above method, the upper slag temperature boundary of the dam foundation rock mass part can be obtained, and a heat conduction mathematical model of the dam foundation rock mass part can be established according to the upper slag temperature boundary of the dam foundation rock mass part and the heat conduction equation. The thermal conductivity coefficient of the dam foundation rock mass part can be parameter inverted to obtain the temperature field distribution of the dam foundation rock mass part.

[0017] Optionally, the thermal diffusivity of the predicted part is obtained by performing parameter inversion on the thermal diffusivity according to the measured temperature of the characteristic point, including: obtaining the trial temperature of the characteristic point according to the trial calculated thermal diffusivity of the predicted part and the mathematical model of heat conduction of the predicted part; wherein, the initial value of the trial calculated thermal diffusivity of the predicted part is the thermal diffusivity of a general rock mass; judging whether the trial calculated temperature satisfies a judgment condition |TS|<0.2; wherein T is the measured temperature of the characteristic point, and S is the trial calculated temperature; if not, adjusting the value of the trial calculated thermal diffusivity of the predicted part until the judgment condition is met; wherein, the trial calculated thermal diffusivity of the predicted part that meets the judgment condition is the thermal diffusivity of the predicted part.

[0018] In the above method, the judgment conditions can be adjusted according to the actual accuracy requirements. The thermal conductivity of the predicted part is obtained through parameter inversion, making the calculation process of the temperature field simpler.

[0019] Optionally, obtaining the temperature field distribution of the predicted part according to the thermal conductivity of the predicted part and the heat conduction mathematical model of the predicted part includes: obtaining the initial conditions and initial boundaries of the heat conduction mathematical model of the predicted part under general conditions; wherein, the initial condition is when t=0, T=0, 0≤x≤∞, and the initial boundary is when x=0, t>0, When x=∞, t>0, T=0, wherein A is the temperature variation amplitude of one side of the predicted part, and P is the temperature variation period of one side of the predicted part; the general conditions include: the calculation area of ​​the temperature field of the predicted part is a semi-infinite object, and the initial temperature of the semi-infinite object is 0°C; the solution of the heat conduction mathematical model of the predicted part under the general conditions is obtained according to Fourier transform Among them, ξ is an integral variable; the second term of the solution of the mathematical model of heat conduction of the predicted part eventually decays to 0 with the increase of time; the temperature field distribution of the predicted part is obtained according to the thermal conductivity of the predicted part and the first term of the solution of the mathematical model of heat conduction of the predicted part.

[0020] In the above method, the solution of the mathematical model of heat conduction is obtained by Fourier transform and the composition of the solution is analyzed. Since the second term of the solution eventually decays to 0 with the increase of time, the temperature field distribution of the predicted part is obtained by substituting the thermal conductivity coefficient of the predicted part into the first term of the solution.

[0021] Optionally, the temperature field calculation method also includes: obtaining the elevations of multiple fitting points in the elevation direction of the slag pile part; obtaining the temperature value corresponding to the elevation of the fitting point according to the temperature field distribution of the slag pile part; obtaining a temperature distribution fitting formula of the slag pile part according to the elevation of the fitting point and the temperature value corresponding to the elevation; wherein the temperature distribution fitting formula of the slag pile part uses the temperature value corresponding to the elevation as the horizontal coordinate and the elevation of the fitting point as the vertical coordinate.

[0022] In the above method, the temperature distribution fitting formula of the slag pile part is obtained according to the elevations of multiple fitting points and the temperatures at the fitting points, and the temperature value at a certain point of the slag pile part can be obtained more conveniently by using the formula. The temperature distribution fitting formula of the slag pile part can also be used to obtain the temperature field distribution of the slag pile of other dams.

[0023] Optionally, the temperature field calculation method also includes: obtaining air temperature boundaries, adiabatic boundaries, fixed temperature boundaries, thermal diffusivity of the dam body, thermal diffusivity of the water body, geothermal gradient of the dam foundation rock mass along the elevation direction, thermal diffusivity of the slag pile part, and thermal diffusivity of the dam foundation rock mass part; wherein the air temperature boundaries include the upper water temperature boundary of the slag pile part, the upper slag pile temperature boundary of the dam foundation rock mass part, the air temperature boundary above the water surface of the upstream surface of the dam body, the dam top air temperature boundary, and the air temperature boundary above the water surface of the downstream surface of the dam body; the adiabatic boundaries include the upstream adiabatic temperature boundary and the downstream adiabatic temperature boundary; the fixed temperature boundary refers to the bottom geothermal temperature boundary; the thermal diffusivity of the dam body is the thermal diffusivity of a general rock mass; and the temperature field of the slag pile-dam foundation rock mass-dam body-water body is obtained according to the air temperature boundary, the adiabatic boundary, the fixed temperature boundary, the thermal diffusivity of the slag pile part, the thermal diffusivity of the dam foundation rock mass part, the thermal diffusivity of the dam body, the thermal diffusivity of the water body, the geothermal gradient, and the heat conduction equation.

[0024] In the above method, by obtaining the boundary conditions of the overall temperature field of slag pile-dam foundation rock mass-dam body-water body, and building a mathematical model of the overall temperature field according to the boundary conditions of the overall temperature field, the temperature distribution in the slag pile-dam foundation rock mass-dam body-water body temperature field, that is, the overall temperature field of the dam, can be intuitively reflected.

[0025] In a second aspect, an embodiment of the present application provides a temperature field calculation device, the device comprising:

[0026] A first acquisition module is used to acquire a temperature boundary on one side of the predicted part;

[0027] A model building module is used to build a model based on the temperature boundary on one side and the heat conduction equation Establishing a mathematical model of heat conduction of the predicted part; wherein a is the thermal conductivity of the predicted part;

[0028] A second acquisition module is used to acquire the measured temperature of at least one characteristic point; wherein the characteristic point is located at the predicted position;

[0029] A parameter inversion module, used to perform parameter inversion on the thermal conductivity of the predicted part according to the measured temperature of the characteristic point, so as to obtain the thermal conductivity of the predicted part;

[0030] The third acquisition module is used to acquire the temperature field distribution of the predicted part according to the thermal conductivity coefficient of the predicted part and the heat conduction mathematical model of the predicted part.

[0031] In a third aspect, the present application provides an electronic device, comprising: a processor and a memory, wherein the memory stores machine-readable instructions executable by the processor, and when the machine-readable instructions are executed by the processor, a method as described in any of the above-mentioned temperature field calculation methods is performed.

[0032] In a fourth aspect, the present application provides a storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the method described in any one of the above-mentioned temperature field calculation methods is executed.

[0033] Other features and advantages of the present application will be described in the following description, and partly become apparent from the description, or be understood by practicing the embodiments of the present application. The purpose and other advantages of the present application can be realized and obtained by the structures specifically pointed out in the written description, claims, and drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the embodiments of the present application will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present application and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0035] Figure 1 A schematic diagram of a flow chart of a temperature field calculation method provided in an embodiment of the present application;

[0036] Figure 2 A schematic diagram of a dam model provided in an embodiment of the present application;

[0037] Figure 3 A schematic diagram of the structure of a temperature field calculation device provided in an embodiment of the present application;

[0038] Figure 4 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0039] The technical solutions in the embodiments of the present application will be described below in conjunction with the drawings in the embodiments of the present application. It should be noted that similar reference numerals and letters represent similar items in the following drawings, so once an item is defined in one drawing, it does not need to be further defined and explained in the subsequent drawings.

[0040] The terms "comprises," "comprising," or any other variation thereof are intended to encompass non-exclusive inclusion, such that a process, method, article, or apparatus that includes a list of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or apparatus. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article, or apparatus that includes the element.

[0041] The terms "first", "second", etc. are only used to distinguish one entity or operation from another entity or operation, and should not be understood as indicating or implying relative importance, nor should they be understood as requiring or implying any such actual relationship or order between these entities or operations.

[0042] In view of the deficiencies in the prior art, the embodiments of the present application provide a temperature field calculation method, device, electronic device and storage medium for more accurately obtaining the temperature field distribution of the dam.

[0043] Please refer to Figure 1 , Figure 1 A flow chart of a temperature field calculation method provided in an embodiment of the present application, the temperature field calculation method is used to more accurately obtain the temperature field calculation result of the dam, including the following steps:

[0044] Step 101: Obtain the temperature boundary on one side of the predicted location.

[0045] Step 102: According to the temperature boundary on one side and the heat conduction equation A mathematical model of heat conduction at the predicted location is established.

[0046] Step 103: Obtain the measured temperature of at least one characteristic point; wherein the characteristic point is located at the predicted position.

[0047] Step 104: Perform parameter inversion on the thermal conductivity of the predicted location according to the measured temperature of the characteristic point to obtain the thermal conductivity of the predicted location.

[0048] Step 105: Obtain the temperature field distribution of the predicted part according to the thermal conductivity of the predicted part and the heat conduction mathematical model of the predicted part.

[0049] Wherein, in step 101, the predicted part refers to any component part of the dam, the thermal conductivity of the same predicted part is the same, and the temperature boundary on one side refers to the temperature value at the boundary on one side of the predicted part.

[0050] In step 102, a is the thermal conductivity of the predicted part, T represents temperature, t represents time, and x is a spatial variable, which here represents the distance to the boundary of one side of the predicted part.

[0051] Wherein, in step 103, the actual temperature corresponding to at least one characteristic point in the predicted part is obtained.

[0052] Among them, in step 104, the thermal conductivity of the predicted position is obtained through parameter inversion, which simplifies the calculation process.

[0053] In step 105, the temperature field distribution of the predicted part can be determined by substituting the thermal conductivity of the predicted part into the heat conduction mathematical model of the predicted part.

[0054] As can be seen from the above, the embodiment of the present application provides a temperature field calculation method, which establishes a mathematical model of heat conduction at the predicted part; obtains the measured temperature of the characteristic point located at the predicted part; performs parameter inversion on the thermal conductivity of the predicted part according to the measured temperature of the characteristic point to obtain the thermal conductivity of the predicted part. The temperature field calculation method can more accurately determine the temperature field distribution of the dam by obtaining accurate temperature field boundary conditions and obtaining the thermal conductivity of the predicted part through parameter inversion.

[0055] In some optional embodiments, the predicted part is the slag pile part of the dam; the temperature boundary on one side of the predicted part is the upper water temperature boundary of the slag pile part; and the thermal conductivity of the predicted part is the thermal conductivity of the slag pile part.

[0056] The upper water temperature boundary of the slag pile part refers to the water temperature at the contact surface between the upper boundary of the slag pile part and the water body; the upper water temperature boundary of the slag pile part can be obtained, according to the upper water temperature boundary of the slag pile part and the heat conduction equation Establish a mathematical model of heat conduction of the slag pile part; wherein a1 is the thermal conductivity of the slag pile part; obtain the measured temperature of at least one first characteristic point; wherein the first characteristic point is located in the slag pile part; perform parameter inversion on the thermal conductivity of the slag pile part according to the measured temperature of the first characteristic point to obtain the thermal conductivity of the slag pile part; obtain the temperature field distribution of the slag pile part according to the thermal conductivity of the slag pile part and the mathematical model of heat conduction of the slag pile part. Since thermometer measuring points are set in the typical dam section of the dam, each measuring point is provided with a thermometer 10 cm away from the upstream dam surface. After the dam is filled with water to the normal water level, the temperature at the thermometer measuring point changes periodically every year, and the temperature change amplitude at the thermometer measuring points at different water depths is also different. The average value of the stable measured values ​​of the thermometers at the same elevation above the slag pile part can be selected as the upper water temperature boundary of the slag pile part.

[0057] In some optional embodiments, the predicted part is the dam foundation rock part; the temperature boundary on one side of the predicted part is the upper slag temperature boundary of the dam foundation rock part; and the thermal conductivity of the predicted part is the thermal conductivity of the dam foundation rock part.

[0058] The upper slag temperature boundary of the dam foundation rock mass part refers to the slag temperature at the contact surface between the upper boundary of the dam foundation rock mass part and the slag part; the upper slag temperature boundary of the dam foundation rock mass part can be obtained, and according to the upper slag temperature boundary of the dam foundation rock mass part and the heat conduction equation Establish a mathematical model of heat conduction of the dam foundation rock mass; wherein a2 is the thermal conductivity of the dam foundation rock mass; obtain the measured temperature of at least one second characteristic point; wherein the second characteristic point is located in the dam foundation rock mass; perform parameter inversion on the thermal conductivity of the dam foundation rock mass according to the measured temperature of the second characteristic point to obtain the thermal conductivity of the dam foundation rock mass; obtain the temperature field distribution of the dam foundation rock mass according to the thermal conductivity of the dam foundation rock mass and the mathematical model of heat conduction of the dam foundation rock mass. A piezometer will be installed about 10m below the heel of the dam to detect changes in osmotic pressure. The piezometer will have a temperature sensitive element. The temperature value can be converted by selecting the temperature value attached to the piezometer. The temperature values ​​converted by multiple thermometers at the same height are averaged. The average value is the slag temperature boundary of the upper part of the dam foundation rock mass.

[0059] In some optional embodiments, step 104 specifically includes: obtaining the trial temperature of the characteristic point according to the trial thermal conductivity of the predicted part and the mathematical model of heat conduction of the predicted part; judging whether the trial temperature satisfies the judgment condition |TS|<0.2; if not, adjusting the value of the trial thermal conductivity of the predicted part until the judgment condition is met; wherein the trial thermal conductivity of the predicted part that meets the judgment condition is the thermal conductivity of the predicted part.

[0060] Among them, the initial value of the calculated thermal conductivity of the predicted part is the thermal conductivity of the general rock mass. In the specific implementation of the present application, the thermal conductivity of the general rock mass is 185 kJ / (m·d·℃). In the absence of engineering site test parameters, it can also be taken as 200 kJ / (m·d·℃) or other thermal conductivity values ​​close to 200 kJ / (m·d·℃); T is the measured temperature of the characteristic point, and S is the calculated temperature; the judgment conditions can also be adjusted so that the obtained thermal conductivity of the predicted part meets different calculation accuracy requirements.

[0061] In some optional embodiments, step 105 specifically includes: obtaining the initial conditions and initial boundaries of the heat conduction mathematical model of the predicted part under general conditions; wherein the initial conditions are when t=0, T=0, 0≤x≤∞, and the initial boundaries are when x=0, t>0, When x=∞, t>0, T=0, wherein A is the temperature variation amplitude of one side of the predicted part, and P is the temperature variation period of one side of the predicted part; the general conditions include: the calculation area of ​​the temperature field of the predicted part is a semi-infinite object, and the initial temperature of the semi-infinite object is 0°C; the solution of the heat conduction mathematical model of the predicted part under the general conditions is obtained according to Fourier transform Among them, ξ is an integral variable; the second term of the solution of the mathematical model of heat conduction of the predicted part eventually decays to 0 with the increase of time; the temperature field distribution of the predicted part is obtained according to the thermal conductivity of the predicted part and the first term of the solution of the mathematical model of heat conduction of the predicted part.

[0062] Among them, since the second term of the solution of the mathematical model of heat conduction in the predicted part contains , will gradually decay to 0 as time goes by, leaving The remaining term is called the stable temperature field. From this stable temperature field term, we can see that the temperature of the predicted part changes periodically, and the change period is the same as the change period of the temperature at the boundary of the predicted part, which is P. When , the temperature variation inside the predicted part can be obtained as The semi-infinite object refers to an object that has the x=0 plane as its only interface and extends infinitely in the x direction.

[0063] In some optional embodiments, the method further includes: obtaining the elevations of multiple fitting points in the elevation direction of the slag pile part; obtaining the temperature value corresponding to the elevation of the fitting point according to the temperature field distribution of the slag pile part; obtaining a temperature distribution fitting formula of the slag pile part according to the elevation of the fitting point and the temperature value corresponding to the elevation; wherein the temperature distribution fitting formula of the slag pile part uses the temperature value corresponding to the elevation as the horizontal coordinate and the elevation of the fitting point as the vertical coordinate.

[0064] Wherein, the multiple is at least two. For similar dams (high arch dams), the temperature field distribution of the slag pile part is similar. The fitting formula can be obtained based on the data in the calculation examples of the slag pile part in multiple projects. After comparing the calculation results of multiple projects, this application has obtained a formula with a relatively high fitting degree: Where ΔT represents the temperature difference between the upper water temperature boundary of the slag pile part and the upper slag temperature boundary of the dam foundation rock mass, that is, the temperature difference between the upper boundary of the slag pile part and the lower boundary of the slag pile part, T1 represents the upper slag temperature boundary of the dam foundation rock mass, y b Indicates the elevation value of the bottom of the dam.

[0065] In some optional embodiments, the method further includes: obtaining air temperature boundaries, adiabatic boundaries, fixed temperature boundaries, thermal conductivity of the dam body, thermal conductivity of the water body, geothermal gradient of the dam foundation rock mass along the elevation direction, thermal conductivity of the slag pile part, and thermal conductivity of the dam foundation rock mass part; wherein the air temperature boundaries include the upper water temperature boundary of the slag pile part, the upper slag pile temperature boundary of the dam foundation rock mass part, the air temperature boundary above the water surface of the upstream surface of the dam body, the dam top air temperature boundary, and the air temperature boundary above the water surface of the downstream surface of the dam body; the adiabatic boundaries include the upstream adiabatic temperature boundary and the downstream adiabatic temperature boundary; the fixed temperature boundary refers to the bottom geothermal temperature boundary; the thermal conductivity of the dam body is the thermal conductivity of a general rock mass; and the temperature field of the slag pile-dam foundation rock mass-dam body-water body is obtained according to the air temperature boundary, the adiabatic boundary, the fixed temperature boundary, the thermal conductivity of the slag pile part, the thermal conductivity of the dam foundation rock mass part, the thermal conductivity of the dam body, the thermal conductivity of the water body, the geothermal gradient, and the heat conduction equation.

[0066] in, Figure 2 This is a schematic diagram of a dam model provided in the embodiment of the present application. For understanding of the relevant parameters for establishing the model, please refer to Figure 2 Since such an idealized object with infinite size does not exist in engineering practice, in actual calculations, adiabatic boundaries can be set at both ends of the predicted part. On this basis, numerical algorithms (finite element method or finite difference method) are used to solve specific practical problems. Among them, the fixed temperature boundary can be obtained by finding the measured temperature analysis results of the ground temperature in the previous hydrogeological data, calculating the geothermal gradient of the dam foundation rock mass along the elevation direction through the previous data, and inferring the deep temperature of the dam foundation rock mass. The deep temperature of the dam foundation rock mass is the fixed temperature boundary.

[0067] Please refer to Figure 3 , Figure 3 A schematic diagram of the structure of a temperature field calculation device provided in an embodiment of the present application, the temperature field calculation device is used to more accurately obtain the temperature field distribution of the dam, including:

[0068] A first acquisition module is used to acquire a temperature boundary on one side of the predicted part;

[0069] A model building module is used to build a model based on the temperature boundary on one side and the heat conduction equation Establishing a mathematical model of heat conduction of the predicted part; wherein a is the thermal conductivity of the predicted part;

[0070] A second acquisition module is used to acquire the measured temperature of at least one characteristic point; wherein the characteristic point is located at the predicted position;

[0071] A parameter inversion module, used to perform parameter inversion on the thermal conductivity of the predicted part according to the measured temperature of the characteristic point, so as to obtain the thermal conductivity of the predicted part;

[0072] The third acquisition module is used to acquire the temperature field distribution of the predicted part according to the thermal conductivity coefficient of the predicted part and the heat conduction mathematical model of the predicted part.

[0073] In some optional embodiments, the first acquisition module is further used to: acquire the upper water temperature boundary of the slag pile part; the model building module is further used to: according to the upper water temperature boundary of the slag pile part and the heat conduction equation A mathematical model of heat conduction of the slag pile part is established; wherein a1 is the thermal conductivity of the slag pile part; the second acquisition module is also used to: obtain the measured temperature of at least one first characteristic point; wherein the first characteristic point is located in the slag pile part; the parameter inversion module is also used to: perform parameter inversion on the thermal conductivity of the slag pile part according to the measured temperature of the first characteristic point, and obtain the thermal conductivity of the slag pile part; the third acquisition module is also used to: obtain the temperature field distribution of the slag pile part according to the thermal conductivity of the slag pile part and the mathematical model of heat conduction of the slag pile part.

[0074] In some other optional embodiments, the first acquisition module is further used to: obtain the upper slag temperature boundary of the dam foundation rock mass part; the model building module is further used to: according to the upper slag temperature boundary of the dam foundation rock mass part and the heat conduction equation A mathematical model of heat conduction of the dam foundation rock mass part is established; wherein a2 is the thermal conductivity of the dam foundation rock mass part; the second acquisition module is also used to: obtain the measured temperature of at least one second characteristic point; wherein the second characteristic point is located in the dam foundation rock mass part; the parameter inversion module is also used to: perform parameter inversion on the thermal conductivity of the dam foundation rock mass part according to the measured temperature of the second characteristic point, and obtain the thermal conductivity of the dam foundation rock mass part; the third acquisition module is also used to: obtain the temperature field distribution of the dam foundation rock mass part according to the thermal conductivity of the dam foundation rock mass part and the mathematical model of heat conduction of the dam foundation rock mass part.

[0075] In some optional embodiments, the above-mentioned parameter inversion module is specifically used to: obtain the trial temperature of the characteristic point according to the trial thermal conductivity of the predicted part and the mathematical model of heat conduction of the predicted part; wherein, the initial value of the trial thermal conductivity of the predicted part is the thermal conductivity of the general rock mass; judge whether the trial temperature satisfies the judgment condition |TS|<0.2; wherein T is the measured temperature of the characteristic point, and S is the trial temperature; if not, adjust the value of the trial thermal conductivity of the predicted part until the judgment condition is met; wherein, the trial thermal conductivity of the predicted part that meets the judgment condition is the thermal conductivity of the predicted part.

[0076] In some optional embodiments, the third acquisition module is specifically used to: obtain the initial conditions and initial boundaries of the heat conduction mathematical model of the predicted part under general conditions; wherein the initial condition is when t=0, T=0, 0≤x≤∞, and the initial boundary is when x=0, t>0, When x=∞, t>0, T=0, wherein A is the temperature variation amplitude of one side of the predicted part, and P is the temperature variation period of one side of the predicted part; the general conditions include: the calculation area of ​​the temperature field of the predicted part is a semi-infinite object, and the initial temperature of the semi-infinite object is 0°C; the solution of the heat conduction mathematical model of the predicted part under the general conditions is obtained according to Fourier transform Among them, ξ is an integral variable; the second term of the solution of the mathematical model of heat conduction of the predicted part eventually decays to 0 with the increase of time; the temperature field distribution of the predicted part is obtained according to the thermal conductivity of the predicted part and the first term of the solution of the mathematical model of heat conduction of the predicted part.

[0077] In some optional embodiments, the temperature field calculation device further includes: a fourth acquisition module, which is used to obtain the temperature boundary, the adiabatic boundary, the fixed temperature boundary, the thermal conductivity of the dam body, the thermal conductivity of the water body, the geothermal gradient of the dam foundation rock mass part along the elevation direction, the thermal conductivity of the slag pile part, and the thermal conductivity of the dam foundation rock mass part; wherein the temperature boundary includes the upper water temperature boundary of the slag pile part, the upper slag pile temperature boundary of the dam foundation rock mass part, the temperature boundary above the water surface of the upstream surface of the dam body, the temperature boundary of the dam top, and the temperature boundary above the water surface of the downstream surface of the dam body The adiabatic boundary includes an upstream adiabatic temperature boundary and a downstream adiabatic temperature boundary, the fixed temperature boundary refers to the bottom geothermal temperature boundary; the thermal conductivity of the dam body is the thermal conductivity of a general rock mass; the overall temperature field acquisition module is used to obtain the temperature field of slag pile-dam foundation rock mass-dam body-water body according to the air temperature boundary, the adiabatic boundary, the fixed temperature boundary, the thermal conductivity of the slag pile part, the thermal conductivity of the dam foundation rock mass part, the thermal conductivity of the dam body, the thermal conductivity of the water body, the geothermal gradient and the heat conduction equation.

[0078] The specific implementation of each module in the above-mentioned temperature field calculation device can refer to the implementation method of the corresponding steps in the temperature field calculation method introduced above.

[0079] Figure 4 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. Figure 4 The electronic device 3 includes: a processor 301 and a memory 302. These components are interconnected and communicate with each other through a communication bus 303 and / or other forms of connection mechanisms (not shown).

[0080] The memory 302 includes one or more (only one is shown in the figure), which may be, but not limited to, a random access memory (RAM), a read only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), etc. The processor 301 and other possible components may access the memory 302 and read and / or write data therein.

[0081] The processor 301 includes one or more (only one is shown in the figure), which can be an integrated circuit chip with signal processing capabilities. The above-mentioned processor 301 can be a general-purpose processor, including a central processing unit (CPU), a micro control unit (MCU), a network processor (NP) or other conventional processors; it can also be a dedicated processor, including a neural network processor (NPU), a graphics processor (GPU), a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components. In addition, when there are multiple processors 301, some of them can be general-purpose processors and the other part can be dedicated processors.

[0082] One or more computer program instructions may be stored in the memory 302 , and the processor 301 may read and execute these computer program instructions to implement a temperature field calculation method provided in an embodiment of the present application.

[0083] Understandably, Figure 4 The structure shown is for illustration only. The electronic device 3 may also include Figure 4 More or fewer components as shown, or with Figure 4 Different structures are shown. Figure 4 The components shown in can be implemented by hardware, software or a combination thereof. The electronic device 3 may be a physical device, such as a PC, a laptop, a tablet computer, a mobile phone, a server, an embedded device, etc., or a virtual device, such as a virtual machine, a virtualized container, etc. Moreover, the electronic device 3 is not limited to a single device, but may also be a combination of multiple devices or a cluster consisting of a large number of devices.

[0084] The present application also provides a computer-readable storage medium, which stores computer program instructions. When the computer program instructions are read and executed by a computer processor, the temperature field calculation method provided in the present application is executed. For example, the computer-readable storage medium can be implemented as Figure 4 The memory 302 in the electronic device 3.

[0085] The above description is only an embodiment of the present application and is not intended to limit the protection scope of the present application. For those skilled in the art, the present application may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the protection scope of the present application.

Claims

1. A temperature field calculation method, characterized in that: include: Obtain the temperature boundary on one side of the predicted part; According to the temperature boundary on one side and the heat conduction equation Establishing a mathematical model of heat conduction of the predicted part; wherein a is the thermal conductivity of the predicted part; Obtaining a measured temperature of at least one characteristic point; wherein the characteristic point is located at the predicted position; Perform parameter inversion on the thermal conductivity of the predicted location according to the measured temperature of the characteristic point to obtain the thermal conductivity of the predicted location; Acquire the temperature field distribution of the predicted part according to the thermal conductivity of the predicted part and the heat conduction mathematical model of the predicted part; The predicted part is the slag dump part of the dam; the method specifically comprises: obtaining the upper water temperature boundary of the slag dump part; selecting the average value of the stable measured values ​​of the thermometers at the same elevation above the slag dump part as the upper water temperature boundary of the slag dump part; according to the upper water temperature boundary of the slag dump part and the heat conduction equation Establish a mathematical model of heat conduction of the slag pile part; wherein a1 is the thermal conductivity of the slag pile part; obtain the measured temperature of at least one first characteristic point; wherein the first characteristic point is located in the slag pile part; perform parameter inversion on the thermal conductivity of the slag pile part according to the measured temperature of the first characteristic point to obtain the thermal conductivity of the slag pile part; obtain the temperature field distribution of the slag pile part according to the thermal conductivity of the slag pile part and the mathematical model of heat conduction of the slag pile part; or, The predicted part is the rock mass of the dam foundation; the method specifically comprises: obtaining the upper slag temperature boundary of the dam foundation rock mass; wherein, selecting the temperature value obtained by converting the attached temperature degree of the piezometer, and averaging the temperature values ​​obtained by converting multiple thermometers at the same height, and taking the average value as the upper slag temperature boundary of the dam foundation rock mass; according to the upper slag temperature boundary of the dam foundation rock mass and the heat conduction equation Establish a mathematical model of heat conduction of the dam foundation rock mass part; wherein a2 is the thermal conductivity of the dam foundation rock mass part; obtain the measured temperature of at least one second characteristic point; wherein the second characteristic point is located in the dam foundation rock mass part; perform parameter inversion on the thermal conductivity of the dam foundation rock mass part according to the measured temperature of the second characteristic point to obtain the thermal conductivity of the dam foundation rock mass part; obtain the temperature field distribution of the dam foundation rock mass part according to the thermal conductivity of the dam foundation rock mass part and the mathematical model of heat conduction of the dam foundation rock mass part.

2. The method according to claim 1, characterized in that The performing parameter inversion on the thermal conductivity of the predicted part according to the measured temperature of the characteristic point to obtain the thermal conductivity of the predicted part includes: The trial temperature of the characteristic point is obtained according to the trial thermal conductivity of the predicted part and the heat conduction mathematical model of the predicted part; wherein the initial value of the trial thermal conductivity of the predicted part is the thermal conductivity of a general rock mass; Determine whether the trial temperature satisfies a determination condition |TS|<0.2; wherein T is the measured temperature of the characteristic point, and S is the trial temperature; If not, adjust the value of the calculated thermal conductivity of the predicted part until the judgment condition is met; Among them, the calculated thermal conductivity of the predicted part that meets the judgment condition is the thermal conductivity of the predicted part.

3. The method according to claim 1, characterized in that: The obtaining of the temperature field distribution of the predicted part according to the thermal conductivity of the predicted part and the heat conduction mathematical model of the predicted part comprises: Obtaining the initial conditions and initial boundaries of the heat conduction mathematical model of the predicted part under general conditions; wherein the initial condition is when t=0, T=0, 0≤x≤∞, and the initial boundary is when x=0, t>0, When x=∞, t>0, T=0, wherein A is the temperature variation amplitude of one side of the predicted part, and P is the temperature variation period of one side of the predicted part; the general conditions include: the calculation area of ​​the temperature field of the predicted part is a semi-infinite object, and the initial temperature of the semi-infinite object is 0°C; The solution of the mathematical model of heat conduction in the predicted position under the general conditions is obtained by Fourier transform. Wherein, ξ is an integral variable; the second term of the solution of the heat conduction mathematical model of the predicted location eventually decays to 0 as time increases; The temperature field distribution of the predicted part is obtained according to the thermal conductivity of the predicted part and the first item of the solution of the mathematical model of heat conduction of the predicted part.

4. The method according to claim 1, characterized in that: The method further comprises: Acquire the elevations of a plurality of fitting points in the elevation direction of the slag pile portion; Obtaining a temperature value corresponding to the elevation of the fitting point according to the temperature field distribution of the slag pile part; A temperature distribution fitting formula for the slag pile part is obtained according to the elevation of the fitting point and the temperature value corresponding to the elevation; wherein the temperature distribution fitting formula for the slag pile part uses the temperature value corresponding to the elevation as the horizontal coordinate and the elevation of the fitting point as the vertical coordinate.

5. The method according to claim 1, characterized in that The method further comprises: Obtain the temperature boundary, adiabatic boundary, fixed temperature boundary, thermal conductivity of the dam body, thermal conductivity of the water body, geothermal gradient of the dam foundation rock mass part along the elevation direction, thermal conductivity of the slag pile part, and thermal conductivity of the dam foundation rock mass part; wherein the temperature boundary includes the upper water temperature boundary of the slag pile part, the upper slag pile temperature boundary of the dam foundation rock mass part, the temperature boundary above the water surface of the upstream face of the dam body, the temperature boundary of the dam top, and the temperature boundary above the water surface of the downstream face of the dam body; the adiabatic boundary includes the upstream adiabatic temperature boundary and the downstream adiabatic temperature boundary; the fixed temperature boundary refers to the bottom geothermal temperature boundary; the thermal conductivity of the dam body is the thermal conductivity of the general rock mass; The temperature field of slag pile-dam foundation rock-dam body-water body is obtained according to the air temperature boundary, the adiabatic boundary, the fixed temperature boundary, the thermal conductivity of the slag pile part, the thermal conductivity of the dam foundation rock part, the thermal conductivity of the dam body, the thermal conductivity of the water body, the geothermal gradient and the heat conduction equation.

6. A temperature field calculation device, characterized in that: The device comprises: A first acquisition module is used to acquire a temperature boundary on one side of the predicted part; A model building module is used to build a model based on the temperature boundary on one side and the heat conduction equation Establishing a mathematical model of heat conduction of the predicted part; wherein a is the thermal conductivity of the predicted part; A second acquisition module is used to acquire the measured temperature of at least one characteristic point; wherein the characteristic point is located at the predicted position; A parameter inversion module, used to perform parameter inversion on the thermal conductivity of the predicted part according to the measured temperature of the characteristic point, so as to obtain the thermal conductivity of the predicted part; A third acquisition module is used to acquire the temperature field distribution of the predicted part according to the thermal conductivity of the predicted part and the heat conduction mathematical model of the predicted part; The predicted part is the slag dump part of the dam; the first acquisition module is specifically used to obtain the upper water temperature boundary of the slag dump part; wherein the average value of the stable measured values ​​of the thermometers at the same elevation above the slag dump part is selected as the upper water temperature boundary of the slag dump part; the model establishment module is specifically used to obtain the upper water temperature boundary of the slag dump part according to the upper water temperature boundary of the slag dump part and the heat conduction equation Establish a mathematical model of heat conduction of the slag pile part; wherein a1 is the thermal conductivity of the slag pile part; the second acquisition module is specifically used to obtain the measured temperature of at least one first characteristic point; wherein the first characteristic point is located in the slag pile part; the parameter inversion module is specifically used to perform parameter inversion on the thermal conductivity of the slag pile part according to the measured temperature of the first characteristic point to obtain the thermal conductivity of the slag pile part; the third acquisition module is specifically used to obtain the temperature field distribution of the slag pile part according to the thermal conductivity of the slag pile part and the mathematical model of heat conduction of the slag pile part; or, The predicted part is the rock mass of the dam foundation; the first acquisition module is specifically used to obtain the upper slag temperature boundary of the dam foundation rock mass; wherein, the temperature value is obtained by converting the attached temperature degree of the piezometer, and the temperature values ​​converted by multiple thermometers at the same height are averaged, and the average value is used as the upper slag temperature boundary of the dam foundation rock mass; the model establishment module is specifically used to obtain the upper slag temperature boundary of the dam foundation rock mass according to the upper slag temperature boundary of the dam foundation rock mass and the heat conduction equation A mathematical model of heat conduction of the dam foundation rock mass is established; wherein a2 is the thermal conductivity of the dam foundation rock mass; the second acquisition module is specifically used to obtain the measured temperature of at least one second characteristic point; wherein the second characteristic point is located in the dam foundation rock mass; the parameter inversion module is specifically used to perform parameter inversion on the thermal conductivity of the dam foundation rock mass according to the measured temperature of the second characteristic point to obtain the thermal conductivity of the dam foundation rock mass; the third acquisition module is specifically used to obtain the temperature field distribution of the dam foundation rock mass according to the thermal conductivity of the dam foundation rock mass and the mathematical model of heat conduction of the dam foundation rock mass.

7. An electronic device, characterized in that: include: A processor and a memory, wherein the memory stores machine-readable instructions executable by the processor, and when the machine-readable instructions are executed by the processor, the method according to any one of claims 1 to 5 is performed.

8. A storage medium, characterized in that: The storage medium stores a computer program, and when the computer program is executed by a processor, the method according to any one of claims 1 to 5 is executed.

Images