CN106768581A - A kind of horizontal stress coefficient of variation measuring method and device - Google Patents
A kind of horizontal stress coefficient of variation measuring method and device Download PDFInfo
- Publication number
- CN106768581A CN106768581A CN201710047366.8A CN201710047366A CN106768581A CN 106768581 A CN106768581 A CN 106768581A CN 201710047366 A CN201710047366 A CN 201710047366A CN 106768581 A CN106768581 A CN 106768581A
- Authority
- CN
- China
- Prior art keywords
- horizontal
- strain value
- core
- value
- measured
- 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.)
- Granted
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 invention provides a kind of horizontal stress coefficient of variation measuring method and device, including：Determine the anelasticity recovery strain value of multiple monitoring positions on columned core curved surface to be measured；As origin, with two vertical lines at midpoint excessively as horizontally and vertically, the axis with midpoint excessively is vertical pivot, sets up coordinate system at bottom surface midpoint with core to be measured；Determine each monitoring position tangent line direction cosines horizontally and vertically relative to coordinate system parallel with the bottom surface of core to be measured, the proportionate relationship between strain value and the horizontal strain component of core to be measured is determined 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 recovered according to known anelastic strain, it is possible to determine the horizontal stress coefficient of variation of shale reservoir, measuring speed is fast, simple to operate, low cost.
Description
Technical field
The present invention relates to unconventionaloil pool development technique field, surveyed in particular to a kind of horizontal stress coefficient of variation
Amount method and device.
Background technology
Unconventionaloil pool development technique is the maximum energy technology revolution in the whole world nearly ten years, with horizontal well multistage fracturing and
" well factory mode " has promoted the successful exploitation of global shale oil gas for the shale oil-gas exploration development technique of core.In shale oil
In gas fracturing developing, horizontal stress coefficient of variation is to determine that can shale reservoir volume fracturing and the crucial ginseng for realizing the economic exploitation
One of number, horizontal stress coefficient of variation refers to the ratio of two-phase horizontal principal stress difference and minimum horizontal principal stress.
In the prior art, it is mostly that the method being combined by small scale fracturing test and shaft lining breakage determines that horizontal stress is poor
Different coefficient, i.e., obtain minimum horizontal principal stress by small scale fracturing test, by hole wall destruction or creep of rock parameter estimation most
The scope of big horizontal principal stress, and then determine horizontal stress coefficient of variation, or directly obtained by water pressure supporing
Take horizontal principal stress information.
But, the method being combined by small scale fracturing test and shaft lining breakage determines horizontal stress coefficient of variation, cost
It is higher, and the coefficient cannot be obtained before FRACTURING DESIGN and assessment；By the method for water pressure supporing, it is necessary to utilize
Drilling platforms is carried out, therefore influence construction operation, and needs are tested in Open-Hole Section, are prone to accidents, and testing cost is high,
Cycle is long.
The content of the invention
In view of this, the purpose of the embodiment of the present invention is to provide a kind of horizontal stress coefficient of variation measuring method and dress
Put, to attempt to solve or alleviate above mentioned problem.
In a first aspect, the embodiment of the invention provides a kind of horizontal stress coefficient of variation measuring method, methods described includes：
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 excessively described midpoint as transverse axis
And the longitudinal axis, the axis using excessively described center sets up coordinate system 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 with the rock to be measured according to the direction cosines
Proportionate relationship between the horizontal strain component of core；
According to the anelasticity recovery strain value and the proportionate relationship that determine, 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；
The level for determining the core to be measured according to the horizontal maximum strain value and the horizontal minimum strain value should
Power coefficient of variation.
With reference in a first aspect, the embodiment of the invention provides the first possible implementation of above-mentioned first aspect, its
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, including：
The anelasticity for determining core to be measured according to the horizontal maximum strain value and the horizontal minimum strain value is recovered
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 horizontal maximum stress value and the ratio of the horizontal minimum stress value, the horizontal stress difference is determined
Coefficient.
With reference to the first possible implementation of first aspect, the of above-mentioned first aspect is the embodiment of the invention provides
Two kinds of possible implementations, wherein, it is described according to the horizontal maximum stress value and the ratio of the horizontal minimum stress value,
Determine the horizontal stress coefficient of variation, including：
Calculate the difference between the ratio and 1 of the horizontal maximum stress value and the horizontal minimum stress value；
The difference is defined as the horizontal stress coefficient of variation.
With reference in a first aspect, the embodiment of the invention provides the third possible implementation of above-mentioned first aspect, its
In, methods described also includes：
Strained according to caused by strain value caused by the horizontal maximum strain value, mineral dehydration and pore pressure loss
The value actual maximum strain value of determination level；
Strained according to caused by strain value caused by the horizontal minimum strain value, mineral dehydration and pore pressure loss
The value actual minimum strain value of determination level；
Determine that the horizontal stress is poor according to the actual maximum strain value of the level and the actual minimum strain value of the level
Different coefficient.
With reference in a first aspect, the embodiment of the invention provides the 4th kind of possible implementation of above-mentioned first aspect, its
In, each described monitoring position is uniformly distributed at sustained height on the columned core curved surface to be measured.
With reference in a first aspect, the embodiment of the invention provides the 5th kind of possible implementation of above-mentioned first aspect, its
In, 45 ° are spaced between the adjacent monitoring position.
With reference in a first aspect, the embodiment of the invention provides the 6th kind of possible implementation of above-mentioned first aspect, its
In, the number of the monitoring position is more than or equal to 3.
With reference in a first aspect, the embodiment of the invention provides the 7th kind of possible implementation of above-mentioned first aspect, its
In, the anelasticity recovery strain value for determining multiple monitoring positions on columned core curved surface to be measured, including：
Each described monitoring position is gathered respectively in the corresponding strain value of multiple predetermined times；
Count the number of times that each corresponding strain value in monitoring position occurs in multiple predetermined times；
Will appear from the most strain value of number of times and be defined as the corresponding anelasticity recovery strain value in the monitoring position.
With reference to the 7th kind of possible implementation of first aspect, the of above-mentioned first aspect is the embodiment of the invention provides
Eight kinds of possible implementations, wherein, strain gauge is provided with each described monitoring position, each is gathered by the strain gauge
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 measurement apparatus, wherein, described device
Including：
First determining module, the anelasticity for determining multiple monitoring positions on columned core curved surface to be measured is recovered
Strain value；
Module is set up, for using the midpoint of the bottom surface of the core to be measured as origin, with two of excessively described midpoint
Used as horizontally and vertically, the axis using excessively described center line sets up coordinate system to vertical line as vertical pivot；
Second determining module, the tangent line parallel with the bottom surface of the core to be measured for determining each monitoring position
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；
3rd determining module, for according to the anelasticity recovery strain value and the proportionate relationship for determining, determining institute
State the horizontal strain component of core to be measured；
4th determining module, for according to the horizontal strain component, determining that the level of the core to be measured is maximum and answering
Variate and horizontal minimum strain value；
5th determining module, for being treated according to the horizontal maximum strain value and the horizontal minimum strain value determination
Measure the horizontal stress coefficient of variation of core.
Horizontal stress coefficient of variation measuring 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 recover parameter, it is possible to determine the water of shale reservoir
Horizontal stress coefficient of variation, measuring speed is fast, simple to operate, low cost.
To enable the above objects, features and advantages of the present invention to become apparent, preferred embodiment cited below particularly, and coordinate
Appended accompanying drawing, is described in detail below.
Brief description of the drawings
Technical scheme in order to illustrate more clearly the embodiments of the present invention, below will be attached to what is used needed for embodiment
Figure is briefly described, and should solve, and the following drawings illustrate only certain embodiments of the present invention, therefore be not construed as to model
The restriction enclosed, for those of ordinary skill in the art, on the premise of not paying creative work, can also be according to these
Accompanying drawing obtains other related accompanying drawings.
Fig. 1 shows the flow chart of the horizontal stress coefficient of variation measuring method that the embodiment of the present invention 1 is provided；
Fig. 2 determines that each is monitored in showing the horizontal stress coefficient of variation measuring method that the embodiment of the present invention 1 is provided
The flow chart of the corresponding anelasticity recovery strain value in position；
Fig. 3 shows the core to be measured in the horizontal stress coefficient of variation measuring method that the embodiment of the present invention 1 is provided
Schematic diagram；
Fig. 4 shows the structural representation of the horizontal stress coefficient of variation measurement apparatus that the embodiment of the present invention 2 is provided.
Specific embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
Middle accompanying drawing, is clearly and completely described to the technical scheme in the embodiment of the present invention, it is clear that described embodiment is only
It is a part of embodiment of the invention, rather than whole embodiments.The present invention generally described and illustrated in accompanying drawing herein is real
The component for applying example can be arranged and designed with a variety of configurations.Therefore, it is of the invention to what is provided in the accompanying drawings below
The detailed description of embodiment is not intended to limit the scope of claimed invention, but is merely representative of selected reality of the invention
Apply example.Based on embodiments of the invention, the institute that those skilled in the art are obtained on the premise of creative work is not made
There is other embodiment, belong to the scope of protection of the invention.
In view of in the prior art, being mostly the method measurement level that is combined by small scale fracturing test and shaft lining breakage
Stress difference coefficient, or directly by water pressure supporing acquisition horizontal principal stress information.But, if by small
The method that type mini-frac and shaft lining breakage are combined determines horizontal stress coefficient of variation, relatively costly, and cannot be set in pressure break
The coefficient is obtained before meter and assessment；By water pressure supporing method, it is necessary to be carried 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 cycle is long.It is of the invention based on this
Embodiment provides a kind of horizontal stress coefficient of variation measuring method and device, is described below by embodiment.
Embodiment 1
A kind of horizontal stress coefficient of variation measuring method is the embodiment of the invention provides, as shown in figure 1, the 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.
, it is necessary to scene before using the horizontal stress coefficient of variation of method provided in an embodiment of the present invention measurement mud shale
It is cylinder to drill through mud shale core, and the mud shale core for drilling through, and the surface of the mud shale core that then will be drilled through is carried out clearly
Wash, obtain core to be measured.
After mud shale core drilling is taken out, at this moment, mud shale core no longer stress, i.e., in unloading condition, the present invention
Embodiment, exactly recovers the corresponding relation 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 position may be located on curved surface at sustained height, it is also possible to 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
Put and the angle interval between adjacent monitoring position is equal, such as, if setting three monitoring positions on core to be measured
Put, then can be that 120 ° are spaced between each monitoring position and its adjacent monitoring position.
Certainly, it is above-mentioned simply by three monitor position as a example by illustrate monitoring position between how to be uniformly distributed, not
There is the specific number for limiting the monitoring position for setting, in addition, can also 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, eight can be at most set on cylindric core curved surface to be measured
Individual monitoring position, therefore, when the horizontal stress coefficient of variation of core to be measured is measured, it may be determined that eight are monitored answering for position
Variate, it is determined that the number of strain value of monitoring position 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, not the number to the above-mentioned monitoring position for using be defined.
Certainly, in order to determine the horizontal stress coefficient of variation of core to be measured, it is thus necessary to determine that core to be measured exists
Anelasticity recovery strain value in horizontal multiple directions, i.e., the number of above-mentioned monitoring position is more than or equal to 3, i.e., in rock to be measured
At least provided with three monitoring positions on core curved surface.
Because the anelasticity recovery strain value of horizontal direction after core off-load to be measured can over time within a period of time
And change, after a period of time, strain value is tended towards stability, and the strain value of stabilization at this moment is defined as into 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, gathers each monitoring position in the corresponding strain value of multiple predetermined times respectively.
S220, counts the number of times that each corresponding strain value in monitoring position occurs in multiple predetermined times.
S230, will appear from the most strain value of number of times and is defined as the corresponding anelasticity recovery strain value in monitoring position.
Wherein, it is above-mentioned when the strain value of each monitoring position is gathered, it is necessary to keep the temperature of core to be measured constant, this
Sample can eliminate the error that core temperature change to be measured is brought.
Wherein, in embodiments of the present invention, each corresponding strain value in monitoring position in a period of time can be measured, and
This of measurement should be more than the time that the anelastic strain value of core to be measured tends towards stability, this section that will be measured for a period of time
Time period is divided into multiple moment in advance, each monitoring position is gathered in the corresponding strain value of predetermined time, when starting one section
Interior, each corresponding strain value in monitoring position can be varied over, when after a period of time, each monitoring position pair
The strain value answered is tended towards stability, i.e., each corresponding strain value in monitoring position is no longer varied over, therefore, can be even
The strain value of the monitoring position of continuous multiple predetermined time measurement is equal, i.e., the corresponding strain value in monitoring position can be more
Individual predetermined time occurs, and therefore, it can count the number of times that each corresponding strain value in monitoring position occurs in multiple predetermined times,
Will appear from the most strain value of number of times and be defined as the corresponding anelasticity recovery strain value in the monitoring position, it is, occurrence number
Most strain values is exactly the strain value of monitoring position stabilization.
Wherein, in order to gather each strain value for monitoring position, strain gauge is provided with above-mentioned each monitoring position, is passed through
Strain gauge gathers 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
To set a strain gauge, it is also possible to set one group of strain gauge, multiple strain gauges can be included in one group of strain gauge, if,
Be respectively provided with one group of strain gauge on each monitoring position, then the strain gauge for including in each corresponding Strain Meter Set in monitoring position
Number is equal.
If setting one group of strain gauge on each monitoring position, by answering that each strain gauge in this group of strain gauge is measured
The average value of variate as the monitoring position strain value.
S120, using the central point of the bottom surface of above-mentioned core to be measured as origin, is made with two vertical lines for crossing the midpoint
For horizontally and vertically, the axis using excessively above-mentioned midpoint sets up coordinate system as vertical pivot.
Wherein, it is above-mentioned can be using the central point of the upper bottom surface of cylindrical core to be measured as origin, it is also possible to cylinder
The central point of the bottom surface of shape core to be measured as origin, in embodiments of the present invention, with the upper bottom surface of core to be measured
Central point is used as origin.
In embodiments of the present invention, the origin of the coordinate system of foundation is represented that transverse axis is represented with X with O, the longitudinal axis is represented with Y,
Vertical pivot represents with Z, sets up the O-XYZ coordinate systems 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 transverse axis of coordinate system and
The direction cosines of the longitudinal axis, the horizontal strain point of above-mentioned anelasticity recovery strain value and core to be measured is determined according to direction cosine
Proportionate relationship between amount.
One strain gauge will be set with each monitoring position below, four monitoring positions will be set altogether, and it is adjacent each
Between monitoring position as a example by 45 ° of interval, the horizontal strain component for how determining above-mentioned strain value and core to be measured is discussed in detail
Between proportionate relationship.
As shown in figure 3, aforementioned four monitoring position is designated as the first monitoring position a respectively1, second monitoring position a2, the 3rd prison
Location puts a3With the 4th monitoring position a4, a strain gauge, the first monitoring position a are respectively provided with each monitoring position1On two
Individual strain is calculated as C1, and the strain value that strain gauge C1 measurements are obtained is defined as into 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 measurements are obtained is defined as the second monitoring position
a2Corresponding strain value, the 3rd monitoring position a3On two strains be calculated as C2, and the strain value that strain gauge C2 measurements are obtained is true
It is set to the 3rd 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 for measuring is defined as the 4th monitoring position a4Corresponding strain value.
First is calculated respectively monitors position a1, second monitoring position a2, the 3rd monitoring position a3With the 4th monitoring position a4With
Parallel direction cosines horizontally and vertically of the tangent line relative to above-mentioned coordinate system in the bottom surface of core to be measured, direction cosine refers to
Be above-mentioned tangent line respectively with the cosine value of 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, therefore, the tangent line is 0 ° with the angle of X-axis, is 90 ° with the angle of Y-axis, thus,
The first monitoring position a can be calculated11 and 0 are respectively with the direction cosines of X-axis and Y-axis respectively, the second monitoring position a2With treat
The tangent line that the bottom surface of measurement core is parallel is parallel with Y-axis, therefore, the tangent line is 90 ° with the angle of X-axis, and the angle with Y-axis is
0 °, thus, it is possible to calculate the second monitoring position is respectively 0 and 1, the 3rd monitoring position with the direction cosines of X-axis and Y-axis respectively
a3The tangent line parallel with the bottom surface of core to be measured is 45 ° with the angle of X-axis, is 135 ° with the angle of Y-axis, thus, it is possible to count
Calculate the 3rd monitoring position a30.707 and -0.707, the 4th monitoring position a are respectively with the direction cosines of X-axis and Y-axis respectively4With
The parallel tangent line in the bottom surface of core to be measured is 135 ° with the angle of X-axis, is 45 ° with the angle of Y-axis, thus, it is possible to calculate
4th monitoring position a4- 0.707 and 0.707, above-mentioned first monitoring position a are respectively with the direction cosines of X-axis and Y-axis respectively1、
Second monitoring position a2, the 3rd monitoring position a3With the 4th monitoring position a4Respectively with the direction cosines such as institute of table 1 of X-axis and Y-axis
Show.
Table 1
Wherein, in upper table, the monitoring position tangent line direction cosines relative to X-axis parallel with the bottom surface of core to be measured
Use liRepresent, monitor parallel direction cosines m of the tangent line relative to Y-axis in the bottom surface of position core to be measurediRepresent, wherein, i
=1,2,3,4.When i values are 1, l1Represent that first monitors the position tangent line parallel with the bottom surface of core to be measured relative to X
The direction cosines of axle, m1Represent first monitor the position tangent line parallel with the bottom surface of core to be measured relative to Y-axis direction more than
String.
When determining each monitoring position tangent line parallel with the bottom surface of core to be measured relative to X-axis and the direction of Y-axis
After cosine, according to equation below determine each corresponding strain value in monitoring position and core to be measured horizontal strain component it
Between proportionality coefficient,
The horizontal component of the corresponding anelasticity recovery strain value in monitoring position 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
It is to be：B=A ε ', wherein, A is the above-mentioned proportionality coefficient for calculating in the formula, and b is corresponding stagnant for each monitoring position
Elastic recovery strain value, ε ' is the horizontal strain component of core to be measured, and ε '=[ε 'x,ε′y,ε′xy]T, wherein, ε 'xIt is X-axis
The components of strain in direction, ε 'yIt is the components of strain of Y direction, ε 'xyIt is shearing strain.
S140, according to the above-mentioned anelasticity recovery strain value and aforementioned proportion relation that determine, determines the water of core to be measured
The flat components of strain.
Can determine that each monitors the corresponding anelasticity recovery strain in position by each strain gauge monitored on position
Value, if setting a strain gauge on each monitoring position, the strain value of the stabilization that the strain gauge is gathered is defined as this
The corresponding strain value in monitoring position, if being provided with one group of strain gauge on each monitoring position, by all strains in the group
The average value for counting the strain value of the stabilization of collection is defined as the corresponding strain value in monitoring position.
The corresponding strain value b in first monitoring position1Represent, the corresponding strain value b in the second monitoring position2Represent, the 3rd
The corresponding strain value b in monitoring position3Represent, the corresponding strain value b in the 4th monitoring position4Represent, then b=[b1,b2,b3,
b4]TIt is the strain value of core to be measured.
According to the formula b=A ε ' and strain value b=[b of determination core to be measured1,b2,b3,b4]T, can calculate
Horizontal strain component, specific calculating process is as follows：
The transposed matrix of the proportionality coefficient between strain value and horizontal strain component is respectively multiplied by the both sides of formula b=A ε ',
I.e.：
ATA ε '=AT b
By solving above-mentioned formula, it can be deduced that ε '=(AT A)-1AT b。
S150, according to above-mentioned horizontal strain component, determines the horizontal maximum strain value of core to be measured and horizontal minimum strain
Value.
The horizontal strain component of core to be measured is defined by the above method, how next introduce should according to level
Become component and determine the horizontal maximum strain value of core to be measured and horizontal minimum strain value.
Strain regime on horizontal plane can represent that the two principal strains are exactly maximum level by two principal strain values
Strain value and horizontal minimum strain value, can be with by horizontal strain component calculated level maximum strain value and horizontal minimum strain value
Obtained by solving equation below group：
Wherein, in above-mentioned formula, λ is unknown parameter to be solved, liIt is monitoring position and the bottom surface of core to be measured
Parallel tangent line relative to X-axis direction cosines, miIt is to monitor the position tangent line parallel with the bottom surface of core to be measured relative to Y
The direction cosines of axle, ε 'xIt is the components of strain of X-direction, ε 'yIt is the components of strain of Y direction, ε 'xyWith ε 'yxIt is to cut and answers
Become, and ε 'xyWith ε 'yxIt is equal in magnitude, but sign is different.
Above-mentioned equation group is system of homogeneous linear equations, coefficient ranks during the NSC of above-mentioned equation group untrivialo solution
Formula is zero, i.e.,：
To be quadratic equation with one unknown formula after above-mentioned determinantal expansion, i.e.,：
(ε′x-λ)(ε′y-λ)-ε′xyε′yx=0
By solving above-mentioned formula, two root λ on unknown parameter λ can be obtained1And λ2, by λ1And λ2In it is larger
One root is defined as the horizontal maximum strain value of core to be measured, by λ1And λ2In a less root be defined as core to be measured
Horizontal minimum strain value.
S160, horizontal stress coefficient of variation is determined according to horizontal maximum strain value and horizontal minimum strain value.
Wherein, it is 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, including：
The anelasticity for determining core to be measured according to horizontal maximum strain value and horizontal minimum strain value recovers mean strain
Value；According to horizontal maximum strain value peace strain value, the horizontal maximum stress value of core to be measured is determined；It is minimum according to level
Strain value peace strain value, determines the horizontal minimum stress value of core to be measured；According to horizontal maximum stress value and level most
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 is determined, according to level
Maximum strain value and horizontal minimum strain value, the mean strain value of core to be measured is determined by equation below；
εIt is flat=(εmax+εmin)/2
Wherein, in above-mentioned formula, εIt is flatFor the anelasticity of core to be measured recovers mean strain value, εmaxIt is rock to be measured
The horizontal maximum strain value of core, εminIt is the horizontal minimum strain value of core to be measured.
According to horizontal strain value peace strain value, by equation below, horizontal maximum stress value and level are determined most
Small stress value；
σi=ei/Jas+em/Jav+P0
Wherein, in above-mentioned formula, the value of i is H and h, represents horizontal maximum stress value when i is H, the table when i is h
Show horizontal minimum stress value, eiRecover the difference of strain, i.e. horizontal strain value and mean strain value partially, mean strain for anelasticity
It is the average value of horizontal maximum strain value and horizontal minimum strain value to be worth, when i is H, eiFor horizontal maximum strain value with it is average
The difference of strain value, when i is h, eiIt is horizontal maximum strain value and the difference of mean strain value, emRecovering body for anelasticity should
Become, i.e., anelasticity recovers 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, P0It is pore pressure, σiIt is horizontal stress value.
Specifically, according to horizontal maximum strain value peace strain value, by equation below, determining the water of core to be measured
Flat maximum stress value is：
σH=(εmax-εIt is flat)/Jas+εIt is flat/Jav+P0
According to horizontal minimum strain value peace strain value, by equation below, determine that the level of core to be measured is minimum
Stress value is：
σh=(εmin-εIt is flat)/Jas+εIt is flat/Jav+P0
By above-mentioned formula, it may be determined that the ratio gone out between horizontal maximum stress value and horizontal minimum stress value is：
Because the value of Jav/Jas is certain, the value of Jas/Jav is 1.56 in some specific embodiments, therefore, level
Ratio between maximum stress value and horizontal minimum stress value can be defined as：
Ratio between the horizontal maximum strain value that will be had confirmed and horizontal minimum strain value substitutes into above-mentioned formula, can
To determine the ratio between horizontal maximum stress value and horizontal minimum stress value, should by horizontal maximum strain value and level minimum
Ratio between variate is designated as εmax/εmin=K, therefore, it can draw：
In above formula, P0It is known, Jav and εminKnown to being, therefore, can determine level most by above-mentioned
Big ratio between stress value and horizontal minimum stress value.
In addition, it can also be seen that ratio between horizontal maximum stress value and horizontal minimum stress value from above-mentioned formula
Increase with the increase of K values.Compare when between the horizontal maximum stress value that core to be measured is determined and horizontal minimum stress value
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：
Difference between the ratio and 1 of calculated level maximum stress value and horizontal minimum stress value；Above-mentioned difference is determined
It is 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, σHIt is horizontal maximum stress value, σhIt is horizontal minimum stress value, above-mentioned formula is carried out into deformation can be with
Obtain,Therefore, by between the ratio and 1 between horizontal maximum stress value and horizontal minimum stress value
Difference is defined as horizontal stress coefficient of variation.
When the horizontal stress coefficient of variation of shale reservoir is measured, because mineral dehydration and pore pressure loss also can bands
To strain, in calculated level stress difference coefficient, in addition it is also necessary to consider strain value and pore pressure ease caused by mineral dehydration
Strain value caused by dissipating, method provided in an embodiment of the present invention also includes：
It is true according to strain value caused by strain value caused by horizontal maximum strain value, mineral dehydration and pore pressure loss
Determine the actual maximum strain value of level；According to strain value and pore pressure loss caused by horizontal minimum strain value, mineral dehydration
Caused strain value determines the actual minimum strain value of level；According to the actual maximum strain value of level and the actual minimum strain value of level
Determine horizontal stress coefficient of variation.
In step S150, it has been determined that go out horizontal maximum strain value for εmax, horizontal minimum strain value is εmin, will be above-mentioned
Strain value is designated as ε caused by mineral dehydrationkw, strain value caused by above-mentioned pore pressure loss is designated as εkx, therefore, it can lead to
Cross formula εmax=εASR max+εkw+εkxThe actual maximum strain value of level is calculated, wherein, ε in above-mentioned formulaASRmaxExpression level
Actual maximum strain value, i.e., the maximum strain value that core to be measured causes by Original strata stress release, according to formula εmin=εASR min
+εkw+εkxThe actual minimum strain value of level is calculated, wherein, εASR minThe actual minimum strain value of expression level, i.e., core to be measured
Because of the minimum strain value that Original strata stress release causes.
As the actual maximum strain value ε of the level for defining core to be measuredASRmaxWith the actual minimum strain value of level
εASR minAfterwards, the actual maximum strain value ε of calculated levelASRmaxMinimum strain value ε actual with levelASR minAverage value εASR puts down, its
Middle εASR puts down=(εASR max+εASR min)/2。
Next, according to formulaCalculated level
The ratio of actual maximum stress value and the actual minimum stress value of level.
In above-mentioned formula, σH realitiesIt is the actual maximum stress value of level, σH realitiesIt is the actual minimum stress value of level, P0It is hole
Pressure, Jav is the non-anelastic strain flexibility of volume.
εASRmax、εASR minAnd εASR puts down、P0Be with Jav it is known, therefore, level reality can be determined by above-mentioned formula
Ratio between border maximum stress value and the actual minimum stress value of level.
After ratio between the actual maximum stress value of the level that determines and the actual minimum stress value of level, according toDetermine the horizontal stress coefficient of variation of core to be measured.
Horizontal stress coefficient of variation measuring method provided in an embodiment of the present invention, only by measuring the stagnant bullet of cylindric core
Property recovery strain value, and according to known anelastic strain recover parameter, it is possible to determine the horizontal stress of shale reservoir
Coefficient of variation, measuring speed is fast, simple to operate, low cost.
Embodiment 2
A kind of horizontal stress coefficient of variation measurement apparatus are the embodiment of the invention provides, as shown in figure 4, the device includes the
One determining module 410, set up module 420, the second determining module 430, the 3rd determining module 440, the 4th determining module 450 and
Five determining modules 460；
Above-mentioned first determining module 410, for determining that multiple monitors the stagnant of positions on columned core curved surface to be measured
Elastic recovery strain value；
It is above-mentioned to set up module 420, for using the midpoint of the bottom surface of core to be measured as origin, to cross two of the midpoint
Vertical line, to cross the axis at the midpoint as vertical pivot, sets up coordinate system as horizontally and vertically；
Above-mentioned second determining module 430, the tangent line parallel with the bottom surface of core to be measured for determining each monitoring position
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 3rd determining module 440, for above-mentioned anelasticity recovery strain value and proportionate relationship according to determination, it is determined that
The horizontal strain component of core to be measured；
Above-mentioned 4th determining module 450, for according to above-mentioned horizontal strain component, determining that the level of core to be measured is maximum
Strain value and horizontal minimum strain value；
Above-mentioned 5th determining module 460, for determining to treat according to above-mentioned horizontal maximum strain value and horizontal minimum strain value
Measure the horizontal stress coefficient of variation of core.
Horizontal stress coefficient of variation measurement apparatus provided in an embodiment of the present invention, only by measuring the stagnant bullet of cylindric core
Property recovery strain value, and according to known anelastic strain recover parameter, it is possible to determine the horizontal stress of shale reservoir
Coefficient of variation, measuring speed is fast, simple to operate, low cost.
The horizontal stress coefficient of variation measurement apparatus that the embodiment of the present invention is provided can be equipment on specific hardware or
Person is installed on software or firmware in equipment etc..The technology of the device that the embodiment of the present invention is provided, its realization principle and generation
Effect is identical with preceding method embodiment, and to briefly describe, device embodiment part does not refer to part, refers to preceding method reality
Apply corresponding contents in example.It is apparent to those skilled in the art that, it is for convenience and simplicity of description, described above
System, device and unit specific work process, may be referred to the corresponding process in above method embodiment, herein no longer
Repeat.
In embodiment provided by the present invention, it should be understood that disclosed apparatus and method, can be by other sides
Formula is realized.Device embodiment described above is only schematical, for example, the division of the unit, only one kind are patrolled
Collect function to divide, there can be other dividing mode when actually realizing, but for example, multiple units or component can combine or can
To be integrated into another system, or some features can be ignored, or not perform.It is another, it is shown or discussed each other
Coupling or direct-coupling or communication connection can be the INDIRECT COUPLING or communication link of device or unit by some communication interfaces
Connect, can be electrical, mechanical or other forms.
The unit that is illustrated as separating component can be or may not be it is physically separate, it is aobvious as unit
The part for showing can be or may not be physical location, you can with positioned at a place, or can also be distributed to multiple
On NE.Some or all of unit therein can be according to the actual needs selected to realize the mesh of this embodiment scheme
's.
In addition, during each functional unit in the embodiment that the present invention is provided can be integrated in a processing unit, also may be used
Being that unit is individually physically present, it is also possible to which two or more units are integrated in a unit.
If the function is to realize in the form of SFU software functional unit and as independent production marketing or when using, can be with
Storage is in a computer read/write memory medium.Based on such understanding, technical scheme is substantially in other words
The part contributed to prior art or the part of the technical scheme can be embodied in the form of software product, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used to so that a computer equipment (can be individual
People's computer, server, or network equipment etc.) perform all or part of step of each embodiment methods described of the invention.
And foregoing storage medium includes：USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only Memory), arbitrary access are deposited
Reservoir (RAM, Random Access Memory), magnetic disc or CD etc. are various can be with the medium of store program codes.
It should be noted that：Similar label and letter represents similar terms in following accompanying drawing, therefore, once a certain Xiang Yi
It is defined in individual accompanying drawing, then it need not be further defined and explained in subsequent accompanying drawing, additionally, term " the
One ", " second ", " the 3rd " etc. are only used for distinguishing description, and it is not intended that indicating or implying relative importance.
Finally it should be noted that：Embodiment described above, specific embodiment only of the invention, is used to illustrate the present invention
Technical scheme, rather than its limitations, protection scope 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, it will be understood by those within the art that：Any one skilled in the art
The invention discloses technical scope in, it can still modify to the technical scheme described in previous embodiment or can be light
Change is readily conceivable that, or equivalent is carried out to which part technical characteristic；And these modifications, change or replacement, do not make
The essence of appropriate technical solution departs from the spirit and scope of embodiment of the present invention technical scheme.Should all cover in protection of the invention
Within the scope of.Therefore, protection scope of the present invention described should be defined by scope of the claims.
Claims (10)
1. a kind of horizontal stress coefficient of variation measuring method, it is characterised in that methods described includes：
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, so that two vertical lines at excessively described midpoint are as transverse axis and indulge
Axle, the axis using excessively described midpoint sets up coordinate system as vertical pivot；
Determine each monitoring position tangent line horizontal stroke relative to the coordinate system parallel with the bottom surface of the core to be measured
The direction cosines of axle and the longitudinal axis, determine the anelasticity recovery strain value with the core to be measured according to the direction cosines
Proportionate relationship between horizontal strain component；
According to the anelasticity recovery strain value and the proportionate relationship that determine, the horizontal strain of the core to be measured is determined
Component；
According to the horizontal strain component, the horizontal maximum strain value and horizontal minimum strain value of the core to be measured are determined；
Determine that the horizontal stress of the core to be measured is poor according to the horizontal maximum strain value and the horizontal minimum strain value
Different coefficient.
2. method according to claim 1, it is characterised in that described according to the horizontal maximum strain value and the level
Minimum strain value determines the horizontal stress coefficient of variation of the core to be measured, including：
Determine that the anelasticity of core to be measured recovers 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 horizontal maximum stress value and the ratio of the horizontal minimum stress value, the horizontal stress difference system is determined
Number.
3. method according to claim 2, it is characterised in that described according to the horizontal maximum stress value and the level
The ratio of minimum stress value, determines the horizontal stress coefficient of variation, including：
Calculate the difference between the ratio and 1 of the horizontal maximum stress value and the horizontal minimum stress value；
The difference is defined as the horizontal stress coefficient of variation.
4. method according to claim 1, it is characterised in that methods described 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
Determine the actual maximum strain value of level；
It is true according to strain value caused by strain value caused by the horizontal minimum strain value, mineral dehydration and pore pressure loss
Determine the actual minimum strain value of level；
The horizontal stress difference system is determined according to the actual maximum strain value of the level and the actual minimum strain value of the level
Number.
5. method according to claim 1, it is characterised in that each described monitoring position is columned described to be measured
It is uniformly distributed on core curved surface.
6. method according to claim 1, it is characterised in that 45 ° are spaced between the adjacent monitoring position.
7. method according to claim 1, it is characterised in that the number of the monitoring position is more than or equal to 3.
8. method according to claim 1, it is characterised in that multiple on the columned core curved surface to be measured of determination
The anelasticity recovery strain value of position is monitored, including：
Each described monitoring position is gathered respectively in the corresponding strain value of multiple predetermined times；
Count the number of times that each corresponding strain value in monitoring position occurs in multiple predetermined times；
Will appear from the most strain value of number of times and be defined as the corresponding anelasticity recovery strain value in the monitoring position.
9. method according to claim 8, it is characterised in that be provided with strain gauge at each described monitoring position, pass through
The strain gauge gathers each described monitoring position in the corresponding strain value of multiple predetermined times.
10. a kind of horizontal stress coefficient of variation measurement apparatus, it is characterised in that described device includes：
First determining module, the anelasticity recovery strain for determining multiple monitoring positions on columned core curved surface to be measured
Value；
Module is set up, it is vertical with two of excessively described midpoint for using the midpoint of the bottom surface of the core to be measured as origin
Used as horizontally and vertically, the axis using excessively described midpoint sets up coordinate system to line as vertical pivot；
Second determining module, it is relative for determining each described monitoring position tangent line parallel with the bottom surface of the core to be measured
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；
3rd determining module, for according to the anelasticity recovery strain value and the proportionate relationship for determining, it is determined that described treat
Measure the horizontal strain component of core；
4th determining module, for according to the horizontal strain component, determining the horizontal maximum strain value of the core to be measured
With horizontal minimum strain value；
5th determining module, it is described to be measured for being determined according to the horizontal maximum strain value and the horizontal minimum strain value
The horizontal stress coefficient of variation of core.
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 true CN106768581A (en) | 2017-05-31 |
| CN106768581B 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) |
Cited By (6)
| 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 |
| CN110243526A (en) * | 2019-07-03 | 2019-09-17 | 东北大学 | A method of it is measured based on flat jack and determines Three-dimensional Rock initial field stress |
| CN111982377A (en) * | 2020-09-01 | 2020-11-24 | 中国地质科学院地质力学研究所 | Deep hole differential stress evaluation method based on secondary sleeve core |
| CN111999165A (en) * | 2020-08-31 | 2020-11-27 | 北京科技大学 | Deep high-stress rock elastic strain recovery monitoring device and method |
| CN112129448A (en) * | 2020-06-07 | 2020-12-25 | 中铁二院工程集团有限责任公司 | Method for measuring ground stress by using elastic recovery deformation of core |
| CN114969906A (en) * | 2022-05-16 | 2022-08-30 | 广西交科集团有限公司 | Method for testing prestress application efficiency of corrugated steel web bridge |
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 | 中国石油天然气股份有限公司 | Oil-gas well sand prevention method and device based on geomechanical model |
-
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 | 中国石油天然气股份有限公司 | Oil-gas well sand prevention method and device based on geomechanical model |
Cited By (9)
| 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 |
| CN110243526A (en) * | 2019-07-03 | 2019-09-17 | 东北大学 | A method of it is measured based on flat jack and determines Three-dimensional Rock initial field stress |
| 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 |
| CN111999165A (en) * | 2020-08-31 | 2020-11-27 | 北京科技大学 | Deep high-stress rock elastic strain recovery monitoring device and method |
| CN111999165B (en) * | 2020-08-31 | 2021-10-22 | 北京科技大学 | Deep high-stress rock elastic strain recovery monitoring device and method |
| CN111982377A (en) * | 2020-09-01 | 2020-11-24 | 中国地质科学院地质力学研究所 | Deep hole differential stress evaluation method based on secondary sleeve core |
| CN111982377B (en) * | 2020-09-01 | 2021-07-23 | 中国地质科学院地质力学研究所 | Deep hole differential stress evaluation method based on secondary sleeve core |
| CN114969906A (en) * | 2022-05-16 | 2022-08-30 | 广西交科集团有限公司 | Method for testing prestress application efficiency of corrugated steel web bridge |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106768581B (en) | 2019-07-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106768581A (en) | A kind of horizontal stress coefficient of variation measuring method and device | |
| US5159833A (en) | Analyzing a hydrocarbon reservoir by determining the response of that reservoir to tidal forces | |
| CN106768580B (en) | A kind of horizontal stress measurement method and device | |
| CN109033504B (en) | Oil-water well casing damage prediction method | |
| CN105160414B (en) | Method and device for predicting type of full oil reservoir production area | |
| US10641923B2 (en) | Method for characterizing and exploiting a subterranean formation comprising a network of fractures | |
| CN105221140A (en) | A kind ofly determine that shale formation can the method for pressure break sex index | |
| CN109241588A (en) | A kind of analogy method of the monolete extension based on quasi-continuous geomechanics model | |
| US20130250726A1 (en) | Quantitative analysis of time-lapse seismic data | |
| CN108073786A (en) | Shale compressibility evaluation method and device | |
| CN106407171B (en) | Tree structure-based security monitoring custom report generation system and method | |
| Neuhaus et al. | Integrated microseismic monitoring for field optimization in the Marcellus Shale-A case study | |
| CN114510854A (en) | Method for evaluating accuracy of acid fracturing numerical simulation result of hole finding in seam following manner | |
| CN112946783A (en) | Method, device and equipment for determining hydrate saturation | |
| CN117852416A (en) | Multimode grouting precontrolled analysis method and system based on digital geological model | |
| CN110863814B (en) | Method and device for determining single-well subsection specific fluid production index of ultra-thick type cuttings limestone oil reservoir | |
| Jianguo Du et al. | Geochemical anomalies connected with great earthquakes in China | |
| CN115270062B (en) | Stress relief method for calculating crustal stress by taking irregular drilling hole shape into consideration | |
| CN107219560B (en) | Mine worked-out section deflection and stability assessment method based on gravity anomaly | |
| CN112241835B (en) | Multi-source information evaluation method for flood disaster in deep shaft engineering | |
| CN107346035A (en) | Method and device for identifying fracture | |
| CN107036575B (en) | Mine worked-out section deformation and Detection of Stability method based on bouguer gravity anomaly | |
| CN110109176B (en) | Rock burst risk assessment method under pathological condition of tunnel microseismic sensor monitoring station network | |
| US12091972B1 (en) | Method and device for determining development status of separation layer of coal mine | |
| US12056422B1 (en) | System for designing perforating gun capable of reducing energy consumption |
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 |