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 PDF

Info

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
formation
shear wave
leakage mode
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
Application number
CN200710090386XA
Other languages
Chinese (zh)
Other versions
CN101285381A (en
Inventor
苏远大
秦瑞宝
孙建孟
邓启才
陈雪莲
汤丽娜
庄春喜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum East China
CNOOC Research Center
Original Assignee
China University of Petroleum East China
CNOOC Research Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum East China, CNOOC Research Center filed Critical China University of Petroleum East China
Priority to CN200710090386XA priority Critical patent/CN101285381B/en
Publication of CN101285381A publication Critical patent/CN101285381A/en
Application granted granted Critical
Publication of CN101285381B publication Critical patent/CN101285381B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种泄漏模式波反演软地层横波速度的方法,属于应用地球物理声波测井技术。具体步骤如下:输入扩径严重的软地层井段的单极或偶极泄漏模式波波形和对应井段地层的常规测井资料,再由阵列波形中提取出泄漏模式波的相速度频散曲线,并确定泄漏模式波的主要能量频段,计算不同地层横波速度下的理论计算的与现场阵列波形提取的泄漏模式波相速度的方差,确定方差为零时对应的地层横波速度,作连续井段处理并输出反演横波速度成果图。本发明解决了现有多极子阵列声波测井仪器,在超大井径、疏松浅层或扩径严重的软地层井段不能直接测量地层横波速度的问题,拓宽了井孔声波测井技术在石油勘探与开发过程中的应用领域和提高了其应用价值。

Figure 200710090386

The invention provides a method for retrieving shear wave velocity in soft strata by leakage mode waves, which belongs to the application of geophysical acoustic wave logging technology. The specific steps are as follows: Input the unipolar or dipole leakage mode wave waveform of the well section in the soft formation with severe diameter expansion and the conventional logging data of the formation in the corresponding well section, and then extract the phase velocity dispersion curve of the leakage mode wave from the array waveform , and determine the main energy frequency band of the leakage mode wave, calculate the variance of the phase velocity of the leakage mode wave extracted from the theoretical calculation and the field array waveform under different formation shear wave velocities, determine the corresponding formation shear wave velocity when the variance is zero, and use it as a continuous well section Process and output the result map of the inversion shear wave velocity. The invention solves the problem that the existing multipole array acoustic logging instrument cannot directly measure the shear wave velocity of the formation in the ultra-large borehole diameter, loose shallow layer or soft formation well section with serious diameter expansion, and broadens the scope of the borehole acoustic logging technology. The application fields in the process of petroleum exploration and development and its application value have been improved.

Figure 200710090386

Description

A kind of method of inversing soft ground horizontal wave velocity by leaky mode waves
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 , k ) = | Σ n = 1 N X n * ( f ) e - ikd ( n - 1 ) | N Σ n = 1 N x n * ( f ) X n ( f )
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.一种泄漏模式波反演软地层横波速度的方法,它是基于内存256M以上的微型计算机的Windows环境,利用Microsoft VC++软件开发工具,通过如下方式实现:1. A method for inversion of soft formation shear wave velocity by leaky mode waves, which is based on the Windows environment of a microcomputer with a memory of more than 256M, utilizes Microsoft VC++ software development tools, and realizes in the following manner: 首先,将扩径严重的软地层井段一个深度点的由多极子阵列声波测井仪器单极或偶极工作模式测量的泄漏模式波的波形数据和对应的该深度点由常规测井测量的地层纵波速度、井径、泥浆声速、泥浆密度参数输入给计算机的几个存储器;Firstly, the waveform data of the leakage mode wave measured by the monopole or dipole working mode of the multipole subarray acoustic logging tool at a depth point in the soft formation well section with severe diameter expansion and the corresponding depth point measured by conventional logging The formation compressional wave velocity, borehole diameter, mud sound velocity, and mud density parameters are input to several memories of the computer; 其次,利用频谱加权相关法、Prony方法、修正的Matrix Pencil法或同态处理方法的阵列波形频散曲线提取技术,由软地层的单极子阵列波形或偶极子阵列波形提取出泄漏模式波的相速度频散曲线,并确定单极泄漏模式波或偶极泄漏模式波的主要能量频段;Secondly, using spectrum weighted correlation method, Prony method, modified Matrix Pencil method or the array waveform dispersion curve extraction technology of homomorphic processing method, the leakage mode wave is extracted from the monopole array waveform or dipole array waveform in soft formation. The phase velocity dispersion curve of , and determine the main energy frequency band of unipolar leakage mode wave or dipole leakage mode wave; 再次,根据Kurkjian提出的充液井孔声场理论,建立井孔多极子声波测井模型,模型中的地层纵波速度、密度、井径、泥浆声速和泥浆密度由传统的测井方法测量已知,改变地层横波速度用迭代法或抛物线法由理论的频散方程D(f,k,Vs,Vc,Vf,ρ,ρf,R)=0解出软地层单极或偶极泄漏模式波的相速度随频率变化的频散曲线,式中,k=2πf/Vp为波数,f-频率,Vp-相速度,Vs-地层横波速度,Vc-地层纵波速度,Vf-井内泥浆声速,ρ-地层密度,ρf-泥浆密度,R-井眼半径;Thirdly, according to the theory of fluid-filled wellbore sound field proposed by Kurkjian, a wellbore multipole acoustic logging model is established. The formation compressional wave velocity, density, borehole diameter, mud sound velocity and mud density in the model are known by traditional logging methods. , to change the formation shear wave velocity, use the iterative method or parabola method to solve the soft formation monopole or dipole from the theoretical dispersion equation D(f, k, V s , V c , V f , ρ, ρ f , R)=0 Dispersion curve of phase velocity of leakage mode wave changing with frequency, where k=2πf/V p is wave number, f-frequency, V p -phase velocity, V s -formation shear wave velocity, V c -formation compressional wave velocity, V f - sound velocity of mud in the well, ρ - formation density, ρ f - mud density, R - borehole radius; 然后,在单极或偶极泄漏模式波主要能量频段内,计算不同地层横波速度下理论计算的泄漏模式波相速度与现场阵列波形资料提取的泄漏模式波相速度的方差之和,得到随地层横波速度变化的方差曲线;利用拉格朗日法对方差曲线精确插值,确定出方差为零时对应的地层横波速度,即为该深度点反演得到的地层横波速度;Then, in the main energy frequency band of the monopole or dipole leakage mode wave, the sum of the variance of the leakage mode wave phase velocity calculated theoretically under different formation shear wave velocities and the leakage mode wave phase velocity extracted from the field array waveform data is calculated to obtain the The variance curve of the shear wave velocity change; use the Lagrangian method to accurately interpolate the variance curve to determine the corresponding formation shear wave velocity when the variance is zero, which is the formation shear wave velocity obtained by inversion at this depth point; 最后,利用上述反演流程对扩径严重的软地层井段进行连续处理,获得软地层井段的地层横波速度成果图,将成果图显示在显示器上,或打印出图。Finally, use the above-mentioned inversion process to continuously process the well sections in soft formations with severe diameter expansion, and obtain the formation shear wave velocity results map of the soft formation well sections, and display the results map on the monitor or print out the map.
CN200710090386XA 2007-04-09 2007-04-09 Process for inversing soft ground horizontal wave velocity by leaky mode waves Expired - Fee Related CN101285381B (en)

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 (10)

* Cited by examiner, † Cited by third party
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 成都理工大学 A method for calculating formation shear wave slowness and attenuation based on dipole acoustic logging
CN119199998B (en) * 2024-11-28 2025-05-23 中国石油大学(华东) Method, device, equipment and storage medium for processing logging while drilling transverse wave speed

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6831876B1 (en) * 2003-07-09 2004-12-14 Goodrich Corporation Acoustic window
CN1601303A (en) * 2004-10-26 2005-03-30 大庆油田有限责任公司 Selection method of sound source frequency in dipole transverse wave well logging of bushing well
CN1601304A (en) * 2004-10-26 2005-03-30 大庆油田有限责任公司 Ground stress multifrequency reverse demonstration method of dipole trans verse wave well logging

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6831876B1 (en) * 2003-07-09 2004-12-14 Goodrich Corporation Acoustic window
CN1601303A (en) * 2004-10-26 2005-03-30 大庆油田有限责任公司 Selection method of sound source frequency in dipole transverse wave well logging of bushing well
CN1601304A (en) * 2004-10-26 2005-03-30 大庆油田有限责任公司 Ground stress multifrequency reverse demonstration method of dipole trans verse wave well logging

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
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
US8379483B2 (en) Method and apparatus for determining radial shear velocity variation from dipole acoustic logging
CN110318740B (en) A method for evaluating formation anisotropy by sonic logging while drilling
CN110456418A (en) Processing and Interpretation Method of Array Acoustic Imaging Logging Data
CN103345001B (en) A kind of method of measuring the fossil lake pool depth of water
CN104407381B (en) Method for acquiring stratum rock brittleness property through radial change of borehole elastic wave velocity
CN103827696A (en) Multi-well anisotropy inversion
MXPA05013185A (en) Methods and systems for acoustic waveform processing.
US9891334B2 (en) System and methodology for determining fracture attributes in a formation
CN110348135A (en) A kind of method of acoustic logging while drilling evaluation in-place permeability
CN114017014A (en) Acoustic wave logging wave acoustics method
Xu et al. Estimation of elastic wave velocity from logging-while-drilling measurement in unconsolidated formation of accretionary wedge in nankai trough
CN112034516B (en) A method for evaluating the three-dimensional anisotropy of shear waves around wells
Che et al. Azimuthal cement evaluation with an acoustic phased-arc array transmitter: numerical simulations and field tests
Sinha et al. Hydraulic fracture characterization using cased-hole sonic data
RU2201606C1 (en) Method of typification and correlation of oil and gas productive rocks by borehole spectral-time parameters
Maalouf et al. Predicting the Elastic Properties of Anisotropic Reservoirs from Dipole Logging While Drilling Using Machine Learning
CN114961700B (en) Inversion and Imaging Method of Radial Profile of P-wave Velocity in Three-dimensional Acoustic Logging
Liu et al. Inversion of gas hydrate saturation and solid frame permeability in a gas hydrate-bearing sediment by Stoneley wave attenuation
Li et al. A Novel Method for Evaluating Formation Fracturing Effect Utilizing Acoustic Logging
Zheng et al. Identifying stress-induced anisotropy and stress direction using cross-dipole acoustic logging
Wang et al. Interpretation of sonic waveforms acquired in high-angle and horizontal wells
Zheng Dispersion Asymptotic Analysis and its Applications in Acoustic Logging

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