CN110529087A - A kind of method and device for evaluating stratum hydraulic fracturing effect - Google Patents
A kind of method and device for evaluating stratum hydraulic fracturing effect Download PDFInfo
- Publication number
- CN110529087A CN110529087A CN201910644976.5A CN201910644976A CN110529087A CN 110529087 A CN110529087 A CN 110529087A CN 201910644976 A CN201910644976 A CN 201910644976A CN 110529087 A CN110529087 A CN 110529087A
- Authority
- CN
- China
- Prior art keywords
- pressure break
- backward
- depth
- log data
- wave
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 230000000694 effects Effects 0.000 title claims abstract description 52
- 230000008901 benefit Effects 0.000 claims abstract description 6
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims description 46
- 239000011435 rock Substances 0.000 claims description 43
- 238000001914 filtration Methods 0.000 claims description 12
- 230000001427 coherent effect Effects 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 9
- 238000005314 correlation function Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 235000013399 edible fruits Nutrition 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 2
- 239000004568 cement Substances 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 23
- 238000005516 engineering process Methods 0.000 abstract description 8
- 238000010606 normalization Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052796 boron Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/14—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Primary Health Care (AREA)
- Economics (AREA)
- Marine Sciences & Fisheries (AREA)
- Animal Husbandry (AREA)
- Agronomy & Crop Science (AREA)
- Acoustics & Sound (AREA)
- Remote Sensing (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention proposes a kind of method and device for evaluating stratum hydraulic fracturing effect, which comprises carries out array sonic log in depth intervals, constructs the dipole acoustic log data of different direction;The dipole acoustic log data of the different direction are filtered, normalized, the forward and backward scattered wave energy envelope of pressure break is calculated, the difference association evaluation well week fracturing effect of the forward and backward dipole acoustic log shear wave velocity difference of pressure break and scattering wave energy is utilized.Compared to existing evaluation of Fracturing Effect on Compact Sandstone technology, this method advantage is that the fracturing effect along pit shaft direction can be evaluated, and can also realize the evaluation of Fracturing Effect on Compact Sandstone within the scope of tens of meters of radial direction of well week.
Description
Technical field
The invention belongs to acoustic logging fields, and in particular to a kind of to evaluate waterpower pressure using dipole shear wave scattering effect
Split the square law device of effect.
Background technique
Unconventional reservoir (such as shale gas) permeability is low, and the major way for improving permeability of reservoir is to generate complexity
Chicken-wire cracking, fracturing yield increasing become the key technology of untraditional reservoir exploitation as the important means for developing fine and close reservoir
One of.After unconventional reservoir implements fracturing reform measure, effective method is needed to determine fracturing work effect, is obtained pressure break and is lured
All multi informations such as fracture condudtiviy, geometric shape, complexity and its orientation are led, shale gas reservoir fracture stimulation operations effect is improved
And gas well deliverability, and improve gas recovery factor.Especially in shale reservoir, the calculating pair of oil reservoir volume (SRV) is transformed
There is important engineering significance in formulating pressing crack construction scheme, evaluation fracturing effect, prediction shale gas yield.
Currently, micro-seismic technology is widely used in hydrofracturing effect assessment, it can be with the shape of dynamic monitoring fracturing fracture
At process.Although micro-seismic technology can effectively in all hundreds of meters of the ranges of detecting shaft fracturing fracture dynamic spread, be
The microseism signal for obtaining high-precision high s/n ratio generally requires to carry out offset well monitoring, this can additionally increase oil-gas exploration
Cost, especially in the higher offshore oil and gas field of exploration cost.In addition to this, evaluation wellbore fracture fracture height is commonly logged well
Method mainly has temperature logging, radioisotope logging, note boron neutron well logging, compensated neutron log, dipole acoustic log etc., thinking
It is to obtain fracturing fracture elevation information according to the difference of the forward and backward well-log information of pressure break, evaluates fracturing effect.Wherein dipole
Acoustic logging mainly evaluates fracture height by the anisotropy of the time difference before and after comparison pressure break, and this method technology maturation is implemented
Convenient, pollution-free, when differentiation, is intuitively quick, is the most common logging method of current evaluation fracturing fracture height.However, this method
The effect that pressure break generates in pit shaft direction can only be assessed, and stratum radial direction fracturing effect cannot be detected, and when well week pressure break is
When netted seam, the technical failure.
Summary of the invention
For the deficiency of existing fracturing effect evaluation method, the present invention is proposed and is disclosed a kind of forward and backward using pressure break
The method that dipole acoustic log scatters the difference evaluation fracturing effect of wave energy.
First aspect present invention proposes a kind of method for evaluating stratum hydraulic fracturing effect, which comprises
S1, array sonic log is carried out in depth intervals, obtain the forward and backward cross-dipole sound of pressure break in depth intervals
Wave well logging time domain four componets data and instrumental azimuth curve;
S2, by before the pressure break four componets data and instrumental azimuth curve terrestrial coordinate system is transformed by instrument coordinates system
Under, construct the dipole acoustic log data of different direction;
S3, the dipole acoustic log data of the different direction are filtered, eliminate well logging random noise and come
Reflection interference from bed boundary;
S4, according to going directly, the filtered dipole acoustic log data are normalized in wave amplitude, using uncommon
Scattered wave energy envelope before depth intervals pressure break described in your Bert transformation calculations;
S5, the post-fracturing scattered wave energy envelope of depth intervals is obtained using step S2 to the identical method of step S4;
S6, the filtered dipole acoustic log data are handled using waveform coherent superposition method, obtains pressure break forward and backwardly
The shear wave velocity curve of layer;
S7, the scattered wave energy envelope forward and backward according to the depth intervals pressure break, calculate a depth of the depth intervals
The discrepancy delta E of the forward and backward scattered wave energy envelope of pressure break at degree point;According to the shear wave velocity curve on the forward and backward stratum of the pressure break, meter
Calculate the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break at one depth point;
S8, it is dissipated according to the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break at one depth point and the pressure break are forward and backward
The discrepancy delta E of ejected wave energy envelope judges the pit shaft direction pressure break state and radial pressure break state of current depth point;
S9, the pressure break state that next depth point is judged using step S7 and the identical method of step S8, it is whole until traversing
A depth intervals;Different depth point pit shaft direction and radial pressure break state are counted, determines the pressure break result in entire depth section.
Optionally, the detailed process of the step S3 are as follows:
S31, it is by depth intervals inner orientationDipole acoustic log dataBandpass filtering is carried out, well logging is eliminated
Random noise obtains filtered dipole acoustic log data v (z, t), whereinValue interval is [0,360 °], and z is depth,
T is the time;
S32, dipole acoustic log data v (z, t) is transformed into frequency wavenumber domain, is eliminated using F-K filtering and comes from stratum boundary
The reflection interference in face obtains filtered dipole acoustic log data w (z, t) by two-dimension fourier inverse transformation.
Optionally, in the step S4, the basis goes directly wave amplitude to the filtered dipole acoustic log data
The formula being normalized are as follows:
G (z, t)=w (z, t)/w0 (1)
Wherein w0For wave amplitude of going directly, w (z, t) is filtered dipole acoustic log data, and g (z, t) is at normalization
Manage result.
It is optionally, described to utilize waveform coherent superposition method formula in the step S6 are as follows:
Wherein, XmIt (t) is N number of acoustic logging instrument z=z0Device receives m-th of reception energy converter in transducer array, and d is
Acoustic logging instrument receives the interval between energy converter, and T is time window TwPosition, v is a certain velocity amplitude in speed interval;
To in entire acoustic waveform or waveform a certain period and given speed interval by formula (2) calculate two-dimensional correlation function
Corr (v, T), when correlation function takes maximum, corresponding v value, that is, found out shear wave velocity vs。
Optionally, in the step S7:
The discrepancy delta E of the forward and backward scattered wave energy envelope of pressure break at one depth point are as follows:
Δ E=B (z0,t)-A(z0,t) (3)
Wherein, B (z0, t) and it is post-fracturing scattered wave energy envelope, A (z0, t) be pressure break before scattered wave energy envelope, t
For the time;
The depth point z=z0Locate the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break are as follows:
Δ s=10^6/VAfter pressure break-10^6/VBefore pressure break (4)
Wherein, VAfter pressure breakFor post-fracturing shear wave velocity, VBefore pressure breakFor the shear wave velocity before pressure break.
Optionally, the specific judgment method of the step S8 are as follows:
If Δ s=0, and Δ E=0 then prove well Zhou Yanshi not by pressure break;If Δ s > 0, and Δ E=0, then prove
Rock along pit shaft direction is by pressure break, and rock radially is not by pressure break;If Δ s=0, and Δ E > 0, then prove along pit shaft
The rock in direction is not by pressure break, and rock radially is by pressure break;If Δ s > 0, and Δ E > 0, then prove along pit shaft direction
Rock is by pressure break, and rock radially is also by pressure break.
Second aspect of the present invention, provides a kind of device for evaluating stratum hydraulic fracturing effect, and described device includes:
Acquisition module: for carrying out array sonic log in depth intervals, it is forward and backward just to obtain pressure break in depth intervals
Hand over dipole acoustic logging time domain four componets data and instrumental azimuth curve;
Conversion module: the four componets data and instrumental azimuth curve are transformed into terrestrial coordinate system by instrument coordinates system
Under, construct the dipole acoustic log data of different direction;
Filter module: being filtered the dipole acoustic log data of the different direction, eliminates well logging and makes an uproar at random
Sound and reflection interference from bed boundary;
Computing module: being normalized the filtered dipole acoustic log data according to through wave amplitude,
Scattered wave energy envelope, post-fracturing scattering wave energy packet before calculating the depth intervals pressure break using Hilbert transform
Network;The filtered dipole acoustic log data are handled using waveform coherent superposition method, obtain the shear wave on the forward and backward stratum of pressure break
Rate curve;
Judgment module: according to the scattered wave energy envelope that the depth intervals pressure break is forward and backward, the depth intervals are calculated
The discrepancy delta E of the forward and backward scattered wave energy envelope of pressure break at one depth point;According to the shear wave velocity on the forward and backward stratum of the pressure break
Curve calculates the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break at one depth point;It is pressed according at one depth point
The discrepancy delta E for splitting the forward and backward scattered wave energy envelope of pressure break described in the discrepancy delta s of forward and backward acoustic wave's slowness judges current depth point
Pit shaft direction pressure break state and radial pressure break state;
Statistical module: using judging that identical method judges the pressure break state of next depth point, until traversing entirely deeply
Spend section;Different depth point pit shaft direction and radial pressure break state are counted, determines the pressure break result in entire depth section.
Optionally, the filter module specifically includes:
Bandpass filtering unit: it is by depth intervals inner orientationDipole acoustic log dataBandpass filtering is carried out,
Well logging random noise is eliminated, filtered dipole acoustic log data v (z, t) is obtained, whereinValue interval is [0,360 °],
Z is depth, and t is the time;
F-K filter unit: transforming to frequency wavenumber domain for dipole acoustic log data v (z, t), is filtered and is eliminated using F-K
Reflection interference from bed boundary obtains filtered dipole acoustic log data w (z, t) by two-dimension fourier inverse transformation.
Optionally, in the computing module, the basis goes directly wave amplitude to the filtered dipole acoustic log number
According to the formula being normalized are as follows:
G (z, t)=w (z, t)/w0
Wherein w0For wave amplitude of going directly, w (z, t) is filtered dipole acoustic log data, and g (z, t) is at normalization
Manage result.
Optionally, in the judgment module, the difference according to the forward and backward acoustic wave's slowness of pressure break at one depth point
The discrepancy delta E of different Δ s and the forward and backward scattered wave energy envelope of the pressure break judge current depth point pit shaft direction pressure break state and
Radial pressure break state specifically:
If Δ s=0, and Δ E=0 then prove well Zhou Yanshi not by pressure break;If Δ s > 0, and Δ E=0, then prove
Rock along pit shaft direction is by pressure break, and rock radially is not by pressure break;If Δ s=0, and Δ E > 0, then prove along pit shaft
The rock in direction is not by pressure break, and rock radially is by pressure break;If Δ s > 0, and Δ E > 0, then prove along pit shaft direction
Rock is by pressure break, and rock radially is also by pressure break.
The beneficial effects of the present invention are:
1. the present invention carries out fracturing effect and comments using the forward and backward shear wave velocity variation of pressure break and scattered wave energy variation simultaneously
Valence, when well week, there are still can be with effective evaluation fracturing effect when netted seam;
2. compared with the routine evaluations method that can only evaluate fracturing height, present invention utilizes the scattered wave after direct wave into
Row fracture evaluation can be evaluated for fracturing effect and be provided with the fracturing effect of rock within the scope of tens of meters of effective evaluation well week
Reliable methods and techniques support.
Detailed description of the invention
It, below will be to needed in the technology of the present invention description in order to illustrate more clearly of technical solution of the present invention
Attached drawing is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for ability
For the those of ordinary skill of domain, without any creative labor, it can also be obtained according to these attached drawings others
Attached drawing.
Fig. 1 is the method flow schematic diagram of evaluation stratum hydraulic fracturing effect provided by the invention;
Fig. 2 is the apparatus structure schematic diagram of evaluation stratum hydraulic fracturing effect provided by the invention;
Fig. 3 is hydraulic fracturing experiments result figure.
Specific embodiment
The present invention proposes that a kind of difference using the forward and backward dipole acoustic log scattering wave energy of pressure break realizes that well week pressure break is imitated
The method of fruit evaluation, compares existing evaluation of Fracturing Effect on Compact Sandstone technology, and this method advantage is that the pressure along pit shaft direction can be evaluated
Effect is split, can also realize the evaluation of Fracturing Effect on Compact Sandstone within the scope of ten meters of number radial direction of well week.
In order to make the invention's purpose, features and advantages of the invention more obvious and easy to understand, below in conjunction with the present invention
Attached drawing in embodiment, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that disclosed below
Embodiment be only a part of the embodiment of the present invention, and not all embodiment.Based on the embodiments of the present invention, this field
Those of ordinary skill's all other embodiment obtained without making creative work, belongs to protection of the present invention
Range.
Referring to Fig. 1, proposing a kind of method for evaluating stratum hydraulic fracturing effect, which comprises
S1, array sonic log is carried out in depth intervals, obtain the forward and backward cross-dipole sound of pressure break in depth intervals
Wave well logging time domain four componets data and instrumental azimuth curve;Specifically, the four componets data can be denoted as respectively XX, XY, YX,
YY。
S2, by before the pressure break four componets data and instrumental azimuth curve terrestrial coordinate system is transformed by instrument coordinates system
Under, construct the dipole acoustic log data of different direction;
S3, the dipole acoustic log data of the different direction are filtered, eliminate well logging random noise and come
Reflection interference from bed boundary;
The detailed process of the step S3 are as follows:
S31, it is by depth intervals inner orientationDipole acoustic log dataBandpass filtering is carried out, well logging is eliminated
Random noise obtains filtered dipole acoustic log data v (z, t), whereinValue interval is [0,360 °], and z is depth,
T is the time;
S32, dipole acoustic log data v (z, t) is transformed into frequency wavenumber domain, is eliminated using F-K filtering and comes from stratum boundary
The reflection interference in face obtains filtered dipole acoustic log data w (z, t) by two-dimension fourier inverse transformation.
S4, according to going directly, the filtered dipole acoustic log data are normalized in wave amplitude, using uncommon
Scattered wave energy envelope before depth intervals pressure break described in your Bert transformation calculations;
In the step S4, the basis goes directly wave amplitude to the filtered dipole acoustic log data progress normalizing
Change the formula of processing are as follows:
G (z, t)=w (z, t)/w0 (1)
Wherein w0For wave amplitude of going directly, w (z, t) is filtered dipole acoustic log data, and g (z, t) is at normalization
Manage result.
S5, the post-fracturing scattered wave energy envelope of depth intervals is obtained using step S2 to the identical method of step S4;
S6, the filtered dipole acoustic log data are handled using waveform coherent superposition method, obtains pressure break forward and backwardly
The shear wave velocity curve of layer;
It is described to utilize waveform coherent superposition method formula are as follows:
Wherein, XmIt (t) is that N number of acoustic logging instrument receives m-th of reception energy converter in transducer array, d is sound wave
Logger receives the interval between energy converter, and T is time window TwPosition, v is a certain velocity amplitude in speed interval;To whole
A certain period and given speed interval in a acoustic waveform or waveform calculate two-dimensional correlation function Corr by formula (2)
(v, T), when correlation function takes maximum, corresponding v value, that is, found out shear wave velocity vs。
S7, the scattered wave energy envelope forward and backward according to the depth intervals pressure break, calculate a depth of the depth intervals
Spend point z=z0Locate the discrepancy delta E of the forward and backward scattered wave energy envelope of pressure break;It is bent according to the shear wave velocity on the forward and backward stratum of the pressure break
Line calculates the depth point z=z0Locate the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break;
One depth point z=z0Locate the discrepancy delta E of the forward and backward scattered wave energy envelope of pressure break are as follows:
Δ E=B (z0,t)-A(z0,t) (3)
Wherein, B (z0, t) and it is post-fracturing energy envelope, A (z0, t) be pressure break before energy envelope, t is the time;
The depth point z=z0Locate the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break are as follows:
Δ s=10^6/VAfter pressure break-10^6/VBefore pressure break (4)
Wherein, VAfter pressure breakFor post-fracturing shear wave velocity, VBefore pressure breakFor the shear wave velocity before pressure break.
S8, according to the depth point z=z0Discrepancy delta s and the pressure break for locating the forward and backward acoustic wave's slowness of pressure break are forward and backward scattered
The discrepancy delta E of ejected wave energy envelope judges the pit shaft direction pressure break state and radial pressure break state of current depth point;When Δ s > 0
When, the sampling interval of acoustic logging instrument is denoted as to the fracturing fracture height h of one depth point;
Specific judgment method are as follows:
If Δ s=0, and Δ E=0 then prove well Zhou Yanshi not by pressure break;If Δ s > 0, and Δ E=0, then prove
Rock along pit shaft direction is by pressure break, and rock radially is not by pressure break;If Δ s=0, and Δ E > 0, then prove along pit shaft
The rock in direction is not by pressure break, and rock radially is by pressure break;If Δ s > 0, and Δ E > 0, then prove along pit shaft direction
Rock is by pressure break, and rock radially is also by pressure break.
Rock near the rock, that is, borehole wall along pit shaft direction, the rock for the separate borehole wall that rock radially is.
Fracturing effect of the present invention according to the forward and backward shear wave velocity difference of pressure break and scattered wave capacity volume variance association evaluation well Zhou Yanshi, root
The fracturing effect that pit shaft direction rock is judged according to Δ s judges radial fracturing effect according to Δ E.
S9, the pressure break state that next depth point is judged using step S7 and the identical method of step S8, it is whole until traversing
A depth intervals;Different depth point pit shaft direction and radial pressure break are counted as a result, obtaining the fracturing height in entire depth section
With radial pressure break extension width.
Specifically, there are the numbers of the depth point of fracturing fracture height in statistics pit shaft direction, so that it may obtain the depth area
Between pit shaft direction fracturing fracture height.If statistics is n, depth there are the number of the depth location of fracturing fracture height
The fracturing fracture height in section is n × h.Different depth point radial direction pressure break state is counted, according to scattered wave earth-layer propagation speed v1
With time t1Conversion obtains radial pressure break extension width d=v1×t1/ 2, since subsequent scattered wave is utilized, belong to backwards to scattered
It penetrates, the distance that scattered wave is walked is twice of actual distance.
Referring to Fig. 2, the present invention provides a kind of device for evaluating stratum hydraulic fracturing effect, described device includes:
Acquisition module 210: for carrying out array sonic log in depth intervals, it is forward and backward to obtain pressure break in depth intervals
Cross-dipole acoustic logging time domain four componets data and instrumental azimuth curve;
Conversion module 220: the four componets data and instrumental azimuth curve are transformed into the earth by instrument coordinates system and sat
Under mark system, the dipole acoustic log data of different direction are constructed;
Filter module 230: being filtered the dipole acoustic log data of the different direction, and it is random to eliminate well logging
Noise and reflection interference from bed boundary;
The filter module 230 specifically includes:
Bandpass filtering unit: it is by depth intervals inner orientationDipole acoustic log dataCarry out band logical filter
Wave eliminates well logging random noise, obtains filtered dipole acoustic log data v (z, t), whereinValue interval be [0,
360 °], z is depth, and t is the time;
F-K filter unit: transforming to frequency wavenumber domain for dipole acoustic log data v (z, t), is filtered and is eliminated using F-K
Reflection interference from bed boundary obtains filtered dipole acoustic log data w (z, t) by two-dimension fourier inverse transformation.
Computing module 240: place is normalized to the filtered dipole acoustic log data according to through wave amplitude
Reason, scattered wave energy envelope, post-fracturing energy envelope before the depth intervals pressure break is calculated using Hilbert transform;Benefit
The filtered dipole acoustic log data are handled with waveform coherent superposition method, obtain the shear wave velocity on the forward and backward stratum of pressure break
Curve;
In the computing module 240, the basis go directly wave amplitude to the filtered dipole acoustic log data into
The formula of row normalized are as follows:
G (z, t)=w (z, t)/w0
Wherein w0For wave amplitude of going directly, w (z, t) is filtered dipole acoustic log data, and g (z, t) is at normalization
Manage result.
Judgment module 250: according to the scattered wave energy envelope that the depth intervals pressure break is forward and backward, the depth area is calculated
Between a depth point z=z0Locate the discrepancy delta E of the forward and backward scattered wave energy envelope of pressure break;According to the forward and backward stratum of the pressure break
Shear wave velocity curve calculates the depth point z=z0Locate the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break;According to the depth point
Z=z0Locate the discrepancy delta E judgement of the discrepancy delta s and the forward and backward scattered wave energy envelope of the pressure break of the forward and backward acoustic wave's slowness of pressure break
The pit shaft direction pressure break state of current depth point and radial pressure break state;
It is described according to the depth point z=z in the judgment module 2500Locate the difference of the forward and backward acoustic wave's slowness of pressure break
The discrepancy delta E of Δ s and the forward and backward energy envelope of the pressure break judges the pit shaft direction pressure break state and radial pressure break of current depth point
State specifically:
If Δ s=0, and Δ E=0 then prove well Zhou Yanshi not by pressure break;If Δ s > 0, and Δ E=0, then prove
Rock along pit shaft direction is by pressure break, and rock radially is not by pressure break;If Δ s=0, and Δ E > 0, then prove along pit shaft
The rock in direction is not by pressure break, and rock radially is by pressure break;If Δ s > 0, and Δ E > 0, then prove along pit shaft direction
Rock is by pressure break, and rock radially is also by pressure break.
Statistical module 260: using judging that identical method judges the pressure break state of next depth point, until traversing entirely
Depth intervals;Different depth point pit shaft direction and radial pressure break state are counted, the pressure break result in entire depth section is obtained.
In order to verify the present invention in the reliability of fracturing effect evaluation, using the process flow of invention to certain mouthful of pressure break
Well is handled.Referring to Fig. 3, Fig. 3 gives fracturing effect evaluation example, give for first natural gamma (GR) and
Hole diameter (CAL) curve, second give the forward and backward shear-wave slowness curve of pressure break.Third road gives depth intervals and pressure break is high
Degree, the 4th gives the forward and backward scattered wave capacity volume variance of pressure break, color brighter (white) bigger, the face that represents capacity volume variance in figure
It is smaller that color darker (black) represents capacity volume variance.It can be seen from the figure that the shear-wave slowness after X957-X975 interval pressure is obvious
It greater than the shear-wave slowness before pressure break, may determine that rock is by pressure break near the interval borehole wall according to step S8, thus may determine that
Fracturing height out.It can be seen that in the section shown in white circle from scattered wave capacity volume variance, post-fracturing scattering wave energy
It is obvious big, it is possible thereby to judge that, in X962-X966 interval, radially extending for pressure break at least reaches 25 meters.
The present invention proposes a kind of difference using the forward and backward dipole acoustic log shear wave velocity difference of pressure break and scattering wave energy
Association evaluation well week fracturing effect, compare existing evaluation of Fracturing Effect on Compact Sandstone technology, this method advantage is to evaluate along pit shaft
The fracturing effect in direction can also realize the evaluation of Fracturing Effect on Compact Sandstone within the scope of tens of meters of well week.
The above, the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations, although referring to before
Stating embodiment, invention is explained in detail, those skilled in the art should understand that: it still can be to preceding
Technical solution documented by each embodiment is stated to modify or equivalent replacement of some of the technical features;And these
It modifies or replaces, the spirit and scope for technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution.
Claims (10)
1. a kind of method for evaluating stratum hydraulic fracturing effect, which is characterized in that the described method includes:
S1, array sonic log is carried out in depth intervals, obtain the forward and backward cross-dipole cement bond logging of pressure break in depth intervals
Well time domain four componets data and instrumental azimuth curve;
S2, by before the pressure break four componets data and instrumental azimuth curve be transformed under terrestrial coordinate system by instrument coordinates system,
Construct the dipole acoustic log data of different direction;
S3, the dipole acoustic log data of the different direction are filtered;
S4, according to going directly, the filtered dipole acoustic log data are normalized in wave amplitude, utilize Martin Hilb
Scattered wave energy envelope before depth intervals pressure break described in special transformation calculations;
S5, the post-fracturing scattered wave energy envelope of depth intervals is obtained using step S2 to the identical method of step S4;
S6, the filtered dipole acoustic log data are handled using waveform coherent superposition method, obtains the forward and backward stratum of pressure break
Shear wave velocity curve;
S7, the energy envelope forward and backward according to the depth intervals pressure break, calculate pressure break at a depth point of the depth intervals
The discrepancy delta E of forward and backward scattered wave energy envelope;According to the shear wave velocity curve on the forward and backward stratum of the pressure break, calculate one
The discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break at depth point;
S8, according to the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break at one depth point and the forward and backward scattered wave of the pressure break
The discrepancy delta E of energy envelope judges the pit shaft direction fracturing effect and radial fracturing effect of current depth point;It, will as Δ s > 0
The sampling interval of acoustic logging instrument is denoted as the fracturing fracture height of one depth point;
S9, the pressure break state that next depth point is judged using step S7 and the identical method of step S8 are entire deep until traversing
Spend section;Different depth point pit shaft direction and radial pressure break are counted as a result, determining the fracturing height and diameter in entire depth section
To pressure break extension width.
2. evaluating the method for stratum hydraulic fracturing effect according to claim 1, which is characterized in that the step S3's is specific
Process are as follows:
S31, it is by depth intervals inner orientationDipole acoustic log dataBandpass filtering is carried out, it is random to eliminate well logging
Noise obtains filtered dipole acoustic log data v (z, t), whereinValue interval is [0,360 °], and z is depth, and t is
Time;
S32, dipole acoustic log data v (z, t) is transformed into frequency wavenumber domain, is eliminated using F-K filtering from bed boundary
Reflection interference obtains filtered dipole acoustic log data w (z, t) by two-dimension fourier inverse transformation.
3. evaluating the method for stratum hydraulic fracturing effect according to claim 2, which is characterized in that in the step S4, institute
State the formula that the filtered dipole acoustic log data are normalized according to wave amplitude of going directly are as follows:
G (z, t)=w (z, t)/w0 (1)
Wherein w0For wave amplitude of going directly, w (z, t) is filtered dipole acoustic log data, and g (z, t) is normalized knot
Fruit.
4. evaluating the method for stratum hydraulic fracturing effect according to claim 2, which is characterized in that in the step S6, institute
It states and utilizes waveform coherent superposition method formula are as follows:
Wherein, XmIt (t) is that N number of acoustic logging instrument receives m-th of reception energy converter in transducer array, d is acoustic logging instrument
Device receives the interval between energy converter, and T is time window TwPosition, v is a certain velocity amplitude in speed interval;To entire sound wave
A certain period and given speed interval in waveform or waveform calculate two-dimensional correlation function Corr (v, T) by formula (2),
When correlation function takes maximum, corresponding v value, that is, found out shear wave velocity vs。
5. evaluating the method for stratum hydraulic fracturing effect according to claim 1, which is characterized in that in the step S7:
In a depth point z=z0Locate the discrepancy delta E of the forward and backward scattered wave energy envelope of pressure break are as follows:
Δ E=B (z0,t)-A(z0,t) (3)
Wherein, B (z0, t) and it is post-fracturing energy envelope, A (z0, t) be pressure break before energy envelope, t is the time;
In depth point z=z0Locate the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break are as follows:
Δ s=10^6/VAfter pressure break-10^6/VBefore pressure break (4)
Wherein, VAfter pressure breakFor post-fracturing shear wave velocity, VBefore pressure breakFor the shear wave velocity before pressure break.
6. evaluating the method for stratum hydraulic fracturing effect according to claim 1, which is characterized in that the step S8's is specific
Judgment method are as follows:
If Δ s=0, and Δ E=0 then prove well Zhou Yanshi not by pressure break;If Δ s > 0, and Δ E=0, then prove along well
The rock in cylinder direction is by pressure break, and rock radially is not by pressure break;If Δ s=0, and Δ E > 0, then prove along pit shaft direction
Rock not by pressure break, rock radially is by pressure break;If Δ s > 0, and Δ E > 0, then the rock along pit shaft direction is proved
Rock by pressure break, and radially is also by pressure break.
7. a kind of device for evaluating stratum hydraulic fracturing effect, which is characterized in that described device includes:
Acquisition module: for carrying out array sonic log in depth intervals, the forward and backward Orthogonal-even of pressure break in depth intervals is obtained
Extremely sub- acoustic logging time domain four componets data and instrumental azimuth curve;
Conversion module: the four componets data and instrumental azimuth curve are transformed under terrestrial coordinate system by instrument coordinates system,
Construct the dipole acoustic log data of different direction;
Filter module: being filtered the dipole acoustic log data of the different direction, eliminate well logging random noise and
Reflection interference from bed boundary;
Computing module: being normalized the filtered dipole acoustic log data according to through wave amplitude, utilizes
Hilbert transform calculates the scattered wave energy envelope before the depth intervals pressure break, post-fracturing scattered wave energy envelope;Benefit
The filtered dipole acoustic log data are handled with waveform coherent superposition method, obtain the shear wave velocity on the forward and backward stratum of pressure break
Curve;
Judgment module: according to the scattered wave energy envelope that the depth intervals pressure break is forward and backward, one of the depth intervals is calculated
The discrepancy delta E of the forward and backward scattered wave energy envelope of pressure break at depth point;According to the shear wave velocity curve on the forward and backward stratum of the pressure break,
Calculate the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break at one depth point;According to before pressure break at one depth point,
The discrepancy delta E of the discrepancy delta s of acoustic wave's slowness and the forward and backward scattered wave energy envelope of the pressure break judges the well of current depth point afterwards
Cylinder direction pressure break result and radial pressure break result;
Statistical module: it uses and judges that identical method judges the pressure break of next depth point as a result, until traversing entire depth area
Between;Different depth point pit shaft direction and radial pressure break are counted as a result, determining the fracturing height in entire depth section and radial pressure
Split extension width.
8. evaluating the device of stratum hydraulic fracturing effect according to claim 7, which is characterized in that the filter module is specific
Include:
Bandpass filtering unit: it is by depth intervals inner orientationDipole acoustic log dataBandpass filtering is carried out, is eliminated
Well logging random noise, obtains filtered dipole acoustic log data v (z, t), whereinValue interval is [0,360 °], and z is
Depth, t are the time;
F-K filter unit: transforming to frequency wavenumber domain for dipole acoustic log data v (z, t), is eliminated and is come from using F-K filtering
The reflection interference of bed boundary obtains filtered dipole acoustic log data w (z, t) by two-dimension fourier inverse transformation.
9. evaluating the device of stratum hydraulic fracturing effect according to claim 7, which is characterized in that in the computing module,
The formula that the filtered dipole acoustic log data are normalized in the through wave amplitude of the basis are as follows:
G (z, t)=w (z, t)/w0
Wherein w0For wave amplitude of going directly, w (z, t) is filtered dipole acoustic log data, and g (z, t) is normalized knot
Fruit.
10. evaluating the device of stratum hydraulic fracturing effect according to claim 7, which is characterized in that in the judgment module,
It is described according to the discrepancy delta s of the forward and backward acoustic wave's slowness of pressure break at one depth point and the forward and backward scattering wave energy of the pressure break
The discrepancy delta E of envelope judges the pit shaft direction pressure break state and radial pressure break state of current depth point specifically:
If Δ s=0, and Δ E=0 then prove well Zhou Yanshi not by pressure break;If Δ s > 0, and Δ E=0, then prove along well
The rock in cylinder direction is by pressure break, and rock radially is not by pressure break;If Δ s=0, and Δ E > 0, then prove along pit shaft direction
Rock not by pressure break, rock radially is by pressure break;If Δ s > 0, and Δ E > 0, then the rock along pit shaft direction is proved
Rock by pressure break, and radially is also by pressure break.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910644976.5A CN110529087B (en) | 2019-07-17 | 2019-07-17 | Method and device for evaluating hydraulic fracturing effect of stratum |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910644976.5A CN110529087B (en) | 2019-07-17 | 2019-07-17 | Method and device for evaluating hydraulic fracturing effect of stratum |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110529087A true CN110529087A (en) | 2019-12-03 |
CN110529087B CN110529087B (en) | 2021-11-05 |
Family
ID=68660682
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910644976.5A Active CN110529087B (en) | 2019-07-17 | 2019-07-17 | Method and device for evaluating hydraulic fracturing effect of stratum |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110529087B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111911146A (en) * | 2020-08-17 | 2020-11-10 | 中国石油大学(华东) | Reservoir fracturing effect evaluation method based on fracture intersection chart |
CN111927428A (en) * | 2020-07-27 | 2020-11-13 | 中海油田服务股份有限公司 | Method and device for determining reservoir transformation effect of synergistic perforation |
CN111980676A (en) * | 2020-09-15 | 2020-11-24 | 中国石油大学(华东) | Method for evaluating well cementation quality through array acoustic logging and processing device |
CN112859164A (en) * | 2021-01-07 | 2021-05-28 | 中国地质科学院 | Pre-judgment method and system for fracturing effect of tight sandstone gas shale |
CN112987094A (en) * | 2021-02-20 | 2021-06-18 | 北京科技大学 | Hydraulic fracturing influence range evaluation method based on coal seam wave velocity field test |
CN114320284A (en) * | 2021-12-29 | 2022-04-12 | 中国石油化工集团有限公司 | Method for evaluating stratum fracturing effect by using dipole acoustic wave time difference correlation matrix |
CN114562244A (en) * | 2020-11-27 | 2022-05-31 | 中国石油天然气股份有限公司 | Method, device, equipment and medium for determining fracturing fracture of reservoir with fracture-cavity beside well |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289881A1 (en) * | 2010-12-21 | 2013-10-31 | Bikash K. Sinha | Hydraulic Fracture Characterization Using Borehole Sonic Data |
CN103775048A (en) * | 2014-01-20 | 2014-05-07 | 中国海洋石油总公司 | Method for determining height of fracturing facture in offshore low-hole-seepage reservoir |
CN106226819A (en) * | 2016-07-11 | 2016-12-14 | 中国石油集团川庆钻探工程有限公司 | Well side fault deep transverse wave reflection imaging identification method |
CN106908837A (en) * | 2017-01-22 | 2017-06-30 | 长江大学 | A kind of fracturing fracture form and fracture height determine method |
US20180283153A1 (en) * | 2014-10-01 | 2018-10-04 | Baker Hughes, A Ge Company, Llc | Methods and materials for evaluating and improving the production of geo-specific shale reservoirs |
-
2019
- 2019-07-17 CN CN201910644976.5A patent/CN110529087B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130289881A1 (en) * | 2010-12-21 | 2013-10-31 | Bikash K. Sinha | Hydraulic Fracture Characterization Using Borehole Sonic Data |
CN103775048A (en) * | 2014-01-20 | 2014-05-07 | 中国海洋石油总公司 | Method for determining height of fracturing facture in offshore low-hole-seepage reservoir |
US20180283153A1 (en) * | 2014-10-01 | 2018-10-04 | Baker Hughes, A Ge Company, Llc | Methods and materials for evaluating and improving the production of geo-specific shale reservoirs |
CN106226819A (en) * | 2016-07-11 | 2016-12-14 | 中国石油集团川庆钻探工程有限公司 | Well side fault deep transverse wave reflection imaging identification method |
CN106908837A (en) * | 2017-01-22 | 2017-06-30 | 长江大学 | A kind of fracturing fracture form and fracture height determine method |
Non-Patent Citations (1)
Title |
---|
王猛 等: ""基于交叉偶极子声波和固井质量测井评价压裂缝方法"", 《石油地质与工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111927428A (en) * | 2020-07-27 | 2020-11-13 | 中海油田服务股份有限公司 | Method and device for determining reservoir transformation effect of synergistic perforation |
CN111927428B (en) * | 2020-07-27 | 2024-02-06 | 中海油田服务股份有限公司 | Method and device for determining reservoir transformation effect by synergistic perforation |
CN111911146A (en) * | 2020-08-17 | 2020-11-10 | 中国石油大学(华东) | Reservoir fracturing effect evaluation method based on fracture intersection chart |
CN111980676A (en) * | 2020-09-15 | 2020-11-24 | 中国石油大学(华东) | Method for evaluating well cementation quality through array acoustic logging and processing device |
CN111980676B (en) * | 2020-09-15 | 2023-11-03 | 中国石油大学(华东) | Method and processing device for evaluating well cementation quality by array acoustic logging |
CN114562244A (en) * | 2020-11-27 | 2022-05-31 | 中国石油天然气股份有限公司 | Method, device, equipment and medium for determining fracturing fracture of reservoir with fracture-cavity beside well |
CN112859164A (en) * | 2021-01-07 | 2021-05-28 | 中国地质科学院 | Pre-judgment method and system for fracturing effect of tight sandstone gas shale |
CN112987094A (en) * | 2021-02-20 | 2021-06-18 | 北京科技大学 | Hydraulic fracturing influence range evaluation method based on coal seam wave velocity field test |
CN112987094B (en) * | 2021-02-20 | 2022-01-11 | 北京科技大学 | Hydraulic fracturing influence range evaluation method based on coal seam wave velocity field test |
CN114320284A (en) * | 2021-12-29 | 2022-04-12 | 中国石油化工集团有限公司 | Method for evaluating stratum fracturing effect by using dipole acoustic wave time difference correlation matrix |
Also Published As
Publication number | Publication date |
---|---|
CN110529087B (en) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110529087A (en) | A kind of method and device for evaluating stratum hydraulic fracturing effect | |
CN104656142B (en) | One kind is using vertical seismic profiling (VSP) and the united seismic layer labeling method of well logging | |
US7660194B2 (en) | Microseismic fracture mapping using seismic source timing measurements for velocity calibration | |
CA2156772C (en) | Single well system for mapping sources of acoustic energy | |
US7872944B2 (en) | Method of monitoring microseismic events | |
CN101630016B (en) | Method for improving imaging quality of vertical seismic profile | |
CN104216008B (en) | Downhole fracturing microseismic event identification method | |
CN107544087B (en) | Method and device for measuring near-surface stratum quality factor | |
CN103592680B (en) | A kind of log data based on FORWARD AND INVERSE PROBLEMS and Depth Domain seismic profile synthetic method | |
CN105259581A (en) | Seismic data time-depth conversion method | |
CN104007462A (en) | Fracture prediction method based on attenuation anisotropy | |
CN1748156A (en) | Permanent downhole resonant source | |
CN104950327A (en) | Method for determining positions of geophones of ground microseismic observation system | |
JPS62269084A (en) | Seismic exporation method having extremely high resolution in horizontal well | |
CN103713312A (en) | Design method of virtual source earthquake observation system | |
CN1245637C (en) | Tube-wave soundex | |
CN105301658A (en) | Method for extracting near-surface deconvolution operator by utilizing record of microlog | |
CN106353807A (en) | Fracture identification method and device | |
CN101609166B (en) | Method for measuring near surface structure of water area | |
CN100552472C (en) | Utilize vertical seismic profiling (VSP) and micro logging to carry out the seismic signal compensation method | |
CN110579799A (en) | Seismic acquisition observation method and system with equal travel time intervals | |
Sicking et al. | Forecasting reservoir performance by mapping seismic emissions | |
CN110579798B (en) | Seismic acquisition observation method and system with equal reflection angle intervals | |
CN106154321A (en) | The detection method of perforation signal and device | |
Murer et al. | Why monitoring with a single downhole microseismic array may not be enough: a case for multiwell monitoring of cyclic steam in diatomite |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |