Fracture inclination angle identification method, device and equipment for horizontal well passing crack and medium

Abstract

The invention discloses a crack inclination angle recognition method, device and equipment for a horizontal well passing crack and a medium, and relates to the technical field of oil-gas exploration, and the method comprises the steps that well hole mode waves of a horizontal well are collected to obtain an array logging result; dipole transverse wave signal extraction is conducted on array waveforms in the array logging result, reflected wave relative amplitude analysis of various transverse waves is conducted on the basis of the corresponding dipole transverse wave signal extraction result, and the target crack dip angle range of the horizontal well passing crack is determined by means of the corresponding relative amplitude analysis result; the well-passing crack is a crack penetrated by a borehole of the horizontal well; the plurality of transverse waves comprise horizontal polarization transverse waves and vertical polarization transverse waves corresponding to the dipole transverse waves; and performing quantitative inversion based on a preset segmentation inversion model, the target crack inclination angle range and the relative amplitude analysis result to obtain a crack inclination angle identification result of the well-passing crack. According to the method, the crack inclination angle of the horizontal well passing crack can be effectively and accurately identified, and the identification cost is reduced.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil and gas exploration, and in particular relates to a fracture dip angle identification method, device, equipment and medium for a horizontal well passing through a well fracture. BACKGROUND

[0002] At present, post-stack and pre-stack technologies constitute the core path of fracture evaluation in seismic. However, due to the restriction of the wavelength scale of seismic waves, in the traditional correlation scheme, only the regional fracture parameters can be statistically evaluated, and the dip angle and opening degree of a single fracture near the horizontal well cannot be quantitatively described.

[0003] Therefore, in some existing correlation schemes, electrical imaging logging can be used to directly obtain in-situ parameters such as fracture dip angle, strike and opening degree around the wellbore. However, the measurement cost of this technology is high, and the push-type wall-attached measurement method is prone to resistance and jamming in the horizontal well, which has the risk of instrument falling, thereby limiting the application range. In view of this problem, some existing correlation schemes use Stoneley waves to evaluate the fractures passing through the well, but these schemes still cannot provide key spatial occurrence information such as the dip angle of the fracture, which will affect the guidance effect of the fracturing development work of the horizontal well in unconventional reservoirs.

[0004] In summary, how to effectively and accurately identify the fracture dip angle of the horizontal well passing through the well fracture and reduce the identification cost is a problem to be solved by those skilled in the art. SUMMARY

[0005] Therefore, the purpose of the present application is to provide a fracture dip angle identification method, device, equipment and medium for a horizontal well passing through a well fracture, which can effectively and accurately identify the fracture dip angle of the horizontal well passing through the well fracture and reduce the identification cost, thereby improving the guidance effect of the fracturing development work of the horizontal well in unconventional reservoirs. The specific scheme is as follows:

[0006] In a first aspect, the present application provides a fracture dip angle identification method for a horizontal well passing through a well fracture, comprising:

[0007] Based on a preset array acoustic logging instrument, borehole mode wave acquisition is performed on the horizontal well to complete array logging operation and obtain array logging results;

[0008] Based on a preset array waveform processing strategy, a dipole shear wave signal extraction result is determined by extracting a dipole shear wave signal from the array waveform in the array logging result;

[0009] Based on the extracted dipole shear wave signal, relative amplitude analysis of the reflected waves of various shear waves is performed, and the target fracture dip angle range corresponding to the wellbore fracture of the horizontal well is determined using the corresponding relative amplitude analysis results; wherein, the wellbore fracture is a fracture that is penetrated by the wellbore of the horizontal well; the various shear waves include horizontally polarized shear waves and vertically polarized shear waves corresponding to dipole shear waves.

[0010] Quantitative inversion is performed based on the preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results to complete the well-pass fracture dip angle identification operation and obtain the fracture dip angle identification result corresponding to the well-pass fracture.

[0011] Optionally, the step of acquiring borehole mode waves in a horizontal well using a preset array acoustic logging tool to complete the array logging operation and obtain the array logging results includes:

[0012] A preset array acoustic logging tool is used to acquire borehole mode waves in fractured formations near horizontal wells to determine the array logging results; the array logging results include array waveforms under the fractured formations; the array waveforms are orthogonal dipole shear wave waveforms.

[0013] Optionally, the step of extracting dipole shear wave signals from the array waveforms in the array logging results based on a preset array waveform processing strategy to determine the dipole shear wave signal extraction result includes:

[0014] Wave signal composition analysis is performed on the orthogonal dipole shear wave waveform in the array logging results to determine the analysis results;

[0015] Based on the analysis results, wavefield separation is performed on the orthogonal dipole shear wave waveform to determine the wavefield separation result;

[0016] Based on the wave field separation results, the direct wave signal is extracted to determine the direct wave signal extraction result;

[0017] Based on the wave field separation results, reflected wave signals are extracted to determine the reflected wave signal extraction results; the reflected wave signal extraction results include reflected wave signals excited by the well-pass fractures corresponding to the horizontal well;

[0018] Based on the results of the direct wave signal extraction and the results of the reflected wave signal extraction, the results of the dipole transverse wave signal extraction are determined.

[0019] Optionally, based on the analysis results, performing wavefield separation on the orthogonal dipole shear wave waveform to determine the wavefield separation result includes:

[0020] Apparent velocity analysis is performed on the direct wave signal and reflected wave signal in the array logging results to determine the apparent velocity analysis results;

[0021] Based on the apparent velocity analysis results and the preset filtering algorithm, wavefield separation is performed on the direct wave signal and reflected wave signal in the array logging results to determine the wavefield separation result.

[0022] Optionally, the step of performing relative amplitude analysis of reflected waves of various shear waves based on the extracted dipole shear wave signal, and determining the target fracture dip angle range corresponding to the wellbore fracture of the horizontal well using the corresponding relative amplitude analysis results, includes:

[0023] Based on the extracted dipole shear wave signal and the preset amplitude statistics rules, amplitude values ​​of shear waves in various modes are statistically analyzed to determine the target amplitude value statistical results; the target amplitude value statistical results include the direct wave amplitude value statistical results and the reflected wave amplitude value statistical results.

[0024] Based on the statistical results of the direct wave amplitude value and the statistical results of the reflected wave amplitude value, a relative amplitude analysis of the reflected wave is performed to determine the first relative amplitude analysis result corresponding to the horizontally polarized transverse wave and the second relative amplitude analysis result corresponding to the vertically polarized transverse wave.

[0025] The first relative amplitude analysis result and the relative amplitude analysis result are respectively determined to be less than a first preset amplitude threshold, so as to determine the first threshold determination result;

[0026] Determine whether the difference between the first relative amplitude analysis result and the relative amplitude analysis result is less than a second preset amplitude threshold, so as to determine the second threshold judgment result;

[0027] Based on the first threshold judgment result, the second threshold judgment result, and the preset well-pass fracture inclination angle range configuration information, the target fracture inclination angle range corresponding to the well-pass fracture of the horizontal well is determined.

[0028] Optionally, based on the extracted dipole shear wave signal and preset amplitude statistical rules, the amplitude values ​​of shear waves in various modes are statistically analyzed to determine the target amplitude value statistical result, including:

[0029] Based on the direct wave signal extraction results and the acoustic system structure parameters and the first preset time window corresponding to the preset array acoustic logging tool, amplitude values ​​are statistically analyzed to determine the statistical results of the first direct wave amplitude value corresponding to the horizontally polarized shear wave.

[0030] Based on the direct wave signal extraction results and the acoustic system structure parameters and the second preset time window corresponding to the preset array acoustic logging tool, amplitude values ​​are statistically analyzed to determine the statistical results of the second direct wave amplitude value corresponding to the vertically polarized shear wave.

[0031] Based on the extracted reflection wave signal and the acoustic system structure parameters and third preset time window corresponding to the preset array acoustic logging tool, amplitude values ​​are statistically analyzed to determine the statistical results of the first reflection wave amplitude value corresponding to the horizontally polarized shear wave.

[0032] Based on the extracted reflected wave signal and the acoustic system structure parameters and fourth preset time window corresponding to the preset array acoustic logging tool, amplitude values ​​are statistically analyzed to determine the statistical results of the second reflected wave amplitude value corresponding to the vertically polarized transverse wave.

[0033] Optionally, the quantitative inversion based on the preset piecewise inversion model, the target fracture dip angle range, and the relative amplitude analysis results to complete the well-through fracture dip angle identification operation includes:

[0034] Obtain the preset segmented inversion model;

[0035] Based on the target fracture dip angle range and the preset segmented inversion model, determine the target inversion model corresponding to the through-well fracture;

[0036] Based on the relative amplitude analysis results and the target inversion model, quantitative inversion is performed to complete the well-pass fracture dip angle identification operation and obtain the fracture dip angle identification result corresponding to the well-pass fracture.

[0037] Secondly, this application provides a fracture inclination angle identification device for horizontal well crossing fractures, comprising:

[0038] The array logging module is used to acquire borehole mode waves in horizontal wells based on a preset array acoustic logging tool, so as to complete the array logging operation and obtain the array logging results.

[0039] The shear wave signal extraction module is used to extract dipole shear wave signals from the array waveform in the array logging results based on a preset array waveform processing strategy, so as to determine the dipole shear wave signal extraction result.

[0040] The relative amplitude analysis module is used to perform relative amplitude analysis of the reflected waves of various shear waves based on the extraction results of the dipole shear wave signal, and to determine the target fracture dip angle range corresponding to the wellbore fracture of the horizontal well using the corresponding relative amplitude analysis results; wherein, the wellbore fracture is a fracture that is penetrated by the wellbore of the horizontal well; the various shear waves include horizontally polarized shear waves and vertically polarized shear waves corresponding to the dipole shear wave;

[0041] The quantitative inversion module is used to perform quantitative inversion based on the preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results, so as to complete the well-pass fracture dip angle identification operation and obtain the fracture dip angle identification result corresponding to the well-pass fracture.

[0042] Thirdly, this application provides an electronic device, comprising:

[0043] Memory, used to store computer programs;

[0044] A processor is used to execute the computer program to implement the steps of the aforementioned method for identifying the inclination angle of fractures in horizontal wells.

[0045] Fourthly, this application provides a computer-readable storage medium for storing a computer program, which, when executed by a processor, implements the steps of the aforementioned method for identifying the inclination angle of horizontal well crossing fractures.

[0046] As can be seen, in this application, based on a preset array acoustic logging tool, borehole mode waves are acquired in a horizontal well to complete the array logging operation and obtain the array logging results; based on a preset array waveform processing strategy, dipole shear wave signals are extracted from the array waveforms in the array logging results to determine the dipole shear wave signal extraction results; based on the dipole shear wave signal extraction results, relative amplitude analysis of the reflected waves of various shear waves is performed, and the target fracture dip angle range corresponding to the wellbore fracture of the horizontal well is determined using the corresponding relative amplitude analysis results; wherein, the wellbore fracture is a fracture that is penetrated by the wellbore of the horizontal well; the various shear waves include horizontally polarized shear waves and vertically polarized shear waves corresponding to the dipole shear waves; based on a preset piecewise inversion model, the target fracture dip angle range, and the relative amplitude analysis results, quantitative inversion is performed to complete the wellbore fracture dip angle identification operation and obtain the fracture dip angle identification result corresponding to the wellbore fracture. In other words, this application first obtains the array logging results of a horizontal well. Then, it extracts the dipole shear wave signal from the array logging results. Using the corresponding dipole shear wave signal extraction results, it performs relative amplitude analysis of the reflected waves of various shear waves to determine the relative amplitude analysis results. Subsequently, using the relative amplitude analysis results, it determines the target fracture dip angle range of the through-well fracture. Combined with a pre-set segmented inversion model and the relative amplitude analysis results, it performs quantitative inversion to obtain the fracture dip angle identification result corresponding to the through-well fracture. This approach effectively and accurately identifies the fracture dip angle of through-well fractures in horizontal wells, reduces identification costs, and thus improves the guidance effect for fracturing development work of horizontal wells in unconventional reservoirs. Attached Figure Description

[0047] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0048] Figure 1 A flowchart of a method for identifying the inclination angle of a horizontal well through a well, provided in this application;

[0049] Figure 2 A schematic diagram of a horizontal well through-hole fracture provided in this application;

[0050] Figure 3(a) is an array waveform diagram of vertically polarized shear waves measured in a fractured formation provided in this application;

[0051] Figure 3(b) is an array waveform diagram of a horizontally polarized shear wave measured in a fractured formation provided in this application;

[0052] Figure 4(a) is a schematic diagram of the array waveform of a vertically polarized transverse wave after wave field separation processing provided in this application;

[0053] Figure 4(b) is a schematic diagram of the array waveform of a horizontally polarized transverse wave after wave field separation processing provided in this application;

[0054] Figure 5 A schematic diagram illustrating the relative amplitude analysis results of reflected waves of horizontally polarized transverse waves and vertically polarized transverse waves provided in this application;

[0055] Figure 6 A schematic diagram illustrating the process for determining the dip angle range of a target crack, as provided in this application;

[0056] Figure 7 A schematic diagram illustrating the quantitative calculation of crack dip angle using a pre-defined piecewise inversion model, provided in this application;

[0057] Figure 8 A schematic diagram of a fracture inclination angle identification device for a horizontal well through-well fracture provided in this application;

[0058] Figure 9 This application provides a structural diagram of an electronic device. Detailed Implementation

[0059] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0060] In some existing solutions, electrical imaging logging can directly obtain in-situ parameters such as the dip angle, orientation, and aperture of fractures around the wellbore. However, this technology is costly, and its push-and-hold method is prone to obstruction and jamming in horizontal wells, posing a risk of instrument drop and limiting its application. To address this issue, some existing solutions use Stoneley waves to evaluate fractures around the well, but these solutions still cannot provide key spatial occurrence information such as fracture dip angle, thus affecting the guidance for fracturing development of horizontal wells in unconventional reservoirs.

[0061] Therefore, this application provides a fracture inclination angle identification scheme for horizontal well through-well fractures, which can effectively and accurately identify the fracture inclination angle of horizontal well through-well fractures and reduce identification costs, thereby improving the guidance effect for fracturing development of horizontal wells in unconventional reservoirs.

[0062] See Figure 1 As shown in the figure, an embodiment of the present invention discloses a method for identifying the fracture dip angle of a horizontal well through a fracture, including:

[0063] Step S11: Based on the preset array acoustic logging tool, the wellbore mode wave is acquired in the horizontal well to complete the array logging operation and obtain the array logging results.

[0064] In this embodiment, the array waveform under fractured formation is first obtained through wellbore mode wave identification. Specifically, a preset array acoustic logging tool is used to acquire wellbore mode waves in fractured formations near horizontal wells to determine the array logging results. The array logging results include the array waveform under the fractured formation; the array waveform is an orthogonal dipole shear wave waveform. The preset array acoustic logging tool can be customized or updated based on actual needs.

[0065] It is understandable that array waveforms refer to orthogonal dipole shear wave waveforms measured in fractured formations. The waveforms mainly include the direct wave signal of the dipole shear wave and the reflected wave signal excited by the well fracture.

[0066] Regarding well-through fractures, combined with Figure 2 As shown, a through-well fracture refers to a fracture that is penetrated by the wellbore. If the fracture is considered as a thin plate structure filled with fluid, the wellbore penetrating the fracture surface forms a "through-well" state. Here, through-well fractures include not only fractures induced during drilling but also pre-existing natural fractures in the formation.

[0067] Step S12: Based on the preset array waveform processing strategy, perform dipole shear wave signal extraction on the array waveform in the array logging results to determine the dipole shear wave signal extraction result.

[0068] In this embodiment, after determining the array logging results, the array waveforms in the array logging results are processed to extract the direct wave and reflected wave signals of the dipole shear wave. Specifically: wave signal composition analysis is performed on the orthogonal dipole shear wave waveforms in the array logging results to determine the analysis results; based on the analysis results, wavefield separation is performed on the orthogonal dipole shear wave waveforms to determine the wavefield separation results; based on the wavefield separation results, direct wave signals are extracted to determine the direct wave signal extraction results; based on the wavefield separation results, reflected wave signals are extracted to determine the reflected wave signal extraction results; the reflected wave signal extraction results include reflected wave signals excited by the through-well fractures corresponding to the horizontal well; based on the direct wave signal extraction results and the reflected wave signal extraction results, the dipole shear wave signal extraction results are determined. In other words, in this embodiment, processing the array waveform refers to performing wavefield separation on the entire dipole shear wave waveform.

[0069] Furthermore, regarding wavefield separation, in this embodiment, the apparent velocity of the direct wave signal and the reflected wave signal in the array logging results is first analyzed to determine the apparent velocity analysis result. Then, based on the apparent velocity analysis result and a preset filtering algorithm, wavefield separation is performed on the direct wave signal and the reflected wave signal in the array logging results to determine the wavefield separation result. The preset filtering algorithm can be customized or updated based on actual needs, including but not limited to FK filtering (Fourier-Krimsky filtering). In other words, in this embodiment, wavefield separation refers to using a corresponding filtering algorithm to separate the direct wave and the reflected wave based on the apparent velocity difference, thereby extracting the direct wave and reflected wave signals from the dipole shear wave.

[0070] Step S13: Based on the extraction results of the dipole shear wave signal, perform relative amplitude analysis of the reflected waves of various shear waves, and use the corresponding relative amplitude analysis results to determine the target fracture dip angle range corresponding to the wellbore fracture of the horizontal well; wherein, the wellbore fracture is a fracture that is penetrated by the wellbore of the horizontal well; the various shear waves include horizontally polarized shear waves and vertically polarized shear waves corresponding to the dipole shear wave.

[0071] In this embodiment, after determining the results of the dipole shear wave signal extraction, the relative amplitudes of the reflected waves in two modes, SV shear wave (Vertically Polarized Shear Wave) and SH shear wave (Horizontally Polarized Shear Wave), are calculated, and the range of crack dip angle is determined using the relative amplitudes of the reflected waves. That is: based on the dipole shear wave signal extraction results and preset amplitude statistical rules, amplitude values ​​of shear waves in various modes are statistically analyzed to determine the target amplitude value statistical results; the target amplitude value statistical results include direct wave amplitude value statistical results and reflected wave amplitude value statistical results; based on the direct wave amplitude value statistical results and the reflected wave amplitude value statistical results, relative amplitude analysis of the reflected waves is performed to determine the first relative amplitude analysis result corresponding to the horizontally polarized shear wave and the second relative amplitude analysis result corresponding to the vertically polarized shear wave; it is determined whether the first relative amplitude analysis result and the relative amplitude analysis result are less than a first preset amplitude threshold to determine the first threshold judgment result; it is determined whether the difference between the first relative amplitude analysis result and the relative amplitude analysis result is less than a second preset amplitude threshold to determine the second threshold judgment result; based on the first threshold judgment result, the second threshold judgment result, and preset wellbore fracture dip angle range configuration information, the target fracture dip angle range corresponding to the wellbore fracture of the horizontal well is determined. The first preset amplitude threshold and the second preset amplitude threshold can be configured or updated based on actual needs.

[0072] Furthermore, regarding the determination of the target amplitude value statistical results, in this embodiment: based on the direct wave signal extraction results and the acoustic system structure parameters corresponding to the preset array acoustic logging tool and the first preset time window, amplitude value statistics are performed to determine the statistical results of the first direct wave amplitude value corresponding to the horizontally polarized shear wave; based on the direct wave signal extraction results and the acoustic system structure parameters corresponding to the preset array acoustic logging tool and the second preset time window, amplitude value statistics are performed to determine the statistical results of the second direct wave amplitude value corresponding to the vertically polarized shear wave; based on the reflected wave signal extraction results and the acoustic system structure parameters corresponding to the preset array acoustic logging tool and the third preset time window, amplitude value statistics are performed to determine the statistical results of the first reflected wave amplitude value corresponding to the horizontally polarized shear wave; based on the reflected wave signal extraction results and the acoustic system structure parameters corresponding to the preset array acoustic logging tool and the fourth preset time window, amplitude value statistics are performed to determine the statistical results of the second reflected wave amplitude value corresponding to the vertically polarized shear wave. The acoustic system structural parameters include, but are not limited to, minimum source distance and receiver distance (analogous to minimum shot-receiver distance and spacing in earthquakes).

[0073] It is understandable that the process of calculating the relative amplitude of the reflected wave mainly includes three steps: direct wave amplitude extraction, reflected wave amplitude extraction, and relative amplitude calculation. Among them, direct wave amplitude extraction refers to setting a specified time window for the pure dipole shear wave direct wave waveform to statistically analyze its amplitude value, and the calculation formula is shown in Equation (1); reflected wave amplitude extraction refers to setting a specified time window for the reflected wave signal after removing the direct wave to perform amplitude statistics, and the formula is the same as Equation (1).

[0074] (1).

[0075] In the formula: The amplitude represents the waveform statistics; express Time waveform amplitude; subscript Indicates the wave type, DSV, DSH, RSV, and RSH can be used. DSV / DSH refers to the direct shear wave of SV / SH, and RSV / RSH refers to the reflected shear wave of SV / SH. Indicates the start time of the calculation. This indicates the end time of the calculation.

[0076] The relative amplitude calculation uses the amplitudes of the direct wave and the reflected wave as inputs, and calculates the relative amplitude of the SV reflected wave according to equation (2). Relative amplitude of SH reflected wave .

[0077] (2).

[0078] Determining the crack dip angle range using the relative amplitude of the reflected wave refers to using... and The differential characteristics of the crack dip angle response are used to classify the range of crack dip angles. In a specific implementation, the specific process for classifying the dip angle range is as follows:

[0079] For unknown tilt angles ( The relative amplitude of reflected waves from well fractures (unit: degrees) is set to a threshold value. , :

[0080] 1) If and All less than If so, the well-pass fracture is determined to be located in the low-angle zone. ≤45), the wellbore fracture is a low-angle fracture, that is, the target fracture dip angle range is 0 to 45 degrees;

[0081] 2) If and All are greater than or equal to , and | - |> If the fracture is located in the high dip angle zone (45 < 0.05°), then it is determined that the fracture is located in the high dip angle zone. <90), the fracture through the well is a high-angle fracture, that is, the target fracture dip angle range is 45 to 90 degrees;

[0082] 3) If and All are greater than or equal to , and | - |≤ If so, the fracture through the well is determined to be a vertical fracture. =90), the well-pass fracture is a vertical fracture, that is, the target fracture dip angle range is 90 degrees.

[0083] Step S14: Based on the preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results, perform quantitative inversion to complete the well-pass fracture dip angle identification operation and obtain the fracture dip angle identification result corresponding to the well-pass fracture.

[0084] In this embodiment, after determining the target fracture dip angle range, the fracture dip angle is calculated using a preset piecewise inversion model based on this dip angle range. That is: obtaining the preset piecewise inversion model; determining the target inversion model corresponding to the well-passing fracture based on the target fracture dip angle range and the preset piecewise inversion model; performing quantitative inversion based on the relative amplitude analysis results and the target inversion model to complete the well-passing fracture dip angle identification operation and obtain the fracture dip angle identification result corresponding to the well-passing fracture.

