CN106768581B - A kind of horizontal stress coefficient of variation measurement method and device - Google Patents
A kind of horizontal stress coefficient of variation measurement method and device Download PDFInfo
- Publication number
- CN106768581B CN106768581B CN201710047366.8A CN201710047366A CN106768581B CN 106768581 B CN106768581 B CN 106768581B CN 201710047366 A CN201710047366 A CN 201710047366A CN 106768581 B CN106768581 B CN 106768581B
- Authority
- CN
- China
- Prior art keywords
- horizontal
- strain
- strain value
- value
- core
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/16—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes for measuring several components of force
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Remote Sensing (AREA)
- Geology (AREA)
- Medicinal Chemistry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The present invention provides a kind of horizontal stress coefficient of variation measurement method and devices, comprising: determines the anelasticity recovery strain value of multiple monitoring positions on columned core curved surface to be measured；It is horizontally and vertically, using the axis for crossing midpoint as vertical pivot, to establish coordinate system to cross two vertical lines at midpoint using the bottom surface midpoint of core to be measured as origin；It determines direction cosines horizontally and vertically of each monitoring position tangent line parallel with the bottom surface of core to be measured relative to coordinate system, determines the proportionate relationship between strain value and the horizontal strain component of core to be measured according to direction cosines；According to above-mentioned strain value and proportionate relationship, horizontal strain component is determined；According to horizontal strain component, horizontal stress coefficient of variation is determined.In the present invention, only by measuring the anelasticity recovery strain value of cylindric core, and parameter is restored according to known anelastic strain, so that it may determine the horizontal stress coefficient of variation of shale reservoir, measuring speed is fast, and it is easy to operate, it is at low cost.
Description
Technical field
The present invention relates to unconventional oil and gas development technique fields, survey in particular to a kind of horizontal stress coefficient of variation
Measure method and device.
Background technique
Unconventional oil and gas development technique is global maximum energy technology revolution nearly ten years, with horizontal well multistage fracturing and
" well factory mode " is that the shale oil-gas exploration development technique of core has promoted the successful exploitation of global shale oil gas.In shale oil
In gas fracturing developing, horizontal stress coefficient of variation is the key that determine that can shale reservoir volume fracturing and realization the economic exploitation ginseng
One of number, horizontal stress coefficient of variation refers to the ratio of two-phase horizontal principal stress difference and minimum horizontal principal stress.
It in the prior art, is that the method combined by small scale fracturing test and shaft lining breakage determines that horizontal stress is poor mostly
Different coefficient, i.e., by small scale fracturing test obtain minimum horizontal principal stress, by hole wall destroy or creep of rock parameter estimation most
The range of big horizontal principal stress, and then determine horizontal stress coefficient of variation, or directly obtain by water pressure supporing
Take horizontal principal stress information.
But horizontal stress coefficient of variation, cost are determined by the method that small scale fracturing test and shaft lining breakage combine
It is higher, and the coefficient can not be obtained before FRACTURING DESIGN and assessment；By the method for water pressure supporing, need to utilize
Drilling platforms carries out, therefore influences construction operation, and needs to be tested in Open-Hole Section, is prone to accidents, and testing cost is high,
Period is long.
Summary of the invention
In view of this, the embodiment of the present invention is designed to provide a kind of horizontal stress coefficient of variation measurement method and dress
It sets, to attempt to solve or alleviate the above problem.
In a first aspect, the embodiment of the invention provides a kind of horizontal stress coefficient of variation measurement methods, which comprises
Determine the anelasticity recovery strain value of multiple monitoring positions on columned core curved surface to be measured；
Using the midpoint of the bottom surface of the core to be measured as origin, using two vertical lines at the excessively described midpoint as horizontal axis
And the longitudinal axis establishes coordinate system using the axis at the excessively described center as vertical pivot；
Determine each monitoring position tangent line parallel with the bottom surface of the core to be measured relative to the coordinate system
Direction cosines horizontally and vertically, determine the anelasticity recovery strain value and the rock to be measured according to the direction cosines
Proportionate relationship between the horizontal strain component of core；
According to the determining anelasticity recovery strain value and the proportionate relationship, the level of the core to be measured is determined
The components of strain；
According to the horizontal strain component, the horizontal maximum strain value and horizontal minimum strain of the core to be measured are determined
Value；
Determine that the level of the core to be measured is answered according to the horizontal maximum strain value and the horizontal minimum strain value
Power coefficient of variation.
With reference to first aspect, the embodiment of the invention provides the first possible implementation of above-mentioned first aspect,
In, the horizontal stress that the core to be measured is determined according to the horizontal maximum strain value and the horizontal minimum strain value
Coefficient of variation, comprising:
Determine that the anelasticity of core to be measured is restored according to the horizontal maximum strain value and the horizontal minimum strain value
Mean strain value；
According to the horizontal maximum strain value and the mean strain value, the horizontal maximum stress of core to be measured is determined
Value；
According to the horizontal minimum strain value and the mean strain value, the horizontal minimum stress of core to be measured is determined
Value；
According to the ratio of the horizontal maximum stress value and the horizontal minimum stress value, the horizontal stress difference is determined
Coefficient.
The possible implementation of with reference to first aspect the first, the embodiment of the invention provides the of above-mentioned first aspect
Two kinds of possible implementations, wherein the ratio according to the horizontal maximum stress value and the horizontal minimum stress value,
Determine the horizontal stress coefficient of variation, comprising:
Calculate the difference between the horizontal maximum stress value and the ratio and 1 of the horizontal minimum stress value；
The difference is determined as the horizontal stress coefficient of variation.
With reference to first aspect, the embodiment of the invention provides the third possible implementation of above-mentioned first aspect,
In, the method also includes:
It is strained according to caused by strain value caused by the horizontal maximum strain value, mineral dehydration and pore pressure loss
Value determines horizontal practical maximum strain value；
It is strained according to caused by strain value caused by the horizontal minimum strain value, mineral dehydration and pore pressure loss
Value determines horizontal practical minimum strain value；
Determine that the horizontal stress is poor according to the practical maximum strain value of the level and the practical minimum strain value of level
Different coefficient.
With reference to first aspect, the embodiment of the invention provides the 4th kind of possible implementation of above-mentioned first aspect,
In, each monitoring position is uniformly distributed at sustained height on the columned core curved surface to be measured.
With reference to first aspect, the embodiment of the invention provides the 5th kind of possible implementation of above-mentioned first aspect,
In, 45 ° are spaced between the adjacent monitoring position.
With reference to first aspect, the embodiment of the invention provides the 6th kind of possible implementation of above-mentioned first aspect,
In, the number of the monitoring position is greater than or equal to 3.
With reference to first aspect, the embodiment of the invention provides the 7th kind of possible implementation of above-mentioned first aspect,
In, the anelasticity recovery strain value of multiple monitoring positions on the columned core curved surface to be measured of determination, comprising:
Each monitoring position is acquired respectively in the corresponding strain value of multiple predetermined times；
Count the number that the corresponding strain value in each monitoring position occurs in multiple predetermined times；
The most strain value of frequency of occurrence is determined as the corresponding anelasticity recovery strain value in the monitoring position.
The 7th kind of possible implementation with reference to first aspect, the embodiment of the invention provides the of above-mentioned first aspect
Eight kinds of possible implementations, wherein be provided with strain gauge at each monitoring position, acquired by the strain gauge each
The monitoring position is in the corresponding strain value of multiple predetermined times.
Second aspect, the embodiment of the invention provides a kind of horizontal stress coefficient of variation measuring devices, wherein described device
Include:
First determining module, for determining that the anelasticity of multiple monitoring positions on columned core curved surface to be measured is restored
Strain value；
Module is established, for using the midpoint of the bottom surface of the core to be measured as origin, with two of the excessively described midpoint
Vertical line, using the axis of the excessively described middle line as vertical pivot, establishes coordinate system as horizontally and vertically；
Second determining module, for determining each monitoring position tangent line parallel with the bottom surface of the core to be measured
Relative to the direction cosines horizontally and vertically of the coordinate system, the anelasticity recovery strain is determined according to the direction cosines
Proportionate relationship between value and the horizontal strain component of the core to be measured；
Third determining module, for determining institute according to the determining anelasticity recovery strain value and the proportionate relationship
State the horizontal strain component of core to be measured；
4th determining module, for determining that the horizontal maximum of the core to be measured is answered according to the horizontal strain component
Variate and horizontal minimum strain value；
5th determining module, for according to the horizontal maximum strain value and the horizontal minimum strain value determine it is described to
Measure the horizontal stress coefficient of variation of core.
Horizontal stress coefficient of variation measurement method provided in an embodiment of the present invention and device, only by measuring cylindric core
Anelasticity recovery strain value, and according to known anelastic strain restore parameter, so that it may determine the water of shale reservoir
Horizontal stress coefficient of variation, measuring speed is fast, easy to operate, at low cost.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment is cited below particularly, and cooperate
Appended attached drawing, is described in detail below.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, and should be solved, and the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as to model
The restriction enclosed for those of ordinary skill in the art without creative efforts, can also be according to these
Attached drawing obtains other relevant attached drawings.
Fig. 1 shows the flow chart of horizontal stress coefficient of variation measurement method provided by the embodiment of the present invention 1；
Fig. 2 shows each monitoring is determined in horizontal stress coefficient of variation measurement method provided by the embodiment of the present invention 1
The flow chart of the corresponding anelasticity recovery strain value in position；
Fig. 3 shows the core to be measured in horizontal stress coefficient of variation measurement method provided by the embodiment of the present invention 1
Schematic diagram；
Fig. 4 shows the structural schematic diagram of horizontal stress coefficient of variation measuring device provided by the embodiment of the present invention 2.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
Middle attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only
It is a part of the embodiment of the present invention, instead of all the embodiments.The present invention being usually described and illustrated herein in the accompanying drawings is real
The component for applying example can be arranged and be designed with a variety of different configurations.Therefore, of the invention to what is provided in the accompanying drawings below
The detailed description of embodiment is not intended to limit the range of claimed invention, but is merely representative of selected reality of the invention
Apply example.Based on the embodiment of the present invention, those skilled in the art institute obtained without making creative work
There are other embodiments, shall fall within the protection scope of the present invention.
In view of in the prior art, being mostly the method measurement level combined by small scale fracturing test and shaft lining breakage
Stress difference coefficient, or horizontal principal stress information is directly obtained by water pressure supporing.But if by small
The method that type mini-frac and shaft lining breakage combine determines horizontal stress coefficient of variation, higher cost, and can not set in pressure break
The coefficient is obtained before meter and assessment；By water pressure supporing method, need to carry out using drilling platforms, therefore influence
Construction operation, and need to be tested in Open-Hole Section, it is prone to accidents, testing cost is high, and the period is long.It is of the invention based on this
Embodiment provides a kind of horizontal stress coefficient of variation measurement method and device, is described below by embodiment.
Embodiment 1
The embodiment of the invention provides a kind of horizontal stress coefficient of variation measurement methods, as shown in Figure 1, this method includes step
Rapid S110-S160, it is specific as follows.
S110 determines the anelasticity recovery strain value of multiple monitoring positions on columned core curved surface to be measured.
Before the horizontal stress coefficient of variation using method provided in an embodiment of the present invention measurement mud shale, scene is needed
The mud shale core for drilling through mud shale core, and drilling through is cylinder, is then carried out the surface of the mud shale core drilled through clear
It washes, obtains core to be measured.
After mud shale core drilling is taken out, at this moment, mud shale core no longer stress is in unloading condition, the present invention
Embodiment exactly restores the corresponding relationship between deformation and horizontal stress according to the anelasticity after core off-load to be measured, comes true
Determine the horizontal stress coefficient of variation of shale reservoir.
Wherein, above-mentioned multiple monitoring positions can be located on curved surface at sustained height, can also be located at different high on curved surface
At degree, the embodiment of the present invention is not to above-mentioned multiple monitoring positions limiting highly specifically on curved surface.
Wherein, above-mentioned each monitoring position is uniformly distributed on columned core curved surface to be measured, i.e., each monitoring position
The angle interval set between adjacent monitoring position is equal, such as, if three monitoring positions are arranged on core to be measured
It sets, then can be and be spaced 120 ° between each monitoring position and its adjacent monitoring position.
Certainly, how above-mentioned only illustrate by taking three monitoring positions as an example is uniformly distributed between monitoring position, not
There is the specific number for limiting the monitoring position being arranged, in addition, can also be with uneven distribution between above-mentioned monitoring position.
In order to more easily calculate the horizontal stress coefficient of variation of core to be measured, it is spaced between adjacent monitoring position
45°。
If being spaced 45 ° between adjacent monitoring position, it at most can be set eight on cylindric core curved surface to be measured
Therefore a monitoring position when measuring the horizontal stress coefficient of variation of core to be measured, can determine answering for eight monitoring positions
The number of variate, the strain value of the monitoring position determined might be less that eight, as long as the level of core to be measured can be measured
Stress difference coefficient, the embodiment of the present invention are not defined the number of the above-mentioned monitoring position used.
Certainly, in order to determining the horizontal stress coefficient of variation of core to be measured, it is thus necessary to determine that core to be measured exists
The number of anelasticity recovery strain value in horizontal multiple directions, i.e., above-mentioned monitoring position is greater than or equal to 3, i.e., in rock to be measured
Three monitoring positions are at least set on core curved surface.
Since the anelasticity recovery strain value of horizontal direction after core off-load to be measured whithin a period of time can be with the time
And change, after a period of time, strain value tends towards stability, and stable strain value is at this moment determined as core to be measured in the direction
Anelasticity recovery strain value, as shown in Fig. 2, determining the anelasticity of multiple monitoring positions on columned core curved surface to be measured
Recovery strain value, including step S210-S230, it is specific as follows.
S210 acquires each monitoring position in the corresponding strain value of multiple predetermined times respectively.
S220 counts the number that the corresponding strain value in each monitoring position occurs in multiple predetermined times.
The most strain value of frequency of occurrence is determined as monitoring the corresponding anelasticity recovery strain value in position by S230.
Wherein, above-mentioned to need to keep the temperature of core to be measured constant when acquiring the strain value of each monitoring position, this
Sample can eliminate the error that core temperature change to be measured is brought.
Wherein, in embodiments of the present invention, the corresponding strain value in each monitoring position in a period of time can be measured, and
This of measurement should be greater than the time that the anelastic strain value of core to be measured tends towards stability for a period of time, by this section of measurement
Period is divided into multiple moment in advance, acquires each monitoring position in the corresponding strain value of predetermined time, when starting one section
In, the corresponding strain value in each monitoring position can vary over, when after a period of time, each monitoring position pair
The strain value answered tends towards stability, i.e., the corresponding strain value in each monitoring position no longer varies over, therefore, can be even
Continue the strain value of the monitoring position of multiple predetermined times measurement be it is equal, i.e., the corresponding strain value in monitoring position can be more
A predetermined time occurs, and therefore, can count the number that the corresponding strain value in each monitoring position occurs in multiple predetermined times,
The most strain value of frequency of occurrence is determined as the corresponding anelasticity recovery strain value in the monitoring position, it is, frequency of occurrence
Most strain values is exactly the stable strain value in the monitoring position.
Wherein, in order to acquire each strain value for monitoring position, it is provided with strain gauge at above-mentioned each monitoring position, is passed through
Strain gauge acquires each monitoring position in the corresponding strain value of multiple predetermined times.
Wherein, above-mentioned strain gauge is attached at the monitoring position of core face to be measured, can at above-mentioned each monitoring position
A strain gauge is arranged, one group of strain gauge also can be set, may include multiple strain gauges in one group of strain gauge, if,
One group of strain gauge is respectively provided on each monitoring position, then the strain gauge for including in the corresponding Strain Meter Set in each monitoring position
Number is equal.
If one group of strain gauge is arranged on each monitoring position, by answering for strain gauge each in this group of strain gauge measurement
Strain value of the average value of variate as the monitoring position.
S120 is made using the central point of the bottom surface of above-mentioned core to be measured as origin with crossing two vertical lines at the midpoint
Horizontally and vertically, using the axis at excessively above-mentioned midpoint as vertical pivot, to establish coordinate system.
Wherein, above-mentioned can be using the central point of the upper bottom surface of cylindrical core to be measured as origin, can also be with cylinder
The central point of the bottom surface of shape core to be measured is as origin, in embodiments of the present invention, with the upper bottom surface of core to be measured
Central point is as origin.
In embodiments of the present invention, the origin of the coordinate system of foundation being indicated that horizontal axis is indicated with X with O, the longitudinal axis is indicated with Y,
Vertical pivot is indicated with Z, establishes the O-XYZ coordinate system of core to be measured.
S130, determine each monitoring position tangent line parallel with the bottom surface of core to be measured relative to the horizontal axis of coordinate system and
The direction cosines of the longitudinal axis determine the horizontal strain point of above-mentioned anelasticity recovery strain value and core to be measured according to direction cosine
Proportionate relationship between amount.
Below will with each monitoring position be arranged a strain gauge, in total be arranged four monitoring positions, and it is adjacent each
It monitors and is spaced for 45 ° between position, the horizontal strain component for how determining above-mentioned strain value Yu core to be measured is discussed in detail
Between proportionate relationship.
As shown in figure 3, aforementioned four monitoring position is denoted as the first monitoring position a respectively1, second monitoring position a2, third prison
Location sets a3With the 4th monitoring position a4, a strain gauge, the first monitoring position a are respectively provided on each monitoring position1On two
A strain is calculated as C1, and the strain value that strain gauge C1 measurement obtains is determined as the first monitoring position a1Corresponding strain value, the
Two monitoring position a2On two strains be calculated as C3, and the strain value that strain gauge C3 measurement obtains is determined as the second monitoring position
a2Corresponding strain value, third monitor position a3On two strains be calculated as C2, and the strain value that strain gauge C2 measurement is obtained is true
It is set to third monitoring position a3Corresponding strain value, the 4th monitoring position a4On two strains be calculated as C4, and strain gauge C4 is surveyed
The strain value measured is determined as the 4th monitoring position a4Corresponding strain value.
Calculate separately the first monitoring position a1, second monitoring position a2, third monitor position a3With the 4th monitoring position a4With
Direction cosines horizontally and vertically of the parallel tangent line in the bottom surface of core to be measured relative to above-mentioned coordinate system, direction cosine refer to
Be cosine value of the above-mentioned tangent line respectively with angle horizontally and vertically, as shown in figure 3, first monitoring position a1With rock to be measured
The tangent line that the bottom surface of core is parallel is parallel with X-axis, and therefore, the angle of the tangent line and X-axis is 0 °, and the angle with Y-axis is 90 °, as a result,
The first monitoring position a can be calculated1It is respectively respectively 1 and 0 with the direction cosines of X-axis and Y-axis, the second monitoring position a2With to
The tangent line that the bottom surface of measurement core is parallel is parallel with Y-axis, and therefore, the angle of the tangent line and X-axis is 90 °, and the angle with Y-axis is
0 °, thus, it is possible to which calculate the second monitoring position is respectively 0 and 1 with the direction cosines of X-axis and Y-axis respectively, third monitors position
a3The angle of the tangent line and X-axis parallel with the bottom surface of core to be measured is 45 °, and the angle with Y-axis is 135 °, thus, it is possible to count
It calculates third and monitors position a3It is respectively respectively 0.707 and -0.707, the 4th monitoring position a with the direction cosines of X-axis and Y-axis4With
The parallel tangent line in the bottom surface of core to be measured and the angle of X-axis are 135 °, and the angle with Y-axis is 45 °, thus, it is possible to calculate
4th monitoring position a4It is respectively respectively -0.707 and 0.707 with the direction cosines of X-axis and Y-axis, above-mentioned first monitoring position a1、
Second monitoring position a2, third monitor position a3With the 4th monitoring position a4Respectively with the direction cosines of X-axis and Y-axis such as 1 institute of table
Show.
Table 1
Wherein, in upper table, direction cosines of the monitoring position tangent line parallel with the bottom surface of core to be measured relative to X-axis
Use liIt indicates, direction cosines m of the parallel tangent line in the bottom surface of monitoring position core to be measured relative to Y-axisiIt indicates, wherein i
=1,2,3,4.When i value is 1, l1Indicate the first monitoring position tangent line parallel with the bottom surface of core to be measured relative to X
The direction cosines of axis, m1Indicate the first monitoring position tangent line parallel with the bottom surface of core to be measured relative to Y-axis direction more than
String.
When determining direction of each monitoring position tangent line parallel with the bottom surface of core to be measured relative to X-axis and Y-axis
After cosine, determine according to the following formula the corresponding strain value in each monitoring position and core to be measured horizontal strain component it
Between proportionality coefficient,
The horizontal component of monitoring position corresponding anelasticity recovery strain value and core to be measured is determined by above-mentioned formula
Between proportionality coefficient after, obtain, the ratio between the corresponding strain value in monitoring position and the horizontal component of core to be measured is closed
System are as follows: b=A ε ', wherein A is the above-mentioned proportionality coefficient being calculated in the formula, and b is that each monitoring position is corresponding stagnant
Elastic recovery strain value, ε ' are the horizontal strain component of core to be measured, and ε '=[ε 'x,ε′y,ε′xy]T, wherein ε 'xFor X-axis
The components of strain in direction, ε 'yFor the components of strain of Y direction, ε 'xyFor shearing strain.
S140 determines the water of core to be measured according to determining above-mentioned anelasticity recovery strain value and aforementioned proportion relationship
The flat components of strain.
The corresponding anelasticity recovery strain in each monitoring position can be determined by the strain gauge on each monitoring position
Value, if a strain gauge is arranged on each monitoring position, the stable strain value which is acquired is determined as this
The corresponding strain value in position is monitored, if one group of strain gauge is provided on each monitoring position, by all strains in the group
The average value for counting the stable strain value of acquisition is determined as the corresponding strain value in monitoring position.
The corresponding strain value b in first monitoring position1It indicates, the corresponding strain value b in the second monitoring position2It indicates, third
Monitor the corresponding strain value b in position3It indicates, the corresponding strain value b in the 4th monitoring position4It indicates, then b=[b1,b2,b3,
b4]TFor the strain value of core to be measured.
According to the formula b=A ε ' and strain value b=[b of determining core to be measured1,b2,b3,b4]T, can be calculated
Horizontal strain component, specific calculating process are as follows:
On the both sides of formula b=A ε ' respectively multiplied by the transposed matrix of the proportionality coefficient between strain value and horizontal strain component,
That is:
ATA ε '=AT b
By solving above-mentioned formula, it can be deduced that ε '=(AT A)-1AT b。
S150 determines the horizontal maximum strain value of core to be measured and horizontal minimum strain according to above-mentioned horizontal strain component
Value.
The horizontal strain component of core to be measured is defined by the above method, is next introduced and how to be answered according to level
Variation, which is measured, determines the horizontal maximum strain value of core to be measured and horizontal minimum strain value.
Strain regime on horizontal plane can indicate that the two principal strains are exactly horizontal maximum by two principal strain values
Strain value and horizontal minimum strain value, calculating horizontal maximum strain value and horizontal minimum strain value by horizontal strain component can be with
It is obtained by solving following equation group:
Wherein, in above-mentioned formula, λ is unknown parameter to be solved, liFor the bottom surface for monitoring position and core to be measured
Direction cosines of the parallel tangent line relative to X-axis, miIt is the monitoring position tangent line parallel with the bottom surface of core to be measured relative to Y
The direction cosines of axis, ε 'xFor the components of strain of X-direction, ε 'yFor the components of strain of Y direction, ε 'xyWith ε 'yxIt is to cut to answer
Become, and ε 'xyWith ε 'yxIt is equal in magnitude, but sign is different.
Above-mentioned equation group is system of homogeneous linear equations, coefficient ranks when the necessary and sufficient conditions of above-mentioned equation group untrivialo solution
Formula is zero, it may be assumed that
It will be quadratic equation with one unknown formula after above-mentioned determinantal expansion, it may be assumed that
(ε′x-λ)(ε′y-λ)-ε′xyε′yx=0
By solving above-mentioned formula, available two root λ about unknown parameter λ1And λ2, by λ1And λ2In it is biggish
One root is determined as the horizontal maximum strain value of core to be measured, by λ1And λ2In a lesser root be determined as core to be measured
Horizontal minimum strain value.
S160 determines horizontal stress coefficient of variation according to horizontal maximum strain value and horizontal minimum strain value.
Wherein, above-mentioned to determine that the horizontal stress of core to be measured is poor according to horizontal maximum strain value and horizontal minimum strain value
Different coefficient, comprising:
Determine that the anelasticity of core to be measured restores mean strain according to horizontal maximum strain value and horizontal minimum strain value
Value；According to the equal strain value of horizontal maximum strain value peace, the horizontal maximum stress value of core to be measured is determined；According to horizontal minimum
The equal strain value of strain value peace, determines the horizontal minimum stress value of core to be measured；Most according to horizontal maximum stress value and level
Ratio between small stress value determines horizontal stress coefficient of variation.
In embodiments of the present invention, after horizontal maximum strain value and horizontal minimum strain value has been determined, according to level
Maximum strain value and horizontal minimum strain value determine the mean strain value of core to be measured by following formula；
εIt is flat=(εmax+εmin)/2
Wherein, in above-mentioned formula, εIt is flatRestore mean strain value, ε for the anelasticity of core to be measuredmaxFor rock to be measured
The horizontal maximum strain value of core, εminFor the horizontal minimum strain value of core to be measured.
Horizontal maximum stress value and level are determined most by following formula according to the equal strain value of horizontal strain value peace
Small stress value；
σi=ei/Jas+em/Jav+P0
Wherein, in above-mentioned formula, the value of i is H and h, indicates horizontal maximum stress value when i is H, the table when i is h
Show horizontal minimum stress value, eiRestore strain, the i.e. difference of horizontal strain value and mean strain value partially, mean strain for anelasticity
Value is the average value of horizontal maximum strain value and horizontal minimum strain value, when i is H, eiFor horizontal maximum strain value and averagely
The difference of strain value, when i is h, eiFor the difference of horizontal maximum strain value and mean strain value, emRestore body for anelasticity to answer
Become, i.e., anelasticity restores mean strain value εIt is flat, Jas is anelasticity recovery strain flexibility, and Jav is that the non-anelastic strain of volume is soft
Degree, P0For pore pressure, σiFor horizontal stress value.
Specifically, determining the water of core to be measured by following formula according to the equal strain value of horizontal maximum strain value peace
Flat maximum stress value are as follows:
σH=(εmax-εIt is flat)/Jas+εIt is flat/Jav+P0
According to the equal strain value of horizontal minimum strain value peace, by following formula, determine that the level of core to be measured is minimum
Stress value are as follows:
σh=(εmin-εIt is flat)/Jas+εIt is flat/Jav+P0
By above-mentioned formula, the ratio between horizontal maximum stress value and horizontal minimum stress value can be determined are as follows:
Since the value of Jav/Jas is certain, the value of Jas/Jav is 1.56 in certain specific embodiments, therefore, horizontal
Ratio between maximum stress value and horizontal minimum stress value can determine are as follows:
Ratio between the horizontal maximum strain value having confirmed and horizontal minimum strain value is substituted into above-mentioned formula, it can
To determine the ratio between horizontal maximum stress value and horizontal minimum stress value, horizontal maximum strain value and horizontal minimum are answered
Ratio between variate is denoted as εmax/εmin=K, therefore, it can be deduced that:
In above formula, P0To be known, Jav and εminKnown to being, therefore, level can be determined most by above-mentioned
Ratio between big stress value and horizontal minimum stress value.
In addition, from above-mentioned formula it can also be seen that ratio between horizontal maximum stress value and horizontal minimum stress value
Increase with the increase of K value.Compare when between the horizontal maximum stress value and horizontal minimum stress value that core to be measured has been determined
After value, according to the ratio between horizontal maximum stress value and horizontal minimum stress value, determine that the horizontal stress of core to be measured is poor
Different coefficient, specifically includes:
Calculate the difference between horizontal maximum stress value and the ratio and 1 of horizontal minimum stress value；Above-mentioned difference is determined
For horizontal stress coefficient of variation.
In embodiments of the present invention, horizontal stress coefficient of variation is between horizontal maximum stress value and horizontal minimum stress value
Difference, with the ratio between horizontal minimum stress value, i.e. DHSR=(σH-σh)/σh, wherein in above-mentioned formula, DHSR is
Horizontal stress coefficient of variation, σHFor horizontal maximum stress value, σhFor horizontal minimum stress value, above-mentioned formula is carried out deformation can be with
It obtains,It therefore, will be between the ratio and 1 between horizontal maximum stress value and horizontal minimum stress value
Difference is determined as horizontal stress coefficient of variation.
When measuring the horizontal stress coefficient of variation of shale reservoir, because mineral dehydration and pore pressure loss also can bands
It strains, when calculating horizontal stress coefficient of variation, it is also necessary to consider strain value caused by mineral dehydration and pore pressure ease
Strain value caused by dissipating, method provided in an embodiment of the present invention further include:
It is true according to strain value caused by strain value caused by horizontal maximum strain value, mineral dehydration and pore pressure loss
Fixed horizontal practical maximum strain value；According to strain value caused by horizontal minimum strain value, mineral dehydration and pore pressure loss
Caused strain value determines horizontal practical minimum strain value；According to horizontal practical maximum strain value and horizontal practical minimum strain value
Determine horizontal stress coefficient of variation.
In step S150, it has been determined that going out horizontal maximum strain value is εmax, horizontal minimum strain value is εmin, will be above-mentioned
Strain value caused by mineral dehydration is denoted as εkw, strain value caused by above-mentioned pore pressure loss is denoted as εkx, therefore, Ke Yitong
Cross formula εmax=εASR max+εkw+εkxCalculate horizontal practical maximum strain value, wherein ε in above-mentioned formulaASRmaxIndicate horizontal
Practical maximum strain value, i.e., maximum strain value caused by core to be measured is discharged because of Original strata stress, according to formula εmin=εASR min
+εkw+εkxCalculate horizontal practical minimum strain value, wherein εASR minIndicate horizontal practical minimum strain value, i.e., core to be measured
Minimum strain value caused by being discharged because of Original strata stress.
As the practical maximum strain value ε of level for defining core to be measuredASRmaxWith horizontal practical minimum strain value
εASR minAfterwards, horizontal practical maximum strain value ε is calculatedASRmaxWith horizontal practical minimum strain value εASR minAverage value εASR is flat,
Middle εASR is flat=(εASR max+εASR min)/2。
Next, according to formulaIt calculates horizontal
The ratio of practical maximum stress value and horizontal practical minimum stress value.
In above-mentioned formula, σH is realFor horizontal practical maximum stress value, σH is realFor horizontal practical minimum stress value, P0For hole
Pressure, Jav are the non-anelastic strain flexibility of volume.
εASRmax、εASR minAnd εASR is flat、P0Be with Jav it is known, therefore, can be determined by above-mentioned formula horizontal real
Ratio between border maximum stress value and horizontal practical minimum stress value.
After the ratio between horizontal practical maximum stress value and horizontal practical minimum stress value has been determined, according toDetermine the horizontal stress coefficient of variation of core to be measured.
Horizontal stress coefficient of variation measurement method provided in an embodiment of the present invention only passes through the stagnant bullet of the cylindric core of measurement
Property recovery strain value, and according to known anelastic strain restore parameter, so that it may determine the horizontal stress of shale reservoir
Coefficient of variation, measuring speed is fast, easy to operate, at low cost.
Embodiment 2
The embodiment of the invention provides a kind of horizontal stress coefficient of variation measuring device, as shown in figure 4, the device includes the
One determining module 410 establishes module 420, the second determining module 430, third determining module 440, the 4th determining module 450 and
Five determining modules 460；
Above-mentioned first determining module 410, for determine on columned core curved surface to be measured it is multiple monitoring positions it is stagnant
Elastic recovery strain value；
It is above-mentioned to establish module 420, for using the midpoint of the bottom surface of core to be measured as origin, to cross two of the midpoint
For vertical line as horizontally and vertically, the axis to cross the midpoint establishes coordinate system as vertical pivot；
Above-mentioned second determining module 430, for determining each monitoring position tangent line parallel with the bottom surface of core to be measured
Relative to the direction cosines horizontally and vertically of above-mentioned coordinate system, according to direction cosines determine above-mentioned anelasticity recovery strain value with
Proportionate relationship between the horizontal component of core to be measured；
Above-mentioned third determining module 440, for determining according to determining above-mentioned anelasticity recovery strain value and proportionate relationship
The horizontal strain component of core to be measured；
Above-mentioned 4th determining module 450, for determining that the level of core to be measured is maximum according to above-mentioned horizontal strain component
Strain value and horizontal minimum strain value；
Above-mentioned 5th determining module 460, for according to above-mentioned horizontal maximum strain value and horizontal minimum strain value determine to
Measure the horizontal stress coefficient of variation of core.
Horizontal stress coefficient of variation measuring device provided in an embodiment of the present invention only passes through the stagnant bullet of the cylindric core of measurement
Property recovery strain value, and according to known anelastic strain restore parameter, so that it may determine the horizontal stress of shale reservoir
Coefficient of variation, measuring speed is fast, easy to operate, at low cost.
Horizontal stress coefficient of variation measuring device provided by the embodiment of the present invention can in equipment specific hardware or
Person is installed on software or firmware in equipment etc..The technology of device provided by the embodiment of the present invention, realization principle and generation
Effect is identical with preceding method embodiment, and to briefly describe, Installation practice part does not refer to place, can refer to preceding method reality
Apply corresponding contents in example.It is apparent to those skilled in the art that for convenience and simplicity of description, foregoing description
System, the specific work process of device and unit, the corresponding process during reference can be made to the above method embodiment, herein no longer
It repeats.
In embodiment provided by the present invention, it should be understood that disclosed device and method, it can be by others side
Formula is realized.The apparatus embodiments described above are merely exemplary, for example, the division of the unit, only one kind are patrolled
Function division is collected, there may be another division manner in actual implementation, in another example, multiple units or components can combine or can
To be integrated into another system, or some features can be ignored or not executed.Another point, shown or discussed is mutual
Coupling, direct-coupling or communication connection can be INDIRECT COUPLING or communication link by some communication interfaces, device or unit
It connects, can be electrical property, mechanical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit
The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple
In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme
's.
In addition, each functional unit in embodiment provided by the invention can integrate in one processing unit, it can also
To be that each unit physically exists alone, can also be integrated in one unit with two or more units.
It, can be with if the function is realized in the form of SFU software functional unit and when sold or used as an independent product
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.
It should also be noted that similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined in a attached drawing, does not then need that it is further defined and explained in subsequent attached drawing, in addition, term " the
One ", " second ", " third " etc. are only used for distinguishing description, are not understood to indicate or imply relative importance.
Finally, it should be noted that embodiment described above, only a specific embodiment of the invention, to illustrate the present invention
Technical solution, rather than its limitations, scope of protection of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, those skilled in the art should understand that: anyone skilled in the art
In the technical scope disclosed by the present invention, it can still modify to technical solution documented by previous embodiment or can be light
It is readily conceivable that variation or equivalent replacement of some of the technical features；And these modifications, variation or replacement, do not make
The essence of corresponding technical solution is detached from the spirit and scope of technical solution of the embodiment of the present invention.Should all it cover in protection of the invention
Within the scope of.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (6)
1. a kind of horizontal stress coefficient of variation measurement method, which is characterized in that the described method includes:
Each preset monitoring position is acquired respectively in the corresponding strain value of multiple predetermined times；
Count the number that the corresponding strain value in each monitoring position occurs in multiple predetermined times；
The most strain value of frequency of occurrence is determined as the corresponding anelasticity recovery strain value in the monitoring position；
Using the midpoint of the bottom surface of core to be measured as origin, using two vertical lines at the excessively described midpoint as horizontally and vertically,
Using the excessively described midpoint and perpendicular to the axis of the horizontal axis and the longitudinal axis as vertical pivot, coordinate system is established；
Determine cross of each monitoring position tangent line parallel with the bottom surface of the core to be measured relative to the coordinate system
The direction cosines of axis and the longitudinal axis determine the anelasticity recovery strain value and the core to be measured according to the direction cosines
Proportionate relationship between horizontal strain component；
According to the determining anelasticity recovery strain value and the proportionate relationship, the horizontal strain of the core to be measured is determined
Component；
It is respectively provided with a strain gauge on each monitoring position, each monitoring is determined by the strain gauge on each monitoring position
The corresponding anelasticity recovery strain value in position is determined according to the determining anelasticity recovery strain value and the proportionate relationship
The horizontal strain component of the core to be measured；
According to the horizontal strain component, two root λ of unknown parameter λ are obtained by formula (1)1And λ2, by λ1And λ2In it is larger
A root be determined as the horizontal maximum strain value of core to be measured, by λ1And λ2In a lesser root be determined as rock to be measured
The horizontal minimum strain value of core；
(ε'x-λ)(ε'y-λ)-ε'xyε'yx=0 (1)
Wherein, ε 'xFor the components of strain of X-direction, ε 'yFor the components of strain of Y direction, ε 'xyWith ε 'yxIt is shearing strain, and
ε'xyWith ε 'yxIt is equal in magnitude, but sign is different；
Determine that the anelasticity of core to be measured restores average according to the horizontal maximum strain value and the horizontal minimum strain value
Strain value；
According to the horizontal maximum strain value and the mean strain value, the horizontal maximum stress value of core to be measured is determined；
According to the horizontal minimum strain value and the mean strain value, the horizontal minimum stress value of core to be measured is determined；
According to the ratio of the horizontal maximum stress value and the horizontal minimum stress value, horizontal stress difference system is determined
Number, wherein calculate the difference between the horizontal maximum stress value and the ratio and 1 of the horizontal minimum stress value；
The difference is determined as the horizontal stress coefficient of variation.
2. the method according to claim 1, wherein the method also includes:
It is true according to strain value caused by strain value caused by the horizontal maximum strain value, mineral dehydration and pore pressure loss
Fixed horizontal practical maximum strain value；
It is true according to strain value caused by strain value caused by the horizontal minimum strain value, mineral dehydration and pore pressure loss
Fixed horizontal practical minimum strain value；
Horizontal stress difference system is determined according to the practical maximum strain value of the level and the practical minimum strain value of level
Number.
3. the method according to claim 1, wherein each monitoring position is columned described to be measured
It is uniformly distributed on core curved surface.
4. the method according to claim 1, wherein being spaced 45 ° between the adjacent monitoring position.
5. the method according to claim 1, wherein the number of the monitoring position is greater than or equal to 3.
6. a kind of horizontal stress coefficient of variation measuring device, which is characterized in that described device includes:
First determining module, for acquiring each preset monitoring position respectively in the corresponding strain value of multiple predetermined times；
Count the number that the corresponding strain value in each monitoring position occurs in multiple predetermined times；
The most strain value of frequency of occurrence is determined as the corresponding anelasticity recovery strain value in the monitoring position；
Module is established, for making with two vertical lines at the excessively described midpoint using the midpoint of the bottom surface of core to be measured as origin
Horizontally and vertically, using the excessively described midpoint and perpendicular to the axis of the horizontal axis and the longitudinal axis as vertical pivot, to establish coordinate system；
Second determining module, for determining that each monitoring position tangent line parallel with the bottom surface of the core to be measured is opposite
In the direction cosines horizontally and vertically of the coordinate system, according to the direction cosines determine the anelasticity recovery strain value with
Proportionate relationship between the horizontal strain component of the core to be measured；
Third determining module, for according to the determining anelasticity recovery strain value and the proportionate relationship, determine it is described to
Measure the horizontal strain component of core；
4th determining module passes through answering on each monitoring position for being respectively provided with a strain gauge on each monitoring position
Become meter and determine the corresponding anelasticity recovery strain value in each monitoring position, according to the determining anelasticity recovery strain value with
The proportionate relationship determines the horizontal strain component of the core to be measured；
According to the horizontal strain component, two root λ of unknown parameter λ are obtained by formula (1)1And λ2, by λ1And λ2In it is larger
A root be determined as the horizontal maximum strain value of core to be measured, by λ1And λ2In a lesser root be determined as rock to be measured
The horizontal minimum strain value of core；
(ε'x-λ)(ε'y-λ)-ε'xyε'yx=0 (1)
Wherein, ε 'xFor the components of strain of X-direction, ε 'yFor the components of strain of Y direction, ε 'xyWith ε 'yxIt is shearing strain, and
ε'xyWith ε 'yxIt is equal in magnitude, but sign is different；
5th determining module, for determining core to be measured according to the horizontal maximum strain value and the horizontal minimum strain value
Anelasticity restore mean strain value；
According to the horizontal maximum strain value and the mean strain value, the horizontal maximum stress value of core to be measured is determined；
According to the horizontal minimum strain value and the mean strain value, the horizontal minimum stress value of core to be measured is determined；
According to the ratio of the horizontal maximum stress value and the horizontal minimum stress value, horizontal stress difference system is determined
Number, wherein calculate the difference between the horizontal maximum stress value and the ratio and 1 of the horizontal minimum stress value；
The difference is determined as the horizontal stress coefficient of variation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710047366.8A CN106768581B (en) | 2017-01-22 | 2017-01-22 | A kind of horizontal stress coefficient of variation measurement method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710047366.8A CN106768581B (en) | 2017-01-22 | 2017-01-22 | A kind of horizontal stress coefficient of variation measurement method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106768581A CN106768581A (en) | 2017-05-31 |
| CN106768581B true CN106768581B (en) | 2019-07-05 |
Family
ID=58945210
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710047366.8A Active CN106768581B (en) | 2017-01-22 | 2017-01-22 | A kind of horizontal stress coefficient of variation measurement method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106768581B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109630098A (en) * | 2019-02-15 | 2019-04-16 | 中国地质科学院地质力学研究所 | A method of poor stress is measured using core |
| CN110243526B (en) * | 2019-07-03 | 2020-11-10 | 东北大学 | Method for determining three-dimensional initial ground stress of rock mass based on flat jack measurement |
| CN112129448A (en) * | 2020-06-07 | 2020-12-25 | 中铁二院工程集团有限责任公司 | Method for measuring ground stress by using elastic recovery deformation of core |
| CN111999165B (en) * | 2020-08-31 | 2021-10-22 | 北京科技大学 | Deep high-stress rock elastic strain recovery monitoring device and method |
| CN111982377B (en) * | 2020-09-01 | 2021-07-23 | 中国地质科学院地质力学研究所 | Deep hole differential stress evaluation method based on secondary sleeve core |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1963423A (en) * | 2006-11-29 | 2007-05-16 | 四川大学 | Method and apparatus for testing stress of cavern wall |
| CN103390112A (en) * | 2013-07-29 | 2013-11-13 | 中国石油大学(华东) | Drilling fluid density design method with consideration of uncertainty of coal seam ground stress |
| CN105317430A (en) * | 2014-07-29 | 2016-02-10 | 中国石油化工股份有限公司 | Method for measuring circumferential stress of well wall surrounding rock of horizontal well in anisotropic formation |
| CN105386756A (en) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | Method for applying dependent variables to calculate porosity of brittle formation |
| CN105784971A (en) * | 2016-03-25 | 2016-07-20 | 大连理工大学 | Experimental facility for utilizing anelastic strain recovery method for in situ analysis of aquo-complex settled layer stress state |
| CN106126866A (en) * | 2016-08-09 | 2016-11-16 | 中国石油天然气股份有限公司 | Sand-proof method for oil-gas well based on geomechanics model and device |
-
2017
- 2017-01-22 CN CN201710047366.8A patent/CN106768581B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1963423A (en) * | 2006-11-29 | 2007-05-16 | 四川大学 | Method and apparatus for testing stress of cavern wall |
| CN103390112A (en) * | 2013-07-29 | 2013-11-13 | 中国石油大学(华东) | Drilling fluid density design method with consideration of uncertainty of coal seam ground stress |
| CN105317430A (en) * | 2014-07-29 | 2016-02-10 | 中国石油化工股份有限公司 | Method for measuring circumferential stress of well wall surrounding rock of horizontal well in anisotropic formation |
| CN105386756A (en) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | Method for applying dependent variables to calculate porosity of brittle formation |
| CN105784971A (en) * | 2016-03-25 | 2016-07-20 | 大连理工大学 | Experimental facility for utilizing anelastic strain recovery method for in situ analysis of aquo-complex settled layer stress state |
| CN106126866A (en) * | 2016-08-09 | 2016-11-16 | 中国石油天然气股份有限公司 | Sand-proof method for oil-gas well based on geomechanics model and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106768581A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106768581B (en) | A kind of horizontal stress coefficient of variation measurement method and device | |
| US5159833A (en) | Analyzing a hydrocarbon reservoir by determining the response of that reservoir to tidal forces | |
| Bonnet et al. | Scaling of fracture systems in geological media | |
| CN109033504B (en) | Oil-water well casing damage prediction method | |
| Hackl et al. | Strain rate patterns from dense GPS networks | |
| CN106709219B (en) | Region initial in-site stress field inversion method and device under complex geological condition | |
| CN108535773A (en) | A kind of micro-seismic event detection method and system | |
| CN106768580B (en) | A kind of horizontal stress measurement method and device | |
| CN112946783A (en) | Method, device and equipment for determining hydrate saturation | |
| CN111815471A (en) | Coal mine resource mining area geological disaster monitoring and early warning management system based on big data | |
| Gross et al. | Rate and state of background stress estimated from the aftershocks of the 1989 Loma Prieta, California, earthquake | |
| CN108180008B (en) | Method and device for predicting collectable reserve, electronic device and storage medium | |
| US20200241157A1 (en) | Real-time surface microseismic monitoring with mobile compact acquisition system | |
| CN110857625B (en) | Method, device and equipment for acquiring reserve of closed water body fracture-cave unit of carbonate rock | |
| CN105354404A (en) | Method for constructing resource quantity parameter distribution model | |
| Heffer et al. | Fracture modelling as part of integrated reservoir characterization | |
| CN115204004A (en) | Three-dimensional ground stress measuring method, device and equipment for carbon dioxide geological sequestration | |
| CN103590813A (en) | Method, tubular column and device for testing dynamic reserves of two-layer commingled producing gas well | |
| CN107014683A (en) | A kind of acquisition methods of earth's surface Rock Displacement Movement parameter | |
| Malthe-Sørenssen et al. | Modeling and characterization of fracture patterns in the Vatnajokull glacier | |
| Dubinya | Analysis of critically stressed fractures behavior at varied depths | |
| CN107036575B (en) | Mine worked-out section deformation and Detection of Stability method based on bouguer gravity anomaly | |
| Liu et al. | Analytical solutions for permeability of a three-dimensional fractal-like tree network model with fractures having variable widths | |
| CN115270062B (en) | Stress relief method for calculating crustal stress by taking irregular drilling hole shape into consideration | |
| CN110851938B (en) | Multi-set reservoir body judgment method and device |
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 |