CN109241209B - Method for rapidly estimating distribution state of ground stress field at any position in China - Google Patents
Method for rapidly estimating distribution state of ground stress field at any position in China Download PDFInfo
- Publication number
- CN109241209B CN109241209B CN201810419188.1A CN201810419188A CN109241209B CN 109241209 B CN109241209 B CN 109241209B CN 201810419188 A CN201810419188 A CN 201810419188A CN 109241209 B CN109241209 B CN 109241209B
- Authority
- CN
- China
- Prior art keywords
- stress
- china
- ground
- magnitude
- distribution state
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for rapidly estimating the distribution state of a ground stress field at any position in China, which comprises the steps of collecting publicly published data related to or reflecting the distribution characteristics of the ground stress field in China; dividing the Chinese whole environment into a plurality of secondary plates according to regional structures in China; inputting the information of 3 aspects into a computer according to the collected data to be used as the basis of quantitative evaluation; assuming that two of the three main stresses are nearly horizontal and one is nearly vertical, and drawing up an estimation method of the magnitude and the direction of the ground stress; according to the estimation method, the computer background calculates and queries the maximum principal stress azimuth and the magnitude of three principal stresses on the basis of relevant actually-measured statistical results according to input parameter conditions and by additionally considering the influence of the earth surface movement speed. According to the present invention, the longitude and latitude of the query position are directly given according to the earth structure characteristics, the earth surface motion and the related stress test results in the existing China, and the state of the earth stress field of the position can be estimated.
Description
Technical Field
The invention relates to a method for quickly estimating the distribution state of a ground stress field at any position in China, belonging to the field of hydrogeology and engineering investigation.
Background
In the rock engineering industry, the distribution state of the initial ground stress field plays a role in the technical links of engineering site selection, underground cavern layout, structural form selection and the like. In the project pre-research stage or on the premise of lacking necessary and available data, if the distribution state of the ground stress field of the project site area can be estimated conveniently, the method has a good guiding effect on the decision of the project.
The existing regional ground stress field estimation method may need engineering geology personnel to consume a large amount of manpower and material resources to look up massive reference documents, collect a large amount of regional geological data, even adopt expensive economic cost to perform some field ground stress test tests, and then summarize the results to estimate the distribution state of the ground stress field of a certain region.
Disclosure of Invention
The ground stress field directly influences the selection of the engineering area and the layout and the model selection of the engineering structure. In view of the characteristics of complicated working process, time consumption and high economic cost of initial crustal stress field distribution of the estimated field area in the past, the invention aims to provide a method for quickly estimating the crustal stress field distribution state of any position in China, collect and integrate the crustal stress test result data in the Chinese disclosed in the past by combining the tectonic features and the surface motion in the China, so that engineering reconnaissance personnel can quickly calculate the corresponding crustal stress distribution state in the early stage of engineering by simply inputting the longitude and latitude of the field area.
The technical scheme of the invention is as follows:
a method for rapidly estimating the distribution state of the ground stress field at any position in China comprises the following steps:
step 1: collecting publicly published data relating to or reflecting the distribution characteristics of the ground stress field in China;
step 2: dividing the Chinese whole environment into a plurality of secondary plates according to regional structures in China;
and step 3: according to the collected data, the following 3 aspects of information are input into a computer as the basis of quantitative evaluation: (1) the magnitude and the direction of the ground surface movement speed; (2) the results of the ground stress test in China; (3) obtaining a statistical relationship according to the distribution of the magnitude of the crustal stress in the rocks with different causes along with the depth;
and 4, step 4: according to the input information, assuming two of the three main stresses are nearly horizontal and one is nearly vertical, and drawing up an estimation method of the magnitude and the direction of the ground stress;
and 5: according to the designed estimation method, according to the input parameter conditions, the computer background takes the relevant actually-measured statistical results as the basis, and simultaneously increases and considers the influence of the ground surface movement speed to calculate and query the maximum principal stress azimuth and the magnitude of three principal stresses;
step 6: when the ground stress of a certain point in China is inquired, the interior of the system is calculated and inquired according to the unit of the tectonic plate of the ground, namely, the source data depended on by the inquiry of the ground stress state of the certain point is only limited in the plate where the certain point is positioned, and the system considers the influence degree of the geographical position of the source data on the inquiry point by adopting the inverse square of the distance.
According to the present invention, the longitude and latitude of the query position are directly given according to the earth structure characteristics, the earth surface motion and the related stress test results in the existing China, and the state of the earth stress field of the position can be estimated. The method is flexible, simple and convenient, can be used for determining the building layout at the initial stage of engineering investigation, can adjust and correct the estimation method according to partial data mastered by a user hand so as to estimate and obtain a more reasonable ground stress field, and is easy to master the analysis method.
Drawings
FIG. 1 illustrates source data for crustal motion velocity in different regions;
FIG. 2 illustrates source data for maximum horizontal stress in different regions;
FIG. 3 illustrates secondary panels to be scribed by regional structures within China;
FIG. 4 shows a diagram of stress query at a location in China.
Detailed Description
The invention is further illustrated by the following figures and examples, which are not to be construed as limiting the invention in any way.
Referring to fig. 1 to 4, the specific implementation steps of the present invention are as follows:
step 1: collecting publicly published data relating to or reflecting the distribution characteristics of the ground stress field in China, comprising: (1) the GPS monitoring results of the earth's surface motion, the nature of which indicates the tectonic plate interaction characteristics under present-day conditions, wherein the motion displacement vector direction generally indicates the direction of maximum principal stress, while the magnitude of velocity is closely related to the degree of tectonic squeezing, i.e. the magnitude of maximum principal stress. Or, the influence of the tectonic plate on the position and the size of the maximum principal stress in the ground stress can be reflected by the position and the size of the ground surface movement speed, so that the quantitative relation between the ground stress state and the ground structure is established, and the method is shown in figure 1; (2) the results of the ground stress test and the seismic source interpretation, the 'world ground stress' research developed by austria organization aims at 800 groups of ground stress test data obtained in China, and refer to fig. 2; (3) zhu Huan Chun et al obtained the relation between ground stress size and rock cause according to the statistics of 300 multiunit actual measurement stress achievements, can help further analysis ground stress state.
Step 2: the whole situation of China is divided into a plurality of secondary plates according to regional structures in China, and the ground stress fields between the secondary plates are connected with each other but are not continuous overall, so that when analyzing regional ground stress field characteristics, the effect of the plates needs to be considered first, and the boundary between adjacent plates is used as the boundary of the ground stress field boundary, see fig. 3.
And step 3: according to the collected data, the following 3 aspects of information are implanted into a computer as the basis of quantitative evaluation: (1) the magnitude and the direction of the ground surface movement speed; (2)800 more groups of Chinese domestic internal ground stress test results (including seismic source mechanism interpretation results); (3) and obtaining statistical relationship along with depth distribution according to 300 groups of crustal stress in rocks with different causes.
And 4, step 4: according to the implantation information, assuming that two of the three main stresses are nearly horizontal and one is nearly vertical, the method for estimating the magnitude and the direction of the ground stress is formulated, and the specific implementation mode is as follows: (1) the maximum horizontal principal stress direction is a weighted average value of the earth surface movement direction and the test result, the speed direction is more representative according to the prior engineering practice experience, and higher weight can be given when the device is used; (2) the three main stress magnitudes are weighted averages of the velocity magnitude and the lithology respectively in consideration of the influences of the velocity magnitude and the lithology.
And 5: according to the designed estimation method, according to the input parameter conditions, the program background takes the relevant actually-measured statistical results as the basis, and simultaneously, the influence of the earth surface movement speed is additionally considered, and the calculation query of the maximum principal stress azimuth and the magnitude of the three principal stresses is carried out. The input parameters comprise longitude and latitude of the query position, rock types (sedimentary rock and negative metamorphic rock, magmatic rock and positive metamorphic rock), uniaxial compressive strength of rock, volume weight of rock mass, surface motion for calculating the maximum principal stress azimuth, and influence weight of test results.
Step 6: when the ground stress of a certain point in China is inquired, the interior of the system is calculated and inquired according to the unit of the tectonic plate of the ground, namely, the source data depended on by the inquiry of the ground stress state of the certain point is only limited in the plate where the certain point is positioned, and the system considers the influence degree of the geographical position of the source data on the inquiry point by adopting the inverse square of the distance.
And 7: query of maximum horizontal principal stress orientation: for any position to be checked, the longitude and latitude of the position are given, the system searches data stored in a background in a plate block according to the distance relation, respectively obtains the earth surface motion direction (indicating the maximum horizontal principal stress direction) of the point and the direction indicated by the measured value, and then calculates the comprehensive result according to a weighting mode. In specific application, when the weight of one factor is zero, the factor is ignored, and only the index value with the weight of 1 is used as the calculation basis.
And 8: three principal stress magnitude queries: according to the data collected in the step 1, the system calculates and provides three distribution characteristics of the magnitude of the main stress according to the input parameter conditions. The system assumes that three principal stresses change linearly with depth, the slope is K, and the intercept is T, namely the principal stress at any depth in the rock mass obeys the definition of the formula (1).
σ＝K·h+T (1)
Where σ and h are a certain principal stress and buried depth, respectively.
And step 9: and (4) calculating according to the formula (2) by taking the maximum horizontal principal stress azimuth at the query point as the weighted value of the earth surface motion azimuth and the azimuth indicated by the measured value at the point according to the step (7).
In the formula (I), the compound is shown in the specification,
the maximum horizontal principal stress azimuth vector is obtained through interpolation calculation according to background earth surface motion and test results, and a and b are corresponding weight coefficients.
Step 10: the distribution rule of the three main stresses according to the step 8 is obeyed to a linear expression with the slope of K and the intercept of T, and K, T is calculated according to the expressions (3) and (4) respectively. The values of the parameters K, T take into account the influence of rock causes and surface motion, and the system takes into account the simplified conditions of two rock causes, namely magma rock, positive metamorphic rock, sedimentary rock and negative metamorphic rock.
K＝A·exp(B·v) (3)
T＝a×ln(v)+b (4)
In the formula, v is the speed of the appointed query part obtained through background earth surface motion data interpolation calculation, and A, a and b are known coefficients and are taken according to the value in the table 1.
Step 11: according to the description of the step 8, the rock mechanical property can influence the maximum principal stress, so that the principal stress distribution parameter K is corrected according to the input rock uniaxial compressive strength according to the formula (4). In estimating the magnitude of the geostress, the system defaults to using the results of step 10, allowing the user to consider the effect of rock strength.
K＝(0.4σc-31)/10000 (4)
In the formula, σcThe uniaxial compressive strength of the rock is in the value range of [25,150]When the average molecular weight is less than 25, 25 is taken, and when the average molecular weight is more than 150, 150 is taken.
Table 1: calculation coefficient value of principal stress magnitude
Wherein σH-maximum horizontal stress; sigmah-minimum horizontal stress; sigmav-vertical stress.
Example 1:
as shown in FIG. 4, in a large-scale hydraulic and hydroelectric project in the southwest area of China, the distribution of the ground stress field in the field area needs to be estimated, and the specific implementation steps are as follows:
1. obtaining rough position information of the area according to the topographic map, and obtaining the longitude and latitude of the position: east longitude 101 deg. 51 'and north latitude 25 deg. 20'.
2. According to the collected past engineering data, the main lithology type is sandstone (according to claim 9, the rock types are 'sedimentary rock and negative metamorphic rock'), the uniaxial compressive strength is 80MPa, and the volume weight of the rock is 2.5KN/m3。
3. Since the data collected earlier is less than perfect, the weights of the measured geostress value and the geostress velocity are assumed to be 0.5 and 0.5, respectively.
4. Inputting the information into a ground stress query system to obtain the ground stress field distribution of the position:
σH＝0.029h+5.591
σv＝0.026h
σh＝0.026h+3.347
maximum horizontal principal stress σHIs N43.5 ° E. In the formula, h is the buried depth.
The above are only specific application examples of the present invention, and other embodiments of the present invention are within the scope of the present invention as claimed by using equivalent alternatives or equivalent variations.
Claims (10)
1. A method for rapidly estimating the distribution state of the ground stress field at any position in China is characterized by comprising the following steps:
step 1: collecting publicly published data relating to or reflecting the distribution characteristics of the ground stress field in China;
step 2: dividing the Chinese whole environment into a plurality of secondary plates according to regional structures in China;
and step 3: according to the collected data, the following 3 aspects of information are input into a computer as the basis of quantitative evaluation: (1) the magnitude and the direction of the ground surface movement speed; (2) the results of the ground stress test in China; (3) obtaining a statistical relationship according to the distribution of the magnitude of the crustal stress in the rocks with different causes along with the depth;
and 4, step 4: according to the input information, assuming two of the three main stresses are nearly horizontal and one is nearly vertical, and drawing up an estimation method of the magnitude and the direction of the ground stress;
and 5: according to the designed estimation method, according to the input parameter conditions, the computer background takes the relevant actually-measured statistical results as the basis, and simultaneously increases and considers the influence of the ground surface movement speed to calculate and query the maximum principal stress azimuth and the magnitude of three principal stresses;
step 6: when the ground stress of a certain point in China is inquired, the interior of the system is calculated and inquired according to the unit of the tectonic plate of the ground, namely, the source data depended on by the inquiry of the ground stress state of the certain point is only limited in the plate where the certain point is positioned, and the system considers the influence degree of the geographical position of the source data on the inquiry point by adopting the inverse square of the distance.
2. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 1, wherein: the query for the maximum horizontal principal stress orientation is such that: and giving longitude and latitude of any position to be checked, searching data stored in a background in the board block by the system according to the distance relation, respectively obtaining the earth surface motion direction of the point and the direction indicated by the measured value, and calculating a comprehensive result according to a weighting mode.
3. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 2, wherein: the query for the magnitude of the three principal stresses is such that: according to the data collected in the step 1 and the input parameter conditions, the system calculates and provides the distribution characteristics of the three main stress magnitudes; the system assumes that three principal stresses are linearly changed along with the depth, the slope is K, and the intercept is T, namely the principal stress at any depth in the rock mass obeys the definition of the formula (1);
σ＝K·h+T (1)
where σ and h are a certain principal stress and buried depth, respectively.
4. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 2, wherein: the maximum horizontal principal stress position at the query point is the weighted value of the earth surface motion position of the point and the position indicated by the measured value, and is obtained by calculation according to the formula (2);
in the formula (I), the compound is shown in the specification,
the maximum horizontal principal stress azimuth vector is obtained through interpolation calculation according to background earth surface motion and test results, and a and b are corresponding weight coefficients.
5. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 3, wherein: the distribution rule of the three main stresses obeys a linear expression with the slope of K and the intercept of T, and K, T is calculated according to the formulas (3) and (4) respectively;
K＝A·exp(B·v) (3)
T＝a·ln(v)+b (4)
in the formula, v is the speed of the appointed query part obtained through background earth surface motion data interpolation calculation, and A, a and b are known coefficients and are taken by table lookup.
6. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 3, wherein: correcting the main stress distribution parameter K according to the input uniaxial compressive strength of the rock and the formula (4);
K＝(0.4σc-31)/10000 (4)
in the formula, σcThe uniaxial compressive strength of the rock is in the value range of [25,150]When the average molecular weight is less than 25, 25 is taken, and when the average molecular weight is more than 150, 150 is taken.
7. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 1, wherein: the data of step 1 includes: (1) GPS monitoring results of surface motion; (2) ground stress test and seismic source interpretation results; (3) and (3) carrying out statistics on relationship data between the obtained ground stress magnitude and the rock cause according to 300 groups of measured stress results.
8. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 1, wherein: the specific implementation manner of the step 4 is as follows: (1) the maximum horizontal principal stress direction is the weighted average of the earth surface motion direction and the test result; (2) the three main stress magnitudes are weighted averages of the velocity magnitude and the lithology respectively in consideration of the influences of the velocity magnitude and the lithology.
9. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 1, wherein: and 5, inputting parameters including longitude and latitude of the inquiry position, rock type, rock uniaxial compressive strength and rock mass volume weight, and surface motion and test results for calculating the maximum principal stress azimuth.
10. The method for rapidly estimating the distribution state of the geostress field at any position in China according to claim 2, wherein: when the query of the maximum horizontal principal stress azimuth is applied specifically, when the weight of one factor is zero, the factor is ignored, and only the index value with the weight of 1 is used as the calculation basis.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810419188.1A CN109241209B (en) | 2018-05-04 | 2018-05-04 | Method for rapidly estimating distribution state of ground stress field at any position in China |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810419188.1A CN109241209B (en) | 2018-05-04 | 2018-05-04 | Method for rapidly estimating distribution state of ground stress field at any position in China |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109241209A CN109241209A (en) | 2019-01-18 |
| CN109241209B true CN109241209B (en) | 2021-08-27 |
Family
ID=65084103
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810419188.1A Active CN109241209B (en) | 2018-05-04 | 2018-05-04 | Method for rapidly estimating distribution state of ground stress field at any position in China |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109241209B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107489424A (en) * | 2017-07-28 | 2017-12-19 | 西安理工大学 | A kind of shield subway work induces stratum deformation and the predictor method influenceed on ancient building |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7882745B2 (en) * | 2006-09-20 | 2011-02-08 | Schlumberger Technology Corporation | Method and system to invert tectonic boundary or rock mass field in in-situ stress computation |
-
2018
- 2018-05-04 CN CN201810419188.1A patent/CN109241209B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107489424A (en) * | 2017-07-28 | 2017-12-19 | 西安理工大学 | A kind of shield subway work induces stratum deformation and the predictor method influenceed on ancient building |
Non-Patent Citations (2)
| Title |
|---|
| "Application of Three-Dimensional Numerical Calculation to Geo-stress Field Analysis of Oil or Gas Reservoir";Hui Deng 等;《2010 International Conference on Computational and Information Sciences》;20101219;第971-974页 * |
| "复杂地应力场反演优化算法及研究新进展";王金安 等;《中国矿业大学学报》;20150312;第189-205页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109241209A (en) | 2019-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hinze et al. | New standards for reducing gravity data: The North American gravity database | |
| Plomerová et al. | Long memory of mantle lithosphere fabric—European LAB constrained from seismic anisotropy | |
| CN100580679C (en) | Civil construction engineering GPS-RTK measurement integrated PDA device | |
| Guo et al. | Generation and verification of three-dimensional network of fractured rock masses stochastic discontinuities based on digitalization | |
| CN105426619A (en) | Analysis method for determining longitudinal bending rigidity of shield tunnel | |
| CN106094022A (en) | A kind of soil layer mean shear wave velocity and the measuring method of soil layer reflecting interface buried depth | |
| CN108824409A (en) | Railway high precipitous rock slope micro seismic monitoring method and system | |
| Xia et al. | Slope stability analysis based on group decision theory and fuzzy comprehensive evaluation | |
| CN102901985A (en) | Depth domain layer speed correcting method suitable for undulating surface | |
| CN104142516A (en) | Method for predicting thickness of thin single sand bed | |
| CN109459787B (en) | coal mine underground structure imaging method and system based on seismic channel wave full-waveform inversion | |
| CN105242328A (en) | Paleo-thermal lithosphere thickness determination method and device | |
| CN114662360A (en) | Inversion method for constructing complex geological three-dimensional ground stress field by using multi-software combination | |
| Li et al. | A quantitative model for the geological strength index based on attribute mathematics and its application | |
| CN109490978A (en) | A kind of frequency domain quick high accuracy forward modeling method on fluctuating stratum | |
| Guo et al. | Topography-dependent eikonal tomography based on the fast-sweeping scheme and the adjoint-state technique | |
| CN103364823A (en) | Vibration source real-time positioning and analyzing system | |
| CN104123455A (en) | Three-dimensional dip inversion method for nonlinear conjugate gradient of magnetotelluric field | |
| CN109241209B (en) | Method for rapidly estimating distribution state of ground stress field at any position in China | |
| Plomerová et al. | Seismic anisotropy—a key to resolve fabrics of mantle lithosphere of Fennoscandia | |
| CN103278852B (en) | Utilize the method for seismic data volume waveform configuration characteristic model predicting oil/gas | |
| CN115586577A (en) | Full-time conversion method for fixed-source transient electromagnetic non-central-point observation data | |
| Hammami et al. | Aquifer structure identification through geostatistical integration of geological parameters: case of the Triassic sandstone aquifer system (SE Tunisia) | |
| CN110008302A (en) | More impact factor precipitation recharge coefficient quantitative evaluation methods | |
| CN113536693B (en) | Aviation-ground-well magnetic anomaly data joint inversion method based on well rock physical property constraint |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |