CN107153096A - Shale reservoir into seam aptitude tests method and device - Google Patents

Shale reservoir into seam aptitude tests method and device Download PDF

Info

Publication number
CN107153096A
CN107153096A CN201710488733.8A CN201710488733A CN107153096A CN 107153096 A CN107153096 A CN 107153096A CN 201710488733 A CN201710488733 A CN 201710488733A CN 107153096 A CN107153096 A CN 107153096A
Authority
CN
China
Prior art keywords
sample
seam
wave
formula
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.)
Pending
Application number
CN201710488733.8A
Other languages
Chinese (zh)
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.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN201710488733.8A priority Critical patent/CN107153096A/en
Publication of CN107153096A publication Critical patent/CN107153096A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/14Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04Analysing solids
    • G01N29/07Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/02Details
    • G01N3/06Special adaptations of indicating or recording means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0019Compressive
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/06Indicating or recording means; Sensing means
    • G01N2203/0658Indicating or recording means; Sensing means using acoustic or ultrasonic detectors

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Acoustics & Sound (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The present invention provides a kind of shale reservoir into seam aptitude tests method and device, including:Each evaluation coefficient of sample is obtained, each evaluating includes Factor of Brittleness, speed stress sensitive sex index and the acoustic emission activity sex index of the sample;By the way that each evaluation coefficient is normalized, obtain the sample into seam merit rating coefficient.The shale reservoir that the present invention is provided into seam aptitude tests method and device, evaluation coefficient using the Factor of Brittleness of sample, speed stress sensitive sex index and acoustic emission activity sex index as evaluation shale reservoir sample into seam ability, so as to during being tested into seam ability shale reservoir, modeling fragility, intrinsic fracture development degree and the rupture activity for forming crack of shale are considered, so that test is more fully accurate, and it is simple to operate, it is adaptable to evaluation of the oil field compact reservoir into seam ability.

Description

Shale reservoir into seam aptitude tests method and device
Technical field
The present invention relates to shale testing field, more particularly to a kind of shale reservoir is into seam aptitude tests method and device.
Background technology
Shale gas exploitation has been included in China's petrol resources exploitation strategy.Due to shale reservoir matrix permeability it is extremely low, it is necessary to Production capacity could be formed by fracturing reform.Different from conventional reservoir pressure break, the pressure break of shale reservoir requires to form disperse co-volume Fracture network (volume fracturing).This depends on the ability and fracturing parameter of shale reservoir formation seam net.But simultaneously not all Reservoir can be transferred through large-scale pressing crack construction and realize the transformation of large-scale volume fracturing, therefore, shale reservoir is into seam ability Evaluate significant for preferred pressure break well section, prediction economic benefit.
Research shows whether reservoir possesses the condition for implementing volume transformation, is sent out depending on the modeling fragility of rock, intrinsic fracture Educate situation and anisotropy.Up to the present, domestic and international shale reservoir is into stitching in aptitude tests scheme, generally by shale reservoir " Factor of Brittleness " is equivalent into seam ability to evaluate.But actually fragility is not equivalent into seam ability, therefore, passes through existing page Rock reservoir is evaluated into the compressibility that seam aptitude tests scheme can not be exactly to shale.
The content of the invention
The present invention provides a kind of shale reservoir into seam aptitude tests method and device, for solve existing shale reservoir into Stitch the technical problem that aptitude tests scheme can not exactly to being evaluated into seam ability.
The first aspect of the invention is to provide a kind of shale reservoir into seam aptitude tests method, including:
Each evaluation coefficient of sample is obtained, it is quick that each evaluating includes the Factor of Brittleness of the sample, speed stress Perceptual index and acoustic emission activity sex index;
By the way that each evaluation coefficient is normalized, obtain the sample into seam merit rating coefficient.
Another aspect of the present invention is to provide a kind of shale reservoir into seam aptitude tests device, including:
Acquisition module, each evaluation coefficient for obtaining sample, each evaluating includes the fragility system of the sample Number, speed stress sensitive sex index and acoustic emission activity sex index;
Processing module, for by the way that each evaluation coefficient is normalized, obtaining being commented into seam ability for the sample Valency coefficient.
The shale reservoir that the present invention is provided is into seam aptitude tests method and device, by the Factor of Brittleness of sample, speed stress Sensitivity indices and acoustic emission activity sex index are as evaluation coefficient of the shale reservoir sample into seam ability is evaluated, so that right During shale reservoir is tested into seam ability, considered the modeling fragility of shale, intrinsic fracture development degree and The rupture activity in crack is formed, by the scheme for evaluating into seam ability jointly based on above-mentioned each evaluation coefficient so that test is more For comprehensively accurately, and simple to operate, it is adaptable to evaluation of the oil field compact reservoir into seam ability.
Brief description of the drawings
Fig. 1 be the embodiment of the present invention one provide shale reservoir into seam aptitude tests method schematic flow sheet;
Fig. 2 be the embodiment of the present invention two provide shale reservoir into seam aptitude tests method schematic flow sheet;
Fig. 3 be the embodiment of the present invention three provide shale reservoir into seam aptitude tests device structural representation.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described.
Schematic flow sheet such as Fig. 1 institute of the shale reservoir that Fig. 1 is provided for the embodiment of the present invention one into seam aptitude tests method Show, methods described includes:
101st, each evaluation coefficient of sample is obtained, each evaluating includes the Factor of Brittleness of the sample, speed should Power sensitivity indices and acoustic emission activity sex index.
In the present embodiment, intrinsic fracture development degree and anisotropy etc. are also into the key factor of seam merit rating. In practical application, experimental provision can be set up in laboratory, for standard sample (one inch or two inches of small core columns), is opened Ultrasonic tesint is studied in exhibition room.Evaluate modeling fragility, the method for cranny development of shale reservoir jointly using ultrasonic wave and sound emission, Carry out shale reservoir into seam Comprehensive Evaluation On Ability.
Wherein, the sample can be the compact reservoir rock core for being derived from underground or appearing, such as shale, tight sand isolith The heart, rock core is cut and 25 millimeters of (mm) × 50mm core columns, 38mm × 76mm core columns or 50*100mm rock core is polished into Post (one inch or two inches of small core columns), the side wall of sample and two end faces keep smooth, and wherein two end faces of sample are parallel Degree is uniform in pressurized condition lower stress to ensure within ± 0.01mm.Accordingly, before 101, methods described can also be wrapped Include:Compact reservoir rock core is prepared into the sample, the sample is cylinder, two faces parallels of the sample ± In the range of 0.01mm.
Optionally, the modeling fragility of shale reservoir can by rock core brittle mineral component ratio, modulus of elasticity and Poisson's ratio, The ratio between compression strength and tensile strength etc. are evaluated.In practical application, it can reflect that the modeling of shale samples is crisp by Factor of Brittleness Property, further, Fragility parameters can be reflected by vertical transverse wave speed.
Specifically, it can be to carry out compressional wave and shear wave under normal temperature and pressure to sample first to obtain the process of Factor of Brittleness Speed measurement, wave speed measurement can use ultrasonic pulse technology, that is, measure sample interior known paths on ultrasonic wave be transmitted through when Between.For example, Fig. 2 is the waveform diagram arrived using ultrasonic pulse technical monitoring in the embodiment of the present invention one, such as Fig. 2 institutes Show, first most compressional wave P ripple arrival times above, first most shear wave S ripple arrival times below, read the arrival of P ripples and S ripples Time, according to known distance, you can calculate the compressional wave and transverse wave speed of standard sample, further utilize the vertical shear wave obtained Speed calculates Young's modulus and Poisson's ratio, have rated the brittleness index of sample.Accordingly, it is described to obtain the crisp of the sample in 101 Property coefficient, can specifically include:
Using ultrasonic pulse technology, the diametric longitudinal wave velocity and transverse wave speed along the sample are gathered;
According to the longitudinal wave velocity, the transverse wave speed and the 3rd formula, calculate and obtain Young's modulus and Poisson's ratio, it is described 3rd formula is:Wherein, VpFor the longitudinal wave velocity, VsFor the shear wave Velocity of wave, E is the Young's modulus, and ν is the Poisson's ratio;
According to the 4th formula, the Young's modulus is normalized, the Young's modulus after being normalized, the described 4th Formula is:Ebrit=(E-10)/(80-10) × 100, wherein, EbritFor the Young's modulus after normalization;
According to the 5th formula, the Poisson's ratio is normalized, the Poisson's ratio after being normalized, the 5th formula For:νbrit=(0.4- ν)/(0.4-0.15) × 100, wherein, νbritFor the Poisson's ratio after normalization;
According to the Young's modulus and Poisson's ratio and the 6th formula after the normalization, calculate and obtain the Factor of Brittleness, 6th formula is:BI=0.5Ebrit+0.5υbrit, wherein, BI is the Factor of Brittleness.
Wherein, the Young's modulus E after normalizationbritUnit be specifically as follows Gpa.In present embodiment, by fragility system Count to reflect the modeling fragility of shale samples, specifically, Factor of Brittleness is bigger, then the fragility of shale is stronger.For example, fragility is worked as Coefficient B I>When 40, rock is brittle;BI>When 60, rock brittleness is then very strong.
Based on above-mentioned embodiment, in practical application, it is possible to use ultrasonic probe monitors velocity of wave.Specifically, by a pair Probe is fixed on the side of standard sample, the diametrically opposite velocity of wave of measurement standard sample.Occur using pulse signal generator Signal, by transmitting probe, ultrasonic wave is propagated in the sample, is received by the corresponding receiving transducer of opposite side, is read arriving for waveform Up to the time, compressional wave and transverse wave speed are obtained.Then corresponding, the utilization ultrasonic pulse technology gathers the diameter along the sample The longitudinal wave velocity and transverse wave speed in direction, can specifically include:
Transmitting probe and receiving transducer are symmetrical arranged on the side wall of the sample, using ultrasonic pulse technology, institute is obtained State the compressional wave and shear wave waveform in sample;
According to the compressional wave and shear wave waveform, the diametric(al) along the sample is obtained, the compressional wave and the shear wave are worn Cross the duration used in the sample;
Diameter difference divided by the compressional wave and the shear wave by calculating the sample are obtained through the duration used in sample Obtain the longitudinal wave velocity and the transverse wave speed.
Again specifically, can be by carrying out uniaxial loading to standard sample, and in loading procedure after obtaining Factor of Brittleness Middle record per unit stress, for example, every MPa (Mpa), change obtains velocity of wave.It is quick that with stress variation its stress is obtained using velocity of wave Perception.Accordingly, the speed stress sensitive sex index of the sample is obtained described in 101, can specifically be included:
By default step-length, the pressure of increase uniaxial loading on to the sample, until the sample burst;
Before the sample burst, by gathering the velocity of wave in sample diameter direction described in each first moment lower edge, institute is obtained Relation of the velocity of wave with stress variation is stated, first moment applies pressure phase on to the sample when being and gathering velocity of wave with last time Than when the current pressure applied on to the sample adds default threshold value at the time of correspondence;
According to relation of the velocity of wave with stress variation, the speed stress sensitive sex index of the sample is obtained.
Again specifically, the acoustic emission activity sex index can specifically include acoustic emission rate, wherein, the acoustic emission rate is The sound emission number occurred in unit interval., can be using ultrasonic probe note during foregoing uniaxial loading in practical application The acoustic emission rate produced by the extension of crack is recorded, to obtain the activity index of reservoir core rupture.Accordingly, the acquisition institute The acoustic emission activity sex index of sample is stated, can specifically be included:
By default step-length, it is all that the pressure of increase uniaxial loading on to the sample, and record occur during this period Sound emission number, until the sample burst;
According to the sound emission number, calculate and obtain the acoustic emission rate.
Optionally, it is possible to use ultrasonic probe monitors sound emission.Accordingly, it is described to press default step-length, increase single shaft The pressure of loading on to the sample, and institute's acoustic emission number that record occurs during this period, until the sample burst, specifically It can include:
By default step-length, the pressure of increase uniaxial loading on to the sample, and utilization is arranged on sample side wall On ultrasonic probe, record institute acoustic emission number, until the sample burst.
102nd, by the way that each evaluation coefficient is normalized, obtain the sample into seam merit rating coefficient.
Specifically, using Factor of Brittleness, speed stress sensitive property coefficient, the activity index ruptured, setting up a synthesis Evaluate compact reservoir sample into seam ability method, this method consider the modeling fragility of rock, intrinsic fracture development degree with And eventually form the rupture activity in crack.Generally speaking, in the present embodiment, reservoir into seam Comprehensive Evaluation On Ability index (Fracbility Index) is made up of three parts, i.e. Factor of Brittleness, stress sensitive property coefficient or fracture density, acoustic emission activity Property.Accordingly, 102 can specifically include:
According to the first formula, each evaluation coefficient is normalized respectively, each evaluation coefficient after being normalized, First formula is:Wherein, XidFor the evaluation coefficient after normalization;
According to each evaluation coefficient and the second formula after normalization, calculate and obtain described into seam merit rating coefficient, it is described Second formula is:Wherein, FI is into seam Comprehensive Evaluation On Ability coefficient, and n is the individual of each evaluation coefficient Number.
Wherein, it is specifically as follows percentages into seam Comprehensive Evaluation On Ability coefficient FI, is respectively evaluated described in the present embodiment The number n of coefficient value is specially 3, XidLater single evaluation coefficient is respectively normalized, that is, represents Factor of Brittleness, stress Sensitivity coefficient, AE activity etc..
Measured specifically, carrying out compressional wave and transverse wave speed under normal temperature and pressure to standard sample first, utilize the speed of acquisition Young's modulus and Poisson's ratio are calculated, the brittleness index of sample is evaluated.And last rupture is forced into core sample progress single shaft, this During sound emission produced by the extension of crack have recorded using ultrasonic probe, and obtain per MPa stress variations velocity of wave.Utilize Velocity of wave obtains its stress sensitivity with stress variation, and the activity index of reservoir core rupture is obtained using sound emission.Using crisp Property coefficient, stress sensitive property coefficient, rupture activity index, overall merit compact reservoir sample into seam ability.
The shale reservoir that the present embodiment is provided is into seam aptitude tests method, by the Factor of Brittleness of sample, speed stress sensitive Sex index and acoustic emission activity sex index are as evaluation coefficient of the shale reservoir sample into seam ability is evaluated, so as to shale During reservoir is tested into seam ability, modeling fragility, intrinsic fracture development degree and the formation of shale have been considered The rupture activity in crack, by the scheme for evaluating into seam ability jointly based on above-mentioned each evaluation coefficient so that test is more complete Face is accurate, and simple to operate, it is adaptable to evaluation of the oil field compact reservoir into seam ability.
Fig. 2 be the embodiment of the present invention two provide shale reservoir into seam aptitude tests method schematic flow sheet, such as Fig. 2 institutes Show, methods described includes:
201st, compact reservoir rock core is prepared into sample;
202nd, at least three ultrasonic probes are arranged on the side wall of the sample, a probe is used for recording sound emission, one Probe is used for monitor along the diametric longitudinal wave velocity of standard sample and transverse wave speed;
203rd, the Young's modulus and Poisson's ratio of sample are obtained using the p-and s-wave velocity of monitoring, and assesses its brittleness index;
204th, to sample uniaxial loading to rupturing, units MPa record compressional wave and transverse wave speed are often increased, and respective record is broken Institute's acoustic emission number of generation before splitting;
205th, to acquisition brittleness index, speed stress sensitivity and rupture sound emission number is normalized, and synthesis is commented Valency compact reservoir into seam ability.
Specifically, the detailed process and method in the present embodiment are referred to the related content in embodiment one.
The shale reservoir that the present embodiment is provided is into seam aptitude tests method, by the Factor of Brittleness of sample, speed stress sensitive Sex index and acoustic emission activity sex index are as evaluation coefficient of the shale reservoir sample into seam ability is evaluated, so as to shale During reservoir is tested into seam ability, modeling fragility, intrinsic fracture development degree and the formation of shale have been considered The rupture activity in crack, by the scheme for evaluating into seam ability jointly based on above-mentioned each evaluation coefficient so that test is more complete Face is accurate, and simple to operate, it is adaptable to evaluation of the oil field compact reservoir into seam ability.
Fig. 3 be the embodiment of the present invention three provide shale reservoir into seam aptitude tests device structural representation, such as Fig. 3 institutes Show, described device includes:
Acquisition module 31, each evaluation coefficient for obtaining sample, each evaluating includes the fragility of the sample Coefficient, speed stress sensitive sex index and acoustic emission activity sex index;
Processing module 32, for by the way that each evaluation coefficient is normalized, obtain the sample into seam ability Evaluation coefficient.
In practical application, experimental provision can be set up in laboratory, for standard sample (one inch or two inches of small rocks Stem), carry out indoor ultrasonic tesint research.Evaluate modeling fragility, the crack hair of shale reservoir jointly using ultrasonic wave and sound emission The method educated, carry out shale reservoir into seam Comprehensive Evaluation On Ability.
Wherein, the sample can be the compact reservoir rock core for being derived from underground or appearing, such as shale, tight sand isolith The heart, rock core is cut and the core column (one of 25mm × 50mm core columns, 38mm × 76mm core columns or 50*100mm is polished into Inch or two inches of small core columns), the side wall of sample and two end faces keep smooth, wherein two faces parallels of sample ± It is uniform in pressurized condition lower stress to ensure within 0.01mm.Accordingly, two faces parallels of the sample ± In the range of 0.01mm.
Optionally, the modeling fragility of shale reservoir can by rock core brittle mineral component ratio, modulus of elasticity and Poisson's ratio, The ratio between compression strength and tensile strength etc. are evaluated.In practical application, it can reflect that the modeling of shale samples is crisp by Factor of Brittleness Property, further, Fragility parameters can be reflected by vertical transverse wave speed.
Specifically, it can be to carry out compressional wave and shear wave under normal temperature and pressure to sample first to obtain the process of Factor of Brittleness Speed measurement, wave speed measurement can use ultrasonic pulse technology, that is, measure sample interior known paths on ultrasonic wave be transmitted through when Between.Accordingly, acquisition module 31 can include:
First collecting unit, for utilizing ultrasonic pulse technology, gathers the diametric longitudinal wave velocity along the sample And transverse wave speed;
First computing unit, for according to the longitudinal wave velocity, the transverse wave speed and the 3rd formula, calculating and obtaining Young Modulus and Poisson's ratio, the 3rd formula is:Wherein, VpFor the compressional wave Velocity of wave, VsFor the transverse wave speed, E is the Young's modulus, and ν is the Poisson's ratio;
First computing unit, is additionally operable to, according to the 4th formula, the Young's modulus be normalized, and obtains normalizing Young's modulus after change, the 4th formula is:Ebrit=(E-10)/(80-10) × 100, wherein, EbritAfter normalization Young's modulus;
First computing unit, is additionally operable to, according to the 5th formula, the Poisson's ratio be normalized, is normalized Poisson's ratio afterwards, the 5th formula is:νbrit=(0.4- ν)/(0.4-0.15) × 100, wherein, νbritAfter normalization Poisson's ratio;
First computing unit, is additionally operable to public according to the Young's modulus and Poisson's ratio after the normalization and the 6th Formula, calculates and obtains the Factor of Brittleness, and the 6th formula is:BI=0.5Ebrit+0.5υbrit, wherein, BI is the fragility system Number.
Wherein, the Young's modulus E after normalizationbritUnit be specifically as follows Gpa.In present embodiment, by fragility system Count to reflect the modeling fragility of shale samples, specifically, Factor of Brittleness is bigger, then the fragility of shale is stronger.
Based on above-mentioned embodiment, in practical application, it is possible to use ultrasonic probe monitors velocity of wave.Specifically, by a pair Probe is fixed on the side of standard sample, the diametrically opposite velocity of wave of measurement standard sample.Occur using pulse signal generator Signal, by transmitting probe, ultrasonic wave is propagated in the sample, is received by the corresponding receiving transducer of opposite side, is read arriving for waveform Up to the time, compressional wave and transverse wave speed are obtained.Then corresponding, first collecting unit includes:
First collection subelement, is visited for the transmitting probe by being symmetricly set on the side wall of the sample and reception Head, using ultrasonic pulse technology, obtains the compressional wave and shear wave waveform in the sample;
First processing subelement, for according to the compressional wave and shear wave waveform, obtaining the diametric(al) along the sample, institute Compressional wave and the shear wave are stated through the duration used in the sample;
First computation subunit, for by calculating the diameter of the sample respectively divided by the compressional wave and the shear wave are worn The duration used in sample is crossed, the longitudinal wave velocity and the transverse wave speed is obtained.
Again specifically, can be by carrying out uniaxial loading to standard sample, and in loading procedure after obtaining Factor of Brittleness Middle record per unit stress, for example, per Mpa, change obtains velocity of wave.Using velocity of wave its stress sensitivity is obtained with stress variation. Accordingly, acquisition module 31 can also include:
First presser unit, for by default step-length, the pressure of increase uniaxial loading on to the sample, until described Sample burst;
Second collecting unit, for before the sample burst, by gathering sample diameter described in each first moment lower edge The velocity of wave in direction, obtains relation of the velocity of wave with stress variation, and first moment is applied to when being and gathering velocity of wave with last time Pressure on the sample is compared, when the current pressure applied on to the sample adds default threshold value at the time of correspondence;
First processing units, for according to relation of the velocity of wave with stress variation, obtaining the speed stress of the sample Sensitivity indices.
Again specifically, the acoustic emission activity sex index can specifically include acoustic emission rate, wherein, the acoustic emission rate is The sound emission number occurred in unit interval., can be using ultrasonic probe note during foregoing uniaxial loading in practical application The acoustic emission rate produced by the extension of crack is recorded, to obtain the activity index of reservoir core rupture.Accordingly, acquisition module 31 It can also include:
Second processing unit, for existing by default step-length, the pressure of increase uniaxial loading on to the sample, and record The institute's acoustic emission number occurred during this, until the sample burst;
Second computing unit, for according to the sound emission number, calculating and obtaining the acoustic emission rate.
Optionally, it is possible to use ultrasonic probe monitors sound emission.Accordingly, the second processing unit, specifically for By default step-length, the pressure of increase uniaxial loading on to the sample, and utilize the ultrasound being arranged on the wall of the sample side Probe, record institute acoustic emission number, until the sample burst.
Specifically, using Factor of Brittleness, speed stress sensitive property coefficient, the activity index ruptured, setting up a synthesis Evaluate compact reservoir sample into seam ability method, this method consider the modeling fragility of rock, intrinsic fracture development degree with And eventually form the rupture activity in crack.Generally speaking, in the present embodiment, reservoir into seam Comprehensive Evaluation On Ability index (Fracbility Index) is made up of three parts, i.e. Factor of Brittleness, stress sensitive property coefficient or fracture density, acoustic emission activity Property.Accordingly, processing module 32 includes:
Normalization unit, for according to the first formula, being normalized, being normalized to each evaluation coefficient respectively Each evaluation coefficient afterwards, first formula is:Wherein, XidFor the evaluation coefficient after normalization;
3rd computing unit, obtains described into seam for according to each evaluation coefficient and the second formula after normalization, calculating Merit rating coefficient, second formula is:Wherein, FI is into seam Comprehensive Evaluation On Ability coefficient, and n is institute State the number of each evaluation coefficient.
Wherein, it is specifically as follows percentages, X into seam Comprehensive Evaluation On Ability coefficient FIidRespectively normalize later Single evaluation coefficient, that is, represent Factor of Brittleness, stress sensitive property coefficient, AE activity etc..
Measured specifically, carrying out compressional wave and transverse wave speed under normal temperature and pressure to standard sample first, utilize the speed of acquisition Young's modulus and Poisson's ratio are calculated, the brittleness index of sample is evaluated.And last rupture is forced into core sample progress single shaft, this During sound emission produced by the extension of crack have recorded using ultrasonic probe, and obtain per MPa stress variations velocity of wave.Utilize Velocity of wave obtains its stress sensitivity with stress variation, and the activity index of reservoir core rupture is obtained using sound emission.Using crisp Property coefficient, stress sensitive property coefficient, rupture activity index, overall merit compact reservoir sample into seam ability.
The shale reservoir that the present embodiment is provided is into seam aptitude tests device, by the Factor of Brittleness of sample, speed stress sensitive Sex index and acoustic emission activity sex index are as evaluation coefficient of the shale reservoir sample into seam ability is evaluated, so as to shale During reservoir is tested into seam ability, modeling fragility, intrinsic fracture development degree and the formation of shale have been considered The rupture activity in crack, by the scheme for evaluating into seam ability jointly based on above-mentioned each evaluation coefficient so that test is more complete Face is accurate, and simple to operate, it is adaptable to evaluation of the oil field compact reservoir into seam ability.
It is apparent to those skilled in the art that, for convenience and simplicity of description, the device of foregoing description Specific work process, may be referred to the corresponding process in preceding method embodiment, will not be repeated here.
It should be noted that each functional unit in various embodiments of the present invention can be integrated in a processing unit, Can be that unit is individually physically present, can also two or more units it is integrated in a unit.It is above-mentioned integrated Unit can both be realized in the form of hardware, it would however also be possible to employ hardware adds the form of SFU software functional unit to realize.
One of ordinary skill in the art will appreciate that:Realizing all or part of step of above-mentioned each method embodiment can lead to The related hardware of programmed instruction is crossed to complete.Foregoing program can be stored in a computer read/write memory medium.The journey Sequence upon execution, performs the step of including above-mentioned each method embodiment;And foregoing storage medium includes:ROM, RAM, magnetic disc or Person's CD etc. is various can be with the medium of store program codes.
Finally it should be noted that:Various embodiments above is merely illustrative of the technical solution of the present invention, rather than its limitations;To the greatest extent The present invention is described in detail with reference to foregoing embodiments for pipe, it will be understood by those within the art that:Its according to The technical scheme described in foregoing embodiments can so be modified, or which part or all technical characteristic are entered Row equivalent substitution;And these modifications or replacement, the essence of appropriate technical solution is departed from various embodiments of the present invention technology The scope of scheme.

Claims (10)

1. a kind of shale reservoir into seam aptitude tests method, it is characterised in that including:
Each evaluation coefficient of sample is obtained, each evaluating includes the Factor of Brittleness of the sample, speed stress sensitivity Index and acoustic emission activity sex index;
By the way that each evaluation coefficient is normalized, obtain the sample into seam merit rating coefficient.
2. according to the method described in claim 1, it is characterised in that described by the way that each evaluation coefficient is normalized, Obtain the sample into seam merit rating coefficient, including:
According to the first formula, each evaluation coefficient is normalized respectively, each evaluation coefficient after being normalized is described First formula is:Wherein, XidFor the evaluation coefficient after normalization;
According to each evaluation coefficient and the second formula after normalization, calculate and obtain described into seam merit rating coefficient, described second Formula is:Wherein, FI is into seam Comprehensive Evaluation On Ability coefficient, and n is the number of each evaluation coefficient.
3. method according to claim 1 or 2, it is characterised in that the Factor of Brittleness of the acquisition sample includes:
Using ultrasonic pulse technology, the diametric longitudinal wave velocity and transverse wave speed along the sample are gathered;
According to the longitudinal wave velocity, the transverse wave speed and the 3rd formula, calculate and obtain Young's modulus and Poisson's ratio, the described 3rd Formula is:Wherein, VpFor the longitudinal wave velocity, VsFor the shear wave Speed, E is the Young's modulus, and ν is the Poisson's ratio;
According to the 4th formula, the Young's modulus is normalized, the Young's modulus after being normalized, the 4th formula For:Ebrit=(E-10)/(80-10) × 100, wherein, EbritFor the Young's modulus after normalization;
According to the 5th formula, the Poisson's ratio is normalized, the Poisson's ratio after being normalized, the 5th formula is: νbrit=(0.4- ν)/(0.4-0.15) × 100, wherein, νbritFor the Poisson's ratio after normalization;
According to the Young's modulus and Poisson's ratio and the 6th formula after the normalization, calculate and obtain the Factor of Brittleness, it is described 6th formula is:BI=0.5Ebrit+0.5υbrit, wherein, BI is the Factor of Brittleness.
4. method according to claim 1 or 2, it is characterised in that the speed stress sensitivity of the acquisition sample Index, including:
By default step-length, the pressure of increase uniaxial loading on to the sample, until the sample burst;
Before the sample burst, by gathering the velocity of wave in sample diameter direction described in each first moment lower edge, the ripple is obtained The fast relation with stress variation, compared with first moment is applies pressure on to the sample when last time gathering velocity of wave, When the current pressure applied on to the sample adds default threshold value at the time of correspondence;
According to relation of the velocity of wave with stress variation, the speed stress sensitive sex index of the sample is obtained.
5. method according to claim 1 or 2, it is characterised in that the acoustic emission activity sex index includes acoustic emission rate, The acoustic emission rate is the sound emission number occurred in the unit time.
6. method according to claim 5, it is characterised in that the acoustic emission activity sex index of the acquisition sample, Including:
By default step-length, the pressure of increase uniaxial loading on to the sample, and record the sound hair of institute occurred during this period Number is penetrated, until the sample burst;
According to the sound emission number, calculate and obtain the acoustic emission rate.
7. method according to claim 3, it is characterised in that the utilization ultrasonic pulse technology, is gathered along the sample Diametric longitudinal wave velocity and transverse wave speed, including:
Transmitting probe and receiving transducer are symmetrical arranged on the side wall of the sample, using ultrasonic pulse technology, the sample is obtained Compressional wave and shear wave waveform in product;
According to the compressional wave and shear wave waveform, the diametric(al) along the sample is obtained, the compressional wave and the shear wave pass through institute State the duration used in sample;
Diameter difference divided by the compressional wave and the shear wave by calculating the sample obtain institute through the duration used in sample State longitudinal wave velocity and the transverse wave speed.
8. method according to claim 6, it is characterised in that described to press default step-length, increase uniaxial loading is described Pressure on sample, and the institute's acoustic emission number occurred during this period is recorded, up to the sample burst, including:
By default step-length, the pressure of increase uniaxial loading on to the sample, and utilization are arranged on the wall of the sample side Ultrasonic probe, record institute acoustic emission number, until the sample burst.
9. according to the method described in claim 1, it is characterised in that before each evaluation coefficient for obtaining sample, in addition to:
Compact reservoir rock core is prepared into the sample, the sample is cylinder, and two faces parallels of the sample exist In the range of ± 0.01mm.
10. a kind of shale reservoir into seam aptitude tests device, it is characterised in that including:
Acquisition module, each evaluation coefficient for obtaining sample, each evaluating includes Factor of Brittleness, the speed of the sample Spend stress sensitive sex index and acoustic emission activity sex index;
Processing module, for by the way that each evaluation coefficient is normalized, obtain the sample into seam merit rating system Number.
CN201710488733.8A 2017-06-23 2017-06-23 Shale reservoir into seam aptitude tests method and device Pending CN107153096A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710488733.8A CN107153096A (en) 2017-06-23 2017-06-23 Shale reservoir into seam aptitude tests method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710488733.8A CN107153096A (en) 2017-06-23 2017-06-23 Shale reservoir into seam aptitude tests method and device

Publications (1)

Publication Number Publication Date
CN107153096A true CN107153096A (en) 2017-09-12

Family

ID=59796409

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710488733.8A Pending CN107153096A (en) 2017-06-23 2017-06-23 Shale reservoir into seam aptitude tests method and device

Country Status (1)

Country Link
CN (1) CN107153096A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108204232A (en) * 2017-12-12 2018-06-26 中国石油天然气股份有限公司 Reservoir fracture development process monitoring method and device and computer storage medium
CN108645596A (en) * 2018-05-11 2018-10-12 北京石油化工学院 A kind of experimental method for evaluating fluid to rock width generation capacity
CN109030205A (en) * 2018-07-06 2018-12-18 沧州子芩信息科技有限公司 Metal material brittleness test method and test device
CN109374757A (en) * 2018-10-30 2019-02-22 中国海洋石油集团有限公司 The method for handling evaluation igneous reservoirs validity using wave amplitude quantification
CN109682676A (en) * 2018-12-29 2019-04-26 上海工程技术大学 A kind of feature extracting method of the acoustic emission signal of fiber tension failure
CN112253101A (en) * 2020-10-09 2021-01-22 中国石油大学(北京) Oil and gas resource exploration method, device, equipment and computer readable storage medium
CN112485401A (en) * 2020-11-10 2021-03-12 成都理工大学 Quantitative experiment method for shale physical property response under influence of multi-factor coupling
CN112649858A (en) * 2019-10-11 2021-04-13 中国石油化工股份有限公司 Shale brittleness prediction method and system based on core test
CN114859435A (en) * 2021-02-04 2022-08-05 中国石油天然气股份有限公司 Shale compressibility evaluation parameter acquisition method and device and computer equipment

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104131803A (en) * 2013-05-10 2014-11-05 中国石油大学(北京) Experiment method for evaluating shale fracturing crack mesh forming capability
CN205353018U (en) * 2016-01-25 2016-06-29 中国石油大学(北京) Utilize ultrasonic wave to survey anisotropy of rock mass's simple and easy experimental apparatus
CN106837322A (en) * 2017-01-23 2017-06-13 中国石油大学(北京) Shale reservoir into seam Capability index acquisition methods and device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104131803A (en) * 2013-05-10 2014-11-05 中国石油大学(北京) Experiment method for evaluating shale fracturing crack mesh forming capability
CN205353018U (en) * 2016-01-25 2016-06-29 中国石油大学(北京) Utilize ultrasonic wave to survey anisotropy of rock mass's simple and easy experimental apparatus
CN106837322A (en) * 2017-01-23 2017-06-13 中国石油大学(北京) Shale reservoir into seam Capability index acquisition methods and device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DAOBING WANG等: "A novel experiment approach for fracability evaluation in tight-gas reservoirs", 《JOURNAL OF NATURAL GAS SCIENCE AND ENGINEERING》 *
郑丽娜: "英台火山岩储层体积压裂适应性评价方法研究", 《2015年全国天然气学术年会》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108204232A (en) * 2017-12-12 2018-06-26 中国石油天然气股份有限公司 Reservoir fracture development process monitoring method and device and computer storage medium
CN108204232B (en) * 2017-12-12 2021-04-30 中国石油天然气股份有限公司 Reservoir fracture development process monitoring method and device and computer storage medium
CN108645596A (en) * 2018-05-11 2018-10-12 北京石油化工学院 A kind of experimental method for evaluating fluid to rock width generation capacity
CN109030205A (en) * 2018-07-06 2018-12-18 沧州子芩信息科技有限公司 Metal material brittleness test method and test device
CN109374757A (en) * 2018-10-30 2019-02-22 中国海洋石油集团有限公司 The method for handling evaluation igneous reservoirs validity using wave amplitude quantification
CN109682676A (en) * 2018-12-29 2019-04-26 上海工程技术大学 A kind of feature extracting method of the acoustic emission signal of fiber tension failure
CN112649858A (en) * 2019-10-11 2021-04-13 中国石油化工股份有限公司 Shale brittleness prediction method and system based on core test
CN112253101A (en) * 2020-10-09 2021-01-22 中国石油大学(北京) Oil and gas resource exploration method, device, equipment and computer readable storage medium
CN112485401A (en) * 2020-11-10 2021-03-12 成都理工大学 Quantitative experiment method for shale physical property response under influence of multi-factor coupling
CN114859435A (en) * 2021-02-04 2022-08-05 中国石油天然气股份有限公司 Shale compressibility evaluation parameter acquisition method and device and computer equipment
CN114859435B (en) * 2021-02-04 2024-07-26 中国石油天然气股份有限公司 Shale compressibility evaluation parameter acquisition method and device and computer equipment

Similar Documents

Publication Publication Date Title
CN107153096A (en) Shale reservoir into seam aptitude tests method and device
CN106053231A (en) Testing device for anisotropism of shale in true-triaxial condition and testing method of testing device
Mehan et al. Analysis of composite failure mechanisms using acoustic emissions
CN103868993B (en) The acoustics method of discrimination of rock three axle list sample fado level yield point and device
US9606036B2 (en) Method for determining geomechanical parameters of a rock sample
CN110864966B (en) Compressibility evaluation method suitable for fractured rock
CN104775810A (en) Method for evaluating compressibility of shale gas reservoir
CN106248494A (en) A kind of method for shale gas well reservoir fragility overall merit
CN105203639B (en) A method of detecting rock brittleness using acoustic-emission
CN104865124A (en) Shale brittleness index determination method based on rock stress-strain curve and ultrasonic longitudinal wave velocity
Carino Impact echo: The fundamentals
CN106896410B (en) The method of the deformation modulus and brittleness index of rock is explained using Sonic Logging Data
CN104747163A (en) Method and device for identifying reservoir fractures in tight sandstone
CN103513270A (en) Gas reservoir identification evaluating method based on acoustic characteristic of rock and device thereof
CN104931363A (en) Jointed rock deformation modulus testing method
Shi et al. A brittleness index evaluation method for weak-brittle rock by acoustic emission technique
CN102183585B (en) Rock core sampling method
CN205879690U (en) Be used for shale anisotropic testing arrangement under true triaxial condition
CN117350145A (en) Rock mechanical parameter while-drilling prediction method based on vibration signals
CN109828031A (en) Rock brittleness evaluation method and device
Abbas et al. Anisotropic index strength behaviour and failure mode validation of weathered shale
CN105068120A (en) Sound wave experiment method and identification method for tight sandstone fracture
CN117148435A (en) Direct inversion method of compact sandstone brittleness index factor
Li et al. Quantitative detection of damage processes in granite by sound signals
CN105352986B (en) Low-field nuclear magnetic resonance rock organic matter detection method and device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170912