CN-122021185-A - Three-dimensional ground stress field inversion method, system and platform based on numerical simulation

CN122021185ACN 122021185 ACN122021185 ACN 122021185ACN-122021185-A

Abstract

The application relates to the technical field of ground stress inversion, in particular to a three-dimensional ground stress field inversion method, system and platform based on numerical simulation, wherein the method comprises the following steps: the earth coordinates of all the earth stress measuring points in the research area and the three-dimensional earth stress data are analyzed to calculate the weight of the residual error corresponding to each earth stress measuring point in the least square method used by the earth stress regression equation to solve the optimized earth stress regression equation, and the earth stress regression equation after the optimization is solved is inverted in the three-dimensional earth stress field to complete the three-dimensional earth stress field inversion of the research area. The application aims to improve the authenticity and reliability of the inversion result of the three-dimensional ground stress field.

Inventors

LI ZHIGUO
LI JIN
YANG MIN
Feng Guipu
YANG MENGQIONG
GUO HONGYANG

Assignees

武汉工程大学
湖北神农磷业科技股份有限公司

Dates

Publication Date: 20260512
Application Date: 20260331

Claims (10)

1. The three-dimensional ground stress field inversion method based on numerical simulation is characterized by comprising the following steps of: Constructing a three-dimensional finite element model of a research area, which is used for carrying out numerical simulation on measured stress values of each ground stress measuring point under a dead weight force field and a preset number of geological structure stress fields so as to obtain a ground stress regression equation, and is used for simulating stress calculation values of all positions in the three-dimensional finite element model to finish inversion of the three-dimensional ground stress field; The optimization solving method of the ground stress regression equation comprises the following steps: S1, collecting three-dimensional ground stress data and ground coordinates of all preset ground stress measuring points in a research area; S2, calculating a first data deviation coefficient of any ground stress measuring point according to the variation difference of the linear distribution rule of the vertical stress along with the depth of the ground in the all ground stress measuring points before and after any ground stress measuring point is eliminated in a research area, and evaluating the degree of data deviation of the vertical stress acquired by any ground stress measuring point due to noise interference; S3, calculating a second data deviation coefficient of any ground stress measuring point according to the variation difference of the nonlinear distribution rule of the difference along with the depth of the ground between the maximum horizontal stress and the minimum horizontal stress in all ground stress measuring points before and after any ground stress measuring point in the research area, and evaluating the degree of data deviation of the horizontal stress acquired by any ground stress measuring point due to noise interference; S4, fusing the first data deviation coefficient and the second data deviation coefficient to determine the confidence coefficient of any ground stress measuring point, and carrying out normalized mapping on the confidence coefficient of all the ground stress measuring points to obtain a weight value, wherein the weight value is used as the weight of the residual error corresponding to each ground stress measuring point in the least square method used for carrying out optimal solution on the ground stress regression equation so as to solve the ground stress regression equation.
2. The numerical simulation-based three-dimensional ground stress field inversion method according to claim 1, wherein the first data deviation coefficient has a positive correlation with the goodness of fit of a linear regression model constructed according to a linear distribution rule after any ground stress measuring point is excluded, and has a negative correlation with the goodness of fit of a linear regression model constructed according to a linear distribution rule before any ground stress measuring point is excluded.
3. The numerical simulation-based three-dimensional ground stress field inversion method according to claim 2, wherein the goodness of fit is determined by the decision coefficients of a linear regression model constructed between the two data of the subsurface depth and the vertical stress of all ground stress measuring points.
4. The numerical simulation-based three-dimensional ground stress field inversion method according to claim 1, wherein the second data deviation coefficient is in positive correlation with the goodness of fit of the nonlinear regression model constructed according to the nonlinear distribution law after any ground stress measuring point is excluded, and in negative correlation with the goodness of fit of the nonlinear regression model constructed according to the nonlinear distribution law before any ground stress measuring point is excluded.
5. The numerical simulation-based three-dimensional ground stress field inversion method according to claim 4, wherein the goodness of fit is determined by the decision coefficients of a nonlinear regression model between two data, namely the difference between the horizontal maximum and minimum principal stresses of all ground stress measuring points and the subsurface depth.
6. The numerical simulation-based three-dimensional ground stress field inversion method according to claim 1, wherein the confidence level is in negative correlation with the first data deviation coefficient and the second data deviation coefficient.
7. The numerical simulation-based three-dimensional ground stress field inversion method according to claim 6, wherein the sum of weight values of all ground stress measuring points obtained by normalized mapping of the confidence is 1.
8. The numerical simulation-based three-dimensional ground stress field inversion method according to any of claims 1-7, wherein the geological formation stress field comprises an x-axis extrusion, a y-axis extrusion, an x-y plane shear, a y-z plane shear, and a z-x plane shear formation stress field.
9. A numerical simulation-based three-dimensional ground stress field inversion system, comprising: the three-dimensional finite element model construction module is used for constructing a three-dimensional finite element model of the research area and transmitting the three-dimensional finite element model to the ground stress numerical simulation module and the three-dimensional ground stress field inversion module; the three-dimensional ground stress data acquisition module is used for realizing the step S1 in the claim 1 and transmitting data to the ground stress numerical simulation module and the three-dimensional ground stress data analysis module; The ground stress numerical simulation module is used for carrying out numerical simulation on measured stress values of all ground stress measuring points under the gravity field and the preset number of geological structure stress fields by utilizing the three-dimensional finite element model so as to acquire a ground stress regression equation, and transmitting the ground stress regression equation to the three-dimensional ground stress data analysis module; The three-dimensional ground stress data analysis module is used for realizing the steps S2-S4 in the claim 1 and transmitting the ground stress regression equation for optimization solution to the three-dimensional ground stress field inversion module; And the three-dimensional ground stress field inversion module is used for carrying out inversion in the three-dimensional finite element model according to a ground stress regression equation of the optimization solution.
10. A numerical simulation based three-dimensional ground stress field inversion platform comprising a memory, a processor and a computer program stored in the memory and running on the processor, characterized in that the processor implements the numerical simulation based three-dimensional ground stress field inversion method according to any of claims 1-7 when executing the computer program.

Description

Three-dimensional ground stress field inversion method, system and platform based on numerical simulation Technical Field The application relates to the technical field of ground stress inversion, in particular to a three-dimensional ground stress field inversion method, system and platform based on numerical simulation. Background The three-dimensional ground stress field inversion is to build a mathematical model through known ground stress field actual measurement data and reversely deduce the initial ground stress distribution state of the underground rock mass in a three-dimensional space. The three-dimensional ground stress inversion result can provide important basis for geological engineering design and geological disaster evaluation, and is widely applied to the fields of geological engineering, oil-gas exploration, geological disaster prevention and control and the like. When inversion of a three-dimensional ground stress field is carried out, the existing method generally carries out numerical simulation on ground stress of ground stress measuring points of a research area under different ground stress influence factors according to a three-dimensional finite element model constructed by geological data of the research area to obtain a stress calculation value, then utilizes a least square method to solve a ground stress regression equation established by the stress calculation value and the three-dimensional ground stress data of the ground stress measuring points of the research area, and finally carries out finite element calculation on the ground stress regression equation obtained by solving to generate the three-dimensional ground stress field in a numerical simulation software platform. However, the method ignores that the three-dimensional ground stress data of each ground stress measuring point can have different degrees of data deviation due to different noise interferences (such as electromagnetic interference, mechanical vibration, instrument self drift and the like) which are included in the ground stress data, and the traditional least square method uses the same weight for the residual errors of all data, so that when the ground stress regression equation constructed by the three-dimensional ground stress data of all the ground stress measuring points is solved by using the least square method, the residual errors calculated by the three-dimensional ground stress data with smaller noise interference in all the three-dimensional ground stress data can not obtain higher weight in the least square method, and the finally solved ground stress regression equation can not reflect the real ground stress distribution situation, thereby influencing the inversion result of the three-dimensional ground stress field. Disclosure of Invention In order to solve the technical problems, the application aims to provide a three-dimensional ground stress field inversion method, system and platform based on numerical simulation, and the adopted technical scheme is as follows: In a first aspect, an embodiment of the present application provides a three-dimensional ground stress field inversion method based on numerical simulation, the method including the steps of: Constructing a three-dimensional finite element model of a research area, which is used for carrying out numerical simulation on measured stress values of each ground stress measuring point under a dead weight force field and a preset number of geological structure stress fields so as to obtain a ground stress regression equation, and is used for simulating stress calculation values of all positions in the three-dimensional finite element model to finish inversion of the three-dimensional ground stress field; The optimization solving method of the ground stress regression equation comprises the following steps: S1, collecting three-dimensional ground stress data and ground coordinates of all preset ground stress measuring points in a research area; S2, calculating a first data deviation coefficient of any ground stress measuring point according to the variation difference of the linear distribution rule of the vertical stress along with the depth of the ground in the all ground stress measuring points before and after any ground stress measuring point is eliminated in a research area, and evaluating the degree of data deviation of the vertical stress acquired by any ground stress measuring point due to noise interference; S3, calculating a second data deviation coefficient of any ground stress measuring point according to the variation difference of the nonlinear distribution rule of the difference along with the depth of the ground between the maximum horizontal stress and the minimum horizontal stress in all ground stress measuring points before and after any ground stress measuring point in the research area, and evaluating the degree of data deviation of the horizontal stress acquired by any ground stress measuring point due to noise inte