CN105370274B - Down-hole formation porosity determines method - Google Patents
Down-hole formation porosity determines method Download PDFInfo
- Publication number
- CN105370274B CN105370274B CN201510930725.5A CN201510930725A CN105370274B CN 105370274 B CN105370274 B CN 105370274B CN 201510930725 A CN201510930725 A CN 201510930725A CN 105370274 B CN105370274 B CN 105370274B
- Authority
- CN
- China
- Prior art keywords
- stratum
- slowness
- formation
- wave
- shear
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 38
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 46
- 239000011707 mineral Substances 0.000 claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 29
- 239000011435 rock Substances 0.000 claims abstract description 29
- 239000011148 porous material Substances 0.000 abstract description 5
- 238000005755 formation reaction Methods 0.000 description 54
- 235000010755 mineral Nutrition 0.000 description 35
- 239000011575 calcium Substances 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 238000012937 correction Methods 0.000 description 5
- 229910052710 silicon Inorganic materials 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 229910052791 calcium Inorganic materials 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000004611 spectroscopical analysis Methods 0.000 description 4
- 230000005251 gamma ray Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 2
- UQSXHKLRYXJYBZ-UHFFFAOYSA-N Iron oxide Chemical compound [Fe]=O UQSXHKLRYXJYBZ-UHFFFAOYSA-N 0.000 description 2
- WCUXLLCKKVVCTQ-UHFFFAOYSA-M Potassium chloride Chemical compound [Cl-].[K+] WCUXLLCKKVVCTQ-UHFFFAOYSA-M 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 238000005056 compaction Methods 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000011777 magnesium Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000010453 quartz Substances 0.000 description 2
- 239000010703 silicon Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 229910052688 Gadolinium Inorganic materials 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- FYYHWMGAXLPEAU-UHFFFAOYSA-N Magnesium Chemical compound [Mg] FYYHWMGAXLPEAU-UHFFFAOYSA-N 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 239000004411 aluminium Substances 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- DLHONNLASJQAHX-UHFFFAOYSA-N aluminum;potassium;oxygen(2-);silicon(4+) Chemical compound [O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[O-2].[Al+3].[Si+4].[Si+4].[Si+4].[K+] DLHONNLASJQAHX-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229910001919 chlorite Inorganic materials 0.000 description 1
- 229910052619 chlorite group Inorganic materials 0.000 description 1
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000001739 density measurement Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000010433 feldspar Substances 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Substances [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 229910052749 magnesium Inorganic materials 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000010445 mica Substances 0.000 description 1
- 229910052618 mica group Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 229910052700 potassium Inorganic materials 0.000 description 1
- 239000001103 potassium chloride Substances 0.000 description 1
- 235000011164 potassium chloride Nutrition 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
- 239000011800 void material Substances 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A kind of down-hole formation porosity of present invention offer determines method, accumulates content according to stratum minerals, determines the rock matrix compressional slowness of underground set depth;According to array sonic log data, the stratum compressional slowness and formation shear slowness of the underground set depth are extracted;Finally according to the rock matrix compressional slowness, stratum compressional slowness and formation shear slowness of the underground set depth, the formation porosity of the underground set depth is determined.Method is determined using down-hole formation porosity provided in an embodiment of the present invention, the porosity value on stratum is determined using stratum compressional slowness and formation shear slowness simultaneously, the porosity of down-hole formation can be accurately obtained, without being corrected to obtained formation pore angle value, the utilization rate for improving Array Sonic Logging Waveformss improves the accuracy of formation pore angle value.
Description
Technical field
The present invention relates to oil and gas exploration technical fields more particularly to a kind of down-hole formation porosity to determine method.
Background technology
Porosity is the ratio of formation void volume and the total volume of rock, it is an important ginseng for describing formation physical property
Number, the accurate formation porosity that obtains is a background task for carrying out underground mine exploration.It can be right using earthquake, logging technique
Subsurface formations porosity is detected.Logging technique mainly carries out the spy of physical signal using various instruments in boring the well opened
It surveys, then underground geological condition and wellbore project situation is detected.
In the prior art, a kind of common porosity determines that method is:Utilize laboratory core velocity of longitudinal wave-porosity, close
Empirical equation is established in degree-porosity experiment, then passes through the physics such as the acoustic speed of logger in-site measurement underground and density
Measurement result is substituted into empirical equation and is calculated, determines the porosity of underground by signal.
But this porosity determines method in the prior art for use, needs to carry out lithology to the porosity value of calculating
Correction and compaction correction etc., since correcting value changes with the variation of formation condition, increase the difficulty of field application.
Invention content
A kind of down-hole formation porosity of offer of the embodiment of the present invention determines method, for solving porosity in the prior art
It determines that method needs to carry out lithology correction and compaction correction to the porosity value of calculating, increases asking for the difficulty of field application
Topic.
Down-hole formation porosity provided in an embodiment of the present invention determines method, including:
Content is accumulated according to stratum minerals, determines the rock matrix compressional slowness of underground set depth;
According to array sonic log data, stratum compressional slowness and the formation shear for extracting the underground set depth are slow
Degree;
According to the rock matrix compressional slowness, stratum compressional slowness and formation shear slowness of the underground set depth, really
The formation porosity of the fixed underground set depth.
In another embodiment, rock matrix compressional slowness, the stratum compressional slowness according to the underground set depth
And formation shear slowness, determine that the formation porosity of the underground set depth includes:
According to
Determine the formation porosity of the underground set depth;
Wherein,Indicate porosity, SmaIndicate rock matrix compressional slowness, SsFor formation shear slowness, SpIndicate stratum
Compressional slowness.
It is described to accumulate content according to stratum minerals in another embodiment, determine the rock matrix longitudinal wave of underground set depth
Before slowness, further include:
Detect the elemental abundance of the stratum element of the underground set depth;
The volume content of subterranean minerals corresponding with the stratum element is determined according to the elemental abundance of the stratum element.
In another embodiment, it is described determined according to the elemental abundance of the stratum element it is corresponding with the stratum element
Layer mineral volume content include:
According to
[M]=[C]-1[E]
Determine the volume content of subterranean minerals corresponding with the element;
Wherein, [M] is the column matrix of the subterranean minerals volume content;[E] is the elemental abundance of the stratum element
Column matrix;[C] is that known stratum element abundance is carried out multinomial recurrence to the contribution of the volume content of corresponding subterranean minerals
Obtained coefficient square formation.
It is described to accumulate content according to stratum minerals in another embodiment, determine the rock matrix longitudinal wave of underground set depth
Slowness includes:
According to
Wherein, viIndicate the volume content of i-th kind of mineral in stratum,For the compressional slowness of i-th kind of mineral in stratum.
It is described according to array sonic log data in another embodiment, extract the stratum longitudinal wave of the underground set depth
Slowness and formation shear slowness include:
According to
Calculate that multiple receivers receive by same emission source the moment t waveform signals sent out related coefficient;Institute
Emission source and the multiple receiver are stated along arranged in a straight line, the spacing between each two adjacent receiver is equal;Wherein, R tables
Show the related coefficient, TWIndicate the time window length of setting;S indicates slowness;Window is received in first receiver when τ is described
Waveform on the windowing moment, first receiver be adjacent with the emission source;fmIndicate what m-th of receiver received
Full waveform signal, at the time of t indicates any point on the full waveform signal, Δ z indicate each two adjacent receiver it
Between spacing;
Determine that the maximum value of the related coefficient obtained in the longitudinal wave range of waveforms of the full waveform signal is corresponding slow
Degree is the compressional slowness on the stratum;
Determine that the maximum value of the related coefficient obtained in the shear wave range of waveforms of the full waveform signal is corresponding slow
Degree is the shear-wave slowness on the stratum.
In another embodiment, the time window length T of the settingWFor three wave periods of the longitudinal wave waveform;Alternatively, described
The time window length T of settingWFor three wave periods of the shear wave waveform.
In another embodiment,
The maximum value for the related coefficient that the determination obtains in the longitudinal wave range of waveforms of the full waveform signal corresponds to
Slowness be the stratum compressional slowness during, the search range of the compressional slowness is 40-140 microseconds/foot;Or
Person,
The maximum value for the related coefficient that the determination obtains in the shear wave range of waveforms of the full waveform signal corresponds to
Slowness be the stratum shear-wave slowness during, the search range of the shear-wave slowness is 60-180 microseconds/foot.
In another embodiment, the waveform sent out in moment t by same emission source for calculating multiple receivers and receiving
Before the related coefficient of signal, further include:
According to the energy of m-th short window in the full waveform signal and first short window to the long window of the m short windows
Normalized energy value ratio maximum value, determine the first arrival time of the longitudinal wave in the full waveform signal.
It is described to be arrived with first short window according to the energy of m-th short window in the full waveform signal in another embodiment
The maximum value R of the ratio of the normalized energy value of the long window of the m short windowsm, determine the longitudinal wave in the full waveform signal
First arrival time include:Determine in the full waveform signal energy of m-th short window and first short window to the m short windows
Long window normalized energy value ratio be 0.2RmaxWhen, the initial time of the corresponding m-th short window is the longitudinal wave
First arrival time.
Down-hole formation porosity provided by the invention determines method, and content is accumulated according to stratum minerals, determines that underground is set
The rock matrix compressional slowness of depth;According to array sonic log data, the stratum longitudinal wave for extracting the underground set depth is slow
Degree and formation shear slowness;Finally according to the rock matrix compressional slowness, stratum compressional slowness and ground of the underground set depth
Layer shear-wave slowness, determines the formation porosity of the underground set depth.Using down-hole formation hole provided in an embodiment of the present invention
Porosity determines method, while the porosity value on stratum is determined using stratum compressional slowness and formation shear slowness, can be accurate
Ground obtains the porosity of down-hole formation, without being corrected to obtained formation pore angle value, improves array sonic log money
The utilization rate of material improves the accuracy of formation pore angle value.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to do one simply to introduce, it should be apparent that, the accompanying drawings in the following description is this hair
Some bright embodiments for those of ordinary skill in the art without having to pay creative labor, can be with
Obtain other attached drawings according to these attached drawings.
Fig. 1 is the flow diagram that down-hole formation porosity of the embodiment of the present invention determines method;
Fig. 2 is the principle signal of the related coefficient for the waveform signal that the embodiment of the present invention determines that multiple receivers receive
Figure;
Fig. 3 is the waveform diagram of full-wave train in the embodiment of the present invention;
Fig. 4 is the schematic diagram of the long short time-window of full waveform in the embodiment of the present invention;
Fig. 5 is the underground porosity and corresponding depth rock core that down-hole formation porosity of the embodiment of the present invention determines that method calculates
Test the effect contrast figure of the porosity obtained.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
The every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
A kind of down-hole formation porosity of present invention offer determines method, and the well logging needed for porosity logging monitoring acquires
The problems such as series is more, and the later stage corrects, the spy detected in well using Element Capture Spectroscopy Logging and two kinds of instruments of array sonic log
Reference number establishes one kind formation porosity easy to operate, with high accuracy and determines method, not only avoids asking for later stage correction
Topic, and improve the utilization rate of Array Sonic Logging Waveformss.
Fig. 1 is the flow diagram that down-hole formation porosity of the embodiment of the present invention determines method.Referring to Fig. 1, of the invention
The down-hole formation porosity that embodiment provides determines that method includes:
S101:Content is accumulated according to stratum minerals, determines the rock matrix compressional slowness of underground set depth;
It is described to accumulate content according to stratum minerals, before the rock matrix compressional slowness for determining underground set depth, also wrap
It includes:
Detect the elemental abundance of the stratum element of the underground set depth;
The volume content of subterranean minerals corresponding with the stratum element is determined according to the elemental abundance of the stratum element.
Specifically, the content for various elements in stratum being detected by Element Capture Spectroscopy Logging instrument realizes subterranean minerals component
Analysis.The Element Capture Spectroscopy Logging instrument not only can be in bore hole borehole measurement, but also can be in casing borehole measurement.The instrument knot
Structure is made of neutron source, crystal counter, photomultiplier, high voltage amplifier electronic circuit, can be in fresh water mud, saturated salt cement
It is used under the conditions of slurry or oil-base mud, potassium chloride mud and gassiness mud etc..In well logging, sent out to stratum by neutron source
Fast neutron is penetrated, stratum element is induced and inelastic scattering reaction occurs, while releasing gamma ray.After Multiple Scattering, in
Son slows down and forms thermal neutron.Thermal neutron is captured, and generates the feature capture gamma ray of element, element is by discharging gamma ray
Return to original state.These inelastic scattering Gamma-ray spectroscopies and capture gamma spectra are detected and recorded with crystal counter.It utilizes
The inelastic gamma spectrum that detector detects, the content that can obtain the elements such as C, O, Si, Ca is handled by spectrum unscrambling;And to wherein
Main capture gamma spectrums can obtain the constituent contents such as Si, Ca, S, Fe, Ti and Gd by spectrum unscrambling processing.According to these elements
Content can obtain the percentage composition that mineral are corresponded in stratum using specific oxides closure modelling technique.Research card
Real, formation elements type is various, but Relatively centralized is in a few element such as oxygen, silicon, aluminium, and these elements composition is each
Kind mineral are even more to concentrate on a few in sedimentary rock, wherein quartz 31.5%, carbonate rock 20%, mica and chlorite
19%, calcedony 9%, feldspar 7.5%, clay mineral 7.5%, iron oxide 3%.It is first in mineral when mineral chemical composition is stablized
Plain degree is held essentially constant.Based on this, Element Capture Spectroscopy Logging is by selected mineral indicator element, by constituent content
It is converted into sedimentary mineral content, meets Comprehensive Evaluation of Reservoir needs.Specifically, the only a few member of mineral characteristic can be characterized by choosing
Indicator element of the element as the mineral.For example, silicon (Si) can be used as the indicator element of quartz, calcium (Ca) can be used as the finger of limestone
Show that element, magnesium (Mg) and calcium (Ca) can be used as dolomitic indicator element.For complicated mineral, then 2-3 kind members can be chosen
Element is used as indicator element, as potassium feldspar can choose Si, Al and K as indicator element, and can by K contents number and Erie
Stone distinguishes.
It is in close relations between mineral and element, it mutually restricts, by between the suitable elemental abundance of determination and mineral content
Relationship, elemental abundance can be converted into mineral content.Specifically, the element according to the stratum element is rich
Degree determines that the volume content of subterranean minerals corresponding with the stratum element includes:
According to
[M]=[C]-1[E] formula (1)
Determine the volume content of subterranean minerals corresponding with the element;
Wherein, [M] is the column matrix of the subterranean minerals volume content;[E] is the elemental abundance of the stratum element
Column matrix;[C] is that known stratum element abundance is carried out multinomial recurrence to the contribution of the volume content of corresponding subterranean minerals
Obtained coefficient square formation.
It is described to accumulate content according to stratum minerals, determine that the rock matrix compressional slowness of underground set depth includes:
According to
Wherein, viIndicate the volume content of i-th kind of mineral in stratum,For the compressional slowness of i-th kind of mineral in stratum,
Unit is microsecond/foot.
S102:According to array sonic log data, stratum compressional slowness and the stratum for extracting the underground set depth are horizontal
Wave slowness;
Specifically, described according to array sonic log data, extract the underground set depth stratum compressional slowness and
Formation shear slowness includes:
According to
Calculate that multiple receivers receive by same emission source the moment t waveform signals sent out related coefficient;Institute
Emission source and the multiple receiver are stated along arranged in a straight line, the spacing between each two adjacent receiver is equal;Wherein, R tables
Show the related coefficient, TWIndicate the time window length of setting;S indicates slowness;Window is received in first receiver when τ is described
Waveform on the windowing moment, first receiver be adjacent with the emission source;fmIndicate what m-th of receiver received
Full waveform signal, at the time of t indicates any point on the full waveform signal, Δ z indicate each two adjacent receiver it
Between spacing;
Determine that the maximum value of the related coefficient obtained in the longitudinal wave range of waveforms of the full waveform signal is corresponding slow
Degree is the compressional slowness on the stratum;
Determine that the maximum value of the related coefficient obtained in the shear wave range of waveforms of the full waveform signal is corresponding slow
Degree is the shear-wave slowness on the stratum.
If array acoustic imaging logging instrument has 6 receivers.The transmitting of the receiver and the setting depth of stratum
Source is equidistantly arranged along straight line, the source between each receiver and emission source away from respectively z0, z1 ... z6.It is each two adjacent
The spacing of receiving transducer is Δ z (such as 6 inches), then away from that can be write as, (m is whole more than or equal to 1 in the source of m-th of receiving transducer
Number):
zm=z0+(m-1)Δz
If the corresponding full waveform signal received of each receiving transducer is respectively f1(t),f2(t),…,f6(t).For
Discrete time series is then at the time of any point i on waveform:
ti=t0+(i-1)ΔT
Wherein, △ T indicate the time interval of sampling, and unit is microsecond;I indicates that sampling number, i are whole more than or equal to 1
Number.
Fig. 2 is the principle signal of the related coefficient for the waveform signal that the embodiment of the present invention determines that multiple receivers receive
Figure.Referring to Fig. 2, the full waveform signal that the same emission source that R1-R6, which is respectively 6 receivers, to be received is sent out.In Fig. 2
The slope of long side of dotted line parallelogram indicate the slowness S, the time window length T of the expression setting of short sideW.Short side
Windowing time instant τ of the window on the waveform that first receiver the receives when position intersected with the waveform of R1 indicates described.Formula 3
It indicates that a plurality of waveform for being included to dotted line parallelogram in Fig. 2 is added, and normalizes.When the phase of these waveforms
Unanimously when (that is, the slowness S is exactly the propagation slowness of corresponding waveform), waveform is strengthened, and as a result (is known as related coefficient)
Close to 1;And when phase is inconsistent, waveform is weakened, and to which amplitude is smaller, is as a result much smaller than 1.Therefore, coefficient R
Value range is 0≤R≤1.As R=0, without any relationship between expression waveform;As R=1, the complete phase of waveform morphology is indicated
Together.
Specifically, by formula 3 from full waveform extracting data longitudinal wave, shear-wave slowness when, need to choose first appropriate
Time window length TW, and determine in rational slowness search range the value of slowness S.When the slowness S value and stratum longitudinal wave or
When the slowness of person's formation shear is close, coefficient R is maximum.That is, the corresponding slowness S of related coefficient maximum value is should
The estimated value of wavelet (stratum longitudinal wave or formation shear) slowness obtains stratum compressional slowness or formation shear slowness.
Optionally, described to set convenient for determining stratum compressional slowness and shear-wave slowness as early as possible in order to improve system operations efficiency
Fixed time window length TWFor three wave periods of the longitudinal wave waveform;Alternatively, the time window length T of the settingWFor the shear wave
Three wave periods of shape.
Optionally, in order to improve the accuracy of stratum compressional slowness and shear-wave slowness, the determination is in the full waveform
The corresponding slowness of maximum value of the related coefficient obtained in the longitudinal wave range of waveforms of signal is the mistake of the compressional slowness on the stratum
The search range of Cheng Zhong, the compressional slowness are 40-140 microseconds/foot;
Alternatively,
The maximum value for the related coefficient that the determination obtains in the shear wave range of waveforms of the full waveform signal corresponds to
Slowness be the stratum shear-wave slowness during, the search range of the shear-wave slowness is 60-180 microseconds/foot.
S103:It is slow according to the rock matrix compressional slowness, stratum compressional slowness and formation shear of the underground set depth
Degree, determines the formation porosity of the underground set depth.
Specifically, the rock matrix compressional slowness, stratum compressional slowness and stratum according to the underground set depth
Shear-wave slowness determines that the formation porosity of the underground set depth includes:
According to
Determine the formation porosity of the underground set depth;
Wherein,Indicate porosity, SmaIndicate rock matrix compressional slowness, SsFor formation shear slowness, SpIndicate stratum
Compressional slowness.
Down-hole formation porosity provided in an embodiment of the present invention determines method, accumulates content according to stratum minerals, determines well
The rock matrix compressional slowness of lower set depth;According to array sonic log data, the stratum of the underground set depth is extracted
Compressional slowness and formation shear slowness;It is finally slow according to the rock matrix compressional slowness of the underground set depth, stratum longitudinal wave
Degree and formation shear slowness, determine the formation porosity of the underground set depth.Using underground provided in an embodiment of the present invention
Formation porosity determines method, while the porosity value on stratum, energy are determined using stratum compressional slowness and formation shear slowness
It is enough accurately obtained the porosity of down-hole formation, without being corrected to obtained formation pore angle value, improves array acoustic
The utilization rate of well-log information improves the accuracy of bottom porosity value.
Further, the waveform signal sent out in moment t by same emission source for calculating multiple receivers and receiving
Related coefficient before, further include:
According to the energy of m-th short window in the full waveform signal and first short window to the long window of the m short windows
Normalized energy value ratio maximum value, determine the first arrival time of the longitudinal wave in the full waveform signal.
Specifically, it is determined that the first arrival time of the longitudinal wave in the full waveform signal may be used but be not limited to following two
Realization method:
In the first realization method, the first arrival time of longitudinal wave can be determined by quotient of amplitudes method.
Irregular due to before longitudinal wave reaches, receiving noise information, Amplitude Ratio longitudinal wave amplitude is much smaller, and
Inside noise signal or inside longitudinal wave signal, their amplitude variations are little, and are made an uproar in longitudinal wave wave to front and back amplitude is put
Amplitude variations are very big between sound and longitudinal wave, therefore can detect discrete wave by the maximum of postpeak in full waveform and leading peak
The longitudinal wave wave of shape is to point.
In second of realization method, can by it is short long when window energy ratio function determine the first arrival time of longitudinal wave.
Fig. 3 is the waveform diagram of full-wave train in the embodiment of the present invention.Referring to Fig. 3, full-wave train include noise, longitudinal wave,
Shear wave and Stoneley wave.Fig. 4 is the schematic diagram of the long short time-window of full waveform in the embodiment of the present invention.Referring to Fig. 4, setting with n-th
A sampling point is that the m-th short window energy of window point is:
Wherein, i indicates that sampling point, i are the integer more than or equal to 1;L indicates the time window length of the short window, AiDescribed in expression
The corresponding waveforms amplitude of i-th of sampling point.
The normalized value of long window energy when the short window of m-th of stepping is
The normalized ratio of long window energy when short with m-th of stepping window of the m-th short window energy is
R is found out according to formula 7mMaximum value Rmax, determining to work as Rm=0.2RmaxWhen, the corresponding m-th short window rises
Begin the first arrival time that the moment is the longitudinal wave.Specifically, the first arrival time of the longitudinal wave is m × △ T.Wherein, △ T are sampling
Time interval.In order to save the time, it can estimate that longitudinal wave first arrival and window are long in advance, estimated value can be selected according to wave shape playback figure,
And ensure first arrival time in estimated range.
Fig. 5 is the underground porosity and corresponding depth rock core that down-hole formation porosity of the embodiment of the present invention determines that method calculates
Test the effect contrast figure of the porosity obtained.Referring to Fig. 5, wherein in rightmost porosity road, dotted line and it is labeled as arrow
It is core experiment porosity, solid wire is that well logging calculates porosity, and the two is numerically substantially close to while size trend one
It causes.Illustrate to determine that method can preferably be detected formation porosity using down-hole formation porosity of the embodiment of the present invention.
Finally it should be noted that:The above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations;To the greatest extent
Present invention has been described in detail with reference to the aforementioned embodiments for pipe, it will be understood by those of ordinary skill in the art that:Its according to
So can with technical scheme described in the above embodiments is modified, either to which part or all technical features into
Row equivalent replacement;And these modifications or replacements, various embodiments of the present invention technology that it does not separate the essence of the corresponding technical solution
The range of scheme.
Claims (9)
1. a kind of down-hole formation porosity determines method, which is characterized in that including:
Content is accumulated according to stratum minerals, determines the rock matrix compressional slowness of underground set depth;
According to array sonic log data, the stratum compressional slowness and formation shear slowness of the underground set depth are extracted;
According to the rock matrix compressional slowness, stratum compressional slowness and formation shear slowness of the underground set depth, institute is determined
State the formation porosity of underground set depth;
The rock matrix compressional slowness, stratum compressional slowness and formation shear slowness according to the underground set depth, really
The formation porosity of the underground set depth includes calmly:
According to
Determine the formation porosity of the underground set depth;
Wherein,Indicate porosity, SmaIndicate rock matrix compressional slowness, SsFor formation shear slowness, SpIndicate that stratum longitudinal wave is slow
Degree.
2. according to the method described in claim 1, it is characterized in that, it is described according to stratum minerals accumulate content, determine that underground is set
Before the rock matrix compressional slowness of depthkeeping degree, further include:
Detect the elemental abundance of the stratum element of the underground set depth;
The volume content of subterranean minerals corresponding with the stratum element is determined according to the elemental abundance of the stratum element.
3. according to the method described in claim 2, it is characterized in that, it is described according to the elemental abundance of the stratum element determine with
The volume content of the corresponding subterranean minerals of the stratum element includes:
According to
[M]=[C]-1[E]
Determine the volume content of subterranean minerals corresponding with the element;
Wherein, [M] is the column matrix of the subterranean minerals volume content;[E] is the row square of the elemental abundance of the stratum element
Battle array;[C] is that known stratum element abundance is carried out multinomial recurrence to the contribution of the volume content of corresponding subterranean minerals to obtain
Coefficient square formation.
4. method according to any one of claims 1 to 3, which is characterized in that it is described to accumulate content according to stratum minerals, really
The rock matrix compressional slowness for determining underground set depth includes:
According to
Wherein, viIndicate the volume content of i-th kind of mineral in stratum,For the compressional slowness of i-th kind of mineral in stratum.
5. according to the method described in claim 1, it is characterized in that, described according to array sonic log data, the well is extracted
The stratum compressional slowness and formation shear slowness of lower set depth include:
According to
Calculate that multiple receivers receive by same emission source the moment t waveform signals sent out related coefficient;The hair
Source and the multiple receiver are penetrated along arranged in a straight line, the spacing between each two adjacent receiver is equal;Wherein, R indicates institute
State related coefficient, TWIndicate the time window length of setting;S indicates slowness;The wave that window is received in first receiver when τ is described
Windowing moment in shape, first receiver are adjacent with the emission source;fmIndicate the all-wave that m-th of receiver receives
Waveform signal, at the time of t indicates any point on the full waveform signal, △ z are indicated between each two adjacent receiver
Spacing;
Determine that the corresponding slowness of maximum value of the related coefficient obtained in the longitudinal wave range of waveforms of the full waveform signal is
The compressional slowness on the stratum;
Determine that the corresponding slowness of maximum value of the related coefficient obtained in the shear wave range of waveforms of the full waveform signal is
The shear-wave slowness on the stratum.
6. according to the method described in claim 5, it is characterized in that, the time window length T of the settingWFor the longitudinal wave waveform
Three wave periods;Alternatively, the time window length T of the settingWFor three wave periods of the shear wave waveform.
7. according to the method described in claim 5, it is characterized in that,
The maximum value for the related coefficient that the determination obtains in the longitudinal wave range of waveforms of the full waveform signal is corresponding slow
During degree is the compressional slowness on the stratum, the search range of the compressional slowness is 40-140 microseconds/foot;Alternatively,
The maximum value for the related coefficient that the determination obtains in the shear wave range of waveforms of the full waveform signal is corresponding slow
During degree is the shear-wave slowness on the stratum, the search range of the shear-wave slowness is 60-180 microseconds/foot.
8. according to claim 5-7 any one of them methods, which is characterized in that it is described calculate that multiple receivers receive by
Same emission source further includes before the related coefficient of the moment t waveform signals sent out:
According to returning for the long window of the energy of m-th short window in the full waveform signal and first short window to the m short windows
One changes the maximum value of the ratio of energy value, determines the first arrival time of the longitudinal wave in the full waveform signal.
9. according to the method described in claim 8, it is characterized in that, described according to m-th short window in the full waveform signal
Energy and first short window to the ratio of the normalized energy value of the long window of the m short windows maximum value Rmax, described in determination
The first arrival time of longitudinal wave in full waveform signal includes:Determine the energy of m-th short window and in the full waveform signal
The ratio of one short window to the normalized energy value of the long window of the m short windows is 0.2RmaxWhen, the corresponding m-th short window
Initial time be the longitudinal wave first arrival time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510930725.5A CN105370274B (en) | 2015-12-14 | 2015-12-14 | Down-hole formation porosity determines method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510930725.5A CN105370274B (en) | 2015-12-14 | 2015-12-14 | Down-hole formation porosity determines method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105370274A CN105370274A (en) | 2016-03-02 |
CN105370274B true CN105370274B (en) | 2018-08-21 |
Family
ID=55372832
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510930725.5A Active CN105370274B (en) | 2015-12-14 | 2015-12-14 | Down-hole formation porosity determines method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105370274B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105866003B (en) * | 2016-05-16 | 2019-05-17 | 中国地质大学(武汉) | A kind of shale gas reservoir total porosity acquisition methods and device |
CN107795320B (en) * | 2016-09-05 | 2021-01-05 | 中国石油化工股份有限公司 | Calculation method for horizontal well carbonate reservoir parameters |
CN108375793A (en) * | 2018-03-01 | 2018-08-07 | 中国石油化工股份有限公司江汉油田分公司物探研究院 | Comprehensive normal compaction speed preparation method and device and pore pressure preparation method and device |
CN115248463B (en) * | 2022-05-24 | 2024-05-24 | 中国石油大学(华东) | Correction method for inelastic scattering influence of D-T source neutron porosity logging |
CN116626779B (en) * | 2023-07-21 | 2023-11-10 | 中海石油(中国)有限公司深圳分公司 | Method and system for calculating porosity of subsurface mountain stratum |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4692910A (en) * | 1985-03-11 | 1987-09-08 | Amoco Corporation | Methods of determining lithological characteristics of an underground formation which utilize compressional velocity and shear velocity data |
CN1015671B (en) * | 1990-06-09 | 1992-02-26 | 潘仁芳 | Oil deposit-described method by earthquake date |
US8630146B2 (en) * | 2009-04-02 | 2014-01-14 | Baker Hughes Incorporated | Method and apparatus for estimating formation permeability and electroacoustic constant of an electrolyte-saturated multi-layered rock taking into account osmosis |
CN103064118A (en) * | 2013-01-08 | 2013-04-24 | 西南石油大学 | Method of acoustic logging and quantifying cavern filling degree |
CN104570127B (en) * | 2013-10-09 | 2017-06-20 | 中国石油化工股份有限公司 | A kind of method of utilization seimic wave velocity Simultaneous Inversion porosity and shale content |
CN104865614B (en) * | 2014-02-20 | 2017-09-29 | 中国石油化工股份有限公司 | Complicated reservoirs Fluid Identification Method based on variable matrix parameter |
CN103984009B (en) * | 2014-04-16 | 2016-04-13 | 刘致水 | Based on three kinds of porosity type quantitative inversion methods of P-wave And S well-log information |
CN104950331B (en) * | 2015-06-30 | 2017-08-11 | 中国石油天然气集团公司 | A kind of porosity of sand mud reservoir and the earthquake prediction method of shale content |
CN105093331B (en) * | 2015-07-13 | 2018-07-10 | 中国石油集团东方地球物理勘探有限责任公司 | The method for obtaining Rock Matrix bulk modulus |
-
2015
- 2015-12-14 CN CN201510930725.5A patent/CN105370274B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN105370274A (en) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105370274B (en) | Down-hole formation porosity determines method | |
CN104453873B (en) | Shale oil-gas economical and effective layer section evaluation method | |
US8573298B2 (en) | Method for petrophysical evaluation of shale gas reservoirs | |
Cannon | Petrophysics: a practical guide | |
US6714480B2 (en) | Determination of anisotropic moduli of earth formations | |
McNally | The prediction of geotechnical rock properties from sonic and neutron logs | |
Franco et al. | Sonic investigation in and around the borehole | |
US20120163123A1 (en) | Stress in formations from azimuthal variation in acoustic and other properties | |
Harrison et al. | Acquisition and analysis of sonic waveforms from a borehole monopole and dipole source for the determination of compressional and shear speeds and their relation to rock mechanical properties and surface seismic data | |
Butel et al. | Prediction of in situ rock strength using sonic velocity | |
Wissler et al. | Calibration of the Early Cretaceous time scale: a combined chemostratigraphic and cyclostratigraphic approach to the Barremian-Aptian interval, Campania Apennines and southern Alps (Italy) | |
CN104484573B (en) | A kind of method of determination stratum stiffness coefficient | |
Timur et al. | Downhole geophysical logging | |
Franquet et al. | Integrated acoustic, mineralogy, and geomechanics characterization of the Huron shale, Southern West Virginia, USA | |
Kowalski | Formation strength parameters from well logs | |
Russell | Seismic reservoir characterization and pre-stack inversion in resource shale plays | |
Wong et al. | Geophysical well logs from a shallow test well and automatic determination of formation velocities from full-waveform sonic logs | |
Opuwari | Petrophysical evaluation of the Albian age gas bearing sandstone reservoirs of the OM field, Orange basin, South Africa | |
Bustos et al. | Cased Hole Solution with Fast Neutrons and Sonic Reflective Waves in Tight Reservoirs | |
Reischman et al. | Horizontal Cased Hole Evaluation Using a New Pulsed Neutron Spectroscopy Tool and Dipole Sonic | |
Musolino et al. | Evaluating uncertainty in empirically derived unconfined compressive strength (UCS) estimates and implications for drilling applications; a case study from the Cooper Basin | |
US11795813B2 (en) | Modified sonic tool for advanced formation evaluation | |
Mwenifumbo et al. | Borehole geophysical logging in the Flin Flon mining camp | |
Mari et al. | Active and Passive acoustic logging applied to the detection of preferential flow in a sedimentary aquifer | |
Murray et al. | Deployment of new sonic-while-drilling technology in a large development project in the Middle East: a case study |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Tang Jun Inventor after: Zhang Chengguang Inventor after: Zheng Gongming Inventor before: Tang Jun |
|
COR | Change of bibliographic data | ||
GR01 | Patent grant | ||
GR01 | Patent grant |