Stone glacier ice content inversion method
By extracting the dynamic model parameters of stone glaciers and combining with the empirical rheology model of stone glaciers, the high cost and low accuracy problems of estimating ice content in stone glaciers are solved, and a higher precision inversion effect is achieved.
Patent Information
- Application Number
- CN202510532393.9
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-04-25
- Publication Date
- 2025-08-08
AI Technical Summary
The prior art has problems of high cost and low accuracy when estimating the ice content of stone glaciers, especially in high altitude areas, which are difficult to apply on a large scale. The observation error and area estimation uncertainty of traditional dynamic models lead to insufficient inversion accuracy.
By extracting the dynamic model parameters of stone glaciers and the average annual deformation rate, combined with the empirical rheology model of stone glaciers, the ice content of stone glaciers was obtained, which overcomes the DInSAR observation error and the uncertainty of traditional dynamic models, and improves the inversion accuracy.
A higher precision ice content inversion of stone glaciers is achieved, reducing the impact of topographic errors and phase disconnection, and improving the accuracy of rheological model parameters.
Smart Images
Figure CN120449456A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of data processing, in particular to a method for inverting ice content in a rock glacier. Background Art
[0002] There are currently three main methods for quantitatively estimating the ice content within rock glaciers: geophysical exploration, empirical formula estimation, and dynamic model inversion. The characteristics of these ice content estimation methods are analyzed as follows:
[0003] Geophysical exploration mainly relies on ground-penetrating radar, electrical resistance tomography, and active source seismic inversion to detect the internal structure and stratification of rock glaciers, and to infer the distribution, thickness, and volume of frozen ice. Although this type of method is relatively accurate in estimating the ice content of rock glaciers using physical exploration, it has the problems of high cost and cannot be applied on a large scale in dangerous and high-altitude areas.
[0004] The ice content of rock glaciers is estimated based on an empirical formula. That is, the thickness of the rock glacier is calculated based on the empirical relationship between its area and thickness. Then, the ice content of the rock glacier is estimated by assuming that the ice content of the rock glacier accounts for a certain proportion of the overall volume (generally 40% to 60%). Although the empirical formula method can efficiently estimate the ice content of a large area of rock glaciers, the estimation of volumetric ice content lacks a physical basis, and the water storage estimated by this method may have large deviations.
[0005] Traditional dynamic modeling methods for inverting rock glacier ice content simplify rock glaciers and, based on physical and mathematical models, input a series of parameters to invert the ice content of rock glaciers. This method can determine the internal ice content of rock glaciers at high altitudes using only remote sensing. Compared with traditional geodetic methods, its high efficiency and large-scale applicability offer significant advantages. However, traditional dynamic models often rely on two-view SAR images to determine rock glacier surface velocity, which often results in significant observational errors. Furthermore, inferring rock glacier thickness based on its area is subject to significant uncertainty, resulting in a need for improved accuracy in inverted ice content. Summary of the Invention
[0006] The purpose of the present invention is to overcome the shortcomings of the prior art and provide a method for inverting the ice content of a rock glacier.
[0007] The object of the present invention is achieved through the following technical solution: A method for inverting the ice content of a rock glacier comprises the following steps:
[0008] S1: Extract rock glacier dynamic model parameters through DEM files;
[0009] S2: Obtain the average annual deformation rate of rock glaciers;
[0010] S3: Import the parameters obtained in step S1 and step S2 into the rock glacier empirical rheological model, and invert to obtain the rock glacier ice content.
[0011] Preferably, step S1 further includes the following steps:
[0012] S11: Extract the area and width of rock glaciers through the rock glacier inventory results;
[0013] S12: Obtain the thickness data of the active layer of the rock glacier, extract the slope angle and altitude of the rock glacier through the DEM file, and calculate the thickness T, air density and permafrost thickness.
[0014]
[0015] Among them, θ slp is the slope angle of the rock glacier, h al is the thickness of the active layer, h core is the thickness of the permafrost layer;
[0016] S13: The air content of the active layer is obtained based on the debris content of the active layer, and the density of the active layer of the rock glacier is calculated by the debris density and the air density.
[0017] ρ al =f d,al ρ d +f a,al ρ a ;
[0018] Among them, ρ d is the debris density, ρ a is the air density, ρ al is the density of the active layer of rock glacier, f d,al is the debris content of the active layer, f a,al is the air content of the active layer;
[0019] S14: The surface factor S of the rock glacier is calculated by the width and thickness parameters of the rock glacier. f ,
[0020]
[0021] Where W is the width of the rock glacier and T is the thickness of the rock glacier.
[0022] Preferably, step S2 further includes the following steps:
[0023] S21: Obtain the LOS deformation rate of rock glaciers by time series analysis;
[0024] S22: Convert the LOS deformation rate into the average annual slope deformation rate, and extract the corresponding rate value for each rock glacier.
[0025]
[0026] Where α is the orientation angle of the rock glacier, is the azimuth of the satellite's flight direction, θ inc is the radar wave incident angle, θ slp is the slope angle of the rock glacier, V LOS is the satellite line-of-sight velocity, V slp is the satellite slope velocity.
[0027] Preferably, step S3 further includes the following steps:
[0028] S31: Convert the debris content of the permafrost layer into a form related to ice content and obtain the equivalent relationship between the density of the permafrost layer and the ice content.
[0029] f w,core =0.4735f d,core +0.0029;
[0030] f d,core =1-f a,core -f w,core -f i ;
[0031] ρ core =f d,core ρ d +f a,core ρ a +f w,core ρ w +f i ρ i ;
[0032] Among them, f w,core is the water content of the permafrost layer, f d,core is the debris content of the permafrost layer, f a,core is the air content of the permafrost layer, ρ core is the density of frozen soil, ρ d is the debris density, ρ a is the air density, ρ w is the water density, f i is the ice content of the rock glacier, ρ i is the density of ice;
[0033] S32: According to the effective viscosity B of the rock glacier and the ice content f i The proportional relationship between them is obtained to obtain the equal relationship between the two.
[0034]
[0035] S33: Based on the proportional relationship between the rock glacier flow law index n and the ice content, the equivalent relationship between the two is obtained.
[0036] n=3f i ;
[0037] S34: Input the parameters and equivalent relationships obtained in steps S31 to S33 into the rheological model, and invert to obtain the ice content of the rock glacier.
[0038]
[0039] Among them, u s is the surface velocity of the rock glacier, and g is the acceleration due to gravity.
[0040] The present invention has the following advantages: the present invention extracts the parameters of the rock glacier dynamic model and obtains the average annual deformation rate of the rock glacier, and then imports the rock glacier empirical rheological model based on the parameters obtained in steps S1 and S2 to invert the ice content of the rock glacier. This method not only overcomes the observation errors caused by terrain errors, atmospheric delays, and phase decoherence during DInSAR observations, but also compensates for the uncertainty problem in the traditional dynamic model when estimating the thickness of the rock glacier based on its area, thereby obtaining more accurate rheological model parameters and further improving the accuracy of the inverted ice content. BRIEF DESCRIPTION OF THE DRAWINGS
[0041] Figure 1 This is a schematic diagram of the process of the inversion method for the ice content of the rock glacier;
[0042] Figure 2 Schematic diagram of the model sensitivity results. DETAILED DESCRIPTION
[0043] To make the objectives, technical solutions, and advantages of the embodiments of the present invention more clear, the technical solutions of the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Generally, the components of the embodiments of the present invention described and shown in the drawings herein can be arranged and designed in various different configurations.
[0044] Therefore, the following detailed description of the embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention as claimed, but rather merely represents selected embodiments of the present invention. All other embodiments derived by persons of ordinary skill in the art based on the embodiments of the present invention without creative effort are also within the scope of protection of the present invention.
[0045] It should be noted that, in the absence of conflict, the embodiments and features of the embodiments of the present invention can be combined with each other.
[0046] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, it does not need to be further defined or explained in subsequent drawings.
[0047] In the description of the present invention, it should be noted that the terms "center," "upper," "lower," "left," "right," "vertical," "horizontal," "inner," "outer," etc., indicating orientations or positional relationships, are based on the orientations or positional relationships shown in the accompanying drawings, or are the orientations or positional relationships in which the inventive product is typically placed when in use, or are the orientations or positional relationships commonly understood by those skilled in the art. These terms are intended only to facilitate the description of the present invention and to simplify the description, and are not intended to indicate or imply that the device or element referred to must have a specific orientation, be constructed, or operate in a specific orientation. Therefore, they should not be construed as limiting the present invention. Furthermore, the terms "first," "second," etc., are used only to distinguish descriptions and should not be construed as indicating or implying relative importance.
[0048] In the description of the present invention, it should also be noted that, unless otherwise expressly specified or limited, the terms "disposed," "installed," "connected," and "connected" should be understood in a broad sense. For example, they may refer to fixed connections, detachable connections, or integral connections; they may refer to mechanical connections or electrical connections; they may refer to direct connections or indirect connections through an intermediate medium; and they may refer to internal communication between two components. Those skilled in the art will understand the specific meanings of the above terms in the present invention based on the specific circumstances.
[0049] In this embodiment, if Figure 1 As shown, a method for inverting the ice content of a rock glacier includes the following steps:
[0050] S1: Extract rock glacier dynamic model parameters through DEM files;
[0051] S2: Obtain the average annual deformation rate of rock glaciers;
[0052] S3: Importing the parameters obtained in steps S1 and S2 into the empirical rheological model of the rock glacier to invert the rock glacier ice content. By extracting rock glacier dynamic model parameters and obtaining the average annual deformation rate of the rock glacier, the parameters obtained in steps S1 and S2 are then imported into the empirical rheological model of the rock glacier to invert the rock glacier ice content. This method not only overcomes observational errors caused by terrain errors, atmospheric delay, and phase decoherence during DInSAR observations, but also compensates for the uncertainty in traditional dynamic models when estimating rock glacier thickness based on its area. This method results in more accurate rheological model parameters and, in turn, improves the accuracy of the inverted ice content.
[0053] Furthermore, step S1 further includes the following steps:
[0054] S11: Extract the area and width of rock glaciers through the rock glacier inventory results;
[0055] S12: Obtain the thickness data of the active layer of the rock glacier, extract the slope angle and altitude of the rock glacier through the DEM file, and calculate the thickness T, air density and permafrost thickness.
[0056]
[0057] Among them, θ slp is the slope angle of the rock glacier, h al is the thickness of the active layer, h core is the thickness of the permafrost layer;
[0058] S13: The air content of the active layer is obtained based on the debris content of the active layer, and the density of the active layer of the rock glacier is calculated by the debris density and the air density.
[0059] ρ al =f d,al ρ d +f a,al ρ a ;
[0060] Among them, ρ d is the debris density, ρ a is the air density, ρ al is the density of the active layer of rock glacier, f d,al is the debris content of the active layer, f a,al is the air content of the active layer;
[0061] S14: The surface factor S of the rock glacier is calculated by the width and thickness parameters of the rock glacier. f ,
[0062]
[0063] Where W is the width of the rock glacier, and t is the thickness of the rock glacier. Specifically, in step S12, the rock glacier active layer thickness data is obtained from the official website of the European Space Agency's Climate Change Initiative. In step S13, the active layer debris content is obtained from existing public documents. The active layer debris content here is 65%, which infers that the active layer air content is 35%.
[0064] Furthermore, step S2 further includes the following steps:
[0065] S21: Obtain the LOS deformation rate of rock glaciers by time series analysis;
[0066] S22: Convert the LOS deformation rate into the average annual slope deformation rate, and extract the corresponding rate value for each rock glacier.
[0067]
[0068] Where α is the orientation angle of the rock glacier, is the azimuth of the satellite's flight direction, θ inc is the radar wave incident angle, θ slp is the slope angle of the rock glacier, V LOS is the satellite line-of-sight velocity, V slp is the satellite slope velocity.
[0069] In this embodiment, step S3 further includes the following steps:
[0070] S31: Convert the debris content of the permafrost layer into a form related to ice content and obtain the equivalent relationship between the density of the permafrost layer and the ice content.
[0071] f w,core =0.4735f d,core +0.0029;
[0072] f d,core =1-f a,core -f w,core -f i ;
[0073] ρ core =f d,core ρ d +f a,core ρ a +f w,core ρ w +f i ρ i ;
[0074] Among them, f w,core is the water content of the permafrost layer, f d,core is the debris content of the permafrost layer, f a,core is the air content of the permafrost layer, ρ core is the density of frozen soil, ρ d is the debris density, ρ a is the air density, ρ w is the water density, f i is the ice content of the rock glacier, ρ i is the density of ice;
[0075] S32: According to the effective viscosity B of the rock glacier and the ice content f i The proportional relationship between them is obtained to obtain the equal relationship between the two.
[0076]
[0077] S33: Based on the proportional relationship between the rock glacier flow law index n and the ice content, the equivalent relationship between the two is obtained.
[0078] n=3f i ;
[0079] S34: Input the parameters and equivalent relationships obtained in steps S31 to S33 into the rheological model, and invert to obtain the ice content of the rock glacier.
[0080]
[0081] Among them, u s is the surface velocity of the rock glacier, and g is the acceleration due to gravity. Specifically, the proportional relationships among air, debris, and water content in the permafrost layer, the proportional relationship between the effective viscosity of the rock glacier and its ice content, and the proportional relationship between the flow regularity index of the rock glacier and its ice content were all obtained from existing public documents. The air content of the permafrost layer here is 7.5%.
[0082] In this embodiment, in order to verify the effect of the method of the present invention on the inversion method of the ice content of the rock glacier, the control variable method was used to explore the relationship between the rock glacier area A and the slope angle θ slp The change of ice content f i and ice volume content V i The influence of the calculation results is shown in Figure 2. Taking a randomly selected rock glacier as an example, the initial values of the area and slope angle are set to 0.514 km 2 and 20.58°, combined into Scn-1, and then generated parameter combinations of 0.4, 0.6, 0.8, 1.2, 1.4 and 1.6 times the initial value according to the 20% change gradient. By observing the f calculated by different parameter combinations i With V i changes, thus obtaining the sensitivity of both to changes in area and slope parameters. The experimental results are as follows Figure 2 shown.
[0083] Furthermore, the results of estimating ice content and ice volume content based on the traditional area dynamic model are shown in Table 1.
[0084] Table 1
[0085] Scheme <![CDATA[A(km 2 )]]> T(m) <![CDATA[h core (m)]]> <![CDATA[f i_area (%)]]> <![CDATA[V i_area (10 -3 km 3 )]]> Scn-0.4 0.205746 36.444860 35.744860 75.55 5.556221 Scn-0.6 0.308619 39.523422 38.823422 76.1 9.118033 Scn-0.8 0.411492 41.864151 41.164151 76.4 12.941182 Scn-1 0.514365 43.774812 43.074812 76.45 16.938397 Scn-1.2 0.617238 45.400490 44.700490 76.5 21.106994 Scn-1.4 0.720111 46.821993 46.121993 76.4 25.374698 Scn-1.6 0.822984 48.089281 47.389281 76.4 29.796475
[0086] The results of the ice content and ice volume content estimation method disclosed in the present invention are shown in Table 2.
[0087] Table 2
[0088] Scheme <![CDATA[θ slp (°)]]> <![CDATA[A(km 2 )]]> T(m) <![CDATA[h core (m)]]> <![CDATA[f i_slp (%)]]> <![CDATA[V i_slp (10 -3 km 3 )]]> Scn-0.4 8.232 0.514365 29.5912 28.891200 67.5 10.030920 Scn-0.6 12.348 0.514365 25.8868 25.186800 69.35 8.984437 Scn-0.8 16.464 0.514365 22.1824 21.482400 69.9 7.723807 Scn-1 20.58 0.514365 18.478 17.778000 69.6 6.364489 Scn-1.2 24.696 0.514365 14.7736 14.073600 68.3 4.944215 Scn-1.4 28.812 0.514365 11.0692 10.369200 64.55 3.442809 Scn-1.6 32.928 0.514365 7.3648 6.664800 62.7 2.149444
[0089] The comparison of ice content and ice volume content parameters of the two schemes is shown in Table 3.
[0090] Table 3
[0091] Scheme <![CDATA[f i_area (%)]]> <![CDATA[f i_slp (%)]]> <![CDATA[V i_area (10 -3 km 3 )]]> <![CDATA[V i_slp (10 -3 km 3 )]]> Scn-0.4 75.55 67.5 5.556221 10.030920 Scn-0.6 76.1 69.35 9.118033 8.984437 Scn-0.8 76.4 69.9 12.941182 7.723807 Scn-1 76.45 69.6 16.938397 6.364489 Scn-1.2 76.5 68.3 21.106994 4.944215 Scn-1.4 76.4 64.55 25.374698 3.442809 Scn-1.6 76.4 62.7 29.796475 2.149444
[0092] As shown in Tables 1 to 3, for different parameter combinations, the ice content obtained by the traditional area-based dynamic model inversion does not change much. The difference in ice content between Scn-0.4 and Scn-1.6 is 1.9%; the difference in ice volume content is 0.024 km 3 The ice content inverted by the method disclosed in the present invention is slightly smaller than that of the area model, in which the ice content difference between Scn-0.4 and Scn-1.6 is 4.8%, and the ice content obtained by the model in the two combinations of Scn-0.4 and Scn-1.6 differs from the distribution of Scn-1.0 by 2.1% and 4.6%. In addition, the ice volume content difference between Scn-0.4 and Scn-1.6 is 0.008 km 3 .
[0093] Another example Figure 2 The ice volume obtained by the two methods shows opposite trends. The ice volume obtained by the area model gradually increases between Scn-0.4 and Scn-1.6, that is, the difference between the minimum and maximum estimated values is 0.024 km 3 , while the results calculated by the method disclosed in the present invention gradually decrease, and the difference between the minimum and maximum estimated values is 0.008km 3 Therefore, compared with traditional dynamic models, this method not only overcomes the observation errors caused by terrain errors, atmospheric delay, and phase incoherence during DInSAR observations, but also makes up for the uncertainty problem in traditional dynamic models when estimating the thickness of rock glaciers based on their area, thereby obtaining more accurate rheological model parameters and improving the accuracy of the inverted ice content.
[0094] Although the present invention has been described in detail with reference to the aforementioned embodiments, it is still possible for those skilled in the art to modify the technical solutions described in the aforementioned embodiments, or to make equivalent substitutions for some of the technical features therein. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.
Claims
1. A method for inverting ice content in rock glaciers, characterized by: The following steps are involved: S1: Extract rock glacier dynamic model parameters through DEM files; S2: Obtain the average annual deformation rate of rock glaciers; S3: Import the parameters obtained in step S1 and step S2 into the rock glacier empirical rheological model, and invert to obtain the rock glacier ice content.
2. The method for inverting ice content of a rock glacier according to claim 1, characterized in that: The step S1 further includes the following steps: S11: Extract the area and width of rock glaciers through the rock glacier inventory results; S12: Obtain the thickness data of the active layer of the rock glacier, extract the slope angle and altitude of the rock glacier through the DEM file, and calculate the thickness T, air density and permafrost thickness. Among them, θ slp is the slope angle of the rock glacier, h al is the thickness of the active layer, h core is the thickness of the permafrost layer; S13: The air content of the active layer is obtained based on the debris content of the active layer, and the density of the active layer of the rock glacier is calculated by the debris density and the air density. r al =f d,al r d +f a,al r d ; Among them, ρ d is the debris density, ρ a is the air density, ρ al is the density of the active layer of rock glacier, f d,al is the debris content of the active layer, f a,al is the air content of the active layer; S14: The surface factor S of the rock glacier is calculated by the width and thickness parameters of the rock glacier. f , Where W is the width of the rock glacier and T is the thickness of the rock glacier.
3. The method for inverting ice content of a rock glacier according to claim 2, characterized in that: The step S2 further includes the following steps: S21: Obtain the LOS deformation rate of rock glaciers by time series analysis; S22: Convert the LOS deformation rate into the average annual slope deformation rate, and extract the rate value corresponding to each rock glacier. Where α is the orientation angle of the rock glacier, is the azimuth of the satellite's flight direction, θ inc is the radar wave incident angle, θ slp is the slope angle of the rock glacier, V LOS is the satellite line-of-sight velocity, V slp is the satellite slope velocity.
4. The method for inverting ice content of a rock glacier according to claim 3, characterized in that: The step S3 further includes the following steps: S31: Convert the debris content of the permafrost layer into a form related to ice content and obtain the equivalent relationship between the density of the permafrost layer and the ice content. <h2 style=";text-align:left;direction:ltr">f<h2 style=";text-align:left;direction:ltr"> w,core <h2 style=";text-align:left;direction:ltr"> =0.4735f<h2 style=";text-align:left;direction:ltr"> d,core <h2 style=";text-align:left;direction:ltr"> +0.0029; f d,core =1-f a,core -f w,core -f i ; r core =f d,core r d +f a,core r a +f w,core r w +f i r i ; Among them, f w,core is the water content of the permafrost layer, f d,core is the debris content of the permafrost layer, f a,core is the air content of the permafrost layer, ρ core is the density of frozen soil, ρ d is the debris density, ρ a is the air density, ρ w is the water density, f i is the ice content of the rock glacier, ρ i is the density of ice; S32: According to the effective viscosity B of the rock glacier and the ice content f i The proportional relationship between them is obtained to obtain the equal relationship between the two. S33: Based on the proportional relationship between the rock glacier flow law index n and the ice content, the equivalent relationship between the two is obtained. <h2 style=";text-align:left;direction:ltr">n=3f<h2 style=";text-align:left;direction:ltr"> i <h2 style=";text-align:left;direction:ltr"> ; S34: Input the parameters and equivalent relationships obtained in steps S31 to S33 into the rheological model, and invert to obtain the ice content of the rock glacier. Among them, u s is the surface velocity of the rock glacier, and g is the acceleration due to gravity.