Method, apparatus, and article of manufacture for identifying fractures near a well based on dipole array waveform data

Abstract

The application discloses a method, device and equipment for identifying a fracture beside a well based on dipole array waveform data, and the method comprises the following steps: preprocessing dipole array waveform data of a target well to obtain sliding transverse wave array waveform data of the target well; extracting a transverse wave time difference curve and a transverse wave arrival time curve from the sliding transverse wave array waveform data; suppressing interference signals of the sliding transverse wave array waveform data based on the transverse wave time difference curve and / or the transverse wave arrival time curve to obtain a sliding reverse transmission transverse wave signal of the target well; separating and imaging the sliding reverse transmission transverse wave signal in a frequency-wave number domain to obtain a sliding reverse transmission transverse wave fracture imaging graph of the target well; and identifying a fracture beside the target well based on the sliding reverse transmission transverse wave fracture imaging graph. The method can improve the resolution and identification accuracy of fracture identification, and provides technical support for subsequent perforation and fracturing of the target well.

Description

Technical Field

[0001] This invention relates to the field of oil and gas reservoir exploration technology, and in particular to a method, apparatus and equipment for identifying well-side fractures based on dipole array waveform data. Background Technology

[0002] Carbonate reservoirs are a major replacement area for onshore oil reserves and a primary target for reserve growth in my country. These reservoirs are characterized by strong heterogeneity, low matrix porosity, and widespread fractures and caverns. Fractures in these reservoirs not only serve as channels for fluid accumulation and migration but also control the development of pores and caverns, the distribution of original fluids within the formation, and the characteristics of mud intrusion, significantly contributing to reservoir permeability. Therefore, accurately identifying and quantitatively evaluating well-side fracture development using array sonic logging data is of great significance for carbonate oil and gas exploration and development.

[0003] Morris et al. (1964) first studied the relationship between P-wave amplitude and fracture width in acoustic logging data, finding that the attenuation of P-waves when passing through horizontal fractures was very small and could be ignored. Furthermore, they found that the transmission coefficient of S-waves was affected by the fracture dip angle, and that S-wave attenuation was more pronounced at horizontal fractures. Paillet (1980) and Zlatev et al. (1988) studied the attenuation phenomenon of P-waves and S-waves at fractures using array acoustic logging data and physical experiments, respectively. Chen (1994) studied the propagation of borehole mode waves in isotropic formations based on a three-dimensional finite difference algorithm of second-order time and fourth-order space. The study showed that when the dipole source was parallel to the formation interface, the formation interface had little influence on the waveform; however, when the dipole was perpendicular to the formation interface, strong dipole S-wave interference existed. Zhu et al. (1994) used aluminum and fluorite borehole models to study the propagation of dipole S-waves in fractures, finding that dipole S-waves were more sensitive to horizontal fractures than vertical fractures. Haldorsen et al. (2006) studied the propagation of P-waves and S-waves in irregular wellbores and anisotropic formations, respectively.

[0004] The above research mainly focuses on how to use wellbore mode waves, such as the velocity and attenuation properties of longitudinal and transverse direct waves, to distinguish well-side fractures. Summary of the Invention

[0005] To improve the accuracy of wellbore fracture identification, thereby enriching process routes and increasing the selection space, this invention provides a method, apparatus, and device for identifying wellbore fractures based on dipole array waveform data.

[0006] In a first aspect, embodiments of the present invention provide a method for identifying well-side fractures based on dipole array waveform data, which may include:

[0007] The dipole array waveform data of the target well is preprocessed to obtain the sliding shear wave array waveform data of the target well;

[0008] Extract the shear wave time difference curve and the shear wave arrival time curve from the waveform data of the gliding shear wave array;

[0009] Based on the shear wave time difference curve and / or the shear wave arrival time curve, interference signal suppression is applied to the shear wave array waveform data to obtain the shear wave signal of the target well.

[0010] Based on the sliding reverse propagation shear wave signal, separation imaging is performed in the frequency-wavenumber domain to obtain the sliding reverse propagation shear wave fracture imaging map of the target well.

[0011] Based on the sliding reverse transverse wave fracture imaging map, the well-side fractures of the target well are identified.

[0012] Optionally, based on the shear wave time difference curve and / or the shear wave arrival time curve, interference signal suppression is applied to the sliding shear wave array waveform data to obtain the sliding reverse propagation shear wave signal of the target well, which may include:

[0013] Based on the shear wave time difference curve and the shear wave arrival time curve, the shear wave direct wave suppression is performed on the gliding shear wave array waveform data.

[0014] Based on the shear wave time difference curve, continuous oscillating wave suppression is applied to the gliding shear wave array waveform data after shear wave direct wave suppression.

[0015] Based on the shear wave time difference curve, scattered noise suppression is performed on the gliding shear wave array waveform data after continuous oscillation wave suppression.

[0016] Based on the shear wave time difference curve, incoherent noise suppression is performed on the scatter noise-suppressed shear wave array waveform data to obtain the sliding reverse propagation shear wave signal of the target well.

[0017] Optionally, methods for suppressing shear wave direct waves may include shear wave variation and Radon domain filtering variation;

[0018] Methods for suppressing continuous oscillating waves may include: two-dimensional filtering, wavelet transform frequency division denoising, or SRME multiple wave attenuation;

[0019] Methods for suppressing scattered noise may include: inverse Q filtering and linear amplitude enhancement;

[0020] Methods for suppressing incoherent noise may include: superposition of common center points and superposition of common reflection points.

[0021] Optionally, the preprocessing of the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well may include:

[0022] Gain recovery is performed on the dipole array waveform data of the target well based on the dipole array acoustic gain curve of the target well.

[0023] For the dipole array waveform data after gain recovery, delay recovery is performed based on the dipole array acoustic delay curve of the target well so that the dipole array waveform data is padded with zeros before time zero.

[0024] Bandpass filtering is performed on the delayed-recovered dipole array waveform data to obtain the sliding shear wave array waveform data of the target well.

[0025] Optionally, the step of extracting the shear wave time difference curve and the shear wave arrival time curve from the gliding shear wave array waveform data may include:

[0026] Based on the correlation coefficient method, the shear wave time difference curve and shear wave arrival time curve are extracted from the shear wave array waveform data.

[0027] Optionally, the method may further include:

[0028] Obtain conventional logging data and array sonic logging data of the target well, and use the conventional logging data as an indicator.

[0029] The conventional logging data may include: natural gamma logging curves, caliper logging curves, or density logging curves;

[0030] Array acoustic logging data can include: original waveform curves of monopole array acoustic waves, gain curves of monopole array acoustic waves, delay curves of monopole array acoustic waves, and original waveform curves of dipole array acoustic waves, gain curves of dipole array acoustic waves, and delay curves of dipole array acoustic waves.

[0031] In a second aspect, embodiments of the present invention provide a method for constructing fracturing technology for oil and gas wells, which may include: designing perforation locations and / or fracturing process schemes based on the wellside fractures of the target well.

[0032] The wellbore fractures of the target well are identified according to the method for identifying wellbore fractures based on dipole array waveform data described in the first aspect.

[0033] Thirdly, embodiments of the present invention provide a device for identifying well-side fractures based on dipole array waveform data, which may include:

[0034] The preprocessing module is used to preprocess the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well.

[0035] The extraction module is used to extract the shear wave time difference curve and the shear wave arrival time curve from the shear wave array waveform data;

[0036] The suppression module is used to suppress interference signals on the shear wave array waveform data based on the shear wave time difference curve and / or the shear wave arrival time curve, so as to obtain the shear wave signal of the target well.

[0037] An imaging module is used to perform separation imaging in the frequency-wavenumber domain based on the sliding back-propagating shear wave signal to obtain a sliding back-propagating shear wave fracture imaging map of the target well.

[0038] The identification module is used to identify the well-side fractures of the target well based on the sliding reverse transverse wave fracture imaging map.

[0039] Fourthly, embodiments of the present invention provide a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for identifying well-side fractures based on dipole array waveform data as described in the first aspect.

[0040] Fourthly, embodiments of the present invention provide a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method for identifying well-side fractures based on dipole array waveform data as described in the first aspect.

[0041] The beneficial effects of the above-described technical solutions provided in the embodiments of the present invention include at least the following:

[0042] This invention provides a method, apparatus, and device for identifying well-side fractures based on dipole array waveform data. The overall concept of this method is to perform glide back-propagation imaging on glide shear wave array waveform data from dipole array waveform data, and then identify well-side fractures based on the glide back-propagation shear wave fracture image. This method can improve the resolution and accuracy of fracture identification, providing technical support for subsequent target well perforation and fracturing.

[0043] Other features and advantages of the invention will be set forth in the following description, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention may be realized and obtained by means of the structures particularly pointed out in the written description and the accompanying drawings.

[0044] The technical solution of the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. Attached Figure Description

[0045] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used in conjunction with embodiments of the invention to explain the invention and do not constitute a limitation thereof. In the drawings:

[0046] Figure 1 This is a flowchart of a method for identifying well-side fractures based on dipole array waveform data provided in an embodiment of the present invention;

[0047] Figure 2 This is a flowchart illustrating a detailed method for identifying well-side fractures based on dipole array waveform data, as provided in this embodiment of the invention.

[0048] Figure 3 This is an example of implementing interference signal suppression processing on the sliding transverse wave waveform provided in this embodiment of the invention;

[0049] Figure 4 This is one of the counterexamples of implementing interference signal suppression processing on the sliding transverse wave waveform provided in the embodiments of the present invention;

[0050] Figure 5 This is a second counterexample of implementing interference signal suppression processing on the sliding transverse wave waveform provided in the embodiments of the present invention;

[0051] Figure 6 This is an example of calculating the taxiing shear wave time difference and shear wave arrival time provided in an embodiment of the present invention;

[0052] Figure 7 This is an example of identifying well-side fractures using a sliding reverse-propagation shear wave fracture imaging map provided in an embodiment of the present invention;

[0053] Figure 8 This is a schematic diagram of the device for identifying well-side cracks based on dipole array waveform data provided in an embodiment of the present invention. Detailed Implementation

[0054] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0055] The inventors discovered that existing research methods essentially rely on the "convexity" and "concavity" variations of curves to identify fractures, resulting in low resolution and accuracy. In view of these problems, this invention is proposed to provide a method, apparatus, and device for identifying well-side fractures based on dipole array waveform data, which overcomes or at least partially solves these problems.

[0056] This invention provides a method for identifying well-side fractures based on dipole array waveform data, referring to... Figure 1 As shown, the method may include the following steps:

[0057] Step S11: Preprocess the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well.

[0058] Step S12: Extract the transverse wave time difference curve and transverse wave arrival time curve from the transverse wave array waveform data.

[0059] Step S13: Based on the shear wave time difference curve and / or shear wave arrival time curve, suppress the interference signal of the sliding shear wave array waveform data to obtain the sliding reverse propagation shear wave signal of the target well.

[0060] Step S14: Based on the sliding back-propagation shear wave signal, perform separation imaging in the frequency-wavenumber domain to obtain the sliding back-propagation shear wave fracture imaging map of the target well.

[0061] Step S15: Identify the well-side fractures of the target well based on the sliding reverse propagation transverse wave fracture imaging map.

[0062] The method provided in this embodiment of the invention involves performing glide back-propagation imaging on glide shear wave array waveform data from dipole array waveform data, and then identifying well-side fractures based on the glide back-propagation shear wave fracture imaging image. This method can improve the resolution and accuracy of fracture identification, providing technical support for subsequent target well perforation and fracturing.

[0063] In a detailed embodiment, refer to Figure 2 As shown, the above method for identifying well-side fractures based on dipole array waveform data may include the following steps:

[0064] Step S20: Obtain conventional logging data and array acoustic logging data for the target well, using conventional logging data as the indicator. Conventional logging data may include: natural gamma ray logging (GR) curves, caliper logging curves, or density logging curves. Array acoustic logging data may include: monopole array acoustic waveform curves, monopole array acoustic gain curves, monopole array acoustic delay curves, and dipole array acoustic waveform curves, dipole array acoustic gain curves, and dipole array acoustic delay curves. The conventional logging curves in this embodiment, such as natural gamma ray logging curves, can be used to corroborate the results and verify the accuracy of the method provided in this embodiment.

[0065] Reference Figure 3 As shown, the first channel is the gamma channel, the second channel is the depth channel, and the third channel displays the raw waveform data of the dipole shear wave (WVXXO).

[0066] Step S21: Preprocess the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well.

[0067] The specific preprocessing process is as follows: First, gain recovery is performed on the dipole array waveform data of the target well based on the dipole array acoustic gain curve of the target well; then, delay recovery is performed on the dipole array waveform data after gain recovery based on the dipole array acoustic delay curve of the target well, so that zeros are padded before the zero time of the dipole array waveform data; finally, bandpass filtering is performed on the dipole array waveform data after delay recovery to obtain the sliding shear wave array waveform data of the target well.

[0068] In practical implementation, the preprocessing of array acoustic logging data mainly includes gain recovery, delay recovery, and bandpass filtering. Gain recovery is performed first, followed by delay recovery, and finally bandpass filtering. When acquiring array acoustic waveforms downhole, automatic gain control is used to ensure the waveform amplitude is at its highest accuracy. This amplifies the waveform amplitude from a decimal to an integer for easier storage and recording. Therefore, for the acquired waveform, gain recovery is performed first. The gain parameter AGN is shown in equation (1) as follows:

[0069] AGN=10 GN*0.05 Equation (1)

[0070] During gain recovery, the array acoustic waveform is divided by AGN. To reduce the amount of data stored during waveform acquisition, waveform data for a period before the arrival of the first wave is often not acquired. Delay recovery aims to fill the array acoustic waveform with zeros before the zero point, thereby obtaining waveform data at accurate times. Bandpass filtering is used to obtain gliding shear wave data; the bandpass filter's frequency range is selected as 3000–5000 Hz.

[0071] Also refer to Figure 3 As shown, Figure 1 The third channel shows the original dipole shear wave waveform (WVXXO), and the fourth channel shows the glide shear wave array waveform data (WVXX). Based on the original dipole shear wave waveform, the glide shear wave waveform can be obtained after gain recovery, delay recovery, and bandpass filtering.

[0072] Step S22: Extract the transverse wave time difference curve and transverse wave arrival time curve from the gliding transverse wave array waveform data.

[0073] This step, based on the correlation coefficient method, extracts the shear wave time difference curve and shear wave arrival time curve from the gliding shear wave array waveform data. (Refer to...) Figure 6 As shown, Figure 6The first track is the gamma track, the second is the depth track, the third is the gliding shear wave waveform, and the fourth is the shear wave time-time-time difference correlation plot obtained using the correlation method. In the correlation plot, the horizontal axis represents the time difference value, with 0 on the left and 300 μs / ft on the right. Lighter colors in the plot represent lower correlation coefficients, and darker colors represent higher correlation coefficients. The darker bands shown in the plot indicate the location of the shear wave time difference. The fifth track is the shear wave time difference curve, and the sixth track is the shear wave arrival time curve.

[0074] In practice, this step calculates the correlation function of the array waveform in two dimensions: time and time difference. The time difference where the maximum function value is located is the transverse wave time difference. Specifically, a time window is first opened for the multi-channel waveforms, and the correlation function value of the multi-channel waveform data within the window is calculated. Then, the slope of the window is gradually changed within the slow search range (DTMIN—DTMAX), and the correlation function value is calculated again. After the entire time difference range is calculated, the time of the window is moved forward by one position, and the same time difference search calculation is performed again until the entire time range (CWBEGIN—CWEND) is calculated. The formula for calculating the correlation function Corr is shown in equation (2) as follows:

[0075]

