Method and system for predicting the boundary of a coal mining subsidence area considering the additional effect of aeolian sand

Abstract

The application discloses a kind of considering wind sand additional effect's coal mining subsidence area boundary prediction method and system, belong to subsidence area boundary prediction field.The prediction method includes: obtaining the synthetic aperture radar image of coal mining area;According to synthetic aperture radar image, the first coal mining subsidence area boundary caused by monitoring coal mining;According to the mapping model of preposition and the wind speed and sand particle size of coal mining area, the wind sand movement boundary caused by wind sand migration is calculated, wherein, mapping model is the model for representing the mapping relationship between wind speed, sand particle size and wind sand movement boundary;The first coal mining subsidence area boundary is corrected using wind sand movement boundary, and the second coal mining subsidence area boundary is obtained.The method and system of the application improve the prediction accuracy of coal mining subsidence area boundary, and have higher dynamic adaptability.

Description

TECHNICAL FIELD

[0001] The present application belongs to the field of subsidence area boundary prediction, and more particularly relates to a coal mining subsidence area boundary prediction method and system considering the additional effect of aeolian sand. BACKGROUND

[0002] In the aeolian sand area, there is a coupling effect between the surface subsidence caused by coal mining and the migration of aeolian sand, and the migration of aeolian sand can cause changes in the boundary of the coal mining subsidence area. The traditional prediction method usually uses a level gauge, a total station, a Global Navigation Satellite System (GNSS) and the like to monitor the boundary of the coal mining subsidence area, and has the problems of low monitoring efficiency, small monitoring range and great influence of climate.

[0003] In order to solve the problems existing in the traditional prediction method, the current research usually uses Interferometric Synthetic Aperture Radar (InSAR) technology to monitor the boundary of the coal mining subsidence area, such as Differential InSAR (D-InSAR) technology, Small Baseline Subset InSAR (SBAS-InSAR) technology and the like. In addition, the existing research mainly focuses on the coal mining subsidence area without the influence of aeolian sand, and the existing prediction theory does not take into account the influence of the coupling effect of surface movement and deformation and aeolian sand movement on the prediction of the boundary of the coal mining subsidence area, thereby causing large prediction deviation. Therefore, how to couple the influence of multiple factors to improve the prediction accuracy of the boundary of the subsidence area has important research significance. SUMMARY

[0004] In view of the deficiencies of the prior art, the purpose of the present application is to provide a coal mining subsidence area boundary prediction method and system considering the additional effect of aeolian sand, so as to solve the problem of low prediction accuracy of the boundary of the subsidence area caused by only considering a single factor in the prior art.

[0005] In order to achieve the above purpose, a first aspect of the present application provides a coal mining subsidence area boundary prediction method considering the additional effect of aeolian sand, comprising: acquiring a synthetic aperture radar image of a coal mining area; monitoring a first coal mining subsidence area boundary caused by coal mining according to the synthetic aperture radar image; calculating an aeolian sand movement boundary caused by aeolian sand migration according to a preset mapping model and the wind speed and sand particle size of the coal mining area, wherein the mapping model is a model representing the mapping relationship between the wind speed, the sand particle size and the aeolian sand movement boundary; and correcting the first coal mining subsidence area boundary by using the aeolian sand movement boundary to obtain a second coal mining subsidence area boundary.

[0006] Optionally, the establishing process of the mapping model comprises: simulating aeolian sand experiments under different wind speeds and different sand particle size combinations, fitting the aeolian sand moving boundary under each aeolian sand experiment as a corresponding mathematical model; extracting model parameters of the mathematical model; and establishing the mapping model according to a plurality of groups of corresponding wind speed, sand particle size and model parameters obtained from the aeolian sand experiments.

[0007] Optionally, the mathematical model is a parabolic model, and parameters of the parabolic model comprise a constant term parameter, a first-order term coefficient and a second-order term coefficient; an upper boundary of the aeolian sand moving boundary is fitted as a first parabolic model, and a lower boundary is fitted as a second parabolic model.

[0008] Optionally, the method further comprises: setting a measuring device according to a preset geometric similarity ratio and a preset time similarity ratio, so as to simulate the aeolian sand experiment by using the measuring device.

[0009] Optionally, the correction of the first coal mining subsidence area boundary by using the aeolian sand moving boundary specifically comprises: superimposing the aeolian sand moving boundary and the first coal mining subsidence area boundary, wherein the center of the aeolian sand moving boundary coincides with the center of the first coal mining subsidence area boundary after superimposition, and the aeolian sand moving boundary is located in a first plane, and the first plane is determined by a working face strike direction and a working face dip direction within the first coal mining subsidence area boundary.

[0010] Optionally, the first coal mining subsidence area boundary caused by coal mining is monitored according to the synthetic aperture radar image, specifically comprising: processing the synthetic aperture radar image to generate a time-series subsidence rate map and a cumulative subsidence map; calculating a prediction parameter according to the time-series subsidence rate map and the cumulative subsidence map; combining the prediction parameter, and predicting a surface subsidence amount, a surface horizontal movement amount, a surface inclination deformation amount, a surface curvature deformation amount and a surface horizontal deformation amount by using a probability integral method; and generating the first coal mining subsidence area boundary according to a prediction result.

[0011] Optionally, the surface horizontal movement amount is:

[0012] ;

[0013] ;

[0014] wherein, , are an axial component and an axial component of the surface horizontal movement amount respectively, is a maximum surface horizontal movement amount, is a first integral variable, is a second integral variable, is a main influence radius, is a subsidence amount of the ground surface, is a coal mining influence propagation angle.

[0015] The second aspect of the present application provides a system for predicting a coal mining subsidence area boundary considering the additional effect of aeolian sand, comprising: an acquisition module configured to acquire a synthetic aperture radar image of a coal mining area; a prediction module configured to monitor a first coal mining subsidence area boundary caused by coal mining according to the synthetic aperture radar image; a calculation module configured to calculate an aeolian sand movement boundary caused by aeolian sand movement according to a preset mapping model and a wind speed and a sand particle size of the coal mining area, wherein the mapping model is a model representing a mapping relationship between the wind speed, the sand particle size and the aeolian sand movement boundary; and a correction module configured to correct the first coal mining subsidence area boundary using the aeolian sand movement boundary to obtain a second coal mining subsidence area boundary.

[0016] The third aspect of the present application provides an electronic device, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the method as described above when executing the program.

[0017] The fourth aspect of the present application provides a non-transitory computer readable storage medium, which stores computer instructions for causing a computer to execute the method as described above.

[0018] Compared with the prior art, the present application has the advantages of providing a method and system for predicting a coal mining subsidence area boundary considering the additional effect of aeolian sand, constructing a mapping model to accurately quantify the influence of wind speed and sand particle size on the subsidence range, dynamically predicting the boundary evolution of the aeolian sand movement area in the coal mining area based on the mapping model, and correcting the subsidence area boundary caused by coal mining using the aeolian sand movement boundary, thereby coupling the influence of both aeolian sand and coal mining to predict the subsidence area boundary, improving the prediction accuracy of the subsidence area boundary, being more beneficial to the design of the coal mining working face, having better engineering practicability, being applicable to subsidence area prediction under different climate conditions, and having higher dynamic adaptability. BRIEF DESCRIPTION OF DRAWINGS

[0019] Figure 1 The flowchart of the method for predicting a coal mining subsidence area boundary considering the additional effect of aeolian sand provided by the embodiments of the present application.

[0020] Figure 2 The flowchart of processing the synthetic aperture radar image provided by the embodiments of the present application.

[0021] Figure 3 The deformation rate provided by the embodiments of the present application.

[0022] Figure 4 A timing analysis diagram provided for an embodiment of the present application.

[0023] Figure 5 A surface subsidence result obtained by probability integral method provided for an embodiment of the present application.

[0024] Figure 6 A dune sand moving boundary under simulation of dune sand experiment provided for an embodiment of the present application.

[0025] Figure 7 A dune sand moving boundary fitting result provided for an embodiment of the present application.

[0026] Figure 8 A subsidence area boundary correction result schematic diagram provided for an embodiment of the present application.

[0027] Figure 9 A block diagram of a coal mining subsidence area boundary prediction system considering the additional effect of dune sand provided for an embodiment of the present application.

[0028] Figure 10 A schematic diagram of an electronic device provided for an embodiment of the present application. DETAILED DESCRIPTION

[0029] In view of the deficiencies in the prior art, the present inventors have long studied and practiced to come up with the technical solution of the present application. The technical solution, its implementation process and principles will be further explained as follows.

[0030] In the following description, many specific details are set forth in order to provide a thorough understanding of the present application. However, the present application can be practiced in a variety of ways beyond the specific embodiments described herein without departing from the scope of the present application, which is set forth in the appended claims.

[0031] In addition, in the description of the present application, it should be understood that the terms "upper", "lower", "inner", "outer", "horizontal", "vertical" and the like indicate the orientation or positional relationship based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and do not indicate or imply that the device or element referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore should not be understood as a limitation of the present application.

[0032] In the description of the present specification, the description referring to the terms "one embodiment", "an embodiment", "the embodiment", and the like means that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described can be combined in any appropriate manner in any one or more embodiments or examples.

[0033] The embodiment of the present application provides a coal mining subsidence area boundary prediction method considering the additional effect of aeolian sand. Referring to Figure 1 , the coal mining subsidence area boundary prediction method considering the additional effect of aeolian sand comprises steps S100-S400.

[0034] Step S100, acquiring a synthetic aperture radar image of a coal mining area.

[0035] Specifically, for example, the synthetic aperture radar image of the coal mining area is collected by using the synthetic aperture radar. The synthetic aperture radar has high resolution and can work all day long, and can obtain a synthetic aperture radar image with extremely high resolution.

[0036] Step S200, monitoring a first coal mining subsidence area boundary caused by coal mining according to the synthetic aperture radar image.

[0037] Specifically, for example, the D-InSAR technology or the SBAS-InSAR technology is used to process the synthetic aperture radar image, and the first coal mining subsidence area boundary caused by coal mining is obtained. The first coal mining subsidence area boundary is the subsidence boundary under the action of only the coal mining factor.

[0038] Step S300, calculating an aeolian sand moving boundary caused by aeolian sand migration according to a preset mapping model and the wind speed and sand particle size of the coal mining area, wherein the mapping model is a model representing the mapping relationship between the wind speed, the sand particle size, and the aeolian sand moving boundary.

[0039] The mapping model represents the relationship among the wind speed v , the sand particle size r , and the aeolian sand moving boundary. The aeolian sand moving boundary is the moving boundary under the action of the aeolian sand factor. According to the pre-established mapping model, the wind speed v and the sand particle size r of the coal mining area are known, and the aeolian sand moving boundary caused by the aeolian sand migration of the coal mining area can be calculated.

[0040] Step S400, correcting the first coal mining subsidence area boundary by using the aeolian sand moving boundary to obtain a second coal mining subsidence area boundary.

[0041] The first coal mining subsidence area boundary considering only the effect of coal mining is corrected by using the aeolian sand moving boundary, and the second coal mining subsidence area boundary considering the additional effect of aeolian sand is obtained, that is, the subsidence area boundary considering the coupling effect of aeolian sand and coal mining is obtained, the prediction accuracy of the subsidence area boundary is improved, the coal mining working face design is more favorable, the engineering practicability is better, and the subsidence area prediction under different climate conditions is applicable, and the dynamic adaptability is higher.

[0042] Before step S300 is performed, the method further includes: establishing a mapping model representing the mapping relationship between the wind speed, the sand particle size and the aeolian sand moving boundary. Preferably, in an embodiment, the establishment process of the mapping model includes the following steps S1-S3.

[0043] Step S1, simulate aeolian sand experiments under different wind speed and different sand particle size combinations, and fit the aeolian sand moving boundary under each aeolian sand experiment into a corresponding mathematical model.

[0044] Preferably, the measuring device is set according to a preset geometric similarity ratio and a preset time similarity ratio to simulate the aeolian sand experiment by using the measuring device. The preset geometric similarity ratio is, for example, 1:4000, 1:5000 or 1:6000, etc.; the preset time similarity ratio is, for example, 1:6, 1:7 or 1:8, etc.; the different wind speeds are, for example, a plurality of wind speeds in the range of 4.4-5.7 m / s; and the different sand particle sizes are, for example, a plurality of sand particle sizes in the range of 0.18-1.25 mm.

[0045] Step S2, extract the model parameters of the mathematical model.

[0046] Preferably, the mathematical model is a parabolic model, and the parameters of the parabolic model include a constant term parameter, a first-order term coefficient and a second-order term coefficient; the upper boundary of the aeolian sand moving boundary is fitted into a first parabolic model, and the lower boundary is fitted into a second parabolic model.

[0047] Taking the parabolic model as an example, the model parameters include: the first parabolic model parameters (constant term parameter u 1, first-order term coefficient u 2 and second-order term coefficient u 3) representing the upper boundary of the aeolian sand moving boundary, and the second parabolic model parameters (constant term parameter d 1, first-order term coefficient d 2 and second-order term coefficient d 3) representing the lower boundary of the aeolian sand moving boundary.

[0048] It should be noted that the aeolian sand moving boundary can also be fitted into other mathematical models, such as a linear regression model.

[0049] Step S3, according to the wind sand experiment obtained from several groups of corresponding wind speed, sand particle size and model parameters, the mapping model is established.

[0050] Taking the mathematical model as a parabolic model, the mapping model includes: the functional relationship between wind speed v , sand particle size r and u 1; the functional relationship between wind speed v , sand particle size r and u 2; the functional relationship between wind speed v , sand particle size r and u 3; the functional relationship between wind speed v , sand particle size r and d 1; the functional relationship between wind speed v , sand particle size r and d 2; the functional relationship between wind speed v , sand particle size r and d 3. According to the above functional relationship, the mapping relationship between wind speed v , sand particle size r and the moving boundary of wind sand can be established.

[0051] Preferably, in an embodiment, in step S200, the first coal mining subsidence area boundary caused by coal mining is monitored according to the synthetic aperture radar image, which specifically includes the following steps S210-S240.

[0052] Step S210, the synthetic aperture radar image is processed to generate a time series subsidence rate map and a cumulative subsidence map.

[0053] Specifically, for example, the synthetic aperture radar image is processed by using D-InSAR technology or SBAS-InSAR technology to generate a time series subsidence rate map and a cumulative subsidence map. Preferably, SBAS-InSAR technology is used because it has relatively high time series resolution and can better capture dynamic evolution.

[0054] The time series subsidence rate map is a kind of chart generated by long-term monitoring and data analysis of surface subsidence through time series analysis technology, which is used to show the subsidence rate change of the coal mining area in different time periods. The cumulative subsidence map reflects the cumulative subsidence change of the coal mining area with time, and directly shows the development trend of the subsidence.

[0055] Step S220, calculating the expected parameters according to the time series subsidence rate map and the cumulative subsidence map.

[0056] Expected parameters include, for example, subsidence coefficient, horizontal movement coefficient, tangent of the main influence angle, inflection point offset distance, and mining influence propagation angle.

[0057] Step S230: Based on the predicted parameters, the probability integral method is used to predict the amount of surface subsidence, the amount of surface horizontal movement, the amount of surface tilt deformation, the amount of surface curvature deformation, and the amount of surface horizontal deformation.

[0058] Preferably, in one embodiment, the land subsidence is:

[0059]

[0060] The horizontal movement of the earth's surface is:

[0061] ;

[0062]

[0063] The amount of surface tilt deformation is:

[0064] ;

[0065]

[0066] The amount of surface curvature deformation is:

[0067] ;

[0068]

[0069] The horizontal deformation of the ground surface is:

[0070] ;

[0071]

[0072] and They are respectively:

[0073] ;

[0074] ;

[0075] in, , These are the horizontal movement of the earth's surface. Axial components, Axial components, This represents the maximum horizontal movement on the Earth's surface. Let be the first integral variable. For the second integral variable, As the main influencing radius, is a subsidence amount of the ground surface, is a propagation angle of a coal mining influence; is a maximum subsidence value of the ground surface, , is a point in a coal mining area, , is an axis component of a ground surface tilt deformation amount, is an axis component of a ground surface tilt deformation amount, , is an axis component of a ground surface curvature deformation amount, is an axis component of a ground surface curvature deformation amount, , is an axis component of a ground surface horizontal deformation amount, is an axis component of a ground surface horizontal deformation amount, , is a tilt value of the ground surface, is a mining thickness, is a subsidence coefficient, is a coal seam inclination, is a horizontal movement coefficient.

[0076] Step S240, generating a first coal mining subsidence area boundary according to the prediction result.

[0077] Specifically, according to the ground surface subsidence amount, the ground surface horizontal movement amount, the ground surface tilt deformation amount, the ground surface curvature deformation amount and the ground surface horizontal deformation amount predicted in step S230, a first coal mining subsidence area boundary caused by coal mining is generated.

[0078] Preferably, in an embodiment, in step S400, the first coal mining subsidence area boundary is corrected by using the aeolian sand movement boundary, specifically including: superimposing the aeolian sand movement boundary and the first coal mining subsidence area boundary, after superimposition, the center of the aeolian sand movement boundary coincides with the center of the first coal mining subsidence area boundary, and the aeolian sand movement boundary is located in a first plane, the first plane is determined by the working face strike direction and the working face inclination direction within the first coal mining subsidence area boundary.

[0079] Through the above steps, the subsidence area boundary caused by considering both the aeolian sand and the coal mining is obtained, the prediction accuracy of the subsidence area boundary is improved, which is more conducive to the design of the coal mining working face, has better engineering practicability, and is suitable for subsidence area prediction under different climate conditions, and has higher dynamic adaptability.

[0080] ​​The technical solutions of the present application will be further described in detail below in combination with several preferred embodiments and the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work belong to the protection scope of the present application. The test methods in the following embodiments without specific conditions are usually according to the conventional conditions.

[0081] In combination with Figures 2-8 , the implementation process of the coal mining subsidence area boundary prediction method considering the additional effect of aeolian sand in the embodiments of the present application is specifically described, including the following steps 1-4.

[0082] Step 1: SBAS-InSAR data acquisition and processing.

[0083] Referring to Figure 2 , 78 scenes of Sentinel-1A SAR images (2017.4-2019.12) are taken as data sources, the spatial baseline threshold is ≤2%, and the time baseline is ≤90 days; the data sources are sequentially processed by image registration, differential interference generation, phase unwrapping, atmospheric error elimination, and high-coherence-point extraction to generate the deformation time series and average displacement rate in Figure 2 , and further generate the deformation rate graph shown in Figure 3 and the time series analysis graph shown in Figure 4 .

[0084] In this embodiment, the subsidence boundary evolution of the wind erosion area and the accumulation area is dynamically predicted. The wind erosion area and the accumulation area are divided as follows: according to the local main wind direction and the working face advancing direction, when the two directions are consistent, the wind erosion area is located on one side of the working face stop line, and the accumulation area is on the other side of the working face open cut; when the two directions are opposite, the wind erosion area is located on one side of the working face open cut, and the accumulation area is on the other side of the working face stop line.

[0085] Step 2: aeolian sand physical simulation experiment.

[0086] According to the existing material conditions and site conditions, the geometric similarity ratio of the aeolian sand physical simulation experiment is set to 1:5000, and the time similarity ratio is set to 1:7. Through the combination experiment of wind speed (4.4-5.7 m / s) and sand particle size (0.18-1.25 mm), the migration distance is measured. The aeolian sand migration range measurement experiment table is shown in Table 1.

[0087] Table 1

[0088]

[0089] The aeolian sand moving boundary obtained by the aeolian sand physical simulation experiment under different wind speeds and sand particle sizes is shown in Figure 5 , and the aeolian sand moving boundary is shown in Figure 5The wind-blown sand moving boundary mathematical model is obtained by fitting the shown boundary. Taking 0.18-0.25 and 0.27-0.55 as examples, the corresponding wind-blown sand moving boundary mathematical model is shown in Table 2, and the corresponding fitting curve is shown in Figure 6 .

[0090] Table 2

[0091]

[0092] According to the model parameters (parabolic parameters) in the obtained wind-blown sand moving boundary mathematical model and the corresponding wind speed and sand particle size, the mathematical model between the parabolic parameters and the wind speed and sand particle size can be inversely deduced and fitted, as shown in Table 3.

[0093] Table 3

[0094]

[0095] According to the measured wind speed and sand particle size, the values of the six parameters are calculated according to Table 3, and the wind-blown sand moving boundary caused by the wind-blown sand movement is obtained through the values of the six parameters.

[0096] Step 3: predicting the first coal mining subsidence area boundary caused by coal mining based on the probability integral method.

[0097] For the movement and deformation of any point (x, y) on the surface within the working face range of the horizontal or gently inclined coal seam, the surface subsidence , the surface horizontal movement , , the surface inclination deformation , , the surface curvature deformation , , and the surface horizontal deformation , are predicted. The prediction results are shown in Figure 7 , the 2201 working face is 2500m long and 300m wide, the 2101 working face is 2000m long and 300m wide, the mining depth is 725m, and the mining thickness is 6.5m, the probability integral parameters and the prediction results are shown in Table 4. Among them, q ′ is the surface subsidence rate, tan β is the tangent of the main influence angle, S1 is the offset distance of the strike inflection point, and S3 is the offset distance of the inclination inflection point.

[0098] Table 4

[0099]

[0100] Step 4: correcting the first coal mining subsidence area boundary by using the aeolian sand moving boundary. The center of the result predicted by the probability integral method is taken as the coordinate origin, the working face strike direction is taken as the horizontal axis, and the working face dip direction is taken as the vertical axis to establish a coordinate system. The aeolian sand moving boundary and the first coal mining subsidence area boundary are drawn into the same scene for superposition to correct the boundary, and the correction result is as shown in Figure 8

[0101] The coal mining subsidence area boundary prediction method considering the additional effect of aeolian sand provided in the embodiment can quantize the influence of wind speed and particle size on the boundary, and is suitable for different climate conditions, such as the northwest wind dominant area. The subsidence area boundary considering both the aeolian sand and the coal mining is obtained, the prediction accuracy of the subsidence area boundary is improved, and the coal mining working face design, such as the isolation coal pillar width, is more beneficial, and the ecological damage area is reduced.

[0102] Based on the same inventive concept, the present application also provides a coal mining subsidence area boundary prediction system considering the additional effect of aeolian sand, which corresponds to the method of any of the above embodiments. Referring to Figure 9 , the coal mining subsidence area boundary prediction system 800 considering the additional effect of aeolian sand comprises an acquisition module 810, a monitoring module 820, a calculation module 830 and a correction module 840.

[0103] The acquisition module 810, for example, performs step S100, is configured to acquire a synthetic aperture radar image of a coal mining area.

[0104] The monitoring module 820, for example, performs step S200, is configured to monitor a first coal mining subsidence area boundary caused by coal mining according to the synthetic aperture radar image.

[0105] The calculation module 830, for example, performs step S300, is configured to calculate an aeolian sand moving boundary caused by aeolian sand movement according to a preset mapping model and wind speed and sand particle size of the coal mining area, wherein the mapping model is a model representing the mapping relationship between the wind speed, the sand particle size and the aeolian sand moving boundary.

[0106] The correction module 840, for example, performs step S400, is configured to correct the first coal mining subsidence area boundary by using the aeolian sand moving boundary to obtain a second coal mining subsidence area boundary.

[0107] ​Through the cooperation of the above-mentioned various components, the mining subsidence area boundary prediction system considering the additional effect of aeolian sand provided by the embodiment of the application constructs a mapping model to accurately quantify the influence of wind speed and sand particle size on the subsidence range, dynamically predicts the boundary evolution of the aeolian sand moving area in the coal mining area based on the mapping model, corrects the subsidence area boundary caused by coal mining by using the aeolian sand moving boundary, thereby coupling the influence of the two factors of aeolian sand and coal mining to predict the subsidence area boundary, improving the prediction accuracy of the subsidence area boundary, being more conducive to the design of the coal mining working face, having better engineering practicability, and being applicable to subsidence area prediction under different climate conditions, and having higher dynamic adaptability.

[0108] Further function description of the above-mentioned various modules is the same as that of the above-mentioned corresponding method embodiment, and will not be described here.

[0109] Based on the same inventive concept, the present application also provides an electronic device corresponding to any of the above-mentioned embodiment methods, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to realize the mining subsidence area boundary prediction method considering the additional effect of aeolian sand according to any one of the above-mentioned embodiments.

