CN104484573B - A kind of method of determination stratum stiffness coefficient - Google Patents
A kind of method of determination stratum stiffness coefficient Download PDFInfo
- Publication number
- CN104484573B CN104484573B CN201410842028.XA CN201410842028A CN104484573B CN 104484573 B CN104484573 B CN 104484573B CN 201410842028 A CN201410842028 A CN 201410842028A CN 104484573 B CN104484573 B CN 104484573B
- Authority
- CN
- China
- Prior art keywords
- ground
- stratum
- cores
- core sample
- coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Geophysics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
- Mining & Mineral Resources (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The method that the present invention proposes a kind of determination stratum stiffness coefficient, including:Acoustic anisotropy measurement is carried out to core sample and bulk density is calculated, stiffness coefficient and compressional wave anisotropy coefficient is obtained;X diffractometries are carried out to core sample, clay content is obtained;Set up the relation of compressional wave anisotropy coefficient and clay content;Set up C11 coresWith C33 cores、C44 coresAnd C66 coresRelation;Using the Conventional Logs of target well destination layer position, its clay content is calculated;Using the P- and S-wave velocity and bulk density log of target well destination layer position, the stiffness coefficient C on stratum is calculated33 groundAnd C44 ground;According to compressional wave anisotropy coefficient and relation, the clay content on stratum and the stiffness coefficient C on stratum of its clay content33 ground, calculate the stiffness coefficient C on stratum11 ground;According to C11 coresWith C33 cores、C44 coresAnd C66 coresRelation, C33 ground、C44 groundAnd C11 ground, calculate the stiffness coefficient C on stratum66 ground.The present invention calculates stiffness coefficient without using the horizontal shear wave velocity being finally inversed by from the Stoneley wave of extraction, avoids horizontal shear wave extraction process, makes calculating simple and effective.
Description
Technical field
The present invention relates to petroleum exploration logging technical field, more particularly to a kind of method of determination stratum stiffness coefficient.
Background technology
Stiffness coefficient is to ensconce rock mechanics logging evaluation mistake using fine and close oil gas, shale oil gas as the unconventionaloil pool of representative
Key parameter necessary in journey.It specifically includes C33、C44、C11、C13And C66Etc. parameter.Wherein C66And C11The two are rigidly
Several well logging computational methods are a difficulties in untraditional reservoir rock mechanics logging evaluation field in recent years.To understand
Certainly this problem, the thinking that forefathers use is:Inverting level is horizontal in the Stoneley wave extracted first from Array Sonic Logging Waveformss
Wave velocity, it then combined to calculate C with bulk density log66;Five rigidity systems are further utilized on this basis
Experimental relationship between number, finally realizes C11Well logging Deng stiffness coefficient is calculated.There is obvious limitation in this way, main
Reason is wanted to be in fast stratum, it is poor to the sensitiveness of horizontal shear wave degree in Stoneley wave, extract error greatly, it is with a low credibility.
Article " Determining formation shear-wave transverse isotropy from borehole
Pointed out in stoneley-wave measurements, Xiaoming Tang, Geophysics.Vol.68, No.1,2003 ":
The way of the horizontal shear wave velocity of extraction is typically only applicable to slowly layer from well logging Stoneley wave.And extract water from Stoneley wave
The process of flat shear wave velocity is also sufficiently complex, realizes that difficulty is big for those skilled in the art.
The content of the invention
The embodiments of the invention provide a kind of stiffness coefficient computational methods, extracted without using from Array Sonic Logging Waveformss
Stoneley wave in the horizontal shear wave velocity that is finally inversed by calculate stiffness coefficient, so as to avoid horizontal shear wave extraction process, make
Simple and effective must be calculated;What is more important this method is applied also for soon in addition to being calculated suitable for the stiffness coefficient of slowly layer
The stiffness coefficient of layer is calculated, with the broader scope of application.This method includes:
Acoustic anisotropy measurement is carried out to core sample and bulk density is calculated, the stiffness coefficient of acquisition core sample,
Compressional wave anisotropy coefficient and shear wave anisotropy coefficient;
X diffractometries are carried out to core sample, its clay content data is obtained;
The compressional wave anisotropy coefficient of core sample and the experimental relationship of its clay content are set up, the shear wave of core sample is each
The experimental relationship of anisotropy coefficient and its clay content;
According to the stiffness coefficient of the core sample of acquisition, stiffness coefficient C is set up11 coresWith C33 cores、C44 coresAnd C66 coresPass
System;
Utilize the Conventional Logs of target well destination layer position, the clay content on continuous depth calculation destination layer position stratum;
Utilize velocity of longitudinal wave, shear wave velocity and the bulk density log of target well destination layer position, continuous depth calculation
The stiffness coefficient C on destination layer position stratum33 groundAnd C44 ground;
According to relation, the clay content on stratum and the stratum of the compressional wave anisotropy coefficient of core sample He its clay content
Stiffness coefficient C33 ground, the stiffness coefficient C on continuous depth calculation stratum11 ground;According to the stiffness coefficient C of core sample11 coresWith
C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient C on stratum33 ground、C44 groundAnd C11 ground, the stiffness coefficient on continuous depth calculation stratum
C66 ground;
Or, the relation of the shear wave anisotropy coefficient according to core sample and its clay content, the clay content on stratum and
The stiffness coefficient C on stratum44 ground, the stiffness coefficient C on continuous depth calculation stratum66 ground;According to the stiffness coefficient C of core sample11 coresWith
C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient C on stratum33 ground、C44 groundAnd C66 ground, the stiffness coefficient on continuous depth calculation stratum
C11 ground。
In one embodiment, the bulk density calculation formula of the core sample is:
Wherein, ρCoreFor the bulk density of core sample;
WSat coresFor the weight of core sample;
VBulk coresFor the volume of core sample.
In one embodiment, the calculation formula of the stiffness coefficient of the core sample is:
Wherein, Vp0For the velocity of longitudinal wave measured by one group of lateral P wave emission of core sample with receiving transducer, the group is visited
Head is vertical with bed plane;
Vp90The velocity of longitudinal wave measured for one group of P wave emission at the top and bottom of core sample and receiving transducer;
Vs190The horizontal stroke measured for the shear wave transmitting of one group of cross polarization at the top and bottom of core sample with receiving transducer
Wave velocity;
Vs290What shear wave transmitting and receiving transducer for another group of cross polarization at the top and bottom of core sample were measured
Velocity of longitudinal wave;
C33 cores、C11 cores、C44 cores、C66 coresIt is the stiffness coefficient of core sample.
In one embodiment, the calculation formula of the compressional wave anisotropy coefficient of the core sample is:
Wherein, εCoreFor the compressional wave anisotropy coefficient of core sample;
The calculation formula of the shear wave anisotropy coefficient of the core sample is:
Wherein, γCoreFor the shear wave anisotropy coefficient of core sample.
In one embodiment, the experimental relationship table of the compressional wave anisotropy coefficient of the core sample and its clay content
It is up to formula:
Wherein, VCL coresFor the clay content of core sample;k1、n1For variable parameter;
The shear wave anisotropy coefficient of the core sample is with the experimental relationship expression formula of its clay content:
Wherein, k2、n2For variable parameter.
In one embodiment, the stiffness coefficient C of the core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelationship expression
Formula is:
Wherein, m, r are variable parameter.
In one embodiment, the Conventional Logs of target well destination layer position include natural gamma curve, sound wave,
Neutron, density, the uranium in natural gamma energy spectrum data, thorium, potassium curve and deep and shallow resistivity curve.
In one embodiment, the stiffness coefficient C on the stratum33 groundCalculation formula be:
Wherein, ρGroundFor the bulk density on stratum, VpFor the velocity of longitudinal wave measured;
The stiffness coefficient C on the stratum44 groundCalculation formula be:
Wherein, Vs is the shear wave velocity measured.
In one embodiment, the relation of the compressional wave anisotropy coefficient according to core sample and its clay content,
The clay content on stratum and the stiffness coefficient C on stratum33 ground, the stiffness coefficient C on continuous depth calculation stratum11 ground, including by following public affairs
Formula is calculated:
The stiffness coefficient C according to core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient on stratum
C33 ground、C44 groundAnd C11 ground, the stiffness coefficient C on continuous depth calculation stratum66 ground, including be calculated as follows:
Wherein, VCLFor the clay content on stratum.
In one embodiment, the relation of the shear wave anisotropy coefficient according to core sample and its clay content,
The clay content on stratum and the stiffness coefficient C on stratum44 ground, the stiffness coefficient C on continuous depth calculation stratum66 ground, including by following public affairs
Formula is calculated:
The stiffness coefficient C according to core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient on stratum
C33 ground、C44 groundAnd C66 ground, the stiffness coefficient C on continuous depth calculation stratum11 ground, including be calculated as follows:
In embodiments of the present invention, without using the water being finally inversed by from the Stoneley wave that Array Sonic Logging Waveformss are extracted
Flat shear wave velocity calculates stiffness coefficient, so as to avoid horizontal shear wave extraction process so that calculate simple and effective;It is even more important
Be this method except calculating outer suitable for the stiffness coefficient of slowly layer, the stiffness coefficient calculating on fast stratum is applied also for, with more
The wide scope of application.
Brief description of the drawings
Accompanying drawing described herein is used for providing a further understanding of the present invention, constitutes the part of the application, not
Constitute limitation of the invention.In the accompanying drawings:
Fig. 1 is a kind of flow chart of the method for determination stratum stiffness coefficient in the embodiment of the present invention;
Fig. 2 is a kind of schematic diagram of single core method of testing progress acoustic anisotropy measurement in the embodiment of the present invention;
Fig. 3 is clay content and stiffness coefficient C in the embodiment of the present invention33And C44Schematic diagram of calculation result;
Fig. 4, Fig. 7 are stiffness coefficient C in the embodiment of the present invention11Schematic diagram of calculation result;
Fig. 5, Fig. 6 are stiffness coefficient C in the embodiment of the present invention66Schematic diagram of calculation result;
Fig. 8 is the Comprehensive Correlation schematic diagram of distinct methods calculating stiffness coefficient result in the embodiment of the present invention.
Embodiment
It is right with reference to embodiment and accompanying drawing for the object, technical solutions and advantages of the present invention are more clearly understood
The present invention is described in further details.Here, the exemplary embodiment of the present invention and its illustrating to be used to explain the present invention, but simultaneously
It is not as a limitation of the invention.
The method of existing calculating stratum stiffness coefficient is needed using the horizontal shear wave velocity obtained from Stoneley wave inverting,
But this method has significant limitations, slowly layer is typically only applicable to, and extract from Stoneley wave the process of horizontal shear wave velocity
Sufficiently complex, so that the process for calculating stratum stiffness coefficient also becomes complicated, general well log interpretation personnel are difficult to grasp.Hair
A person of good sense proposes a kind of new method, this method when calculating stratum stiffness coefficient without using horizontal shear wave velocity, and then can be with
Horizontal shear wave extraction process is avoided, and this method has the broader scope of application, is applicable not only to slowly layer, is also equally applicable to
Fast stratum, can solve the problem that the above-mentioned problems in the prior art.It is specifically described below.
Fig. 1 is a kind of flow chart of the method for determination stratum stiffness coefficient in the embodiment of the present invention, as shown in figure 1, the party
Method includes:
Step 101:Acoustic anisotropy measurement is carried out to core sample and bulk density is calculated, the firm of core sample is obtained
Property coefficient, compressional wave anisotropy coefficient and shear wave anisotropy coefficient;
Step 102:X diffractometries are carried out to core sample, its clay content data is obtained;
Step 103:Set up the compressional wave anisotropy coefficient of core sample and the experimental relationship of its clay content, core sample
Shear wave anisotropy coefficient and its clay content experimental relationship;
Step 104:According to the stiffness coefficient of the core sample of acquisition, stiffness coefficient C is set up11 coresWith C33 cores、C44 coresWith
C66 coresRelation;
Step 105:Using the Conventional Logs of target well destination layer position, continuous depth calculation destination layer position stratum it is viscous
Native content;
Step 106:Using velocity of longitudinal wave, shear wave velocity and the bulk density log of target well destination layer position, continuously
The stiffness coefficient C on depth calculation destination layer position stratum33 groundAnd C44 ground;
Step 107:Contained according to the clay on the relation of the compressional wave anisotropy coefficient of core sample and its clay content, stratum
Amount and the stiffness coefficient C on stratum33 ground, the stiffness coefficient C on continuous depth calculation stratum11 ground;According to the stiffness coefficient of core sample
C11 coresWith C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient C on stratum33 ground、C44 groundAnd C11 ground, continuous depth calculation stratum it is firm
Property coefficient C66 ground;
Or, the relation of the shear wave anisotropy coefficient according to core sample and its clay content, the clay content on stratum and
The stiffness coefficient C on stratum44 ground, the stiffness coefficient C on continuous depth calculation stratum66 ground;According to the stiffness coefficient of water saturation core sample
C11 coresWith C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient C on stratum33 ground、C44 groundAnd C66 ground, continuous depth calculation stratum it is firm
Property coefficient C11 ground。
When it is implemented, the preparation method of core sample is:From the full diameter core that emphasis prospect pit drilling extracting core is obtained
Coring is carried out along the direction parallel to stratum stratification, the horizontal plug sample core of acquisition is carried out obtaining experiment institute after water saturation
Core sample (could also say that water saturation core sample).
To core sample using single core method of testing progress acoustic anisotropy measurement, principle as shown in Fig. 2 described list
Core method of testing may be referred to file " Seismic anisotropy in sedimentary rocks, Part 1:A
Single-plug laboratory method ", Zhijing Wang, Geophysics.Vol.67, No.5,2002.
In Fig. 2, there are 3 groups of probes at the top and bottom of core sample, wherein 1 group is P wave emission and receiving transducer, for surveying
Measure velocity of longitudinal wave Vp90, the stiffness coefficient C for calculating core sample11 cores;The other 2 groups shear waves for cross polarization are launched with connecing
Probe is received, for measuring shear wave velocity Vs190And Vs290, the stiffness coefficient C for calculating core sample44 coresAnd C66 cores;In core
Sample laterally has 2 groups of P wave emissions and receiving transducer, and 1 group vertical with bed plane, measurement velocity of longitudinal wave Vp0, for calculating core sample
The stiffness coefficient C of product33 cores;1 group with bed plane into 45 degree angles, measure velocity of longitudinal wave Vp45, the rigidity system for calculating core sample
Number C13 cores。
It is also desirable to calculate the bulk density of core sample, its calculation formula is:
In formula, ρCoreFor the bulk density of core sample;WSat coresFor the weight of core sample;VBulk coresFor core sample
Volume.
The velocity of wave and bulk density result of calculation obtained according to acoustic anisotropy measurement, calculates the rigidity system of core sample
Number, specific formula is as follows:
Compressional wave anisotropy coefficient, shear wave anisotropy coefficient refer to Thomsen in 1986 propose be used for describe to indulge
The parameter of wave velocity anisotropic character, shear wave velocity anisotropic character, specifically refers to article " Weak elastic
anisotropy:Geophysics.Vol.51,1986 ".
The expression formula of the compressional wave anisotropy coefficient of core sample is:
In formula, εCoreFor the compressional wave anisotropy coefficient of core sample.
The expression formula of the shear wave anisotropy coefficient of core sample is:
In formula, γCoreFor the shear wave anisotropy coefficient of core sample.
Method used in X diffractometries is carried out to core and refers to professional standard SY/T 5163-2010.It is real by X diffraction
Test, obtain the clay content information of core sample.
According to the P-wave And S anisotropy coefficient and its clay content data of core sample, two groups of experiment passes are set up respectively
System:First group of clay content using core sample is independent variable, and the compressional wave anisotropy coefficient using core sample is dependent variable, table
It is up to formula;
In formula, VCL coresFor the clay content of core sample;k1、n1For variable parameter.
Second group of clay content using core sample is independent variable, using the shear wave anisotropy coefficient of core sample as because becoming
Measure, expression formula is;
In formula, k2、n2For variable parameter.
Set up the stiffness coefficient C of core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelation, be specially:With C33 cores、
C44 coresAnd C66 coresCombination parameter based on these three parametersFor independent variable, with C11 coresFor dependent variable, set up
Therebetween relation, specific expression formula is:
Wherein, m, r are variable parameter.
When calculating the clay content on stratum, it is possible to use the natural gamma ray logging curve of target well destination layer position is calculated,
Formula is as follows:
In formula, VCLFor the clay content on the stratum of current depth point;
GR is the natural gamma ray logging curve values of current depth point;
GRminFor the GR characteristic value of destination layer position clean sandstone section;
GRmaxFor the GR characteristic value of destination layer position pure shale section.
Can also using optimization algorithm using target well destination layer position more well-log informations (GR, sound wave,
Neutron, density, the uranium in natural gamma energy spectrum data, thorium, potassium curve and deep and shallow resistivity curve) cutd open to calculate polymineralic rock
Face, so as to obtain the clay content on stratum, or utilizes one or more of surveys with target well destination layer position using other method
Well data calculates the clay content on stratum.
The stiffness coefficient C on stratum is calculated using the velocity of longitudinal wave and bulk density log of target well destination layer position33 ground's
Formula is:
In formula, ρGroundFor the bulk density on stratum, VpFor the velocity of longitudinal wave on stratum;
The stiffness coefficient C on stratum is calculated using the shear wave velocity and bulk density log of target well destination layer position44 ground's
Formula is:
In formula, Vs is the shear wave velocity on stratum.
Calculate the stiffness coefficient C on stratum11 groundAnd C66 groundThere are two methods:The first is first to calculate C11 ground, C is calculated afterwards66 ground.Root
According to the compressional wave anisotropy coefficient and relation, the clay content on stratum and the stiffness coefficient on stratum of its clay content of core sample
C33 ground, the stiffness coefficient C on continuous depth calculation stratum11 ground;According to the stiffness coefficient C of core sample11 coresWith C33 cores、C44 coresWith
C66 coresRelation, the stiffness coefficient C on stratum33 ground、C44 groundAnd C11 ground, the stiffness coefficient C on continuous depth calculation stratum66 ground。
Specifically, according to the relational expression (8) of compressional wave anisotropy coefficient and its clay content, compressional wave anisotropy coefficient
Calculation formula (6) and the stiffness coefficient C on stratum33 groundCalculation formula (12), finally give the stiffness coefficient C on stratum11 groundCalculating
Formula is as follows:
According to the stiffness coefficient C of core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelational expression (10), stratum rigidity
Coefficient C33 groundCalculation formula (12), C44 groundCalculation formula (13) and C11 groundCalculation formula (14), finally give the firm of stratum
Property coefficient C66 groundCalculation formula it is as follows:
Wherein, VCLFor the clay content on stratum.
Second is first to calculate C66 ground, C is calculated afterwards11 ground.Specially:According to the shear wave anisotropy coefficient of core sample and
The stiffness coefficient C of the relation of its clay content, the clay content on stratum and stratum44 ground, the stiffness coefficient on continuous depth calculation stratum
C66 ground;According to the stiffness coefficient C of water saturation core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient on stratum
C33 ground、C44 groundAnd C66 ground, the stiffness coefficient C on continuous depth calculation stratum11 ground。
Specifically, according to the relational expression (9) of shear wave anisotropy coefficient and its clay content, shear wave anisotropy coefficient
Calculation formula (7) and stratum stiffness coefficient C44 groundCalculation formula (13), finally give the stiffness coefficient C on stratum66 groundCalculating it is public
Formula is as follows:
According to the stiffness coefficient C of core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelational expression (10), stratum rigidity
Coefficient C33 groundCalculation formula (12) and stratum stiffness coefficient C66 groundCalculation formula (16), finally give the stiffness coefficient on stratum
C11 groundCalculation formula it is as follows:
This hair is illustrated by taking stratum (slowly layer) near fine and close oil-gas reservoir a bite emphasis prospect pit hydrocarbon source rock interval as an example
Bright technical scheme.
The first step, plunger sample core is drilled through from the full diameter drilling extracting core of the emphasis well along the direction parallel to stratum stratification
Sample, single core acoustic anisotropy measurement and bulk density are carried out to the horizontal plug sample core sample of acquisition after water saturation
Calculate, the stiffness coefficient and P-wave And S anisotropy coefficient for obtaining horizontal plug sample core sample are as shown in table 1.
The stiffness coefficient and P-wave And S anisotropy coefficient of the horizontal plug sample core sample of table 1
Second step, carries out X diffractometries to the horizontal plug sample core sample obtained by the first step, obtains every piece of horizontal columns
The clay content of sample core sample is filled in, as shown in table 2.
The clay content of the horizontal plug sample core sample of table 2
3rd step, the relational expression of compressional wave anisotropy coefficient and clay content is set up according to Tables 1 and 2, is specially:
The relational expression of shear wave anisotropy coefficient and clay content is set up according to Tables 1 and 2, is specially:
4th step, measures and calculates according to the first step stiffness coefficient of obtained horizontal plug sample core sample, set up
C11 coresWith C33 cores、C44 coresAnd C66 coresExperimental relationship, it is as follows:
I.e. in this example, slope m=1, intercept r=0.
5th step, the clay for calculating stratum using the natural gamma ray logging curve and formula (11) of target well destination layer position contains
Amount;
6th step, using P- and S-wave velocity and bulk density curve according to formula (12) and the continuous depth gauge of formula (13)
Calculate the stiffness coefficient C on stratum33 groundAnd C44 ground, as a result as shown in Figure 3.Wherein, the GR in first is natural gamma curve value, single
Position is gAPI;Second is depth track, and the 3rd road shows interval transit time, and wherein DTSM is shear wave slowness, and DTCO is compressional wave time difference,
Unit is us/m;V in 4th roadCLFor the clay content on stratum;5th and the 6th road be shown respectively by interval transit time and
Density curve calculates the stiffness coefficient C on obtained stratum33 groundAnd C44 ground, unit is Gpa.
7th step, C is first calculated using first method11 groundC is calculated again66 ground:The clay content on the stratum of acquisition is substituted into public
Formula (18), obtains the compressional wave anisotropy coefficient on stratum, then according to formula (14) and the stiffness coefficient C on stratum33 ground, final
To the stiffness coefficient C on stratum11 ground, it is designated as C11a, result of calculation is as shown in Figure 4.
According to formula (15) and the stiffness coefficient C on stratum33 ground、C44 groundAnd C11 ground, finally give the stiffness coefficient C on stratum66 ground,
It is designated as C66a, result of calculation is as shown in Figure 5.
Or, C is first calculated using second method66 groundC is calculated again11 ground:The clay content on the stratum of acquisition is substituted into formula
(19) the shear wave anisotropy coefficient on stratum, is obtained, then according to formula (16) and the stiffness coefficient C on stratum44 ground, finally give
The stiffness coefficient C on stratum66 ground, it is designated as C66b, as a result as shown in Figure 6.
According to formula (17) and the stiffness coefficient C on stratum33 ground、C44 groundAnd C66 ground, finally give the stiffness coefficient C on stratum11 ground,
It is designated as C11b, result of calculation is as shown in Figure 7.
In interval of interest, using the method for calculating stiffness coefficient in the prior art, i.e., realized using horizontal shear wave velocity
Stiffness coefficient C66Well logging characterize, formula is:
In formula, ρ is bulk density log;VshFor the horizontal shear wave velocity extracted from Stoneley wave well-log information.
Compare the C obtained based on horizontal shear wave velocity and experimental relationship11And C66The C obtained with first method11a、C66a
And the C that second method is obtained11b、C66b, comparative result is as shown in Figure 8.As seen from Figure 8, two provided using the present invention
The C that the method for kind is obtained11aWith C11bIt is very identical, C66aWith C66bAlso it is very identical.By Fig. 8 it can also be seen that using from stone profit
The stiffness coefficient C that the horizontal shear wave velocity extracted in ripple well-log information is obtained66With C66a、C66bAlso it coincide substantially, it is horizontal based on level
The C that wave velocity and experimental relationship are obtained11With C11a、C11bAlso it is consistent substantially, this illustrates in the case of unhorizontal shear wave velocity,
The stiffness coefficient C that can equally obtain being consistent with actual conditions using the present invention66And C11, and the present invention is due to avoiding level
The extraction process of shear wave velocity, therefore than the method for forefathers there is the broader scope of application (forefathers' method can be only applied to slowly
Layer).
In summary, calculated using the method for the invention provided during stiffness coefficient without using the inverting from Stoneley wave
The horizontal shear wave velocity gone out, and then avoid horizontal shear wave extraction process so that calculate simple and effective;The method that the present invention is provided
In addition to being calculated suitable for the stiffness coefficient of slowly layer, the stiffness coefficient for applying also for fast stratum is calculated, with broader applicable model
Enclose.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention, for the skill of this area
For art personnel, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles of the invention, made
Any modification, equivalent substitution and improvements etc., should be included in the scope of the protection.
Claims (9)
1. a kind of method of determination stratum stiffness coefficient, it is characterised in that including:
Acoustic anisotropy measurement is carried out to core sample and bulk density is calculated, stiffness coefficient, the compressional wave of core sample is obtained
Anisotropy coefficient and shear wave anisotropy coefficient;
X diffractometries are carried out to core sample, its clay content data is obtained;
The compressional wave anisotropy coefficient of core sample and the experimental relationship of its clay content are set up, the shear wave of core sample is each to different
The experimental relationship of property coefficient and its clay content;
According to the stiffness coefficient of the core sample of acquisition, stiffness coefficient C is set up11 coresWith C33 cores、C44 coresAnd C66 coresRelation;
Utilize the Conventional Logs of target well destination layer position, the clay content on continuous depth calculation destination layer position stratum;
Utilize velocity of longitudinal wave, shear wave velocity and the bulk density log of target well destination layer position, continuous depth calculation target
The stiffness coefficient C on layer position stratum33 groundAnd C44 ground;
According to the firm of the relation of the compressional wave anisotropy coefficient of core sample and its clay content, the clay content on stratum and stratum
Property coefficient C33 ground, the stiffness coefficient C on continuous depth calculation stratum11 ground;According to the stiffness coefficient C of core sample11 coresWith C33 cores、
C44 coresAnd C66 coresRelation, the stiffness coefficient C on stratum33 ground、C44 groundAnd C11 ground, the stiffness coefficient C on continuous depth calculation stratum66 ground;
Or, relation, the clay content on stratum and the stratum of the shear wave anisotropy coefficient according to core sample and its clay content
Stiffness coefficient C44 ground, the stiffness coefficient C on continuous depth calculation stratum66 ground;According to the stiffness coefficient C of water saturation core sample11 cores
With C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient C on stratum33 ground、C44 groundAnd C66 ground, the rigidity system on continuous depth calculation stratum
Number C11 ground;
The compressional wave anisotropy coefficient of the core sample is with the experimental relationship expression formula of its clay content:
Wherein, VCL coresFor the clay content of core sample;k1、n1For variable parameter;
The shear wave anisotropy coefficient of the core sample is with the experimental relationship expression formula of its clay content:
Wherein, k2、n2For variable parameter.
2. according to the method described in claim 1, it is characterised in that the bulk density calculation formula of the core sample is:
Wherein, ρCoreFor the bulk density of core sample;
WSat coresFor the weight of core sample;
VBulk coresFor the volume of core sample.
3. method according to claim 2, it is characterised in that the calculation formula of the stiffness coefficient of the core sample is:
Wherein, Vp0For the velocity of longitudinal wave measured by one group of lateral P wave emission of core sample with receiving transducer, group probe with
Bed plane is vertical;
Vp90The velocity of longitudinal wave measured for one group of P wave emission at the top and bottom of core sample and receiving transducer;
Vs190The shear wave speed measured for the shear wave transmitting of one group of cross polarization at the top and bottom of core sample with receiving transducer
Degree;
Vs290The compressional wave speed measured for the shear wave transmitting of another group of cross polarization at the top and bottom of core sample with receiving transducer
Degree;
C33 cores、C11 cores、C44 cores、C66 coresIt is the stiffness coefficient of core sample.
4. method according to claim 3, it is characterised in that the calculating of the compressional wave anisotropy coefficient of the core sample
Formula is:
Wherein, εCoreFor the compressional wave anisotropy coefficient of core sample;
The calculation formula of the shear wave anisotropy coefficient of the core sample is:
Wherein, γCoreFor the shear wave anisotropy coefficient of core sample.
5. according to the method described in claim 1, it is characterised in that the stiffness coefficient C of the core sample11 coresWith C33 cores、
C44 coresAnd C66 coresRelational expression be:
Wherein, m, r are variable parameter.
6. method according to claim 5, it is characterised in that the Conventional Logs of the target well destination layer position include
Natural gamma curve, sound wave, neutron, density, the uranium in natural gamma energy spectrum data, thorium, potassium curve and deep and shallow resistivity curve.
7. method according to claim 6, it is characterised in that the stiffness coefficient C on the stratum33 groundCalculation formula be:
Wherein, ρGroundFor the bulk density on stratum, VpFor the velocity of longitudinal wave measured;
The stiffness coefficient C on the stratum44 groundCalculation formula be:
Wherein, Vs is the shear wave velocity measured.
8. method according to claim 7, it is characterised in that the compressional wave anisotropy coefficient according to core sample and
The stiffness coefficient C of the relation of its clay content, the clay content on stratum and stratum33 ground, the stiffness coefficient on continuous depth calculation stratum
C11 ground, including be calculated as follows:
The stiffness coefficient C according to core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient C on stratum33 ground、
C44 groundAnd C11 ground, the stiffness coefficient C on continuous depth calculation stratum66 ground, including be calculated as follows:
Wherein, VCLFor the clay content on stratum;
ε is the compressional wave anisotropy coefficient on stratum, with the stiffness coefficient C on stratum33 groundAnd C11 groundTo represent:
9. method according to claim 8, it is characterised in that the shear wave anisotropy coefficient according to core sample and
The stiffness coefficient C of the relation of its clay content, the clay content on stratum and stratum44 ground, the stiffness coefficient on continuous depth calculation stratum
C66 ground, including be calculated as follows:
The stiffness coefficient C according to core sample11 coresWith C33 cores、C44 coresAnd C66 coresRelation, the stiffness coefficient C on stratum33 ground、
C44 groundAnd C66 ground, the stiffness coefficient C on continuous depth calculation stratum11 ground, including be calculated as follows:
Wherein, γ is the shear wave anisotropy coefficient on stratum, with the stiffness coefficient C on stratum44 groundAnd C66 groundTo represent:
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410842028.XA CN104484573B (en) | 2014-12-30 | 2014-12-30 | A kind of method of determination stratum stiffness coefficient |
AU2015275302A AU2015275302B2 (en) | 2014-12-30 | 2015-12-23 | Method and apparatus of determining stiffness coefficients of formation |
US14/984,013 US10031250B2 (en) | 2014-12-30 | 2015-12-30 | Method and apparatus of determining stiffness coefficients of formation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410842028.XA CN104484573B (en) | 2014-12-30 | 2014-12-30 | A kind of method of determination stratum stiffness coefficient |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104484573A CN104484573A (en) | 2015-04-01 |
CN104484573B true CN104484573B (en) | 2017-09-15 |
Family
ID=52759114
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410842028.XA Active CN104484573B (en) | 2014-12-30 | 2014-12-30 | A kind of method of determination stratum stiffness coefficient |
Country Status (3)
Country | Link |
---|---|
US (1) | US10031250B2 (en) |
CN (1) | CN104484573B (en) |
AU (1) | AU2015275302B2 (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104965227B (en) * | 2015-06-25 | 2018-04-03 | 中国石油天然气股份有限公司 | A kind of computational methods and device of compact reservoir well logging stiffness coefficient |
CN105093295B (en) * | 2015-06-26 | 2017-03-08 | 中国石油天然气股份有限公司 | A kind of continuous deep treatment method of horizontal shear wave sensitivity and device |
CN108843312B (en) * | 2018-06-04 | 2019-04-30 | 西南石油大学 | Oilfield reservoir in-layer heterogeneity integrated evaluating method |
CN109684764B (en) * | 2019-01-02 | 2023-06-02 | 中铁第六勘察设计院集团有限公司 | Method for calculating reasonable depth of pile foundation embedded in steel pipe column |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630014B (en) * | 2008-07-16 | 2011-10-05 | 中国石油天然气集团公司 | Method for imaging anisotropic medium through utilization of vertical seismic profile data |
CN102466818A (en) * | 2010-11-11 | 2012-05-23 | 中国石油集团东方地球物理勘探有限责任公司 | Method for imaging anisotropic medium by utilizing cross-well seismic data |
CN103323575A (en) * | 2013-05-20 | 2013-09-25 | 中国石油天然气股份有限公司 | Sedimentary rock reconstruction method and system |
CN102706965B (en) * | 2012-06-12 | 2014-08-27 | 西安邮电大学 | Novel accurate measuring method of anisotropic parameter of rock |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6351991B1 (en) * | 2000-06-05 | 2002-03-05 | Schlumberger Technology Corporation | Determining stress parameters of formations from multi-mode velocity data |
US6714873B2 (en) * | 2001-12-17 | 2004-03-30 | Schlumberger Technology Corporation | System and method for estimating subsurface principal stresses from seismic reflection data |
JP4250945B2 (en) | 2002-10-18 | 2009-04-08 | 株式会社Ihi | Cell foundation structure and cell foundation construction method |
US7751980B2 (en) * | 2006-12-22 | 2010-07-06 | Schlumberger Technology Corporation | Method and apparatus for evaluating elastic mechanical properties of a transversely isotropic formation |
US8175807B2 (en) * | 2008-02-20 | 2012-05-08 | Schlumberger Technology Corporation | Estimating horizontal stress from three-dimensional anisotropy |
US8184502B2 (en) * | 2008-02-28 | 2012-05-22 | Exxonmobil Upstream Research Company | Rock physics model for simulating seismic response in layered fractured rocks |
RU2602409C2 (en) * | 2011-07-27 | 2016-11-20 | Шлюмбергер Текнолоджи Б.В. | Inversion anisotropy multi-well system |
CN103628866B (en) | 2013-11-22 | 2017-02-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for obtaining stratum pressure coefficients of stratum |
CN103792581B (en) | 2013-12-31 | 2015-06-03 | 孙赞东 | Shale gas reservoir crustal stress logging prediction method based on rock physics model |
MX2016014193A (en) * | 2014-04-30 | 2017-05-03 | Halliburton Energy Services Inc | Characterizing a downhole environment using stiffness coefficients. |
-
2014
- 2014-12-30 CN CN201410842028.XA patent/CN104484573B/en active Active
-
2015
- 2015-12-23 AU AU2015275302A patent/AU2015275302B2/en active Active
- 2015-12-30 US US14/984,013 patent/US10031250B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630014B (en) * | 2008-07-16 | 2011-10-05 | 中国石油天然气集团公司 | Method for imaging anisotropic medium through utilization of vertical seismic profile data |
CN102466818A (en) * | 2010-11-11 | 2012-05-23 | 中国石油集团东方地球物理勘探有限责任公司 | Method for imaging anisotropic medium by utilizing cross-well seismic data |
CN102706965B (en) * | 2012-06-12 | 2014-08-27 | 西安邮电大学 | Novel accurate measuring method of anisotropic parameter of rock |
CN103323575A (en) * | 2013-05-20 | 2013-09-25 | 中国石油天然气股份有限公司 | Sedimentary rock reconstruction method and system |
Non-Patent Citations (8)
Title |
---|
Measurement of stiffness coefficients of anisotropic materials from pointlike generation and detection of acoustic waves;Audoin B, et al.,;《Journal of Applied Physics》;19961231;第80卷(第7期);3760-3771 * |
Seismic anisotropy in sedimentary rocks, part 1: A single‐plug laboratory method;Wang Z.;《Geophysics》;20021031;第67卷(第5期);1415~1422 * |
Weak Elastic Anisotropy;Thomsen L.;《 Geophysics》;19861031;第51卷(第10期);1954-1966 * |
从声波测井数据中反演地层各向异性参数的改进方法;王瑞甲,等;《声学技术》;20090430;第28卷(第2期);17-18 * |
储层沉积岩地震波速度及衰减与岩石物性间的关系;A.I.Best,等著,唐朗青,译;《天然气勘探与开发》;19950331;第18卷(第1期);27~42 * |
沉积储集岩的速度、衰减与岩石物理性质间的关系;A.I.Best,等;《国外油气勘探》;19950731;第7卷(第7期);425~441 * |
沉积岩中黏土矿物和常见非黏土矿物X射线衍射分析方法;国家能源局;《中华人民共和国石油天然气行业标准SY/T5163-2010》;20100501;1-47 * |
空隙度-渗透率-粘土含量对纵波速度的影响;Theodoros Klimentos;《国外油气勘探》;19930531;第5卷(第3期);307~317 * |
Also Published As
Publication number | Publication date |
---|---|
CN104484573A (en) | 2015-04-01 |
US10031250B2 (en) | 2018-07-24 |
US20160186557A1 (en) | 2016-06-30 |
AU2015275302A1 (en) | 2016-07-14 |
AU2015275302B2 (en) | 2017-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Paillet et al. | Acoustic waves in boreholes | |
Goertz-Allmann et al. | Combining microseismic and geomechanical observations to interpret storage integrity at the In Salah CCS site | |
CN105277982B (en) | A kind of mud shale total content of organic carbon earthquake prediction method | |
CN105988137B (en) | Logging character curve Reconstruction method of the sand-conglomerate body based on Core-Calibrated Logging | |
CN103867197A (en) | Complex lithology natural gas reservoir interval transit time discriminating method | |
CN104005760B (en) | Based on the crack detection method of azimuthal anisotropy elastic impedance | |
US20160320504A1 (en) | System and method for processing microseismic data | |
Suzuki et al. | Quantifying uncertainties in attenuation estimation at methane-hydrate-bearing zones using sonic waveform logs | |
Song et al. | Full-waveform based microseismic source mechanism studies in the Barnett Shale: Linking microseismicity to reservoir geomechanics | |
CN104155693A (en) | Angle gather seismic response numerical computation method of reservoir fluid fluidity | |
CN104484573B (en) | A kind of method of determination stratum stiffness coefficient | |
AU2014407527B2 (en) | Integrating vertical seismic profile data for microseismic anisotropy velocity analysis | |
CN110058323A (en) | A kind of tight sand formation brittleness index calculation method | |
CN107728205B (en) | A kind of Formation pressure prediction method | |
Su et al. | Elastic-wave evaluation of downhole hydraulic fracturing: Modeling and field applications | |
Cai et al. | Intelligent calculation method of relative sonic attenuation and its application to fracture evaluation in tight sandstone reservoir | |
CN107422384B (en) | A method of vertical seismic profile data is imitated using well-log information | |
WU et al. | Application status and research progress of dipole acoustic well logging | |
CN102809761A (en) | Control method for acquiring speed characteristics of coal rock | |
Xu et al. | Estimation of elastic wave velocity from logging-while-drilling measurement in unconsolidated formation of accretionary wedge in nankai trough | |
Tao et al. | Application of multipole array sonic logging to acid hydralic fracturing | |
Knutsen | On determination of Gmax by bender element and cross-hole testing | |
Sinha et al. | Hydraulic fracture characterization using cased-hole sonic data | |
Guo et al. | Multimodal dispersion curves extracted from deep seismic sounding profile to characterize the sedimentary structure | |
CN104516021A (en) | Ray elastic parameter inversion method capable of improving stability and precision of analysis formula |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |