CN103605900A - Cross-scale complicated geologic body ground stress field identification method and device - Google Patents

Abstract

The invention relates to a cross-scale complicated geologic body ground stress field identification method, and belongs to the technical field of geological engineering, rock mass engineering and numerical analysis. The cross-scale complicated geologic body ground stress field identification method comprises the following steps: an optimization model capable of identifying ground stress is established; a genetic algorithm is adopted for solving the optimization model so as to obtain protolith parameters and side pressure coefficients of a geologic body; the protolith parameters of the geologic body are substituted into a three-dimensional numerical model to obtain a stress field. A device used for implementing the method comprises an optimization model establishing module, a genetic algorithm calculation module and a ground stress field obtaining module, wherein the optimization model establishing module is used for establishing the optimization model capable of identifying the ground stress; the genetic algorithm is used for solving the optimization model to obtain the protolith parameters and the side pressure coefficients of the geologic body; the ground stress field obtaining module is used for substituting the protolith parameters of the geologic body into the three-dimensional numerical model so as to obtain the ground stress field. According to the cross-scale complicated geologic body ground stress field identification method and device, the initial ground stress field in the geologic body can be obtained, so that the problem that existing ground stress inversion theories and methods are limited by ground stress inversion of a cross-scale complicated geologic body.

Description

Across the stress field recognition methods of scale complex geologic body and device

Technical field

The present invention relates to a kind of Geological Engineering, rock mass engineering project and numerical analysis techniques field, especially a kind of across the stress field recognition methods of scale complex geologic body and device.

Background technology

Terrestrial stress is the natural stress being present in stratum, is the internally-powered that causes rock mass engineering project deformation failure and rock burst dynamic disaster, is the key factor that rock mass engineering project stability analysis and damage control must be considered.

Geologic body has experienced very long tectonic movement and repeatedly geology transformation effect, seals the residual structure acting force of different times in stratum up for safekeeping.Differently the difference of the suffered geologic function type of plastid, effect degree, action time and experience number of times, causes the terrestrial stress size and Orientation of geologic body on room and time, to have variability.Especially the larger structural feature such as tomography, fold and shear zone, may cause the terrestrial stress distribution in geologic body to undergo mutation or reverse, across inverting and the Study on regularity of yardstick stress field, brought very big difficulty thus, become Geotechnical Engineering field always and study for a long period of time and insoluble technical barrier.

Geostress survey is still the most important means of current Study on Rock Stress.According to works scope, significance level and complexity, carry out the geostress survey quantity adapting with it, carry out on this basis regretional analysis and Study on Rock Stress.Owing to being subject to the restriction of time and funds, the geostress survey quantity of most of engineering is very limited, and measured the impact of the factors such as measure, toward contact, there is discreteness significantly in the terrestrial stress test result that obtained thus, thereby bring very large difficulty to stress field Study on regularity, be also difficult to one of three key factors of quantitative Application to thus the quantitative test of Geological Engineering.

In order to improve the reliability of stress field research, at present people studys and exploringly INVERSION OF STRESS FIELD theory and analogy method always.Eighties of last century has proposed Back Analysis of Viscoelastic method the eighties first, and the deformation monitoring by excavation project carries out Mechanics Parameters of Rock Mass and terrestrial stress inverting [12].Along with artificial technology's development, people have carried out a large amount of terrestrial stress inverting combining with numerical evaluation based on neural network research.In recent years, the size and Orientation of the deformation failure characteristic estimation stress field of people's application tectonic structure person's movements and expression and country rock.

The inverting of summing up the stress field of current geologic body has obtained remarkable progress, especially adopt artificial neural network to combine with numerical model, the calculated value of terrestrial stress of take approaches measured value as inverting foundation, for the inverting of complicated field stress field provides good research approach.But existing terrestrial stress inverting exists limitation to be, terrestrial stress inverting based on neural network or to limitedly training and the prediction of actual Stress Measurement result, or by means of expanding terrestrial stress sample, improves numerical model the precision of prediction of terrestrial stress, therefore also there is certain limitation for the terrestrial stress inverting across yardstick geologic body of complexity in existing terrestrial stress inversion theory and method.

Summary of the invention

For the weak point existing in the problems referred to above, the invention provides a kind of can obtain initially stress field in geologic body across the stress field recognition methods of scale complex geologic body and device.

For achieving the above object, the invention provides a kind ofly across the recognition methods of scale complex geologic body stress field, comprise the following steps:

The Optimized model that S100, foundation can be identified terrestrial stress;

S200, employing genetic algorithm solve Optimized model, to obtain protolith parameter and the coefficient of horizontal pressure of geologic body;

S300, by the protolith parameter substitution Three-dimension Numerical Model of geologic body, to obtain stress field.

Above-mentioned across the recognition methods of scale complex geologic body stress field, wherein, in step S100, comprising:

S101, employing genetic planning, according to known terrestrial stress measuring point, set up the funtcional relationship of terrestrial stress monitoring point terrestrial stress and protolith parameter and coefficient of horizontal pressure;

It is optimization aim that the quadratic sum of the terrestrial stress calculating value of S102, base area stress monitoring point and the difference of measured value reaches minimum, and sets up the Optimized model of the initial stress and coefficient of horizontal pressure.

Above-mentioned across the recognition methods of scale complex geologic body stress field, wherein, in step S101, j terrestrial stress measuring point (j=1,15) terrestrial stress (size and Orientation) and protolith parameter (density γ, elastic modulus E, Poisson ratio μ) and coefficient of horizontal pressure (x direction coefficient of horizontal pressure λ are set up in ground ₁with y direction coefficient of horizontal pressure λ ₂) between funtcional relationship as follows:

σ' _j1=f _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j2=f _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j3=f _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j1=g _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j2=g _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j3=g _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)。

Above-mentioned across the recognition methods of scale complex geologic body stress field, wherein, in step S102, the objective function of Optimized model is as follows:

$\mathrm{Min}\underset{j=1}{\overset{15}{\mathrm{\Σ}}}(\underset{i=1}{\overset{3}{\mathrm{\Σ}}}{({\mathrm{\σ}}_{\mathrm{ji}}^{\′}-{\mathrm{\σ}}_{\mathrm{ji}})}^2+\underset{i=1}{\overset{3}{\mathrm{\Σ}}}{({\mathrm{\α}}_{\mathrm{ji}}^{\′}-{\mathrm{\α}}_{\mathrm{ji}})}^2),$

Its constraint condition is as follows:

2.4<γ ₁<2.8，2.8<γ ₂<3.6，2.5<γ ₃<2.9，1.8<γ ₄<2.2

6<E ₁<20，20<E ₂<34，6<E ₃<20，0.8<E ₄<1.3

0.20<μ ₁<0.26，0.20<μ ₂<0.24，0.20<μ ₃<0.26，0.35<μ ₄<0.45

0.8<λ ₁<1.2，1.1<λ ₂<1.5。

Above-mentioned across the recognition methods of scale complex geologic body stress field, wherein, in step S200, protolith parameter and the coefficient of horizontal pressure of geologic body are as follows:

Upper panel protolith parameter:

The protolith parameter of ore body:

Lower panel protolith parameter:

Tomography protolith parameter:

Mining area coefficient of horizontal pressure:

Above-mentioned across the recognition methods of scale complex geologic body stress field, wherein, in step S300, by obtained geologic body protolith parameter

substitution Three-dimension Numerical Model calculates, thus obtained stress field.

It is a kind of across scale complex geologic body stress field recognition device that the present invention also provides, and comprising:

Optimize module and set up module, for setting up the Optimized model that can identify terrestrial stress;

Genetic algorithm computing module, for Optimized model is solved, to obtain protolith parameter and the coefficient of horizontal pressure of geologic body;

Stress field obtains module, for by the protolith parameter substitution Three-dimension Numerical Model of geologic body, to obtain stress field.

Above-mentioned device, wherein, sets up in module in described optimization module, comprising:

Adopt genetic planning, according to known terrestrial stress measuring point, set up the funtcional relationship of terrestrial stress monitoring point terrestrial stress and protolith parameter and coefficient of horizontal pressure;

It is optimization aim that the quadratic sum of the terrestrial stress calculating value of base area stress monitoring point and the difference of measured value reaches minimum, and sets up the Optimized model of the initial stress and coefficient of horizontal pressure.

Above-mentioned device, wherein, after genetic algorithm computing module is processed, protolith parameter and the coefficient of horizontal pressure of geologic body are as follows:

Upper panel protolith parameter:

The protolith parameter of ore body:

Lower panel protolith parameter:

Tomography protolith parameter:

Mining area is stress coefficient of horizontal pressure initially:

Above-mentioned device, wherein, obtains module by stress field, by obtained geologic body protolith parameter substitution Three-dimension Numerical Model calculates, thus obtained stress field.

Compared with prior art, the present invention has the following advantages:

1, utilize the function of numerical modeling, different types of lithology of geologic body and rift structure are carried out to subregion processing, realize the impact of geologic body Spatial Variability on terrestrial stress;

2, the model of mind of employing genetic planning is set up the relation of terrestrial stress (size and Orientation) with protolith parameter and the coefficient of horizontal pressure of terrestrial stress measuring point, not only solve neural network model and be difficult to provide explicit function, and can simulate nonlinearity function therebetween with the given function of tree structure;

3, the minimum of the quadratic sum based on the calculated value of all terrestrial stress measuring points in geologic body and the difference of measured value is optimization aim, adopts genetic algorithm to carry out the optimizing of protolith parameter and coefficient of horizontal pressure, can obtain total optimization parameter;

4, the protolith parameter solving and coefficient of horizontal pressure substitution Three-dimension Numerical Model are carried out to Stress calculation, directly obtain thus the primary stress field of geologic body and boundary condition as the initial stress of engineering numerical analysis; Give the protolith parameter of geologic body simultaneously.

Accompanying drawing explanation

Fig. 1 is the process flow diagram of method part in the present invention;

Fig. 2 is the structural drawing of device part in the present invention;

Fig. 3 is the typical geology sectional view of the large-scale Copper-nickel Deposits in Jinchuan.

Main element symbol description is as follows:

1-optimizes module and sets up module 2-genetic algorithm computing module

3-stress field obtains module

Embodiment

As shown in Figure 1, the invention provides a kind ofly across the recognition methods of scale complex geologic body stress field, comprise the following steps:

The Optimized model that S100, foundation can be identified terrestrial stress.

Concrete, in step S100, comprising:

S101, employing genetic planning, according to known terrestrial stress measuring point, set up the funtcional relationship of terrestrial stress monitoring point terrestrial stress and protolith parameter and coefficient of horizontal pressure, and its funtcional relationship is as follows:

σ' _j1=f _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j2=f _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j3=f _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j1=g _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j2=g _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j3=g _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)。

It is optimization aim that the quadratic sum of the terrestrial stress calculating value of S102, base area stress monitoring point and the difference of measured value reaches minimum, and sets up the Optimized model of the initial stress and coefficient of horizontal pressure, and in step S102, the objective function of Optimized model is as follows:

$\mathrm{Min}\underset{j=1}{\overset{15}{\mathrm{\&Sigma;}}}(\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{({\mathrm{\&sigma;}}_{\mathrm{ji}}^{\&prime;}-{\mathrm{\&sigma;}}_{\mathrm{ji}})}^2+\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{({\mathrm{\&alpha;}}_{\mathrm{ji}}^{\&prime;}-{\mathrm{\&alpha;}}_{\mathrm{ji}})}^2),$

Its constraint condition is as follows:

2.4<γ ₁<2.8，2.8<γ ₂<3.6，2.5<γ ₃<2.9，1.8<γ ₄<2.2

6<E ₁<20，20<E ₂<34，6<E ₃<20，0.8<E ₄<1.3

0.20<μ ₁<0.26，0.20<μ ₂<0.24，0.20<μ ₃<0.26，0.35<μ ₄<0.45

0.8<λ ₁<1.2，1.1<λ ₂<1.5。

S200, employing genetic algorithm solve Optimized model, to obtain protolith parameter and the coefficient of horizontal pressure of geologic body.

Concrete, in step S200, protolith parameter and the coefficient of horizontal pressure of geologic body are as follows:

Upper panel protolith parameter:

The protolith parameter of ore body:

Lower panel protolith parameter:

Tomography protolith parameter:

Mining area coefficient of horizontal pressure:

S300, by the protolith parameter substitution Three-dimension Numerical Model of geologic body, to obtain stress field.

Concrete, in step S300, by obtained geologic body protolith parameter

substitution Three-dimension Numerical Model calculates, thus obtained stress field.

As shown in Figure 2, the invention provides a kind ofly across scale complex geologic body stress field recognition device, comprise and optimize that module is set up module 1, genetic algorithm computing module 2 obtains module 3 with stress field.

Optimize module and set up module 1 for setting up the Optimized model that can identify terrestrial stress.

Specifically comprise:

Adopt genetic planning, according to known terrestrial stress measuring point, set up the funtcional relationship of terrestrial stress monitoring point terrestrial stress and protolith parameter and coefficient of horizontal pressure, its funtcional relationship is as follows:

σ' _j1=f _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j2=f _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j3=f _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j1=g _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j2=g _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j3=g _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)。

It is optimization aim that the quadratic sum of the terrestrial stress calculating value of base area stress monitoring point and the difference of measured value reaches minimum, and sets up the Optimized model of the initial stress and coefficient of horizontal pressure, and the objective function of Optimized model is as follows:

$\mathrm{Min}\underset{j=1}{\overset{15}{\mathrm{\&Sigma;}}}(\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{({\mathrm{\&sigma;}}_{\mathrm{ji}}^{\&prime;}-{\mathrm{\&sigma;}}_{\mathrm{ji}})}^2+\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{({\mathrm{\&alpha;}}_{\mathrm{ji}}^{\&prime;}-{\mathrm{\&alpha;}}_{\mathrm{ji}})}^2),$

Its constraint condition is as follows:

2.4<γ ₁<2.8，2.8<γ ₂<3.6，2.5<γ ₃<2.9，1.8<γ ₄<2.2

6<E ₁<20，20<E ₂<34，6<E ₃<20，0.8<E ₄<1.3

0.20<μ ₁<0.26，0.20<μ ₂<0.24，0.20<μ ₃<0.26，0.35<μ ₄<0.45

0.8<λ ₁<1.2，1.1<λ ₂<1.5。

Genetic algorithm computing module 2 is for Optimized model is solved, to obtain protolith parameter and the coefficient of horizontal pressure of geologic body.After genetic algorithm computing module is processed, protolith parameter and the coefficient of horizontal pressure of geologic body are as follows:

Upper panel protolith parameter:

The protolith parameter of ore body:

Lower panel protolith parameter:

Tomography protolith parameter:

Mining area is stress coefficient of horizontal pressure initially:

Stress field obtains module 3, for by the protolith parameter substitution Three-dimension Numerical Model of geologic body, to obtain stress field.By stress field, obtain module, by obtained geologic body protolith parameter

substitution Three-dimension Numerical Model calculates, thus obtained stress field.

Experiment case study:

The present invention be take Jinchuan large ore deposit as case history, and this working of an invention process is described.Figure 3 shows that the typical log sheet of the large-scale Copper-nickel Deposits in Jinchuan.Terrestrial stress identification step is as follows thus:

1, set up geologic model

Level by means of mineral deposit is carried out geologic body subregion and rift structure identification (as Fig. 1) with vertical geologic section, according to the lithology of geology of mineral deposit body and tectonic structure, geologic body is divided into upper dish, ore body and three geologic division of lower wall, and considers the fault tectonic of lower wall.

2, set up mechanical model

1) set up the coordinate system of numerical model

It is x coordinate that surface level points to east, and energized north y coordinate, vertically upward for z coordinate is set up coordinate system.

2) determine geologic division rock mass parameter

According to three of geology of mineral deposit model subregions and a tomography, determine corresponding rock mass parameter.Be not subject to rock mass that engineering excavation affects in elastic stage, determine that thus the geologic body rock mass parameter of each subregion is as follows:

(1) upper dish rock mass unit weight, elastic modulus and Poisson ratio are γ ₁, E ₁, μ ₁;

(2) ore body dish rock mass unit weight, elastic modulus and Poisson ratio are γ ₂, E ₂, μ ₂;

(3) lower wall rock mass unit weight, elastic modulus and Poisson ratio are γ ₃, E ₃, μ ₃;

(4) shear stiffness of lower wall tomography and normal stiffness γ ₄, E ₄, μ ₄.

3, determine numerical procedure

Complicated geological mass ground stress state is conventionally by gravity and the tectonic stress λ of geologic body ₁, λ ₂(coefficient of horizontal pressure of x direction and y direction) totally 14 factors forms.Therefore, adopt the orthogonal numerical analysis of 14 factor 2 levels to simulate the stress field under various combination condition.Rule of thumb, with mining area geostress survey value, determine that 2 levels of 14 factors are in Table 1.

14 factors and 2 level values of table 1 Jinchuan copper-nickel mine bed quadrature analysis

3, actual measurement terrestrial stress measuring point is selected

1) terrestrial stress measuring point selection principle

The geologic body of different engineerings is owing to being subject to the restriction of engineering type, significance level and time and funds, and method and the quantity of carrying out on-the-spot geostress survey there are differences.When the quantity of in-site measurement is not more than 15 measuring points, will all adopt the test result of measuring point to carry out Intelligent Recognition.When on-the-spot geostress survey adopts measuring point that different measuring technologies obtain over 15 at different times, consider the workload of terrestrial stress identification, should be taken into account the factors such as point position, the degree of depth, lithology and tomography, selectively adopt actual measurement terrestrial stress data.Should make each subregion and be which has terrestrial stress measuring point in the different degree of depth, and the reliable monitoring result of stress monitoring means as far as possible selectively.

2) determine point position and result

According to selecting point position and quantity, provide measuring point coordinate and terrestrial stress test result.

4, measuring point terrestrial stress calculates result

Adopt the Three-dimension Numerical Model of setting up, according to determined orthogonal numerical analytical plan, analyze successively, carry out altogether Three-dimensional numerical calculation 16 times.For each calculating, extract terrestrial stress monitoring and point out the components of stress, and calculate principle stress and principal direction of stress.

5, terrestrial stress intelligent recognition model is set up

1) set up the relation of measuring point terrestrial stress and protolith parameter and coefficient of horizontal pressure

According to the result of calculation of calculating the terrestrial stress measuring point obtaining, adopt genetic planning to set up j terrestrial stress measuring point (j=1,15) terrestrial stress (size and Orientation) and protolith parameter (density γ, elastic modulus E, Poisson ratio μ) and coefficient of horizontal pressure (x direction coefficient of horizontal pressure λ ₁with y direction coefficient of horizontal pressure λ ₂) between funtcional relationship as follows:

σ' _j1=f _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j2=f _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j3=f _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j1=g _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j2=g _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j3=g _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

2) set up the Optimized model of protolith parameter and coefficient of horizontal pressure

In the plastid of base area, the quadratic sum of the terrestrial stress calculating value of 15 measuring points and the difference of measured value reaches minimum for optimization aim, sets up the Optimized model of the initial stress and coefficient of horizontal pressure.

(1) objective function

$\mathrm{Min}\underset{j=1}{\overset{15}{\mathrm{\&Sigma;}}}(\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{({\mathrm{\&sigma;}}_{\mathrm{ji}}^{\&prime;}-{\mathrm{\&sigma;}}_{\mathrm{ji}})}^2+\underset{i=1}{\overset{3}{\mathrm{\&Sigma;}}}{({\mathrm{\&alpha;}}_{\mathrm{ji}}^{\&prime;}-{\mathrm{\&alpha;}}_{\mathrm{ji}})}^2)$

(2) constraint condition

2.4<γ ₁<2.8，2.8<γ ₂<3.6，2.5<γ ₃<2.9，1.8<γ ₄<2.2

6<E ₁<20，20<E ₂<34，6<E ₃<20，0.8<E ₄<1.3

0.20<μ ₁<0.26，0.20<μ ₂<0.24，0.20<μ ₃<0.26，0.35<μ ₄<0.45

0.8<λ ₁<1.2，1.1<λ ₂<1.5

6, geologic body protolith parameter and coefficient of horizontal pressure Intelligent Recognition

Adopt genetic algorithm to solve the Optimized model of setting up, obtain thus protolith parameter and the coefficient of horizontal pressure of geologic body,

Upper panel protolith parameter:

The protolith parameter of ore body:

Lower panel protolith parameter:

Tomography protolith parameter:

Mining area coefficient of horizontal pressure:

7, geologic body ground stress field inversion

The geologic body protolith parameter that identification obtains according to function

substitution Three-dimension Numerical Model calculates, and thus obtained stress field is initially stress field of mining area.

Only as described above, be only preferred embodiment of the present invention, professional who are familiar with this art such as, and after understanding technological means of the present invention, natural energy, according to actual needs, is changed under instruction of the present invention.Therefore all equal variation and modifications of doing according to the present patent application the scope of the claims, once should still remain within the scope of the patent.

Claims (10)

1. across the recognition methods of scale complex geologic body stress field, comprise the following steps:

The Optimized model that S100, foundation can be identified terrestrial stress;

S200, employing genetic algorithm solve Optimized model, to obtain protolith parameter and the coefficient of horizontal pressure of geologic body;

S300, by the protolith parameter substitution Three-dimension Numerical Model of geologic body, to obtain stress field.

2. according to claim 1ly across the recognition methods of scale complex geologic body stress field, it is characterized in that, in step S100, comprising:

S101, employing genetic planning, according to known terrestrial stress measuring point, set up the funtcional relationship of terrestrial stress monitoring point terrestrial stress and protolith parameter and coefficient of horizontal pressure;

It is optimization aim that the quadratic sum of the terrestrial stress calculating value of S102, base area stress monitoring point and the difference of measured value reaches minimum, and sets up the Optimized model of the initial stress and coefficient of horizontal pressure.

3. according to claim 2 across the recognition methods of scale complex geologic body stress field, it is characterized in that, in step S101, set up j terrestrial stress measuring point (j=1,15) terrestrial stress (size and Orientation) and protolith parameter (density γ, elastic modulus E, Poisson ratio μ) and coefficient of horizontal pressure (x direction coefficient of horizontal pressure λ ₁with y direction coefficient of horizontal pressure λ ₂) between funtcional relationship as follows:

σ' _j1=f _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j2=f _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

σ' _j3=f _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j1=g _1i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j2=g _2i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)

α' _j3=g _3i(γ ₁,E ₁,μ ₁,γ ₂,E ₂,μ ₂,γ ₃,E ₃,μ ₃,γ ₄,E ₄,μ ₄,λ ₁,λ ₂)。

4. according to claim 2ly across the recognition methods of scale complex geologic body stress field, it is characterized in that, in step S102, the objective function of Optimized model is as follows:

Its constraint condition is as follows:

2.4<γ ₁<2.8，2.8<γ ₂<3.6，2.5<γ ₃<2.9，1.8<γ ₄<2.2

6<E ₁<20，20<E ₂<34，6<E ₃<20，0.8<E ₄<1.3

0.20<μ ₁<0.26，0.20<μ ₂<0.24，0.20<μ ₃<0.26，0.35<μ ₄<0.45

0.8<λ ₁<1.2，1.1<λ ₂<1.5?。

5. according to claim 1ly across the recognition methods of scale complex geologic body stress field, it is characterized in that, in step S200, protolith parameter and the coefficient of horizontal pressure of geologic body are as follows:

Upper panel protolith parameter:

The protolith parameter of ore body:

Lower panel protolith parameter:

Tomography protolith parameter:

Mining area is stress coefficient of horizontal pressure initially:

6. according to claim 1ly across the recognition methods of scale complex geologic body stress field, it is characterized in that, in step S300, by obtained geologic body protolith parameter

substitution Three-dimension Numerical Model calculates, thus obtained stress field.

7. implement the claims described in 1 across a device for scale complex geologic body stress field recognition methods, it is characterized in that, comprising:

Optimize module and set up module, for setting up the Optimized model that can identify terrestrial stress;

Genetic algorithm computing module, for Optimized model is solved, to obtain protolith parameter and the coefficient of horizontal pressure of geologic body;

Stress field obtains module, for by the protolith parameter substitution Three-dimension Numerical Model of geologic body, to obtain stress field.

8. device according to claim 7, is characterized in that, in described optimization module, sets up in module, comprising:

Adopt genetic planning, according to known terrestrial stress measuring point, set up the funtcional relationship of terrestrial stress monitoring point terrestrial stress and protolith parameter and coefficient of horizontal pressure;

It is optimization aim that the quadratic sum of the terrestrial stress calculating value of base area stress monitoring point and the difference of measured value reaches minimum, and sets up the Optimized model of the initial stress and coefficient of horizontal pressure.

9. device according to claim 7, is characterized in that, after genetic algorithm computing module is processed, protolith parameter and the coefficient of horizontal pressure of geologic body are as follows:

Upper panel protolith parameter:

The protolith parameter of ore body:

Lower panel protolith parameter:

Tomography protolith parameter:

Mining area is stress coefficient of horizontal pressure initially:

10. device according to claim 7, is characterized in that, by stress field, obtains module, by obtained geologic body protolith parameter

substitution Three-dimension Numerical Model calculates, thus obtained stress field.

Images