[0110] Figure 10 A more specific electronic device hardware structure schematic diagram provided by the embodiment is shown, which can include a processor 910, a memory 920, an input / output interface 930, a communication interface 940 and a bus 950. The processor 910, the memory 920, the input / output interface 930 and the communication interface 940 are connected to each other through the bus 950 for communication within the device.

[0111] The processor 910 can be implemented by a general CPU (Central Processing Unit, central processor), a microprocessor, an application specific integrated circuit (Application Specific Integrated Circuit, ASIC) or one or more integrated circuits, etc., for executing related programs to realize the technical solutions provided by the embodiments of the present application.

[0112] The memory 920 can be implemented in the form of ROM (Read Only Memory, read-only memory), RAM (Random Access Memory, random access memory), static storage device, dynamic storage device, etc. The memory 920 can store an operating system and other application programs, and when the technical solutions provided by the embodiments of the present application are implemented by software or firmware, the related program codes are saved in the memory 920 and executed by the processor 910.

[0113] The input / output interface 930 is configured to connect an input / output module to realize information input and output. The input / output module can be configured in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. The input device can include a keyboard, a mouse, a touch screen, a microphone, various sensors, etc., and the output device can include a display, a speaker, a vibrator, an indicator light, etc.

[0114] The communication interface 940 is configured to connect a communication module (not shown in the figure) to realize communication interaction between the device and other devices. The communication module can realize communication through a wired manner (such as a USB, a network cable, etc.) or a wireless manner (such as a mobile network, WIFI, Bluetooth, etc.).

[0115] The bus 950 includes a channel to transmit information between various components (such as the processor 910, the memory 920, the input / output interface 930, and the communication interface 940) of the device.

[0116] It should be noted that although the above device only shows the processor 910, the memory 920, the input / output interface 930, the communication interface 940, and the bus 950, in the specific implementation process, the device can also include other components necessary for normal operation. In addition, those skilled in the art can understand that the above device can also only include components necessary for implementing the embodiments of the present specification, and does not necessarily include all components shown in the figure.

[0117] The electronic device of the above embodiments is used to implement the corresponding coal mining subsidence area boundary prediction method considering the additional effect of wind-blown sand in any of the above embodiments, and has the beneficial effects of the corresponding method embodiments, which are not repeated here.

[0118] Based on the same inventive concept, the present application also provides a non-transitory computer readable storage medium, which stores computer instructions for causing the computer to execute the coal mining subsidence area boundary prediction method considering the additional effect of wind-blown sand as described in any of the above embodiments.

[0119] The computer readable media of the embodiments can include permanent and non-permanent, removable and non-removable media, which can be implemented by any method or technology to store information. The information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette, magnetic tape, magnetic disk storage or other magnetic storage device, or any other non-transmission medium that can be used to store information accessible by a computing device.

[0120] The storage medium of the above embodiments stores computer instructions for causing the computer to execute the coal mining subsidence area boundary prediction method considering the additional action of sand dune, as described in any of the above embodiments, and has the beneficial effects of the corresponding method embodiments, which are not described here.

[0121] Those skilled in the art should understand that the above discussion of any of the embodiments is only exemplary; under the idea of the present application, the technical features of the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other changes of different aspects of the embodiments of the present application as described above. In order to be brief, they are not provided in detail.

[0122] In addition, in order to simplify the description and discussion, and so as not to make the embodiments of the present application difficult to understand, the well-known power / ground connections of integrated circuit (IC) chips and other components can or can not be shown in the provided drawings. In addition, the devices can be shown in the form of block diagrams in order to avoid making the embodiments of the present application difficult to understand, and this also takes into account the fact that the details of the implementation of these block diagram devices are highly dependent on the platform to be implemented the embodiments of the present application (i.e. these details should be fully within the understanding of those skilled in the art). Where specific details (e.g. circuits) are set forth in order to describe an exemplary embodiment of the present application, it will be apparent to those skilled in the art that the embodiments of the present application can be implemented without these specific details or with variations on these specific details. Therefore, these descriptions should be considered illustrative rather than limiting.