[0076] In the formula, DT represents the time difference, Time represents the start time of the current time window, CWLENTH represents the duration of the time window, N represents the number of channels in the array waveform, and RRSP represents the distance between two adjacent receivers. Based on the shear wave time difference, and considering the wellbore diameter and logging instrument length, the shear wave arrival time curve can be directly calculated.

[0077] Step S23: Based on the shear wave time difference curve and the shear wave arrival time curve, suppress the shear wave direct wave of the gliding shear wave array waveform data.

[0078] Methods for suppressing shear waves and direct waves may include shear wave transformation and Radon domain filtering transformation.

[0079] Step S24: Based on the transverse wave time difference curve, perform continuous oscillating wave suppression on the gliding transverse wave array waveform data after the transverse wave direct wave suppression.

[0080] Methods for suppressing continuous oscillating waves may include: two-dimensional filtering, wavelet transform frequency division denoising, or SRME multiple wave attenuation.

[0081] Step S25: Based on the transverse wave time difference curve, perform scattered noise suppression on the gliding transverse wave array waveform data after continuous oscillation wave suppression.

[0082] Methods for suppressing scattered noise may include: inverse Q filtering and linear amplitude enhancement.

[0083] Step S26: Based on the shear wave time difference curve, perform incoherent noise suppression on the scatter noise-suppressed sliding shear wave array waveform data to obtain the sliding back-propagation shear wave signal of the target well. The incoherent noise suppression methods may include: common center point superposition and common reflection point superposition.

[0084] Steps S23 to S26 above are for suppressing interference signals on the shear wave array waveform data to obtain the shear wave signal of the target well. This method is actually a tiered suppression process for various types of interference signals. The four interference signal suppression techniques—shear wave direct wave suppression, continuous oscillation wave suppression, scattered noise suppression, and incoherent noise suppression—are executed sequentially to achieve the purpose of suppressing interference signals. Specifically, shear wave direct wave suppression can be achieved using signal processing methods such as shear wave transform and Radon domain filtering; continuous oscillation wave suppression can be achieved using signal processing methods such as two-dimensional filtering, wavelet transform frequency division denoising, or SRME multiple wave attenuation; scattered noise suppression can be achieved using signal processing methods such as inverse Q filtering and linear amplitude enhancement; and incoherent noise suppression is typically achieved using signal processing methods such as common center point superposition and common reflection point superposition.

[0085] Reference Figure 3 As shown, the fourth track displays the gliding shear wave waveform (WVXX); the fifth track is the waveform after filtering out the direct signal of the gliding shear wave using shear wave transform (SLXX), at which point the gliding shear wave backpropagation wave is highlighted; the sixth track is the waveform after compensating the amplitude of the far end of the gliding shear wave backpropagation wave using the amplitude linear enhancement method (GNXX), from which the far end waveform of the gliding shear wave backpropagation wave can be observed, that is, the amplitude of the waveform with a longer recording time is significantly enhanced, and the scattered noise is effectively suppressed; the seventh track is the waveform after attenuating the gliding shear wave backpropagation wave multiple times using SRME (PDXX), compared with the sixth track, the oscillation period of the gliding shear wave backpropagation wave in this track is reduced, which is manifested in the image as a thinner waveform "strip"; the eighth track is the waveform after superimposing the waveforms received by multiple receivers (CXX), at which point the uncorrelated noise is further suppressed, and the gliding shear wave backpropagation wave signal is more clearly displayed.

[0086] Step S27: Based on the sliding back-propagation shear wave signal, perform separation imaging in the frequency-wavenumber domain to obtain the sliding back-propagation shear wave fracture imaging map of the target well.

[0087] This step in the sliding back-propagating shear wave imaging primarily involves using a frequency-wavenumber domain separation method to process the sliding back-propagating shear wave, achieving separation of the ascending and descending waves. The final result is a sliding back-propagating shear wave fracture image with the ascending reflected wave on the left, the borehole in the middle, and the descending reflected wave on the right. (Refer to...) Figure 3As shown, the ninth channel is the imaging result of the glide reverse propagation shear wave crack. This is an image image after processing the up and down reflected waves, based on the glide reverse propagation shear wave waveform of the eighth channel. The left side of the image shows the up wave of the glide reverse propagation shear wave, and the right side shows the down wave of the glide reverse propagation shear wave.

[0088] Step S28: Identify the well-side fractures of the target well based on the sliding reverse propagation shear wave fracture imaging map.

[0089] In this step, refer to the imaging map of the sliding reverse transverse wave crack. Figure 3 and Figure 7 As shown, the inclined phase axis from the upper left to the lower right can be observed and labeled with a straight line. The intersection of the straight line and the well axis is the location where the fracture passes through the wellbore, and the amplitude of the inclined phase axis represents fracture parameters such as fracture opening.

[0090] Also refer to Figure 7 As shown, Figure 3 The first track is the gamma curve track, the second is the depth track, and the third is the imaging result of the sliding back-propagating shear wave through-well fracture. In the sliding back-propagating shear wave through-well fracture imaging image, the left side shows the upward sliding back-propagating shear wave, the middle vertical line represents the wellbore, and the right side shows the downward sliding back-propagating shear wave. Multiple sets of tilted, in-phase axes can be observed, all representing through-well fractures, intersecting the wellbore at depths of 5215m, 5233m, and 5243m. These three depths represent the depth at which the fracture penetrates the wellbore. The location of the fracture depth typically represents the reservoir interval, which will provide a reference for the subsequent design of perforation and fracturing schemes for this well.

[0091] It should also be noted that the tiered suppression method for multiple types of interference signals described in steps S23 to S26 of this invention is the key to interference signal processing in this embodiment of the invention. The ability to suppress multiple types of interference signals and the tiered processing are two crucial points. To illustrate the advantages of this method, examples are provided to demonstrate the role of these two key points.

[0092] Reference Figure 4 As shown, the scattered noise suppression step was missing, and subsequent processing was performed directly on the waveform SLXX (fifth channel). In the final sliding back-propagating transverse wave fracture imaging image IMAGEXX (eighth channel), there is indeed a large amount of unfiltered scattered noise around the wellbore, which seriously affects the clear display of the effective back-propagating wave signal.

[0093] Reference Figure 5As shown, the suppression sequence of the interference signal does not follow the key point of the stepwise processing method. Instead, it moves the suppression of multiple oscillating interference signals in the subsequent processing to the first step. Because the amplitude of the direct glide shear wave signal is much higher than that of the glide reverse propagation wave at this point, the main purpose of this step becomes to suppress the multiple oscillations of the direct signal, and the glide reverse propagation wave is also suppressed as noise in the process. This results in the effective glide shear wave reverse propagation wave signal being difficult to identify in the final glide reverse propagation shear wave crack imaging image IMAGEXX (ninth channel).

[0094] The above-mentioned method for identifying well-side fractures based on dipole array waveform data provided in this embodiment of the invention can effectively suppress multiple types of interference signals, thereby making the sliding transverse wave back propagation signal clearer and thus accurately identifying well-side fractures.

[0095] Based on the same inventive concept, this embodiment of the invention also provides a method for constructing oil and gas well fracturing technology. The method may include: designing perforation locations and / or fracturing process schemes based on the well-side fractures of the target well; wherein the well-side fractures of the target well are identified according to the above-mentioned method for identifying well-side fractures based on dipole array waveform data.

[0096] Based on the same inventive concept, this invention also provides a device for identifying well-side fractures based on dipole array waveform data, referring to... Figure 8 As shown, the device may include:

[0097] The preprocessing module 11 is used to preprocess the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well.

[0098] Extraction module 12 is used to extract the shear wave time difference curve and the shear wave arrival time curve from the shear wave array waveform data;

[0099] The suppression module 13 is used to suppress interference signals on the shear wave array waveform data based on the shear wave time difference curve and / or shear wave arrival time curve, so as to obtain the shear wave signal of the target well.

[0100] The imaging module 14 is used to perform separation imaging in the frequency-wavenumber domain based on the sliding back propagation shear wave signal to obtain the sliding back propagation shear wave fracture imaging map of the target well.

[0101] The identification module 15 is used to identify the well-side fractures of the target well based on the sliding reverse propagation transverse wave fracture imaging map.

[0102] In an optional embodiment, the suppression module 13 is specifically used for:

[0103] Based on the shear wave time difference curve and the shear wave arrival time curve, the shear wave direct wave suppression is performed on the gliding shear wave array waveform data.

[0104] Based on the shear wave time difference curve, continuous oscillating wave suppression is applied to the gliding shear wave array waveform data after shear wave direct wave suppression.

[0105] Based on the shear wave time difference curve, scattered noise suppression is performed on the gliding shear wave array waveform data after continuous oscillation wave suppression.

[0106] Based on the shear wave time difference curve, incoherent noise suppression is performed on the scatter noise-suppressed shear wave array waveform data to obtain the sliding reverse propagation shear wave signal of the target well.

[0107] Among them, the methods for suppressing shear wave direct waves include shear wave transformation and Radon domain filtering transformation;

[0108] Methods for suppressing continuous oscillating waves include: two-dimensional filtering, wavelet transform frequency division denoising, or SRME multiple wave attenuation;

[0109] Methods for suppressing scattered noise include: inverse Q filtering and linear amplitude enhancement;

[0110] Methods for suppressing incoherent noise include: superposition of common center points and superposition of common reflection points.

[0111] In another alternative embodiment, the preprocessing module 11 is specifically used for:

[0112] Gain recovery is performed on the dipole array waveform data of the target well based on the dipole array acoustic gain curve of the target well.

[0113] For the dipole array waveform data after gain recovery, delay recovery is performed based on the dipole array acoustic delay curve of the target well so that the dipole array waveform data is padded with zeros before time zero.

[0114] Bandpass filtering is performed on the delayed-recovered dipole array waveform data to obtain the sliding shear wave array waveform data of the target well.

[0115] In another optional embodiment, the extraction module 12 is specifically used to: extract the shear wave time difference curve and the shear wave arrival time curve from the gliding shear wave array waveform data based on the correlation coefficient method.

[0116] In another alternative embodiment, refer to Figure 8As shown, the device may further include: an acquisition module 10, which is used to acquire conventional logging data and array acoustic logging data of the target well, so as to use the conventional logging data as an indicator; wherein, the conventional logging data may include: natural gamma logging curves, caliper logging curves, or density logging curves; the array acoustic logging data may include: monopole array acoustic raw waveform curves, monopole array acoustic gain curves, monopole array acoustic delay curves, and dipole array acoustic raw waveform curves, dipole array acoustic gain curves, and dipole array acoustic delay curves.

[0117] Based on the same inventive concept, this embodiment of the invention also provides a computer-readable storage medium storing a computer program thereon, which, when executed by a processor, implements the above-described method for identifying well-side fractures based on dipole array waveform data.

[0118] Based on the same inventive concept, this embodiment of the invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, it implements the above-described method for identifying wellside fractures based on dipole array waveform data.

[0119] The principles underlying the problems solved by the above-mentioned devices, media, related equipment, and oil and gas well fracturing process construction methods in the embodiments of the present invention are similar to those of the aforementioned methods. Therefore, their implementation can refer to the implementation of the aforementioned methods, and repeated details will not be repeated.

[0120] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage and optical storage) containing computer-usable program code.

[0121] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0122] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0123] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0124] Obviously, those skilled in the art can make various modifications and variations to this invention without departing from its spirit and scope. Therefore, if these modifications and variations fall within the scope of the claims of this invention and their equivalents, this invention also intends to include these modifications and variations.

Claims

1. A method for identifying well-side fractures based on dipole array waveform data, characterized in that, include: The dipole array waveform data of the target well is preprocessed to obtain the sliding shear wave array waveform data of the target well; Extract the shear wave time difference curve and the shear wave arrival time curve from the waveform data of the gliding shear wave array; Based on the shear wave time difference curve and the shear wave arrival time curve, the shear wave array waveform data is subjected to shear wave direct wave suppression; based on the shear wave time difference curve, the shear wave array waveform data after shear wave direct wave suppression is subjected to continuous oscillation wave suppression; based on the shear wave time difference curve, the shear wave array waveform data after continuous oscillation wave suppression is subjected to scattered noise suppression; based on the shear wave time difference curve, the shear wave array waveform data after scattered noise suppression is subjected to incoherent noise suppression, so as to obtain the shear wave signal of the target well. Based on the sliding reverse propagation shear wave signal, separation imaging is performed in the frequency-wavenumber domain to obtain the sliding reverse propagation shear wave fracture imaging map of the target well. Based on the sliding reverse transverse wave fracture imaging map, the well-side fractures of the target well are identified.

2. The method according to claim 1, characterized in that, Methods for suppressing shear waves include shear wave transformation and Radon domain filtering transformation; Methods for suppressing continuous oscillating waves include: two-dimensional filtering, wavelet transform frequency division denoising, or SRME multiple wave attenuation; Methods for suppressing scattered noise include: inverse Q filtering and linear amplitude enhancement; Methods for suppressing incoherent noise include: superposition of common center points and superposition of common reflection points.

3. The method according to claim 1, characterized in that, The preprocessing of the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well includes: Gain recovery is performed on the dipole array waveform data of the target well based on the dipole array acoustic gain curve of the target well; For the dipole array waveform data after gain recovery, delay recovery is performed based on the dipole array acoustic delay curve of the target well so that the dipole array waveform data is padded with zeros before time zero. Bandpass filtering is performed on the delayed-recovered dipole array waveform data to obtain the sliding shear wave array waveform data of the target well.

4. The method according to claim 1, characterized in that, The extraction of the shear wave time difference curve and shear wave arrival time curve from the gliding shear wave array waveform data includes: Based on the correlation coefficient method, the shear wave time difference curve and shear wave arrival time curve are extracted from the shear wave array waveform data.

5. The method according to any one of claims 1 to 4, characterized in that, Also includes: Obtain conventional logging data and array sonic logging data of the target well, and use the conventional logging data as an indicator. The conventional logging data includes: natural gamma logging curves, caliper logging curves, or density logging curves; The array acoustic logging data includes: the original waveform curve of monopole array acoustic wave, the gain curve of monopole array acoustic wave, the delay curve of monopole array acoustic wave, as well as the original waveform curve of dipole array acoustic wave, the gain curve of dipole array acoustic wave, and the delay curve of dipole array acoustic wave.

6. A method for constructing fracturing technology in oil and gas wells, characterized in that, include: Design perforation locations and / or fracturing process schemes based on the wellside fractures of the target well; The wellbore fractures of the target well are identified by the method for identifying wellbore fractures based on dipole array waveform data according to any one of claims 1 to 5.

7. A device for identifying well-side fractures based on dipole array waveform data, characterized in that, include: The preprocessing module is used to preprocess the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well. The extraction module is used to extract the shear wave time difference curve and the shear wave arrival time curve from the shear wave array waveform data; The suppression module is used to suppress the shear wave array waveform data based on the shear wave time difference curve and the shear wave arrival time curve, to suppress the shear wave direct wave, to suppress the shear wave array waveform data after shear wave direct wave suppression, to suppress the continuous oscillation wave, to suppress the shear wave array waveform data after continuous oscillation wave suppression, and to suppress the incoherent noise, to obtain the shear wave signal of the target well. An imaging module is used to perform separation imaging in the frequency-wavenumber domain based on the sliding back-propagating shear wave signal to obtain a sliding back-propagating shear wave fracture imaging map of the target well. The identification module is used to identify the well-side fractures of the target well based on the sliding reverse transverse wave fracture imaging map.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by the processor, the program implements the method for identifying well-side fractures based on dipole array waveform data as described in any one of claims 1 to 5.

9. A computer device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method for identifying well-side fractures based on dipole array waveform data as described in any one of claims 1 to 5.

Images