CN103257081A - Method and device for recovering oil and gas reservoir rock mechanics underground in-situ model - Google Patents

Method and device for recovering oil and gas reservoir rock mechanics underground in-situ model Download PDF

Info

Publication number
CN103257081A
CN103257081A CN201310157416XA CN201310157416A CN103257081A CN 103257081 A CN103257081 A CN 103257081A CN 201310157416X A CN201310157416X A CN 201310157416XA CN 201310157416 A CN201310157416 A CN 201310157416A CN 103257081 A CN103257081 A CN 103257081A
Authority
CN
China
Prior art keywords
rock
parameter
mechanics
poisson ratio
anisotropy
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
Application number
CN201310157416XA
Other languages
Chinese (zh)
Other versions
CN103257081B (en
Inventor
师永民
张玉广
王磊
吴文娟
柴智
秦小双
师锋
李晓敏
郭馨蔚
吴洛菲
熊文涛
徐蕾
师春爱
方媛媛
师翔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking University
Original Assignee
Peking University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Peking University filed Critical Peking University
Priority to CN201310157416.XA priority Critical patent/CN103257081B/en
Publication of CN103257081A publication Critical patent/CN103257081A/en
Application granted granted Critical
Publication of CN103257081B publication Critical patent/CN103257081B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention provides a method and a device for recovering an oil and gas reservoir rock mechanics underground in-situ model. The method comprises the following steps of: measuring rock mechanics parameters of an oil and gas reservoir; counting rock mechanics property influence parameters of the oil and gas reservoir; calculating the anisotropy strength of the rock mechanics parameters; obtaining a correlation between the rock mechanics parameters and the rock mechanics property influence parameters of the oil and gas reservoir; obtaining a correlation between the rock mechanics property influence parameters and the anisotropy strength; and recovering the oil and gas reservoir rock mechanics underground in-situ model by adopting the rock mechanics parameters, the anisotropy strength of the rock mechanics parameters, the correlation between the rock mechanics parameters and the rock mechanics property influence parameters of the oil and gas reservoir and the correlation of the rock mechanics property influence parameters and the anisotropy strength. By the method, the effectiveness and the accuracy for recovering the oil and gas reservoir rock mechanics underground in-situ model can be improved.

Description

Method and device that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics recovers
Technical field
The present invention relates to hydrocarbon-bearing pool development technique field, particularly relate to the method that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics recovers, and, the device that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics recovers.
Background technology
The hydrocarbon-bearing pool rock mass be preserve the rock mass porous media material and up and down country rock under certain depth, be under uniform temperature, the pressure, saturated some fluids the geology entity of being formed by a series of structural planes and structure.Compare with the near surface engineering geology, generally bury dark (the location layer depth is generally thousands of rice), be under higher three complicated confined pressures (can reach 200MPa), higher temperature (can reach 200 ° of C) and the higher pore pressure (can reach 200Mpa) and heterogeneous fluid saturated (oil, gas, water) condition.The residing environment of rock mass at this and shallow surface and lower crust upper mantle place has a great difference, and rock mass also has the unique mechanical characteristic.
For the complicated confined pressure in this underground thousands of meters depths, high temperature, the hydrocarbon-bearing pool rock mass of high pore pressure and heterogeneous fluid, traditional static mechanical property research be with drilling and coring delivery under surface condition, utilize triaxial stress tester or sonic test instrument that the rock core that takes out is carried out Experiments of Machanics, as drilling and coring delivery being taken normal temperature in the ground experiment chamber, normal pressure, dry sample, or consideration applies different confined pressures, or consider to fill different saturated fluids and carry out the rock mechanics experiment, obviously, the rock mechanics parameter that adopts this prior art to obtain can not represent fully and be in underground high temperature, high pressure, high pore pressure, hydrocarbon-bearing pool rock mass mechanics character under the heterogeneous saturated flow concrete conditions in the establishment of a specific crime namely can't guarantee the complicated confined pressure in underground thousands of meters depths, high temperature, validity and the accuracy of the hydrocarbon-bearing pool rock mass mechanics parameter of high pore pressure and heterogeneous fluid.
Summary of the invention
Technical matters to be solved by this invention provides the method and apparatus that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics recovers, in order to improve validity and the accuracy that the underground former bit model of hydrocarbon-bearing pool rock mass mechanics recovers.
In order to address the above problem, the invention discloses the method that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics recovers, comprising:
Measure the rock mechanics parameter of hydrocarbon-bearing pool;
Add up the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool;
Calculate the strength of anisotropy of described rock mechanics parameter;
Obtain the incidence relation of rock mechanics parameter and the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool;
Obtain the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy;
Adopt described rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and, the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
Preferably, the rock mechanics parameter comprises Young modulus, Poisson ratio and Biot coefficient, and described Biot coefficient is used for characterizing the poroelasticity parameter; The strength of anisotropy of described rock mechanics parameter comprises Young modulus strength of anisotropy, Poisson ratio strength of anisotropy and Biot coefficient strength of anisotropy.
Preferably, described rock mass mechanics property effect parameter comprises the external action parameter, and the described step of obtaining the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy is,
Obtain the incidence relation of described external action parameter and strength of anisotropy.
Preferably, the strength of anisotropy of described rock mechanics parameter is calculated by following formula and is obtained:
E k = E h - E v E h ;
μ k = μ h - μ v μ h ;
Wherein, E kBe the Young modulus strength of anisotropy; E hBe the horizontal direction Young modulus; E vBe the vertical direction Young modulus; μ kBe the Poisson ratio strength of anisotropy; μ hBe the horizontal direction Poisson ratio; μ vBe the vertical direction Poisson ratio.
Preferably, described external action parameter comprises the temperature profile parameter, confined pressure characteristic parameter, pore fluid characteristic parameter, pore pressure characteristic parameter; The incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy comprises:
The strength of anisotropy of described rock mechanics parameter becomes big along with the increase of temperature profile parameter;
And/or,
The strength of anisotropy of described rock mechanics parameter reduces along with the increase of confined pressure characteristic parameter;
And/or,
The strength of anisotropy of described rock mechanics parameter reduces along with the increase of pore fluid characteristic parameter; Wherein, described pore fluid characteristic parameter is fluid saturation;
And/or,
The strength of anisotropy of described rock mechanics parameter becomes big along with the increase of pore pressure characteristic parameter.
Preferably, the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool comprises the internal influence parameter, and described internal influence parameter comprises the rock composition characteristic parameter, rock texture characteristic parameter, rock structure characteristic parameter, blowhole FRACTURE CHARACTERISTICS parameter;
Described rock mechanics parameter comprises rock elasticity parameter and rock strength parameter, and described rock elasticity parameter comprises Poisson ratio, Young modulus;
Described rock composition characteristic parameter comprises rock particles characteristic parameter and chink characteristic parameter, and the incidence relation of the rock mechanics parameter of described hydrocarbon-bearing pool and rock mass mechanics property effect parameter comprises:
If quartz content increases in the described rock particles characteristic parameter, then Young modulus rises;
And/or,
If described rock particles characteristic parameter andesine content increases, then Young modulus reduces;
And/or,
If described rock particles characteristic parameter CLAY MINERALS AND THEIR SIGNIFICANCE content increases, then Poisson ratio reduces;
And/or,
If described chink characteristic parameter is the calcitic cementation rock, then its Young modulus is greater than the Young modulus of other type of rock;
And/or,
If the rock that described rock texture characteristic parameter is the crystal class origin cause of formation, then its Young modulus and Poisson ratio are higher than the rock that the rock texture characteristic parameter is the cementing matter diagenesis;
And/or,
If described rock structure characteristic parameter bored sample by horizontal direction, its Young modulus and Poisson ratio Young modulus and the Poisson ratio of boring sample greater than vertical direction then;
And/or,
If blowhole FRACTURE CHARACTERISTICS parameter is the hole fracture spaces, then along with the increase of hole fracture spaces, the Young modulus of rock and the value of Poisson ratio reduce.
Preferably, the incidence relation of the rock mechanics parameter of described hydrocarbon-bearing pool and rock mass mechanics property effect parameter also comprises:
If described temperature profile parameter increases, then the hole fracture density of rock increases, and the Young modulus of rock reduces, and the Poisson ratio of rock reduces, and the poroelasticity coefficient of rock increases;
And/or,
If described confined pressure characteristic parameter increases, then the hole fracture density of rock reduces, and the Young modulus of rock and Poisson ratio increase, and the poroelasticity coefficient of rock reduces;
And/or,
If described pore fluid characteristic parameter is fluid saturation, then along with the increase of fluid saturation, the Young modulus of rock, Poisson ratio and poroelasticity coefficient all increase; And, the Poisson ratio the when Poisson ratio in rock core oil-containing period is moisture greater than rock core, along with the increase of fluid saturation, its Poisson ratio difference increases;
And/or,
If described pore pressure characteristic parameter increases, then the Young modulus of rock, Poisson ratio reduce the increase of poroelasticity coefficient.
Preferably, the step of the rock mechanics parameter of described measurement hydrocarbon-bearing pool comprises:
Adopt static measurement method to obtain static rock mechanics parameter;
Adopt dynamic measurement method to obtain dynamic rock mechanics parameter;
Wherein, described static rock mechanics parameter and dynamic rock mechanics parameter obtain through synchro measure;
Described method also comprises:
Adopt described static rock mechanics parameter and dynamic rock mechanics parameter, the sound diversity ratio of computing rock mechanics parameter;
Described step of carrying out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics according to the rock mechanics parameter after adjusting is,
Adopt described static rock mechanics parameter, dynamic rock mechanics parameter, the sound diversity ratio of rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
Preferably, described static rock mechanics parameter comprises static rock elasticity parameter and static rock strength parameter, and described dynamic rock mechanics parameter comprises dynamic rock elasticity parameter and dynamic rock strength parameter; Described static rock elasticity parameter comprises static Poisson ratio, static Young modulus and static poroelasticity coefficient, and described dynamic rock elasticity parameter comprises dynamic Poisson ratio, kinetic Youngs modulus and dynamic hole elasticity coefficient.
Preferably, described static Poisson ratio μ calculates by following formula and obtains:
μ = Δϵ 1 Δϵ 2
Described dynamic Poisson ratio μ calculates by following formula and obtains:
μ = Δt s 2 - 2 Δt p 2 2 ( Δt s 2 - Δt p 2 ) = V p 2 - 2 V s 2 2 ( V p 2 - V s 2 )
Wherein, described Δ ε 1Axial strain increment for rock sample; Δ ε 2Transverse strain increment for rock sample; Δ t sBe the shear wave slowness of rock sample, μ s/ m; Δ t pBe the compressional wave time difference of rock sample, μ s/ m; V sBe the shear wave velocity of rock sample, m/s; V pBe the velocity of longitudinal wave of rock sample, m/s;
Described static Young modulus E calculates by following formula and obtains:
E = F / S ΔL / L
Described kinetic Youngs modulus E calculates by following formula and obtains:
E = ρ V s 2 3 ( V p V s ) 2 - 4 ( V p V s ) 2 - 1
Wherein, described ρ is the density of rock, g/cm 3F is the suffered power of test rock sample; S is the sectional area of rock sample; Δ L is that rock sample is being subjected to the elongation of force direction (or shortening amount); L is that rock sample is in the former length that is subjected on the force direction; V SBe the shear wave velocity of rock sample, m/s; Vp is the velocity of longitudinal wave of rock sample, m/s;
Described static poroelasticity coefficient B is calculated by following formula and is obtained:
B = 1 - K dry K s
Wherein, described K DryBe the bulk modulus of dry rock, Ks is the bulk modulus of the mineral of composition rock;
Described dynamic hole elasticity coefficient B calculates by following formula and obtains:
B = 1 - K dry K s = 1 - 1 ΔV s 1 ΔV s - ΔV f = 1 - ΔV s - ΔV f ΔV s = ΔV f ΔV s
Wherein, described K DryBe the bulk modulus of dry rock, Ks is the bulk modulus of the mineral of composition rock, Δ V fBe pore fluid fluctuations in discharge amount, Δ V sBe the rock skeleton volume change, Δ V mBe rock forming mineral particle volume variable quantity;
The sound diversity ratio of described rock mechanics parameter comprises Young modulus sound diversity ratio and Poisson ratio sound diversity ratio, and described Young modulus sound diversity ratio is calculated by following formula and obtained:
E c = E d - E s E s
Wherein, described E dBe kinetic Youngs modulus; E sBe static Young modulus;
Described Poisson ratio sound diversity ratio is calculated by following formula and is obtained:
μ c = μ d - μ s μ s ;
Wherein, described μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio.
Preferably, satisfy the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the temperature profile parameter:
E d - E s E s = aT 2 + bT + C ;
μ d - μ s μ s = mT 2 + nT + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; T is temperature; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant;
Satisfy the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the confined pressure characteristic parameter:
E d - E s E s = aP c 2 + bP c + C ;
μ d - μ s μ s = m P c 2 + nP c + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; P cBe confined pressure; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant;
Satisfy the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the fluid saturation:
E d - E s E s = a 1 S w 2 + b 1 S w + C 1 ;
μ d - μ s μ s = m 1 S w 2 + n 1 S w + k 1 ;
E d - E s E s = a 2 S o 2 + b 2 S o + C 2 ;
μ d - μ s μ s = m 2 S o 2 + n 2 S o + k 2 ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; S wBe water saturation; S oBe oil saturation; a 1, a 2Be the quadratic term coefficient; b 1, b 2Be coefficient once; C 1, C 2Be constant; m 1, m 2Be the quadratic term coefficient; n 1, n 2Be coefficient once; k 1, k 2Be constant;
Satisfy binomial relation shown in following formula between the sound diversity ratio of described rock mechanics parameter and the pore pressure:
E d - E s E s = a P o 2 + bP o + C ;
μ d - μ s μ s = m P o 2 + nP o + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; P oBe pore pressure; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant.
The embodiment of the invention also discloses the device that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics recovers, comprising:
Rock mechanics parameter measurement module is used for measuring the rock mechanics parameter of hydrocarbon-bearing pool;
Influence the parametric statistics module, be used for the rock mass mechanics property effect parameter of the described hydrocarbon-bearing pool of statistics;
The strength of anisotropy computing module is for the strength of anisotropy of calculating described rock mechanics parameter;
First concerns acquisition module, for the incidence relation of rock mechanics parameter and the rock mass mechanics property effect parameter of obtaining described hydrocarbon-bearing pool;
Second concerns acquisition module, is used for obtaining the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy;
Former bit model recovers module, be used for adopting described rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and, the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
Compared with prior art, the present invention includes following advantage:
The present invention amplifies concept, theory and the method for rock mass mechanics in the near surface engineering geology in the hydrocarbon-bearing pool rock mass under the complicated confined pressure in underground thousands of meters depths, high temperature, high pore pressure and the heterogeneous fluid environment, has formed the concept of hydrocarbon-bearing pool rock mass mechanics.Mechanical property that to be the underground reservoir rock mass show under the specific degree of depth, temperature, three confined pressures, pore pressure and heterogeneous fluid saturation conditionss that the hydrocarbon-bearing pool rock mass mechanics is described.The influence of external environment conditions such as build-in attributes such as condition and the degree of depth, temperature, confined pressure, pore pressure, heterogeneous fluid is oozed in rock composition, structure, structure, the hole that the hydrocarbon-bearing pool rock mass mechanics mainly is subjected to structural plane and structure inside.
Be drilling and coring delivery to be taken the ground experiment chamber apply different confined pressures or consider that saturated dissimilar fluid carries out rock mechanics experiment, seldom considers the influence of pore pressure to the traditional static mechanical property research of hydrocarbon-bearing pool rock mass.Consider that simultaneously the rock mass mechanics test under buried depth, temperature, confined pressure, pore pressure and the saturated flow concrete conditions in the establishment of a specific crime is not open in the prior art.The underground original position static mechanical parameter field of hydrocarbon-bearing pool requirement of engineering generally is to obtain dynamic mechanics parameter by well logging, earthquake, obtains indirectly through dynamic and static correction again.Underground former bit model static mechanical character and dynamic and static difference correction problem are blank of studying both at home and abroad.To this, the present invention has creatively proposed to load temperature, three confined pressures, pore pressure and saturated fluids simultaneously under laboratory condition and has carried out the dynamic and static mechanical property test of rock mass, set up dynamic and static correction relationship formula, and then obtain the scheme of the underground former bit model static mechanical character of hydrocarbon-bearing pool rock mass mechanics.
Innovative point of the present invention relates to:
1) under laboratory condition, loads simultaneously under temperature, three confined pressures, pore pressure and saturated flow concrete conditions in the establishment of a specific crime and carry out the rock mass mechanics property detection, and then recover the mechanical property of hydrocarbon-bearing pool rock mass under underground former bit model, analyze different factors to the influence of mechanical property;
2) rock under the different experimental conditions has been carried out dynamically and the synchronous mensuration of static mechanical parameter, analyzed different factors to the influence of sound difference, set up dynamic and static correction relationship formula;
3) analyzed and researched rock composition, structure, construct and preserve rock build-in attributes such as characteristic to the influence of mechanical property;
4) qualitative analysis the influence to rock mechanics parameter strength of anisotropy of confined pressure, fluid, pore pressure and temperature.
The former bit model of hydrocarbon-bearing pool rock mass mechanics proposed by the invention recovers research, can guarantee validity and the accuracy of the hydrocarbon-bearing pool rock mass mechanics parameter of the complicated confined pressure in underground thousands of meters depths, high temperature, high pore pressure and heterogeneous fluid, and help applied geophysics (well logging, the earthquake) raising of theory and technology, the optimization of petroleum engineering scheme (drilling well, pressure break, oil recovery etc.) and the raising of natural fracture precision of prediction etc., can also improve oil-gas exploration and development efficiency.
Description of drawings
Fig. 1 is the flow chart of steps of the method embodiment 1 of the underground former bit model recovery of a kind of hydrocarbon-bearing pool rock mass mechanics of the present invention;
Fig. 2 is the flow chart of steps of the method embodiment 2 of the underground former bit model recovery of a kind of hydrocarbon-bearing pool rock mass mechanics of the present invention;
Fig. 3 is the stress-strain-time relation curve map of rock in a kind of example of the present invention;
Fig. 4 is Biot coefficient test experiments synoptic diagram in a kind of example of the present invention;
Fig. 5 is core sampling synoptic diagram in a kind of example of the present invention;
Fig. 6 is mechanical properties of rock and the mineral content curve synoptic diagram that crosses in a kind of example of the present invention;
Fig. 7 is rock sample Young modulus and quartz content histogram in a kind of example of the present invention;
Fig. 8 is that the acoustie emission event number varies with temperature histogram in a kind of example of the present invention;
Fig. 9 be in a kind of example of the present invention Poisson ratio with different fluid type change curve synoptic diagram;
Figure 10 is the structured flowchart of the device embodiment of the underground former bit model recovery of a kind of hydrocarbon-bearing pool rock mass mechanics of the present invention.
Embodiment
For above-mentioned purpose of the present invention, feature and advantage can be become apparent more, the present invention is further detailed explanation below in conjunction with the drawings and specific embodiments.
One of core idea of the present invention is, be research object with rock elasticity mechanics parameter commonly used in the hydrocarbon-bearing pool engineering, under the condition that the hydrocarbon-bearing pool prototype is recovered, from the funtcional relationship of rock mass mechanics parameter and confined pressure, pore pressure, temperature and pore fluid, difference and two aspects of influence factor thereof of rock mass mechanics parameter dynamic and static measurement result, studied the influence to the rock mechanics parameter of The structure of rock, structure, composition and hole fissured structure; The variation of hydrocarbon-bearing pool rock mass power mechanics elastic parameter and the relation between confined pressure, pore pressure, temperature and the pore fluid; Difference and influence factor thereof between triaxial stress test (static test) result and sonic test (dynamic test) result, the recurrence of the relationship between the sound parameter; Influence the rock mass mechanics elastic parameter each to microscopic mechanism of abnormal sound power, and and confined pressure, temperature, pore pressure and fluid between relation.Thereby macroscopical experimental phenomena and micro-analysis are combined, analyze composition, structure and the hole FRACTURE CHARACTERISTICS of rock to the influence of rock mechanics parameter, obtain influencing the microscopic mechanism of rock mechanics parameter character; Set up the method for the former bit model of indoor recovery hydrocarbon-bearing pool rock mass, by the influence to the rock mass mechanics parameter of multiple spot experiment quantitative test temperature, confined pressure, pore pressure and pore fluid; And by the result difference between quantitative test rock mass mechanics parameter dynamic measurement method and the static measurement method and the reason that makes a difference, and provide conversion formula between the sound attitude parameter, and the further anisotropic influence factor of qualitative analysis rock mass, provide the anisotropic Changing Pattern of rock mass mechanics parameter under different temperatures, the confined pressure, thereby guarantee validity and the accuracy of the hydrocarbon-bearing pool rock mass mechanics parameter of the complicated confined pressure in underground thousands of meters depths, high temperature, high pore pressure and heterogeneous fluid.
With reference to figure 1, show the flow chart of steps of the method embodiment 1 that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics of the application recovers, specifically can comprise the steps:
Step 101 is measured the rock mechanics parameter of hydrocarbon-bearing pool;
In a preferred embodiment of the present invention, described rock mechanics parameter can comprise: rock elasticity parameter and rock strength parameter, and wherein, described rock elasticity parameter can comprise Poisson ratio, Young modulus and be used for to characterize the Biot coefficient of poroelasticity parameter; Described rock composition characteristic parameter can comprise rock particles characteristic parameter and chink characteristic parameter etc.
Step 102, the rock mass mechanics property effect parameter of adding up described hydrocarbon-bearing pool;
In a preferred embodiment of the present invention, the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool can comprise the internal influence parameter, and described internal influence parameter can comprise: rock composition characteristic parameter, rock texture characteristic parameter, the rock structure characteristic parameter, blowhole FRACTURE CHARACTERISTICS parameter.
More preferably, the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool can also comprise the external action parameter; Described external action parameter can comprise: temperature profile parameter, confined pressure characteristic parameter, pore fluid characteristic parameter, pore pressure characteristic parameter.
Step 103 is calculated the strength of anisotropy of described rock mechanics parameter;
The anisotropy of mechanical properties of rock is one of important physical property of rock.Velocity anisotropy in the seismic wave propagation process injects the non-homogeneous of water and involves in the oilfield development process, the stability of the borehole wall, the success ratio of pressing crack construction, and the successful Application of highly difficult orientation well etc. all is subjected to the anisotropic influence of mechanical properties of rock.
The rock that reservoir exploration and exploitation are studied mostly is sea, river facies sedimentogeneous rock.Because the stratum is formed by deposition repeatedly often, therefore, its bedding on longitudinal profile is very tangible.This also be cause sedimentogeneous rock in the vertical direction with horizontal direction on mechanical property have the major reason of notable difference.What it should be noted that sea, river mobilely has a stronger directivity.Make irregular rock particles in deposition process, arrange with distribution and also have certain directivity.Velocity anisotropy in the seismic wave propagation process injects the non-homogeneous of water and involves in the oilfield development process, the stability of the borehole wall, the success ratio of pressing crack construction, and the successful Application of highly difficult orientation well etc. all is subjected to the anisotropic influence of mechanical properties of rock.
Anisotropy is the important physical parameter of rock.In the various factors of petrologic make-up, the lattice advantage is orientated and hangs down the microfissure that aligns under the confined pressure is the topmost factor that determines the rock macroscopic anisotropy.Anisotropy mainly is lattice advantage orientation (preferential direction) decision by mineral composition under the high confining pressure.Usually the preferential direction of mineral crystal is formed by the arrangement of anisometric crustal particle in local rule growth, the mobile mediophyric wire arrangement of magma and the sedimentogeneous rock in the rock.
For a long time, for the Study on influencing factors deficiency of rock anisotropy intensity, inventor herein's method has by experiment been carried out Quantitative Study to confined pressure, pore pressure, fluid, temperature to the influence of rock mechanics strength of anisotropy.The sandstone that prolongs oil field prolongation group is chosen in experiment, and drilling through diameter at vertical stratification and parallel bedding both direction respectively is 19.7mm, and length is that 50mm and diameter are 25mm, and length is two kinds of cylindrical samples of 62.5mm.Earlier the sample that processes is put into drying box, oven dry 140h is as dry sample about 48 ° of C.Rock core is divided into four groups to experimentize: first group is the experiment of alternating temperature degree, and second group is to become the confined pressure experiment, and the 3rd group is the experiment of fluid saturation degree, and the 4th group is the experiment of change pore pressure.Its key step is as follows: the collection of (1) on-the-spot big rock core, adopt from prolonging the big tube rock core that prolongation group in oil field drills through; (2) big rock core axially and radially core respectively; (3) title of mark rock core and the direction of coring.
In a preferred embodiment of the present invention, the strength of anisotropy of described rock mechanics parameter comprises Young modulus strength of anisotropy, Poisson ratio strength of anisotropy and Biot coefficient strength of anisotropy, specifically can adopt following formula to calculate:
E k = E h - E v E h ;
μ k = μ h - μ v μ h ;
Wherein, E kBe the Young modulus strength of anisotropy; E hBe the horizontal direction Young modulus; E vBe the vertical direction Young modulus; μ kBe the Poisson ratio strength of anisotropy; μ hBe the horizontal direction Poisson ratio; μ vBe the vertical direction Poisson ratio.
Step 104 is obtained the incidence relation of rock mechanics parameter and the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool;
Rock mass mechanics property effect parameter in described hydrocarbon-bearing pool comprises under the situation of internal influence parameter that the incidence relation of the rock mechanics parameter of described hydrocarbon-bearing pool and rock mass mechanics property effect parameter can comprise:
If quartz content increases in the described rock particles characteristic parameter, then Young modulus rises;
And/or,
If described rock particles characteristic parameter andesine content increases, then Young modulus reduces;
And/or,
If described rock particles characteristic parameter CLAY MINERALS AND THEIR SIGNIFICANCE content increases, then Poisson ratio reduces;
And/or,
If described chink characteristic parameter is the calcitic cementation rock, then its Young modulus is greater than the Young modulus of other type of rock;
And/or,
If the rock that described rock texture characteristic parameter is the crystal class origin cause of formation, then its Young modulus and Poisson ratio are higher than the rock that the rock texture characteristic parameter is the cementing matter diagenesis;
And/or,
If described rock structure characteristic parameter bored sample by horizontal direction, its Young modulus and Poisson ratio Young modulus and the Poisson ratio of boring sample greater than vertical direction then;
And/or,
If blowhole FRACTURE CHARACTERISTICS parameter is the hole fracture spaces, then along with the increase of hole fracture spaces, the Young modulus of rock and the value of Poisson ratio reduce.
Rock mass mechanics property effect parameter in described hydrocarbon-bearing pool comprises under the situation of external action parameter that the incidence relation of the rock mechanics parameter of described hydrocarbon-bearing pool and rock mass mechanics property effect parameter can comprise:
If described temperature profile parameter increases, then the hole fracture density of rock increases, and the Young modulus of rock reduces, and the Poisson ratio of rock reduces, and the poroelasticity coefficient of rock increases;
And/or,
If described confined pressure characteristic parameter increases, then the hole fracture density of rock reduces, and the Young modulus of rock and Poisson ratio increase, and the poroelasticity coefficient of rock reduces;
And/or,
If described pore fluid characteristic parameter is fluid saturation, then along with the increase of fluid saturation, the Young modulus of rock, Poisson ratio and poroelasticity coefficient all increase; And, the Poisson ratio the when Poisson ratio in rock core oil-containing period is moisture greater than rock core, along with the increase of fluid saturation, its Poisson ratio difference increases;
And/or,
If described pore pressure characteristic parameter increases, then the Young modulus of rock, Poisson ratio reduce the increase of poroelasticity coefficient.
Step 105 is obtained the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy;
In embodiments of the present invention, the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy comprises the incidence relation of described external action parameter and strength of anisotropy.As previously mentioned, described external action parameter can comprise the temperature profile parameter, confined pressure characteristic parameter, pore fluid characteristic parameter, pore pressure characteristic parameter; In a kind of preferred exemplary of the embodiment of the invention, the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy can comprise:
The strength of anisotropy of described rock mechanics parameter becomes big along with the increase of temperature profile parameter;
And/or,
The strength of anisotropy of described rock mechanics parameter reduces along with the increase of confined pressure characteristic parameter;
And/or,
The strength of anisotropy of described rock mechanics parameter reduces along with the increase of pore fluid characteristic parameter; Wherein, described pore fluid characteristic parameter is fluid saturation;
And/or,
The strength of anisotropy of described rock mechanics parameter becomes big along with the increase of pore pressure characteristic parameter.
Rock mechanics parameter strength of anisotropy is influenced the aspect, and the inventor herein draws through experimental analysis, and the strength of anisotropy of rock mechanics parameter becomes big with temperature along with the increase of temperature.The strength of anisotropy of Biot coefficient is negative, and by the definition of strength of anisotropy as can be known, the Biot coefficient of vertical direction is greater than the Biot coefficient of horizontal direction, and this is because the hole fracture density of vertical direction causes greater than horizontal direction; The strength of anisotropy of the mechanical properties of rock of different lithology increases the variation tendency unanimity with temperature, but the amplitude difference.Mineralogical composition, structure, cementing matter and the hole fracture density of the temperature variant amplitude of various mechanical properties of rock and rock are contrasted, comparing result shows, different lithology strength of anisotropy amplitude of variation and hole fracture density have certain correlativity, but are not complete correspondences.Vertically go up the hole crack and grow than horizontal direction, rock is more loose, when heating, more is easy to generate new crack, and therefore, when temperature raise, the differences of mechanical properties between different directions can increase gradually.Anisotropy is incomplete corresponding with the hole fracture density, and it is a plurality of to show that anisotropic influence factor should have; According to the petrographic thin section observations, the difference that puts in order of the mineralogical composition of different directions and mineral is obvious, and from the influence of mineral to mechanical properties of rock, the composition of mineral also is a factor that influences rock anisotropy intensity with aligning.
According to rock anisotropy intensity with the change curve of confined pressure as can be known, the strength of anisotropy of rock mechanics parameter reduces with the increase of confined pressure.The strength of anisotropy of Biot coefficient is negative, by the definition of strength of anisotropy as can be known, the Biot coefficient of vertical direction is greater than the Biot coefficient of horizontal direction, according to the hole fissured structure of level and vertical direction as can be known, the vertical hole fracture density that goes up is greatly to cause the vertical big main cause of Biot coefficient that goes up.The strength of anisotropy of the mechanical properties of rock of different lithology increases the variation tendency unanimity with confined pressure, but the amplitude difference.Various mechanical properties of rock are contrasted with the amplitude of confined pressure variation and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows, different lithology strength of anisotropy amplitude of variation and hole fracture density have certain correlativity, but are not complete correspondences.Vertically go up the hole crack and grow than horizontal direction, rock is more loose, and when confined pressure increased, the crack was than the easier closure of horizontal direction, and therefore, when confined pressure increased, the differences of mechanical properties between different directions can reduce gradually.With anisotropy vary with temperature trend different be under the less situation of Young modulus amplitude of variation, bigger Poisson ratio and Biot amplitude of variation may be arranged with a kind of lithology.Anisotropy is incomplete corresponding with the hole fracture density, and the difference that shows different directions hole fracture density is not to cause anisotropic unique factor.According to the petrographic thin section observations, the difference that puts in order of the mineralogical composition of different directions and mineral is obvious, and from the influence of mineral to mechanical properties of rock, the composition of mineral also is a factor that influences rock anisotropy intensity with aligning.
According to rock mechanics parameter strength of anisotropy with the change curve of fluid saturation as can be known, the strength of anisotropy of rock mechanics parameter reduces with the increase of fluid saturation.The strength of anisotropy of Biot coefficient is negative, by the definition of strength of anisotropy as can be known, the Biot coefficient of vertical direction is greater than the Biot coefficient of horizontal direction, according to the hole fissured structure of level and vertical direction as can be known, the vertical hole fracture density that goes up is greatly to cause the vertical big main cause of Biot coefficient that goes up.The strength of anisotropy of the mechanical properties of rock of different lithology increases the variation tendency unanimity with fluid saturation, but the amplitude difference; Lithology of the same race, fluid type not simultaneously, its amplitude of variation is also inequality.Various mechanical properties of rock are contrasted with the amplitude of fluid saturation variation and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows, different lithology strength of anisotropy and hole fracture density have certain correlativity, but are not complete correspondences.Vertically going up the hole crack grows than horizontal direction, rock is more loose, when increasing identical fluid saturation, the fluid of vertical direction filling is more, the range of decrease of the amplification of its Young modulus Poisson ratio and Biot coefficient is greater than horizontal direction, the mechanics difference of different directions is dwindled, and strength of anisotropy diminishes, and it is identical with the variation tendency of confined pressure with mechanical properties of rock.Under the less situation of Young modulus amplitude of variation, bigger Poisson ratio and Biot amplitude of variation may be arranged with a kind of lithology.Anisotropy is incomplete corresponding with the hole fracture density, and the difference that shows different directions hole fracture density is not to cause anisotropic unique factor; According to the petrographic thin section observations, the difference that puts in order of the mineralogical composition of different directions and mineral is obvious, and from the influence of mineral to mechanical properties of rock, the composition of mineral also is a factor that influences rock anisotropy intensity with aligning.
According to the change curve of rock mechanics parameter strength of anisotropy with pore pressure, the strength of anisotropy of rock mechanics parameter increases with the increase of pore pressure.The strength of anisotropy of Biot coefficient is negative, by the definition of strength of anisotropy as can be known, the Biot coefficient of vertical direction is greater than the Biot coefficient of horizontal direction, according to the hole fissured structure of level and vertical direction as can be known, the vertical hole fracture density that goes up is greatly to cause the vertical big main cause of Biot coefficient that goes up.The strength of anisotropy of the mechanical properties of rock of different lithology increases the variation tendency unanimity with pore pressure, but the amplitude difference.Various mechanical properties of rock are contrasted with the amplitude of pore pressure variation and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows, different lithology strength of anisotropy amplitude of variation and hole fracture density have certain correlativity, but are not complete correspondences.Vertically go up the hole crack and grow than horizontal direction, rock is more loose, when pore pressure increases, the crack is than the easier unlatching of horizontal direction, therefore, and when confined pressure increases, differences of mechanical properties between different directions can increase gradually, and it is identical that its trend and anisotropy vary with temperature trend.Under the less situation of Young modulus amplitude of variation, bigger Poisson ratio and Biot amplitude of variation may be arranged with a kind of lithology.Anisotropy is incomplete corresponding with the hole fracture density, and the difference that shows different directions hole fracture density is not to cause anisotropic unique factor.According to the petrographic thin section observations, the difference that puts in order of the mineralogical composition of different directions and mineral is obvious, and from the influence of mineral to mechanical properties of rock, the composition of mineral also is a factor that influences rock anisotropy intensity with aligning.
Anisotropy is the important physical parameter of rock.The microfissure that aligns under lattice advantage orientation, the low confined pressure and the pore texture of rock are the topmost factors that determines the rock macroscopic anisotropy.Along with the increase of temperature, the orientation of mineral does not change, and vertically goes up fracture pore and pushes the speed and be higher than horizontal direction, and rock mechanics parameter strength of anisotropy strengthens; Along with the increase of confined pressure, the orientation of mineral does not change, but vertical speed that goes up microfracture close down is higher than the speed of crack closure on the horizontal direction, so the strength of anisotropy of rock mechanics parameter reduces gradually; When fluid saturation increased, vertically the Fluid Volume of filling was more in the rock sample, and the hole fracture density that vertical rock sample and horizontal rock sample are not filled with is tending towards approaching, and therefore, the strength of anisotropy of rock mechanics parameter reduces; When pore pressure increased, the clean confined pressure that acts on the rock skeleton reduced, and the fracture pore of vertical direction is all opened with the fracture pore vertically, but the vertical direction opening speed is higher than horizontal direction, and the anisotropy of rock mechanics parameter increases gradually.
Step 106, adopt described rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and, the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
With reference to figure 2, show the flow chart of steps of the method embodiment 2 that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics of the application recovers, specifically can comprise the steps:
Step 201 is measured the rock mechanics parameter of hydrocarbon-bearing pool;
In specific implementation, described step 201 can comprise following substep:
Substep S11 adopts static measurement method to obtain static rock mechanics parameter;
Substep S12 adopts dynamic measurement method to obtain dynamic rock mechanics parameter.
In embodiments of the present invention, described static rock mechanics parameter and dynamic rock mechanics parameter obtain through synchro measure, described static rock mechanics parameter can comprise static rock elasticity parameter and static rock strength parameter, and described dynamic rock mechanics parameter can comprise dynamic rock elasticity parameter and dynamic rock strength parameter; Described static rock elasticity parameter can comprise static Poisson ratio, static Young modulus and static poroelasticity coefficient etc., and described dynamic rock elasticity parameter can comprise dynamic Poisson ratio, kinetic Youngs modulus and dynamic hole elasticity coefficient etc.
In specific implementation, can adopt static measurement method to obtain static rock mechanics parameter by following substep:
Substep S11 carries out static state to rock sample and loads;
Substep S12 measures the stress-strain curve of described rock sample;
Substep S13, the stress-strain curve of the described rock sample of foundation calculates static rock mechanics parameter;
In specific implementation, can adopt dynamic measurement method to obtain dynamic rock mechanics parameter by following substep:
Substep S21 measures the velocity of propagation of sound wave in rock sample;
Substep S22 calculates dynamic rock mechanics parameter according to described velocity of propagation conversion.
Certainly, the measuring method of above-mentioned static rock mechanics parameter and dynamic rock mechanics parameter is only as example, and it all is feasible that those skilled in the art adopt any measuring method according to actual conditions, and the present invention need not this to be limited.
Step 202 adopts described static rock mechanics parameter and dynamic rock mechanics parameter, the sound diversity ratio of computing rock mechanics parameter;
Step 203, the rock mass mechanics property effect parameter of adding up described hydrocarbon-bearing pool;
Step 204 is calculated the strength of anisotropy of described rock mechanics parameter;
Step 205 is obtained the incidence relation of rock mechanics parameter and the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool;
Step 206, adopt described static rock mechanics parameter, dynamic rock mechanics parameter, the sound diversity ratio of rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and, the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
The embodiment of the invention is amplified rock mass mechanics concept, theory and method in the near surface engineering geology in the hydrocarbon-bearing pool rock mass of the complicated confined pressure in underground thousands of meters depths, high temperature, high pore pressure and heterogeneous fluid, form the concept of hydrocarbon-bearing pool rock mass mechanics, the various mechanical properties of thinking it to be underground reservoir rock mass porous media material showing under the specific degree of depth, temperature, three confined pressures, pore pressure and heterogeneous fluid saturation conditionss.Main oozed external environment condition restrictions such as build-in attribute such as condition and the degree of depth, temperature, confined pressure, pore pressure, heterogeneous fluid by rock composition, structure, structure, the hole of structural plane and structure inside.
The hydrocarbon-bearing pool rock mass mechanics is that a research hydrocarbon-bearing pool rock mass be out of shape under various force fields and the subject of failure law.The hydrocarbon-bearing pool rock mass mechanics is different from traditional rock mass mechanics, its study to as if be in the higher confined pressure of underground deep layer (can reach 200MPa), higher temperature (can reach 200 ° of C), therefore the heterogeneous body complicated rock mass of higher pore pressure (can reach 100Mpa) and saturated heterogeneous fluid has different research characteristics.The hydrocarbon-bearing pool rock mass mechanics is the basis of carrying out reservoir exploration exploitation measures such as terrestrial stress test, earthquake fluid prediction, time shift earthquake and reservoir waterfrac treatment transformation, and the Changing Pattern of furtheing investigate hydrocarbon-bearing pool rock mass under the underground prototype condition is significant.
The rock mass mechanics elastic parameter mainly comprises Poisson ratio, Young modulus and Biot coefficient etc.Poisson ratio is the performance of exosyndrome material transversely deforming when tension or pressurized, and it is exactly the ratio of transverse strain and longitudinal strain.Young modulus is to describe the physical quantity of solid material opposing deformability, and its physical significance is object when being subjected to stress generation elastic deformation, the ratio of stress and strain.The Biot coefficient is the ratio of the bulk modulus of the bulk modulus of dry rock and mineral grain, has reflected the contribution of pore space to the rock bulk property, is one of most important parameter in the poroelasticity The Characteristics.Poisson ratio, Young modulus and Biot coefficient are in drilling well, completion, borehole wall stability analysis, sanding prediction, formation drillability prediction and drill bit is preferred, have in the sleeve pipe damage mechanism, land subsidence design important use is all arranged; Simultaneously, because pore fluids such as oil gas can make the Poisson ratio of rock and Young modulus value that bigger change takes place, therefore, Poisson ratio and Young modulus also are the important parameters of oil gas geophysics (particularly exploration geophysics) research.Aspect terrestrial stress calculating, the Young modulus in the rock mass mechanics parameter and Poisson ratio equivalence are to carry out the important parameter that terrestrial stress calculates, and stress influence is fairly obvious over the ground; The occurrence in the crack of reservoir fracturing transformation process and bearing of trend are also determined by the rock mass mechanics parameter.
Hydrocarbon-bearing pool rock mass mechanics character is subjected to the influence of the factors such as composition, pore texture, confined pressure, pore pressure, saturated fluid and temperature of rock mass.If can understand the rock mechanics parameter with the Changing Pattern of these influence factors, just the seismic event that can observe according to ground and indoor rock mechanics experiment are finally inversed by the micro-structure (hole, crackle etc.) of elastic characteristic (as elastic modulus, Poisson ratio), surroundings (temperature and confined pressure), rock of deep under ground hydrocarbon-bearing pool and fluid state (fluid saturation, pore pressure) etc., thereby provide important basic data for efficient exploration and the reasonable development of hydrocarbon-bearing pool.
The mechanical property of hydrocarbon-bearing pool rock mass is subjected to the influence of many-sided factor, sums up to get up roughly can be divided into two classes: the one, and the build-in attribute of rock is as the mineralogical composition of rock, structure, structure, factor of porosity etc.; The 2nd, the residing external environment condition of rock is as heterogeneous fluid saturated in temperature, confined pressure, pore pressure and the hole etc.
Mechanical properties of rock, velocity of wave characteristic are the internal factors such as composition, structure, porous dehiscence crack quantity and distribution thereof that are controlled by rock.But for the rock under certain mode of occurence, its residing stress, temperature conditions also will inevitably produce a very large impact its mechanical property; Conventionally test can not reflect this environment to the influence of mechanical properties of rock, the mechanics parameter of under the condition of rock sample stratum of living in confined pressure, temperature, pore pressure, heterogeneous fluid, measuring only, and representative rock that could be real is in underground Mechanical Characters of Composite Ground.Simultaneously, the dynamic and static mechanics parameter of rock and its stress recovery, thermal equilibrium environment all have relation, therefore during concerning between research Rocks at Dynamic or Static Loading mathematic(al) parameter, should under identical condition, carry out, could guarantee precision and accuracy that the Rocks at Dynamic or Static Loading mathematic(al) parameter contrasts like this.
Below the basic mechanical character of hydrocarbon-bearing pool rock mass is described:
The hydrocarbon-bearing pool rock mass is under certain depth, the geology entity of be under uniform temperature, the pressure, saturated some fluids being made up of a series of structural planes and structure.It is the object of petroleum engineering geological research, geophysical survey, drilling well, completion, oil-gas mining and hydrocarbon zone transformation at target all be that rock mass is preserved in hydrocarbon-bearing pool.It is more and more that the difficulty that petroleum exploration and development now faces is adopted reserves, and object becomes increasingly complex, and the contradiction that the problem that the maturing field development late stage faces and the exploitation of traditional extensive style expose is more and more.Therefore these problems all need the theory and technology progress of rock mass mechanics aspect to solve, and the mechanical property of hydrocarbon-bearing pool rock mass are carried out detail knowledge and scrutiny has very important significance.
Because hydrocarbon-bearing pool rock mass external environment condition of living in is different with face of land rock mass with inner structure, its mechanical property more is subjected to the influence of ground environment, and its environment of living in has following characteristics:
1) the hydrocarbon-bearing pool rock mass is subjected to the influence of overlying strata pressure and horizontal direction terrestrial stress, and stress is very complicated, so its distortion is also very complicated and changeable with destruction.
2) it is higher to be in the hydrocarbon-bearing pool rock mass environment temperature of deep under ground.
3) being distributed with a large amount of holes, crack, crack and dissolution pore and solution cavity in the hydrocarbon-bearing pool rock mass, is the strong porous medium of heterogeneous body, differs widely such as metal material with common material.
4) not only be distributed with all kinds of spaces such as a large amount of holes, crack, solution cavity in the hydrocarbon-bearing pool rock mass, and be full of oil, gas, water multiphase fluid therein.Rock mechanics properties is subjected to the double influence of rock skeleton and fluid after the saturated fluid.
According to the stress-strain-time relation of hydrocarbon-bearing pool rock mass, its mechanical attribute can be divided into elasticity, plasticity and viscosity.Stress-strain-time relation curve map below in conjunction with rock shown in Figure 3 further specifies elasticity, plasticity and viscosity.
Elasticity refers to that in certain range of stress object is subjected to the external force effect to produce distortion, and can recover the character of its original shape and size size after the removal external force (load) immediately.The distortion of its generation is called elastic deformation, generally the material with elastic property is called elastic medium.Elasticity can be divided into two types again by its strain-stress relation: ess-strain linear elasticity in linear relation (or Hooke type elasticity, ideal elasticity) (referring to (a) among Fig. 3), and ess-strain is non-directional nonlinear elasticity.
After plasticity refers to object receiving force, when stress surpasses yield stress, still can continue distortion and not fracture at once, remove external force (load) after, the character that distortion can not recover again fully.Irreclaimable that part of distortion is called plastic yield, or claims permanent strain, residual deformation.The essence of object plastic yield is the slip that takes place under the shear stress effect between the interior of articles lattice.Plastic yield only taking place under external force, or the object of plastic yield only takes place in certain range of stress, is called plastic medium.The desirable stress-strain curve of plastic material is shown in (b) among Fig. 3, when stress is lower than yield stress σ yThe time, material character is elasticity; Stress reaches after the yield stress, and distortion constantly increases and stress is constant, and stress-strain curve is horizontal linear.After yield point, stress-strain curve is the rising curve, illustrates that crystal grain slides into after the reposition, causes phase embedding between well, jam-packed and crystal grain to increase, as makes it to continue to slide, and corresponding increase stress.This phenomenon is called as strain hardening, as (c) among Fig. 3.
Viscosity refers to be out of shape and can not finish instantaneous behind the object receiving force, and strain rate is called viscosity with the character that the size of stress changes.Desirable cohesive material (as Newtonian fluid), it is zeroaxial straight line that its stress-strain speed is closed, as (d) among Fig. 3.Strain rate is called flow deformation with the distortion of STRESS VARIATION.
Distortion and destruction relation according to rock can also be divided into rock property fragility and ductility.After fragility refers to object receiving force, be out of shape the character of just breaking when very little.Ductility refers to that object can bear moderate finite deformation and do not lose the character of its bearing capacity.The ductility of material and fragility are to divide according to the gradient size on negative slope on the overall strain before its stress destruction and the total stress-strain curve.It is little to destroy preceding overall strain, and the steep person in negative slope is fragility, otherwise is ductility.Be that standard is divided with 5% generally on the engineering, overall strain is ductile material greater than 5%, otherwise is hard brittle material.He De (Heard, 1963) is boundary with 3% and 5%, and rock is divided three classes: overall strain is brittle rock less than 3%; Overall strain is half brittle rock 3%~5%; Overall strain is the ductility rock greater than 5%.By above standard, most of surface rock all is fragility or half fragility under low confined pressure condition.Certainly the ductility of rock is relative with fragility, can transform mutually under certain conditions, and under at high-temperature and high-pressure conditions, the brittle rock under the normal temperature and pressure can show very high ductility.
Rock is the aggregate of mineral, have complicated composition and structure, so its mechanical attribute also is complicated.Which kind of proterties rock presents under external force, depends on its composition and structure on the one hand; On the other hand also and its stress condition, closely related as size and combined situation, load mode and speed, the stress path etc. of load.For example, at normal temperatures and pressures, rock is neither desirable resilient material, simply plasticity and cohesive material, and often show elastic-plastic, bullet-glue-mould or sticking-character such as elasticity.In addition, rock is composed the condition of depositing, and is also very big to its affect trait as temperature, water, terrestrial stress etc.
The rock mechanics parameter generally refers to elastic parameter (as Young's modulus of elasticity, modulus of shearing, bulk modulus and Poisson ratio etc.) and the intensive parameter (as uniaxial compressive strength, shear of rock intensity, rock tensile strength, cohesion and angle of internal friction etc.) of rock.The method of asking for the rock mechanics parameter mainly contains two kinds: the one, in the laboratory rock sample is surveyed, and the 2nd, ask for the rock mechanics parameter with the geophysics well-log information.
(1) Poisson ratio:
The Poisson ratio of rock refers to that rock is subjected to normal stress to do the time spent, and laterally the ratio of shortening and longitudinal tensile strain claims the horizontal compression coefficient again, is an important parameter of expression rock mechanics.The rock Poisson ratio is one of important parameter of estimating deformation behaviour of rock, in drilling well, completion, borehole wall stability analysis and FRACTURING DESIGN etc. important use is arranged all.For common rock, Poisson ratio is between 0~0.5.The Poisson ratio of rock is more big, represents that its elasticity is more little, and plasticity is more big, the more easy fracture of rock or pressure break.In the deposition rock, the Poisson ratio of mud stone is bigger, illustrates that its plasticity is big and yielding.
The computing method of Poisson ratio are divided into static and dynamic two kinds, and its computing formula is as follows respectively:
μ = Δϵ 1 Δϵ 2
μ = Δt s 2 - 2 Δt p 2 2 ( Δt s 2 - Δt p 2 ) = V p 2 - 2 V s 2 2 ( V p 2 - V s 2 )
Wherein, μ is the Poisson ratio of rock sample; Δ ε 1Axial strain increment for rock sample; Δ ε 2Transverse strain increment for rock sample; Δ t sBe the shear wave slowness of rock sample, μ s/ m; Δ t pBe the compressional wave time difference of rock sample, μ s/ m; V sBe the shear wave velocity of rock sample, m/s; V pBe the velocity of longitudinal wave of rock sample, m/s.
(2) Young modulus
Young modulus is rock in regime of elastic deformation the time, and rock has reflected that at the suffered stress of normal orientation and along the ratio of the elongation of answering forced direction to cause rock is subjected to external force to make the non-deformability of time spent, is indispensable physical quantity in the Mechanics Calculation; The Young modulus of rock is the basic mechanical character of rock, is research rock mass deformation, the requisite important parameter of a series of problems such as blast wave and seismic wave propagation in the rock.The Young modulus measuring method of rock has two kinds in experimental rock mechanics.A kind of is the load deflection experiment of rock, measures the relation between strain and the stress in the rock sample loading procedure, obtains the Young modulus of rock from the stress-strain curve of rock.Another kind method is to measure the velocity of propagation of elastic wave in rock, thereby calculates the Young modulus of rock again by the density of rock.The Young modulus that is obtained by first method is called static Young modulus, and the Young modulus that is obtained by second method is called kinetic Youngs modulus.The computing formula of static Young modulus and kinetic Youngs modulus is as follows respectively:
E = F / S ΔL / L
E = ρ V s 2 3 ( V p V s ) 2 - 4 ( V p V s ) 2 - 1
Wherein, E is the Young modulus of rock; Described ρ is the density of rock, g/cm 3F is the suffered power of test rock sample; S is the sectional area of rock sample; Δ L is that rock sample is being subjected to the elongation of force direction (or shortening amount); L is that rock sample is in the former length that is subjected on the force direction; V SBe the shear wave velocity of rock sample, m/s; Vp is the velocity of longitudinal wave of rock sample, m/s.
(3) Biot coefficient
The Biot coefficient is the ratio of the bulk modulus of rock volume modulus and the mineral grain of forming rock, and the Biot coefficient has reflected the contribution of pore space to the rock bulk property, and it is one of most important parameter in the poroelasticity feature.
The formation rock that contacts in the reservoir exploration exploitation all has hole, and being full of the fluid that contains pressure in the hole, pore pressure is the basis of carrying out reservoir engineerings such as oil gas well production forecast, terrestrial stress calculating, waterfrac treatment design, reservoir compacting, surface subsidence problem and Environmental Studies to the influence of rock elasticity strain.The Biot coefficient that has only definite rock could accurately calculate hole to the influence of the whole mechanical property of rock.
B = 1 - K dry K s
Wherein, B is the poroelasticity coefficient, abbreviates the Biot coefficient as, and Kdry is the bulk modulus (coefficient of volume compressibility) of dry rock, and Ks is the bulk modulus (the particle compressibility coefficient of solid) of the mineral of composition rock.
With reference to Biot coefficient test experiments synoptic diagram shown in Figure 4, adding confined pressure, hole pressure, temperature to experiment condition, add axial load to the confined pressure size, fixed orifice is pressed piston, the volume V of record rock 1With pore fluid flow F 1, confined pressure and axial load are increased 10MPa respectively, record rock volume V at this moment 2With pore fluid flow F 2, then can get:
ΔV f=F 2-F 1
ΔV s=V 2-V 1
The difference of rock deformation amount and hole deflection is rock skeleton particle deformation amount, then
ΔV m=ΔV s-ΔV f
By the definition of Biot coefficient as can be known:
B = 1 - K dry K s = 1 - 1 ΔV s 1 ΔV s - ΔV f = 1 - ΔV s - ΔV f ΔV s = ΔV f ΔV s
Wherein, B is the poroelasticity coefficient, abbreviates the Biot coefficient as, and Kdry is the bulk modulus (coefficient of volume compressibility) of dry rock, and Ks is the bulk modulus (the particle compressibility coefficient of solid) of the mineral of composition rock, Δ V fBe pore fluid fluctuations in discharge amount, Δ V sBe the rock skeleton volume change, Δ V mBe rock forming mineral particle volume variable quantity.
The rock mechanics parameter is the basic data that solves many petroleum engineering technology, and stable to solving drilling shaft lining, hydraulic fracture crack initiation and propagation law have significance.In specific implementation, the intensity of rock and static elasticity parameter can directly be measured by rock core test.But the rock core of obtaining in the laboratory might not represent the character of the rock stratum that will study fully, and these character core and process in may be changed.This just requires when carrying out shop experiment, needs to pay attention to taking and making of test specimen, makes it to have enough representativenesses, and keeps the natural structure state, is not subjected to unnecessary influence as far as possible; Making will be satisfied accuracy requirement to test specimen; Pay attention to the structure description of test specimen, as joint fissure development degree, distribution situation and orientation thereof etc.; Be also noted that sample dimensions (shape, size and ratio of height to diameter etc.), test condition and environment are to the influence of test.Though by taking suitable sample preparation method can solve subproblem, but still some problem needs to be resolved hurrily.This is because rock mechanics properties depends on the interaction of forming between crystal, particle and the consolidating material, and such as crack, joint, aspect with than the existence of minor fault.Be difficult on the one hand according to its character of composition particle this rock mechanics properties is described, particularly its intensity; On the other hand, the distribution of crack, joint, aspect and tomography is so changeable, to such an extent as to be subjected to the mechanical property of the bulk rock mass that this separation influences.
The experiment work of rock sample is the physics of understanding rock under varying environment, the important channel of mechanical property, also is to carry out the previous work that rock engineering should carry out.The purpose of shop experiment is to understand the character of rock, except carrying out physical property tests such as density, water percentage, porosity, water-intake rate and perviousness, should carry out the strength characteristics of rock and the test of deformation characteristic emphatically.Correct analysis theories, true and reliable experimental result can help to understand existence condition and the evolutionary process of the crust deep part that also can't observe directly so far.The rock that to fetch from the field or down-hole rock core drill through the rock core of normal diameter with the rotary rock drill model machine, cut suitable length with cutting machine again, with two-sided stone mill the rock section are polished at last, satisfy the accuracy requirement of sample.
Core sampling synoptic diagram with reference to shown in Figure 5 in taking from the on-the-spot general out-of-shape of core, can not be directly used in experiment.Need before the experiment on-the-spot core is processed.The process of indoor processing core is: extract earlier the cylindrical sample of a φ 25mm at the scene on the core with diamond core bit, and flat, the polishing of the two ends car with cylindrical sample then, the basal plane deviation makes length-diameter ratio 〉=1.5 of rock sample in 2.5% scope.
The measuring method of rock mechanics parameter is divided into two kinds of static test and dynamic tests.The assay method commonly used of rock mechanics parameter has two kinds: the one, and static measurement method calculates after rock sample being carried out measuring its stress-strain curve after static state loads in the laboratory; The 2nd, dynamic measurement method namely utilizes geophysical method, and by measuring the velocity of propagation of sound wave in rock sample, conversion calculates dynamic rock mechanics parameter again.
Static measurement method and dynamic measurement method respectively have its advantage and deficiency.On the one hand, dynamic measurement results is compared with the static measurement result, all has clear superiority at the aspects such as real-time of the quantity of data, procurement cost, data; What use in actual engineering on the other hand is the static elastic parameter of rock, and simultaneously existing mechanics constitutive relation is based on all generally that static parameter sets up, and therefore The result of dynamic test need be converted into static result and just can use.Therefore the influence factor of difference is carried out qualitative and quantitative description between rock mass mechanics parameter sound being measured, and can help us to utilize this instrument simply and easily of ultrasound wave to obtain the static parameter of rock.The research of the influence factor of sound difference has important practical significance.
Below laboratory dynamic mechanics parameter test is described:
The rock velocity of wave is the basic parameter that various fields such as oil-gas exploration and exploitation, petroleum drilling are paid close attention to.The velocity of propagation of sound wave in rock depends on the elastic constant of rocks such as the density, Young modulus, Poisson ratio, modulus of shearing, bulk modulus of rock, in the laboratory, the longitudinal and transverse wave velocity that utilization records just can be obtained the elastic modulus of rock, as Poisson ratio, compressibility coefficient, side cut modulus and Young modulus etc., the mechanics parameter of rock is the basis of petroleum drilling, fracturing engineering design and FRACTURE PREDICTION.
In the experiment of general acoustic wave test of rocks institute's externally applied forces normally small, instantaneous, therefore what think that rock shows is elastic characteristic, thereby the propagation of sound wave in rock also be to meet the propagation law of elastic wave, the fluctuation form of utilizing elastic wave wave equation and elastic wave calculate and the inverting rock mass in wave structure feature and distribution situation be feasible.During the walking of sound wave, the internal state of amplitude, frequency spectrum and rock has close contact, namely along with the degree of crushing increase of rock mass (rock), springy texture, stress reduce, sound wave can corresponding generation velocity of wave reduces, amplitude reduces, waveform degenerates, the dominant frequency in the frequency spectrum moves to low side.
Relation between the propagation law of sound wave and the rock mass mechanics index mainly shows:
(1) elastic modulus of rock mass and sound wave velocity of propagation therein has close mathematical relation, and the stretching of rock mass and compression deformation mainly are to rely on the compressional wave performance, can objectively reflect the physical mechanics characteristics such as intensity, distortion of rock mass;
(2) because mechanism's unevenness of rock mass makes sound wave produce very complicated variation in communication process, wave process is also ever-changing.Sound wave is propagated in rock mass and is caused the fissure-plane effect, can produce serious sound wave diffraction phenomenon what is more; Phenomenons such as the structural plane in the rock mass can make the velocity of wave of sound wave slack-off, and waveform generation is unusual.The factor that influences rock mass engineering project characteristic essence mainly is the variation that wherein has structural plane;
(3) velocity of propagation and its confined pressure of sound wave in rock mass has certain relation, this is by the resulting result of a large amount of shop experiments, increase along with the suffered pressure of rock, the velocity of propagation of sound wave also can increase thereupon, but rock mass reaches capacity pressure and after destroying, its velocity of wave also can correspondingly reduce;
(4) along with the increase of the compressive strength of rock mass, the velocity of propagation of sound wave also can increase, and good correlativity is arranged between the two.
In practice, the method for test rock and mineral elastic wave velocity mainly contains two kinds of sound wave pulse transmission beam method and reflectometrys.
(1) pulse transmission method:
The ultimate principle of pulse transmission method measurement elastic wave velocity is that length and the ultrasound wave of experiments of measuring sample passes the used time of laboratory sample, obtains the ultrasound wave elastic wave velocity of sample by experiment with laboratory sample length divided by the time.
When deducting the walking of top anvil up and down when during ultrasound wave walking in sample being total walking when base (system), and before experiment, can measure during the walking of top anvil; The initial length that the length of sample equals sample deducts compressed length, and compressed length is measured by HEL-electric pole displacement meter, and its linear measure longimetry precision is 0.01mm, and the velocity of wave of rock can be expressed as:
V=L/t。
(2) pulse reflection method:
The ultimate principle of ultrasonic pulse reflection method for measuring elastic wave velocity is by ultrasonic probe transponder pulse ripple, when running into two borders of barrier, pulsating wave reflects, reflection wave is received by same transducer, by the mistiming of twice reflection wave of oscillograph recording, obtain the ultrasound wave elastic wave velocity of sample by experiment with laboratory sample length divided by the time:
V=2L/t
Wherein, t is the mistiming that ultrasound wave passes through the rock upper and lower interface.
Because natural rock is made up of several mineral materials, often contains microfissure, acoustic attenuation is bigger.Pulse transmission method therefore commonly used is measured the elastic wave velocity of rock.
Below test describes to the dynamic modulus of elasticity, Poisson ratio:
The experimental principle of the described dynamic modulus of elasticity, Poisson ratio test can adopt following formula to characterize:
μ s = Δt s 2 - 2 Δt p 2 2 ( Δt s 2 - Δt p 2 ) = v p 2 - 2 v s 2 2 ( v p 2 - v s 2 )
E = ρ v s 2 3 ( v p v s ) 2 - 4 ( v p v s ) 2 - 1
Wherein, E is the Young modulus of rock; ρ is the density of rock, g/cm 3F is the suffered power of test rock sample; S is the sectional area of rock sample; Δ L is that rock sample is being subjected to the elongation of force direction (or shortening amount); L is that rock sample is in the former length that is subjected on the force direction; V SBe the shear wave velocity of rock sample, m/s; Vp is the velocity of longitudinal wave of rock sample, m/s; μ s is the Poisson ratio of rock sample; Δ ε 1Axial strain increment for rock sample; Δ ε 2Transverse strain increment for rock sample; Δ t SBe the shear wave slowness of rock sample, μ S/m; Δ t PCompressional wave time difference for rock sample.The longitudinal wave propagation speed definition in solid, liquid, gas owing to rising-elastic wave propagation speed that compression deformation produces; The shear wave propagation speed definition is in fixing, the elasticity velocity of wave propagation that produces owing to shearing deformation;
As from the foregoing, only need are measured rock density, velocity of wave can calculate the dynamic modulus of elasticity, Poisson ratio in length and breadth.
Below laboratory static mechanical parameter test is described:
Static rock test result is the basis of multiple rock test.At present the most static rock test method of long usefulness has rock uniaxial compression test, three compression tests of rock etc.
One, rock uniaxial compression test
Uniaxial compression test also claims no side direction confined pressure compression test, is the most frequently used rock strength test.The sample of test is generally right cylinder, and in order to reduce the influence of end effect, the ratio of length and diameter is generally 2~3.The pressing machine of test usefulness is wanted can load continuously and not impact, can test in the scope of gross tonnage (GT) 10%~90%.Measure sample deformation, the measurement instrument that all can meet the demands with precision and span.In the indoor compression test of standard, rock core normally through processing and place the crosshead of testing machine and worktable between carry out compression test.The suffered confined pressure of test specimen is zero, and the loading velocity pressurization of 5~8 kilograms every square centimeter of per second is adopted in test, till destroying.Axially loading is to load continuously, by measuring axial stress and axial and radial deformation, the research rock mechanics properties.
Two, three compression tests of rock
The rock of deep layer is in the anisotropic stress field, namely is subjected to the triaxial stress effect.In petroleum drilling or recovery process, near the shear stress values well or the production perforation increases.If the borehole wall is permeable, then having a principal effective stress component in the principal effective stress on well or pore wall is zero.And when mud cake obstruction fluid made it to pass through well or pore wall, this value can be higher slightly.Uniaxial compressive strength σ C is the important intensive parameter of describing wellbore stability.For the better stability of assessment well or pore texture, how must understand the mechanical property of rock and strength characteristics and be along with external load variation and variation.
According to the test of above-mentioned dynamic mechanics parameter and static mechanical parameter test can obtain sound difference to surround and watch mechanism as follows:
Dynamic and static load strain amplitude is different with the load frequency, and static state belongs to unlimited low-frequency big strain (10 -5~10 -3) load; Dynamic test is high-frequency small strain, and sound wave is small strain load, its typical frequency range: about tens hertz of seismic prospecting, acoustic logging 10~20KHz, indoor ultrasonic investigation 0.5~1.0MHz.The representative value of the strain that earthquake produces is 10 -6Just because of this difference of dynamic and static load, just make rock to the response difference of dynamic and static load, its dynamic and static elastic parameter is also different.
The dynamic and static elastic parameter of average materials such as iron and steel, aluminium, organic glass is basic identical.Why the dynamic and static elastic parameter of rock exists difference, is because the inherent feature of rock material itself causes.Rock is heterogeneous complex media, the microfissure that distributing in the rock, and there is fluid in microfissure inside, and this is the immanent cause of the dynamic and static elastic parameter difference of rock just.
To this, the inventor herein has carried out sound and has measured the difference experimental design.By original experimental provision is reequiped, set up the dynamic and static elastic parameter synchronism detection of cover rock device.The test rock sample is contained in the autoclave, and with the sealing of oil resistant rubber cover, applies oil pressure with pressing machine, provides confined pressure with wobble pump, and the stabilizing hydraulic pressure source provides stable pore pressure to rock sample.Measure the velocity of propagation of longitudinal and transverse ripple in rock sample with the ultrasonic transducer of the axial both ends of the surface of rock sample, obtain kinetic Youngs modulus and Poisson ratio after the conversion.The compressional wave resonance frequency of ultrasonic transducer is 850KHz, and shear wave is 450KHz, is coupled with masking foil between transducer and rock sample.The resistance strain gage on rock sample surface is used for measuring the distortion of rock sample loading procedure, by calculating static Young modulus and the Poisson ratio that may be able to obtain rock.Use HP3054A data acquisition and control system, whole mensuration process is finished under real-time monitoring of computer fully.The factor that influences the rock mechanics characteristic is more, and experimental result certain deviation can occur usually.The experimental result at fewer strong point may be caused bigger fluctuation for experimental analysis.Therefore, should adopt great deal of experiment data to carry out the more representational result that regretional analysis is lived.
Testing used rock sample is that the sandstone rock core that prolongs the oil field drills through and processes at lathe indoor with diamond head, its diameter is 53mm, its both ends of the surface polish, polish at grinding machine and sand paper, make the smooth finish of smooth finish, both ends of the surface of size, outside surface of rock sample and the depth of parallelism all reach the world and cover up society of theoretical and applied mechanics (ISRM) to rock static state and dynamic test request.
Before experiment, measure the length L (mm) of test core, diameter D 2Quality W(g) and longitudinal and transverse ripple T time delay of sound wave test macro (mm), Po, T So(μ s), and with the organic glass sample as standard component, dynamic test system has been carried out once test.Velocity of longitudinal wave by test organic glass sample is 2760m/s, and shear wave velocity is 1390m/s, and the deviation of the experimental data that these results and rock acoustical testing provide has proved that the dynamic test data are reliable within 5%.Before rock sample is carried out official testing, first precompressed three cycles of stress, eliminated the influence that stress lags behind basically.
In order to represent the sound method of testing to the influence of rock mechanics parameter, introduce the concept of rock mechanics parameter sound diversity ratio.Rock mechanics parameter sound diversity ratio:
E c = E d - E s E s
μ c = μ d - μ s μ s
Wherein, E cBe Young modulus sound diversity ratio; E dBe kinetic Youngs modulus; E sBe static Young modulus; μ cBe Poisson ratio sound diversity ratio; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio.
Below the influence factor of rock mass mechanics parameter is described:
The mechanical property of hydrocarbon-bearing pool rock mass is subjected to the influence of many-sided factor, sums up to get up roughly can be divided into two classes: the one, and the build-in attribute of rock is as the mineralogical composition of rock, structure, structure, factor of porosity etc.; The 2nd, the residing external environment condition of rock is as heterogeneous fluid saturated in temperature, confined pressure, pore pressure and the hole etc.
One, rock mass mechanics character internal influence factor
The mechanical property of hydrocarbon-bearing pool rock mass is subjected to the influence of many-sided factor, sums up to get up roughly can be divided into two classes: the one, and the build-in attribute of rock is as the mineralogical composition of rock, structure, structure, factor of porosity etc.; The 2nd, the residing external environment condition of rock is as heterogeneous fluid saturated in temperature, confined pressure, pore pressure and the hole etc.By analyzing the rock mechanics parameter with the Changing Pattern of these influence factors, just the seismic event that can observe according to ground and indoor rock mechanics are tested the micro-structure (hole, crackle etc.) of the elastic characteristic (as elastic modulus, Poisson ratio) that is finally inversed by deep under ground hydrocarbon-bearing pool, surroundings (temperature and confined pressure), rock and fluid state (fluid saturation, pore pressure) etc.
(1) rock composition characteristic parameter:
Rock composition comprises rock particles and chink.Rock composition is to embody the body internal cause of mechanical properties of rock, is the main body of different materials.Rock is by mineral composition, and the mineralogical composition of rock produces directly influence to the physico-mechanical properties of rock.Proportion such as the brightness feldspar is bigger than grouan, because the essential mineral composition of brightness feldspar is pyroxene and hornblende, its proportion is bigger than constituting granitic quartz and orthoclase; Quartzitic compressive strength rate marble is much higher, is because the quartzy high cause of strength ratio kalzit.As seen, although the rock class is identical, structure is also identical with structure, if the mineralogical composition difference, the physico-mechanical properties of rock also can present tangible difference, but can not think simply, and the rock strength that contains the high strength mineral is just certain high.After rock is stressed, if wherein the higher mineral of intensity do not contact in rock mutually, the influence of low-intensity mineral in the middle of then the stress transmission will be subjected to, rock just differs and demonstrates high intensity surely, have only the mineral of working as to be evenly distributed, when the high strength mineral form firm skeleton in The structure of rock, just can play the effect that improves rock strength.
(1) rock particles:
The sedimentogeneous rock particle mainly comprises quartz, feldspar and landwaste etc.Quartzy modal in petroclastic rock, most important detrital component, silica sand is a kind of hard, wear-resisting silicate mineral, its essential mineral composition is SiO 2Its physical property and chemical property are all very stable.Compositional maturity and textural maturity be all than higher, quartzyly in diagenetic process various variations can take place, and comprises quartzy secondary enlargement, corrosion, explanation and dissolution phenomena and stress deformation-wavy extinction and non-plastic fracture.
Feldspar is the aluminum silicate rock forming mineral of calcic, sodium and potassium, is the general name of feldspar family mineral.Feldspar ratio in the earth's crust all can occur in pyrogenic rock, metamorphosed rock, sedimentogeneous rock up to 60%.Feldspar is the essential mineral composition of nearly all pyrogenic rock, and the hardness fluctuations of feldspar is in 6-6.5, and proportion fluctuates in 2-2.5, and property is crisp, and high compression strength is arranged, and acid is had stronger chemical stability.Mainly contain potassium feldspar, soda feldspar, lime feldspar, baryta fledspar, baryta feldspar, microcline, orthoclase, sanidine etc.They all have vitreous luster, and color is varied, have colourless, adularescent, yellow, pink, green, grey, black etc.The feldspar component mainly is potassium feldspar and acid plagioclase, basic plagioclase is rarely found, and the compositional maturity of arkosic arenite and textural maturity are all lower, but feldspar can produce a large amount of secondary poress at diagenetic stage, therefore stable influential to rock, yielding and destruction.Secondary enlargement, corrosion, explanation dissolving and plasticity brittle deformation also can take place at the diagenetic process andesine.
Landwaste is the fragment of parent rock rock, is the mineral aggregate that keeps the parent rock structure.The landwaste of silexite, middle persilicic rock distributes the widest, but sedimentogeneous rock, chemical rock also can form landwaste.Landwaste composition difference, compositional maturity has marked difference.Its stability is poor more than quartz, and wherein the volcanics landwaste decomposes easily; The schist that is rich in schistosity is that the phyllite landwaste is easily deformable through diagenetic compaction; Some basin intraclasts as soft landwaste such as carbonate sand formation cuttings, oolith and spares, can be bent in diagenetic process, flatten, break, even are crushed to form false assorted base.The carbonate sand formation cuttings are easy to dissolving under the weak acidic medium effect, form a large amount of secondary poress, so its mechanical property is dynamic change in diagenetic process.
The mineral composition of rock, grain graininess, sorting, arrangement, psephicity etc. have also determined intensity and the stress characteristics of rock mass simultaneously because rerum natura and oiliness to the hydrocarbon-bearing pool reservoir have a direct impact.For specific petroliferous basin, under the close condition of Diagn, pore texture, chink and cementing matter, cementation type, degree of consolidation and fluid properties, rock skeleton mineralogical composition has bigger influence to mechanical property.The inventor herein chooses the part rock core that prolongs oil field prolongation group, and mineralogical composition and the mechanical property thereof of different rock cores have been carried out Experiments of Machanics and the analysis and research of XRD total rock, has analyzed the relation between mineral content and the mechanical properties of rock.XRD is the abbreviation of X-Ray Diffraction, is by rock sample being carried out X-ray diffraction, analyze its diffracting spectrum, the research means of information such as the structure of composition, inner atom or the molecule of acquisition rock sample or form.
Choose 4 and prolong oil field prolongation group with four blocks of rock cores of one deck position, four blocks of rock cores are argillaceous agglutination, all from the horizontal direction sampling, guarantee that its hole bedding architecture is similar.Rock is become to be grouped into by multi mineral, and the content of every kind of mineral all can influence rock mechanics properties to a certain extent.The mineral that this XRD identifies have 20 kinds nearly, only feldspar just comprises four types (lime feldspar, soda feldspar, potassium feldspar and microclines), for the ease of analyzing, adopted a short-cut method, the mineral of similar performance are sorted out, thereby mineral type is divided into quartz, feldspar, clay mineral and other mineral (content is less, mainly is kalzit, black mica etc.).
For example, a kind of example of testing rock sample rock mechanics parameter and mineral composition component list is as follows:
Figure BDA00003130394900341
Figure BDA00003130394900351
On the basis of tables of data, rock forming mineral is sorted out, and done the curve that crosses of mineral content and mechanical properties of rock, specifically can be with reference to mechanical properties of rock shown in Figure 6 and the mineral content curve synoptic diagram that crosses, as can be seen from Figure 6, along with the increase of quartz content, the Young modulus value of rock has the trend of rising; And along with the increase of feldspar content, Young modulus descends gradually.There is not tangible correlationship between clay mineral content and other mineral content and the Young modulus; The value of Poisson ratio mainly is subjected to clay mineral content, and data analysis shows that along with the increase of clay mineral content, Poisson ratio reduces gradually.The mineral content of Biot coefficient and rock relation not obvious (correlativity the highest be Biot coefficient and clay mineral content), the curve that crosses of Biot coefficient and mineral content shows that Biot coefficient and clay mineral content are proportionate basically, but local abnormity point is arranged also, it is a plurality of that analysis result shows that the factor that influences rock Biot coefficient may have, and causes itself and mineral content correlativity relatively poor.
(2) chink:
Chink mainly comprises clay mineral and cementing matter.Clay mineral is the important component of petroclastic rock.Clay mineral is the essential mineral of forming tonstein and soil.They are that some contain the hydrosilicate mineral that aluminium, magnesium etc. are the master.Except sepiolite, polygorskite tool chain layer structure, all the other equal tool layer structures.Particle is superfine, is generally less than 0.01 millimeter.Add the plasticity that has behind the water in various degree.Clay mineral can be divided into him by the origin cause of formation and give birth to (Lu Yuan) and two spontaneous classes, and the Lu Yuan clay mineral is the inheritance component that derives from parent rock; The authigenic clay mineral refer to the clay cement that forms in the chemical precipitation mode at deposition-diagenetic stage, and common have smalite, smectite, illite and chlorite.On structure, the authigenic clay mineral can be the film shape or the shell shape surrounds the petroclastic rock particle, wherein smectite and illite/smectite mixed layer clay and hydroexpansivity are strong, changeableness, poor stability, especially in drilling process, influence the stability of the borehole wall, can make permeability plugging in the expansion of waterflooding process medium clay soil.Clay depends on the position of clay particle in rock to the influence of seismic character, and the rock that has higher clay content on the statistics has lower seismic character and higher Vp/Vs.
Cementing matter refers to that the diagenesis phase is reaching the chemical precipitates that plays cohesive action outside detrital grain and the matrix between the rock particles.The composition of cementing matter is very complicated, main cementing matter is siliceous (quartz, calcedony etc.), carbonate mineral (kalzit, rauhkalk etc.), next is irony (haematite, limonite etc.), sometimes visible sulfate mineral (gypsum, muriacite etc.), zeolite minerals (euthalite, laumonite etc.), clay mineral (smalite, hydromica, chlorite etc.).Cementing matter often occupies pore space with the form of hole sark, and content generally is no more than 50% in petroclastic rock, and it plays cementation to detrital grain, makes it become hard rock.
Cementing matter is by the fixing rock particles of cementation.Cementation refers to be settled out mineral matter (cementing matter) from pore solution, with the loose fixed effect of sediment, it plays control action to the Elasticity feature of rock.The type of cementing matter directly influences degree of consolidation, thereby influences mechanical properties of rock, silico-calcium matter cementing matter consolidation strength height, and rock strength is big.The inventor herein chooses the part rock core that prolongs oil field prolongation group, and different cementing origin cause of formation rock mechanics properties has been carried out experimental study.Specifically can be with reference to rock sample Young modulus and quartz content histogram as shown in Figure 7, as shown in Figure 7, though the quartz content of calcitic cementation rock lower (9.25%), but because its cementation type is calcitic cementation, cementing matter intensity height, therefore the Young modulus of this rock sample is greater than other four kinds of rock samples (argillaceous agglutination, consolidation strength is low).
(2) rock texture characteristic parameter:
The The structure of rock feature is the key factor that influences physical-mechanical properties of rock.According to the The structure of rock feature, rock can be divided into two classes: the CRYSTALLINE ROCKS class, as most magmatite, metamorphosed rock and carbonate rock; Another kind of is the rock of cementing matter diagenesis, as petroclastic rock in the sedimentogeneous rock etc.
The crystallization diagenesis is by crystallization in magma or the solution, and recrystallization forms.Mineral crystal is bound up securely by the power of directly contact generation, and adhesion is strong.The rock porosity of crystallization diagenesis is little, compact structure, water-intake rate change for a short time, has higher intensity and stability than the rock of glued diagenesis.But with regard to the crystallization diagenesis, the size of crystallization crystal grain has a significant effect to the intensity of rock, and between 118~137MPa, what granitello had then can reach 196~245MPa as the compressive strength of coarse-grained granite; And for example the compressive strength of marble is generally between 79~118MPa, but the firmest ls then about 196MPa, even can reach 255MPa.The rock that mineralogical composition is identical with structure type is described, the size of crystallizable mineral crystal grain is fairly obvious to the influence of its intensity.
For example, a kind of example of rock mechanics parameter list of different structure rock is as follows:
Figure BDA00003130394900371
Figure BDA00003130394900381
As seen from the above table, the value of the Young modulus of the rock of the crystal class origin cause of formation and Poisson ratio generally is higher than the rock of cementing matter diagenesis.
(3) rock structure characteristic parameter:
In the rock between the different minerals aggregate, between each ingredient of rock or the mutual relationship between mineral aggregate and other ingredients of rock, be called rock structure.Also think that by the people structure of rock should be mutual relationship and the filling method in shape, size and the space of forming the mineral aggregate of rock, i.e. the geometric feature of the combination of these mineral aggregates.For example, sheet fiber crops structure, clumpy structure, fluxion structure, pillow structure, pore columnar structure, miarolitic structure etc.
The structure of rock mainly is to be determined by the space distribution of each ingredient of rock and mutual Rankine-Hugoniot relations thereof to the influence of its physico-mechanical properties.As having schistose structure, tabular structure, phyllitic structure, gneissose structure when rock, and during rhyotaxitic structure, its mineralogical composition distributed pole in rock is inhomogeneous.Some intensity are low, the mineral of easy-weathering, and how in a certain direction enrichment slivering zonal arrangement perhaps becomes local aggregation, and the physico-mechanical properties that makes rock in a certain direction or local great changes will take place.Rock stress destruction and rock suffer weathering, at first all are that these defectives from rock begin to take place.Another kind of situation is, though the distribution of different mineralogical composition in rock is uniformly, owing to exist the hole of bedding, crack and the various origin causes of formation, and make rock mechanics properties present evident difference in different directions.
The inventor herein uses the method for CT scan that bedding architecture and the hole FRACTURE CHARACTERISTICS of rock sample are studied.CT is the abbreviation of computer tomography technology (Computed Tomography), is to be proposed the earliest and designed by Britain Hounsfield in 1972.To the later stage eighties, the CT technology begins to be applied to hydrocarbon-bearing pool research, and develops into the important tool of research reservoir porous characteristic.Use the CT scan imaging technique, can under the prerequisite that does not change rock core formalness and inner structure, observe the internal void structure of rock core, core porosity is distributed carry out quantitatively characterizing, for the seepage flow mechanism of studying low permeability reservoir provides foundation.
The inventor herein chooses the vertical stratification direction and the horizontal stratification direction is taken a sample respectively, and the rerum natura mechanical property of the rock sample of two orthogonal directions is tested.The mechanics parameter table that obtains different directions sampling rock core is as follows:
Figure BDA00003130394900391
As can be known from the above table, the factor of porosity that horizontal direction is bored sample is greater than the factor of porosity of vertical direction rock sample, and the value of Young modulus and Poisson ratio is also greater than the value of vertical direction.
(4) blowhole FRACTURE CHARACTERISTICS parameter:
Have hole in the rock, hole can tear seam/crack (crack) and hole (pore) by its shape; Be divided into by the origin cause of formation and time sequencing: primary pore crack and secondary pores crack.The hole FRACTURE CHARACTERISTICS of rock is subjected to shouting the influence of factors such as feature, inner wall coarse degree, rock particles size, grain sorting and compacting.
Hole shared volume in rock is very little, the but influence to rock property very big (when existing crack and crack to be compacted, the bulk modulus of rock can differ 5 times).The method that the inventor herein utilizes constant speed to press mercury and rock mechanics test to combine, relation between rock mechanics parameter and the space FRACTURE CHARACTERISTICS is studied, the result shows that along with the increase of hole fracture spaces, the Young modulus of rock and the value of Poisson ratio reduce gradually.
Two, rock mass mechanics character external influence factors:
The external influence factors that the mechanical property of hydrocarbon-bearing pool rock mass is subjected to comprises heterogeneous fluid saturated in temperature, confined pressure, pore pressure and the hole etc.If can understand the rock mechanics parameter with the Changing Pattern of these influence factors, just the seismic event that can observe according to ground and indoor rock mechanics experiment are finally inversed by elastic characteristic (as elastic modulus, Poisson ratio), surroundings (temperature and confined pressure) and the fluid state (fluid saturation, pore pressure) etc. of deep under ground hydrocarbon-bearing pool, thereby provide important basic data for efficient exploration and the reasonable development of hydrocarbon-bearing pool.
The inventor herein adopts the method for single argument research, changes the value of unitary variant, simultaneously the value of other variable is remained on the scope of oil-gas stratigraphic deposits, by the Changing Pattern of record rock mechanics parameter with unitary variant; Simultaneously experimental apparatus is reequiped, when guaranteeing the rock mechanics experiment, the sound test is carried out synchronously; To a kind of rock sample, carried out the sampling respectively of vertical direction and horizontal direction in addition, the anisotropy of rock mechanics parameter has been analyzed.On this basis, this chapter to influence the rock mechanics parameter with and sound difference, anisotropic external factor carried out quantitative test, and provided the relational expression of static and dynamic correction.
(1) temperature profile parameter:
Hydrocarbon-bearing pool is under the hot environment of depths, stratum, owing to can't directly observe the material in earth deep, mainly rely on the structure in geophysical method inverting earth deep at present, and the not uniqueness of inversion result has caused very big inconvenience to explaining underground actual configuration.Therefore, by under condition of different temperatures to the rock mechanics properties research that experimentizes, obtain the rock mechanics parameter with the variation of temperature rule, reliable constraint condition and theoretical foundation can be provided for the reasonable dismissal of geophysics inversion result.
Simultaneously; in the waterflooding of conventional oil reservoir and the gas injection development process of heavy crude reservoir; because the temperature of injected media and the temperature contrast of oil reservoir are very big; after exploitation for many years; change has taken place in the temperature of oil reservoir; this has also caused the change of oil reservoir rock mass mechanics character, and the research temperature is conducive to casing protection and the remaining oil research technology of latter period of oilfield development to the influence of rock mass mechanics character.
To this, the inventor herein chooses 5 groups of rock samples and experimentizes, and every group of rock sample comprises two parallel sample, and one is the horizontal direction sampling, and one is the vertical direction sampling.In the test, sample is carried out sound measure simultaneously, the Young modulus of rock, Poisson ratio and Biot coefficient are to keep confined pressure, pore pressure and fluid saturation constant, and measurement calculates under the change temperature conditions; Confined pressure, pore pressure and oil-containing water saturation are chosen the oil reservoir real data, are respectively 20Mpa, 10Mpa, 50%, 50%, and range of temperature is 30-165 ° of C, measure a rock mechanics parameter every 15 ° of C.
Obtained following experimental result: along with the increase of temperature, the value of Young modulus and Poisson ratio reduces gradually; Along with the increase of temperature, the Biot coefficient raises gradually; The mechanical properties of rock of different lithology increases the variation tendency unanimity with temperature, but the amplitude difference.Mineralogical composition, structure, cementing matter and the hole fracture density of the temperature variant amplitude of various mechanical properties of rock and rock are contrasted, comparing result shows, it is littler than the rock of argillaceous agglutination that the Young modulus of the rock of calcitic cementation varies with temperature amplitude, but its Poisson ratio and Biot index variation amplitude and argillaceous agglutination rock type are seemingly; Various lithological change amplitude trend and its hole fracture density good relationship, but neither be corresponding fully, showing that the change of mechanical properties of rock was subjected to the influence of many-sided factor when temperature raise, the hole fracture density may be major influence factors.On this analysis foundation, the theory relation between mechanical properties of rock and the hole fracture density and blowhole fracture density have been carried out observation analysis with the variation that the temperature increase takes place.
A. temperature is to the influence of rock fracture void density
Rock sample is packed into behind the triaxial cell, apply axial compression and confined pressure to setting value.Heat gradually then, install 4 calibrate AE sensors in the upper and lower, triaxial cell respectively.The acoustic emission registering instrument detects the rock burst acoustic emission signal of each temperature section in real time.Acoustic emission signal comprises parameters such as acoustie emission event number, energy, duration, amplitude.300 acoustic emission signals of per second record, its data volume is huge, analyzes for convenient, and acoustic emission signal is average by second, draws temperature-acoustic emission time diagram again.
Fig. 8 is from 30-150 degree temperature section packsand count of acoustic emission event rate under the high temperature condition of triaxial stress, as shown in Figure 8: (1) the acoustics record occurs, but number of times is less at 35 ° of C-40 ° of C places, illustrates that other thermal rupture event has appearred in this section rock; The acoustics record occurs continuously near 60 ° of C, illustrates that this section rock burst is more active; (2) when the temperature range of 110 ° of C-135 ° of C, the acoustie emission event relatively placidity; (3) when the temperature range of 135 ° of C-150 ° of C, the acoustics record occurs continuously, illustrates that more violent thermal rupture has taken place this section rock.
Rock has experienced the violent phase of a plurality of thermal ruptures in temperature-rise period, its acoustie emission event number, duration, energy etc. present several peak regions.During heating-up temperature, the sample internal moisture is overflowed to the infiltration of rock mass outside along the crack in a large amount of evaporation processes, thereby makes the hole crack widen lengthening.Simultaneously, action of high temperature makes and the volumetric expansion of sample endoparticle also makes the sample internal crack increase.Along with the rising of temperature, the void density of rock increases gradually.
B. the rock fracture void density is to the influence of mechanical properties of rock
According to the theory of damage mechanics, from the angle of temperature to the influence of blowhole and fracture density, temperature is analyzed the influence of mechanical properties of rock.
(b1) temperature is to the influence of Young modulus
According to the theory of damage mechanics, can draw:
E , E = 1 + 16 45 ( 1 - μ 2 ) ( 10 - 3 μ ) 2 - μ f
In the formula, E ' is the effective modulus of elasticity that contains the rock particles of crackle; E is the Young modulus of rock; μ is the Poisson ratio of rock particles; F is micro-crack density.
Can be got by following formula, along with the increase of crack density, the Young modulus value of rock reduces; Otherwise along with reducing of crack density, the effective modulus of rock improves.Along with the rising of temperature, the hole fracture density of rock increases, and the Young modulus of rock reduces gradually.
(b2) temperature is to the influence of Poisson ratio
Contain the Poisson ratio μ ' of crackle rock and the crack density ε of rock following relation arranged:
ϵ = 16 45 μ ′ + 2 ( 1 - ξ ) ( 1 + 3 μ ) - 1 2 - 4 μ
In the formula, ε is the density of crackle, and μ is the Poisson ratio that does not contain the crystalline matrix of crackle, ξ be crackle contain the liquid saturation degree, μ ' is the Poisson ratio that contains the crackle rock.
What this experiment was adopted is the single argument research method, so in the humid test, the liquid saturation degree that contains of rock sample remains unchanged, so ξ is definite value, and by following formula as can be known, Poisson ratio and the crack density of rock are directly proportional.Along with the rising of temperature, the hole fracture density of rock increases, and the Poisson ratio of rock reduces gradually.
(b3) temperature is to the influence of Biot coefficient
Figure BDA00003130394900432
In the formula, B is the poroelasticity coefficient, abbreviates the Biot coefficient as, K DryBe the bulk modulus (coefficient of volume compressibility) of dry rock, Ks is the bulk modulus (the particle compressibility coefficient of solid) of the mineral of composition rock.
Along with reducing of crack density, the value of the bulk modulus of dry rock is more and more approaching, and then the value of Biot coefficient is more and more littler.Along with the rising of temperature, the hole fracture density of rock increases, and the Biot coefficient of rock increases gradually.
Comprehensively to the heat observation of phenomenon in the process of rock, the inventor herein thinks that breaking appears in rock in the process of heating and causes that it is the basic reason that causes mechanical properties of rock to change that the rock fracture volume of voids increases.Along with the increase of temperature, the sample internal moisture is sought the crack and is overflowed to the infiltration of rock mass outside in a large amount of evaporation processes, thereby makes hole widen lengthening.Simultaneously, action of high temperature makes and the volumetric expansion of sample endoparticle also makes the sample internal crack increase.Along with the rising of temperature, the void density of rock increases gradually.Along with the increase of temperature, the hole fracture density of rock increases, thereby causes the Young modulus of rock and the value of Poisson ratio to diminish, and the value of Biot coefficient increases.
According to rock sound difference with the variation of temperature curve as can be known, the sound diversity ratio of the Young modulus of rock and Poisson ratio increases gradually; The mechanical properties of rock of different lithology increases the variation tendency unanimity with temperature, but the amplitude difference.Mineralogical composition, structure, cementing matter and the hole fracture density of the temperature variant amplitude of various mechanical properties of rock and rock are contrasted, comparing result shows that the correlativity of different lithology sound difference amplitude and mineral content, cementing matter and hole fracture density, structure is all relatively poor.Microscopic mechanism according to sound difference (dynamically is high frequency-little displacement, static state is low frequency-big displacement), the hole crack should be the main cause that influences sound difference, but the result from this test, the sound diversity ratio increases (increase of hole fracture density) with temperature and increases, this point matches with the mechanism of surrounding and watching of sound difference, but the amplitude of variation of different lithology is bad with corresponding hole fissured structure match result, this explanation influences the reason of rock sound difference amplitude may be more complicated (in temperature changing process, character such as the velocity of wave of rock all can change etc.).
Satisfy binomial relation between rock mechanics parameter sound diversity ratio and the temperature, namely
E d - E s E s = aT 2 + bT + C
μ d - μ s μ s = mT 2 + nT + k
In the formula, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; T is temperature; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant.The regression formula of five kinds of rock samples of inventor herein's statistics is as shown in the table:
Wherein, R 2Be related coefficient.
(2) confined pressure characteristic parameter:
The hydrocarbon-bearing pool rock mass is under the higher confined pressure, measured Poisson ratio is not suitable for the situation of subterranean strata under atmospheric conditions, rock mechanics properties has greatly changed facing to the change of confined pressure, this phenomenon is found by many experimental rock mechanics workers, but along with increase or the reduction of confined pressure, the Changing Pattern of mechanical properties of rock is gone back neither one at present and is set forth more clearly.
The inventor herein experimentizes by choosing 5 groups of rock samples, and every group of rock sample comprises two parallel sample, and one is the horizontal direction sampling, and one is the vertical direction sampling.In the test, sample is carried out sound measure simultaneously, the Young modulus of rock, Poisson ratio and Biot coefficient are to keep temperature, pore pressure and oil-containing water saturation constant, and measurement calculates under the change confined pressure condition; Temperature, pore pressure and fluid saturation are chosen the oil reservoir real data, are respectively 40 ° of C, 10Mpa, 50%, 50%, and the confined pressure variation range is 20-90Mpa, and each 10Mpa measures a rock mechanics parameter.
Obtain following experimental result: along with the increase of confined pressure, the value of Young modulus and Poisson ratio increases gradually, and the Biot coefficient reduces gradually; The mechanical properties of rock of different lithology increases the variation tendency unanimity with temperature, but the amplitude difference.Various mechanical properties of rock are contrasted with the amplitude of confined pressure variation and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows that the Young modulus of the rock of calcitic cementation, Poisson ratio and Biot coefficient are littler than the rock of argillaceous agglutination with the confined pressure amplitude of variation.Various lithological change amplitude trend and its hole fracture density correlativity are good, and variation tendency is almost coincide; Relatively poor with the correlativity of mineralogical composition.Comparative analysis is the result show, when confined pressure raise, the change of mechanical properties of rock mainly was subjected to the influence of hole fracture density.On this analysis foundation, in confined pressure increase process, the variation of blowhole fracture density is observed.
Along with the continuous increase of confined pressure, the microfracture in the rock is constantly closed.Cause in the process of crack closure in the confined pressure increase, the uneven micro-crack that makes of surface of contact is closed brokenly.Cause in the process of crack closure in the confined pressure increase, the uneven micro-crack that makes of surface of contact is closed brokenly, and a micro-crack can become several, and the length of micro-crack reduces tempestuously, and its width slowly reduces.Crackle is closed gradually under low confined pressure effect, but the variation of its geometric configuration is little.After confined pressure surpassed certain value, the micro-crack change of shape was apparent in view, become gradually flat, until closure.Along with the increase of confined pressure, fracture width diminishes even is closed, and the fracture density of rock diminishes, and it is finer and close that rock becomes.
Along with reducing of fracture density, the Young modulus of rock, Poisson ratio increase gradually, and the Biot coefficient value reduces.Along with the increase of confined pressure, the hole crack of rock dwindles gradually or is closed, and therefore, for rock mechanics parameter influence aspect, along with the increase of confined pressure, the Young modulus of rock, the value of Poisson ratio increase, and the value of Biot coefficient diminishes.
Along with the increase of confined pressure, the sound diversity ratio of the Young modulus of rock and Poisson ratio reduces gradually, and the mechanical properties of rock of different lithology increases the variation tendency unanimity with confined pressure, but the amplitude difference.Mineralogical composition, structure, cementing matter and the hole fracture density of the temperature variant amplitude of various mechanical properties of rock and rock are contrasted, comparing result shows that the correlativity of different lithology sound difference amplitude and mineral content, cementing matter and hole fracture density, structure is all relatively poor; Under the less situation of Young modulus amplitude of variation, bigger Poisson ratio amplitude of variation may be arranged with a kind of lithology.Microscopic mechanism according to sound difference (dynamically is high frequency-little displacement, static state is low frequency-big displacement), the hole crack should be the main cause that influences sound difference, but the result from this test, the sound diversity ratio increases (the hole fracture density reduces) with confined pressure and reduces, this point matches with the mechanism of surrounding and watching of sound difference, but the amplitude of variation of different lithology is bad with corresponding hole fissured structure match result, and the reason that this explanation influences rock sound difference amplitude is comparatively complicated.Satisfy binomial relation between rock mechanics parameter sound diversity ratio and the confined pressure, that is:
E d - E s E s = aP c 2 + bP c + C
μ d - μ s μ s = mP c 2 + nP c + k
In the formula, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; P cBe confined pressure; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant.The regression formula of five kinds of rock samples of inventor herein's statistics is as shown in the table:
Figure BDA00003130394900481
Wherein, R 2Be related coefficient.
(3) pore fluid characteristic parameter:
The hydrocarbon-bearing pool rock mass is porous medium, but filling water, gas or wet goods fluid therebetween.The existence of fluid can influence the seismologic parameter characteristic of rock, and wave field characteristics such as speed, decay are changed along with the variation of rock fluid saturation.The extraction of oil and the injection of water also can cause mechanical properties of rock to change in the oilfield process simultaneously, cause occurring a large amount of sleeve pipes and damage.Therefore study the variation of blowhole fluid content to the influence of physical properties of rock, the understanding of reservoir characteristics in the oil-gas exploration and the analysis of seismic data are all had great importance.
The inventor herein experimentizes by choosing 5 groups of rock samples, and every group of rock sample comprises two parallel sample, and one is the horizontal direction sampling, and one is the vertical direction sampling.In the test, sample is carried out sound measure simultaneously, the Young modulus of rock, Poisson ratio and Biot coefficient are to keep temperature, pore pressure and oil-containing water saturation constant, and measurement calculates under the change confined pressure condition; Temperature, pore pressure and fluid saturation are chosen the oil reservoir real data, are respectively 40 ° of C, 10Mpa, 50%, 50%, and the confined pressure variation range is 20-90Mpa, and each 10Mpa measures a rock mechanics parameter.
Obtain following experimental result: along with the increase of fluid saturation (oil-containing or moisture), Young modulus, Poisson ratio and Biot coefficient all increase gradually, and the mechanical properties of rock of different lithology increases the variation tendency unanimity with temperature, but the amplitude difference; Simultaneously, different fluid type (oil-containing and moisture), its amplitude of variation is also variant.Mechanical properties of rock is all different with confined pressure, temperature and pore pressure with the variation tendency of fluid saturation.Various mechanical properties of rock are contrasted with the amplitude of variation of fluid and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows, the correlativity of various lithological change amplitude trend and these influence factors is all relatively poor, and what degree of agreement was the highest is the hole fissured structure.Analysis result shows, fluid saturation influences mechanism than more complicated other three kinds of factors to mechanical properties of rock.On this analysis foundation, in the convection cell saturation degree increase process, rock mechanics properties changes to be analyzed.
For the rock of drying, rock is in " not discharge opeing state ", and according to the Hadizadeh.J theory, compressibility coefficient β and the shear modulus G that then can get rock this moment are respectively:
β=β eff=β s-dη/dP
G=G eff
In the formula, β EffCompressibility coefficient for dry rock; G EffModulus of shearing for dry rock; β sBe rock matrix compressibility; η is rock porosity; P is blowhole hydraulic pressure.
For the rock of complete saturated fluid, rock is in " discharge opeing state ", and compressibility coefficient β and shear modulus G under the draining situation are respectively:
β = β eff + ( β eff - β s ) 2 η ( β P - β s ) 2 + ( 1 - η ) ( β eff - β s ) 2
G=G eff
For natural rock, the fluid saturation of rock may be between the 0-1, and just rock " discharge opeing " and " not discharge opeing " two kinds of situations may occur simultaneously when being under pressure.Therefore, for the influence of true description fluid to the rock deformation characteristic, introduce the average discharge opeing coefficient of blowhole, its implication is exactly the cumulative volume of rock discharge opeing hole and the ratio of blowhole cumulative volume, and its value equals the pick up of rock, namely contains the liquid saturation degree.Then the compressibility coefficient β of natural rock and shear modulus G are:
β = aβ eff + ( 1 - α ) · ( β eff + ( β eff - β s ) 2 η ( β P - β s ) 2 + ( 1 - η ) ( β eff - β s ) 2 )
G=G eff
Can be got by following formula:
d&beta; d&alpha; = - ( &beta; eff - &beta; s ) 2 &eta; ( &beta; P - &beta; s ) 2 + ( 1 - &eta; ) ( &beta; eff - &beta; s ) 2 ) < 0
So rock compressibility β reduces along with the increase that contains the liquid saturation degree.Close between the Young modulus of rock, Poisson ratio and compressibility coefficient, the modulus of shearing and be:
μ=(3-2βG)/(6+2βG)=9/(6+2βG)-1
E=9G/(3+Gβ)
Can be got by following formula, μ, E and β are inversely proportional to, and namely are directly proportional with containing liquid saturation degree α, and namely along with the increase that contains the liquid saturation degree, μ, E also increase thereupon.As from the foregoing, the existence of pore fluid has increased the effective bulk modulus of rock, because the modulus of shearing of fluid is 0, so after being full of fluid, the modulus of shearing of rock remains unchanged, so the Poisson ratio of rock and Young modulus value become big.
Contain the liquid saturation degree to the influence of Biot coefficient:
B = 1 - K dry K s = 1 - 1 &Delta;V s 1 &Delta;V s - &Delta;V f = 1 - &Delta;V s - &Delta;V f &Delta;V s = &Delta;V f &Delta;V s
&Delta;V s = &Delta;P K
K = 1 &beta;
Can be learnt by above three formulas
B = &Delta;V f &Delta;P &CenterDot; &beta;
By following formula as can be known, thus B and β value be inversely proportional to because β and fluid saturation are inversely proportional to, so B and contain the liquid saturation degree and be directly proportional.Along with the increase that contains the liquid saturation degree, the Biot coefficient of rock increases.
Analyze theoretically, the increase of fluid saturation directly influences rock mechanics properties, do not need to change mechanical properties of rock by other state-owned inherent attribute, the inventor herein thinks, after the saturated with fluid rock, rock has become a kind of " new material " that flows solid coupling by original solid phase, and fluid has become the part of " new material ", so its influence to mechanical properties of rock is more direct.
By rock sound difference with the variation of fluid saturation as can be known, along with the increase of fluid saturation, the sound diversity ratio of the Young modulus of rock and Poisson ratio reduces gradually; The mechanical properties of rock of different lithology increases the variation tendency unanimity with confined pressure, but the amplitude difference.Same lithology but fluid type not simultaneously, its amplitude of variation is also variant.Various rock mechanics sound differences are contrasted with the amplitude of fluid saturation variation and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows that the correlativity of different lithology sound difference amplitude and hole fracture density, mineralogical composition and cementing matter is all relatively poor; Under the less situation of Young modulus amplitude of variation, bigger Poisson ratio amplitude of variation may be arranged with a kind of lithology.Microscopic mechanism according to sound difference (dynamically is high frequency-little displacement, static state is low frequency-big displacement), the hole crack should be the main cause that influences sound difference, result from this test, the sound diversity ratio reduces with the increase (the hole fracture density reduces) of fluid saturation, and this point matches with the mechanism of surrounding and watching of sound difference; But the amplitude of variation of different lithology is relatively poor with corresponding hole fissured structure correlativity, the reason that this explanation fluid saturation influences rock sound difference amplitude is comparatively complicated, fill the crack except fluid and cause the dwindling of factor of porosity, fluid also can exert an influence to rock mechanics properties to the dissolving of clay mineral etc.Satisfy binomial relation between rock mechanics parameter sound diversity ratio and the fluid saturation, that is:
E d - E s E s = a 1 S w 2 + b 1 S w + C 1
&mu; d - &mu; s &mu; s = m 1 S w 2 + n 1 S w + k 1
E d - E s E s = a 2 S o 2 + b 2 S o + C 2
&mu; d - &mu; s &mu; s = m 2 S o 2 + n 2 S o + k 2
In the formula, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; S wBe water saturation; S oBe oil saturation; a 1, a 2Be the quadratic term coefficient; b 1, b 2Be coefficient once; C 1, C 2Be constant; m 1, m 2Be the quadratic term coefficient; n 1, n 2Be coefficient once; k 1, k 2Be constant.
The regression formula of five kinds of rock samples of inventor herein's statistics is seen shown in following two tables:
Figure BDA00003130394900525
Figure BDA00003130394900531
Wherein, R 2Be related coefficient.
By rock mechanics parameter strength of anisotropy with the variation of fluid saturation as can be known, the strength of anisotropy of rock mechanics parameter reduces with the increase of fluid saturation.The strength of anisotropy of Biot coefficient is negative, by the definition of strength of anisotropy as can be known, the Biot coefficient of vertical direction is greater than the Biot coefficient of horizontal direction, according to the hole fissured structure of level and vertical direction as can be known, the vertical hole fracture density that goes up is greatly to cause the vertical big main cause of Biot coefficient that goes up.The strength of anisotropy of the mechanical properties of rock of different lithology increases the variation tendency unanimity with fluid saturation, but the amplitude difference; Lithology of the same race, fluid type not simultaneously, its amplitude of variation is also inequality.Various mechanical properties of rock are contrasted with the amplitude of fluid saturation variation and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows, different lithology strength of anisotropy and hole fracture density have certain correlativity, but are not complete correspondences.Vertically going up the hole crack grows than horizontal direction, rock is more loose, when increasing identical fluid saturation, the fluid of vertical direction filling is more, the range of decrease of the amplification of its Young modulus Poisson ratio and Biot coefficient is greater than horizontal direction, the mechanics difference of different directions is dwindled, and strength of anisotropy diminishes, and it is identical with the variation tendency of confined pressure with mechanical properties of rock.Under the less situation of Young modulus amplitude of variation, bigger Poisson ratio and Biot amplitude of variation may be arranged with a kind of lithology.Anisotropy is incomplete corresponding with the hole fracture density, and the difference that shows different directions hole fracture density is not to cause anisotropic unique factor; According to the petrographic thin section observations, the difference that puts in order of the mineralogical composition of different directions and mineral is obvious, and from the influence of mineral to mechanical properties of rock, the composition of mineral also is a factor that influences rock anisotropy intensity with aligning.
(4) pore fluid characteristic parameter
Reservoir rock is porous medium, but filling water, gas or wet goods fluid therebetween.The existence of fluid can make the mechanics parameter of rock change.Research blowhole fluid content changes the influence to the rock mechanics parameter, and the explanation of remaining oil saturation in the identification of oil-gas-water layer in the exploration, the four-dimensional earthquake is all had great importance.The inventor herein studies with the Changing Pattern of fluid type the mechanics parameter that prolongs the dissimilar rocks of oil field prolongation group reservoir.Test is found by experiment, and Poisson ratio convection cell type changes the most responsive, and the result as shown in Figure 9.
As shown in Figure 9, the Poisson ratio of the value that prolongs oil field prolongation group rock core oil-containing Poisson ratio in period when moisture, along with the increase of fluid saturation, it is big that its difference becomes gradually.This analysis result can provide support for the prediction of oil-water-layer in the exploration.
(5) pore pressure characteristic parameter
The hydrocarbon-bearing pool rock mass is porous medium, therebetween the filling pore fluid.Pore fluid applies pore pressure to rock skeleton.Pore pressure can a part of confined pressure of negative function on rock skeleton.The variation of pore pressure can cause the distortion of rock, causes the mechanics parameter of rock to change.The existence of pore pressure plays decisive action to mechanical property and the distortion of rock mass, such as mining area, oil field problems such as land subsidence on stream.
The inventor herein experimentizes by choosing 5 groups of rock samples, and every group of rock sample comprises two parallel sample, and one is the horizontal direction sampling, and one is the vertical direction sampling.In the test, sample is carried out sound measure simultaneously, the Young modulus of rock, Poisson ratio and Biot coefficient are to keep confined pressure, fluid saturation and temperature-resistant, and measurement calculates under the change pore pressure condition; Confined pressure, fluid saturation and temperature are chosen the oil reservoir real data, are respectively 20Mpa, 50% saturation water, 50% saturated oil, 50% saturation water, 40 ° of C.The pore pressure variation range is 0-18Mpa, every point of 3Mpa test.
Obtain following experimental result: along with the increase of pore pressure, the value of Young modulus and Poisson ratio reduces gradually, and the Biot coefficient increases gradually; The mechanical properties of rock of different lithology increases the variation tendency unanimity with temperature, but the amplitude difference.Mechanical property is opposite with its variation tendency with confined pressure with the variation tendency of pore pressure; Various mechanical properties of rock are contrasted with the amplitude of pore pressure variation and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows that the Young modulus of the rock of calcitic cementation, Poisson ratio and Biot coefficient are littler than the rock of argillaceous agglutination with the pore pressure amplitude of variation.Various lithological change amplitude trend and its hole fracture density correlativity are good, and be relatively poor with the correlativity of mineralogical composition.Comparative analysis is the result show, when pore pressure raise, the change of mechanical properties of rock mainly was subjected to the influence of hole fracture density.On this analysis foundation, in pore pressure increase process, the force-bearing situation of rock is analyzed.
The total stress σ that subsurface rock is subjected to is divided into two parts, and a part is born by rock skeleton, is called effective stress σ eA part is born by pore fluid, is called pore pressure P f, i.e. σ=σ e+ P f
According to the Terzaghi theorem, as can be known: (1) has all increased same amount when confined pressure and pore pressure, and then the volume of rock is almost constant; (2) in shear failure, if normal stress and pore pressure all increase same amount, shear resistance does not increase.Under the constant situation of confined pressure, along with the increase of pore pressure, the pressure that rock skeleton is subjected to is more and more littler, and then the crack closure degree of rock is more and more lower.By confined pressure to the influence of rock mechanics parameter as can be known, along with the increase of crack closure degree, the Young modulus of rock and Poisson ratio increase, the Biot coefficient reduces, and then can get, along with the increase of pore pressure, the Young modulus of rock and Poisson ratio reduce, and the Biot coefficient increases.
By rock sound difference with the change curve of pore pressure as can be known, along with the increase of pore pressure, the sound diversity ratio of the Young modulus of rock and Poisson ratio increases gradually; The mechanical properties of rock of different lithology increases the variation tendency unanimity with confined pressure, but the amplitude difference.Mineralogical composition, structure, cementing matter and the hole fracture density of the temperature variant amplitude of various mechanical properties of rock and rock are contrasted, and comparing result shows that different lithology sound difference amplitude and hole fracture density have certain correlativity; All relatively poor with the correlativity of mineralogical composition and cementing matter.Under the less situation of Young modulus amplitude of variation, bigger Poisson ratio amplitude of variation may be arranged with a kind of lithology.Microscopic mechanism according to sound difference (dynamically is high frequency-little displacement, static state is low frequency-big displacement), the hole crack should be the main cause that influences sound difference, result from this test, the sound diversity ratio increases (increase of hole fracture density) with pore pressure and increases, and this point matches with the mechanism of surrounding and watching of sound difference.The amplitude of variation of different lithology has certain correlativity with corresponding hole fissured structure, but not exclusively coincide, the reason that this explanation influences rock sound difference amplitude is comparatively complicated, except the hole fracture density is the principal element, also has other influence factors (mineralogical composition etc.).Satisfy binomial relation between rock mechanics parameter sound diversity ratio and the pore pressure, that is:
E d - E s E s = aP o 2 + bP o + C
&mu; d - &mu; s &mu; s = m P o 2 + np o + k
In the formula, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; P oBe pore pressure; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant.The regression formula of five kinds of rock samples of inventor herein's statistics is as shown in the table:
Figure BDA00003130394900573
Wherein, R 2Be related coefficient.
By rock mechanics parameter strength of anisotropy with the variation of pore pressure as can be known, the strength of anisotropy of rock mechanics parameter increases with the increase of pore pressure.The strength of anisotropy of Biot coefficient is negative, by the definition of strength of anisotropy as can be known, the Biot coefficient of vertical direction is greater than the Biot coefficient of horizontal direction, according to the hole fissured structure of level and vertical direction as can be known, the vertical hole fracture density that goes up is greatly to cause the vertical big main cause of Biot coefficient that goes up.The strength of anisotropy of the mechanical properties of rock of different lithology increases the variation tendency unanimity with pore pressure, but the amplitude difference.Various mechanical properties of rock are contrasted with the amplitude of pore pressure variation and mineralogical composition, structure, cementing matter and the hole fracture density of rock, comparing result shows, different lithology strength of anisotropy amplitude of variation and hole fracture density have certain correlativity, but are not complete correspondences.Vertically go up the hole crack and grow than horizontal direction, rock is more loose, when pore pressure increases, the crack is than the easier unlatching of horizontal direction, therefore, and when confined pressure increases, differences of mechanical properties between different directions can increase gradually, and it is identical that its trend and anisotropy vary with temperature trend.Under the less situation of Young modulus amplitude of variation, bigger Poisson ratio and Biot amplitude of variation may be arranged with a kind of lithology.Anisotropy is incomplete corresponding with the hole fracture density, and the difference that shows different directions hole fracture density is not to cause anisotropic unique factor.According to the petrographic thin section observations, the difference that puts in order of the mineralogical composition of different directions and mineral is obvious, and from the influence of mineral to mechanical properties of rock, the composition of mineral also is a factor that influences rock anisotropy intensity with aligning.
To sum up, hydrocarbon-bearing pool rock mass mechanics character is to be determined by the build-in attribute of rock (rock composition, structure, structure, hole FRACTURE CHARACTERISTICS) and external environment condition (confined pressure, temperature, fluid and pore pressure).The Young modulus of hydrocarbon-bearing pool rock mass and the value of Poisson ratio increase along with the increase of quartz content and cementing matter intensity; Rock mass structure is more loose, and educate all the more in the hole crack, and the value of its Young modulus and Poisson ratio is more little, and the value of Biot coefficient is more big.Along with the increase of temperature, it is big that the hole in the rock and crevice volume become, and causes the Young modulus of rock mass and the value of Poisson ratio to diminish, and it is big that the value of Biot coefficient becomes; Along with the increase of confined pressure, the hole in the rock and part crack pressurized closure cause the hole crack density of rock to diminish, and it is big that the Young modulus of rock and Poisson ratio become, and the Biot coefficient value dwindles; During fill fluid, along with the increase of fluid saturation, the coupling between rock particles and the hole improves in the rock, and fluid has been born the pressure of partial action on rock skeleton simultaneously, causes Young modulus, Poisson ratio and the Biot coefficient value of rock to become big; The reacting condition of Poisson ratio convection cell type is the most responsive, and the value of the Poisson ratio of the value of the Poisson ratio during the rock saturated oil during than saturation water is big; Along with the increase of pore pressure, the net pressure that acts on the rock skeleton reduces, thereby causes the Young modulus of rock, Poisson ratio to diminish, and it is big that the Biot coefficient value becomes.
Rock is porous medium, and the hole that exists in the rock and crack are the basic reasons that causes rock sound rock mechanics parameter difference.Increase along with temperature, the sample internal moisture is sought the crack and is overflowed to the infiltration of rock mass outside in a large amount of evaporation processes, thereby make hole widen lengthening, simultaneously, action of high temperature makes the volumetric expansion of sample endoparticle, also make the sample internal crack increase, therefore along with the increase of temperature, the sound diversity ratio of rock mechanics parameter is increasing; When confined pressure increased, it was closed gradually to get hole and microfissure pressurized in the rock, and the difference of the dynamic and static elastic parameter of rock reduces; Contain the increase of liquid saturation degree along with rock, increasing by fluid-filled volume of voids, the hole that ultrasound wave can be walked around is also just more and more littler, the effect of rock interior microfissure system weakens, coupling between rock particles improves, therefore along with the increase that contains the liquid saturation degree, the sound diversity ratio of Young modulus and Poisson ratio is more and more littler; When pore pressure increased, the quiet confined pressure that acts on the rock skeleton reduced gradually, caused blowhole and microfracture to be opened, thereby caused the sound difference of rock mechanics parameter to increase gradually.
With reference to Figure 10, show the structured flowchart of the device embodiment of the underground former bit model recovery of a kind of hydrocarbon-bearing pool rock mass mechanics of the present invention, specifically can comprise as lower module:
Rock mechanics parameter measurement module 111 is used for measuring the rock mechanics parameter of hydrocarbon-bearing pool;
Influence parametric statistics module 112, be used for the rock mass mechanics property effect parameter of the described hydrocarbon-bearing pool of statistics;
Strength of anisotropy computing module 113 is for the strength of anisotropy of calculating described rock mechanics parameter;
First concerns acquisition module 114, for the incidence relation of rock mechanics parameter and the rock mass mechanics property effect parameter of obtaining described hydrocarbon-bearing pool;
Second concerns acquisition module 115, is used for obtaining the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy;
Former bit model recovers module 116, be used for adopting described rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and, the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
In specific implementation, described rock mechanics parameter can comprise Young modulus, Poisson ratio and be used for characterizing the Biot coefficient of poroelasticity parameter, in this case, the strength of anisotropy of described rock mechanics parameter can comprise Young modulus strength of anisotropy, Poisson ratio strength of anisotropy and Biot coefficient strength of anisotropy.
In a preferred embodiment of the present invention, described rock mass mechanics property effect parameter can comprise the external action parameter, in this case, the step of the described incidence relation that obtains rock mass mechanics property effect parameter and strength of anisotropy can for,
Obtain the incidence relation of described external action parameter and strength of anisotropy.
As the concrete a kind of example used of the embodiment of the invention, the strength of anisotropy of described rock mechanics parameter is calculated by following formula and is obtained:
E k = E h - E v E h ;
&mu; k = &mu; h - &mu; v &mu; h ;
Wherein, E kBe the Young modulus strength of anisotropy; E hBe the horizontal direction Young modulus; E vBe the vertical direction Young modulus; μ kBe the Poisson ratio strength of anisotropy; μ hBe the horizontal direction Poisson ratio; μ vBe the vertical direction Poisson ratio.
In specific implementation, described external action parameter can comprise the temperature profile parameter, confined pressure characteristic parameter, pore fluid characteristic parameter, pore pressure characteristic parameter; In this case, the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy can comprise:
The strength of anisotropy of described rock mechanics parameter becomes big along with the increase of temperature profile parameter;
And/or,
The strength of anisotropy of described rock mechanics parameter reduces along with the increase of confined pressure characteristic parameter;
And/or,
The strength of anisotropy of described rock mechanics parameter reduces along with the increase of pore fluid characteristic parameter; Wherein, described pore fluid characteristic parameter is fluid saturation;
And/or,
The strength of anisotropy of described rock mechanics parameter becomes big along with the increase of pore pressure characteristic parameter.
In a preferred embodiment of the present invention, the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool can comprise the internal influence parameter, and described internal influence parameter can comprise: rock composition characteristic parameter, rock texture characteristic parameter, the rock structure characteristic parameter, blowhole FRACTURE CHARACTERISTICS parameter.
In specific implementation, described rock mechanics parameter can comprise rock elasticity parameter and rock strength parameter, and described rock elasticity parameter can comprise Poisson ratio, Young modulus and poroelasticity coefficient; Described rock composition characteristic parameter can comprise rock particles characteristic parameter and chink characteristic parameter, and in this case, the incidence relation of the rock mechanics parameter of described hydrocarbon-bearing pool and rock mass mechanics property effect parameter can comprise:
If quartz content increases in the described rock particles characteristic parameter, then Young modulus rises;
And/or,
If described rock particles characteristic parameter andesine content increases, then Young modulus reduces;
And/or,
If described rock particles characteristic parameter CLAY MINERALS AND THEIR SIGNIFICANCE content increases, then Poisson ratio reduces;
And/or,
If described chink characteristic parameter is the calcitic cementation rock, then its Young modulus is greater than the Young modulus of other type of rock;
And/or,
If the rock that described rock texture characteristic parameter is the crystal class origin cause of formation, then its Young modulus and Poisson ratio are higher than the rock that the rock texture characteristic parameter is the cementing matter diagenesis;
And/or,
If described rock structure characteristic parameter bored sample by horizontal direction, its Young modulus and Poisson ratio Young modulus and the Poisson ratio of boring sample greater than vertical direction then;
And/or,
If blowhole FRACTURE CHARACTERISTICS parameter is the hole fracture spaces, then along with the increase of hole fracture spaces, the Young modulus of rock and the value of Poisson ratio reduce.
In another kind of preferred embodiment of the present invention, the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool can also comprise the external action parameter; Described external action parameter can comprise: temperature profile parameter, confined pressure characteristic parameter, pore fluid characteristic parameter, pore pressure characteristic parameter.
In this case, the incidence relation of the rock mechanics parameter of described hydrocarbon-bearing pool and rock mass mechanics property effect parameter can comprise:
If described temperature profile parameter increases, then the hole fracture density of rock increases, and the Young modulus of rock reduces, and the Poisson ratio of rock reduces, and the poroelasticity coefficient of rock increases;
And/or,
If described confined pressure characteristic parameter increases, then the hole fracture density of rock reduces, and the Young modulus of rock and Poisson ratio increase, and the poroelasticity coefficient of rock reduces;
And/or,
If described pore fluid characteristic parameter is fluid saturation, then along with the increase of fluid saturation, the Young modulus of rock, Poisson ratio and poroelasticity coefficient all increase; And, the Poisson ratio the when Poisson ratio in rock core oil-containing period is moisture greater than rock core, along with the increase of fluid saturation, its Poisson ratio difference increases;
And/or,
If described pore pressure characteristic parameter increases, then the Young modulus of rock, Poisson ratio reduce the increase of poroelasticity coefficient.
In specific implementation, described rock mechanics parameter measurement module 111 can comprise following submodule:
The static measurement submodule is used for adopting static measurement method to obtain static rock mechanics parameter;
The kinetic measurement submodule is used for adopting dynamic measurement method to obtain dynamic rock mechanics parameter;
In this case, described device also comprises:
Sound diversity ratio computing module is used for adopting described static rock mechanics parameter and dynamic rock mechanics parameter, the sound diversity ratio of computing rock mechanics parameter;
Described former bit model recovers module 114 and also is used for adopting described static rock mechanics parameter, dynamic rock mechanics parameter, the sound diversity ratio of rock mechanics parameter, and, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
In specific implementation, described static rock mechanics parameter and dynamic rock mechanics parameter obtain through synchro measure, described static rock mechanics parameter can comprise static rock elasticity parameter and static rock strength parameter, and described dynamic rock mechanics parameter can comprise dynamic rock elasticity parameter and dynamic rock strength parameter; Described static rock elasticity parameter can comprise static Poisson ratio, static Young modulus and static poroelasticity coefficient etc., and described dynamic rock elasticity parameter can comprise dynamic Poisson ratio, kinetic Youngs modulus and dynamic hole elasticity coefficient etc.
As the concrete example of using of the embodiment of the invention, described static Poisson ratio μ can calculate by following formula and obtain:
&mu; = &Delta;&epsiv; 1 &Delta;&epsiv; 2
Described dynamic Poisson ratio μ can calculate by following formula and obtain:
&mu; = &Delta;t s 2 - 2 &Delta;t p 2 2 ( &Delta;t s 2 - &Delta;t p 2 ) = V p 2 - 2 V s 2 2 ( V p 2 - V s 2 )
Wherein, described Δ ε 1Axial strain increment for rock sample; Δ ε 2Transverse strain increment for rock sample; Δ t sBe the shear wave slowness of rock sample, μ s/ m; Δ t pBe the compressional wave time difference of rock sample, μ s/ m; V sBe the shear wave velocity of rock sample, m/s; V pBe the velocity of longitudinal wave of rock sample, m/s;
Described static Young modulus E can calculate by following formula and obtain:
E = F / S &Delta;L / L
Described kinetic Youngs modulus E can calculate by following formula and obtain:
E = &rho; V s 2 3 ( V p V s ) 2 - 4 ( V p V s ) 2 - 1
Wherein, described ρ is the density of rock, g/cm 3F is the suffered power of test rock sample; S is the sectional area of rock sample; Δ L is that rock sample is being subjected to the elongation of force direction (or shortening amount); L is that rock sample is in the former length that is subjected on the force direction; V SBe the shear wave velocity of rock sample, m/s; Vp is the velocity of longitudinal wave of rock sample, m/s;
Described static poroelasticity coefficient B can be calculated by following formula and obtain:
B = 1 - K dry K s
Wherein, described K DryBe the bulk modulus of dry rock, Ks is the bulk modulus of the mineral of composition rock;
Described dynamic hole elasticity coefficient B can calculate by following formula and obtain:
B = 1 - K dry K s = 1 - 1 &Delta;V s 1 &Delta;V s - &Delta;V f = 1 - &Delta;V s - &Delta;V f &Delta;V s = &Delta;V f &Delta;V s
Wherein, described K DryBe the bulk modulus of dry rock, Ks is the bulk modulus of the mineral of composition rock, Δ V fBe pore fluid fluctuations in discharge amount, Δ V sBe the rock skeleton volume change, Δ V mBe rock forming mineral particle volume variable quantity;
The sound diversity ratio of described rock mechanics parameter can comprise Young modulus sound diversity ratio and Poisson ratio sound diversity ratio, and described Young modulus sound diversity ratio can be calculated by following formula and obtain:
E c = E d - E s E s
Wherein, described E dBe kinetic Youngs modulus; E sBe static Young modulus;
Described Poisson ratio sound diversity ratio can be calculated by following formula and obtain:
&mu; c = &mu; d - &mu; s &mu; s ;
Wherein, described μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio.
In a kind of preferred exemplary of practical application of the present invention, satisfy the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the temperature profile parameter:
E d - E s E s = aT 2 + bT + C ;
&mu; d - &mu; s &mu; s = mT 2 + nT + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; T is temperature; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant;
Satisfy the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the confined pressure characteristic parameter:
E d - E s E s = aP c 2 + bP c + C ;
&mu; d - &mu; s &mu; s = m P c 2 + nP c + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; P cBe confined pressure; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant;
Satisfy the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the fluid saturation:
E d - E s E s = a 1 S w 2 + b 1 S w + C 1 ;
&mu; d - &mu; s &mu; s = m 1 S w 2 + n 1 S w + k 1 ;
E d - E s E s = a 2 S o 2 + b 2 S o + C 2 ;
&mu; d - &mu; s &mu; s = m 2 S o 2 + n 2 S o + k 2 ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; S wBe water saturation; S oBe oil saturation; a 1, a 2Be the quadratic term coefficient; b 1, b 2Be coefficient once; C 1, C 2Be constant; m 1, m 2Be the quadratic term coefficient; n 1, n 2Be coefficient once; k 1, k 2Be constant;
Satisfy binomial relation shown in following formula between the sound diversity ratio of described rock mechanics parameter and the pore pressure:
E d - E s E s = a P o 2 + bP o + C ;
&mu; d - &mu; s &mu; s = m P o 2 + nP o + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; P oBe pore pressure; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant.
For device embodiment, because it is similar substantially to method embodiment, so description is fairly simple, relevant part gets final product referring to the part explanation of method embodiment.
More than the device that recovers of method that the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics provided by the present invention is recovered and the underground former bit model of a kind of hydrocarbon-bearing pool rock mass mechanics, be described in detail, used specific case among the inventor herein principle of the present invention and embodiment are set forth, the explanation of above embodiment just is used for helping to understand method of the present invention and core concept thereof; Simultaneously, for one of ordinary skill in the art, according to thought of the present invention, the part that all can change in specific embodiments and applications, in sum, this description should not be construed as limitation of the present invention.

Claims (12)

1. the method that the underground former bit model of hydrocarbon-bearing pool rock mass mechanics recovers is characterized in that, comprising:
Measure the rock mechanics parameter of hydrocarbon-bearing pool;
Add up the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool;
Calculate the strength of anisotropy of described rock mechanics parameter;
Obtain the incidence relation of rock mechanics parameter and the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool;
Obtain the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy;
Adopt described rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and, the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
2. method according to claim 1 is characterized in that, described rock mechanics parameter comprises Young modulus, Poisson ratio and Biot coefficient, and described Biot coefficient is used for characterizing the poroelasticity parameter; The strength of anisotropy of described rock mechanics parameter comprises Young modulus strength of anisotropy, Poisson ratio strength of anisotropy and Biot coefficient strength of anisotropy.
3. method according to claim 1 and 2 is characterized in that, described rock mass mechanics property effect parameter comprises the external action parameter, and the described step of obtaining the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy is,
Obtain the incidence relation of described external action parameter and strength of anisotropy.
4. method according to claim 2 is characterized in that, the strength of anisotropy of described rock mechanics parameter is calculated by following formula and obtained:
E k = E h - E v E h ;
&mu; k = &mu; h - &mu; v &mu; h ;
Wherein, E kBe the Young modulus strength of anisotropy; E hBe the horizontal direction Young modulus; E vBe the vertical direction Young modulus; μ kBe the Poisson ratio strength of anisotropy; μ hBe the horizontal direction Poisson ratio; μ vBe the vertical direction Poisson ratio.
5. method according to claim 3 is characterized in that, described external action parameter comprises the temperature profile parameter, confined pressure characteristic parameter, pore fluid characteristic parameter, pore pressure characteristic parameter; The incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy comprises:
The strength of anisotropy of described rock mechanics parameter becomes big along with the increase of temperature profile parameter;
And/or,
The strength of anisotropy of described rock mechanics parameter reduces along with the increase of confined pressure characteristic parameter;
And/or,
The strength of anisotropy of described rock mechanics parameter reduces along with the increase of pore fluid characteristic parameter; Wherein, described pore fluid characteristic parameter is fluid saturation;
And/or,
The strength of anisotropy of described rock mechanics parameter becomes big along with the increase of pore pressure characteristic parameter.
6. method according to claim 1 is characterized in that, the rock mass mechanics property effect parameter of described hydrocarbon-bearing pool comprises the internal influence parameter, described internal influence parameter comprises the rock composition characteristic parameter, the rock texture characteristic parameter, rock structure characteristic parameter, blowhole FRACTURE CHARACTERISTICS parameter;
Described rock mechanics parameter comprises rock elasticity parameter and rock strength parameter, and described rock elasticity parameter comprises Poisson ratio, Young modulus;
Described rock composition characteristic parameter comprises rock particles characteristic parameter and chink characteristic parameter, and the incidence relation of the rock mechanics parameter of described hydrocarbon-bearing pool and rock mass mechanics property effect parameter comprises:
If quartz content increases in the described rock particles characteristic parameter, then Young modulus rises;
And/or,
If described rock particles characteristic parameter andesine content increases, then Young modulus reduces;
And/or,
If described rock particles characteristic parameter CLAY MINERALS AND THEIR SIGNIFICANCE content increases, then Poisson ratio reduces;
And/or,
If described chink characteristic parameter is the calcitic cementation rock, then its Young modulus is greater than the Young modulus of other type of rock;
And/or,
If the rock that described rock texture characteristic parameter is the crystal class origin cause of formation, then its Young modulus and Poisson ratio are higher than the rock that the rock texture characteristic parameter is the cementing matter diagenesis;
And/or,
If described rock structure characteristic parameter bored sample by horizontal direction, its Young modulus and Poisson ratio Young modulus and the Poisson ratio of boring sample greater than vertical direction then;
And/or,
If blowhole FRACTURE CHARACTERISTICS parameter is the hole fracture spaces, then along with the increase of hole fracture spaces, the Young modulus of rock and the value of Poisson ratio reduce.
7. method according to claim 5 is characterized in that, the incidence relation of the rock mechanics parameter of described hydrocarbon-bearing pool and rock mass mechanics property effect parameter also comprises:
If described temperature profile parameter increases, then the hole fracture density of rock increases, and the Young modulus of rock reduces, and the Poisson ratio of rock reduces, and the poroelasticity coefficient of rock increases;
And/or,
If described confined pressure characteristic parameter increases, then the hole fracture density of rock reduces, and the Young modulus of rock and Poisson ratio increase, and the poroelasticity coefficient of rock reduces;
And/or,
If described pore fluid characteristic parameter is fluid saturation, then along with the increase of fluid saturation, the Young modulus of rock, Poisson ratio and poroelasticity coefficient all increase; And, the Poisson ratio the when Poisson ratio in rock core oil-containing period is moisture greater than rock core, along with the increase of fluid saturation, its Poisson ratio difference increases;
And/or,
If described pore pressure characteristic parameter increases, then the Young modulus of rock, Poisson ratio reduce the increase of poroelasticity coefficient.
8. method according to claim 1 is characterized in that, the step of the rock mechanics parameter of described measurement hydrocarbon-bearing pool comprises:
Adopt static measurement method to obtain static rock mechanics parameter;
Adopt dynamic measurement method to obtain dynamic rock mechanics parameter;
Wherein, described static rock mechanics parameter and dynamic rock mechanics parameter obtain through synchro measure;
Described method also comprises:
Adopt described static rock mechanics parameter and dynamic rock mechanics parameter, the sound diversity ratio of computing rock mechanics parameter;
Described step of carrying out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics according to the rock mechanics parameter after adjusting is,
Adopt described static rock mechanics parameter, dynamic rock mechanics parameter, the sound diversity ratio of rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
9. method according to claim 8 is characterized in that, described static rock mechanics parameter comprises static rock elasticity parameter and static rock strength parameter, and described dynamic rock mechanics parameter comprises dynamic rock elasticity parameter and dynamic rock strength parameter; Described static rock elasticity parameter comprises static Poisson ratio, static Young modulus and static poroelasticity coefficient, and described dynamic rock elasticity parameter comprises dynamic Poisson ratio, kinetic Youngs modulus and dynamic hole elasticity coefficient.
10. method according to claim 9 is characterized in that, described static Poisson ratio μ calculates by following formula and obtains:
&mu; = &Delta;&epsiv; 1 &Delta;&epsiv; 2
Described dynamic Poisson ratio μ calculates by following formula and obtains:
&mu; = &Delta;t s 2 - 2 &Delta;t p 2 2 ( &Delta;t s 2 - &Delta;t p 2 ) = V p 2 - 2 V s 2 2 ( V p 2 - V s 2 )
Wherein, described Δ ε 1Axial strain increment for rock sample; Δ ε 2Transverse strain increment for rock sample; Δ t sBe the shear wave slowness of rock sample, μ s/ m; Δ t pBe the compressional wave time difference of rock sample, μ s/ m; V sBe the shear wave velocity of rock sample, m/s; V pBe the velocity of longitudinal wave of rock sample, m/s;
Described static Young modulus E calculates by following formula and obtains:
E = F / S &Delta;L / L
Described kinetic Youngs modulus E calculates by following formula and obtains:
E = &rho;V s 2 3 ( V p V s ) 2 - 4 ( V p V s ) 2 - 1
Wherein, described ρBe the density of rock, g/cm 3F is the suffered power of test rock sample; S is the sectional area of rock sample; △ L is that rock sample is being subjected to the elongation of force direction (or shortening amount); L is that rock sample is in the former length that is subjected on the force direction; V SBe the shear wave velocity of rock sample, m/s; Vp is the velocity of longitudinal wave of rock sample, m/s;
Described static poroelasticity coefficient B is calculated by following formula and is obtained:
B = 1 - K dry K s
Wherein, described K DryBe the bulk modulus of dry rock, Ks is the bulk modulus of the mineral of composition rock;
Described dynamic hole elasticity coefficient B calculates by following formula and obtains:
B = 1 - K dry K s = 1 - 1 &Delta;V s 1 &Delta;V s - &Delta;V f = 1 - &Delta;V s - &Delta;V f &Delta;V s = &Delta;V f &Delta;V s
Wherein, described K DryBe the bulk modulus of dry rock, Ks is the bulk modulus of the mineral of composition rock, Δ V fBe pore fluid fluctuations in discharge amount, Δ V sBe the rock skeleton volume change, Δ V mBe rock forming mineral particle volume variable quantity;
The sound diversity ratio of described rock mechanics parameter comprises Young modulus sound diversity ratio and Poisson ratio sound diversity ratio, and described Young modulus sound diversity ratio is calculated by following formula and obtained:
E c = E d - E s E s
Wherein, described E dBe kinetic Youngs modulus; E sBe static Young modulus;
Described Poisson ratio sound diversity ratio is calculated by following formula and is obtained:
&mu; c = &mu; d - &mu; s &mu; s ;
Wherein, described μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio.
11. method according to claim 10 is characterized in that, satisfies the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the temperature profile parameter:
E d - E s E s = aT 2 + bT + C ;
&mu; d - &mu; s &mu; s = mT 2 + nT + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; T is temperature; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant;
Satisfy the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the confined pressure characteristic parameter:
E d - E s E s = aP c 2 + bP c + C ;
&mu; d - &mu; s &mu; s = mP c 2 + nP c + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; P cBe confined pressure; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant;
Satisfy the binomial relation shown in the following formula between the sound diversity ratio of described rock mechanics parameter and the fluid saturation:
E d - E s E s = a 1 S w 2 + b 1 S w + C 1 ;
&mu; d - &mu; s &mu; s = m 1 S w 2 + n 1 S w + k 1 ;
E d - E s E s = a 2 S o 2 + b 2 S o + C 2 ;
&mu; d - &mu; s &mu; s = m 2 S o 2 + n 2 S o + k 2 ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; S wBe water saturation; S oBe oil saturation; a 1, a 2Be the quadratic term coefficient; b 1, b 2Be coefficient once; C 1, C 2Be constant; m 1, m 2Be the quadratic term coefficient; n 1, n 2Be coefficient once; k 1, k 2Be constant;
Satisfy binomial relation shown in following formula between the sound diversity ratio of described rock mechanics parameter and the pore pressure:
E d - E s E s = aP o 2 + bP o + C ;
&mu; d - &mu; s &mu; s = mP o 2 + nP o + k ;
Wherein, E dBe kinetic Youngs modulus; E sBe static Young modulus; μ dBe dynamic Poisson ratio; μ sBe static Poisson ratio; P oBe pore pressure; A is the quadratic term coefficient; B is coefficient once; C is constant; M is the quadratic term coefficient; N is coefficient once; K is constant.
12. the device that the underground former bit model of hydrocarbon-bearing pool rock mass mechanics recovers is characterized in that,
Comprise:
Rock mechanics parameter measurement module is used for measuring the rock mechanics parameter of hydrocarbon-bearing pool;
Influence the parametric statistics module, be used for the rock mass mechanics property effect parameter of the described hydrocarbon-bearing pool of statistics;
The strength of anisotropy computing module is for the strength of anisotropy of calculating described rock mechanics parameter;
First concerns acquisition module, for the incidence relation of rock mechanics parameter and the rock mass mechanics property effect parameter of obtaining described hydrocarbon-bearing pool;
Second concerns acquisition module, is used for obtaining the incidence relation of described rock mass mechanics property effect parameter and strength of anisotropy;
Former bit model recovers module, be used for adopting described rock mechanics parameter, the strength of anisotropy of rock mechanics parameter, the incidence relation of the rock mechanics parameter of hydrocarbon-bearing pool and rock mass mechanics property effect parameter, and, the incidence relation of rock mass mechanics property effect parameter and strength of anisotropy carries out the recovery of the underground former bit model of hydrocarbon-bearing pool rock mass mechanics.
CN201310157416.XA 2013-04-28 2013-04-28 A kind of method that hydrocarbon-bearing pool rock mass mechanics ground in-situ model recovers and device Expired - Fee Related CN103257081B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310157416.XA CN103257081B (en) 2013-04-28 2013-04-28 A kind of method that hydrocarbon-bearing pool rock mass mechanics ground in-situ model recovers and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310157416.XA CN103257081B (en) 2013-04-28 2013-04-28 A kind of method that hydrocarbon-bearing pool rock mass mechanics ground in-situ model recovers and device

Publications (2)

Publication Number Publication Date
CN103257081A true CN103257081A (en) 2013-08-21
CN103257081B CN103257081B (en) 2015-09-02

Family

ID=48961131

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310157416.XA Expired - Fee Related CN103257081B (en) 2013-04-28 2013-04-28 A kind of method that hydrocarbon-bearing pool rock mass mechanics ground in-situ model recovers and device

Country Status (1)

Country Link
CN (1) CN103257081B (en)

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200039A (en) * 2014-09-17 2014-12-10 中国石油大学(华东) Quantitative forecasting method of tectonic fissure occurrence
CN104297783A (en) * 2014-06-30 2015-01-21 中国石油天然气集团公司 Interpretation method and system for hydrofracture micro-seismic event
CN104535450A (en) * 2014-12-31 2015-04-22 中国石油天然气股份有限公司 Well drilling site marine facies carbonate rock stratum system lithology analysis method and device
CN104677734A (en) * 2015-02-12 2015-06-03 中国石油大学(华东) Multi-scale information integration tight formation brittleness index measuring method
CN104865370A (en) * 2015-04-22 2015-08-26 中国矿业大学 Dual-purpose gripper for transfusion under ultrasonic load and desorption experiment
CN105158341A (en) * 2015-09-18 2015-12-16 重庆大学 Sound-speed-based method for identifying material of transformer winding
CN105424495A (en) * 2015-11-06 2016-03-23 太原理工大学 Device and method for testing plane stress state of coal rock sheet under thermal-mechanical coupling condition
CN105761158A (en) * 2016-03-03 2016-07-13 长江大学 Analysis method for rock salt underground gas storage limit operation pressure
CN106289585A (en) * 2015-05-19 2017-01-04 中国石油大学(华东) A kind of Visco elastic strain surveys the measuring method of the test device of on-the-spot crustal stress
CN106383172A (en) * 2016-10-10 2017-02-08 三峡大学 Surrounding rock damage prediction method based on energy release coefficient
CN106682757A (en) * 2015-12-31 2017-05-17 中国石油天然气股份有限公司 Method and device for determining collapse pressure of borehole wall
CN106772673A (en) * 2016-11-29 2017-05-31 西南石油大学 A kind of shale reservoir crustal stress predicts modeling method
CN107016156A (en) * 2016-01-28 2017-08-04 波音公司 The method that test data for failure mode in composite layer simplifies
CN107036907A (en) * 2017-04-07 2017-08-11 北京科技大学 A kind of fissure rock intensive analysis Forecasting Methodology
CN107748205A (en) * 2017-12-05 2018-03-02 中国空气动力研究与发展中心计算空气动力研究所 A kind of elastic constant measurement method varied with temperature
CN107917865A (en) * 2016-10-11 2018-04-17 中国石油化工股份有限公司 A kind of tight sandstone reservoir multi-parameter Permeability Prediction method
CN108195669A (en) * 2017-10-09 2018-06-22 中国石油化工股份有限公司 The correction of oil reservoir Under Confining Pressure static mechanical parameter and Forecasting Methodology
CN108519317A (en) * 2018-04-24 2018-09-11 钦州学院 Rock stress-seepage coupling test device under uniaxial direct tensile load
CN108732048A (en) * 2018-05-31 2018-11-02 长安大学 A kind of elastoplasticity yield point stress of graded broken stone repeated-load test determines method
CN108843303A (en) * 2018-07-19 2018-11-20 西南石油大学 A kind of casing damage in oil-water well prediction technique based on mud stone creep model
CN109033672A (en) * 2018-08-09 2018-12-18 中国石油天然气股份有限公司 The dynamic crack determination method and device of pore constriction network model in displacement simulation
CN109447944A (en) * 2018-09-21 2019-03-08 中国石油天然气股份有限公司 The Lithofacies Identification method and system of carbonate rock
CN109828031A (en) * 2019-02-15 2019-05-31 西南石油大学 Rock brittleness evaluation method and device
CN110006568A (en) * 2019-04-12 2019-07-12 中国地质科学院地质力学研究所 It is a kind of to obtain the method for three-dimensional ground stress using core and its obtain system
CN110133724A (en) * 2019-05-06 2019-08-16 中国石油天然气集团有限公司 A kind of application array sonic log data progress gas-bearing formation knowledge method for distinguishing
CN110836929A (en) * 2019-11-26 2020-02-25 中国科学院武汉岩土力学研究所 Jointed rock mass vibration attenuation physical simulation test method and device
CN111089904A (en) * 2019-12-20 2020-05-01 西南石油大学 Indoor measurement wave velocity confining pressure and dispersion correction method considering formation characteristics
CN111175128A (en) * 2020-01-13 2020-05-19 西南石油大学 Compact sandstone tensile strength prediction method based on spectral analysis
CN111208565A (en) * 2020-03-04 2020-05-29 中国石油大学(北京) KT model-based hole seam parameter inversion method and device and storage medium
CN112001055A (en) * 2019-11-07 2020-11-27 中海石油(中国)有限公司 Low-amplitude thin oil reservoir water content prediction method based on microstructure
CN112324429A (en) * 2020-08-12 2021-02-05 中国石油天然气股份有限公司 Method for determining reasonable production pressure difference of water-gas reservoir
CN112673147A (en) * 2018-09-13 2021-04-16 沙特阿拉伯石油公司 Downhole tool for fracturing a hydrocarbon containing formation
CN113361131A (en) * 2021-06-25 2021-09-07 中国矿业大学 Static mechanical parameter calculation method for primary structure coal
US11175207B2 (en) 2016-12-20 2021-11-16 Petrochina Company Limited Method for obtaining conversion relationship between dynamic and static elastic parameters
CN114113083A (en) * 2021-11-10 2022-03-01 中国石油大学(北京) Method, device, medium, apparatus and program product for evaluating heterogeneity of conglomerate
CN114526066A (en) * 2022-01-24 2022-05-24 西南石油大学 Method for compacting and recovering continental facies shale oil reservoir
CN115266514A (en) * 2022-05-11 2022-11-01 中国石油大学(华东) Dynamic evaluation device and method for rock mechanical parameters in high-pressure fluid injection process
CN116066092A (en) * 2023-02-03 2023-05-05 西南石油大学 Shale formation collapse pressure increment prediction method based on hardness experiment

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1034064A (en) * 1987-10-14 1989-07-19 阿莫科公司 Be used for detecting the method and apparatus of the elastic anisotropy of measuring the geology rock sample
CN101886996A (en) * 2010-05-31 2010-11-17 付志亮 Triaxial compression rheological test system capable of simulating engineering geological environment
US20110015907A1 (en) * 2009-07-20 2011-01-20 Crawford Brian R Petrophysical Method For Predicting Plastic Mechanical Properties In Rock Formations
CN102011580A (en) * 2010-11-08 2011-04-13 西南石油大学 Method for predicting failure pressure of reservoir with acid damage
CN102607959A (en) * 2012-03-28 2012-07-25 中国石油大学(华东) Experimental device and method for measuring rock mechanics parameters under action of ultrasonic wave and chemistry

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1034064A (en) * 1987-10-14 1989-07-19 阿莫科公司 Be used for detecting the method and apparatus of the elastic anisotropy of measuring the geology rock sample
US20110015907A1 (en) * 2009-07-20 2011-01-20 Crawford Brian R Petrophysical Method For Predicting Plastic Mechanical Properties In Rock Formations
CN101886996A (en) * 2010-05-31 2010-11-17 付志亮 Triaxial compression rheological test system capable of simulating engineering geological environment
CN102011580A (en) * 2010-11-08 2011-04-13 西南石油大学 Method for predicting failure pressure of reservoir with acid damage
CN102607959A (en) * 2012-03-28 2012-07-25 中国石油大学(华东) Experimental device and method for measuring rock mechanics parameters under action of ultrasonic wave and chemistry

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘斌等: "不同围压下岩石中泊松比的各向异性", 《地球物理学报》, vol. 45, no. 6, 30 November 2002 (2002-11-30) *
吴文娟等: "超低渗油气藏非对称压裂数值模拟理论及应用", 《北京大学学报(自然科学版)》, vol. 48, no. 6, 30 November 2012 (2012-11-30), pages 895 - 901 *
韩林等: "高压油气藏对砂岩力学特性影响的试验研究", 《岩石力学与工程学报》, vol. 28, no. 5, 31 May 2009 (2009-05-31), pages 965 - 972 *

Cited By (61)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104297783A (en) * 2014-06-30 2015-01-21 中国石油天然气集团公司 Interpretation method and system for hydrofracture micro-seismic event
CN104200039A (en) * 2014-09-17 2014-12-10 中国石油大学(华东) Quantitative forecasting method of tectonic fissure occurrence
CN104535450A (en) * 2014-12-31 2015-04-22 中国石油天然气股份有限公司 Well drilling site marine facies carbonate rock stratum system lithology analysis method and device
CN104677734A (en) * 2015-02-12 2015-06-03 中国石油大学(华东) Multi-scale information integration tight formation brittleness index measuring method
CN104865370A (en) * 2015-04-22 2015-08-26 中国矿业大学 Dual-purpose gripper for transfusion under ultrasonic load and desorption experiment
CN104865370B (en) * 2015-04-22 2016-06-15 中国矿业大学 Seepage flow, desorption experiment Two-purpose paper clip holder under a kind of ultrasonic pressing
CN106289585B (en) * 2015-05-19 2019-04-09 中国石油大学(华东) A kind of measurement method of the test device of the live crustal stress of Visco elastic strain survey
CN106289585A (en) * 2015-05-19 2017-01-04 中国石油大学(华东) A kind of Visco elastic strain surveys the measuring method of the test device of on-the-spot crustal stress
CN105158341A (en) * 2015-09-18 2015-12-16 重庆大学 Sound-speed-based method for identifying material of transformer winding
CN105158341B (en) * 2015-09-18 2018-05-25 重庆大学 Transformer winding material discrimination method based on the velocity of sound
CN105424495A (en) * 2015-11-06 2016-03-23 太原理工大学 Device and method for testing plane stress state of coal rock sheet under thermal-mechanical coupling condition
CN106682757B (en) * 2015-12-31 2020-04-10 中国石油天然气股份有限公司 Method and device for determining borehole wall collapse pressure
CN106682757A (en) * 2015-12-31 2017-05-17 中国石油天然气股份有限公司 Method and device for determining collapse pressure of borehole wall
CN107016156A (en) * 2016-01-28 2017-08-04 波音公司 The method that test data for failure mode in composite layer simplifies
CN105761158A (en) * 2016-03-03 2016-07-13 长江大学 Analysis method for rock salt underground gas storage limit operation pressure
CN106383172A (en) * 2016-10-10 2017-02-08 三峡大学 Surrounding rock damage prediction method based on energy release coefficient
CN106383172B (en) * 2016-10-10 2018-09-28 三峡大学 A kind of rock convergence measure prediction technique based on energy release coefficient
CN107917865A (en) * 2016-10-11 2018-04-17 中国石油化工股份有限公司 A kind of tight sandstone reservoir multi-parameter Permeability Prediction method
CN107917865B (en) * 2016-10-11 2020-01-31 中国石油化工股份有限公司 compact sandstone reservoir multi-parameter permeability prediction method
CN106772673A (en) * 2016-11-29 2017-05-31 西南石油大学 A kind of shale reservoir crustal stress predicts modeling method
US11175207B2 (en) 2016-12-20 2021-11-16 Petrochina Company Limited Method for obtaining conversion relationship between dynamic and static elastic parameters
CN107036907B (en) * 2017-04-07 2019-10-22 北京科技大学 A kind of fissure rock intensive analysis prediction technique
CN107036907A (en) * 2017-04-07 2017-08-11 北京科技大学 A kind of fissure rock intensive analysis Forecasting Methodology
CN108195669A (en) * 2017-10-09 2018-06-22 中国石油化工股份有限公司 The correction of oil reservoir Under Confining Pressure static mechanical parameter and Forecasting Methodology
CN108195669B (en) * 2017-10-09 2020-04-17 中国石油化工股份有限公司 Method for correcting and predicting static mechanical parameters of rock under oil reservoir confining pressure condition
CN107748205A (en) * 2017-12-05 2018-03-02 中国空气动力研究与发展中心计算空气动力研究所 A kind of elastic constant measurement method varied with temperature
CN108519317A (en) * 2018-04-24 2018-09-11 钦州学院 Rock stress-seepage coupling test device under uniaxial direct tensile load
CN108519317B (en) * 2018-04-24 2023-10-31 北部湾大学 Rock stress-seepage coupling test device under direct tensile load
CN108732048B (en) * 2018-05-31 2020-06-05 长安大学 Method for determining elastic-plastic yield point stress of graded broken stone repeated load test
CN108732048A (en) * 2018-05-31 2018-11-02 长安大学 A kind of elastoplasticity yield point stress of graded broken stone repeated-load test determines method
CN108843303B (en) * 2018-07-19 2021-08-10 西南石油大学 Oil-water well casing damage prediction method based on mudstone creep model
CN108843303A (en) * 2018-07-19 2018-11-20 西南石油大学 A kind of casing damage in oil-water well prediction technique based on mud stone creep model
CN109033672B (en) * 2018-08-09 2022-01-04 中国石油天然气股份有限公司 Dynamic crack determination method and device for displacement simulation mesoporous throat network model
CN109033672A (en) * 2018-08-09 2018-12-18 中国石油天然气股份有限公司 The dynamic crack determination method and device of pore constriction network model in displacement simulation
CN112673147A (en) * 2018-09-13 2021-04-16 沙特阿拉伯石油公司 Downhole tool for fracturing a hydrocarbon containing formation
CN109447944A (en) * 2018-09-21 2019-03-08 中国石油天然气股份有限公司 The Lithofacies Identification method and system of carbonate rock
CN109447944B (en) * 2018-09-21 2020-08-11 中国石油天然气股份有限公司 Lithofacies identification method and system for carbonate rock
CN109828031A (en) * 2019-02-15 2019-05-31 西南石油大学 Rock brittleness evaluation method and device
CN110006568A (en) * 2019-04-12 2019-07-12 中国地质科学院地质力学研究所 It is a kind of to obtain the method for three-dimensional ground stress using core and its obtain system
CN110006568B (en) * 2019-04-12 2020-01-07 中国地质科学院地质力学研究所 Method and system for acquiring three-dimensional ground stress by using rock core
CN110133724B (en) * 2019-05-06 2021-01-29 中国石油天然气集团有限公司 Method for identifying gas layer by applying array acoustic logging data
CN110133724A (en) * 2019-05-06 2019-08-16 中国石油天然气集团有限公司 A kind of application array sonic log data progress gas-bearing formation knowledge method for distinguishing
CN112001055B (en) * 2019-11-07 2024-04-09 中海石油(中国)有限公司 Low-amplitude thin oil reservoir water content prediction method based on microstructure
CN112001055A (en) * 2019-11-07 2020-11-27 中海石油(中国)有限公司 Low-amplitude thin oil reservoir water content prediction method based on microstructure
CN110836929B (en) * 2019-11-26 2021-01-19 中国科学院武汉岩土力学研究所 Jointed rock mass vibration attenuation physical simulation test method and device
CN110836929A (en) * 2019-11-26 2020-02-25 中国科学院武汉岩土力学研究所 Jointed rock mass vibration attenuation physical simulation test method and device
CN111089904B (en) * 2019-12-20 2021-05-28 西南石油大学 Indoor measurement wave velocity confining pressure and dispersion correction method considering formation characteristics
CN111089904A (en) * 2019-12-20 2020-05-01 西南石油大学 Indoor measurement wave velocity confining pressure and dispersion correction method considering formation characteristics
CN111175128A (en) * 2020-01-13 2020-05-19 西南石油大学 Compact sandstone tensile strength prediction method based on spectral analysis
CN111208565B (en) * 2020-03-04 2021-08-27 中国石油大学(北京) KT model-based hole seam parameter inversion method and device and storage medium
CN111208565A (en) * 2020-03-04 2020-05-29 中国石油大学(北京) KT model-based hole seam parameter inversion method and device and storage medium
CN112324429A (en) * 2020-08-12 2021-02-05 中国石油天然气股份有限公司 Method for determining reasonable production pressure difference of water-gas reservoir
CN112324429B (en) * 2020-08-12 2023-08-22 中国石油天然气股份有限公司 Method for determining reasonable production pressure difference of water-gas reservoir
CN113361131A (en) * 2021-06-25 2021-09-07 中国矿业大学 Static mechanical parameter calculation method for primary structure coal
CN114113083A (en) * 2021-11-10 2022-03-01 中国石油大学(北京) Method, device, medium, apparatus and program product for evaluating heterogeneity of conglomerate
CN114113083B (en) * 2021-11-10 2024-01-19 中国石油大学(北京) Method, device, medium, device, and program product for evaluating conglomerate heterogeneity
CN114526066A (en) * 2022-01-24 2022-05-24 西南石油大学 Method for compacting and recovering continental facies shale oil reservoir
CN114526066B (en) * 2022-01-24 2024-06-11 西南石油大学 Method for recovering compaction of land shale oil reservoir
CN115266514A (en) * 2022-05-11 2022-11-01 中国石油大学(华东) Dynamic evaluation device and method for rock mechanical parameters in high-pressure fluid injection process
CN116066092A (en) * 2023-02-03 2023-05-05 西南石油大学 Shale formation collapse pressure increment prediction method based on hardness experiment
CN116066092B (en) * 2023-02-03 2024-04-26 西南石油大学 Shale formation collapse pressure increment prediction method based on hardness experiment

