CN105093295B - Continuous deep processing method and device for horizontal transverse wave sensitivity - Google Patents

Abstract

The invention provides a continuous depth processing method and a device for horizontal transverse wave sensitivity, wherein the processing method comprises the following steps: acquiring logging data, calculating a first rigidity coefficient and a second rigidity coefficient according to the logging data, and calculating a third rigidity coefficient and a fourth rigidity coefficient according to the first rigidity coefficient and the second rigidity coefficient; substituting the first rigidity coefficient, the second rigidity coefficient, the third rigidity coefficient and the fourth rigidity coefficient into a Stoneley wave frequency dispersion equation, and calculating an equivalent modulus of the instrument; inverting the Stoneley wave according to the equivalent modulus of the instrument to generate a horizontal transverse wave speed; and calculating the horizontal transverse wave sensitivity according to the horizontal transverse wave speed. By the method, the reliability of the Stoneley wave for inverting the horizontal transverse wave can be inspected, so that the elastic property of the unconventional oil and gas stratum can be more accurately characterized, and a good foundation is laid for subsequent evaluation of minimum horizontal principal stress, fracture pressure and the like.

Description

Continuous deep processing method and device for horizontal transverse wave sensitivity

Technical Field

The invention relates to an application technology of acoustic logging in petroleum exploration and development, in particular to a continuous depth processing method and a continuous depth processing device for horizontal transverse wave sensitivity.

Background

Along with the continuous deep development of oil gas exploration and development, unconventional oil gas such as dense gas, shale gas, coal bed gas, dense oil and the like shows huge potential under economic and technical conditions, and global oil gas resources are expanded for a second time. Experiments have shown that unconventional hydrocarbon reservoirs are typically transversely isotropic formations. Evaluating elastic properties of transversely isotropic formations5 rigidity coefficients, i.e. C₁₁,C₁₃,C₃₃,C₄₄,C₆₆According to the Annie model, C₁₃＝f(C₃₃,C₄₄)＝C₁₂＝f(C₁₁,C₆₆) In which C is₃₃And C₄₄Can be obtained from conventional acoustic logging data, if C can be obtained₆₆The elastic properties of unconventional hydrocarbon reservoirs may then be evaluated. In patent US6920082B2, a method for inverting horizontal shear waves using Stoneley waves is provided, but only a single point inversion method for horizontal shear waves is provided in this patent, and it is explicitly stated that this method is only applicable to slow formations.

In fact, the method for inverting the horizontal transverse wave by the stoneley wave is not only suitable for the slow stratum, but also has higher sensitivity of the horizontal transverse wave to the stoneley wave in the slow stratum, namely the inverted horizontal transverse wave is more reliable. However, in a well with a depth of six to seven kilometers, both a slow stratum and a fast stratum may exist, and it is an urgent problem to determine the sensitivity of horizontal transverse waves to stoneley waves continuously and deeply and further investigate the credibility of the stoneley waves for inverting the horizontal transverse waves.

Disclosure of Invention

The main purpose of the embodiments of the present invention is to provide a method and an apparatus for continuous depth processing of horizontal transverse wave sensitivity, which determine the sensitivity of a horizontal transverse wave to a Stoneley wave by continuous depth, and further examine the reliability of inversion of the horizontal transverse wave by the Stoneley wave.

In order to achieve the above object, an embodiment of the present invention provides a continuous depth processing method for horizontal shear wave sensitivity, where the method includes: obtaining logging data, and calculating a first rigidity coefficient C according to the logging data₃₃And a second rigidity coefficient C₄₄And according to said first rigidity coefficient C₃₃And a second rigidity coefficient C₄₄Calculating a third stiffness coefficient C₁₂And a fourth rigidity coefficient C₁₃(ii) a Wherein,C₁₂＝C₁₃＝C₃₃-2C₄₄rho is density, V_PIs the velocity of longitudinal wave, V_SIs the transverse wave velocity; the first rigidity coefficient C₃₃A second rigidity coefficient C₄₄Third rigidity coefficient C₁₂And a fourth rigidity coefficient C₁₃Substituting the Stoneley wave frequency dispersion equation and calculating the equivalent modulus of the instrument; inverting the Stoneley wave according to the equivalent modulus of the instrument to generate a horizontal transverse wave speed; and calculating the horizontal transverse wave sensitivity according to the horizontal transverse wave speed.

In one embodiment, after calculating the instrument equivalent modulus, before inverting the stoneley wave according to the instrument equivalent modulus to generate the horizontal shear wave velocity, the method further comprises: and judging whether the Stoneley wave parameter meets the frequency spectrum weighted average slowness theorem, and if so, outputting the equivalent modulus of the instrument.

In one embodiment, the spectrum weighted average slowness theorem is as follows:

wherein,is the slowness of the Stoneley wave obtained from a non-dispersive array waveform processing method; s_ST(ω,M_T) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves.

In one embodiment, the stoneley wave dispersion equation is:

D(k,ω,V_SH,C₄₄,M_T,R,a,ρ,V_f,ρ_f,C₁₁,C₁₃,C₃₃) 0, where k is wavenumber, V_SHFor horizontal transverse wave, R is the borehole radius, a is the instrument radius, V_fAnd ρ_fFluid velocity and density, respectively，C₁₁Is a coefficient of rigidity, and C₁₁＝C₁₂+2C₆₆，

In one embodiment, the horizontal shear wave velocity is calculated by the following formula:

wherein,is the slowness of the Stoneley wave obtained from a non-dispersive array waveform processing method; s_ST(ω,V_SH) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves.

In one embodiment, the horizontal shear wave sensitivity is calculated by the following formula:

wherein, Sensitivity is Sensitivity; v_SHHorizontal transverse waves; omega is angular frequency; v_phase(ω) is the Stoneley wave phase velocity spectrum.

The embodiment of the invention also provides a continuous depth processing device for the horizontal transverse wave sensitivity, which comprises: a rigidity coefficient calculation unit for acquiring logging data and calculating a first rigidity coefficient C according to the logging data₃₃And a second rigidity coefficient C₄₄And according to said first rigidity coefficient C₃₃And a second rigidity coefficient C₄₄Calculating a third stiffness coefficient C₁₂And a fourth rigidity coefficient C₁₃(ii) a Wherein,C₁₂＝C₁₃＝C₃₃-2C₄₄rho is density, V_PIs the velocity of longitudinal wave, V_SIs the transverse wave velocity; an instrument equivalent modulus calculation unit for calculating the first rigidity coefficient C₃₃A second rigidity coefficient C₄₄Third rigidity coefficient C₁₂And a fourth rigidity coefficient C₁₃Substituting the Stoneley wave frequency dispersion equation and calculating the equivalent modulus of the instrument; the horizontal transverse wave velocity generating unit is used for inverting the Stoneley wave according to the equivalent modulus of the instrument to generate a horizontal transverse wave velocity; and the horizontal transverse wave sensitivity calculating unit is used for calculating the horizontal transverse wave sensitivity according to the horizontal transverse wave speed.

In an embodiment, the apparatus further includes: and the instrument equivalent modulus output unit is used for judging whether the Stoneley wave parameter meets the frequency spectrum weighted average slowness theorem or not, and if so, outputting the instrument equivalent modulus.

In one embodiment, the spectrum weighted average slowness theorem is as follows:

wherein,is the slowness of the Stoneley wave obtained from the non-dispersive array waveform processing device; s_ST(ω,M_T) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves.

In one embodiment, the stoneley wave dispersion equation is:

D(k,ω,V_SH,C₄₄,M_T,R,a,ρ,V_f,ρ_f,C₁₁,C₁₃,C₃₃) 0, where k is wavenumber, V_SHFor horizontal transverse wave, R is the borehole radius, a is the instrument radius, V_fAnd ρ_fRespectively fluid velocity and density, C₁₁Is a coefficient of rigidity, and C₁₁＝C₁₂+2C₆₆，

In one embodiment, the horizontal shear wave velocity is calculated by the following formula:

wherein,is the slowness of the Stoneley wave obtained from the non-dispersive array waveform processing device; s_ST(ω,V_SH) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves.

In one embodiment, the horizontal shear wave sensitivity is calculated by the following formula:

wherein, Sensitivity is Sensitivity; v_SHHorizontal transverse waves; omega is angular frequency; v_phase(ω) is the Stoneley wave phase velocity spectrum.

The method and the device have the advantages that the credibility of the Stoneley wave for inverting the horizontal transverse wave can be inspected, so that the elastic property of the unconventional oil and gas stratum can be more accurately characterized, and a good foundation is laid for subsequent evaluation of minimum horizontal principal stress, fracture pressure and the like.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without inventive labor.

FIGS. 1 and 2 are flow charts of a method for continuous depth processing of horizontal shear wave sensitivity according to an embodiment of the present invention;

FIG. 3 is a graph illustrating the results of continuous depth processing of horizontal shear sensitivity using well log data according to an embodiment of the present invention;

fig. 4 and 5 are schematic structural views of a continuous depth processing apparatus for horizontal shear wave sensitivity according to an embodiment of the present invention.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

The embodiment of the invention provides a method and a device for continuously and deeply processing horizontal transverse wave sensitivity. The present invention will be described in detail below with reference to the accompanying drawings.

The embodiment of the invention provides a continuous depth processing method of horizontal transverse wave sensitivity, as shown in fig. 1, the processing method mainly comprises the following steps:

step S101: acquiring logging data, calculating a first rigidity coefficient and a second rigidity coefficient according to the logging data, and calculating a third rigidity coefficient and a fourth rigidity coefficient according to the first rigidity coefficient and the second rigidity coefficient;

step S102: substituting the first rigidity coefficient, the second rigidity coefficient, the third rigidity coefficient and the fourth rigidity coefficient into a Stoneley wave frequency dispersion equation, and calculating an equivalent modulus of the instrument;

step S103: inverting the Stoneley wave according to the equivalent modulus of the instrument to generate a horizontal transverse wave speed;

step S104: and calculating the horizontal transverse wave sensitivity according to the horizontal transverse wave speed.

Through the steps S101 to S104, the sensitivity of the horizontal transverse wave to the Stoneley wave can be continuously and deeply judged by combining the logging data and the frequency dispersion characteristic of the Stoneley wave, the reliability of the Stoneley wave for inverting the horizontal transverse wave can be inspected, the elastic property of the unconventional oil and gas stratum can be accurately represented, and a good foundation is laid for subsequent evaluation of minimum horizontal main stress, fracture pressure and the like.

The following describes in detail the method for continuous depth processing of horizontal shear wave sensitivity according to the embodiment of the present invention with reference to the above steps.

In the step S101, logging data is obtained, a first stiffness coefficient and a second stiffness coefficient are calculated according to the logging data, and a third stiffness coefficient and a fourth stiffness coefficient are calculated according to the first stiffness coefficient and the second stiffness coefficient.

Wherein, the logging data comprises: the data logging method comprises density logging data and acoustic logging data, wherein the acoustic logging data comprise longitudinal wave data and transverse wave data.

Firstly, the rigidity coefficient C is calculated according to the logging data₃₃、C₄₄Specifically, the formula is combined by using acoustic logging data and density logging dataSeparately calculating the rigidity coefficients C₃₃、C₄₄Where ρ is density, V_PIs the velocity of longitudinal wave, V_SIs the shear wave velocity.

Obtaining the rigidity coefficient C in the calculation₃₃、C₄₄Then, the rigidity coefficient C can be calculated₁₂And C₁₃. Specifically, the above is calculatedCoefficient of rigidity C₁₃Is according to equation C in the Annie model₁₃＝C₃₃-2C₄₄Is calculated, wherein C₃₃And C₄₄The stiffness coefficient is the one determined above. It should be noted that the anie model adopted in the embodiment of the present invention is applicable to shale reservoirs, and for other unconventional oil and gas reservoirs, corresponding relations and calculation methods may be adopted, which is not limited to the present invention.

Calculating the stiffness coefficient C as described above₁₂Then it is according to equation C in the anie model₁₂＝C₁₃To calculate. Similarly, the equation is applicable to shale reservoirs, and for other unconventional hydrocarbon reservoirs, corresponding relational equations and calculation methods can be adopted, and the invention is not limited thereto.

After each stiffness coefficient is calculated in step S101, step S102 may be executed to calculate the equivalent modulus of the instrument by substituting the first stiffness coefficient, the second stiffness coefficient, the third stiffness coefficient, and the fourth stiffness coefficient into the stoneley wave dispersion equation. The equivalent modulus of the instrument is that under the condition that the wavelength of sound waves is far larger than the radius of a logging instrument, the instrument can be equivalently simulated into an elastic rod, and the elastic parameter of the elastic rod can be replaced by the equivalent modulus of elasticity.

The stoneley wave frequency dispersion equation is a curve equation of the variation of the slowness of the stoneley wave along with the frequency, and in the embodiment of the invention, the stoneley wave frequency dispersion equation is as follows:

D(k,ω,V_SH,C₄₄,M_T,R,a,ρ,V_f,ρ_f,C₁₁,C₁₃,C₃₃)＝0，

wherein k is the wave number, V_SHHorizontal shear wave velocity, R is borehole radius, a is instrument radius, V_fAnd ρ_fFluid velocity and density, respectively, which are obtained from well log data and drilling mud; c₁₁Is a coefficient of rigidity, and C₁₁＝C₁₂+2C₆₆，

Under the condition that other parameters are known or measurable, the equivalent modulus M of the instrument can be obtained by inputting each rigidity coefficient into the Stoneley wave frequency dispersion equation_T。

After the instrument equivalent modulus is obtained by calculation, in step S103, the stoneley wave is inverted according to the instrument equivalent modulus to generate the horizontal shear wave velocity. Horizontal shear waves are waves that are parallel to the isotropic plane for the transversely isotropic formation vibration direction and propagation direction. Inversion of horizontal shear waves from stoneley waves is the inversion of horizontal shear waves to determine their magnitude using their sensitivity to stoneley waves. Specifically, in step S103, the horizontal shear wave velocity is calculated by the following formula:

wherein,is the slowness of the Stoneley wave obtained from a non-dispersive array waveform processing method; s_ST(ω,V_SH) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves. And, by solving the Stoneley wave dispersion equation described above, the wave number k, S, corresponding to each ω can be found_ST(ω)＝k/ω。

In practical applications, as shown in fig. 2, between the step S102 and the step S103, the method for continuous depth processing of horizontal shear wave sensitivity according to the embodiment of the present invention may further include a determining step S105: and judging whether the Stoneley wave parameter meets the frequency spectrum weighted average slowness theorem, and if so, outputting the equivalent modulus of the instrument.

Specifically, the spectrum weighted average slowness theorem is as follows:

wherein,is the slowness of the Stoneley wave obtained from a non-dispersive array waveform processing method; s_ST(ω,M_T) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves.

The slowness theorem of the frequency spectrum weighted average shows that the slowness obtained by the time domain non-dispersive array processing method is the weighted average of the slowness dispersion curve of the frequency domain, wherein the weight function is omega²A²(ω)。

Further, after the horizontal shear wave velocity is generated in step S103, step S104 is executed to calculate the horizontal shear wave sensitivity from the horizontal shear wave velocity. In step S104, the horizontal shear wave sensitivity is calculated by the following formula:

wherein, Sensitivity is Sensitivity; v_SHHorizontal transverse waves; omega is angular frequency; v_phase(ω) is the Stoneley wave phase velocity spectrum.

After the above steps, the result of the horizontal shear wave sensitivity is shown in fig. 3, wherein the first path is a depth path and represents the distance from the measurement well section (i.e. the target layer) to the wellhead; the second path is a natural gamma curve (GR) and a well diameter curve, wherein the solid line is the natural gamma curve and represents the change of lithology; the dotted line is a well diameter curve and represents the quality of the well; the third is the density and longitudinal wave time difference curve, where the solid line is the longitudinal wave time difference curve and the dotted line is the density curve, usually used to calculate the porosity, here used to calculate the stiffness coefficient; fourth passThe time difference curve of the transverse wave and the time difference curve of the Stoneley wave are shown, wherein the solid line is the time difference curve of the Stoneley wave and is used for horizontal transverse wave inversion; the dotted line is a transverse wave time difference curve which is used for calculating a rigidity coefficient; the fifth and sixth curves are stiffness coefficient curves, wherein the fifth curve is C₃₃,C₄₄The solid line is C₃₃Curve, dashed line C₄₄A curve; the sixth channel is respectively C₆₆,C₁₃The solid line is C₆₆Curve, dotted line is C₁₃The curve is used for representing the elastic property of the stratum; the seventh trace is a sensitivity curve used to characterize the horizontal shear wave sensitivity.