[0123] While the application has been described in terms of particular embodiments, many alternatives, modifications and variations will be apparent to those skilled in the art from the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) can use the embodiments discussed.

[0124] Embodiments of the application are intended to cover all such alternatives, modifications, and variations as come within the scope of the appended claims. Accordingly, any and all such alternatives, modifications, equivalents, improvements and the like are intended to be encompassed by the present application.

Claims

1. A method for predicting the boundary of a coal mining subsidence area considering the additional effect of wind-blown sand, characterized by, The method comprises: acquiring synthetic aperture radar images of a coal mining area; monitoring a first coal mining subsidence area boundary caused by coal mining according to the synthetic aperture radar images; calculating a sand drift moving boundary caused by sand drift migration according to a preset mapping model and a wind speed and a sand particle size of the coal mining area, wherein the mapping model is a model representing a mapping relationship between the wind speed, the sand particle size and the sand drift moving boundary; correcting the first coal mining subsidence area boundary by using the sand drift moving boundary to obtain a second coal mining subsidence area boundary; the establishment process of the mapping model comprises: simulating sand drift experiments under different wind speeds and different sand particle size combinations, fitting the sand drift moving boundaries in each sand drift experiment into corresponding mathematical models; extracting model parameters of the mathematical models; and establishing the mapping model according to a plurality of groups of wind speed, sand particle size and model parameters corresponding to each other obtained through the sand drift experiments; the mathematical model is a parabolic model, and parameters of the parabolic model comprise a constant term parameter, a first term coefficient and a second term coefficient; an upper boundary of the sand drift moving boundary is fitted into a first parabolic model, and a lower boundary of the sand drift moving boundary is fitted into a second parabolic model; the correction of the first coal mining subsidence area boundary by using the sand drift moving boundary specifically comprises: superimposing the sand drift moving boundary and the first coal mining subsidence area boundary, wherein a center of the sand drift moving boundary coincides with a center of the first coal mining subsidence area boundary after superimposition, and the sand drift moving boundary is located in a first plane, and the first plane is determined by a working face strike direction and a working face dip direction within the first coal mining subsidence area boundary.

2. The method of claim 1, wherein, The method further comprises: setting a measuring device according to a preset geometric similarity ratio and a preset time similarity ratio to simulate the sand drift experiments by using the measuring device.

3. The method of claim 1, wherein, The monitoring of the first coal mining subsidence area boundary caused by coal mining according to the synthetic aperture radar images specifically comprises: processing the synthetic aperture radar images to generate a time-series subsidence rate graph and a cumulative subsidence graph; calculating prediction parameters according to the time-series subsidence rate graph and the cumulative subsidence graph; predicting a surface subsidence amount, a surface horizontal movement amount, a surface inclination deformation amount, a surface curvature deformation amount and a surface horizontal deformation amount by using a probability integral method in combination with the prediction parameters; generating the first coal mining subsidence area boundary according to a prediction result.

4. The method according to claim 3, wherein, The surface horizontal movement amount is: ； wherein, , are the horizontal movement of the surface, axis component, axis component, is the maximum horizontal movement of the surface, is the first integration variable, is the second integration variable, is the main influence radius, is the subsidence of the surface, is the coal mining influence propagation angle.

5. A system for predicting the boundary of a coal mining subsidence area taking into account the additional effect of wind-blown sand, characterized in that, The method comprises: an acquisition module configured to acquire synthetic aperture radar images of a coal mining area; a monitoring module configured to monitor a first coal mining subsidence area boundary caused by coal mining according to the synthetic aperture radar images; a calculation module configured to calculate a sand drift moving boundary caused by sand drift migration according to a preset mapping model and a wind speed and a sand particle size of the coal mining area, wherein the mapping model is a model representing a mapping relationship between the wind speed, the sand particle size and the sand drift moving boundary; a correction module configured to correct the first coal mining subsidence area boundary by using the sand drift moving boundary to obtain a second coal mining subsidence area boundary; the system is configured to perform the method according to any one of claims 1-4.

6. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor implements the method of any one of claims 1-4 when executing the program.

7. A non-transitory computer-readable storage medium storing computer instructions, wherein, The computer instructions are for causing a computer to perform the method of any one of claims 1-4.

Images