Also Published As

Publication number Publication date
CN103257081B (en) 2015-09-02

Similar Documents

Publication Publication Date Title
CN103257081B (en) A kind of method that hydrocarbon-bearing pool rock mass mechanics ground in-situ model recovers and device
CN103278614B (en) Method and device for correcting dynamic and static rock mechanical parameters
CN103267678B (en) A kind of method that hydrocarbon-bearing pool rock mass mechanics ground in-situ model recovers and device
CN103256046B (en) Unconventionaloil pool hides method and the device that horizontal well stitches the simulation of long fracturing parameter entirely
Liang et al. Wellbore stability model for shale gas reservoir considering the coupling of multi-weakness planes and porous flow
CN103278389A (en) Method for synchronous measurements on dynamic and static elastic parameters of rocks
WO2017204689A1 (en) Method for determining the mechanical properties of reservoir rock
Walker et al. Fault zone architecture and fluid flow in interlayered basaltic volcaniclastic-crystalline sequences
Nauroy Geomechanics applied to the petroleum industry
Damani et al. Mapping of hydraulic fractures under triaxial stress conditions in laboratory experiments using acoustic emissions
Hamza et al. Characterization of anisotropic elastic moduli and stress for unconventional reservoirs using laboratory static and dynamic geomechanical data
Ng et al. Small strain stiffness of natural granitic saprolite in Hong Kong
Yue et al. Investigation of acoustic emission response and fracture morphology of rock hydraulic fracturing under true triaxial stress
Du et al. New expression of the changing stress field in low-permeability reservoir and its application in secondary exploitation
Yang et al. Study on evaluation method of fracture forming ability of shale oil reservoirs in Fengcheng Formation, Mahu sag
Shi et al. Characterization of hydraulic fracture configuration based on complex in situ stress field of a tight oil reservoir in Junggar Basin, Northwest China
Adero Experimental investigations of mechanical anisotropy of Freiberg gneiss: implications for hydraulic stimulation
Abbas et al. Practical approach for sand-production prediction during production
Han et al. Stress-Dependent flow, mechanical, and acoustic properties for an unconventional oil reservoir rock
Halleck et al. Reduction of jet perforator penetration in rock under stress
Al-Malikee et al. Indirect prediction of rock elasticity and compressibility strength using well log data at selected sites within Rumaila Oilfield, Southern Iraq
Zhang et al. Geomechanical evaluation enabled successful stimulation of a HPHT tight gas reservoir in western China
Munir Development of correlation between rock classification system and modulus of deformation
Hill et al. Techniques for determining subsurface stress direction and assessing hydraulic fracture azimuth
HASAN et al. Correlation of Shear Wave Velocity with SPT-N for a Tower-Building Site at Erbil City

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150902

Termination date: 20190428

CF01 Termination of patent right due to non-payment of annual fee