By the continuous depth processing method for the horizontal transverse wave sensitivity, the reliability of the Stoneley wave for inverting the horizontal transverse wave can be inspected, so that the elastic property of an unconventional oil and gas stratum can be more accurately characterized, and a good foundation is laid for subsequent evaluation of minimum horizontal principal stress, fracture pressure and the like.

An embodiment of the present invention further provides a continuous depth processing apparatus for horizontal transverse wave sensitivity, as shown in fig. 4, the processing apparatus mainly includes: a stiffness coefficient calculation unit 1, an instrument equivalent modulus calculation unit 2, a horizontal shear wave velocity generation unit 3, a horizontal shear wave sensitivity calculation unit 4, and the like.

The stiffness coefficient calculation unit 1 is configured to obtain logging data, calculate a first stiffness coefficient and a second stiffness coefficient according to the logging data, and calculate a third stiffness coefficient and a fourth stiffness coefficient according to the first stiffness coefficient and the second stiffness coefficient.

Wherein, the logging data comprises: the data logging method comprises density logging data and acoustic logging data, wherein the acoustic logging data comprise longitudinal wave data and transverse wave data.

First, the stiffness coefficient calculation unit 1 calculates a stiffness coefficient C from the above-mentioned logging data₃₃、C₄₄Specifically, the formula is combined by using acoustic logging data and density logging dataSeparately calculating the rigidity coefficients C₃₃、C₄₄Where ρ is density, V_PIs the velocity of longitudinal wave, V_SIs the shear wave velocity.

Obtaining the rigidity coefficient C in the calculation₃₃、C₄₄Then, the rigidity coefficient C can be calculated₁₂And C₁₃. Specifically, the above-mentioned rigidity coefficient C is calculated₁₃Is according to equation C in the Annie model₁₃＝C₃₃-2C₄₄Is calculated, wherein C₃₃And C₄₄The stiffness coefficient is the one determined above. It should be noted that the anie model adopted in the embodiment of the present invention is applicable to shale reservoirs, and for other unconventional oil and gas reservoirs, corresponding relations and calculation methods may be adopted, which is not limited to the present invention.

Calculating the stiffness coefficient C as described above₁₂Then it is according to equation C in the anie model₁₂＝C₁₃To calculate. Similarly, the equation is applicable to shale reservoirs, and for other unconventional hydrocarbon reservoirs, corresponding relational equations and calculation methods can be adopted, and the invention is not limited thereto.

After the stiffness coefficients are calculated by the stiffness coefficient calculation unit 1, the instrument equivalent modulus calculation unit 2 can be triggered, and the first stiffness coefficient, the second stiffness coefficient, the third stiffness coefficient and the fourth stiffness coefficient are substituted into the stoneley wave frequency dispersion equation to calculate the instrument equivalent modulus. The equivalent modulus of the instrument is that under the condition that the wavelength of sound waves is far larger than the radius of a logging instrument, the instrument can be equivalently simulated into an elastic rod, and the elastic parameter of the elastic rod can be replaced by the equivalent modulus of elasticity.

The stoneley wave frequency dispersion equation is a curve equation of the variation of the slowness of the stoneley wave along with the frequency, and in the embodiment of the invention, the stoneley wave frequency dispersion equation is as follows:

D(k,ω,V_SH,C₄₄,M_T,R,a,ρ,V_f,ρ_f,C₁₁,C₁₃,C₃₃)＝0，

wherein k is the wave number, V_SHHorizontal shear wave velocity, R is borehole radius, a is instrument radius, V_fAnd ρ_fFluid velocity and density, respectively, which are obtained from well log data and drilling mud; c₁₁Is a coefficient of rigidity, and C₁₁＝C₁₂+2C₆₆，

Under the condition that other parameters are known or measurable, the equivalent modulus calculation unit 2 of the instrument inputs each rigidity coefficient into the Stoneley wave frequency dispersion equation to obtain the equivalent modulus M of the instrument_T。

After the instrument equivalent modulus is obtained through calculation, the stoneley wave is inverted according to the instrument equivalent modulus through the horizontal transverse wave velocity generation unit 3, and the horizontal transverse wave velocity is generated. Horizontal shear waves are waves that are parallel to the isotropic plane for the transversely isotropic formation vibration direction and propagation direction. Inversion of horizontal shear waves from stoneley waves is the inversion of horizontal shear waves to determine their magnitude using their sensitivity to stoneley waves. Specifically, the horizontal shear wave velocity generation unit 3 calculates the horizontal shear wave velocity by the following formula:

wherein,is the slowness of the Stoneley wave obtained from a non-dispersive array waveform processing method; s_ST(ω,V_SH) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves. Moreover, by solving the Stoneley wave frequency dispersion equation, the wave number k corresponding to each omega can be found, and S is_ST(ω)＝k/ω。

In practical applications, as shown in fig. 5, the apparatus for continuous depth processing of horizontal shear wave sensitivity according to the embodiment of the present invention may further include an instrument equivalent modulus output unit 5, configured to determine whether the stoneley wave parameter satisfies the spectrum weighted average slowness theorem, and if so, output the instrument equivalent modulus.

Specifically, the spectrum weighted average slowness theorem is as follows:

wherein,is the slowness of the Stoneley wave obtained from a non-dispersive array waveform processing method; s_ST(ω,M_T) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves.

The slowness theorem of the frequency spectrum weighted average shows that the slowness obtained by the time domain non-dispersive array processing method is the weighted average of the slowness dispersion curve of the frequency domain, wherein the weight function is omega²A²(ω)。

Further, after the horizontal shear wave velocity generating means 3 generates the horizontal shear wave velocity, the horizontal shear wave sensitivity can be calculated from the horizontal shear wave velocity by the above-described horizontal shear wave sensitivity calculating means 4. Specifically, the horizontal shear wave sensitivity calculating unit 4 calculates the horizontal shear wave sensitivity by the following formula:

wherein, Sensitivity is Sensitivity; v_SHHorizontal transverse waves; omega is angular frequency; v_phase(ω) is the Stoneley wave phase velocity spectrum.

After the processes performed by the units, the result of the horizontal shear wave sensitivity is shown in fig. 3, where the first path is a depth path and represents the distance from the measurement well section (i.e., the target layer) to the wellhead; the second path is a natural gamma curve (GR) and a well diameter curve, wherein the solid line is the natural gamma curve and represents the change of lithology; the dotted line is a well diameter curve and represents the quality of the well; the third is the density and longitudinal wave time difference curve, where the solid line is the longitudinal wave time difference curve and the dotted line is the density curve, usually used to calculate the porosity, here used to calculate the stiffness coefficient; the fourth path is a transverse wave time difference curve and a Stoneley wave time difference curve, wherein the solid line is the Stoneley wave time difference curve and is used for horizontal transverse wave inversion; the dotted line is a transverse wave time difference curve which is used for calculating a rigidity coefficient; the fifth and sixth curves are stiffness coefficient curves, wherein the fifth curve is C₃₃,C₄₄The solid line is C₃₃Curve, dashed line C₄₄A curve; the sixth channel is respectively C₆₆,C₁₃The solid line is C₆₆Curve, dotted line is C₁₃The curve is used for representing the elastic property of the stratum; the seventh trace is a sensitivity curve used to characterize the horizontal shear wave sensitivity.

Through the continuous depth processing device for the horizontal transverse wave sensitivity, provided by the embodiment of the invention, the credibility of the Stoneley wave for inverting the horizontal transverse wave can be inspected, so that the elastic property of an unconventional oil and gas stratum can be more accurately characterized, and a good foundation is laid for subsequent evaluation of minimum horizontal principal stress, fracture pressure and the like.

It will be understood by those skilled in the art that all or part of the steps in the method for implementing the above embodiments may be implemented by relevant hardware instructed by a program, and the program may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, optical disk, etc.

The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (12)

1. A continuous depth processing method for horizontal shear wave sensitivity is characterized by comprising the following steps:

obtaining logging data, and calculating a first rigidity coefficient C according to the logging data₃₃And a second rigidity coefficient C₄₄And according to said first rigidity coefficient C₃₃And a second rigidity coefficient C₄₄Calculating a third stiffness coefficient C₁₂And a fourth rigidity coefficient C₁₃(ii) a Wherein,C₁₂＝C₁₃＝C₃₃-2C₄₄rho is density, V_PIs the velocity of longitudinal wave, V_SIs the transverse wave velocity;

the first rigidity coefficient C₃₃A second rigidity coefficient C₄₄Third rigidity coefficient C₁₂And a fourth rigidity coefficient C₁₃Substituting the Stoneley wave frequency dispersion equation and calculating the equivalent modulus of the instrument;

inverting the Stoneley wave according to the equivalent modulus of the instrument to generate a horizontal transverse wave speed;

and calculating the horizontal transverse wave sensitivity according to the horizontal transverse wave speed.

2. The method for continuous depth processing of horizontal shear wave sensitivity according to claim 1, wherein after calculating the instrument equivalent modulus and before inverting the stoneley wave according to the instrument equivalent modulus to generate the horizontal shear wave velocity, the method further comprises:

and judging whether the Stoneley wave parameter meets the frequency spectrum weighted average slowness theorem, and if so, outputting the equivalent modulus of the instrument.

3. The method of continuous depth processing of horizontal shear sensitivity of claim 2, wherein the spectral weighted average slowness theorem is:

$S_{ST}^{*}=\frac{\underset{\mathrm{\Ω}}{\∫}{S}_{ST}(\mathrm{\ω},M_T){\mathrm{\ω}}^2{A}^2\left(\mathrm{\ω}\right)d\mathrm{\ω}}{\underset{\mathrm{\Ω}}{\∫}{\mathrm{\ω}}^2{A}^2\left(\mathrm{\ω}\right)d\mathrm{\ω}},$

wherein,is the slowness of the Stoneley wave obtained from a non-dispersive array waveform processing method; s_ST(ω,M_T) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves; m_TIs the instrument equivalent modulus.

4. The continuous depth processing method for horizontal shear wave sensitivity according to claim 1, wherein the Stoneley wave frequency dispersion equation is as follows:

D(k,ω,V_SH,C₄₄,M_T,R,a,ρ,V_f,ρ_f,C₁₁,C₁₃,C₃₃)＝0，

wherein k is the wave number, V_SHHorizontal shear wave velocity, R is borehole radius, a is instrument radius, V_fAnd ρ_fRespectively fluid velocity and density, C₁₁Is a coefficient of rigidity, and C₁₁＝C₁₂+2C₆₆，M_TIs the instrument equivalent modulus; ω is the angular frequency.

5. The continuous depth processing method of horizontal shear wave sensitivity according to claim 1, wherein the horizontal shear wave velocity is calculated by the following formula:

$S_{ST}^{*}=\frac{\underset{\mathrm{\Ω}}{\∫}{S}_{ST}(\mathrm{\ω},V_{SH}){\mathrm{\ω}}^2{A}^2\left(\mathrm{\ω}\right)d\mathrm{\ω}}{\underset{\mathrm{\Ω}}{\∫}{\mathrm{\ω}}^2{A}^2\left(\mathrm{\ω}\right)d\mathrm{\ω}},$

wherein,is the slowness of the Stoneley wave obtained from a non-dispersive array waveform processing method; s_ST(ω,V_SH) Stoneley wave slowness, V, obtained for a dispersion curve_SHHorizontal shear wave velocity; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves.

6. The continuous depth processing method of horizontal shear wave sensitivity according to claim 1, wherein the horizontal shear wave sensitivity is calculated by the following formula:

$Sensitivity=\left(\frac{V_{SH}}{V_{phase}\left(\mathrm{\ω}\right)}\right)\frac{\∂V_{phase}\left(\mathrm{\ω}\right)}{\∂V_{SH}},$

wherein, Sensitivity is Sensitivity; v_SHHorizontal shear wave velocity; omega is angular frequency; v_phase(ω) is the Stoneley wave phase velocity spectrum.

7. A continuous depth processing device for horizontal shear wave sensitivity, comprising:

