CN107829727A - A kind of method and device for determining crustal stress - Google Patents
A kind of method and device for determining crustal stress Download PDFInfo
- Publication number
- CN107829727A CN107829727A CN201710817354.9A CN201710817354A CN107829727A CN 107829727 A CN107829727 A CN 107829727A CN 201710817354 A CN201710817354 A CN 201710817354A CN 107829727 A CN107829727 A CN 107829727A
- Authority
- CN
- China
- Prior art keywords
- interval
- interest
- opening position
- svel
- drilling well
- 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.)
- Pending
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/14—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
Abstract
The embodiment of the present application discloses a kind of method and device for determining crustal stress.Methods described is provided with the SVEL of core sample in interval of interest, and methods described includes:SVEL based on core sample, establish incidence relation of the interval of interest in the SVEL and the SVEL of non-perpendicular bed plane in the vertical stratification face of drilling well opening position;Based on incidence relation and interval of interest in the SVEL in the vertical stratification face of drilling well opening position, SVEL of the interval of interest in the non-perpendicular bed plane of drilling well opening position is determined;According to the density log curve of interval of interest, and interval of interest determines minimum horizontal principal stress and maximum horizontal principal stress of the interval of interest in drilling well opening position in the SVEL in the vertical stratification face of drilling well opening position and the SVEL of non-perpendicular bed plane.The technical scheme that the embodiment of the present application provides, the degree of accuracy of identified crustal stress can be improved.
Description
Technical field
The application is related to oil gas development technology field, more particularly to a kind of method and device for determining crustal stress.
Background technology
Crustal stress is that earth solid dielectric is added by gravity, a variety of earth structure power and astronomical factor and spy pick engineering
The effect of power, cause the mechanics parameter of response deformation in media interior unit, there is important influence to oil-gas exploration and development
And meaning.
It is commonly used in the prior art determine crustal stress the main process of method be:Based on transverse isotropy model, adopt
The shear wave velocity V vertical with bed plane of stratum opening position is determined with the method for acoustic loggingSHIt is (0 °), vertical with bed plane vertical
Wave velocity VP(0 °) and Stoneley wave speed, and according to shear wave velocity VSH(0 °), velocity of longitudinal wave VP(0 °) and Stoneley wave
Speed calculates the stiffness coefficient C at the position of stratum11And C13, then calculated under certain approximate condition at the position of stratum
Stiffness coefficient C33、C44And C66, it is last according to the stiffness coefficient C being calculated11、C13、C33、C44And C66, calculate the stratum position
Put the horizontal minimum principal stress at place and horizontal maximum principal stress.However, be difficult often to obtain Stoneley wave by conventional logging, and
And acquired Stoneley wave is easily influenceed by factors such as in-place permeability, crack and pit shaft degree of irregularities, may cause really
The degree of accuracy of fixed Stoneley wave speed is relatively low, so as to influence the result of calculation of final crustal stress.Moreover, prior art is worked as
Employed in approximate condition it is inapplicable when, may cause to calculate the stiffness coefficient C at the position of stratum33、C44And C66It is accurate
Spend it is relatively low, so as to cause the degree of accuracy of horizontal minimum principal stress and horizontal maximum principal stress that is calculated relatively low.
The content of the invention
The purpose of the embodiment of the present application is to provide a kind of method and device for determining crustal stress, with determined by raising should
The degree of accuracy of power.
In order to solve the above technical problems, it is such reality that the embodiment of the present application, which provides a kind of method and device for determining crustal stress,
Existing:
A kind of method for determining crustal stress, there is provided the SVEL of core sample, methods described include in purposeful interval:
Based on the SVEL of the core sample, the interval of interest is established in the vertical stratification face of drilling well opening position
The incidence relation of SVEL and the SVEL of non-perpendicular bed plane;
SVEL based on the incidence relation and the interval of interest in the vertical stratification face of drilling well opening position, it is determined that
SVEL of the interval of interest in the non-perpendicular bed plane of drilling well opening position;
According to the density log curve of the interval of interest, and the interval of interest is in the vertical stratification of drilling well opening position
The SVEL of the SVEL in face and non-perpendicular bed plane, determine minimum level master of the interval of interest in drilling well opening position
Stress and maximum horizontal principal stress.
In preferred scheme, the SVEL of the core sample includes the SVEL in vertical stratification face and non-perpendicular stratification
The SVEL in face;
Correspondingly, the SVEL based on the core sample, the interval of interest is established in drilling well opening position
The incidence relation of the SVEL in vertical stratification face and the SVEL of non-perpendicular bed plane, including:
The SVEL in the vertical stratification face based on the core sample and the SVEL of non-perpendicular bed plane, using most
The method that a young waiter in a wineshop or an inn multiplies fitting establishes SVEL and non-perpendicular layer of the interval of interest in the vertical stratification face of drilling well opening position
The incidence relation of the SVEL in reason face.
In preferred scheme, the SVEL in the vertical stratification face includes:First shear wave velocity and the first velocity of longitudinal wave;Its
In, first shear wave velocity represents the shear wave velocity in vertical stratification face, and first velocity of longitudinal wave represents the vertical stratification
The velocity of longitudinal wave in face;The SVEL of the non-perpendicular bed plane includes:Second shear wave velocity, the second velocity of longitudinal wave and the 3rd are vertical
Wave velocity;Wherein, the velocity of longitudinal wave that it is 45 degree with the angle of bed plane that second velocity of longitudinal wave, which represents, the 3rd compressional wave speed
The velocity of longitudinal wave that it is 0 degree with the angle of bed plane that degree, which represents,;Second shear wave velocity represents that with the angle of bed plane be 0 degree
Shear wave velocity;
The first shear wave velocity and the second shear wave velocity based on the core sample, are built using the method for least square fitting
The interval of interest is found in the second shear wave velocity of drilling well opening position and the incidence relation of the first shear wave velocity;
The first velocity of longitudinal wave and the second velocity of longitudinal wave based on the core sample, are built using the method for least square fitting
The interval of interest is found in the second velocity of longitudinal wave of drilling well opening position and the incidence relation of the first velocity of longitudinal wave;
The first velocity of longitudinal wave and the 3rd velocity of longitudinal wave based on the core sample, are built using the method for least square fitting
The interval of interest is found in the 3rd velocity of longitudinal wave of drilling well opening position and the incidence relation of the first velocity of longitudinal wave.
In preferred scheme, the density log curve according to the interval of interest, and the interval of interest is in drilling well
The SVEL in the vertical stratification face of opening position and the SVEL of non-perpendicular bed plane, determine the interval of interest in well drilling positions
The minimum horizontal principal stress and maximum horizontal principal stress at place are put, including:
According to the density log curve, density of the interval of interest in drilling well opening position is determined;
According to the interval of interest in the SVEL in the vertical stratification face of drilling well opening position, the sound wave of non-perpendicular bed plane
Speed, determine at least one stiffness coefficient of the interval of interest in drilling well opening position;Wherein, described in the stiffness coefficient represents
Stiffness coefficient in the stiffness matrix of transverse isotropy petrophysical model corresponding to interval of interest;
According to the interval of interest in the density and stiffness coefficient of drilling well opening position, determine the interval of interest in well drilling positions
Put the minimum horizontal principal stress and maximum horizontal principal stress at place.
In preferred scheme, it is described according to the interval of interest drilling well opening position density and stiffness coefficient, it is determined that described
Interval of interest drilling well opening position minimum horizontal principal stress and maximum horizontal principal stress, including:
According to the interval of interest in the stiffness coefficient of drilling well opening position, determine the interval of interest in drilling well opening position
Kinetic Youngs modulus and dynamic Poisson's ratio;Wherein, the kinetic Youngs modulus includes dynamic level Young's modulus and dynamic is vertical
Young's modulus, the dynamic Poisson's ratio include dynamic level Poisson's ratio and dynamic vertical Poisson's ratio;
Determine static corresponding to kinetic Youngs modulus and dynamic Poisson's ratio difference of the interval of interest in drilling well opening position
Young's modulus and static Poisson's ratio;Wherein, the static Young's modulus includes quiescent levels Young's modulus and static vertical Young
Modulus, the static Poisson's ratio include quiescent levels Poisson's ratio and static vertical Poisson's ratio;
Determine Biot spy coefficient of the interval of interest in drilling well opening position;
According to the interval of interest in the special system of the density of drilling well opening position, static Young's modulus, static Poisson's ratio and Biot
Number, determines minimum horizontal principal stress and maximum horizontal principal stress of the interval of interest in drilling well opening position.
In preferred scheme, the kinetic Youngs modulus and dynamic Poisson's ratio for determining the interval of interest in drilling well opening position
Static Young's modulus and static Poisson's ratio corresponding to respectively, including:
According to the SVEL in the vertical stratification face of the core sample and the SVEL of non-perpendicular bed plane, institute is established
Interval of interest is stated in the kinetic Youngs modulus of drilling well opening position and the incidence relation of static Young's modulus, and establishes the purpose
Interval is in the dynamic Poisson's ratio of drilling well opening position and the incidence relation of static Poisson's ratio;
According to the interval of interest in the kinetic Youngs modulus of drilling well opening position and the incidence relation of static Young's modulus, with
And dynamic Poisson's ratio and the incidence relation of static Poisson's ratio, determine kinetic Youngs modulus of the interval of interest in drilling well opening position
With dynamic Poisson's ratio respectively corresponding to static Young's modulus and static Poisson's ratio.
In preferred scheme, the SVEL in the vertical stratification face according to the core sample and non-perpendicular bed plane
SVEL, the interval of interest is established in the kinetic Youngs modulus of drilling well opening position and the incidence relation of static Young's modulus,
And the interval of interest is established in the dynamic Poisson's ratio of drilling well opening position and the incidence relation of static Poisson's ratio, including:
According to the SVEL in the vertical stratification face of the core sample and the SVEL of non-perpendicular bed plane, institute is determined
State the stiffness coefficient of core sample;
According to the stiffness coefficient of the core sample, the kinetic Youngs modulus and dynamic Poisson of the core sample are determined
Than;
Determine the static Young's modulus of the core sample and static Poisson's ratio;
According to the static Young's modulus and kinetic Youngs modulus of the core sample, the interval of interest is established in well drilling positions
Put the kinetic Youngs modulus at place and the incidence relation of static Young's modulus, and according to the static Poisson's ratio of the core sample and
Dynamic Poisson's ratio, the interval of interest is established in the dynamic Poisson's ratio of drilling well opening position and the incidence relation of static Poisson's ratio.
In preferred scheme, the Biot spy's coefficient for determining the interval of interest in drilling well opening position, including:
Determine Biot spy's coefficient of the core sample;
Using the average value of Biot spy's coefficient of multiple core samples as the interval of interest in drilling well opening position
Biot spy's coefficient.
In preferred scheme, the core sample includes the core sample that geometry is cylinder;Wherein, the axle of the cylinder
Line is parallel or vertical with bed plane.
A kind of device for determining crustal stress, described device provide the SVEL of core sample in interval of interest;The dress
Put including:Incidence relation determining module, interval SVEL determining module and crustal stress determining module;Wherein,
The incidence relation determining module, for the SVEL based on the core sample, establish the interval of interest
In the incidence relation of the SVEL and the SVEL of non-perpendicular bed plane in the vertical stratification face of drilling well opening position;
The interval SVEL determining module, for based on the incidence relation and the interval of interest in drilling well position
The SVEL in the vertical stratification face at place, determine sound wave speed of the interval of interest in the non-perpendicular bed plane of drilling well opening position
Degree;
The crustal stress determining module, for the density log curve according to the interval of interest, and the target zone
Section determines the interval of interest in the SVEL in the vertical stratification face of drilling well opening position and the SVEL of non-perpendicular bed plane
In the minimum horizontal principal stress and maximum horizontal principal stress of drilling well opening position.
A kind of method and device for determining crustal stress of the embodiment of the present application, can based on the SVEL of the core sample
To establish the interval of interest in the SVEL in the vertical stratification face of drilling well opening position and the SVEL of non-perpendicular bed plane
Incidence relation;SVEL based on the incidence relation and the interval of interest in the vertical stratification face of drilling well opening position,
SVEL of the interval of interest in the non-perpendicular bed plane of drilling well opening position can be determined;Can be according to the interval of interest
Density log curve, and SVEL and non-perpendicular stratification of the interval of interest in the vertical stratification face of drilling well opening position
The SVEL in face, determine minimum horizontal principal stress and maximum horizontal principal stress of the interval of interest in drilling well opening position.By
It is near without introducing without obtaining Stoneley wave in it is determined that during minimum horizontal principal stress and maximum horizontal principal stress
Like condition, so higher using the degree of accuracy of crustal stress determined by the present processes.
Brief description of the drawings
, below will be to embodiment or existing in order to illustrate more clearly of the embodiment of the present application or technical scheme of the prior art
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments described in application, for those of ordinary skill in the art, do not paying the premise of creative labor
Under, other accompanying drawings can also be obtained according to these accompanying drawings.
Fig. 1 is the schematic diagram of core sample in the embodiment of the present application;
Fig. 2 is a kind of flow chart for the embodiment of the method for determining crustal stress of the application;
Fig. 3 is the composition structure chart for the device embodiment that the application determines crustal stress.
Embodiment
The embodiment of the present application provides a kind of method and device for determining crustal stress.
In order that those skilled in the art more fully understand the technical scheme in the application, it is real below in conjunction with the application
The accompanying drawing in example is applied, the technical scheme in the embodiment of the present application is clearly and completely described, it is clear that described implementation
Example only some embodiments of the present application, rather than whole embodiments.It is common based on the embodiment in the application, this area
The every other embodiment that technical staff is obtained under the premise of creative work is not made, it should all belong to the application protection
Scope.
The embodiment of the present application provides a kind of method for determining crustal stress.The method for determining crustal stress provides purposeful
The SVEL of core sample in interval.
In the present embodiment, the interval of interest can be the still undetermined interval of crustal stress.
In the present embodiment, multiple core samples can be obtained in the interval of interest.Acoustic wave of rock can be passed through
Measuring instrument measures the attribute information of each core sample, to provide data basis for subsequent step.
In the present embodiment, the SVEL of the core sample can include:The SVEL in vertical stratification face and
The SVEL of non-perpendicular bed plane.Wherein, the SVEL in the vertical stratification face can include:First shear wave velocity and
One velocity of longitudinal wave;First shear wave velocity can represent the shear wave velocity V in vertical stratification faceSH(0 °), the first compressional wave speed
Degree can represent the velocity of longitudinal wave V in the vertical stratification faceP(0°).The SVEL of the non-perpendicular bed plane can include:
Second shear wave velocity, the second velocity of longitudinal wave and the 3rd velocity of longitudinal wave;Second velocity of longitudinal wave can represent the folder with bed plane
Angle is 45 degree of velocity of longitudinal wave VP(45 °), the compressional wave speed that the 3rd velocity of longitudinal wave can represent with the angle of bed plane is 0 degree
Spend VP(90°);The shear wave velocity V that second shear wave velocity can represent with the angle of bed plane is 0 degreeSH(90°).The layer
Reason face is the divisional plane between each rock stratum layer position in the interval of interest.
In the present embodiment, the core sample can include the core sample that geometry is cylinder.Wherein, it is described
The axis of cylinder can be parallel or vertical with bed plane.The diameter of the cylinder can be 25 millimeters, and height can be 50 millimeters.Example
Such as, Fig. 1 is the schematic diagram of core sample in the embodiment of the present application.(a) and (b) is respectively that the axis hangs down with bed plane in Fig. 1
Perpendicular and parallel rock core cylindrical sample.Angle in Fig. 1 is the angle of cylinder axis and bed plane.As shown in figure 1, it is directed to axle
The line rock core cylindrical sample vertical with bed plane, can be measured by acoustic wave of rock measuring instrument the rock core cylindrical sample it is each to
Different in nature acoustic velocity, respectively shear wave velocity VSH(0°)、VSV(0°)、VSV(90 °) and VSH(90 °), and velocity of longitudinal wave VP
(0°)、VP(45 °) and VP(90°).Under normal circumstances, can be by shear wave velocity due to the anisotropic character of rock core cylindrical sample
VSH(0°)、VSV(0 °) and VSVThe average value of (90 °), the final shear wave velocity V as the rock core cylindrical sampleSH(0°).It is similar
Ground, for the axis rock core cylindrical sample parallel with bed plane, rock core circle can also be measured by acoustic wave of rock measuring instrument
The anisotropy acoustic velocity of post sample, respectively shear wave velocity VSH(90°)、VSV(90°)、VSV(0 °) and VSH(0 °), and it is vertical
Wave velocity VP(90°)、VP(45 °) and VP(0°).Can also be by shear wave velocity VSH(90°)、VSV(0 °) and VSV(90 °) are averaged
Value, the final shear wave velocity V as the rock core cylindrical sampleSH(90°)., can be with the subsequent step of the embodiment of the present application
Illustrated exemplified by the rock core cylindrical sample vertical with bed plane by axis.
Fig. 2 is a kind of flow chart for the embodiment of the method for determining crustal stress of the application.As shown in Fig. 2 described definitely should
The method of power, comprises the following steps.
Step S101:Based on the SVEL of the core sample, the interval of interest hanging down in drilling well opening position is established
The incidence relation of the SVEL and the SVEL of non-perpendicular bed plane of straight bed plane.
In the present embodiment, the drilling well position in the interval of interest can refer to by the interval of interest
The drilling well opened up, it is able to detect that the position of acoustic velocity.
In the present embodiment, the SVEL based on the core sample, the interval of interest is established in drilling well position
The incidence relation of the SVEL in the vertical stratification face at place and the SVEL of non-perpendicular bed plane, can be included, based on described
The SVEL in the vertical stratification face of core sample and the SVEL of non-perpendicular bed plane, least square fitting can be used
Method establishes the interval of interest in the SVEL in the vertical stratification face of drilling well opening position and the sound wave speed of non-perpendicular bed plane
The incidence relation of degree.Specifically, the first shear wave velocity and the second shear wave velocity based on the core sample, minimum can be used
Two methods for multiplying fitting establish the interval of interest in the second shear wave velocity of drilling well opening position and associating for the first shear wave velocity
Relation.The first velocity of longitudinal wave and the second velocity of longitudinal wave based on the core sample, the method that least square fitting can be used
The interval of interest is established in the second velocity of longitudinal wave of drilling well opening position and the incidence relation of the first velocity of longitudinal wave.Based on the rock
The first velocity of longitudinal wave and the 3rd velocity of longitudinal wave of heart sample, the method for least square fitting can be used to establish the interval of interest
In the 3rd velocity of longitudinal wave of drilling well opening position and the incidence relation of the first velocity of longitudinal wave.
Step S102:Sound wave based on the incidence relation and the interval of interest in the vertical stratification face of drilling well opening position
Speed, determine SVEL of the interval of interest in the non-perpendicular bed plane of drilling well opening position.
In the present embodiment, the method that can use acoustic logging, obtains the interval of interest in drilling well opening position
The SVEL in vertical stratification face, i.e., first shear wave velocity V of the described interval of interest in drilling well opening positionSH(0 °) and the first compressional wave
Speed VP(0°).In this way, based on the interval of interest in the second shear wave velocity of drilling well opening position and the pass of the first shear wave velocity
Connection relation, it may be determined that second shear wave velocity V of the interval of interest in drilling well opening positionSH(90°).Based on the interval of interest
In the second velocity of longitudinal wave of drilling well opening position and the incidence relation of the first velocity of longitudinal wave, it may be determined that the interval of interest is in drilling well
Second velocity of longitudinal wave V of opening positionP(45°).Based on the interval of interest in the 3rd velocity of longitudinal wave of drilling well opening position and first
The incidence relation of velocity of longitudinal wave, it may be determined that threeth velocity of longitudinal wave V of the interval of interest in drilling well opening positionP(90°)。
Step S103:According to the density log curve of the interval of interest, and the interval of interest is in drilling well opening position
The SVEL in vertical stratification face and the SVEL of non-perpendicular bed plane, determine the interval of interest in drilling well opening position
Minimum horizontal principal stress and maximum horizontal principal stress.
In the present embodiment, according to the density log curve of the interval of interest, and the interval of interest is in drilling well
The SVEL in the vertical stratification face of opening position and the SVEL of non-perpendicular bed plane, determine the interval of interest in well drilling positions
The minimum horizontal principal stress and maximum horizontal principal stress at place are put, can specifically be included, can be according to the density log curve, really
Density of the fixed interval of interest in drilling well opening position.Can be according to the interval of interest in the vertical stratification face of drilling well opening position
SVEL, the SVEL of non-perpendicular bed plane, determine at least one rigidity of the interval of interest in drilling well opening position
Coefficient.Wherein, the stiffness coefficient can represent the firm of transverse isotropy petrophysical model corresponding to the interval of interest
The stiffness coefficient spent in matrix.Can according to the interval of interest drilling well opening position density and stiffness coefficient, it is determined that described
Minimum horizontal principal stress and maximum horizontal principal stress of the interval of interest in drilling well opening position.
Specifically, expressions below can be used to characterize the stiffness matrix of the transverse isotropy petrophysical model:
Wherein, stiffness coefficient C at least one stiffness coefficient11、C33、C44、C66And C13Between be separate;
Stiffness coefficient C11Represent element, stiffness coefficient C corresponding to the first row and first row in the stiffness matrix12Represent the rigidity
Element, stiffness coefficient C corresponding to the first row and secondary series in matrix13Represent the first row and tertial member in the stiffness matrix
Element, stiffness coefficient C33Represent in the stiffness matrix element corresponding to the third line and the 3rd row, stiffness coefficient C44Represent described firm
Spend matrix in fourth line and the 4th row corresponding to element, stiffness coefficient C66Represent the 6th row and the 6th row pair in the stiffness matrix
The element answered.Stiffness coefficient C12It is the coefficient of a dependent, can be determined by following formula:C12=C11-2C66。
In the present embodiment, stiffness coefficient of the interval of interest in drilling well opening position can be determined using following formula
C11、C33、C44、C66And C13:
C11=VP(90°)2ρ
C33=VP(0°)2ρ
C44=VSH(0°)2ρ
C66=VSH(90°)2ρ
Wherein, ρ represents density of the interval of interest in drilling well opening position.
In the present embodiment, according to the interval of interest drilling well opening position density and stiffness coefficient, it is determined that described
Interval of interest can specifically include in the minimum horizontal principal stress and maximum horizontal principal stress of drilling well opening position, can be according to institute
State stiffness coefficient of the interval of interest in drilling well opening position, determine the interval of interest drilling well opening position kinetic Youngs modulus and
Dynamic Poisson's ratio;Wherein, the kinetic Youngs modulus includes dynamic level Young's modulus and the vertical Young's modulus of dynamic, described dynamic
State Poisson's ratio includes dynamic level Poisson's ratio and dynamic vertical Poisson's ratio.The interval of interest can be determined in drilling well opening position
Static Young's modulus and static Poisson's ratio corresponding to kinetic Youngs modulus and dynamic Poisson's ratio difference.Wherein, the static Young
Modulus includes quiescent levels Young's modulus and static vertical Young's modulus, the static Poisson's ratio include quiescent levels Poisson's ratio and
Static vertical Poisson's ratio.Biot spy coefficient of the interval of interest in drilling well opening position can be determined.Can be according to the purpose
Interval determines that the interval of interest exists in the density of drilling well opening position, static Young's modulus, static Poisson's ratio and Biot spy's coefficient
The minimum horizontal principal stress and maximum horizontal principal stress of drilling well opening position.
In the present embodiment, dynamic Young of the interval of interest in drilling well opening position can be determined using following formula
Modulus and dynamic Poisson's ratio:
Wherein, EhAnd EvRepresent that the interval of interest hangs down in the dynamic level Young's modulus of drilling well opening position and dynamic respectively
To Young's modulus;νhAnd νvRepresent that the interval of interest is moored in the dynamic level Poisson's ratio and dynamic vertical of drilling well opening position respectively
Loose ratio.
In the present embodiment, kinetic Youngs modulus and dynamic Poisson's ratio of the interval of interest in drilling well opening position are determined
Static Young's modulus and static Poisson's ratio corresponding to respectively, can specifically include, can be according to the perpendicular layers of the core sample
The SVEL of the SVEL in reason face and non-perpendicular bed plane, establish dynamic Young of the interval of interest in drilling well opening position
The incidence relation of modulus and static Young's modulus, and the interval of interest is established in the dynamic Poisson's ratio of drilling well opening position and quiet
The incidence relation of state Poisson's ratio.Can be according to kinetic Youngs modulus of the interval of interest in drilling well opening position and static Young mould
The incidence relation of amount, and dynamic Poisson's ratio and the incidence relation of static Poisson's ratio, determine the interval of interest in drilling well position
Static Young's modulus and static Poisson's ratio corresponding to kinetic Youngs modulus and the dynamic Poisson's ratio difference at place.
In the present embodiment, according to the SVEL in the vertical stratification face of the core sample and non-perpendicular bed plane
SVEL, the interval of interest is established in the kinetic Youngs modulus of drilling well opening position and the incidence relation of static Young's modulus,
And the interval of interest is established in the dynamic Poisson's ratio of drilling well opening position and the incidence relation of static Poisson's ratio, it can specifically wrap
Include, institute can be determined according to the SVEL in the vertical stratification face of the core sample and the SVEL of non-perpendicular bed plane
State the stiffness coefficient of core sample.According to the stiffness coefficient of the core sample, the dynamic poplar of the core sample can be determined
Family name's modulus and dynamic Poisson's ratio.The static Young's modulus of the core sample and static Poisson's ratio can be determined.Can according to institute
The static Young's modulus and kinetic Youngs modulus of core sample are stated, establishes dynamic Young of the interval of interest in drilling well opening position
The incidence relation of modulus and static Young's modulus, and static Poisson's ratio and dynamic Poisson's ratio according to the core sample, build
The interval of interest is found in the dynamic Poisson's ratio of drilling well opening position and the incidence relation of static Poisson's ratio.
In the present embodiment, triaxial compression test method can be used to measure the static Young's modulus of the core sample
With static Poisson's ratio.
In the present embodiment, Biot spy coefficient of the interval of interest in drilling well opening position is determined, can specifically be included,
Biot spy's coefficient of the core sample can be determined.The average value of Biot spy's coefficient of multiple core samples can be made
For the interval of interest drilling well opening position Biot spy's coefficient.
In the present embodiment, triaxial compression test method can also be used to measure the special system of Biot of the core sample
Number.
In the present embodiment, minimum level of the interval of interest in drilling well opening position can be determined using following formula
Principal stress and maximum horizontal principal stress:
Wherein, σhAnd σHMinimum horizontal principal stress and maximum horizontal of the interval of interest in drilling well opening position are represented respectively
Principal stress, unit are megapascal (MPa) (MPa);σvThe pressure of overlying rock is represented, density log curve and vertical depth meter can be passed through
Obtain, unit MPa;α represents the interval of interest in Biot spy's coefficient of drilling well opening position, dimensionless;σpRepresent hole
Fluid pressure, unit MPa;εhAnd εHThe horizontal minimal construction coefficient of strain and horizontal maximum tectonic strain coefficient are represented respectively,
It can be determined by mini frac test and borehole wall stability analysis method, dimensionless;ξ represents hole bullet coefficient, generally takes
It is worth for 0, dimensionless.
The embodiment of the method for determining crustal stress, based on the attribute information, can establish the interval of interest and bore
The incidence relation of the SVEL in the vertical stratification face of well opening position and the SVEL of non-perpendicular bed plane;Based on the association
The SVEL of relation and the interval of interest in the vertical stratification face of drilling well opening position, it may be determined that the interval of interest is boring
The SVEL of the non-perpendicular bed plane of well opening position;Can be and described according to the density log curve of the interval of interest
Interval of interest determines the mesh in the SVEL in the vertical stratification face of drilling well opening position and the SVEL of non-perpendicular bed plane
Interval drilling well opening position minimum horizontal principal stress and maximum horizontal principal stress.Due to it is determined that minimum horizontal principal stress and
During maximum horizontal principal stress, without obtaining Stoneley wave, without approximate condition is introduced, so the side using the application
The degree of accuracy of crustal stress determined by method is higher.
Fig. 3 is the composition structure chart for the device embodiment that the application determines crustal stress.The device for determining crustal stress carries
For the attribute information of core sample in interval of interest, wherein, the attribute information can include:The SVEL in vertical stratification face
With the SVEL of non-perpendicular bed plane.As shown in figure 3, the device for determining crustal stress can include:Incidence relation determines
Module 100, interval SVEL determining module 200 and crustal stress determining module 300.
The incidence relation determining module 100, it can be used for being based on the attribute information, establish the interval of interest and boring
The incidence relation of the SVEL in the vertical stratification face of well opening position and the SVEL of non-perpendicular bed plane.
The interval SVEL determining module 200, it can be used for existing based on the incidence relation and the interval of interest
The SVEL in the vertical stratification face of drilling well opening position, determine non-perpendicular bed plane of the interval of interest in drilling well opening position
SVEL.
The crustal stress determining module 300, it can be used for the density log curve according to the interval of interest, and it is described
Interval of interest determines the mesh in the SVEL in the vertical stratification face of drilling well opening position and the SVEL of non-perpendicular bed plane
Interval drilling well opening position minimum horizontal principal stress and maximum horizontal principal stress.
The device embodiment for determining crustal stress is corresponding with the embodiment of the method for the determination crustal stress, it is possible to achieve
Determine the technical scheme of the embodiment of the method for crustal stress, and the technique effect of adquisitiones embodiment.
In the 1990s, the improvement for a technology can clearly distinguish be on hardware improvement (for example,
Improvement to circuit structures such as diode, transistor, switches) or software on improvement (improvement for method flow).So
And as the development of technology, the improvement of current many method flows can be considered as directly improving for hardware circuit.
Designer nearly all obtains corresponding hardware circuit by the way that improved method flow is programmed into hardware circuit.Cause
This, it cannot be said that the improvement of a method flow cannot be realized with hardware entities module.For example, PLD
(Programmable Logic Device, PLD) (such as field programmable gate array (Field Programmable Gate
Array, FPGA)) it is exactly such a integrated circuit, its logic function is determined by user to device programming.By designer
Voluntarily programming comes a digital display circuit " integrated " on a piece of PLD, without asking chip maker to design and make
Special IC chip.Moreover, nowadays, substitution manually makes IC chip, this programming is also used instead mostly " patrols
Volume compiler (logic compiler) " software realizes that software compiler used is similar when it writes with program development,
And the source code before compiling also write by handy specific programming language, this is referred to as hardware description language
(Hardware Description Language, HDL), and HDL is also not only a kind of, but have many kinds, such as ABEL
(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description
Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL
(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby
Hardware Description Language) etc., VHDL (Very-High-Speed are most generally used at present
Integrated Circuit Hardware Description Language) and Verilog2.Those skilled in the art
It will be apparent to the skilled artisan that only need method flow slightly programming in logic and being programmed into integrated circuit with above-mentioned several hardware description languages
In, it is possible to it is readily available the hardware circuit for realizing the logical method flow.
Controller can be implemented in any suitable manner, for example, controller can take such as microprocessor or processing
Device and storage can by the computer of the computer readable program code (such as software or firmware) of (micro-) computing device
Read medium, gate, switch, application specific integrated circuit (Application Specific Integrated Circuit,
ASIC), the form of programmable logic controller (PLC) and embedded microcontroller, the example of controller include but is not limited to following microcontroller
Device:ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20 and Silicone Labs C8051F320, are deposited
Memory controller is also implemented as a part for the control logic of memory.
It is also known in the art that in addition to realizing controller in a manner of pure computer readable program code, it is complete
Entirely can by by method and step carry out programming in logic come controller with gate, switch, application specific integrated circuit, may be programmed
The form of logic controller and embedded microcontroller etc. realizes identical function.Therefore this controller is considered one kind
Hardware component, and it is used to realize that the device of various functions can also to be considered as the structure in hardware component to what is included in it.Or
Even, it not only can be able to will be the software module of implementation method for realizing that the device of various functions is considered as but also can be Hardware Subdivision
Structure in part.
Device that above-described embodiment illustrates, module, it can specifically be realized by computer chip or entity, or by with certain
The product of kind of function is realized.
For convenience of description, it is divided into various modules during description apparatus above with function to describe respectively.Certainly, this is being implemented
The function of each module can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application
On the part that is contributed in other words to prior art can be embodied in the form of software product, in a typical configuration
In, computing device includes one or more processors (CPU), input/output interface, network interface and internal memory.The computer is soft
Part product can include some instructions make it that a computer equipment (can be personal computer, server, or network
Equipment etc.) perform method described in some parts of each embodiment of the application or embodiment.The computer software product can
To be stored in internal memory, internal memory may include the volatile memory in computer-readable medium, random access memory
(RAM) and/or the form such as Nonvolatile memory, such as read-only storage (ROM) or flash memory (flash RAM).Internal memory is computer
The example of computer-readable recording medium.Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by
Any method or technique come realize information store.Information can be computer-readable instruction, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only storage
(ROM), Electrically Erasable Read Only Memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic cassette tape, tape magnetic rigid disk storage or
Other magnetic storage apparatus or any other non-transmission medium, the information that can be accessed by a computing device available for storage.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation
Data-signal and carrier wave.
Each embodiment in this specification is described by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment stressed is the difference with other embodiment.It is real especially for device
For applying example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part is referring to embodiment of the method
Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as:Personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer including any of the above system or equipment
DCE etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by
Task is performed and connected remote processing devices by communication network.In a distributed computing environment, program module can be with
In the local and remote computer-readable storage medium including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without departing from spirit herein, it is desirable to which appended claim includes these deformations and changed without departing from the application's
Spirit.
Claims (10)
- A kind of 1. method for determining crustal stress, it is characterised in that the SVEL of core sample in interval of interest is provided with, it is described Method includes:Based on the SVEL of the core sample, sound wave of the interval of interest in the vertical stratification face of drilling well opening position is established Speed and the incidence relation of the SVEL of non-perpendicular bed plane;SVEL based on the incidence relation and the interval of interest in the vertical stratification face of drilling well opening position, it is determined that described SVEL of the interval of interest in the non-perpendicular bed plane of drilling well opening position;According to the density log curve of the interval of interest, and the interval of interest is in the vertical stratification face of drilling well opening position The SVEL of SVEL and non-perpendicular bed plane, determine minimum horizontal principal stress of the interval of interest in drilling well opening position And maximum horizontal principal stress.
- A kind of 2. method for determining crustal stress according to claim 1, it is characterised in that the sound wave speed of the core sample Degree includes the SVEL in vertical stratification face and the SVEL of non-perpendicular bed plane;Correspondingly, the SVEL based on the core sample, the interval of interest is established in the vertical of drilling well opening position The incidence relation of the SVEL of bed plane and the SVEL of non-perpendicular bed plane, including:The SVEL in the vertical stratification face based on the core sample and the SVEL of non-perpendicular bed plane, using a most young waiter in a wineshop or an inn The method for multiplying fitting establishes SVEL and non-perpendicular bed plane of the interval of interest in the vertical stratification face of drilling well opening position SVEL incidence relation.
- A kind of 3. method for determining crustal stress according to claim 2, it is characterised in that the sound wave in the vertical stratification face Speed includes:First shear wave velocity and the first velocity of longitudinal wave;Wherein, first shear wave velocity represents the shear wave in vertical stratification face Speed, first velocity of longitudinal wave represent the velocity of longitudinal wave in the vertical stratification face;The SVEL of the non-perpendicular bed plane Including:Second shear wave velocity, the second velocity of longitudinal wave and the 3rd velocity of longitudinal wave;Wherein, the second velocity of longitudinal wave expression and stratification The angle in face is 45 degree of velocity of longitudinal wave, the velocity of longitudinal wave that it is 0 degree with the angle of bed plane that the 3rd velocity of longitudinal wave, which represents,;Institute State the second shear wave velocity to represent the shear wave velocity for being 0 degree with the angle of bed plane;The first shear wave velocity and the second shear wave velocity based on the core sample, institute is established using the method for least square fitting Interval of interest is stated in the second shear wave velocity of drilling well opening position and the incidence relation of the first shear wave velocity;The first velocity of longitudinal wave and the second velocity of longitudinal wave based on the core sample, institute is established using the method for least square fitting Interval of interest is stated in the second velocity of longitudinal wave of drilling well opening position and the incidence relation of the first velocity of longitudinal wave;The first velocity of longitudinal wave and the 3rd velocity of longitudinal wave based on the core sample, institute is established using the method for least square fitting Interval of interest is stated in the 3rd velocity of longitudinal wave of drilling well opening position and the incidence relation of the first velocity of longitudinal wave.
- 4. a kind of method for determining crustal stress according to claim 1, it is characterised in that described according to the interval of interest Density log curve, and SVEL and non-perpendicular stratification of the interval of interest in the vertical stratification face of drilling well opening position The SVEL in face, minimum horizontal principal stress and maximum horizontal principal stress of the interval of interest in drilling well opening position are determined, wrapped Include:According to the density log curve, density of the interval of interest in drilling well opening position is determined;According to the interval of interest in the SVEL in the vertical stratification face of drilling well opening position, the sound wave speed of non-perpendicular bed plane Degree, determines at least one stiffness coefficient of the interval of interest in drilling well opening position;Wherein, the stiffness coefficient represents the mesh Interval corresponding to transverse isotropy petrophysical model stiffness matrix in stiffness coefficient;According to the interval of interest in the density and stiffness coefficient of drilling well opening position, determine the interval of interest in drilling well opening position Minimum horizontal principal stress and maximum horizontal principal stress.
- 5. a kind of method for determining crustal stress according to claim 4, it is characterised in that described according to the interval of interest In the density and stiffness coefficient of drilling well opening position, minimum horizontal principal stress and most of the interval of interest in drilling well opening position is determined Big horizontal principal stress, including:According to the interval of interest in the stiffness coefficient of drilling well opening position, dynamic of the interval of interest in drilling well opening position is determined Young's modulus and dynamic Poisson's ratio;Wherein, the kinetic Youngs modulus includes dynamic level Young's modulus and the vertical Young of dynamic Modulus, the dynamic Poisson's ratio include dynamic level Poisson's ratio and dynamic vertical Poisson's ratio;Determine static Young corresponding to kinetic Youngs modulus and dynamic Poisson's ratio difference of the interval of interest in drilling well opening position Modulus and static Poisson's ratio;Wherein, the static Young's modulus includes quiescent levels Young's modulus and static vertical Young's modulus, The static Poisson's ratio includes quiescent levels Poisson's ratio and static vertical Poisson's ratio;Determine Biot spy coefficient of the interval of interest in drilling well opening position;According to the interval of interest in the density of drilling well opening position, static Young's modulus, static Poisson's ratio and Biot spy's coefficient, really Minimum horizontal principal stress and maximum horizontal principal stress of the fixed interval of interest in drilling well opening position.
- 6. a kind of method for determining crustal stress according to claim 4, it is characterised in that described to determine the interval of interest Static Young's modulus and static Poisson's ratio corresponding to kinetic Youngs modulus and dynamic Poisson's ratio difference in drilling well opening position, bag Include:According to the SVEL in the vertical stratification face of the core sample and the SVEL of non-perpendicular bed plane, the mesh is established Interval in the kinetic Youngs modulus of drilling well opening position and the incidence relation of static Young's modulus, and establish the interval of interest In the dynamic Poisson's ratio of drilling well opening position and the incidence relation of static Poisson's ratio;According to the interval of interest in the kinetic Youngs modulus of drilling well opening position and the incidence relation of static Young's modulus, Yi Jidong The incidence relation of state Poisson's ratio and static Poisson's ratio, the interval of interest is determined in the kinetic Youngs modulus of drilling well opening position and is moved Static Young's modulus and static Poisson's ratio corresponding to state Poisson's ratio difference.
- 7. a kind of method for determining crustal stress according to claim 6, it is characterised in that described according to the core sample The SVEL in vertical stratification face and the SVEL of non-perpendicular bed plane, establish the interval of interest in drilling well opening position The incidence relation of kinetic Youngs modulus and static Young's modulus, and establish the interval of interest and moored in the dynamic of drilling well opening position The loose incidence relation than with static Poisson's ratio, including:According to the SVEL in the vertical stratification face of the core sample and the SVEL of non-perpendicular bed plane, the rock is determined The stiffness coefficient of heart sample;According to the stiffness coefficient of the core sample, the kinetic Youngs modulus and dynamic Poisson's ratio of the core sample are determined;Determine the static Young's modulus of the core sample and static Poisson's ratio;According to the static Young's modulus and kinetic Youngs modulus of the core sample, the interval of interest is established in drilling well opening position Kinetic Youngs modulus and static Young's modulus incidence relation, and static Poisson's ratio and dynamic according to the core sample Poisson's ratio, the interval of interest is established in the dynamic Poisson's ratio of drilling well opening position and the incidence relation of static Poisson's ratio.
- 8. a kind of method for determining crustal stress according to claim 4, it is characterised in that described to determine the interval of interest In Biot spy's coefficient of drilling well opening position, including:Determine Biot spy's coefficient of the core sample;Using the average value of Biot spy's coefficient of multiple core samples as the interval of interest drilling well opening position Biot Special coefficient.
- 9. a kind of method for determining crustal stress according to claim 1, it is characterised in that the core sample includes geometry It is shaped as the core sample of cylinder;Wherein, the axis of the cylinder is parallel or vertical with bed plane.
- 10. a kind of device for determining crustal stress, it is characterised in that described device provides the sound wave speed of core sample in interval of interest Degree;Described device includes:Incidence relation determining module, interval SVEL determining module and crustal stress determining module;Wherein,The incidence relation determining module, for the SVEL based on the core sample, establish the interval of interest and boring The incidence relation of the SVEL in the vertical stratification face of well opening position and the SVEL of non-perpendicular bed plane;The interval SVEL determining module, for based on the incidence relation and the interval of interest in drilling well opening position The SVEL in vertical stratification face, determine SVEL of the interval of interest in the non-perpendicular bed plane of drilling well opening position;The crustal stress determining module, exist for the density log curve according to the interval of interest, and the interval of interest The SVEL in the vertical stratification face of drilling well opening position and the SVEL of non-perpendicular bed plane, determine that the interval of interest is boring The minimum horizontal principal stress and maximum horizontal principal stress of well opening position.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710817354.9A CN107829727A (en) | 2017-09-12 | 2017-09-12 | A kind of method and device for determining crustal stress |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710817354.9A CN107829727A (en) | 2017-09-12 | 2017-09-12 | A kind of method and device for determining crustal stress |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107829727A true CN107829727A (en) | 2018-03-23 |
Family
ID=61643810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710817354.9A Pending CN107829727A (en) | 2017-09-12 | 2017-09-12 | A kind of method and device for determining crustal stress |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107829727A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110320277A (en) * | 2018-03-30 | 2019-10-11 | 中国石油天然气股份有限公司 | In-situ stress state test method |
CN111460602A (en) * | 2019-12-06 | 2020-07-28 | 西南石油大学 | Transverse isotropic stratum ground stress prediction method based on rock physics modeling |
CN111579405A (en) * | 2020-05-20 | 2020-08-25 | 河南工程学院 | Multifunctional rock test system |
CN114200538A (en) * | 2020-09-17 | 2022-03-18 | 中国石油化工股份有限公司 | Ground stress direction determining method and device, storage medium and computer equipment |
CN117189092A (en) * | 2023-08-16 | 2023-12-08 | 中国矿业大学 | Soft rock ground stress testing method based on drilling cuttings particle size distribution |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1242746A1 (en) * | 1984-12-10 | 1986-07-07 | Военная Ордена Ленина И Ордена Октябрьской Революции Краснознаменная Академия Бронетанковых Войск Им.Маршала Советского Союза Малиновского Р.Я. | Installation for studying stress and shifts of ground under vehicle supports |
WO2009035869A1 (en) * | 2007-09-11 | 2009-03-19 | Schlumberger Canada Limited | Joint inversion of borehole acoustic radial profiles for in situ stresses as well as third-order nonlinear dynamic moduli, linear dynamic elastic moduli, and static elastic moduli in an isotropically stressed reference state |
CN103792581A (en) * | 2013-12-31 | 2014-05-14 | 刘致水 | Shale gas reservoir crustal stress logging prediction method based on rock physics model |
CN104975836A (en) * | 2015-06-30 | 2015-10-14 | 中国石油天然气股份有限公司 | Rock sample hydraulic fracture morphology acoustic emission diagnosis experimental method and device |
CN106053231A (en) * | 2016-07-18 | 2016-10-26 | 西南石油大学 | Testing device for anisotropism of shale in true-triaxial condition and testing method of testing device |
CN106295119A (en) * | 2016-07-26 | 2017-01-04 | 中国石油大学(华东) | A kind of shale gas stratum crustal stress computational methods |
CN106353822A (en) * | 2016-08-05 | 2017-01-25 | 成都理工大学 | Shale mechanical parameter logging method based on VTI transverse isotropy model |
CN106501872A (en) * | 2016-09-26 | 2017-03-15 | 中国石油天然气股份有限公司 | A kind of computational methods of fracture reservoir ground stress characteristics and device |
CN206177657U (en) * | 2016-08-05 | 2017-05-17 | 成都理工大学 | Shale mechanics parameter testing sample |
CN106855897A (en) * | 2016-11-25 | 2017-06-16 | 长江大学 | Suitable for the research method of the wellbore stability on pressure exhaustion stratum |
CN106855485A (en) * | 2016-12-20 | 2017-06-16 | 中国石油天然气股份有限公司 | A kind of conversion method of sound state elastic parameter |
-
2017
- 2017-09-12 CN CN201710817354.9A patent/CN107829727A/en active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1242746A1 (en) * | 1984-12-10 | 1986-07-07 | Военная Ордена Ленина И Ордена Октябрьской Революции Краснознаменная Академия Бронетанковых Войск Им.Маршала Советского Союза Малиновского Р.Я. | Installation for studying stress and shifts of ground under vehicle supports |
WO2009035869A1 (en) * | 2007-09-11 | 2009-03-19 | Schlumberger Canada Limited | Joint inversion of borehole acoustic radial profiles for in situ stresses as well as third-order nonlinear dynamic moduli, linear dynamic elastic moduli, and static elastic moduli in an isotropically stressed reference state |
CN103792581A (en) * | 2013-12-31 | 2014-05-14 | 刘致水 | Shale gas reservoir crustal stress logging prediction method based on rock physics model |
CN104975836A (en) * | 2015-06-30 | 2015-10-14 | 中国石油天然气股份有限公司 | Rock sample hydraulic fracture morphology acoustic emission diagnosis experimental method and device |
CN106053231A (en) * | 2016-07-18 | 2016-10-26 | 西南石油大学 | Testing device for anisotropism of shale in true-triaxial condition and testing method of testing device |
CN106295119A (en) * | 2016-07-26 | 2017-01-04 | 中国石油大学(华东) | A kind of shale gas stratum crustal stress computational methods |
CN106353822A (en) * | 2016-08-05 | 2017-01-25 | 成都理工大学 | Shale mechanical parameter logging method based on VTI transverse isotropy model |
CN206177657U (en) * | 2016-08-05 | 2017-05-17 | 成都理工大学 | Shale mechanics parameter testing sample |
CN106501872A (en) * | 2016-09-26 | 2017-03-15 | 中国石油天然气股份有限公司 | A kind of computational methods of fracture reservoir ground stress characteristics and device |
CN106855897A (en) * | 2016-11-25 | 2017-06-16 | 长江大学 | Suitable for the research method of the wellbore stability on pressure exhaustion stratum |
CN106855485A (en) * | 2016-12-20 | 2017-06-16 | 中国石油天然气股份有限公司 | A kind of conversion method of sound state elastic parameter |
Non-Patent Citations (1)
Title |
---|
翁定为等: "致密油水平井分段多簇优化设计方法", 《中国石油大学学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110320277A (en) * | 2018-03-30 | 2019-10-11 | 中国石油天然气股份有限公司 | In-situ stress state test method |
CN110320277B (en) * | 2018-03-30 | 2021-08-03 | 中国石油天然气股份有限公司 | Ground stress field state testing method |
CN111460602A (en) * | 2019-12-06 | 2020-07-28 | 西南石油大学 | Transverse isotropic stratum ground stress prediction method based on rock physics modeling |
CN111579405A (en) * | 2020-05-20 | 2020-08-25 | 河南工程学院 | Multifunctional rock test system |
CN114200538A (en) * | 2020-09-17 | 2022-03-18 | 中国石油化工股份有限公司 | Ground stress direction determining method and device, storage medium and computer equipment |
CN117189092A (en) * | 2023-08-16 | 2023-12-08 | 中国矿业大学 | Soft rock ground stress testing method based on drilling cuttings particle size distribution |
CN117189092B (en) * | 2023-08-16 | 2024-04-09 | 中国矿业大学 | Soft rock ground stress testing method based on drilling cuttings particle size distribution |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107829727A (en) | A kind of method and device for determining crustal stress | |
CN111665560B (en) | Hydrocarbon reservoir identification method, apparatus, computer device and readable storage medium | |
US20210304060A1 (en) | Method and system for enhancing artificial intelligence predictions using well data augmentation | |
AU2014384715B2 (en) | Facies definition using unsupervised classification procedures | |
CN108150161A (en) | Shale gassiness evaluation square law device | |
CN106153495B (en) | Rock gas adsorbs phase volume measurement method and device | |
CN106896406A (en) | The method and apparatus that High-quality Reservoir is predicted based on impedance Domain Properties | |
CN108019205A (en) | A kind of method and device of definite reservoir compressibility index | |
JP2020527725A (en) | Methods and systems for determining bulk density, porosity, and porosity distribution of surface formations | |
Liu et al. | Asymmetric propagation mechanism of hydraulic fracture networks in continental reservoirs | |
CN107367762A (en) | A kind of method and device for determining reservoir parameter | |
CN107817520A (en) | The pressure coefficient Forecasting Methodology and system of marine facies mud shale stratum | |
Soares et al. | An experimental study for mechanical formation damage | |
CN107942404A (en) | A kind of method and device in definite crack | |
CN108035709A (en) | A kind of definite method and device of shale reservoir quality | |
CN108387931B (en) | The determination method and apparatus of reservoir pore pressure | |
Askari et al. | A fully integrated method for dynamic rock type characterization development in one of Iranian off-shore oil reservoir | |
CN106873030A (en) | The determination method and apparatus of Fault Seal | |
CN106202813A (en) | A kind of method and device calculating oil viscosity | |
US11525352B2 (en) | Method and system to automate formation top selection using well logs | |
CN106772599B (en) | A kind of method and device calculating formation shear speed | |
CN109630090A (en) | Store up lid combination logging evaluation method and device | |
US11391683B2 (en) | Method to determine pore size distribution of rocks with rough surface from capillary pressure by nuclear magnetic resonance | |
CN106682757A (en) | Method and device for determining collapse pressure of borehole wall | |
CN108387711A (en) | The determination method and apparatus of TOC characterization parameters, brittleness characterization parameter |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180323 |