A surface crack temperature field determination method, system and electronic equipment

Abstract

The application provides a surface crack temperature field determination method and system and electronic equipment, relates to the ground crack monitoring technical field, constructs a crack surface three-dimensional geometric model; according to the crack wind speed, determines the crack heat convection coefficient; according to the rock-soil conduction data and the rock-soil physics data, determines the rock-soil heat conduction amount; according to the solar radiation data, determines the solar radiation intensity; according to the crack surface three-dimensional geometric model, the crack heat convection coefficient, the rock-soil heat conduction amount and the solar radiation intensity, constructs a surface crack temperature field calculation model; according to the surface crack temperature field calculation model, determines the surface crack temperature field data in the target mining area. The application can complete the monitoring of the surface crack temperature field data in the target mining area by constructing the surface crack temperature field calculation model, provides data support for the unmanned aerial vehicle infrared monitoring equipment, and provides a basis for the surface crack evolution.

Description

Technical Field

[0001] This invention relates to the field of ground fissure monitoring technology, and in particular to a method, system and electronic equipment for determining the temperature field of ground fissures. Background Technology

[0002] Surface fissures caused by coal mining have always been a major research focus in mine subsidence management. Efficient, intelligent, real-time, and accurate monitoring of surface fissure characteristics and morphological evolution is a prerequisite for managing mining subsidence and surface fissures. Current research on surface fissures in coal mine subsidence areas mainly focuses on their formation mechanisms and the damage they cause to the surface. Monitoring surface fissures is an effective means of studying them. Currently, researchers worldwide primarily focus on the shape, size, and orientation of surface fissures, but studies on their temperature fields are scarce.

[0003] For example, Chinese patent document CN113295843A discloses "A Monitoring Device for Ground Fissures Caused by Mining in Shallow Coal Seams." This device monitors the ground surface around the ground fissures in real time by installing ground fissure monitoring devices around the monitoring area and transmitting the data to a monitoring platform, thereby identifying the spatial orientation, diffusion points, and size of the ground fissures. Chinese patent document CN114155378A discloses "Automatic Extraction Method, Device, and Storage Medium for Ground Fissures in Coal Mining Subsidence." This method accurately obtains the distribution and morphological characteristics of ground fissures by deploying ground image control points and performing multiple measurements of the ground fissure monitoring area using UAV photography for correction. The monitoring devices in these two patents are relatively primitive, with low measurement accuracy and inability to study the temperature field of surface fissures. Changes in the temperature field of surface fissures are also a significant factor affecting coal mine subsidence. During the day, the temperature of the mine surface rises under sunlight. Fissures generated during coal mining connect the surface and underground mine tunnels, causing some of the ventilation air to leak to the surface through the fissures. The heat carried away by the leaking air makes the temperature around the fissures lower than the average temperature of the surface directly exposed to sunlight. As temperatures drop at night, the surface temperature decreases while heat from deeper soil is transferred to the surface through cracks, causing the crack temperature to be higher than the average surface temperature. Monitoring changes in the temperature field of surface cracks is a pressing issue that needs to be addressed. Summary of the Invention

[0004] The purpose of this invention is to provide a method, system, and electronic device for determining the surface fracture temperature field. By constructing a calculation model of the surface fracture temperature field, it is possible to monitor the surface fracture temperature field data in the target mining area, provide data support for UAV infrared monitoring equipment, and provide a basis for the evolution of surface fractures.

[0005] To achieve the above objectives, the present invention provides the following solution:

[0006] A method for determining the temperature field of surface fissures includes:

[0007] Acquire fracture data and solar radiation data within the target mining area; the fracture data includes surface fracture geometric parameters, fracture wind speed, soil-rock conduction data, and soil-rock physics data.

[0008] A three-dimensional geometric model of the fractured surface is constructed based on the geometric parameters of the surface fractures; the three-dimensional geometric model of the fractured surface is used to describe the surface soil and fracture geological structure within the target mining area;

[0009] The fracture thermal convection coefficient is determined based on the fracture wind velocity.

[0010] The thermal conductivity of the soil and rock is determined based on the soil and rock conductivity data and the soil and rock physics data.

[0011] Determine the solar radiation intensity based on the solar radiation data;

[0012] Based on the three-dimensional geometric model of the fractured surface, the thermal convection coefficient of the fracture, the thermal conductivity of the soil and rock, and the solar radiation intensity, a calculation model for the surface fracture temperature field is constructed.

[0013] Based on the surface fracture temperature field calculation model, the surface fracture temperature field data within the target mining area are determined.

[0014] Optionally, the surface fracture temperature field calculation model is constructed using Comsol with MATLAB software.

[0015] Optionally, the fissure heat convection coefficient is: h0 = 5.8 + 3.7v a ;

[0016] Where h0 is the fracture heat convection coefficient; v a The wind speed at the crack.

[0017] Optionally, the thermal conductivity of the soil and rock is:

[0018] Where ρ is the soil density; C g V is the volumetric heat capacity of the soil; k is the thermal conductivity of the soil; v a Where is the wind speed through the crack; T is the temperature; t is the time. For divergence; Q0 represents the temperature increment; Q0 represents the thermal conductivity of the soil and rock.

[0019] Optionally, the solar radiation intensity is: Q θ =I Dθ +I dθ ;

[0020] In the formula, Q θ Solar radiation intensity; IDθ I represents the intensity of direct solar radiation absorbed by the soil. Dθ =I0exp(-γl) (1 / sinβ) cosθ τ I0 is the normal surface solar radiation intensity at the upper boundary of the atmosphere, γ is the proportionality constant, l is the solar radiation distance from the upper boundary of the atmosphere to the ground when the sun is at the top of the sky, β is the solar altitude angle, and θ is the solar radiation intensity at the upper boundary of the atmosphere. τ The angle of incidence of the sun; I dθ I is the total scattered radiation intensity received by the soil. dθ =I ds +I dg I ds I represents the intensity of atmospheric scattering in the sky. ds =CI0exp(-γl) (1 / sinβ) (1-F g C is the ratio of the intensity of scattered radiation from the sky to the intensity of direct radiation reaching the ground, and F is the intensity of scattered radiation from the sky to the ground. g F is the angle coefficient between the ground and the horizontal plane. g =0.5(1-cosθ), where θ is the angle between the ground and the horizontal plane, I dg I represents the intensity of radiation reflected by the surface soil. dg =ρ g exp(-γl) (1 / sinβ) (C+sinβ)F g , ρ g This refers to the ground reflectivity.

[0021] A system for determining the temperature field of surface fractures, comprising:

[0022] The data acquisition module is used to acquire crack data and solar radiation data within the target mining area; the crack data includes surface crack geometric parameters, crack wind speed, soil-rock conduction data, and soil-rock physics data.

[0023] The fracture surface three-dimensional geometric model construction module is used to construct a fracture surface three-dimensional geometric model based on the surface fracture geometric parameters; the fracture surface three-dimensional geometric model is used to describe the surface soil and fracture geological structure within the target mining area;

[0024] A fracture thermal convection coefficient determination module is used to determine the fracture thermal convection coefficient based on the fracture wind velocity.

[0025] The soil and rock thermal conductivity determination module is used to determine the soil and rock thermal conductivity based on the soil and rock conductivity data and the soil and rock physics data.

[0026] A solar radiation intensity determination module is used to determine the solar radiation intensity based on the solar radiation data.

[0027] The surface fracture temperature field calculation model construction module is used to construct a surface fracture temperature field calculation model based on the three-dimensional geometric model of the fracture surface, the fracture thermal convection coefficient, the thermal conductivity of the soil and rock, and the solar radiation intensity.

[0028] The surface fracture temperature field data determination module is used to determine the surface fracture temperature field data in the target mining area based on the surface fracture temperature field calculation model.

[0029] An electronic device includes a memory and a processor, the memory storing a computer program, and the processor running the computer program to cause the electronic device to perform the method for determining the temperature field of a surface crack.

[0030] Optionally, the memory is a readable storage medium.

[0031] According to specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0032] This invention provides a method, system, and electronic device for determining the surface fracture temperature field. The method includes: acquiring fracture data and solar radiation data within a target mining area; the fracture data includes surface fracture geometric parameters, fracture wind speed, soil-rock conductivity data, and soil-rock physics data; constructing a three-dimensional geometric model of the fracture surface based on the surface fracture geometric parameters; the three-dimensional geometric model of the fracture surface is used to describe the surface soil and fracture geological structure within the target mining area; determining the fracture thermal convection coefficient based on the fracture wind speed; determining the soil-rock thermal conductivity based on the soil-rock conductivity data and the soil-rock physics data; determining the solar radiation intensity based on the solar radiation data; constructing a surface fracture temperature field calculation model based on the three-dimensional geometric model of the fracture surface, the fracture thermal convection coefficient, the soil-rock thermal conductivity, and the solar radiation intensity; and determining the surface fracture temperature field data within the target mining area based on the surface fracture temperature field calculation model. This invention, by constructing a surface fracture temperature field calculation model, can monitor the surface fracture temperature field data within a target mining area, providing data support for UAV infrared monitoring equipment and providing a basis for the evolution of surface fractures. Attached Figure Description

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

[0034] Figure 1 This is a flowchart of the method for determining the surface fracture temperature field in Embodiment 1 of the present invention;

[0035] Figure 2 This is a first principle diagram of the method for determining the surface fracture temperature field in Embodiment 1 of the present invention;

[0036] Figure 3 This is a second principle diagram of the method for determining the surface fracture temperature field in Embodiment 1 of the present invention;

[0037] Figure 4 This is a structural diagram of the model for calculating the surface fracture temperature field in Embodiment 1 of the present invention;

[0038] Figure 5 This is a schematic diagram of the ground fissure obtained in Embodiment 1 of the present invention;

[0039] Figure 6 This is a schematic diagram of the ground fissure features extracted in Embodiment 1 of the present invention;

[0040] Figure 7 This is a schematic diagram of the temperature field results of the ground fissures calculated in Embodiment 1 of the present invention. Detailed Implementation

[0041] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0042] The purpose of this invention is to provide a method, system, and electronic device for determining the surface fracture temperature field. By constructing a surface fracture temperature field calculation model, it is possible to monitor the surface fracture temperature field data in the target mining area, provide data support for UAV infrared monitoring equipment, and provide a basis for the evolution of surface fractures.

[0043] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0044] Example 1

[0045] like Figure 1 As shown in the figure, this embodiment provides a method for determining the temperature field of surface fractures, the method including:

[0046] Step 101: Obtain fracture data and solar radiation data within the target mining area; fracture data includes surface fracture geometric parameters, fracture wind speed, soil-rock conduction data, and soil-rock physics data;

[0047] Step 102: Construct a three-dimensional geometric model of the fractured surface based on the geometric parameters of the surface fractures; the three-dimensional geometric model of the fractured surface is used to describe the surface soil and fracture geological structure within the target mining area;

[0048] Step 103: Determine the fracture heat convection coefficient based on the fracture wind velocity; the fracture heat convection coefficient is: h0 = 5.8 + 3.7v a Where h0 is the fracture thermal convection coefficient; v a The wind speed at the crack.

[0049] Step 104: Determine the thermal conductivity of the soil and rock based on the soil conductivity data and soil physics data; the thermal conductivity of the soil and rock is: Where ρ is the soil density; C g V is the volumetric heat capacity of the soil; k is the thermal conductivity of the soil; v a Where is the wind speed through the crack; T is the temperature; t is the time. For divergence; Q0 represents the temperature increment; Q0 represents the thermal conductivity of the soil and rock.

[0050] Step 105: Determine the solar radiation intensity based on the solar radiation data; the solar radiation intensity is: Q θ =I Dθ +I dθ ;

[0051] In the formula, Q θ Solar radiation intensity; I Dθ I represents the intensity of direct solar radiation absorbed by the soil. Dθ =I0exp(-γl) (1 / sinβ) cosθ τ I0 is the normal surface solar radiation intensity at the upper boundary of the atmosphere, γ is the proportionality constant, l is the solar radiation distance from the upper boundary of the atmosphere to the ground when the sun is at the top of the sky, β is the solar altitude angle, and θ is the solar radiation intensity at the upper boundary of the atmosphere. τ The angle of incidence of the sun; I dθ I is the total scattered radiation intensity received by the soil. dθ =I ds +I dg I ds I represents the intensity of atmospheric scattering in the sky. ds =CI0exp(-γl) (1 / sinβ) (1-F g C is the ratio of the intensity of scattered radiation from the sky to the intensity of direct radiation reaching the ground, and F is the intensity of scattered radiation from the sky to the ground. g F is the angle coefficient between the ground and the horizontal plane. g =0.5(1-cosθ), where θ is the angle between the ground and the horizontal plane, I dg I represents the intensity of radiation reflected by the surface soil. dg =ρg exp(-γl) (1 / sinβ) (C+sinβ)F g , ρ g This refers to the ground reflectivity.

[0052] Step 106: Based on the three-dimensional geometric model of the fractured surface, the thermal convection coefficient of the fracture, the thermal conductivity of the soil and rock, and the solar radiation intensity, a calculation model of the surface fracture temperature field is constructed using Comsol with MATLAB software.

[0053] like Figures 2-3 Soil was collected uniformly from the shallow topsoil layer of the mining area, and the volumetric heat capacity C of the shallow topsoil layer was determined. g Laboratory tests were conducted on soil thermal conductivity (k), soil heat absorption and dissipation, and other thermal conductivity properties, as well as soil density (ρ), soil structure, and porosity, to obtain corresponding soil thermal and physical property parameters. Due to the presence of ground fissures, some of the underground ventilation airflow was transferred to the surface through these fissures; the wind speed (v) leaking through the fissures was measured using a wind speed monitoring device. a The ground fissure monitoring device was used to scan the morphology of surface fissures and measure characteristic parameters such as the length and width of the surface fissures; the solar altitude angle β and solar incidence angle θ were collected at any time at the monitoring location. τ The angle θ between the ground and the horizontal plane, and the ground reflectivity ρ g Wait for solar radiation data.

[0054] Based on the measured air velocity leaking through the crack, the thermal convection coefficient is calculated using the formula h0 = 5.8 + 3.7v. a The thermal convection coefficient of the crack was calculated; the solar radiation intensity was calculated based on the measured solar radiation data. The heat of solar radiation is usually expressed as radiation intensity, and the total solar radiation intensity Q at any solar incidence angle is calculated. θ =I Dθ +I dθ Substitute the test and calculation results into the soil heat conduction equation. The heat conduction equation of soil and rock was further optimized and derived, and a coupled mathematical model of surface crack temperature and air leakage was established with the calculation of soil and rock heat conduction, the calculation of the thermal convection coefficient of cracks and solar radiation intensity as input conditions.

[0055] Based on the measured surface crack morphology information, a three-dimensional geometric model of the cracked surface was established in AutoCAD to reconstruct the monitored surface soil and crack geological structure. Using the Comsol with MATLAB platform, a surface crack temperature field calculation module was developed based on the theory of partial differential equation solving and the established surface crack temperature-leakage coupled mathematical model. The established three-dimensional model was imported into Comsol with MATLAB to establish a numerical model of the soil geological structure. Using the calculated solar radiation intensity data, the surface-to-surface radiation module of Comsol with MATLAB was called to simulate the heat transferred to the surface by solar radiation, and the latitude and longitude, date and time, and solar radiation intensity of the measured crack locations were set, with the date and time accurate to the second. The solid heat transfer module was called to set the heat transfer conditions of other soil surfaces. Since the soil is in contact with other soils on all sides, the heat outflow and inflow conditions are the same; therefore, the soil boundaries were approximately considered thermally insulating. The bottom of the soil is affected by deep geothermal activity; a constant temperature was set at the bottom of the soil during model construction to simulate the impact of deep geothermal activity on the soil. When the Earth's surface is exposed to sunlight, diffuse reflection occurs. Therefore, a diffuse reflection condition is added to the soil surface. The heat transfer module is coupled with the surface-to-surface radiation module to construct a model of the material's actual heating condition.

[0056] The measured thermal conductivity, density, constant-pressure heat capacity, porosity, and other data of the surface soil in the mining area, as well as the constant-pressure heat capacity, density, and thermal conductivity of air in the soil pores, are assigned to the geological structure model material. The energy conservation module is then used for calculation. After the model is built and its components and surrounding environmental conditions are assigned values ​​according to actual conditions, the solver function in Comsol with MATLAB is used to calculate the temperature change of the geological structure model over a day. The structure of the model for calculating the surface fracture temperature field is as follows: Figure 4 As shown, real-time testing and quantification of solar irradiation conditions were performed. Using Comsol with MATLAB software, the temperature field distribution around the studied surface fissures was calculated at different times of day and under different light intensities. Based on the data obtained from the above model simulation, the relationship between surface fissure characteristic parameters, solar irradiance, and temperature distribution was statistically analyzed. Temperature field data are shown below. Figures 5-7 As shown.

[0057] Since the air volume and velocity of ventilation in underground mines remain constant during normal mine operation, the measured wind velocity v a The data is also constant and can be monitored and updated periodically. Solar irradiation conditions, however, change in real time. By testing these conditions in real time, solar radiation is further quantified, thus achieving accurate calculation of the real-time temperature field of the ground fissures.

[0058] Step 107: Determine the surface fracture temperature field data within the target mining area based on the surface fracture temperature field calculation model.

[0059] The technical solution provided by the present invention will be described in detail below.

[0060] Step 1: Collect soil samples evenly from the shallow topsoil layer of the mining area. Conduct laboratory tests on the thermal conductivity and physical properties of the shallow topsoil layer to determine the soil heat capacity C. g =2.0 J·cm -3 ·K -1 The thermal conductivity of the soil is k = 1.2 W·m. - 1 K -1 Soil thermal properties parameters such as soil heat absorption and dissipation, and soil density ρ = 2.6 g·cm³ -3 Physical parameters, etc.

[0061] Step 2: Due to the presence of ground fissures, some of the underground ventilation air is transmitted to the surface through these fissures. A wind speed monitoring device is used to measure the wind speed v leaking from the ground fissures. a =1.2m / s.

[0062] Step 3, Figure 4 It is a ground fissure image obtained using a ground fissure monitoring device. After obtaining the ground fissure image, the image is analyzed to obtain characteristic parameters such as the length and width of the ground fissure.

[0063] Step 4: Collect solar radiation data at a specific moment at the monitoring location, such as solar altitude angle β = 30° and solar incidence angle θ. τ =45°, the angle between the ground and the horizontal plane θ = 15°, and the ground reflectivity ρ g =20%, etc.

[0064] The solar radiation data in step 4 can provide real-time solar radiation parameters for the surface fracture temperature field calculation model through real-time monitoring and calculation, thereby realizing dynamic and intelligent monitoring of the surface fracture temperature field.

[0065] Step 5: Based on the wind velocity of the leaking air in the crack measured in Step 2, use the thermal convection coefficient calculation formula h0 = 5.8 + 3.7v a The thermal convection coefficient of the crack was calculated.

[0066] Step 6: Calculate the solar radiation intensity based on the solar radiation data measured in Step 4, using the formula Q. θ =I Dθ +I dθ The total solar radiation intensity Q at any solar incidence angle is calculated. θ =924W·m 2 In the formula, IDθ =I0exp(-γl) (1 / sinβ) cosθ τ ;I dθ =I ds +I dg In the formula, I ds =CI0exp(-γl) (1 / sinβ) (1-F g ), I dg =ρ g exp(-γl) (1 / sinβ) (C+sinβ)F g .

[0067] Step 7: Substitute the test and calculation results from steps 1-6 into the soil heat conduction equation. The heat conduction equation of soil and rock was further optimized and derived, and a coupled mathematical model of surface crack temperature and air leakage was established with the calculation of soil and rock heat conduction, the calculation of the thermal convection coefficient of cracks and solar radiation intensity as input conditions.

[0068] Step 8, as follows Figure 6 As shown, based on the surface crack morphology information measured in step 3, a two-dimensional geometric model of the cracked surface is established in AutoCAD to restore the monitored surface soil and crack geological structure.

[0069] Step 9: Using the Comsol with MATLAB platform, based on the theory of solving partial differential equations, develop a surface fracture temperature field calculation module using the surface fracture temperature-air leakage coupled mathematical model established in Step 7.

[0070] The two-dimensional models established in steps 10 and 9 are imported into Comsol with MATLAB software to create a numerical model of soil geological structure.

[0071] Step 11: Using the solar radiation intensity data calculated in Step 6, call the Comsol with MATLAB surface-to-surface radiation module to simulate the heat transferred to the Earth's surface by the sun's radiation, and set the latitude and longitude, date and time, and solar radiation intensity of the measured crack location, with the date and time set to the second.

[0072] Step 12: Use the solid heat transfer module to set the heat transfer conditions for other soil surfaces. The soil is in contact with other soil surfaces, and the heat outflow and inflow conditions are the same; therefore, the soil boundaries are approximated as thermal insulation. The bottom of the soil is affected by deep geothermal activity; a constant temperature is set at the bottom of the soil during model building to simulate the impact of deep geothermal activity on the soil. The surface experiences diffuse reflection when exposed to sunlight; therefore, a diffuse reflection condition is added to the soil surface.

[0073] Step 13: Couple the heat transfer module with the surface-to-surface radiation module to construct a model of the material's actual heating condition.

[0074] Step 14: Assign the data obtained in Step 1, such as the thermal conductivity, density, constant pressure heat capacity, porosity, etc. of the surface soil in the mining area, as well as the constant pressure heat capacity, density, thermal conductivity, etc. of the air in the soil pores, to the geological structure numerical model material established in Step 10.

[0075] Step 15: After the model is built and the model itself and its surrounding environment are assigned values ​​according to the actual conditions, the temperature distribution of the geological structure model is calculated using the solver function in Comsol with MATLAB.

[0076] Step 16, as follows Figure 7 As shown, the temperature field distribution around the surface fissures under study was obtained through calculations using Comsol with MATLAB software.

[0077] Example 2

[0078] In order to implement the method corresponding to Embodiment 1 above and achieve the corresponding functions and technical effects, a surface fracture temperature field determination system is provided below, including:

[0079] The data acquisition module is used to acquire crack data and solar radiation data within the target mining area; the crack data includes surface crack geometric parameters, crack wind speed, soil-rock conduction data, and soil-rock physics data.

[0080] The fracture surface 3D geometric model construction module is used to construct a fracture surface 3D geometric model based on the surface fracture geometric parameters; the fracture surface 3D geometric model is used to describe the surface soil and fracture geological structure within the target mining area.

[0081] The fracture thermal convection coefficient determination module is used to determine the fracture thermal convection coefficient based on the fracture wind velocity.

[0082] The module for determining the thermal conductivity of soil and rock is used to determine the thermal conductivity of soil and rock based on soil and rock conductivity data and soil and rock physics data.

[0083] The solar radiation intensity determination module is used to determine the solar radiation intensity based on solar radiation data.

[0084] The surface fracture temperature field calculation model construction module is used to construct a surface fracture temperature field calculation model based on the three-dimensional geometric model of the fracture surface, the fracture thermal convection coefficient, the thermal conductivity of the soil and rock, and the solar radiation intensity.

[0085] The surface fracture temperature field data determination module is used to determine the surface fracture temperature field data in the target mining area based on the surface fracture temperature field calculation model.

[0086] Example 3

[0087] This embodiment provides an electronic device, including a memory and a processor. The memory stores a computer program, and the processor runs the computer program to enable the electronic device to execute the method for determining the surface fracture temperature field described in Embodiment 1. The memory is a readable storage medium.

[0088] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the systems disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple; relevant parts can be referred to the method section.

[0089] This document uses specific examples to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. Furthermore, those skilled in the art will recognize that, based on the ideas of the present invention, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for determining the temperature field of surface fissures, characterized in that, include: Acquire fracture data and solar radiation data within the target mining area; the fracture data includes surface fracture geometric parameters, fracture wind speed, soil-rock conduction data, and soil-rock physics data. A three-dimensional geometric model of the fractured surface is constructed based on the geometric parameters of the surface fractures; the three-dimensional geometric model of the fractured surface is used to describe the surface soil and fracture geological structure within the target mining area; The fracture thermal convection coefficient is determined based on the fracture wind velocity. Based on the soil and rock conductivity data and the soil and rock physics data, the soil and rock thermal conductivity is determined; the soil and rock thermal conductivity is: ; Where ρ is the soil density; is the volumetric heat capacity of the soil; k is the thermal conductivity of the soil. Where is the wind speed through the crack; T is the temperature; t is the time. For divergence; This represents the temperature increment. This refers to the thermal conductivity of soil and rock. Based on the solar radiation data, the solar radiation intensity is determined; the solar radiation intensity is: ; In the formula, Solar radiation intensity; The intensity of direct solar radiation absorbed by the soil. , γ represents the intensity of direct solar radiation at the normal surface of the upper atmosphere boundary, where γ is a proportionality constant. This is the distance of solar radiation from the upper boundary of the atmosphere to the ground when the sun is at the top of the sky. The solar altitude angle, The angle of incidence of the sun; The total intensity of scattered radiation received by the soil. , The intensity of atmospheric scattering in the sky. C is the ratio of the intensity of scattered radiation from the sky to the intensity of direct radiation reaching the ground. This is the angle coefficient between the ground and the horizontal plane. , The angle between the ground and the horizontal plane. The intensity of radiation reflected by the surface soil; , Ground reflectivity; Based on the three-dimensional geometric model of the fractured surface, the thermal convection coefficient of the fracture, the thermal conductivity of the soil and rock, and the solar radiation intensity, a calculation model for the surface fracture temperature field is constructed. Based on the surface fracture temperature field calculation model, the surface fracture temperature field data within the target mining area are determined.

2. The method for determining the surface fracture temperature field according to claim 1, characterized in that, The surface fracture temperature field calculation model was constructed using Comsol with MATLAB software.

3. The method for determining the surface fracture temperature field according to claim 1, characterized in that, The fissure thermal convection coefficient is: ; in, The fissure heat convection coefficient; The wind speed at the crack.

4. A system for determining the temperature field of surface fissures, characterized in that, include: The data acquisition module is used to acquire crack data and solar radiation data within the target mining area; the crack data includes surface crack geometric parameters, crack wind speed, soil-rock conduction data, and soil-rock physics data. The fracture surface three-dimensional geometric model construction module is used to construct a fracture surface three-dimensional geometric model based on the surface fracture geometric parameters; the fracture surface three-dimensional geometric model is used to describe the surface soil and fracture geological structure within the target mining area; A fracture thermal convection coefficient determination module is used to determine the fracture thermal convection coefficient based on the fracture wind velocity. The soil and rock thermal conductivity determination module is used to determine the soil and rock thermal conductivity based on the soil and rock conductivity data and the soil and rock physics data; the soil and rock thermal conductivity is: ; Where ρ is the soil density; is the volumetric heat capacity of the soil; k is the thermal conductivity of the soil. Where is the wind speed through the crack; T is the temperature; t is the time. For divergence; This represents the temperature increment. This refers to the thermal conductivity of soil and rock. A solar radiation intensity determination module is used to determine the solar radiation intensity based on the solar radiation data; the solar radiation intensity is: ; In the formula, Solar radiation intensity; The intensity of direct solar radiation absorbed by the soil. , γ represents the intensity of direct solar radiation at the normal surface of the upper atmosphere boundary, where γ is a proportionality constant. This is the distance of solar radiation from the upper boundary of the atmosphere to the ground when the sun is at the top of the sky. The solar altitude angle, The angle of incidence of the sun; The total intensity of scattered radiation received by the soil. , The intensity of atmospheric scattering in the sky. C is the ratio of the intensity of scattered radiation from the sky to the intensity of direct radiation reaching the ground. This is the angle coefficient between the ground and the horizontal plane. , The angle between the ground and the horizontal plane. The intensity of radiation reflected by the surface soil; , Ground reflectivity; The surface fracture temperature field calculation model construction module is used to construct a surface fracture temperature field calculation model based on the three-dimensional geometric model of the fracture surface, the fracture thermal convection coefficient, the thermal conductivity of the soil and rock, and the solar radiation intensity. The surface fracture temperature field data determination module is used to determine the surface fracture temperature field data in the target mining area based on the surface fracture temperature field calculation model.

5. An electronic device, characterized in that, The device includes a memory and a processor, the memory being used to store a computer program, and the processor running the computer program to cause the electronic device to perform a method for determining the surface fracture temperature field according to any one of claims 1 to 3.

6. An electronic device according to claim 5, characterized in that, The memory is a readable storage medium.

Images