a rigidity coefficient calculation unit for acquiring logging data and calculating a first rigidity coefficient C according to the logging data₃₃And a second rigidity coefficient C₄₄And according to said first rigidity coefficient C₃₃And a second rigidity coefficient C₄₄Calculating a third stiffness coefficient C₁₂And a fourth rigidity coefficient C₁₃(ii) a Wherein,C₁₂＝C₁₃＝C₃₃-2C₄₄rho is density, V_PIs the velocity of longitudinal wave, V_SIs the transverse wave velocity;

an instrument equivalent modulus calculation unit for calculating the first rigidity coefficient C₃₃A second rigidity coefficient C₄₄Third rigidity coefficient C₁₂And a fourth rigidity coefficient C₁₃Substituting the Stoneley wave frequency dispersion equation and calculating the equivalent modulus of the instrument;

the horizontal transverse wave velocity generating unit is used for inverting the Stoneley wave according to the equivalent modulus of the instrument to generate a horizontal transverse wave velocity;

and the horizontal transverse wave sensitivity calculating unit is used for calculating the horizontal transverse wave sensitivity according to the horizontal transverse wave speed.

8. The apparatus for continuous depth processing of horizontal shear wave sensitivity according to claim 7, further comprising:

and the instrument equivalent modulus output unit is used for judging whether the Stoneley wave parameter meets the frequency spectrum weighted average slowness theorem or not, and if so, outputting the instrument equivalent modulus.

9. The apparatus according to claim 8, wherein the spectral weighted mean slowness theorem is:

$S_{ST}^{*}=\frac{\underset{\mathrm{\Ω}}{\∫}{S}_{ST}(\mathrm{\ω},M_T){\mathrm{\ω}}^2{A}^2\left(\mathrm{\ω}\right)d\mathrm{\ω}}{\underset{\mathrm{\Ω}}{\∫}{\mathrm{\ω}}^2{A}^2\left(\mathrm{\ω}\right)d\mathrm{\ω}},$

wherein,is the slowness of the Stoneley wave obtained from the non-dispersive array waveform processing device; s_ST(ω,M_T) (ii) the stoneley wave slowness obtained for the dispersion curve; omega is angular frequency; a (omega) is StoneleyAn amplitude spectrum of the wave; m_TIs the instrument equivalent modulus.

10. The apparatus for continuous depth processing of horizontal shear wave sensitivity according to claim 7, wherein the Stoneley wave dispersion equation is:

D(k,ω,V_SH,C₄₄,M_T,R,a,ρ,V_f,ρ_f,C₁₁,C₁₃,C₃₃)＝0，

wherein k is the wave number, V_SHHorizontal shear wave velocity, R is borehole radius, a is instrument radius, V_fAnd ρ_fRespectively fluid velocity and density, C₁₁Is a coefficient of rigidity, and C₁₁＝C₁₂+2C₆₆，M_TIs the instrument equivalent modulus; ω is the angular frequency.

11. The apparatus according to claim 7, wherein the horizontal shear wave velocity is calculated by the following formula:

$S_{ST}^{*}=\frac{\underset{\mathrm{\Ω}}{\∫}{S}_{ST}(\mathrm{\ω},V_{SH}){\mathrm{\ω}}^2{A}^2\left(\mathrm{\ω}\right)d\mathrm{\ω}}{\underset{\mathrm{\Ω}}{\∫}{\mathrm{\ω}}^2{A}^2\left(\mathrm{\ω}\right)d\mathrm{\ω}},$

wherein,is the slowness of the Stoneley wave obtained from the non-dispersive array waveform processing device; s_ST(ω,V_SH) Stoneley wave slowness, V, obtained for a dispersion curve_SHHorizontal shear wave velocity; omega is angular frequency; a (ω) is the amplitude spectrum of Stoneley waves.

12. The apparatus according to claim 7, wherein the horizontal shear wave sensitivity is calculated by the following formula:

$Sensitivity=\left(\frac{V_{SH}}{V_{phase}\left(\mathrm{\ω}\right)}\right)\frac{\∂V_{phase}\left(\mathrm{\ω}\right)}{\∂V_{SH}},$

wherein, Sensitivity is Sensitivity; v_SHHorizontal shear wave velocity; omega is angular frequency; v_phase(ω) is the Stoneley wave phase velocity spectrum.