CN101285381B - Process for inversing soft ground horizontal wave velocity by leaky mode waves - Google Patents
Process for inversing soft ground horizontal wave velocity by leaky mode waves Download PDFInfo
- Publication number
- CN101285381B CN101285381B CN200710090386XA CN200710090386A CN101285381B CN 101285381 B CN101285381 B CN 101285381B CN 200710090386X A CN200710090386X A CN 200710090386XA CN 200710090386 A CN200710090386 A CN 200710090386A CN 101285381 B CN101285381 B CN 101285381B
- Authority
- CN
- China
- Prior art keywords
- velocity
- wave
- stratum
- dipole
- shear 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for inversion of speed of soft formation transverse waves by leaky mode waves, belonging to a logging technology applying physical geophysical sound waves. The specific steps are as follows: inputting wave shape of monopole or dipole leaky mode waves of a soft formation well section with serious diameter extension and conventional well measurement data of corresponding well section formation; extracting a phase velocity dispersion curve of the leaky mode waves from array wave shape and determining the main energy frequency range of the leaky mode waves; at various speeds of formation transverse waves, calculating variance between the leaky mode waves speed calculated by theoretical arithmetic and the leaky mode waves phase velocity extracted from field array wave shape; determining the corresponding formation transverse wave speed in a state in which the variance is zero so as to carry out successive well section processing and output an inversion transverse speed result map. The invention solves the problem of the prior multipole subarray acoustic logging apparatus that the formation transverse wave speed can not be directly measured in the soft formation well sections with ultra-large well diameters, loose superficial layers or serious diameter extension. The method expands the application fields of the borehole acoustic logging technology in oil exploration process and oil development process, and improves the application value of the borehole acoustic logging technology.
Description
Technical field
The invention belongs to the applied geophysics acoustic logging, when particularly relating to the serious soft formation exploration of super large hole diameter, loose shallow-layer or hole enlargement, the method for utilizing the leakage mode ripple of frequency dispersion to come inverting stratum shear wave velocity.
Background technology
Development along with oil exploration and development technique, the equal an urgent demand acoustic logging in field such as seismic prospecting and drilling engineering provides stratum shear wave velocity section accurately, as the inverting of earthquake elastic wave impedance, the processing of earthquake AVO attribute, drilling stratum pressure prediction and borehole wall stability evaluation.At present, the wellhole acoustic logging provides the effective means of continuous stratum shear wave velocity section.But the stratum shear wave velocity less than well in the soft formation of the mud velocity of sound, the symmetrical monopole sound source of traditional acoustic logging can not excite the shear wave on stratum, but the Stoneley wave inverting stratum shear wave velocity that utilizes one pole to excite.In recent years, developed asymmetrical sound source of the dipole, realized measuring the dipole acoustic logging technique of any stratum shear wave velocity by exciting wellhole beam mode ripple.During actual the well logging, the excitation intensity of one pole Stoneley wave and DIPOLE FLEXURAL WAVES and Dispersion all are subjected to parameter influences such as formation properties and hole diameter serious, even at big hole diameter, loose shallow-layer or the more serious soft formation of hole enlargement, doublet source often can not excite bending wave, the one pole sound source can not excite Stoneley wave, shown in 430～435 meters and 446～452 meters well sections among Fig. 1.Therefore, at big hole diameter, loose shallow-layer or the more serious soft formation of hole enlargement,, all can't guarantee to provide the shear wave velocity on stratum from the angle of existing dipole acoustic log commercial measurement and the angle of one pole Stoneley wave inverting.
From a large amount of on-the-spot multipoles (one pole and dipole) array sonic log data, find the shear wave velocity that the stratum can't be directly measured in big hole diameter, loose shallow-layer or the more serious soft formation dipole well logging of hole enlargement.Set up wellhole multipole array acoustic log theoretical model, the one pole of numerical simulation study soft formation, dipole acoustic log response characteristic show, serious soft formation one pole, the dipole of hole enlargement all excites the bigger leakage mode ripple of amplitude in theory, and the one pole Stoneley wave and the DIPOLE FLEXURAL WAVES amplitude that are used to measure shear wave velocity are very little, even complete obiteration, as shown in Figures 2 and 3.The leakage mode ripple of soft formation one pole, dipole excitation is a kind of frequency dispersion mode wave, and promptly phase velocity changes with frequency change, as shown in Figure 4.The Dispersion of leakage mode ripple mainly is sensitive to the mud velocity of sound and stratum shear wave velocity in stratum velocity of longitudinal wave, the well, and as shown in Figure 5, wherein, this mode wave reaches 10% to the sensitivity of stratum shear wave velocity.Therefore, in the soft formation situation that one pole Stoneley wave and DIPOLE FLEXURAL WAVES all do not excite, can consider shear wave velocity with the Dispersion inverting soft formation of the leakage mode ripple of relatively stable one pole or dipole.
Summary of the invention
The object of the present invention is to provide a kind of method of inversing soft ground horizontal wave velocity by leaky mode waves, the technical problem that the leakage mode ripple that utilizes frequency dispersion with solution carries out inverting to the shear wave velocity of the more serious soft formation of big hole diameter, loose shallow-layer or hole enlargement.
The objective of the invention is Windows environment, utilize Microsoft VC++ SDK, realize in the following way based on the microcomputer more than the internal memory 256M:
The Wave data of the leakage mode ripple of being measured by multipole array sonic log instrument one pole or dipole mode of operation of at first, that hole enlargement is a serious depth point of soft formation well section and stratum velocity of longitudinal wave, hole diameter, the mud velocity of sound, the mud density parameter that this corresponding depth point is measured by conventional logging input to several memories of computer.
Secondly, utilize the Matrix Pencil method of frequency spectrum weighting correlation method, Prony method, correction or the array waveform dispersion curve extractive technique of homomorphism processing method, monopole array waveform or dipole array waveform extracting by soft formation go out leakage mode phase velocity of wave dispersion curve, and the main energy frequency range of definite one pole leakage mode ripple or dipole leakage mode ripple.
Once more, liquid-filled well bore sound field theory according to people such as Kurkjian proposition, set up wellhole multipole array acoustic log model, stratum velocity of longitudinal wave in the model, density, hole diameter, the mud velocity of sound and mud density are measured known by traditional logging method, change the stratum shear wave velocity with iterative method or parabolic method dispersion equation D (f by theory, k, V
s, V
c, V
f, ρ, ρ
f, R)=0 solve soft formation one pole or dipole leakage mode phase velocity of wave dispersion curve with frequency change, in the formula, k=
2π f/V
pBe wave number, f-frequency, V
p-phase velocity, V
s-stratum shear wave velocity, V
c-stratum velocity of longitudinal wave, V
fThe mud velocity of sound in the-well, ρ-density of earth formations, ρ
f-mud density, R-well radius.
Then, in one pole or the main energy frequency range of dipole leakage mode ripple, calculate the variance sum of the leakage mode phase velocity of wave of theoretical down leakage mode phase velocity of wave of calculating of Different Strata shear wave velocity and the extraction of on-the-spot array waveform information, obtain the variance curve that changes with the stratum shear wave velocity.Utilize Lagrangian method to the variance curve precise interpolation, the stratum shear wave velocity of correspondence is the stratum shear wave velocity that this depth point inverting obtains when determining variance and be zero.
At last, utilize above-mentioned inverting flow process that the serious soft formation well section of hole enlargement is handled continuously, obtain the stratum shear wave velocity result map of soft formation well section, result map is presented on the display, or print and publish picture.
The present invention has overcome existing multipole array sonic log instrument can not directly be measured the stratum shear wave velocity in the serious loose soft formation well section of hole enlargement shortcoming, utilize the leakage mode ripple of frequency dispersion in one pole or the dipole array waveform to come inverting stratum shear wave velocity, this method has been widened the application of wellhole acoustic logging in oil exploration and development process, has improved its using value.
Description of drawings
Fig. 1 is the waveform information and the contrast schematic diagram of loose shallow-layer well section one pole and dipole acoustic log.
Fig. 2 is the array time domain waveform schematic diagram of the theoretical monopole acoustic logging of synthesizing of soft formation.Fig. 3 is the array time domain waveform schematic diagram that the theoretical synthetic dipole acoustic of soft formation is logged well.
Fig. 4 is the leakage mode phase velocity of wave dispersion curve figure that soft formation monopole sound source excites.
Fig. 5 is the sensitivity curves figure of soft formation monopole leakage mode ripple to the mud velocity of sound in stratum shear wave velocity, velocity of longitudinal wave and the well.
Fig. 6 is the FB(flow block) according to the leakage mode ripple inverting stratum shear wave velocity method of a kind of frequency dispersion provided by the present invention.
Fig. 7 is the monopole array time domain waveform schematic diagram of the theoretical soft formation that synthesizes.
Fig. 8 is the spectrum diagram of soft formation monopole array waveform correspondence.
Fig. 9 utilizes the frequency spectrum weighted, coherent figure of the leakage mode ripple that frequency spectrum weighting correlation method obtains from monopole array waveform and the dispersion curve (shown in the heavy line) that extracts.
Figure 10 is a wellhole monopole acoustic logging theoretical model schematic diagram.
Figure 11 is leakage mode phase velocity of wave dispersion curve (small circle) and theoretical leakage mode phase velocity of wave dispersion curve (solid line) figure that calculates of adjustment stratum shear wave velocity that extracts from the array waveform according to the method for the invention.
Figure 12 be according to the method for the invention in the main energy frequency range of one pole leakage mode ripple, the Different Strata shear wave velocity is the leakage mode wave dispersion curve and the variance sum of extracting dispersion curve of calculating down, is finally inversed by the schematic diagram of soft formation shear wave velocity when being equalled zero by both variances.
Figure 13 is the stratum shear wave velocity result map that the one pole leakage mode ripple inverting of on-the-spot shallow-layer well section (velocity of longitudinal wave of stratum about 2500 meter per seconds, stratum shear wave velocity about 1200) multipole array sonic log data is drawn according to the method for the invention.
The specific embodiment
The present invention is as follows about the description of the drawings:
Fig. 1 has reflected the waveform information and the contrast of loose shallow-layer well section one pole and dipole acoustic log, and this well section is a soft formation well section, and the 1st road provides natural gamma, bit size and CAL among the figure, and the 2nd, 3 roads provide dipole and one pole waveform respectively.Can find out from the one pole and the dipole waveform of 430～435 meters and 446～452 meters two well sections:, do not excited one pole Stoneley wave and DIPOLE FLEXURAL WAVES because hole enlargement seriously causes one pole and dipole acoustic log data to have only a kind of leakage mode ripple.
Fig. 2 has reflected the array time domain waveform of the theoretical monopole acoustic logging of synthesizing of soft formation, and as can be seen from Figure, this soft formation situation one pole sound source has only excited the leakage mode ripple, does not have Stoneley wave.
Fig. 3 has reflected the array time domain waveform that the theoretical synthetic dipole acoustic of soft formation is logged well, and as can be seen from Figure, this soft formation situation doublet source has only excited the leakage mode ripple, does not have bending wave.
Fig. 4 has reflected the leakage mode phase velocity of wave dispersion curve that soft formation monopole sound source excites, and wherein abscissa is a frequency, and ordinate is a phase velocity.
Fig. 5 reflects the sensitivity curves of soft formation monopole leakage mode ripple to the mud velocity of sound in stratum shear wave velocity, velocity of longitudinal wave and the well, the Dispersion of leakage mode ripple mainly is sensitive to the stratum velocity of longitudinal wave in low-frequency range as can be seen from Figure, and high band mainly is sensitive to the mud velocity of sound in the well.Find simultaneously at the Dispersion of main this mode wave of energy frequency band the stratum shear wave velocity also to be compared sensitivity, sensitivity coefficient reaches about 10%.
Fig. 7 has reflected the monopole array time domain waveform of theoretical synthetic soft formation, wherein has only a kind of leakage mode ripple, does not have Stoneley wave.
Fig. 8 has reflected the frequency spectrum of soft formation monopole array waveform correspondence, and main concentration of energy is at 5 KHzs～10 KHzs.
Describe the present invention in detail below in conjunction with drawings and Examples.
At the serious soft formation of hole enlargement, one pole sound source and doublet source have all excited the leakage mode compressional wave of the very strong frequency dispersion of amplitude, and this mode wave is relatively stable, and are sensitive to very much the stratum shear wave velocity.Therefore, the present invention greater than 30 centimetres, the stratum shear wave velocity soft formation less than 1400 meter per seconds, can adopt one pole leakage mode ripple inverting stratum shear wave velocity and two kinds of specific embodiments of dipole leakage mode ripple inverting stratum shear wave velocity at hole diameter.
Scheme one, the concrete steps of one pole leakage mode ripple inverting stratum shear wave velocity are:
(1). waveform (as shown in Figure 7) data of the leakage mode ripple of being measured by multipole array sonic log instrument monopolar operation pattern of the depth point of soft formation well section that hole enlargement is serious and stratum velocity of longitudinal wave, hole diameter, the mud velocity of sound, the mud density parameter that this corresponding depth point is measured by conventional logging input to several memories of computer.
(2). utilize FFT that monopole array waveform is as shown in Figure 7 done spectrum analysis, obtain the frequency spectrum of per pass waveform, as shown in Figure 8.On this basis, utilize following formula to obtain frequency spectrum weighting correlation diagram, and extract one pole leakage mode phase velocity of wave dispersion curve by index of correlation maximum by the array waveform, as shown in Figure 9,
F is a frequency in the formula, k=2 π f/V
pBe wave number, V
pBe phase velocity; N is the number of receiver, and d is the spacing of adjacent two receivers; X
n(f) be the frequency spectrum of the waveform correspondence of n receiver, complex conjugate is got in " * " expression.
(3). set up wellhole monopole acoustic logging theoretical model, as shown in figure 10, derive theoretical dispersion equation D (f, k, the V of wellhole monopole mode wave according to the inside and outside wave equation of wellhole and fringe conditions
s, V
c, V
f, ρ, ρ
f, R)=0, wherein, k=2 π f/V
pBe wave number, f-frequency, V
p-phase velocity, V
s-stratum shear wave velocity, V
c-stratum velocity of longitudinal wave, V
fThe mud velocity of sound in the-well, ρ-density of earth formations, ρ
f-mud density, R-well radius.Stratum velocity of longitudinal wave, density, hole diameter, the mud velocity of sound and mud density are measured known by other logging methods in the model, change the stratum shear wave velocity and find the solution one pole leakage mode phase velocity of wave dispersion curve by dispersion equation, shown in solid line among Figure 11 with Newton iteration method.
(4). in the main energy frequency range of one pole leakage mode ripple, i.e. 5 KHzs～10 KHzs, calculate the leakage mode phase velocity of wave (solid line among Figure 11) of theoretical calculating under the Different Strata shear wave velocity and the variance sum of the leakage mode phase velocity of wave (small circle among Figure 11) that the array waveform information extracts, obtain variance curve, as shown in figure 12 with the variation of stratum shear wave velocity.
(5). by to the variance curve interpolation, determine the stratum shear wave velocity of correspondence when variance is zero, be the stratum shear wave velocity that this depth point inverting obtains, as shown in figure 12.This stratum shear wave velocity that obtains by the method for the invention is 507.32m/s, and actual stratum shear wave velocity is 508m/s.
Scheme two, the concrete steps of dipole leakage mode ripple inverting stratum shear wave velocity are:
(1). the Wave data of the leakage mode ripple of being measured by multipole array sonic log instrument dipole mode of operation of the depth point of soft formation well section that hole enlargement is serious and stratum velocity of longitudinal wave, hole diameter, the mud velocity of sound, the mud density parameter that this corresponding depth point is measured by conventional logging input to several memories of computer.
(2). utilize FFT electrode couple subarray waveform to do spectrum analysis, obtain the frequency spectrum of per pass waveform.On this basis, utilize frequency spectrum weighting correlation method to extract dipole leakage mode phase velocity of wave dispersion curve.
(3). set up wellhole dipole acoustic logging theory model, derive the theoretical dispersion equation of wellhole dipole mode wave according to the inside and outside wave equation of wellhole and fringe conditions.Stratum velocity of longitudinal wave, density, hole diameter, the mud velocity of sound and mud density are known by other logging methods measurements in the model, change the stratum shear wave velocity and find the solution dipole leakage mode phase velocity of wave dispersion curve with Newton iteration method by dispersion equation.
(4). in the main energy frequency range of dipole leakage mode ripple, i.e. 2 KHzs～5 KHzs, calculate the variance sum of the leakage mode phase velocity of wave of theoretical down leakage mode phase velocity of wave of calculating of Different Strata shear wave velocity and the extraction of on-the-spot dipole array waveform information, obtain the variance curve that changes with the stratum shear wave velocity.
(5). by to the variance curve interpolation, determine the stratum shear wave velocity of correspondence when variance is zero, be the stratum shear wave velocity that this depth point inverting obtains.
In sum, with the multipole array acoustic log data of the method for the invention processing shallow-layer well section, export continuous stratum shear wave velocity section, as shown in figure 13.The stratum velocity of longitudinal wave of this shallow-layer well section is about 2500 meter per seconds, and the stratum shear wave velocity belongs to typical soft formation about 1200.Shown in the 1st road bit size, hole diameter and gamma ray curve among the figure, this well uses 9.5 inches drill bits well-diggings serious hole enlargement phenomenon to occur 430～435 meters and 446～453 meters well sections, and hole diameter reaches 50 centimetres.The 2nd and 3 roads provide the one pole of multipole array sonic log data and the full-wave train waveform of dipole respectively among the figure, can find out from the full-wave train waveform, two well section one pole and the dipole serious in hole enlargement all excite the very strong leakage mode ripple of amplitude, are not used for one pole Stoneley wave and the DIPOLE FLEXURAL WAVES that soft formation is measured shear wave velocity and excite.Handle the dipole array waveform information with traditional slowness time matching method, the stratum shear wave velocity that obtains index of correlation plane and correspondence is shown in the 4th road among the figure.Find out from the index of correlation plan view, only the index of correlation 435～445 meters regular slim hole section DIPOLE FLEXURAL WAVES is bigger, can extract reliable stratum shear wave velocity by DIPOLE FLEXURAL WAVES, and very poor even do not have in the correlation of the serious well section DIPOLE FLEXURAL WAVES of other two hole enlargements, can't extract the stratum shear wave velocity this moment by the dipole data.The 5th road provides shear wave velocity, the contrast of stratum velocity of longitudinal wave that the stratum shear wave velocity that utilizes the inventive method to utilize the inverting of one pole leakage mode ripple Dispersion and regular well section dipole extract among the figure.Comparing result shows that the inventive method can obtain continuous stratum shear wave velocity section by the inverting of leakage mode ripple at the loose shallow-layer of hole condition complexity.In the soft formation well section of 435～445 meters hole diameter rules, the DIPOLE FLEXURAL WAVES waveform is better, identical fine of the shear wave velocity result (small circle) that the inversion result of the method for the invention (heavy line) and dipole extract; 430～435 meters and 446～453 meters serious well sections of hole enlargement, do not excited DIPOLE FLEXURAL WAVES, can't extract the stratum shear wave velocity from the dipole data, the method of the invention is stable, reliable by the stratum shear wave velocity that the one pole leakage mode ripple of frequency dispersion is finally inversed by, consistent with stratum velocity of longitudinal wave trend and the engineering that can be used for oil exploration and exploitation is used than rationally the ripple time difference in length and breadth.
1, frequency spectrum weighting correlation method list of references: Tang Xiaoming, Zheng Chuanhan work,, the acoustics of quantitatively logging well, petroleum industry publishing house, p52-54 in 2004.
2, Prony method list of references: Lang, S.W.Kurkjian, J.H. etc., 1987, Estimatingslowness dispersion from arrays of sonic logging waveforms, Geophysics, Vol.52, p530-544.
3, Matrix Pencil method list of references: Ekstrom, M.P., 1995, Dispersion estimationfrom borehole acoustic arrays using a modified matrix pencil algorithm, 29thAsilomar Conference on Signals Systems and Computers, p449-453.
4, homomorphism processing method list of references: Ellefsen, K.J.Cheng, C.H.Tubman, K.M., 1989, Estimating phase velocity and attenuation of guided waves inacoustic logging data, Geophysics, VOL, 54, p1054-1059.
5, wellhole sound field Modeling Theory list of references: Kurkjian, J.H.and Chang, S.K., 1986, Acoustic multipole sources in fuid-filled boreholes, Geophysics.51 (1) 148-163.
Claims (1)
1. the method for an inversing soft ground horizontal wave velocity by leaky mode waves, it is based on the Windows environment of the above microcomputer of internal memory 256M, utilizes Microsoft VC++ SDK, realizes in the following way:
The Wave data of the leakage mode ripple of being measured by multipole array sonic log instrument one pole or dipole mode of operation of at first, that hole enlargement is a serious depth point of soft formation well section and stratum velocity of longitudinal wave, hole diameter, the mud velocity of sound, the mud density parameter that this corresponding depth point is measured by conventional logging input to several memories of computer;
Secondly, utilize the Matrix Pencil method of frequency spectrum weighting correlation method, Prony method, correction or the array waveform dispersion curve extractive technique of homomorphism processing method, monopole array waveform or dipole array waveform extracting by soft formation go out leakage mode phase velocity of wave dispersion curve, and the main energy frequency range of definite one pole leakage mode ripple or dipole leakage mode ripple;
Once more, liquid-filled well bore sound field theory according to the Kurkjian proposition, set up wellhole multipole array acoustic log model, stratum velocity of longitudinal wave in the model, density, hole diameter, the mud velocity of sound and mud density are measured known by traditional logging method, change the stratum shear wave velocity with iterative method or parabolic method dispersion equation D (f by theory, k, V
s, V
c, V
f, ρ, ρ
f, R)=0 solve soft formation one pole or dipole leakage mode phase velocity of wave dispersion curve with frequency change, in the formula, k=2 π f/V
pBe wave number, f-frequency, V
p-phase velocity, V
s-stratum shear wave velocity, V
c-stratum velocity of longitudinal wave, V
fThe mud velocity of sound in the-well, ρ-density of earth formations, ρ
f-mud density, R-well radius;
Then, in one pole or the main energy frequency range of dipole leakage mode ripple, calculate the variance sum of the leakage mode phase velocity of wave of theoretical down leakage mode phase velocity of wave of calculating of Different Strata shear wave velocity and the extraction of on-the-spot array waveform information, obtain the variance curve that changes with the stratum shear wave velocity; Utilize Lagrangian method to the variance curve precise interpolation, the stratum shear wave velocity of correspondence is the stratum shear wave velocity that this depth point inverting obtains when determining variance and be zero;
At last, utilize above-mentioned inverting flow process that the serious soft formation well section of hole enlargement is handled continuously, obtain the stratum shear wave velocity result map of soft formation well section, result map is presented on the display, or print and publish picture.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200710090386XA CN101285381B (en) | 2007-04-09 | 2007-04-09 | Process for inversing soft ground horizontal wave velocity by leaky mode waves |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200710090386XA CN101285381B (en) | 2007-04-09 | 2007-04-09 | Process for inversing soft ground horizontal wave velocity by leaky mode waves |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101285381A CN101285381A (en) | 2008-10-15 |
| CN101285381B true CN101285381B (en) | 2011-06-15 |
Family
ID=40057763
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200710090386XA Expired - Fee Related CN101285381B (en) | 2007-04-09 | 2007-04-09 | Process for inversing soft ground horizontal wave velocity by leaky mode waves |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101285381B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103364822B (en) * | 2012-04-05 | 2016-08-24 | 核工业北京地质研究院 | A kind of layered approach being applicable to sedimentary basin shallow layer stratum |
| CN103233727B (en) * | 2013-05-13 | 2015-04-08 | 中国石油大学(华东) | Inversion method of stratum shear wave velocity radial sections |
| CN107165623B (en) * | 2017-06-07 | 2019-03-19 | 南方科技大学 | Monopole while-drilling acoustic logging instrument matched with bottom drilling tool combination for use and method for measuring slow formation transverse wave speed |
| CN109025979A (en) * | 2018-04-24 | 2018-12-18 | 长江大学 | A kind of forward simulation alteration stratum is with the method for boring multipole sound wave dispersion curve |
| CN110173259B (en) * | 2019-05-14 | 2022-08-30 | 中国石油化工集团有限公司 | Stoneley wave energy compensation correction method |
| CN112901156B (en) * | 2021-01-14 | 2022-11-04 | 思凡(上海)石油设备有限公司 | Rock debris logging detection system |
| CN114017014B (en) * | 2021-12-01 | 2023-10-24 | 北京华晖探测科技股份有限公司 | Acoustic logging wave acoustic method |
| CN116009085B (en) * | 2023-02-02 | 2024-03-12 | 哈尔滨工业大学 | Soft stratum transverse wave speed measurement method and device based on full waveform inversion |
| CN117270053B (en) * | 2023-09-28 | 2024-04-16 | 成都理工大学 | Stratum transverse wave slowness and attenuation calculation method based on dipole acoustic logging |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6831876B1 (en) * | 2003-07-09 | 2004-12-14 | Goodrich Corporation | Acoustic window |
| CN1601304A (en) * | 2004-10-26 | 2005-03-30 | 大庆油田有限责任公司 | Ground stress multifrequency reverse demonstration method of dipole trans verse wave well logging |
| CN1601303A (en) * | 2004-10-26 | 2005-03-30 | 大庆油田有限责任公司 | Selection method of sound source frequency in dipole transverse wave well logging of bushing well |
-
2007
- 2007-04-09 CN CN200710090386XA patent/CN101285381B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6831876B1 (en) * | 2003-07-09 | 2004-12-14 | Goodrich Corporation | Acoustic window |
| CN1601304A (en) * | 2004-10-26 | 2005-03-30 | 大庆油田有限责任公司 | Ground stress multifrequency reverse demonstration method of dipole trans verse wave well logging |
| CN1601303A (en) * | 2004-10-26 | 2005-03-30 | 大庆油田有限责任公司 | Selection method of sound source frequency in dipole transverse wave well logging of bushing well |
Non-Patent Citations (2)
| Title |
|---|
| 伍先运等.利用斯通利波速度反演地层横波速度的方法与应用研究.《地球物理学报》.1995,第38卷第216-223页. * |
| 赖富强等.利用多极子陈列声波测井预测地层破裂压力.《勘探地球物理进展》.2007,第30卷(第1期),第39-42页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101285381A (en) | 2008-10-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101285381B (en) | Process for inversing soft ground horizontal wave velocity by leaky mode waves | |
| CN103233727B (en) | Inversion method of stratum shear wave velocity radial sections | |
| US6098021A (en) | Estimating formation stress using borehole monopole and cross-dipole acoustic measurements: theory and method | |
| US7649804B2 (en) | In-situ determination of yield stress state of earth formations | |
| US8379483B2 (en) | Method and apparatus for determining radial shear velocity variation from dipole acoustic logging | |
| CN105277982B (en) | A kind of mud shale total content of organic carbon earthquake prediction method | |
| CN110348135B (en) | Method for evaluating stratum permeability by acoustic logging while drilling | |
| CN109563736A (en) | The estimation of horizontal stress and nonlinear constant in the anisotropic formation of interlayer carbonate horizon in such as organic shale reservoir | |
| CA2610455A1 (en) | Radial profiling of formation mobility using horizontal and vertical shear slowness profiles | |
| MXPA05013185A (en) | Methods and systems for acoustic waveform processing. | |
| Prioul et al. | Forward modeling of fracture-induced sonic anisotropy using a combination of borehole image and sonic logs | |
| Xu et al. | Estimation of elastic wave velocity from logging-while-drilling measurement in unconsolidated formation of accretionary wedge in nankai trough | |
| Lee et al. | Near-borehole characteristics of hydraulic fractures and fracturing-induced sonic-wave attenuation | |
| CN114017014A (en) | Acoustic wave logging wave acoustics method | |
| CN112034516B (en) | Method for evaluating three-dimensional anisotropy of well-cycle transverse wave | |
| CN114961700B (en) | Three-dimensional acoustic logging longitudinal wave velocity radial profile inversion imaging method | |
| Li et al. | A Novel Method for Evaluating Formation Fracturing Effect Utilizing Acoustic Logging | |
| RU2201606C1 (en) | Method of typification and correlation of oil and gas productive rocks by borehole spectral-time parameters | |
| Zheng | Dispersion Asymptotic Analysis and its Applications in Acoustic Logging | |
| Liu et al. | Inversion of gas hydrate saturation and solid frame permeability in a gas hydrate-bearing sediment by Stoneley wave attenuation | |
| Sinha et al. | Applications of sonics and ultrasonics in geophysical prospecting | |
| Maalouf et al. | Predicting the Elastic Properties of Anisotropic Reservoirs from Dipole Logging While Drilling Using Machine Learning | |
| Zheng et al. | Identifying stress-induced anisotropy and stress direction using cross-dipole acoustic logging | |
| Velez et al. | Core data integration/validation of sonic derived anisotropic mechanical properties to expedite well decisions in unconventional reservoirs | |
| Li et al. | Effect of perforation on borehole acoustic wave logging: Modeling and application |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110615 Termination date: 20120409 |