Method for determining reclamation elevation of coal mining collapse cultivated land in mining area

Abstract

The invention discloses a mining area coal mining collapse cultivated land reclamation elevation determination method, and relates to the technical field of mining area coal mining collapse cultivated land reclamation, and the method comprises the steps: obtaining a mining area collapse region, a coal mining plan, an underground coal occurrence condition, a cultivated land damage boundary and a land use type, and carrying out the preprocessing; according to the pre-processed mining area collapse area, obtaining original surface elevation data before mining of the mining area through satellite image analysis; based on the preprocessed coal mining plan and underground coal occurrence conditions, predicting surface subsidence amount data at different stages by adopting a probability integral method, and superposing the predicted subsidence amount data and the original terrain elevation data to generate subsidence surface elevation data; through satellite image analysis and a probability integral method, the accurate reclamation elevation determination method for the mining area coal mining collapsed cultivated land is provided, firstly, original surface elevation data before mining of the mining area is obtained through satellite images, and the data obtaining precision and efficiency are remarkably improved;

Description

Technical Field

[0001] The invention relates to the technical field of reclamation of farmland subsided by coal mining in mining areas, in particular to a method for determining the elevation of reclamation of farmland subsided by coal mining in mining areas. Background Art

[0002] As global energy demand continues to grow, coal, as a traditional energy source, still occupies an important position in the energy structure of many countries. With the advancement of coal mining technology, open-pit and underground coal mining have become the mainstream, but coal mining activities often lead to land subsidence, changes in groundwater levels and environmental damage in mining areas, seriously affecting the quality of farmland and land use. Coal mining subsidence is caused by the collapse of the strata due to underground coal mining, resulting in surface subsidence and affecting the surrounding ecological environment. Therefore, how to effectively reclaim coal mining subsidence areas has become a key environmental issue in mining area development.

[0003] In recent years, land reclamation technology in mining areas has been continuously developing. The traditional reclamation method that relies on manual measurement and empirical design is inefficient and prone to errors, making it difficult to ensure the quality of reclamation. With the advancement of remote sensing technology, lidar and drone technology, the use of high-precision satellite images to obtain surface elevation data has become an important tool for reclamation design. These technologies make it possible to accurately assess the subsidence area, thereby providing a reliable basis for subsequent earthwork calculations and reclamation projects. At the same time, the application of modern computer modeling and numerical simulation technology has made the dynamic prediction of coal mining subsidence a reality, thereby improving the scientificity and accuracy of reclamation projects and greatly improving reclamation efficiency and construction quality. Summary of the invention

[0004] In view of the above existing problems, the present invention is proposed.

[0005] Therefore, the present invention provides a method for determining the reclamation elevation of cultivated land subsided by coal mining in mining areas to solve the problems in the prior art of being unable to accurately and dynamically adjust the reclamation elevation and ignoring the poor reclamation effect caused by changes in the subsidence process.

[0006] In order to solve the above technical problems, the present invention provides the following technical solutions: In a first aspect, the present invention provides a method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area, which comprises: Obtain mining subsidence areas, coal mining plans, underground coal occurrence conditions, cultivated land damage boundaries and land use types, and perform preprocessing; According to the pre-processed mining subsidence area, the original surface elevation data of the mining area before mining is obtained through satellite image analysis; Based on the pre-processed coal mining plan and underground coal occurrence conditions, the probability integral method is used to predict the surface subsidence data at different stages, and the predicted subsidence data is superimposed with the original terrain elevation data to generate the subsidence surface elevation data; According to the subsidence surface elevation data and the pre-processed cultivated land damage boundary, the reclamation area is divided into several terraced sub-areas; According to the terrace sub-areas, calculate the earthwork volume of the terrace sub-areas and generate earthwork balance data; According to the earthwork balance data and the pre-processed land use type, the reclamation construction elevation is dynamically adjusted; Carry out the reclamation project according to the adjusted reclamation construction elevation.

[0007] As a preferred solution of the method for determining the reclamation elevation of cultivated land in coal mining subsidence in mining areas of the present invention, the method comprises: obtaining the mining area subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundary and land use type, and performing preprocessing; the specific steps are: Obtain mining subsidence areas and farmland damage boundaries through drone aerial surveys; Obtain coal mining plans through the mine management department and historical data; Obtain underground coal occurrence conditions through geological exploration reports and drilling data provided by exploration units; Obtain land use types through land use planning maps provided by government departments and field surveys; Use GIS software to coordinate and time standardize the mining area subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundaries and land use types; Interpolation methods are used to fill in missing data on the mining subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundaries and land use types after coordinate standardization and time standardization. The mining subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundary and land use type after filling the data will be normalized.

[0008] As a preferred solution of the method for determining the elevation of cultivated land reclamation after coal mining in a mining area of ​​the present invention, the original surface elevation data before mining in the mining area is obtained by satellite image analysis based on the pre-processed mining area subsidence area, and the specific steps are: According to the pre-processed mining subsidence area, the WorldView-3 satellite was selected to obtain two images of different perspectives before mining, which were used as the left image and the right image. Performing geometric correction processing on the acquired left image and right image; According to the geometrically corrected left and right images, the focal length of the satellite sensor and the baseline distance during image acquisition are obtained; Use the left and right images after geometric correction to perform stereo image matching and calculate the disparity data of each pixel ; Based on the parallax data, the focal length of the satellite sensor and the baseline distance when the image was collected, the surface elevation of each location is calculated to generate a digital elevation model DEM, which is expressed as: ; in, Indicated in Surface elevation data for the location, Indicates the baseline distance when the satellite image was taken. represents the focal length of the satellite sensor, represents the horizontal coordinate in the space coordinate, Represents the ordinate in the spatial coordinate; The median filter method is used to denoise and filter the generated digital elevation model DEM; For the digital elevation model DEM after denoising and filtering, the missing data is filled by interpolation method; The digital elevation model DEM after filling in the missing data is exported to a standard geographic information format as the original surface elevation data before mining.

[0009] As a preferred solution of the method for determining the elevation of cultivated land reclamation in mining areas subsided by coal mining in the present invention, the method adopts the probability integral method to predict the surface subsidence data at different stages based on the pre-processed coal mining plan and underground coal occurrence conditions, and superimposes the predicted subsidence data with the original terrain elevation data to generate the subsidence surface elevation data. The specific steps are as follows: Based on the pre-processed coal mining plan and underground coal occurrence conditions, the probability integral method is used to predict the surface subsidence data at different stages. The expression is: ; in, Indicates at time Time, Location Surface subsidence data at It indicates the conversion ratio of subsidence caused by coal mining. Indicates the coal mining speed, Indicates the thickness of the coal seam, Indicates the depth of coal mining, Indicates the coal mining area, Indicates the depth of the coal seam, represents the spatial distribution function of subsidence, represents the sedimentation diffusion coefficient, Mathematical notation for exponential functions, represents the base of natural logarithms, represents the sinking time function, represents the sinking rate coefficient, Indicates the time of mining. represents the geological influence factor function, represents the cosine function, Indicates the coal seam inclination, Indicates the inclination of coal seam The cosine value of represents the geological influence factor; Superimpose the predicted subsidence data with the original surface elevation data to generate subsidence surface elevation data .

[0010] As a preferred solution of the method for determining the elevation of reclamation of subsided farmland in a mining area according to the present invention, the reclamation area is divided into a number of terraced sub-areas according to the subsided surface elevation data and the pre-processed farmland damage boundary, and the specific steps are as follows: Obtaining the damage boundary of cultivated land through drone aerial survey ; The subsidence surface elevation data was rasterized by GIS software, and the rasterized subsidence surface elevation data and the boundary of the damaged cultivated land area were spatially superimposed to obtain the rasterized subsidence surface elevation data in the reclamation area. ; Based on the rasterized subsidence surface elevation data in the reclamation area, the slope value of each grid cell is calculated; ; in, Indicates the time Time, Location The slope value at represents the symbol of partial derivative, and Respectively expressed in The elevation at the location is and rate of change of direction; Set the slope threshold range to ; in, Indicates the lower limit of the slope range, Indicates the upper limit of the slope range; Through the reclassification function of GIS software, the slope value of each grid cell in the reclamation area is divided into Reclassify the grid cells, and the classified grid cells will be marked with the terraced sub-area to which each grid cell belongs; when < When ; when ≤ ≤ When ; when > When ; Set the side length of the grid cell to , then the area of ​​the grid cell is .

[0011] As a preferred solution of the method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area according to the present invention, the specific steps of calculating the earthwork volume of the terraced sub-area and the earthwork balance data are as follows: Set the reclamation construction elevation to ; For each grid unit, when the gridded subsidence surface elevation data in the reclamation area is greater than the reclamation construction elevation, the grid unit belongs to the height-cutting area, and the height-cutting earthwork volume is calculated. ; For each grid unit, when the gridded subsidence surface elevation data in the reclamation area is lower than the reclamation construction elevation, the grid unit belongs to the padding area, and the padding earthwork volume is calculated. ; The earthwork balance data is generated by accumulating the earthwork volume of all grid cells, and the expression is: ; in, Indicates at time , earthwork balance data for the entire reclamation area.

[0012] As a preferred solution of the method for determining the reclamation elevation of coal mining subsided cultivated land in the mining area of ​​the present invention, the reclamation construction elevation is dynamically adjusted according to the earthwork balance data and the pre-processed land use type, and the expression is: According to the earthwork balance data and the pre-processed land use type, the reclamation construction elevation is dynamically adjusted. The expression is: ; in, Indicates at time When the reclamation construction elevation is adjusted, represents the earthwork volume influencing factor, Indicates the adjustment factor for land use type.

[0013] As a preferred solution of the method for determining the reclamation elevation of coal mining subsided farmland in the mining area of ​​the present invention, the reclamation project is carried out according to the adjusted reclamation construction elevation, and the specific steps are as follows: According to the adjusted reclamation construction elevation, the elevation of each sub-area is cut and lowered; After the height reduction and lowering operations are completed, a laser land leveler is used to measure the surface flatness of each terrace sub-area; When the flatness is qualified, the reclaimed soil is loosened, drained and the pH value is adjusted, and vegetation is selected for cultivation according to the land use type; Install land subsidence monitors, weather monitors and soil moisture sensors in the reclaimed areas to monitor surface subsidence, soil quality and crop growth; Through long-term monitoring of surface subsidence, soil quality and crop growth, the quality and flatness of reclaimed arable land are evaluated and an annual assessment report is generated.

[0014] In a second aspect, the present invention provides a computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: when the computer program is executed by the processor, it implements any step of the method for determining the elevation of reclamation of cultivated land subsided by coal mining in a mining area as described in the first aspect of the present invention.

[0015] In a third aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein: when the computer program is executed by a processor, any step of the method for determining the elevation of reclamation of cultivated land caused by coal mining subsidence in a mining area as described in the first aspect of the present invention is implemented.

[0016] The beneficial effect of the present invention is that the present invention proposes an accurate method for determining the reclamation elevation of cultivated land subsided by coal mining in mining areas through satellite image analysis and probability integration method. Firstly, the original surface elevation data of the mining area before mining is obtained by using satellite images, which significantly improves the accuracy and efficiency of data acquisition, avoids the errors of traditional manual measurement, and is suitable for rapid assessment of large-scale mining areas. Secondly, through the probability integration method based on coal mining plan and underground coal occurrence conditions, the surface subsidence at different stages is accurately predicted, which provides a scientific basis for the dynamic adjustment of reclamation elevation. This innovative step effectively improves the accuracy of subsidence prediction, ensures the timely adjustment of construction elevation during the reclamation process, avoids engineering errors caused by subsidence changes, and reduces resource waste. Generally speaking, the present invention optimizes the reclamation design and construction process, improves the reclamation effect and engineering quality, and has significant economic benefits and environmental protection significance. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0018] Figure 1 This is a flow chart of the method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area in Example 1.

[0019] Figure 2 This is a schematic diagram of generating subsidence surface elevation data in Example 1. DETAILED DESCRIPTION

[0020] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are described in detail below in conjunction with the accompanying drawings.

[0021] In the following description, many specific details are set forth to facilitate a full understanding of the present invention, but the present invention may also be implemented in other ways different from those described herein, and those skilled in the art may make similar generalizations without violating the connotation of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below.

[0022] Secondly, the term "one embodiment" or "embodiment" as used herein refers to a specific feature, structure, or characteristic that may be included in at least one implementation of the present invention. The term "in one embodiment" that appears in different places in this specification does not necessarily refer to the same embodiment, nor does it refer to a separate or selective embodiment that is mutually exclusive with other embodiments.

[0023] Example 1, reference Figure 1 and Figure 2 , which is the first embodiment of the present invention, provides a method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area, comprising the following steps: S1. Obtain mining subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundary and land use type, and perform preprocessing; Obtain mining subsidence areas and farmland damage boundaries through drone aerial surveys; Obtain coal mining plans through the mine management department and historical data; The coal mining plan refers to the coal mining speed, coal mining depth and coal mining area; Obtain underground coal occurrence conditions through geological exploration reports and drilling data provided by exploration units; The underground coal occurrence conditions refer to the thickness, inclination and depth of the coal seam; Obtain land use types through land use planning maps provided by government departments and field surveys; The land use types mentioned are agricultural land, industrial land and ecological restoration land; Use GIS software to coordinate and time standardize the mining area subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundaries and land use types; Interpolation methods are used to fill in missing data on the mining subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundaries and land use types after coordinate standardization and time standardization. Normalize the mining subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundary and land use type after filling the data; By integrating the mining area subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundaries and land use types, it provides comprehensive data required for reclamation projects; The use of interpolation and normalization methods can improve the integrity and consistency of data, effectively reduce manual intervention, and improve work efficiency and data credibility.

[0024] S2. According to the pre-processed mining area subsidence area, the original surface elevation data of the mining area before mining is obtained through satellite image analysis; According to the pre-processed mining subsidence area, the WorldView-3 satellite was selected to obtain two images of different perspectives before mining, which were used as the left image and the right image. Performing geometric correction processing on the acquired left image and right image; According to the geometrically corrected left and right images, the focal length of the satellite sensor and the baseline distance during image acquisition are obtained; Use the left and right images after geometric correction to perform stereo image matching and calculate the disparity data of each pixel , the expression is:

[0025] in, Indicates location The parallax at represents the left image, Indicates the position in the left image The horizontal coordinate at represents the right image, Indicates the position in the right image The horizontal coordinate at .

[0026] Based on the parallax data, the focal length of the satellite sensor and the baseline distance when the image was collected, the surface elevation of each location is calculated to generate a digital elevation model DEM, which is expressed as: ; in, Indicated in Surface elevation data for the location, Indicates the baseline distance when the satellite image was taken. represents the focal length of the satellite sensor, represents the horizontal coordinate in the space coordinate, Represents the ordinate in the spatial coordinate; The median filter method is used to denoise and filter the generated digital elevation model DEM; For the digital elevation model DEM after denoising and filtering, the missing data is filled by interpolation method; Export the digital elevation model DEM after filling in the missing data into a standard geographic information format as the original surface elevation data before mining; Through satellite images and geometric correction methods, high-precision surface data of large-scale mining areas can be quickly obtained, avoiding the high cost and low efficiency of traditional manual measurement; The surface elevation data obtained from satellite images does not require on-site measurement, which effectively avoids the entry of personnel in complex or dangerous environments and ensures the safety and reliability of data collection; Through stereo image matching and parallax calculation, detailed surface elevation change data can be provided, thereby more accurately reflecting the original surface morphology of the mining area and providing strong support for subsequent subsidence prediction.

[0027] S3. Based on the pre-processed coal mining plan and underground coal occurrence conditions, the probability integral method is used to predict the surface subsidence data at different stages, and the predicted subsidence data is superimposed with the original terrain elevation data to generate the subsidence surface elevation data; Based on the pre-processed coal mining plan and underground coal occurrence conditions, the coal mining speed, coal mining depth, coal mining area, coal seam thickness, coal seam inclination and coal seam depth are obtained. Based on the obtained data, the probability integral method is used to predict the surface subsidence data at different stages. The expression is: ; in, Indicates at time Time, Location Surface subsidence data at It indicates the conversion ratio of subsidence caused by coal mining. Indicates the coal mining speed, Indicates the thickness of the coal seam, Indicates the depth of coal mining, Indicates the coal mining area, Indicates the depth of the coal seam, represents the spatial distribution function of subsidence, represents the sedimentation diffusion coefficient, represents the exponential function, represents the base of natural logarithms, represents the sinking time function, represents the sinking rate coefficient, Indicates the time of mining. represents the geological influence factor function, represents the cosine function, Indicates the coal seam inclination, Indicates the inclination of coal seam The cosine value of represents the geological influence factor; Superimpose the predicted subsidence data with the original surface elevation data to generate subsidence surface elevation data , the expression is: ; in, Indicates at time ,Location The subsidence surface elevation data at The probability integral method can accurately predict the surface subsidence at different stages of the coal mining process through coal seam thickness, coal seam inclination, coal seam depth, coal mining speed, coal mining depth and coal mining area, providing a scientific basis for reclamation projects; By superimposing the predicted subsidence data with the original terrain elevation data, more accurate subsidence surface elevation data can be obtained, providing a precise basis for the calculation of subsequent reclamation construction elevations and earthwork volumes, and reducing re-repairs caused by errors.

[0028] S4. Based on the subsidence surface elevation data and the pre-processed cultivated land damage boundary, the reclamation area is divided into several terraced sub-areas; Obtaining the damage boundary of cultivated land through drone aerial survey ; The subsidence surface elevation data was rasterized by GIS software, and the rasterized subsidence surface elevation data and the boundary of the damaged cultivated land area were spatially superimposed to obtain the rasterized subsidence surface elevation data in the reclamation area. The expression is: ; Based on the rasterized subsidence surface elevation data in the reclamation area, the slope value of each grid cell is calculated; ; in, Indicates the time Time, Location The slope value at represents the symbol of partial derivative, and Respectively expressed in The elevation at the location is and rate of change of direction; Set the slope threshold range to ; in, Indicates the lower limit of the slope range, Indicates the upper limit of the slope range; Through the reclassification function of GIS software, the slope value of each grid cell in the reclamation area is divided into Reclassify the grid cells, and the classified grid cells will be marked with the terraced sub-area to which each grid cell belongs; when < When ; when ≤ ≤ When ; when > When ; Set the side length of the grid cell to , then the area of ​​the grid cell is ; By combining the subsidence surface elevation data with the damaged boundaries of cultivated land and using GIS technology for spatial superposition and analysis, the reclaimed area can be automatically and efficiently divided into terraced sub-areas, avoiding the deviations and errors that occur during manual division.

[0029] S5. Calculate the earthwork volume of the terraced sub-areas according to the terraced sub-areas, and generate earthwork balance data; Set the reclamation construction elevation to ; For each grid unit, when the gridded subsidence surface elevation data in the reclamation area is greater than the reclamation construction elevation, the grid unit belongs to the height-cutting area, and the height-cutting earthwork volume is calculated. , expressed as: ; in, Represents a grid cell In time The amount of earthwork cut at that time; For each grid unit, when the gridded subsidence surface elevation data in the reclamation area is lower than the reclamation construction elevation, the grid unit belongs to the padding area, and the padding earthwork volume is calculated. , the expression is: ; in, Represents a grid cell In time The amount of earthwork lowered when The earthwork balance data is generated by accumulating the earthwork volume of all grid cells, and the expression is: ; in, Indicates at time , earthwork balance data for the entire reclamation area; By calculating the amount of earthwork to be cut and filled for each terraced sub-area, the required earthwork volume can be accurately assessed to avoid excessive construction, thus providing a scientific basis for the rational planning of reclamation construction. By generating earthwork balance data, the earthwork operation plan can be adjusted in real time to ensure maximum efficiency in earthwork use and avoid uneven earthwork allocation.

[0030] S6. Dynamically adjust the reclamation construction elevation based on earthwork balance data and the pre-processed land use type; According to the earthwork balance data and the pre-processed land use type, the reclamation construction elevation is dynamically adjusted. The expression is: ; in, Indicates at time When the reclamation construction elevation is adjusted, represents the earthwork volume influencing factor, An adjustment factor indicating the land use type; Land use types refer to agricultural land, industrial land, and ecological restoration land; Depending on the type of land, Then take different values; For agricultural land, , indicating that the adjustment of agricultural land has a small impact; For industrial land, , indicating that the adjustment of industrial land has a greater impact; For ecological restoration land, It indicates that the adjustment of ecological restoration land has the greatest impact; Combine earthwork balance data and land use type to dynamically adjust the reclamation construction elevation. The construction plan can be adjusted in real time according to the changes in earthwork construction progress and land use needs, thereby improving construction flexibility and accuracy. By dynamically adjusting the reclamation construction elevation, we can maximize the adaptation to the needs of different land types, optimize the function of land reclamation, and ensure that the reclaimed land can meet subsequent agricultural and vegetation requirements.

[0031] S7. Carry out the reclamation project according to the adjusted reclamation construction elevation; According to the adjusted reclamation construction elevation, the elevation of each sub-area is cut and lowered; After the height reduction and lowering operations are completed, a laser land leveler is used to measure the surface flatness of each terrace sub-area; When the flatness is qualified, the reclaimed soil is loosened, drained and the pH value is adjusted, and vegetation is selected for cultivation according to the land use type; Install land subsidence monitors, weather monitors and soil moisture sensors in the reclaimed areas to monitor surface subsidence, soil quality and crop growth; Evaluate the quality and flatness of reclaimed farmland through long-term monitoring of surface subsidence, soil quality and crop growth, and generate annual assessment reports; According to the adjusted reclamation construction elevation, accurate height reduction and lowering operations were carried out to ensure that the reclamation elevation of each sub-area met the planning requirements and to avoid deviations during the construction process to the greatest extent; By installing equipment such as ground subsidence monitors, meteorological monitors and soil moisture sensors in the reclamation area, it is possible to monitor the surface subsidence, soil quality and crop growth conditions in the reclamation area in real time, identify problems in a timely manner and make adjustments to ensure the quality of reclamation.

[0032] This embodiment also provides a computer device, which is suitable for the method of determining the elevation of reclamation of cultivated land subsided by coal mining in mining areas, including: a memory and a processor; the memory is used to store computer executable instructions, and the processor is used to execute computer executable instructions to implement the method of determining the elevation of reclamation of cultivated land subsided by coal mining in mining areas as proposed in the above embodiment. The computer device may be a terminal, and the computer device includes a processor, a memory, a communication interface, a display screen and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The communication interface of the computer device is used to communicate with an external terminal in a wired or wireless manner, and the wireless manner can be achieved through WIFI, an operator network, NFC (near field communication) or other technologies. The display screen of the computer device may be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer device may be a touch layer covered on the display screen, or a key, trackball or touchpad provided on the housing of the computer device, or an external keyboard, touchpad or mouse, etc.

[0033] This embodiment also provides a storage medium on which a computer program is stored. When the program is executed by a processor, the method for determining the elevation of cultivated land reclamation in a mining area due to coal mining subsidence is implemented as proposed in the above embodiment; the storage medium can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as a static random access memory (SRAM), an electrically erasable programmable read-only memory (EEPROM), an erasable programmable read-only memory (EPROM), a programmable read-only memory (PROM), a read-only memory (ROM), a magnetic storage device, a flash memory, a disk or an optical disk.

[0034] In summary, the present invention proposes a precise method for determining the reclamation elevation of cultivated land subsided by coal mining in mining areas through satellite image analysis and probability integration method. Firstly, the original surface elevation data of the mining area before mining is obtained by using satellite images, which significantly improves the accuracy and efficiency of data acquisition and avoids the errors of traditional manual measurement. It is suitable for rapid assessment of large-scale mining areas. Secondly, through the probability integration method based on coal mining plan and underground coal occurrence conditions, the surface subsidence at different stages is accurately predicted, which provides a scientific basis for the dynamic adjustment of the reclamation elevation. This innovative step effectively improves the accuracy of subsidence prediction, ensures the timely adjustment of construction elevation during the reclamation process, avoids engineering errors caused by subsidence changes, and reduces resource waste. In general, the present invention optimizes the reclamation design and construction process, improves the reclamation effect and engineering quality, and has significant economic benefits and environmental protection significance.

[0035] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.

Claims

1. A method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area, characterized by: include, Obtain mining subsidence areas, coal mining plans, underground coal occurrence conditions, cultivated land damage boundaries and land use types, and perform preprocessing; According to the pre-processed mining subsidence area, the original surface elevation data of the mining area before mining is obtained through satellite image analysis; Based on the pre-processed coal mining plan and underground coal occurrence conditions, the probability integral method is used to predict the surface subsidence data at different stages, and the predicted subsidence data is superimposed with the original terrain elevation data to generate the subsidence surface elevation data; According to the subsidence surface elevation data and the pre-processed cultivated land damage boundary, the reclamation area is divided into several terraced sub-areas; According to the terrace sub-areas, calculate the earthwork volume of the terrace sub-areas and generate earthwork balance data; According to the earthwork balance data and the pre-processed land use type, the reclamation construction elevation is dynamically adjusted; Carry out the reclamation project according to the adjusted reclamation construction elevation.

2. The method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area according to claim 1, characterized in that: The specific steps of obtaining the mining area subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundary and land use type, and performing preprocessing are as follows: Obtain mining subsidence areas and farmland damage boundaries through drone aerial surveys; Obtain coal mining plans through the mine area management department and historical data review; Obtain underground coal occurrence conditions through geological exploration reports and drilling data provided by exploration units; Obtain land use types through land use planning maps provided by government departments and field surveys; Use GIS software to coordinate and time standardize the mining area subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundaries and land use types; Interpolation methods are used to fill in missing data on the mining subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundaries and land use types after coordinate standardization and time standardization. The mining subsidence area, coal mining plan, underground coal occurrence conditions, cultivated land damage boundary and land use type after filling the data will be normalized.

3. The method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area according to claim 2, characterized in that: The original surface elevation data before mining of the mining area is obtained by satellite image analysis based on the pre-processed mining area subsidence area. The specific steps are as follows: According to the pre-processed mining subsidence area, the WorldView-3 satellite was selected to obtain two images of different perspectives before mining, which were used as the left image and the right image. Performing geometric correction processing on the acquired left image and right image; According to the geometrically corrected left and right images, the focal length of the satellite sensor and the baseline distance during image acquisition are obtained; Use the left and right images after geometric correction to perform stereo image matching and calculate the disparity data of each pixel ; Based on the parallax data, the focal length of the satellite sensor and the baseline distance when the image was collected, the surface elevation of each location is calculated to generate a digital elevation model DEM, which is expressed as: ； in, Indicated in Surface elevation data for the location, Indicates the baseline distance when the satellite image was taken. represents the focal length of the satellite sensor, represents the horizontal coordinate in the space coordinate, Represents the ordinate in the spatial coordinate; The median filter method is used to denoise and filter the generated digital elevation model DEM; For the digital elevation model DEM after denoising and filtering, the missing data is filled by interpolation method; The digital elevation model DEM after filling in the missing data is exported to a standard geographic information format as the original surface elevation data before mining.

4. The method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area according to claim 3, characterized in that: The method uses the probability integral method to predict the surface subsidence data at different stages based on the pre-processed coal mining plan and underground coal occurrence conditions, and superimposes the predicted subsidence data with the original terrain elevation data to generate the subsidence surface elevation data. The specific steps are as follows: Based on the pre-processed coal mining plan and underground coal occurrence conditions, the probability integral method is used to predict the surface subsidence data at different stages. The expression is: ； in, Indicates at time Time, Location Surface subsidence data at It indicates the conversion ratio of subsidence caused by coal mining. Indicates the coal mining speed, Indicates the thickness of the coal seam, Indicates the depth of coal mining, Indicates the coal mining area, Indicates the depth of the coal seam, represents the spatial distribution function of subsidence, represents the sedimentation diffusion coefficient, Mathematical notation for exponential functions, represents the base of natural logarithms, represents the sinking time function, represents the sinking rate coefficient, Indicates the time of mining. represents the geological influence factor function, represents the cosine function, Indicates the coal seam inclination, Indicates the inclination of coal seam The cosine value of represents the geological influence factor; Superimpose the predicted subsidence data with the original surface elevation data to generate subsidence surface elevation data .

5. The method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area according to claim 4, characterized in that: According to the subsidence surface elevation data and the pre-processed cultivated land damage boundary, the reclamation area is divided into several terraced sub-areas. The specific steps are as follows: Obtaining the damage boundary of cultivated land through drone aerial survey ; The subsidence surface elevation data was rasterized by GIS software, and the rasterized subsidence surface elevation data and the boundary of the damaged cultivated land area were spatially superimposed to obtain the rasterized subsidence surface elevation data in the reclamation area. ; Based on the rasterized subsidence surface elevation data in the reclamation area, the slope value of each grid cell is calculated; ； in, Indicates the time Time, Location The slope value at represents the symbol of partial derivative, and Respectively expressed in The elevation at the location is and rate of change of direction; Set the slope threshold range to ; in, Indicates the lower limit of the slope range, Indicates the upper limit of the slope range; Through the reclassification function of GIS software, the slope value of each grid cell in the reclamation area is divided into Reclassify the grid cells, and the classified grid cells will be marked with the terraced sub-area to which each grid cell belongs; when < When ; when ≤ ≤ When ; when > When ; Set the side length of the grid cell to , then the area of ​​the grid cell is .

6. The method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area according to claim 5, characterized in that: The specific steps of calculating the earthwork volume of the terraced sub-area and the earthwork balance data according to the terraced sub-area are as follows: Set the reclamation construction elevation to ; For each grid unit, when the gridded subsidence surface elevation data in the reclamation area is greater than the reclamation construction elevation, the grid unit belongs to the height-cutting area, and the height-cutting earthwork volume is calculated. ; For each grid unit, when the gridded subsidence surface elevation data in the reclamation area is lower than the reclamation construction elevation, the grid unit belongs to the padding area, and the padding earthwork volume is calculated. ; The earthwork balance data is generated by accumulating the earthwork volume of all grid cells, and the expression is: ； in, Indicates at time , earthwork balance data for the entire reclamation area.

7. The method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area according to claim 6, characterized in that: The reclamation construction elevation is dynamically adjusted based on the earthwork balance data and the pre-processed land use type. The expression is: According to the earthwork balance data and the pre-processed land use type, the reclamation construction elevation is dynamically adjusted. The expression is: ； in, Indicates at time When the reclamation construction elevation is adjusted, represents the earthwork volume influencing factor, Indicates the adjustment factor for land use type.

8. The method for determining the elevation of reclamation of farmland subsided by coal mining in a mining area according to claim 7, characterized in that: The reclamation project is carried out according to the adjusted reclamation construction elevation, and the specific steps are as follows: According to the adjusted reclamation construction elevation, the elevation of each sub-area is cut and lowered; After the height reduction and lowering operations are completed, a laser land leveler is used to measure the surface flatness of each terrace sub-area; When the flatness is qualified, the reclaimed soil is loosened, drained and the pH value is adjusted, and vegetation is selected for cultivation according to the land use type; Install land subsidence monitors, weather monitors and soil moisture sensors in the reclaimed areas to monitor surface subsidence, soil quality and crop growth; Through long-term monitoring of surface subsidence, soil quality and crop growth, the quality and flatness of reclaimed arable land are evaluated and an annual assessment report is generated.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method for determining the elevation of cultivated land reclamation subsided by coal mining in a mining area as described in any one of claims 1 to 8 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for determining the elevation of cultivated land reclamation subsided by coal mining in a mining area as described in any one of claims 1 to 8 are implemented.

Images