[0085] It should be understood that calculating the crack dip angle using the preset segmented inversion model means that different inversion models are used according to the target crack dip angle range predicted in step S13. The preset segmented inversion model can be shown in the following equation (3):

[0086] (3).

[0087] In the formula, This represents a low-angle inversion model, specifically in the form of a logarithmic equation. , among them and These are the equation parameters; This represents a high-angle inversion model, specifically in the form of linear equations. , among them and These are the equation parameters. It's understandable that if the target crack dip angle ranges from 0 to 45 degrees, then [the appropriate parameter] is selected. The target inversion model is selected, and quantitative inversion is performed; if the target crack dip angle ranges from 45 to 90 degrees, then... The target inversion model is used, and quantitative inversion is performed; if the target fracture dip angle range is 90 degrees, then the fracture dip angle of the corresponding well-passing fracture is determined to be 90 degrees.

[0088] In summary, this embodiment considers that borehole mode waves are extremely sensitive to the response of fractures near the well, especially dipole shear waves. The SH and SV waves of dipole shear waves follow different reflection coefficient equations when incident on the inclined solid-liquid interface, and their responses to changes in the inclination angle of the fracture passing through the well differ significantly, particularly in horizontal wells. Therefore, to address these differences in response characteristics, a scheme for identifying the inclination angle of fractures passing through horizontal wells based on borehole mode waves is proposed, aiming to achieve quantitative inversion of the fracture inclination angle. In this scheme, the difference in the response characteristics of SH and SV reflected waves to fractures passing through the well is used as the input for inverting the fracture inclination angle, and a piecewise inversion model based on the relative amplitudes of SH and SV reflected waves is proposed. This improves the accuracy and reliability of identifying the fracture inclination angle of fractures passing through horizontal wells and reduces the identification cost.

[0089] Therefore, this application first obtains the array logging results of a horizontal well, then extracts the dipole shear wave signal from the array logging results, and uses the corresponding dipole shear wave signal extraction results to perform relative amplitude analysis of various shear wave reflections to determine the relative amplitude analysis results. Subsequently, using the relative amplitude analysis results, the target fracture dip angle range of the through-well fracture is determined, and quantitative inversion is performed by combining a pre-set segmented inversion model and the relative amplitude analysis results to obtain the fracture dip angle identification result corresponding to the through-well fracture. This approach effectively and accurately identifies the fracture dip angle of through-well fractures in horizontal wells, reduces identification costs, and thus improves the guidance effect for fracturing development work of horizontal wells in unconventional reservoirs.

[0090] The following is combined with Figures 2 to 7 The schematic diagram disclosed herein provides a detailed description of the technical solutions of the embodiments of this application.

[0091] In one specific implementation, array waveforms were first acquired in the fractured formation, as shown in Figures 3(a) and (b), illustrating the orthogonal dipole shear wave waveforms measured in the fractured formation. Figure 3(a) shows the SV shear wave array received waveform, with the horizontal axis representing time in milliseconds (ms) and the vertical axis representing the distance between the source and receiver in meters (m). The horizontal solid line indicates the location of the fracture. It can be observed that the SV shear wave exhibits significant reflection at the fracture location. The arrival time of the SV reflected wave increases with decreasing source distance (marked by a dashed rectangle), while the arrival time of the SV direct wave decreases with decreasing source distance (marked by a solid rectangle). Figure 3(b) shows the SH shear wave array received waveform. Since the formation outside the well is a homogeneous isotropic medium, the travel time patterns of the SH shear wave's direct and reflected waves are consistent with those of the SV shear wave. Further observation of the reflected wave amplitude reveals that, under the same source distance receiving conditions, the amplitude of the SV reflected wave is lower than that of the SH reflected wave.

[0092] Next, the direct and reflected wave signals of the dipole shear wave are extracted from the processed array waveform. Figures 4(a) and (b) show the dipole shear wave array waveforms after wavefield separation processing. As shown, the FK filtering method effectively separates the direct and reflected waves. Compared with the original full waveform in Figure 3, only the weaker reflected wave signal is retained in Figure 4, and the direct wave component is significantly suppressed.

[0093] Next, the relative amplitudes of the reflected waves in the SV and SH modes were calculated. Figure 5 The results of relative amplitude analysis of the dipole shear wave signal after wavefield separation and the reflected waves of the two modes are presented. Since the formation outside the well is a homogeneous isotropic medium, the direct wave amplitudes of the SV and SH shear waves are shown. They are completely identical, both at 54.824 mV, but the amplitude of the reflected wave of the SV transverse wave ( =9.275mV), lower than the amplitude of the reflected wave of the SH transverse wave ( =19.901mV), that is, 9.275 < 19.901, which in turn leads to the relative reflection amplitude of the reflected wave of the SV transverse wave ( =0.171), which is also lower than the relative reflection amplitude of the reflected wave of the SH reflected transverse wave ( =0.363), that is, 0.171 < 0.363.

[0094] Furthermore, the range of crack dip angles is determined by utilizing the relative amplitude of the reflected waves. Figure 6 This paper demonstrates the specific process for determining the fracture dip angle range using the relative amplitudes of SV and SH reflected waves. An amplitude threshold is set based on the elastic parameters of the rock in this block. =0.03、 =0.01. As shown in the figure, firstly, because and All greater than Therefore, it can be determined that the wellbore fracture is not a low-angle fracture; secondly, because | - |=0.192≥ Therefore, it was determined that the dip angle of the well-passing fracture was located in the high dip angle region and belonged to the high dip angle fracture.

[0095] Finally, the crack dip angle was calculated using a piecewise inversion model. Figure 7 The diagram shows the inversion chart used to quantitatively calculate the crack dip angle using a pre-defined piecewise inversion model. The horizontal axis in the chart represents the relative reflection amplitude of the crack's SH-reflected wave. The vertical axis represents the crack dip angle, in degrees. The dashed line in the figure represents the low-angle inversion model. Specifically, it takes the form of a logarithmic equation. Furthermore, the solid line in the figure represents the high-angle inversion model. Specifically, it takes the form of a linear equation. Based on the calibration results of the rock elastic parameters in this block, the specific parameter values ​​for the inversion model are as follows: Low dip angle zone High-angle region , .

[0096] By substituting the corresponding values ​​into the inversion model The value can be used to determine the fracture dip angle of the corresponding wellbore fracture. Based on the aforementioned steps, it is known that the fracture dip angle to be inverted is located in a high dip angle region; therefore, a high-angle inversion model should be used, such as... Figure 7 As shown by the arrow in the image, substitute... =0.363, the calculated inclination angle of the fracture through the well is 0.363. .

[0097] See Figure 8 As shown in the figure, this application also discloses a fracture inclination angle identification device for horizontal well crossing fractures, including:

[0098] The array logging module 11 is used to acquire borehole mode waves in a horizontal well based on a preset array acoustic logging instrument, so as to complete the array logging operation and obtain the array logging results.

[0099] The shear wave signal extraction module 12 is used to extract dipole shear wave signals from the array waveform in the array logging results based on a preset array waveform processing strategy, so as to determine the dipole shear wave signal extraction result.

[0100] The relative amplitude analysis module 13 is used to perform relative amplitude analysis of the reflected waves of various shear waves based on the extraction results of the dipole shear wave signal, and to determine the target fracture dip angle range corresponding to the wellbore fracture of the horizontal well using the corresponding relative amplitude analysis results; wherein, the wellbore fracture is a fracture that is penetrated by the wellbore of the horizontal well; the various shear waves include horizontally polarized shear waves and vertically polarized shear waves corresponding to the dipole shear wave;

[0101] The quantitative inversion module 14 is used to perform quantitative inversion based on the preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results, so as to complete the well-pass fracture dip angle identification operation and obtain the fracture dip angle identification result corresponding to the well-pass fracture.

[0102] In some specific embodiments, the array logging module 11 may specifically include:

[0103] The wellbore pattern wave identification unit is used to acquire a preset array acoustic logging tool to collect wellbore pattern waves in fractured formations near horizontal wells in order to determine the array logging results; the array logging results include array waveforms under the fractured formations; the array waveforms are orthogonal dipole shear wave waveforms.

[0104] In some specific embodiments, the shear wave signal extraction module 12 may specifically include:

[0105] The signal analysis unit is used to perform wave signal composition analysis on the orthogonal dipole shear wave waveform in the array logging results to determine the analysis results;

[0106] The wave field separation unit is used to perform wave field separation on the orthogonal dipole shear wave waveform based on the analysis results, so as to determine the wave field separation result;

[0107] The direct wave extraction unit is used to extract the direct wave signal based on the wave field separation result, so as to determine the direct wave signal extraction result;

[0108] The reflected wave extraction unit is used to extract the reflected wave signal based on the wave field separation result, so as to determine the reflected wave signal extraction result; the reflected wave signal extraction result includes the reflected wave signal excited by the well-pass fracture corresponding to the horizontal well;

[0109] The extraction result determination unit is used to determine the dipole transverse wave signal extraction result based on the direct wave signal extraction result and the reflected wave signal extraction result.

[0110] In some specific embodiments, the wavefield separation unit may specifically include:

[0111] The apparent velocity analysis subunit is used to perform apparent velocity analysis on the direct wave signal and reflected wave signal in the array logging results to determine the apparent velocity analysis results;

[0112] Based on the apparent velocity analysis results and a preset filtering algorithm, the wavefield separation subunit performs wavefield separation on the direct wave signal and reflected wave signal in the array logging results to determine the wavefield separation result.

[0113] In some specific embodiments, the relative amplitude analysis module 13 may specifically include:

[0114] An amplitude value statistics unit is used to perform amplitude value statistics for various modes of shear waves based on the dipole shear wave signal extraction results and preset amplitude statistics rules, so as to determine the target amplitude value statistics results; the target amplitude value statistics results include direct wave amplitude value statistics results and reflected wave amplitude value statistics results;

[0115] The relative amplitude analysis unit is used to perform relative amplitude analysis of the reflected wave based on the statistical results of the direct wave amplitude value and the statistical results of the reflected wave amplitude value, so as to determine the first relative amplitude analysis result corresponding to the horizontally polarized transverse wave and the second relative amplitude analysis result corresponding to the vertically polarized transverse wave.

[0116] The first threshold judgment unit is used to determine whether the first relative amplitude analysis result and the relative amplitude analysis result are less than a first preset amplitude threshold, so as to determine the first threshold judgment result;

[0117] The second threshold judgment unit is used to determine whether the difference between the first relative amplitude analysis result and the relative amplitude analysis result is less than the second preset amplitude threshold, so as to determine the second threshold judgment result.

[0118] The dip angle range determination unit is used to determine the target fracture dip angle range corresponding to the cross-well fracture of the horizontal well based on the first threshold judgment result, the second threshold judgment result, and the preset cross-well fracture dip angle range configuration information.

[0119] In some specific embodiments, the amplitude value statistics unit may specifically include:

[0120] The first amplitude value statistics subunit is used to perform amplitude value statistics based on the direct wave signal extraction result, the acoustic system structure parameters corresponding to the preset array acoustic logging tool, and the first preset time window, so as to determine the first direct wave amplitude value statistics result corresponding to the horizontally polarized shear wave.

[0121] The second amplitude value statistics subunit is used to perform amplitude value statistics based on the direct wave signal extraction result and the acoustic system structure parameters and the second preset time window corresponding to the preset array acoustic logging tool, so as to determine the second direct wave amplitude value statistics result corresponding to the vertically polarized shear wave.

[0122] The third amplitude value statistics subunit is used to perform amplitude value statistics based on the reflected wave signal extraction result, the acoustic system structure parameters corresponding to the preset array acoustic logging tool, and the third preset time window, so as to determine the first reflected wave amplitude value statistics result corresponding to the horizontally polarized shear wave.

[0123] The fourth amplitude value statistics unit is used to perform amplitude value statistics based on the reflected wave signal extraction result, the acoustic system structure parameters corresponding to the preset array acoustic logging tool, and the fourth preset time window, so as to determine the second reflected wave amplitude value statistics result corresponding to the vertically polarized transverse wave.

[0124] In some specific embodiments, the quantitative inversion module 14 may specifically include:

[0125] The model acquisition unit is used to acquire a preset segmented inversion model;

[0126] The target model determination unit is used to determine the target inversion model corresponding to the well-through fracture based on the target fracture dip angle range and the preset segmented inversion model;

[0127] The quantitative inversion unit is used to perform quantitative inversion based on the relative amplitude analysis results and the target inversion model to complete the well-pass fracture dip angle identification operation and obtain the fracture dip angle identification result corresponding to the well-pass fracture.

[0128] Furthermore, embodiments of this application also disclose an electronic device, Figure 9 This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content of the diagram should not be construed as limiting the scope of this application.

[0129] Figure 9 This is a schematic diagram of the structure of an electronic device 20 provided in an embodiment of this application. Specifically, the electronic device 20 may include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the fracture inclination angle identification method for horizontal well through-well fractures disclosed in any of the foregoing embodiments. Alternatively, the electronic device 20 in this embodiment may specifically be an electronic computer.

[0130] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.

[0131] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored thereon can include operating system 221, computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0132] The operating system 221 is used to manage and control the various hardware devices on the electronic device 20 and the computer program 222, which may be Windows Server, Netware, Unix, Linux, etc. In addition to including a computer program capable of performing the fracture inclination angle identification method for horizontal well crossing fractures disclosed in any of the foregoing embodiments, the computer program 222 may further include computer programs capable of performing other specific tasks.

[0133] Furthermore, this application also discloses a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, it implements the aforementioned method for identifying the fracture inclination angle of horizontal well crossing fractures. Specific steps of this method can be found in the corresponding content disclosed in the foregoing embodiments, and will not be repeated here.

[0134] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from other embodiments. Similar or identical parts between embodiments can be referred to interchangeably. For the apparatus disclosed in the embodiments, since it corresponds to the method disclosed in the embodiments, the description is relatively simple; relevant parts can be referred to in the method section.

[0135] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0136] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

[0137] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0138] The technical solutions provided in this application have been described in detail above. Specific examples have been used to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for identifying fracture dip of a fracture crossing a horizontal well, the method comprising: determining a first fracture dip angle; determining a second fracture dip angle; and determining a third fracture dip angle. The method comprises the following steps: Based on a preset array acoustic logging instrument, borehole mode waves of a horizontal well are collected to complete array logging operation and obtain array logging results; Based on a preset array waveform processing strategy, dipole shear wave signals are extracted from array waveforms in the array logging results to determine dipole shear wave signal extraction results; Based on the dipole shear wave signal extraction results, reflection wave relative amplitude analysis of multiple shear waves is performed, and corresponding relative amplitude analysis results are used to determine a target fracture dip angle range corresponding to a through-hole fracture of the horizontal well; the through-hole fracture is a fracture penetrated by a borehole of the horizontal well; the multiple shear waves include horizontal polarization shear waves and vertical polarization shear waves corresponding to dipole shear waves; Based on a preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results, quantitative inversion is performed to complete through-hole fracture dip angle identification operation and obtain fracture dip angle identification results corresponding to the through-hole fracture.

2. The method of claim 1, wherein, The method comprises the following steps: Based on a preset array acoustic logging instrument, borehole mode waves of a horizontal well are collected to complete array logging operation and obtain array logging results; 3. The method of claim 2, wherein, Based on a preset array waveform processing strategy, dipole shear wave signals are extracted from array waveforms in the array logging results to determine dipole shear wave signal extraction results; Based on the dipole shear wave signal extraction results, reflection wave relative amplitude analysis of multiple shear waves is performed, and corresponding relative amplitude analysis results are used to determine a target fracture dip angle range corresponding to a through-hole fracture of the horizontal well; the through-hole fracture is a fracture penetrated by a borehole of the horizontal well; the multiple shear waves include horizontal polarization shear waves and vertical polarization shear waves corresponding to dipole shear waves; Based on a preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results, quantitative inversion is performed to complete through-hole fracture dip angle identification operation and obtain fracture dip angle identification results corresponding to the through-hole fracture. The method comprises the following steps: Based on a preset array acoustic logging instrument, borehole mode waves of a horizontal well are collected to complete array logging operation and obtain array logging results; Based on a preset array waveform processing strategy, dipole shear wave signals are extracted from array waveforms in the array logging results to determine dipole shear wave signal extraction results; 4. The method of claim 3, wherein, Based on the dipole shear wave signal extraction results, reflection wave relative amplitude analysis of multiple shear waves is performed, and corresponding relative amplitude analysis results are used to determine a target fracture dip angle range corresponding to a through-hole fracture of the horizontal well; the through-hole fracture is a fracture penetrated by a borehole of the horizontal well; the multiple shear waves include horizontal polarization shear waves and vertical polarization shear waves corresponding to dipole shear waves; Based on a preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results, quantitative inversion is performed to complete through-hole fracture dip angle identification operation and obtain fracture dip angle identification results corresponding to the through-hole fracture. The method comprises the following steps:

5. The method of claim 3, wherein, Based on a preset array acoustic logging instrument, borehole mode waves of a horizontal well are collected to complete array logging operation and obtain array logging results; Based on a preset array waveform processing strategy, dipole shear wave signals are extracted from array waveforms in the array logging results to determine dipole shear wave signal extraction results; Based on the dipole shear wave signal extraction results, reflection wave relative amplitude analysis of multiple shear waves is performed, and corresponding relative amplitude analysis results are used to determine a target fracture dip angle range corresponding to a through-hole fracture of the horizontal well; the through-hole fracture is a fracture penetrated by a borehole of the horizontal well; the multiple shear waves include horizontal polarization shear waves and vertical polarization shear waves corresponding to dipole shear waves; Based on a preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results, quantitative inversion is performed to complete through-hole fracture dip angle identification operation and obtain fracture dip angle identification results corresponding to the through-hole fracture. The method comprises the following steps: Based on a preset array acoustic logging instrument, borehole mode waves of a horizontal well are collected to complete array logging operation and obtain array logging results; Based on a preset array waveform processing strategy, dipole shear wave signals are extracted from array waveforms in the array logging results to determine dipole shear wave signal extraction results; Based on the dipole shear wave signal extraction results, reflection wave relative amplitude analysis of multiple shear waves is performed, and corresponding relative amplitude analysis results are used to determine a target fracture dip angle range corresponding to a through-hole fracture of the horizontal well; the through-hole fracture is a fracture penetrated by a borehole of the horizontal well; the multiple shear waves include horizontal polarization shear waves and vertical polarization shear waves corresponding to dipole shear waves; Based on a preset segmented inversion model, the target fracture dip angle range, and the relative amplitude analysis results, quantitative inversion is performed to complete through-hole fracture dip angle identification operation and obtain fracture dip angle identification results corresponding to the through-hole fracture. Based on the dipole shear wave signal extraction result and a preset amplitude statistical rule, amplitude value statistics of multiple modes of shear waves are respectively performed to determine a target amplitude value statistical result; the target amplitude value statistical result includes a direct wave amplitude value statistical result and a reflected wave amplitude value statistical result; Based on the direct wave amplitude value statistical result and the reflected wave amplitude value statistical result, relative amplitude analysis of the reflected wave is performed to determine a first relative amplitude analysis result corresponding to a horizontally polarized shear wave and a second relative amplitude analysis result corresponding to a vertically polarized shear wave; Whether the first relative amplitude analysis result and the relative amplitude analysis result are smaller than a first preset amplitude threshold is respectively judged to determine a first threshold judgment result; Whether a difference between the first relative amplitude analysis result and the relative amplitude analysis result is smaller than a second preset amplitude threshold is judged to determine a second threshold judgment result; Based on the first threshold judgment result, the second threshold judgment result and preset through-hole fracture dip angle range configuration information, a target fracture dip angle range corresponding to a through-hole fracture of the horizontal well is determined.

6. The method of fracture dip identification for a horizontal well crossing a fracture according to claim 5, wherein, The based on the dipole shear wave signal extraction result and a preset amplitude statistical rule, amplitude value statistics of multiple modes of shear waves are respectively performed to determine a target amplitude value statistical result, includes: Based on the direct wave signal extraction result and a first preset time window corresponding to a sound system structure parameter of the preset array acoustic logging instrument, amplitude value statistics are performed to determine a first direct wave amplitude value statistical result corresponding to a horizontally polarized shear wave; Based on the direct wave signal extraction result and a second preset time window corresponding to the sound system structure parameter of the preset array acoustic logging instrument, amplitude value statistics are performed to determine a second direct wave amplitude value statistical result corresponding to a vertically polarized shear wave; Based on the reflected wave signal extraction result and a third preset time window corresponding to the sound system structure parameter of the preset array acoustic logging instrument, amplitude value statistics are performed to determine a first reflected wave amplitude value statistical result corresponding to the horizontally polarized shear wave; Based on the reflected wave signal extraction result and a fourth preset time window corresponding to the sound system structure parameter of the preset array acoustic logging instrument, amplitude value statistics are performed to determine a second reflected wave amplitude value statistical result corresponding to the vertically polarized shear wave.

7. The method of claim 1 to 6, wherein, The based on the preset segmented inversion model, the target fracture dip angle range and the relative amplitude analysis result to perform quantitative inversion to complete the through-hole fracture dip angle identification operation, includes: A preset segmented inversion model is acquired; Based on the target fracture dip angle range and the preset segmented inversion model, a target inversion model corresponding to the through-hole fracture is determined; Based on the relative amplitude analysis result and the target inversion model, quantitative inversion is performed to complete the through-hole fracture dip angle identification operation and obtain a fracture dip angle identification result corresponding to the through-hole fracture.

8. A fracture dip identification device for horizontal well crossing fractures, characterized in that, It includes: An array logging module is configured to perform borehole mode wave acquisition on a horizontal well based on a preset array acoustic logging instrument to complete array logging operation and obtain array logging results. The transverse wave signal extraction module is configured to perform dipole transverse wave signal extraction on the array waveforms in the array logging result based on a preset array waveform processing strategy, so as to determine a dipole transverse wave signal extraction result. The relative amplitude analysis module is configured to perform relative amplitude analysis on reflected waves of multiple transverse waves based on the dipole transverse wave signal extraction result, and determine a target fracture dip angle range corresponding to a fracture penetrated by the horizontal well based on corresponding relative amplitude analysis results. The quantitative inversion module is configured to perform quantitative inversion based on a preset piecewise inversion model, the target fracture dip angle range, and the relative amplitude analysis result, so as to complete a fracture dip angle identification operation and obtain a fracture dip angle identification result corresponding to the fracture penetrated by the horizontal well.

9. An electronic device, comprising: The memory is configured to store a computer program. The processor is configured to execute the computer program to implement the fracture dip angle identification method for a fracture penetrated by a horizontal well according to any one of claims 1 to 7. The memory is configured to store a computer program.

10. A computer-readable storage medium, characterized in that, The memory is configured to store a computer program.

Images