[go: up one dir, main page]

CN116050657B - Surface elevation prediction method, device, terminal equipment and medium for closed mine - Google Patents

Surface elevation prediction method, device, terminal equipment and medium for closed mine

Info

Publication number
CN116050657B
CN116050657B CN202310171366.4A CN202310171366A CN116050657B CN 116050657 B CN116050657 B CN 116050657B CN 202310171366 A CN202310171366 A CN 202310171366A CN 116050657 B CN116050657 B CN 116050657B
Authority
CN
China
Prior art keywords
ground surface
mine
lifting
surface point
hossfeld
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202310171366.4A
Other languages
Chinese (zh)
Other versions
CN116050657A (en
Inventor
杨泽发
柴嘉遥
吴立新
李志伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central South University
Original Assignee
Central South University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Central South University filed Critical Central South University
Priority to CN202310171366.4A priority Critical patent/CN116050657B/en
Publication of CN116050657A publication Critical patent/CN116050657A/en
Application granted granted Critical
Publication of CN116050657B publication Critical patent/CN116050657B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B15/00—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
    • G01B15/06—Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons for measuring the deformation in a solid
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Human Resources & Organizations (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Electromagnetism (AREA)
  • Development Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Mathematical Physics (AREA)
  • Quality & Reliability (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Animal Husbandry (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Mining & Mineral Resources (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Agronomy & Crop Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)

Abstract

本申请适用于地表形变监测技术领域，提供了一种关闭矿井的地表抬升预测方法、装置、终端设备及介质，通过获取关闭矿井的雷达观测数据，并利用时序InSAR方法对雷达观测数据进行处理，得到关闭矿井的LOS向形变数据；利用Hossfeld函数模型对每个地表点位的地表抬升进行拟合，并建立该地表点位对应的LOS向形变数据与Hossfeld模型的关系表达式；利用遗传算法对关系表达式的参数进行反演，得到该地表点位对应的参数值；利用参数值对Hossfeld函数模型进行更新，得到该地表点位的地表抬升预测模型；利用地表抬升预测模型对该地表点位，在待测时刻的地表抬升进行预测。本申请能提高关闭矿井的地表抬升预测的准确性。

The present application is applicable to the field of surface deformation monitoring technology, and provides a surface uplift prediction method, apparatus, terminal equipment, and medium for closed mines. By acquiring radar observation data of the closed mine and processing the radar observation data using a time-series InSAR method, the LOS deformation data of the closed mine is obtained; the surface uplift of each surface point is fitted using a Hossfeld function model, and a relationship expression between the LOS deformation data corresponding to the surface point and the Hossfeld model is established; the parameters of the relationship expression are inverted using a genetic algorithm to obtain the parameter values corresponding to the surface point; the Hossfeld function model is updated using the parameter values to obtain a surface uplift prediction model for the surface point; and the surface uplift prediction model is used to predict the surface uplift of the surface point at the time to be measured. The present application can improve the accuracy of surface uplift prediction for closed mines.

Description

Surface elevation prediction method, device, terminal equipment and medium for closed mine
Technical Field
The application belongs to the technical field of surface deformation monitoring, and particularly relates to a surface lifting prediction method, device, terminal equipment and medium for a closed mine.
Background
Currently, conventional measurement methods such as leveling, global navigation satellite system (GNSS, global Navigation SATELLITE SYSTEM) and the like are generally used for monitoring the surface elevation of a closed mine, but the conventional measurement methods have some limitations such as high cost, limitation of measurement range, risk of field measurement in mining areas and the like. The synthetic aperture radar interferometry (InSAR, interferometric Synthetic Aperture Rader) has the advantages of wide coverage, all-weather all-day time, high resolution, low cost, historical backtracking and the like, and is an effective monitoring technology for observing the surface elevation of a closed mine.
In recent years, students at home and abroad try to monitor the surface time sequence lifting of a closed mine for a long period (years or even decades) by using satellite-borne InSAR data, and the evolution rule from residual settlement to slow rebound exists on the surface of the closed mine with most of water system supply enrichment. The learner observes continuous rebound of underground water after the mine is closed, and introduces InSAR technology to combine with a physical mechanical model to predict the surface elevation of the closed mine caused by the rebound of underground water. The existing prediction models for the surface lifting are mainly based on soil mechanics models such as an effective stress principle, and the models have the advantages that the surface lifting of the closed mine is predicted by introducing a physical mechanics model and combining parameters such as the underground water level, and the surface lifting situation predicted based on the physical mechanics rule is more fit with reality theoretically, but the prediction method is highly dependent on observation parameters such as the underground water level of the mine, and the parameters can be obtained only by on-site measurement, and most of the mine is not observed after being closed, so that the prediction of the surface lifting is inaccurate or can not be performed due to the lack of the parameters such as the underground water level required by the prediction model of the closed mine.
Disclosure of Invention
The embodiment of the application provides a method, a device, terminal equipment and a medium for predicting the surface elevation of a closed mine, which can solve the problem of inaccurate surface elevation prediction of the closed mine at present.
In a first aspect, an embodiment of the present application provides a method for predicting surface elevation of a closed mine, including:
Acquiring radar observation data of a target closed mine in a preset time period, and processing the radar observation data by using a time sequence InSAR method to obtain LOS deformation data of the target closed mine, wherein the LOS deformation data comprises earth surface elevation values of a plurality of earth surface points of the target closed mine;
for any one of a plurality of surface points, performing the steps of:
Fitting the surface elevation of the surface point by using Hossfeld function models, and establishing a relational expression of LOS deformation data corresponding to the surface point and the Hossfeld model;
Inverting the parameters of the relational expression by using a genetic algorithm to obtain a parameter value corresponding to the earth surface point location;
updating Hossfeld the function model by using the parameter value to obtain a ground surface elevation prediction model of the ground surface point location;
And predicting the earth surface lifting of the earth surface point position at the moment to be detected by using an earth surface lifting prediction model. Optionally, fitting the surface elevation of the surface point location using a Hossfeld function model includes:
By calculation formula
Obtaining a ground surface lifting fitting value D H (t) of the ground surface point, wherein a represents a ground surface lifting value influence parameter, b represents a ground surface lifting time influence parameter, c represents a ground surface lifting speed influence parameter, and t represents unit time.
Optionally, the established relational expression is as follows:
Wherein, the Indicating the unwrapping phase of the target shut-down mine at the surface point (x, y), indicating the wavelength of the radar,Represents the surface elevation value of the (x, y) point in the time period from t n-1 to t n, and n represents the maximum value of the preset time period.
Optionally, the expression of the surface lift prediction model is as follows:
Wherein t k represents the k time, k=1, 2,..n, D (t k) represents the predicted value of the surface elevation of the surface point at time t k, and all represent the parameter values corresponding to the surface point obtained by inversion of the genetic algorithm.
In a second aspect, an embodiment of the present application provides a surface lift prediction apparatus for closing a mine, including:
The system comprises a data acquisition module, a target closed mine and a target closed mine, wherein the data acquisition module is used for acquiring radar observation data of a preset time period of the target closed mine and processing the radar observation data by utilizing a time sequence InSAR method to obtain LOS deformation data of the target closed mine;
The relational expression construction module is used for fitting the surface elevation of the surface point by utilizing the Hossfeld function model and establishing a relational expression of LOS deformation data corresponding to the surface point and the Hossfeld model;
The parameter inversion module is used for inverting the parameters of the relational expression by utilizing a genetic algorithm to obtain a parameter value corresponding to the earth surface point position;
The prediction model construction module is used for updating Hossfeld function models by utilizing parameter values to obtain a ground surface elevation prediction model of the ground surface point position;
The ground surface lifting prediction module is used for predicting the ground surface lifting of the ground surface point position at the moment to be detected by using the ground surface lifting prediction model.
In a third aspect, an embodiment of the present application provides a terminal device, including a memory, a processor, and a computer program stored in the memory and capable of running on the processor, where the processor implements the above-mentioned method for predicting surface elevation of a closed mine when executing the computer program.
In a fourth aspect, an embodiment of the present application provides a computer readable storage medium, where a computer program is stored, where the computer program when executed by a processor implements the above-mentioned method for predicting surface elevation of a closed mine.
The scheme of the application has the following beneficial effects:
In some embodiments of the application, radar observation data of a target closed mine in a preset time period is obtained, the radar observation data is processed by using a time sequence InSAR method to obtain LOS deformation data of the target closed mine, then a Hossfeld function model is used for fitting the earth surface elevation of the earth surface point, a relational expression of the LOS deformation data corresponding to the earth surface point and a Hossfeld model is established, then a genetic algorithm is used for inverting parameters of the relational expression to obtain parameter values corresponding to the earth surface point, then a Hossfeld function model is updated by using the parameter values to obtain an earth surface elevation prediction model of the earth surface point, finally the earth surface elevation of the earth surface point is predicted by using the earth surface elevation prediction model, and the earth surface elevation at the moment to be measured is predicted. The parameters of the relational expression are inverted by utilizing a genetic algorithm, and the Hossfeld function model is updated by the inverted parameter values to obtain a ground surface lifting prediction model, so that the ground surface lifting can be predicted by only referring to InSAR data, the influence of complex terrain environmental factors is avoided, and the accuracy of the ground surface lifting prediction of a closed mine is improved.
Other advantageous effects of the present application will be described in detail in the detailed description section which follows.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments or the description of the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for predicting surface elevation of a closed mine according to an embodiment of the present application;
FIG. 2 is a schematic diagram showing a comparison of ground surface elevation prediction effects of a ground surface elevation prediction model according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a device for predicting the surface elevation of a closed mine according to an embodiment of the present application;
Fig. 4 is a schematic structural diagram of a terminal device according to an embodiment of the present application.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth such as the particular system architecture, techniques, etc., in order to provide a thorough understanding of the embodiments of the present application. It will be apparent, however, to one skilled in the art that the present application may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.
It should be understood that the terms "comprises" and/or "comprising," when used in this specification and the appended claims, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
It should also be understood that the term "and/or" as used in the present specification and the appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes such combinations.
As used in the present description and the appended claims, the term "if" may be interpreted as "when..once" or "in response to a determination" or "in response to detection" depending on the context. Similarly, the phrase "if a determination" or "if a [ described condition or event ] is detected" may be interpreted in the context of meaning "upon determination" or "in response to determination" or "upon detection of a [ described condition or event ]" or "in response to detection of a [ described condition or event ]".
Furthermore, the terms "first," "second," "third," and the like in the description of the present specification and in the appended claims, are used for distinguishing between descriptions and not necessarily for indicating or implying a relative importance.
Reference in the specification to "one embodiment" or "some embodiments" or the like means that a particular feature, structure, or characteristic described in connection with the embodiment is included in one or more embodiments of the application. Thus, appearances of the phrases "in one embodiment," "in some embodiments," "in other embodiments," and the like in the specification are not necessarily all referring to the same embodiment, but mean "one or more but not all embodiments" unless expressly specified otherwise. The terms "comprising," "including," "having," and variations thereof mean "including but not limited to," unless expressly specified otherwise.
Aiming at the problem of inaccurate surface elevation prediction of the closed mine at present, the application provides a surface elevation prediction method, a device, terminal equipment and medium of the closed mine, wherein radar observation data of a target closed mine in a preset time period are obtained, the radar observation data are processed by utilizing a time sequence InSAR method to obtain LOS directional deformation data of the target closed mine, then a Hossfeld function model is utilized to fit the surface elevation of the surface point, a relational expression of the LOS directional deformation data corresponding to the surface point and a Hossfeld model is established, then a genetic algorithm is utilized to invert parameters of the relational expression to obtain parameter values corresponding to the surface point, the Hossfeld function model is updated by utilizing the parameter values to obtain a surface elevation prediction model of the surface point, and finally the surface elevation of the surface point at the moment to be detected is predicted by utilizing the surface elevation prediction model. The parameters of the relational expression are inverted by utilizing a genetic algorithm, and the Hossfeld function model is updated by the inverted parameter values to obtain a ground surface lifting prediction model, so that the ground surface lifting can be predicted by only referring to InSAR data, the influence of complex terrain environmental factors is avoided, and the accuracy of the ground surface lifting prediction of a closed mine is improved.
As shown in fig. 1, the method for predicting the surface elevation of the closed mine mainly comprises the following steps:
And 11, acquiring radar observation data of a target closed mine in a preset time period, and processing the radar observation data by using a time sequence InSAR method to obtain LOS deformation data of the target closed mine.
The radar observation data includes Sentinel (an observation satellite designed and transmitted by the european aerospace agency, which has high resolution, long-term continuous consistent data acquisition capability and system observation scheme design capability), down-track observation data, SRTM data (Shuttle RadarTopography Mission, a data jointly measured by the united states aerospace agency and national survey and drawing agency of the national defense department). Wherein the SRTM data is used as an external reference DEM (Digital Elevation Model ), which is an entity floor model that represents floor elevation in the form of a set of ordered arrays of values.
The time sequence InSAR method comprises the steps of performing geometric registration by an external DEM, performing relevant registration based on amplitude and intensity information of images, performing image registration such as fine registration by using precise orbit data, generating a differential interference pair by using a time base line and a space base line formed by a main image and a sub image, calculating a coherence coefficient of a phase point in time by using a differential phase, a filtering phase and a terrain residual phase after differential interference, and then obtaining lifting information of the earth surface through phase unwrapping, geocoding and the like.
The LOS deformation data comprises surface elevation values of a plurality of surface points of the target closed mine.
The radar observation data is processed by using a time sequence InSAR method to obtain LOS deformation data of a target closed mine, so that limitation of regional topography on surface deformation prediction can be avoided, and the method is safer and simpler and more convenient without depending on geodetic data (GPS measurement and the like).
The ground surface lifting belongs to ground surface deformation, and particularly refers to ground surface deformation caused by mine groundwater rebound.
For any surface point position of a closed mine, the following steps are executed:
And step 12, fitting the surface elevation of the surface point by using a Hossfeld function model, and establishing a relational expression of LOS deformation data corresponding to the surface point and the Hossfeld model.
It should be noted that under the general condition, the trend of the surface elevation shows an S-shaped increase, and the Hossfeld function model is very excellent in describing the trend curve, so that the Hossfeld function model is utilized to fit the surface elevation, and meanwhile, the parameters of the Hossfeld function model are simpler, so that complicated parameter observation and calculation processes are omitted, and the simplicity of the method is reflected.
In the step, the specific process of fitting the surface elevation of the surface point by using Hossfeld function model is as follows:
specifically, by a calculation formula
Obtaining a ground surface lifting fitting value D H (t) of the ground surface point, wherein a represents a ground surface lifting value influence parameter, b represents a ground surface lifting time influence parameter, c represents a ground surface lifting speed influence parameter, and t represents unit time.
In the above formula, the actual values of the influencing parameters (b, c) are unknown, so that the actual values need to be solved, and in order to solve the influencing parameters, in some embodiments of the application, a relational expression of the LOS-direction deformation data and the Hossfeld model is constructed.
In the following step, a specific process of establishing a relational expression between LOS-direction deformation data corresponding to the earth surface point location and Hossfeld model is illustrated as an example:
Specifically, the established relational expression is as follows:
Wherein, the Indicating the unwrapping phase of the target shut-down mine at the surface point (x, y), indicating the wavelength of the radar,Represents the surface elevation value of the (x, y) point in the time period from t n-1 to t n, and n represents the maximum value of the preset time period.
And 13, inverting the parameters of the relational expression by utilizing a genetic algorithm to obtain the parameter value corresponding to the earth surface point location.
The genetic algorithm can automatically acquire and accumulate knowledge about the search space in the data search process, and adaptively control the search process to obtain an optimal solution, so that the genetic algorithm is adopted to invert the influence parameter parameters in the step 12 until the model is fitted, and a group of optimal influence parameter values corresponding to the surface point position can be obtained.
Illustratively, in some embodiments of the present application, the maximum population number of the genetic algorithm may be set to 20000, the maximum number of iterations of the population to 5000, and the population crossover probability to 0.7.
It should be noted that the process of solving the optimal solution by using the genetic algorithm belongs to common general knowledge, and detailed solving process thereof is not described herein.
And 14, updating Hossfeld the function model by using the parameter values to obtain a ground surface lifting prediction model of the ground surface point.
Specifically, the expression of the constructed ground surface elevation prediction model is as follows:
Wherein t k represents the k time, k=1, 2,..n, D (t k) represents the predicted value of the surface elevation of the surface point at time t k, and all represent the parameter values corresponding to the surface point obtained by inversion of the genetic algorithm. And 15, predicting the earth surface lifting at the moment to be detected by using an earth surface lifting prediction model to the earth surface point position.
Illustratively, in some embodiments of the present application, at a surface point of a closed-loop mine, a set of impact parameter values for the surface point obtained by inversion of a genetic algorithm are a=0.091, b=216000, and c=3.61, respectively. The corresponding surface deformation prediction model for the surface point is as follows:
When the ground surface elevation value of the ground surface point position at a certain moment (t k) needs to be predicted, substituting the corresponding t k into the above formula to obtain the ground surface elevation predicted value at the moment.
The method comprises the steps of obtaining radar observation data of a target closed mine in a preset time period, processing the radar observation data by utilizing a time sequence InSAR method to obtain LOS deformation data of the target closed mine, fitting the earth surface elevation of the earth surface point by utilizing a Hossfeld function model, establishing a relational expression of the LOS deformation data corresponding to the earth surface point and a Hossfeld model, inverting parameters of the relational expression by utilizing a genetic algorithm to obtain parameter values corresponding to the earth surface point, updating a Hossfeld function model by utilizing the parameter values to obtain an earth surface elevation prediction model of the earth surface point, and predicting the earth surface elevation of the earth surface point at a moment to be detected by utilizing the earth surface elevation prediction model. The parameters of the relational expression are inverted by utilizing a genetic algorithm, and the Hossfeld function model is updated by the inverted parameter values to obtain a ground surface lifting prediction model, so that the ground surface lifting can be predicted by only referring to InSAR data, the influence of complex terrain environmental factors is avoided, and the accuracy of the ground surface lifting prediction of a closed mine is improved.
In order to better understand the method for predicting the surface elevation of the closed mine, which is provided by the application, an example is described below with reference to a specific embodiment.
In the embodiment of the application, 133 scene Sentinel derated observation data including target closed mines of 2016, 10 and 2021, 12 are selected, and the acquired SRTM data with the spatial resolution of 90 meters is used as external reference DEM.
And processing the acquired observation data by using a time sequence InSAR technology to obtain the LOS directional deformation result of the region. Data analysis shows that the accumulated deformation of the earth surface of the research area reaches 0.15m in the observation stage, namely the earth surface of the research area is lifted, and the phenomenon is caused by groundwater rebound through reference literature analysis.
After the InSAR observation value (LOS deformation) of the target closed mine is calculated, any earth surface point location is selected according to the growth rule (generally presenting S-shaped growth) presented by earth surface deformation, and a Hossfeld model (capable of better reflecting an S-shaped curve) is used as an earth surface deformation fitting model to obtain an earth surface lifting fitting model of the earth surface point location as follows:
Wherein D H (t) represents a surface elevation fitting value of the surface point, a represents a surface elevation value influence parameter, b represents a surface elevation time influence parameter, c represents a surface elevation speed influence parameter, and t represents unit time.
And then establishing a relational expression of LOS deformation data corresponding to the surface point location and Hossfeld model, which is specifically as follows:
Wherein, the Indicating the unwrapping phase of the target shut-down mine at the surface point (x, y), indicating the wavelength of the radar,Represents the surface elevation value of the (x, y) point in the time period from t n-1 to t n, and n represents the maximum value of the preset time period.
If the deformation condition of the earth surface single point is to be accurately described, accurate estimation is required to be performed on the parameters in the relational expression, and in consideration of the fact that the genetic algorithm can automatically acquire and accumulate knowledge about the search space in the data search process and adaptively control the search process to obtain the optimal solution, in the embodiment of the application, the genetic algorithm is adopted to perform inversion updating on the parameters (a, b and c) in the relational expression.
Through the steps, a group of parameter values corresponding to the surface point of the target closed mine are obtained, wherein a=0.091, b=216000 and c=3.61.
Correspondingly, the expression of the ground surface elevation prediction model of the ground surface point is as follows:
The ground surface elevation is predicted based on the above, and the ground surface elevation prediction is compared with the result shown in fig. 2, wherein the abscissa of fig. 2 represents time (year and month) and the ordinate represents elevation value (millimeter).
As can be seen from fig. 2, the lifting curve (dotted line) of the ground lifting prediction model is very consistent with the lifting trend of the LOS to the deformation curve (solid line) observed by the InSAR, which indicates that the ground lifting prediction model can well describe the LOS to dynamic displacement of the ground lifting. By calculation, the fitted RMSE (Root Mean Squared Error, root mean square error) is 0.0018m respectively, which is about 1.2% of the maximum LOS deformation value, so that the dynamic change of the surface elevation of the closed mine can be predicted by using the surface elevation prediction model.
The following illustrates an exemplary device for predicting the surface elevation of a closed mine.
As shown in fig. 3, an embodiment of the present application provides a surface lift prediction apparatus for a closed mine, the surface lift prediction apparatus 300 for a closed mine including:
the data acquisition module 301 is configured to acquire radar observation data of a target closed mine in a preset period of time, and process the radar observation data by using a time sequence InSAR method to obtain LOS-direction deformation data of the target closed mine, where the LOS-direction deformation data includes surface elevation values of a plurality of surface points of the target closed mine.
The relational expression construction module 302 is configured to fit the surface elevation of the surface point location by using a Hossfeld function model, and establish a relational expression of LOS-direction deformation data corresponding to the surface point location and the Hossfeld model.
And the parameter inversion module 303 is configured to invert the parameters of the relational expression by using a genetic algorithm to obtain a parameter value corresponding to the earth surface point location.
The prediction model construction module 304 is configured to update Hossfeld the function model with the parameter value, and obtain a surface elevation prediction model of the surface point.
The ground surface elevation prediction module 305 is configured to predict the ground surface elevation of the ground surface point location at the time to be measured by using a ground surface elevation prediction model.
It should be noted that, because the content of information interaction and execution process between the above devices/units is based on the same concept as the method embodiment of the present application, specific functions and technical effects thereof may be referred to in the method embodiment section, and will not be described herein.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions. The functional units and modules in the embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit, where the integrated units may be implemented in a form of hardware or a form of a software functional unit. In addition, the specific names of the functional units and modules are only for distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working process of the units and modules in the above system may refer to the corresponding process in the foregoing method embodiment, which is not described herein again.
As shown in fig. 4, an embodiment of the present application provides a terminal device, as shown in fig. 4, the terminal device D10 of this embodiment includes at least one processor D100 (only one processor is shown in fig. 4), a memory D101, and a computer program D102 stored in the memory D101 and executable on the at least one processor D100, the processor D100 implementing the steps in any of the respective method embodiments described above when executing the computer program D102.
Specifically, when the processor D100 executes the computer program D102, radar observation data of a target closed mine in a preset time period is obtained, the radar observation data is processed by using a time sequence InSAR method to obtain LOS deformation data of the target closed mine, then a Hossfeld function model is used for fitting the earth surface elevation of the earth surface point, a relational expression of LOS deformation data corresponding to the earth surface point and a Hossfeld model is established, then a genetic algorithm is used for inverting parameters of the relational expression to obtain parameter values corresponding to the earth surface point, a Hossfeld function model is updated by using the parameter values to obtain an earth surface elevation prediction model of the earth surface point, and finally the earth surface elevation of the earth surface point at the moment to be measured is predicted by using the earth surface elevation prediction model. The parameters of the relational expression are inverted by utilizing a genetic algorithm, and the Hossfeld function model is updated by the inverted parameter values to obtain a ground surface lifting prediction model, so that the ground surface lifting can be predicted by only referring to InSAR data, the influence of complex terrain environmental factors is avoided, and the accuracy of the ground surface lifting prediction of a closed mine is improved.
The Processor D100 may be a central processing Unit (CPU, centralProcessing Unit), and the Processor D100 may also be other general purpose processors, digital signal processors (DSP, digital Signal Processor), application SPECIFIC INTEGRATED Circuits (ASIC), off-the-shelf programmable gate arrays (FPGA, field-Programmable GateArray) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory D101 may in some embodiments be an internal storage unit of the terminal device D10, for example a hard disk or a memory of the terminal device D10. The memory D101 may also be an external storage device of the terminal device D10 in other embodiments, for example, a plug-in hard disk, a smart memory card (SMC, smartMedia Card), a Secure Digital (SD) card, a flash memory card (FLASH CARD) or the like, which are provided on the terminal device D10. Further, the memory D101 may also include both an internal storage unit and an external storage device of the terminal device D10. The memory D101 is used for storing an operating system, an application program, a boot loader (BootLoader), data, other programs, etc., such as program codes of the computer program. The memory D101 may also be used to temporarily store data that has been output or is to be output.
Embodiments of the present application also provide a computer readable storage medium storing a computer program which, when executed by a processor, implements steps for implementing the various method embodiments described above.
Embodiments of the present application provide a computer program product enabling a terminal device to carry out the steps of the method embodiments described above when the computer program product is run on the terminal device.
The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the present application may implement all or part of the flow of the method of the above embodiments, and may be implemented by a computer program to instruct related hardware, where the computer program may be stored in a computer readable storage medium, and when the computer program is executed by a processor, the computer program may implement the steps of each of the method embodiments described above. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium can include at least any entity or device capable of carrying computer program code to a surface lift prediction device/terminal equipment that shuts down a mine, a recording medium, a computer Memory, a Read-only Memory (ROM), a random access Memory (RAM, randomAccess Memory), an electrical carrier signal, a telecommunications signal, and a software distribution medium. Such as a U-disk, removable hard disk, magnetic or optical disk, etc. In some jurisdictions, computer readable media may not be electrical carrier signals and telecommunications signals in accordance with legislation and patent practice.
In the foregoing embodiments, the descriptions of the embodiments are emphasized, and in part, not described or illustrated in any particular embodiment, reference is made to the related descriptions of other embodiments.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus/network device and method may be implemented in other manners. For example, the apparatus/network device embodiments described above are merely illustrative, e.g., the division of the modules or units is merely a logical functional division, and there may be additional divisions in actual implementation, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection via interfaces, devices or units, which may be in electrical, mechanical or other forms.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
The technical scheme provided by the application has the advantages that the dynamic change process of the surface lifting of the closed mine along with time can be well described and reasonably predicted, and the method is important for the heavy demands of closed mine reuse, mine geological disaster prevention and control and the like.
While the foregoing is directed to the preferred embodiments of the present application, it will be appreciated by those skilled in the art that various modifications and adaptations can be made without departing from the principles of the present application, and such modifications and adaptations are intended to be comprehended within the scope of the present application.

Claims (5)

1. A method for predicting surface elevation of a closed mine, comprising:
Acquiring radar observation data of a target closed mine in a preset time period, and processing the radar observation data by using a time sequence InSAR method to obtain LOS deformation data of the target closed mine, wherein the LOS deformation data comprises earth surface elevation values of a plurality of earth surface points of the target closed mine;
For any one of the plurality of surface points, performing the steps of:
fitting the surface elevation of the surface point location by using Hossfeld function models, wherein the fitting comprises the following steps:
By calculation formula
Obtaining a ground surface lifting fitting value D H (t) of the ground surface point, wherein a represents a ground surface lifting value influence parameter, b represents a ground surface lifting time influence parameter, c represents a ground surface lifting speed influence parameter, t represents unit time, and establishing a relational expression of LOS deformation data corresponding to the ground surface point and the Hossfeld function model as follows:
Wherein, the Indicating the unwrapping phase of the target shut-down mine at the surface point (x, y), lambda indicating the wavelength of the radar,Representing the ground surface elevation value of the (x, y) point in the time period from t n-1 to t n, and n represents the maximum value of the preset time period;
Inverting the parameters of the relational expression by utilizing a genetic algorithm to obtain a parameter value corresponding to the earth surface point location;
updating the Hossfeld function model by using the parameter value to obtain a ground surface lifting prediction model of the ground surface point position;
And predicting the earth surface lifting of the earth surface point position at the moment to be detected by using the earth surface lifting prediction model.
2. The prediction method according to claim 1, wherein the expression of the surface elevation prediction model is as follows:
Wherein t k represents a k time, k=1, 2,..n, D (t k) represents a predicted value of the surface elevation of the surface point at time t k, and a ', b ', c ' each represent a parameter value corresponding to the surface point obtained by inversion of a genetic algorithm.
3. A surface lift prediction apparatus for a closed mine, comprising:
The system comprises a data acquisition module, a target closed mine and a target acquisition module, wherein the data acquisition module is used for acquiring radar observation data of a target closed mine in a preset time period, and processing the radar observation data by utilizing a time sequence InSAR method to obtain LOS deformation data of the target closed mine;
the relational expression construction module is used for fitting the surface elevation of the surface point by using Hossfeld function models, and comprises the following steps:
By calculation formula
Obtaining a ground surface lifting fitting value D H (t) of the ground surface point, wherein a represents a ground surface lifting value influence parameter, b represents a ground surface lifting time influence parameter, c represents a ground surface lifting speed influence parameter, t represents unit time, and establishing a relational expression of LOS deformation data corresponding to the ground surface point and the Hossfeld function model as follows:
Wherein, the Indicating the unwrapping phase of the target shut-down mine at the surface point (x, y), lambda indicating the wavelength of the radar,Representing the ground surface elevation value of the (x, y) point in the time period from t n-1 to t n, and n represents the maximum value of the preset time period;
The parameter inversion module is used for inverting the parameters of the relational expression by utilizing a genetic algorithm to obtain a parameter value corresponding to the earth surface point position;
the prediction model construction module is used for updating the Hossfeld function model by utilizing the parameter value to obtain a ground surface lifting prediction model of the ground surface point position;
The ground surface lifting prediction module is used for predicting the ground surface lifting of the ground surface point position at the moment to be detected by using the ground surface lifting prediction model.
4. A terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements a method of predicting surface lift to shut down a mine as claimed in any one of claims 1 to 2.
5. A computer readable storage medium storing a computer program, wherein the computer program when executed by a processor implements the method of surface lift prediction for closing a mine as claimed in any one of claims 1 to 2.
CN202310171366.4A 2023-02-28 2023-02-28 Surface elevation prediction method, device, terminal equipment and medium for closed mine Active CN116050657B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310171366.4A CN116050657B (en) 2023-02-28 2023-02-28 Surface elevation prediction method, device, terminal equipment and medium for closed mine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310171366.4A CN116050657B (en) 2023-02-28 2023-02-28 Surface elevation prediction method, device, terminal equipment and medium for closed mine

Publications (2)

Publication Number Publication Date
CN116050657A CN116050657A (en) 2023-05-02
CN116050657B true CN116050657B (en) 2025-09-09

Family

ID=86118343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310171366.4A Active CN116050657B (en) 2023-02-28 2023-02-28 Surface elevation prediction method, device, terminal equipment and medium for closed mine

Country Status (1)

Country Link
CN (1) CN116050657B (en)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110058236B (en) * 2019-05-21 2023-04-07 中南大学 InSAR and GNSS weighting method oriented to three-dimensional surface deformation estimation
CN110991048B (en) * 2019-12-04 2021-08-27 中国矿业大学 Prediction method for surface subsidence of closed well industrial and mining
CN111323776B (en) * 2020-02-27 2021-04-13 长沙理工大学 Method for monitoring deformation of mining area
CN113253270A (en) * 2021-06-11 2021-08-13 中国测绘科学研究院 Method and system for inverting underground mining parameters based on InSAR and Okada models

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于InSAR的废弃煤矿地表抬升建模与预测;柴嘉遥;中国优秀硕士学位论文全文数据库;20250215;全文 *

Also Published As

Publication number Publication date
CN116050657A (en) 2023-05-02

Similar Documents

Publication Publication Date Title
Kalenjuk et al. Processing of mobile laser scanning data for large‐scale deformation monitoring of anchored retaining structures along highways
Matwij et al. Determination of underground mining-induced displacement field using multi-temporal TLS point cloud registration
Tapete et al. Localising deformation along the elevation of linear structures: An experiment with space-borne InSAR and RTK GPS on the Roman Aqueducts in Rome, Italy
CN117648873B (en) Land subsidence prediction method, training method, device, equipment
Shi et al. Early soil consolidation from magnetic extensometers and full resolution SAR interferometry over highly decorrelated reclaimed lands
Wang et al. Interferometric synthetic aperture radar statistical inference in deformation measurement and geophysical inversion: A review
CN118259280B (en) Method, system and terminal for deformation assessment of airport in reclamation area by combining InSAR and GNSS
Wang et al. Analysis and prediction of regional land subsidence with InSAR technology and machine learning algorithm
Armaş et al. InSAR validation based on GNSS measurements in Bucharest
Huang et al. An Improved Adaptive Template Size Pixel‐Tracking Method for Monitoring Large‐Gradient Mining Subsidence
CN118566920B (en) Method, equipment and medium for winding same-vibration interference pattern
Yao et al. Research on Surface Deformation of Ordos Coal Mining Area by Integrating Multitemporal D‐InSAR and Offset Tracking Technology
CN114114257B (en) A method and device for detecting the correlation between dam area deformation and water level
Brighenti et al. UAV survey method to monitor and analyze geological hazards: the case study of the mud volcano of Villaggio Santa Barbara, Caltanissetta (Sicily)
Ren et al. Monitoring Yanwan deep-seated toppling deformation with the impact of water-level fluctuation by SAR observations
Wang et al. Monitoring, Analyzing, and Modeling for Single Subsidence Basin in Coal Mining Areas Based on SAR Interferometry with L‐Band Data
CN116299466B (en) Geological deformation monitoring method and device for power transmission channel
Tan et al. Deformation Monitoring and Spatiotemporal Evolution of Mining Area with Unmanned Aerial Vehicle and D‐InSAR Technology
Shahbazi et al. Detection of buildings with potential damage using differential deformation maps
Schmid et al. Georeferencing of terrestrial radar images in geomonitoring using kernel correlation
CN116050657B (en) Surface elevation prediction method, device, terminal equipment and medium for closed mine
Manconi et al. Surface displacements following the Mw 6.3 L’Aquila earthquake: One year of continuous monitoring via Robotized Total Station
Becek et al. Identifying land subsidence using global digital elevation models
Beshr et al. Using modified inverse distance weight and principal component analysis for spatial interpolation of foundation settlement based on geodetic observations
Hassan et al. Integration of GNSS observations with volunteered geographic information for improved navigation performance

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant