CN116952119A - Deformation monitoring model building method and device, electronic equipment and storage medium - Google Patents
Deformation monitoring model building method and device, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN116952119A CN116952119A CN202310919294.7A CN202310919294A CN116952119A CN 116952119 A CN116952119 A CN 116952119A CN 202310919294 A CN202310919294 A CN 202310919294A CN 116952119 A CN116952119 A CN 116952119A
- Authority
- CN
- China
- Prior art keywords
- deformation monitoring
- monitoring model
- model
- boundary
- deformation
- 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.)
- Pending
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 194
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000000463 material Substances 0.000 claims abstract description 47
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 23
- 239000011707 mineral Substances 0.000 claims abstract description 23
- 238000006073 displacement reaction Methods 0.000 claims description 54
- 229920001971 elastomer Polymers 0.000 claims description 32
- 239000000806 elastomer Substances 0.000 claims description 32
- 238000004590 computer program Methods 0.000 claims description 22
- 239000003245 coal Substances 0.000 claims description 18
- 238000005457 optimization Methods 0.000 claims description 18
- 238000010206 sensitivity analysis Methods 0.000 claims description 18
- 239000002184 metal Substances 0.000 claims description 5
- 229910052751 metal Inorganic materials 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000005065 mining Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000007405 data analysis Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
-
- 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/904—SAR modes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application provides a method and a device for establishing a deformation monitoring model, electronic equipment and a storage medium. The method comprises the following steps: determining the type of a deformation monitoring model corresponding to the current region to be detected according to mineral resources in the current region to be detected; establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model; acquiring material parameters of a current region to be detected at a current time, inputting the material parameters into a deformation monitoring model, and outputting a deformation monitoring result based on a current boundary condition and an initial condition; and obtaining deformation field data of the current region to be detected at the current time, and optimizing the deformation monitoring model according to the deformation field data and the deformation monitoring result to obtain an optimized deformation monitoring model. The method and the device can effectively ensure the accuracy of the deformation monitoring result, and are not limited by the professional level of the user.
Description
Technical Field
The present application relates to the field of earth surface deformation monitoring technologies, and in particular, to a method and an apparatus for establishing a deformation monitoring model, an electronic device, and a storage medium.
Background
Along with the high-speed development of the economy in China, the demands for mineral resources such as coal, iron and the like are also growing year by year. In recent years, the problem of surface deformation is made worse by mining disturbances of mineral resources. Therefore, it is necessary to monitor the surface deformation of the mining area and its surroundings in time.
The synthetic aperture radar interference (Corner Reflectors-Interferometric Synthetic Aperture Radar, CR-InSAR) technology of the corner reflector can measure the deformation of the micrometer or even nanometer level, and is widely applied to the field of monitoring the surface deformation. When the CR-InSAR technology is adopted to monitor the surface deformation, a large amount of data processing and analysis are required, and in addition, the radar signal is considered to be influenced by atmospheric refraction, and real-time data correction is required during data processing. This makes the CR-InSAR technique necessary under the direction of a highly skilled technician to generate accurate measurements.
In summary, when the CR-InSAR technology is adopted to monitor the surface deformation, the surface deformation is limited by the professional level of a user, and the accuracy of a measurement result cannot be effectively ensured.
Disclosure of Invention
The embodiment of the application provides a method, a device, electronic equipment and a storage medium for establishing a deformation monitoring model, which are used for solving the problems that the surface deformation monitoring is limited by the professional level of a user and the accuracy of a measurement result cannot be effectively ensured when a CR-InSAR technology is adopted.
In a first aspect, an embodiment of the present application provides a method for establishing a deformation monitoring model, including:
determining the type of a deformation monitoring model corresponding to a current region to be detected according to mineral resources in the current region to be detected;
establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model;
acquiring material parameters of a current region to be detected at a current time, inputting the material parameters into the deformation monitoring model, and outputting a deformation monitoring result based on a current boundary condition and an initial condition;
and obtaining deformation field data of the current region to be detected at the current time, and optimizing the deformation monitoring model according to the deformation field data and the deformation monitoring result to obtain an optimized deformation monitoring model.
In a possible implementation manner, the determining, according to mineral resources in a current area to be measured, a type of a deformation monitoring model corresponding to the current area to be measured includes:
when the current region to be detected is a coal mine, determining a deformation monitoring model corresponding to the current region to be detected as an elastomer model;
when the current area to be measured is nonferrous metal ore, determining the deformation monitoring model corresponding to the current area to be measured as a porous medium model.
In one possible implementation, when the deformation monitoring model is an elastomer model, the boundary condition includes: force boundary conditions, displacement boundary conditions, equilibrium boundary conditions, and material property boundary conditions;
and establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model, wherein the boundary conditions and initial conditions comprise:
according toEstablishing a force boundary condition of the elastomer model;
wherein sigma n Represents boundary normal stress, u represents displacement, n represents normal vector, f n Representing the force exerted by the boundary;
according toEstablishing displacement boundary conditions of the elastomer model; wherein (1)>Represents the boundary of the current region to be measured omega, u 0 Representing the displacement of the boundary fixation;
according toEstablishing equilibrium boundary conditions of the elastomer model; wherein ε represents strain and σ represents stress;
according toEstablishing a material property boundary condition of the elastomer model; where g represents the material properties at the boundary.
In one possible implementation manner, the establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model further includes:
according to
Establishing initial conditions of the elastomer model;
where u (x, y, z, t=0) represents the time t=0, the displacement at the position with coordinates (x, y, z), u_0 (x, y, z) represents the initial displacement at the position, v_0 (x, y, z) represents the initial velocity at the position, a represents the first amplitude of vibration, lx represents the length of the coal mine, ly represents the width of the coal mine, lz represents the height of the coal mine, ω represents the circular frequency of vibration, B represents the second amplitude of vibration, σ (x, y, z, t=0) represents the time t=0, the stress at the position with coordinates (x, y, z), σ_0 (x, y, z) represents the initial stress at the position, C represents the third amplitude of vibration, and D represents the fourth amplitude of vibration.
In one possible implementation, when the deformation monitoring model is a porous medium model, the boundary condition includes: displacement boundary conditions, force boundary conditions, and material property boundary conditions;
and establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model, wherein the boundary conditions and initial conditions comprise:
according toEstablishing a displacement boundary condition of the porous medium model;
wherein u represents displacement, n represents normal vector, and g_n represents normal displacement boundary condition;
establishing a force boundary condition of the porous medium model according to σn=f_n;
wherein sigma represents stress, n represents normal vector, and f_n represents normal stress boundary condition;
establishing a material property boundary condition of the porous medium model according to epsilon=epsilon sigma;
where ε represents strain.
In one possible implementation manner, the establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model further includes:
establishing an initial displacement condition of the porous medium model according to u (x, y, z, t=0) =0;
where u (x, y, z, t=0) represents the displacement at the position of the coordinates (x, y, z) at time t=0;
according to sigma ij (x,y,z,t＝0)＝σ 0 δ ij Establishing an initial stress condition of the porous medium model;
wherein sigma ij (x, y, z, t=0) represents the stress tensor at the position of coordinates (x, y, z), σ, at time t=0 0 Represents an initial stress constant, δ when i=j ij =1, otherwise, δ ij ＝0。
In one possible implementation manner, according to the deformation field data and the deformation monitoring result, optimizing the deformation monitoring model to obtain an optimized deformation monitoring model, including:
performing sensitivity analysis on boundary conditions of the deformation monitoring model to obtain sensitivity analysis results of all the boundary conditions;
determining an output error of the deformation monitoring model according to the deformation field data and the deformation monitoring result;
and when the output error is larger than an error threshold, optimizing and adjusting each boundary condition according to the sensitivity analysis result of each boundary condition and the output error until the output error of the deformation monitoring model after optimizing and adjusting is smaller than the error threshold, and obtaining the deformation monitoring model after optimizing and finishing.
In a second aspect, an embodiment of the present application provides a device for establishing a deformation monitoring model, including:
the modeling module is used for determining the type of the deformation monitoring model corresponding to the current region to be detected according to mineral resources in the current region to be detected;
the modeling module is further used for establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model;
the optimization module is used for acquiring material parameters of the current region to be detected at the current time, inputting the material parameters into the deformation monitoring model, and outputting a deformation monitoring result based on the current boundary condition and the initial condition;
the optimization module is used for acquiring deformation field data of the current region to be detected at the current time, and optimizing the deformation monitoring model according to the deformation field data and the deformation monitoring result to obtain an optimized deformation monitoring model.
In a third aspect, an embodiment 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, the processor implementing the steps of the method according to the first aspect or any one of the possible implementations of the first aspect, when the computer program is executed by the processor.
In a fourth aspect, embodiments of the present application provide a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the method as described above in the first aspect or any one of the possible implementations of the first aspect.
The embodiment of the application provides a method, a device, electronic equipment and a storage medium for establishing a deformation monitoring model, which are used for carrying out surface deformation monitoring on a current region to be detected, and a deformation monitoring result can be obtained by inputting material parameters at the current moment into the optimized deformation monitoring model without carrying out a large amount of data analysis and data correction work. According to the embodiment of the application, the deformation monitoring model is optimized and adjusted by adopting actually measured deformation field data in advance, so that the optimized and adjusted deformation monitoring model is not limited by the professional level of a user, and a deformation monitoring result with higher accuracy is correspondingly output in real time according to the input material parameters. In addition, considering that the mining activities of different mineral resources have different influence degrees on the surface deformation, the embodiment of the application correspondingly establishes different types of deformation monitoring models aiming at different mineral resources, thereby being convenient for monitoring the surface deformation caused by different mineral resources in a targeted manner and further improving the accuracy of the deformation monitoring result.
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 described below, it being obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flowchart of an implementation of a method for establishing a deformation monitoring model according to an embodiment of the present application;
fig. 2 is a schematic structural diagram of a device for establishing a deformation monitoring model according to an embodiment of the present application;
fig. 3 is a schematic diagram of an electronic 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.
For the purpose of making the objects, technical solutions and advantages of the present application more apparent, the following description will be made by way of specific embodiments with reference to the accompanying drawings.
Fig. 1 is a flowchart of an implementation of a method for establishing a deformation monitoring model according to an embodiment of the present application, which is described in detail below:
and step 101, determining the type of the deformation monitoring model corresponding to the current region to be detected according to mineral resources in the current region to be detected.
Surface deformation problems are typically caused by disturbances in the exploitation of mineral resources. The influence degree of different mineral resources on the surface deformation is considered to be different. In the embodiment of the application, different types of deformation monitoring models are correspondingly established for different mineral resources.
In some embodiments, step 101 may comprise:
when the current region to be detected is a coal mine, determining a deformation monitoring model corresponding to the current region to be detected as an elastomer model;
when the current area to be measured is nonferrous metal ore, determining the deformation monitoring model corresponding to the current area to be measured as a porous medium model.
In an elastomer model, the earth is considered to be composed of a plurality of elastic layers, each having different physical properties and densities. The elastomer model may simulate the propagation of seismic waves in different geological media through elastic wave propagation. Thus, the elastomer model is more suitable for studying the elastic properties of the earth's internal mass.
In a porous media model, the subsurface medium is considered to be composed of a plurality of pores and solids, and the process of transporting fluid movement in the porous medium can be described by Darcy's law. Therefore, the porous medium model is more suitable for researching problems of underground hydrology, oil and gas reservoirs and the like.
Based on the characteristics of the elastomer model and the porous medium model, when the area to be measured is a coal mine, the elastomer model is correspondingly established, and when the area to be measured is nonferrous metal ore, such as iron ore, the porous medium model is correspondingly established. When building an elastomer model or a porous medium model, this can be achieved by SciPr library.
And 102, establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model.
In some embodiments, when the deformation monitoring model is an elastomeric model, the boundary conditions include: force boundary conditions, displacement boundary conditions, equilibrium boundary conditions, and material property boundary conditions.
When establishing force boundary conditions, it is often specified that certain parts of the object are subjected to external forces. In particular according toForce boundary conditions of the elastomeric model are established.
Wherein sigma n Represents boundary normal stress, u represents displacement, n represents normal vector, f n Representing the force exerted by the boundary.
When establishing displacement boundary conditions, variations or deformations of certain parts of the object are typically specified. In particular according toEstablishing displacement edges of elastomeric modelsBoundary conditions.
Wherein,,represents the boundary of the current region to be measured omega, u 0 Representing a displacement of the boundary fixation.
In establishing the equilibrium boundary condition, the sum of the forces and moments to which the surface of the object is subjected is typically specified to be 0. In particular according toEstablishing balance boundary conditions of the elastomer model; where ε represents strain and σ represents stress.
The material properties of an object are typically specified when establishing a material property boundary. In particular according toEstablishing a material property boundary condition of the elastomer model; where g represents the material properties at the boundary.
The initial condition refers to the state of the elastomer at time t=0. The method specifically comprises the following steps: initial displacement, initial velocity, and initial stress.
In some embodiments, it may be according to
Establishing initial conditions of an elastomer model;
where u (x, y, z, t=0) represents the time t=0, the displacement at the position of coordinates (x, y, z), u_0 (x, y, z) represents the initial displacement at the position, v_0 (x, y, z) represents the initial velocity at the position, a represents the first amplitude of vibration, lx represents the length of the coal mine, ly represents the width of the coal mine, lz represents the height of the coal mine, ω represents the circular frequency of vibration, B represents the second amplitude of vibration, σ (x, y, z, t=0) represents the time t=0, the stress at the position of coordinates (x, y, z), σ_0 (x, y, z) represents the initial stress at the position, C represents the third amplitude of vibration, and D represents the fourth amplitude of vibration.
In some embodiments, when the deformation monitoring model is a porous media model, the boundary conditions include: displacement boundary conditions, force boundary conditions, and material property boundary conditions.
The displacement boundary conditions of the porous media model generally refer to displacement values or derivative values of displacement specified at certain boundaries (e.g., surfaces). In particular according toEstablishing displacement boundary conditions of the porous medium model;
where u represents displacement, n represents normal vector, and g_n represents normal displacement boundary condition.
The force boundary conditions of a porous media model generally refer to values of force or values of force normal guides specified on certain boundaries. Specifically, a force boundary condition of the porous medium model can be established according to σn=f_n;
where σ represents stress, n represents normal vector, and f_n represents normal stress boundary condition.
The material property boundary conditions of a porous media model generally refer to the values or rules of variation of material parameters specified at certain boundaries. Specifically, a material property boundary condition of the porous medium model can be established according to epsilon=epsilon sigma;
where ε represents strain.
The initial conditions include: initial displacement conditions and initial stress conditions. In some embodiments, the initial displacement condition of the porous media model may be established from u (x, y, z, t=0) =0;
where u (x, y, z, t=0) represents the displacement at the time t=0 at the position of coordinates (x, y, z).
In some embodiments, one may rely on σ ij (x,y,z,t＝0)＝σ 0 δ ij Establishing an initial stress condition of a porous medium model;
wherein sigma ij (x, y, z, t=0) represents the stress tensor at the position of coordinates (x, y, z), σ, at time t=0 0 Represents an initial stress constant, δ when i=j ij =1, otherwise, δ ij ＝0。
According to the initial stress conditions, the initial stress in each direction in the porous medium model is the same and does not change with the position.
And step 103, acquiring material parameters of the current region to be detected at the current time, inputting the material parameters into the deformation monitoring model, and outputting a deformation monitoring result based on the current boundary condition and the initial condition.
The material parameters mainly comprise elastic modulus, poisson ratio, shear strength and the like. And inputting the material parameters into a deformation monitoring model. Based on the boundary condition and the initial condition, the deformation monitoring model can be converted into a numerical model by adopting a finite element method, a boundary element method, a finite difference method or the like, and the numerical model is solved to obtain surface deformation data, namely a deformation monitoring result, and the surface deformation data is output.
And 104, acquiring deformation field data of the current region to be detected at the current time, and optimizing the deformation monitoring model according to the deformation field data and the deformation monitoring result to obtain an optimized deformation monitoring model.
When acquiring deformation field data, CR-InSAR technology can be adopted. When the CR-InSAR technology is adopted to acquire deformation field data, SAR image data are acquired by utilizing satellites or aircrafts; the interference pattern of the current region to be detected can be obtained by carrying out phase unwrapping, filtering, registering, difference and other treatments on SAR image data; and calculating deformation field data according to the interference pattern.
The CR-InSAR technology is adopted to acquire more accurate deformation field data. And optimizing the deformation monitoring model by taking the deformation field data as a reference until the deformation monitoring model with the optimized deformation is obtained.
In some embodiments, step 104 may include:
and carrying out sensitivity analysis on boundary conditions of the deformation monitoring model to obtain sensitivity analysis results of all the boundary conditions.
And determining the output error of the deformation monitoring model according to the deformation field data and the deformation monitoring result.
When the output error is larger than the error threshold, optimizing and adjusting each boundary condition according to the sensitivity analysis result and the output error of each boundary condition until the output error of the deformation monitoring model after optimizing and adjusting is smaller than the error threshold, and obtaining the deformation monitoring model after optimizing and finishing.
The main idea of the sensitivity analysis is to vary each input parameter over a range of input parameters for a given mathematical model and then evaluate the extent of influence of each input parameter by observing the variation of the output results. According to the embodiment of the application, the influence degree of each boundary condition on the deformation monitoring result can be obtained by carrying out sensitivity analysis on the boundary condition of the deformation monitoring model, so that the deformation monitoring model is conveniently optimized.
When sensitivity analysis is performed on boundary conditions, a reasonable boundary range is set for each boundary condition. The boundary range should cover the possible actual boundary conditions and the interactions between the boundary conditions need to be taken into account. And then, sequentially carrying out simulation calculation on all boundary conditions in the boundary range to obtain a model calculation result. In the simulation calculation, a numerical simulation method such as a finite element method may be used. And sequentially analyzing the deviation between each simulation calculation result and the actually measured deformation field data, and determining the sensitivity of each boundary condition to the deformation monitoring result to obtain a sensitivity analysis result. In determining the sensitivity analysis results for each boundary condition, analysis may be performed by sensitivity analysis software.
According to the sensitivity analysis result, determining which boundary conditions need to be optimally adjusted and the optimal adjustment direction of each boundary condition; respectively carrying out optimization adjustment according to the optimization adjustment direction of each boundary condition to obtain an optimized deformation monitoring model; re-outputting a deformation monitoring result by adopting the optimized and adjusted deformation monitoring model, and calculating the deviation between the deformation monitoring result and deformation field data, namely outputting an error; and when the output error is greater than the error threshold, continuing to perform optimization adjustment on each boundary condition of the deformation monitoring model until the output error is less than or equal to the error threshold, and completing optimization adjustment to obtain the optimized deformation monitoring model.
The optimized deformation monitoring model can be used for monitoring the surface deformation of the current region to be detected. And inputting the material parameters of the current region to be detected at the current time into the optimized deformation monitoring model, and correspondingly outputting the deformation monitoring result of the current region to be detected by the deformation monitoring model.
Compared with the prior art, the embodiment of the application has the beneficial effects that:
according to the embodiment of the application, the deformation monitoring model is established for monitoring the surface deformation of the current region to be detected, a large amount of data analysis and data correction work are not needed, and the deformation monitoring result can be obtained by inputting the material parameters at the current moment into the optimized deformation monitoring model. According to the embodiment of the application, the deformation monitoring model is optimized and adjusted by adopting actually measured deformation field data in advance, so that the optimized and adjusted deformation monitoring model is not limited by the professional level of a user, and a deformation monitoring result with higher accuracy is correspondingly output in real time according to the input material parameters. In addition, considering that the mining activities of different mineral resources have different influence degrees on the surface deformation, the embodiment of the application correspondingly establishes different types of deformation monitoring models aiming at different mineral resources, thereby being convenient for monitoring the surface deformation caused by different mineral resources in a targeted manner and further improving the accuracy of the deformation monitoring result.
It should be understood that the sequence number of each step in the foregoing embodiment does not mean that the execution sequence of each process should be determined by the function and the internal logic, and should not limit the implementation process of the embodiment of the present application.
The following are device embodiments of the application, for details not described in detail therein, reference may be made to the corresponding method embodiments described above.
Fig. 2 is a schematic structural diagram of a deformation monitoring model building device according to an embodiment of the present application, and for convenience of explanation, only a portion related to the embodiment of the present application is shown, which is described in detail below:
as shown in fig. 2, the deformation monitoring model building apparatus 2 includes: a modeling module 21 and an optimization module 22.
The modeling module 21 is configured to determine a type of deformation monitoring model corresponding to the current area to be measured according to mineral resources in the current area to be measured.
The modeling module 21 is further configured to establish a boundary condition and an initial condition of the deformation monitoring model according to the type of the deformation monitoring model.
The optimizing module 22 is configured to obtain a material parameter of the current area to be measured at the current time, input the material parameter to the deformation monitoring model, and output a deformation monitoring result based on the current boundary condition and the initial condition.
The optimizing module 22 is further configured to obtain deformation field data of the current area to be measured at the current time, and optimize the deformation monitoring model according to the deformation field data and the deformation monitoring result, so as to obtain an optimized deformation monitoring model.
In one possible implementation, the modeling module 21 is configured to determine, when the current area to be measured is a coal mine, that the deformation monitoring model corresponding to the current area to be measured is an elastomer model.
The modeling module 21 is further configured to determine, when the current area to be measured is a nonferrous metal ore, that the deformation monitoring model corresponding to the current area to be measured is a porous medium model.
In one possible implementation, when the deformation monitoring model is an elastomer model, the boundary conditions include: force boundary conditions, displacement boundary conditions, equilibrium boundary conditions, and material property boundary conditions;
a modeling module 21 for according toEstablishing a force boundary condition of the elastomer model;
wherein sigma n Represents boundary normal stress, u represents displacement, n represents normal vector, f n Representing the force exerted by the boundary.
The modeling module 21 is also used for the following stepsEstablishing an elastomer moldA displacement boundary condition of the model; wherein (1)>Represents the boundary of the current region to be measured omega, u 0 Representing a displacement of the boundary fixation.
The modeling module 21 is also used for the following stepsEstablishing balance boundary conditions of the elastomer model; where ε represents strain and σ represents stress.
The modeling module 21 is also used for the following stepsEstablishing a material property boundary condition of the elastomer model; where g represents the material properties at the boundary.
In one possible implementation, the modeling module 21 is configured to, according to
Establishing initial conditions of an elastomer model;
where u (x, y, z, t=0) represents the time t=0, the displacement at the position of coordinates (x, y, z), u_0 (x, y, z) represents the initial displacement at the position, v_0 (x, y, z) represents the initial velocity at the position, a represents the first amplitude of vibration, lx represents the length of the coal mine, ly represents the width of the coal mine, lz represents the height of the coal mine, ω represents the circular frequency of vibration, B represents the second amplitude of vibration, σ (x, y, z, t=0) represents the time t=0, the stress at the position of coordinates (x, y, z), σ_0 (x, y, z) represents the initial stress at the position, C represents the third amplitude of vibration, and D represents the fourth amplitude of vibration.
In one possible implementation, when the deformation monitoring model is a porous medium model, the boundary conditions include: displacement boundary conditions, force boundary conditions, and material property boundary conditions.
A modeling module 21 for according toEstablishing displacement boundary conditions of the porous medium model;
where u represents displacement, n represents normal vector, and g_n represents normal displacement boundary condition.
A modeling module 21, further configured to establish a force boundary condition of the porous medium model according to σn=f_n;
where σ represents stress, n represents normal vector, and f_n represents normal stress boundary condition.
The modeling module 21 is further configured to establish a material property boundary condition of the porous medium model according to epsilon=epsilon sigma;
where ε represents strain.
In one possible implementation, the modeling module 21 is configured to establish an initial displacement condition of the porous medium model according to u (x, y, z, t=0) =0;
where u (x, y, z, t=0) represents the displacement at the time t=0 at the position of coordinates (x, y, z).
The modeling module 21 is also configured to generate a model according to σ ij (x,y,z,t＝0)＝σ 0 δ ij Establishing an initial stress condition of a porous medium model;
wherein sigma ij (x, y, z, t=0) represents the stress tensor at the position of coordinates (x, y, z), σ, at time t=0 0 Represents an initial stress constant, δ when i=j ij =1, otherwise, δ ij ＝0。
In one possible implementation, the optimization module 22 is configured to perform sensitivity analysis on boundary conditions of the deformation monitoring model, so as to obtain sensitivity analysis results of each boundary condition.
The optimization module 22 is further configured to determine an output error of the deformation monitoring model according to the deformation field data and the deformation monitoring result.
And the optimization module 22 is further configured to perform optimization adjustment on each boundary condition according to the sensitivity analysis result and the output error of each boundary condition when the output error is greater than the error threshold, until the output error of the deformation monitoring model after optimization adjustment is less than the error threshold, and obtain the deformation monitoring model after optimization.
Compared with the prior art, the embodiment of the application has the beneficial effects that:
the modeling module 21 is used for performing surface deformation monitoring on the current region to be detected by establishing a deformation monitoring model, and can obtain a deformation monitoring result by inputting material parameters at the current moment into the optimized deformation monitoring model without performing a large amount of data analysis and data correction work. The optimization module 22 performs optimization adjustment on the deformation monitoring model by adopting actually measured deformation field data in advance, so that the deformation monitoring model after optimization adjustment can be not limited by the professional level of a user, and correspondingly outputs a deformation monitoring result with higher accuracy according to the input material parameters in real time. Moreover, considering that the mining activities of different mineral resources have different influence degrees on the surface deformation, the modeling module 21 correspondingly establishes different types of deformation monitoring models aiming at the different mineral resources, so that the surface deformation caused by the different mineral resources can be monitored in a targeted manner, and the accuracy of the deformation monitoring result is further improved.
Fig. 3 is a schematic diagram of an electronic device according to an embodiment of the present application. As shown in fig. 3, the electronic apparatus 3 of this embodiment includes: a processor 30, a memory 31 and a computer program 32 stored in said memory 31 and executable on said processor 30. The steps of the above-described embodiments of the method for creating the deformation monitoring model, such as steps 101 to 104 shown in fig. 1, are implemented by the processor 30 when executing the computer program 32. Alternatively, the processor 30 may perform the functions of the modules/units of the apparatus embodiments described above, such as the functions of the modules 21 to 23 shown in fig. 2, when executing the computer program 32.
Illustratively, the computer program 32 may be partitioned into one or more modules/units that are stored in the memory 31 and executed by the processor 30 to complete the present application. The one or more modules/units may be a series of computer program instruction segments capable of performing the specified functions for describing the execution of the computer program 32 in the electronic device 3. For example, the computer program 32 may be divided into modules 21 to 23 shown in fig. 2.
The electronic device 3 may be a computing device such as a desktop computer, a notebook computer, a palm computer, a cloud server, etc. The electronic device 3 may include, but is not limited to, a processor 30, a memory 31. It will be appreciated by those skilled in the art that fig. 3 is merely an example of the electronic device 3 and does not constitute a limitation of the electronic device 3, and may include more or fewer components than shown, or may combine certain components, or different components, e.g., the electronic device may further include an input-output device, a network access device, a bus, etc.
The processor 30 may be a central processing unit (Central Processing Unit, CPU), other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), field-programmable gate arrays (Field-Programmable Gate Array, FPGA) 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 31 may be an internal storage unit of the electronic device 3, such as a hard disk or a memory of the electronic device 3. The memory 31 may be an external storage device of the electronic device 3, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card) or the like, which are provided on the electronic device 3. Further, the memory 31 may also include both an internal storage unit and an external storage device of the electronic device 3. The memory 31 is used for storing the computer program and other programs and data required by the electronic device. The memory 31 may also be used for temporarily storing data that has been output or is to be output.
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.
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/electronic device and method may be implemented in other manners. For example, the apparatus/electronic device embodiments described above are merely illustrative, e.g., the division of the modules or units is merely a logical function 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.
In addition, each functional unit in the embodiments of the present application 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. The integrated units may be implemented in hardware or in software functional units.
The integrated modules/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 embodiment, or 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 the computer program may implement the steps of the method embodiment of establishing the deformation monitoring model when the computer program is executed by a processor. 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 may include: any entity or device capable of carrying the computer program code, a recording medium, a U disk, a removable hard disk, a magnetic disk, an optical disk, a computer Memory, a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), an electrical carrier signal, a telecommunications signal, a software distribution medium, and so forth. The above embodiments are only for illustrating the technical solution of the present application, and not for limiting the same; although the application has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present application, and are intended to be included in the scope of the present application.
Claims (10)
1. The method for establishing the deformation monitoring model is characterized by comprising the following steps of:
determining the type of a deformation monitoring model corresponding to a current region to be detected according to mineral resources in the current region to be detected;
establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model;
acquiring material parameters of a current region to be detected at a current time, inputting the material parameters into the deformation monitoring model, and outputting a deformation monitoring result based on a current boundary condition and an initial condition;
and obtaining deformation field data of the current region to be detected at the current time, and optimizing the deformation monitoring model according to the deformation field data and the deformation monitoring result to obtain an optimized deformation monitoring model.
2. The method for building a deformation monitoring model according to claim 1, wherein the determining the type of the deformation monitoring model corresponding to the current area to be measured according to mineral resources in the current area to be measured comprises:
when the current region to be detected is a coal mine, determining a deformation monitoring model corresponding to the current region to be detected as an elastomer model;
when the current area to be measured is nonferrous metal ore, determining the deformation monitoring model corresponding to the current area to be measured as a porous medium model.
3. The method for building a deformation monitoring model according to claim 2, wherein when the deformation monitoring model is an elastomer model, the boundary condition includes: force boundary conditions, displacement boundary conditions, equilibrium boundary conditions, and material property boundary conditions;
and establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model, wherein the boundary conditions and initial conditions comprise:
according toEstablishing a force boundary condition of the elastomer model;
wherein sigma n Represents boundary normal stress, u represents displacement, n represents normal vector, f n Representing the force exerted by the boundary;
according toEstablishing displacement boundary conditions of the elastomer model; wherein (1)>Represents the boundary of the current region to be measured omega, u 0 Representing the displacement of the boundary fixation;
according toEstablishing equilibrium boundary conditions of the elastomer model; wherein ε represents strain and σ represents stress;
according toEstablishing a material property boundary condition of the elastomer model; where g represents the material properties at the boundary.
4. A method of constructing a deformation monitoring model according to claim 3, wherein the constructing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model further comprises:
according to
Establishing initial conditions of the elastomer model;
where u (x, y, z, t=0) represents the time t=0, the displacement at the position with coordinates (x, y, z), u_0 (x, y, z) represents the initial displacement at the position, v_0 (x, y, z) represents the initial velocity at the position, a represents the first amplitude of vibration, lx represents the length of the coal mine, ly represents the width of the coal mine, lz represents the height of the coal mine, ω represents the circular frequency of vibration, B represents the second amplitude of vibration, σ (x, y, z, t=0) represents the time t=0, the stress at the position with coordinates (x, y, z), σ_0 (x, y, z) represents the initial stress at the position, C represents the third amplitude of vibration, and D represents the fourth amplitude of vibration.
5. The method for building a deformation monitoring model according to claim 2, wherein when the deformation monitoring model is a porous medium model, the boundary condition includes: displacement boundary conditions, force boundary conditions, and material property boundary conditions;
and establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model, wherein the boundary conditions and initial conditions comprise:
according toEstablishing a displacement boundary condition of the porous medium model;
wherein u represents displacement, n represents normal vector, and g_n represents normal displacement boundary condition;
establishing a force boundary condition of the porous medium model according to σn=f_n;
wherein sigma represents stress, n represents normal vector, and f_n represents normal stress boundary condition;
establishing a material property boundary condition of the porous medium model according to epsilon=epsilon sigma;
where ε represents strain.
6. The method for building a deformation monitoring model according to claim 5, wherein the building boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model further comprises:
establishing an initial displacement condition of the porous medium model according to u (x, y, z, t=0) =0;
where u (x, y, z, t=0) represents the displacement at the position of the coordinates (x, y, z) at time t=0;
according to sigma ij (x,y,z,t＝0)＝σ 0 δ ij Establishing an initial stress condition of the porous medium model;
wherein sigma ij (x, y, z, t=0) represents the stress tensor at the position of coordinates (x, y, z), σ, at time t=0 0 Represents an initial stress constant, δ when i=j ij =1, otherwise, δ ij ＝0。
7. The method for building a deformation monitoring model according to any one of claims 1 to 6, wherein optimizing the deformation monitoring model according to the deformation field data and the deformation monitoring result to obtain an optimized deformation monitoring model comprises:
performing sensitivity analysis on boundary conditions of the deformation monitoring model to obtain sensitivity analysis results of all the boundary conditions;
determining an output error of the deformation monitoring model according to the deformation field data and the deformation monitoring result;
and when the output error is larger than an error threshold, optimizing and adjusting each boundary condition according to the sensitivity analysis result of each boundary condition and the output error until the output error of the deformation monitoring model after optimizing and adjusting is smaller than the error threshold, and obtaining the deformation monitoring model after optimizing and finishing.
8. A deformation monitoring model building device, characterized by comprising:
the modeling module is used for determining the type of the deformation monitoring model corresponding to the current region to be detected according to mineral resources in the current region to be detected;
the modeling module is further used for establishing boundary conditions and initial conditions of the deformation monitoring model according to the type of the deformation monitoring model;
the optimization module is used for acquiring material parameters of the current region to be detected at the current time, inputting the material parameters into the deformation monitoring model, and outputting a deformation monitoring result based on the current boundary condition and the initial condition;
the optimization module is used for acquiring deformation field data of the current region to be detected at the current time, and optimizing the deformation monitoring model according to the deformation field data and the deformation monitoring result to obtain an optimized deformation monitoring model.
9. An electronic device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor, when executing the computer program, realizes the steps of the method of building a deformation monitoring model according to any of the preceding claims 1 to 7.
10. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the method of establishing a deformation monitoring model according to any of the preceding claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310919294.7A CN116952119A (en) | 2023-07-25 | 2023-07-25 | Deformation monitoring model building method and device, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310919294.7A CN116952119A (en) | 2023-07-25 | 2023-07-25 | Deformation monitoring model building method and device, electronic equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116952119A true CN116952119A (en) | 2023-10-27 |
Family
ID=88456123
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310919294.7A Pending CN116952119A (en) | 2023-07-25 | 2023-07-25 | Deformation monitoring model building method and device, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116952119A (en) |
-
2023
- 2023-07-25 CN CN202310919294.7A patent/CN116952119A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Priest | Determination of discontinuity size distributions from scanline data | |
| Xiang et al. | Experimental investigation of frequency-based multi-damage detection for beams using support vector regression | |
| CN110954950B (en) | Underground transverse wave velocity inversion method, device, computing equipment and storage medium | |
| WO2017120088A1 (en) | Effective permeability upscaling for a discrete fracture network | |
| US12073154B2 (en) | Pressure cancelation of unsteady boundary conditions during simulation of fluid flows | |
| Louboutin et al. | Performance prediction of finite-difference solvers for different computer architectures | |
| Liu et al. | Shift boundary material point method: an image-to-simulation workflow for solids of complex geometries undergoing large deformation | |
| CN112578450A (en) | Forward modeling method and device based on finite difference of frequency dispersion medium scalar wave equation | |
| US20150073730A1 (en) | Mechanical strain gauge simulation | |
| Loyola et al. | General statistics-based methodology for the determination of the geometrical and mechanical representative elementary volumes of fractured media | |
| Xu et al. | Modal parameter identification of a quayside container crane based on data-driven stochastic subspace identification | |
| EP3211594B1 (en) | Seismic modeling system providing seismic survey data inpainting based upon suspect region boundary comparisons and related methods | |
| EP3217354B1 (en) | Seismic modeling system providing seismic survey data frequency domain inpainting and related methods | |
| Bathe | The finite element method with “overlapping finite elements” | |
| CN111679318A (en) | Method and system for identifying fracture system in reservoir | |
| CN111665556A (en) | Method for constructing stratum acoustic wave propagation velocity model | |
| CN116952119A (en) | Deformation monitoring model building method and device, electronic equipment and storage medium | |
| Madenci et al. | Peridynamics enabled digital image correlation for tracking crack paths | |
| CN109581494B (en) | Pre-stack migration method and device | |
| Li et al. | A 3D reflection ray-tracing method based on linear traveltime perturbation interpolation | |
| CN108646288B (en) | Method and device for establishing near-surface model | |
| CN112649848B (en) | Method and device for solving earthquake wave impedance by utilizing wave equation | |
| Goldfeld | Identification of the stiffness distribution in statically indeterminate beams | |
| CN111665550A (en) | Underground medium density information inversion method | |
| CN111665549A (en) | Inversion method of stratum acoustic wave attenuation factor |
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 |