Method, device and equipment for identifying near-well crack based on dipole array waveform data

Abstract

The invention discloses a method, device and equipment for identifying a near-well crack based on dipole array waveform data, and the method comprises the steps: carrying out the preprocessing of the dipole array waveform data of a target well, so as to obtain the 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; based on the transverse wave time difference curve and / or the transverse wave arrival time curve, interference signal suppression is carried out on the sliding transverse wave array waveform data to obtain a sliding back propagation transverse wave signal of the target well; performing separation imaging in a frequency-wave number domain based on the gliding back propagation shear wave signal to obtain a gliding back propagation shear wave crack imaging picture of the target well; and identifying the well-side crack of the target well based on the sliding back propagation transverse wave crack imaging graph. According to the method, the resolution ratio and the recognition precision of crack recognition can be improved, and technical support is provided for subsequent perforation and fracturing of a target well.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas reservoir exploration, and in particular to a method, a device and equipment for identifying wellside cracks based on dipole array waveform data. Background Art

[0002] Carbonate reservoirs are the main replacement area for my country's onshore oil reserve resources and the main target for reserve growth. The reservoir rock mass is highly heterogeneous, the matrix porosity is low, and fractures and caves are widely developed. In the reservoir, fractures are not only channels for fluid storage and migration, but also control the development of dissolved pores and caves, the distribution of original fluids in the formation, and the characteristics of mud invasion, etc., which greatly contribute to the permeability of the reservoir. Therefore, accurate identification and quantitative evaluation of the development of wellbore fractures through array acoustic logging data are of great significance to the exploration and development of carbonate oil and gas.

[0003] Morris et al. (1964) first studied the relationship between the amplitude of longitudinal waves and the width of fractures in acoustic logging data, and found that the attenuation of longitudinal waves when passing through horizontal fractures was very small and could be basically ignored. In addition, they also found that the transmission coefficient of shear waves was affected by the fracture inclination, and the shear wave attenuation was more obvious at horizontal fractures. Paillet (1980) and Zlatev et al. (1988) studied the attenuation of longitudinal and shear waves at fractures through 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 in time and fourth order in space. The study showed that when the dipole sound source is parallel to the formation interface, the formation interface has little effect on the waveform. When the dipole is perpendicular to the formation interface, there is strong dipole shear wave interference. Zhu et al. (1994) used aluminum and fluorite borehole models to study the propagation of dipole shear waves in fractures and found that dipole shear waves are more sensitive to horizontal fractures than vertical fractures. Haldorsen et al. (2006) studied the propagation of P- and S-waves in irregular boreholes and anisotropic formations respectively.

[0004] The above research content mainly focuses on how to use borehole mode waves, such as the velocity and attenuation properties of direct longitudinal waves and direct transverse waves, to distinguish wellbore cracks. Summary of the invention

[0005] In order to improve the recognition accuracy of wellbore fractures, thereby enriching process routes and increasing selection space, a method, device and apparatus for identifying wellbore fractures based on dipole array waveform data are provided in an embodiment of the present invention.

[0006] In a first aspect, an embodiment of the present invention provides a method for identifying wellbore fractures based on dipole array waveform data, which may include:

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

[0008] Extracting a shear wave time difference curve and a shear wave arrival time curve from the sliding shear wave array waveform data;

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

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

[0011] The wellbore fractures of the target well are identified based on the sliding back-propagation shear wave fracture imaging image.

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

[0013] Based on the shear wave time difference curve and the shear wave arrival time curve, performing shear wave direct wave suppression on the sliding shear wave array waveform data;

[0014] Based on the shear wave time difference curve, continuous shock wave suppression is performed on the sliding shear wave array waveform data after the shear wave direct wave suppression;

[0015] Based on the shear wave time difference curve, the sliding shear wave array waveform data after continuous shock wave suppression is subjected to scattered noise suppression;

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

[0017] Optionally, the method of suppressing the shear wave direct wave may include shear wave variation, Radon domain filtering variation;

[0018] Methods for suppressing continuous shock 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, linear amplitude enhancement;

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

[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] Performing gain recovery on the dipole array waveform data of the target well based on the dipole array acoustic wave gain curve of the target well;

[0023] Performing delay recovery on the dipole array waveform data after gain recovery based on the dipole array acoustic wave delay curve of the target well, so as to fill zero before the zero time of the dipole array waveform data;

[0024] The dipole array waveform data after delayed recovery is subjected to bandpass filtering to obtain the sliding shear wave array waveform data of the target well.

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

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

[0027] Optionally, the method may further include:

[0028] Acquiring conventional logging data and array acoustic logging data of the target well, and using the conventional logging data as an indicator channel;

[0029] Wherein, the conventional logging data may include: natural gamma logging curve, caliper logging curve or density logging curve;

[0030] The array acoustic logging data may include: a monopole array acoustic wave original waveform curve, a monopole array acoustic wave gain curve, a monopole array acoustic wave delay curve, and a dipole array acoustic wave original waveform curve, a dipole array acoustic wave gain curve, and a dipole array acoustic wave delay curve.

[0031] In a second aspect, an embodiment of the present invention provides an oil and gas well fracturing process construction method, which may include: designing a perforation position and / or a fracturing process scheme based on wellside fractures of a target well;

[0032] Wherein, 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] In a third aspect, an embodiment of the present invention provides a device for identifying wellbore fractures based on dipole array waveform data, which may include:

[0034] A preprocessing module, used for preprocessing the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well;

[0035] An extraction module, used for extracting a shear wave time difference curve and a shear wave arrival time curve from the sliding shear wave array waveform data;

[0036] A suppression module, used for suppressing interference signals on the sliding 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 sliding reverse propagation shear wave signal of the target well;

[0037] An imaging module, used for performing separation imaging in the frequency-wavenumber domain based on the sliding back-propagation shear wave signal, so as to obtain a sliding back-propagation shear wave fracture imaging map of the target well;

[0038] An identification module is used to identify the wellbore fractures of the target well based on the sliding back-propagation shear wave fracture imaging image.

[0039] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method for identifying wellbore fractures based on dipole array waveform data as described in the first aspect.

[0040] In a fourth aspect, an embodiment of the present invention provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method for identifying wellbore fractures based on dipole array waveform data as described in the first aspect is implemented.

[0041] The beneficial effects of the above technical solution provided by the embodiment of the present invention include at least:

[0042] The embodiment of the present invention provides a method, device and apparatus for identifying wellside fractures based on dipole array waveform data. The overall inventive concept of the method is to perform sliding back propagation imaging on the sliding shear wave array waveform data in the dipole array waveform data, and then perform wellside fracture identification on the sliding back propagation shear wave fracture imaging image. The method can improve the resolution and recognition accuracy of fracture identification, and provide technical support for subsequent perforation and fracturing of target wells.

[0043] Other features and advantages of the present invention will be described in the following description, and partly become apparent from the description, or understood by practicing the present invention. The purpose and other advantages of the present invention can be realized and obtained by the structures particularly pointed out in the written description and the accompanying drawings.

[0044] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the present invention and do not constitute a limitation of the present invention. In the accompanying drawings:

[0046] Figure 1 A flow chart of a method for identifying wellbore fractures based on dipole array waveform data provided in an embodiment of the present invention;

[0047] Figure 2 A flowchart of a method for identifying wellbore fractures based on dipole array waveform data in detail provided in an embodiment of the present invention;

[0048] Figure 3 An example of interference signal suppression processing for a sliding shear wave waveform provided in an embodiment of the present invention;

[0049] Figure 4 It is one of the counter-examples of implementing interference signal suppression processing for the sliding shear wave waveform provided in the embodiment of the present invention;

[0050] Figure 5 The second counter-example of interference signal suppression processing implemented for the sliding shear wave waveform provided in the embodiment of the present invention;

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

[0052] Figure 7 An example of identifying wellbore fractures using sliding back-propagation shear wave fracture imaging provided in an embodiment of the present invention;

[0053] Figure 8 It is a schematic diagram of the structure of a device for identifying wellbore fractures based on dipole array waveform data provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0054] The exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present disclosure are shown in the accompanying drawings, it should be understood that the present disclosure can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided to enable a more thorough understanding of the present disclosure and to fully convey the scope of the present disclosure to those skilled in the art.

[0055] The inventors found that the existing research methods essentially use the "concave and convex" changes of the curve to identify fractures, and their resolution and recognition accuracy are both low. In view of the above problems, the present invention is proposed to provide a method, device and equipment for identifying wellside fractures based on dipole array waveform data to overcome the above problems or at least partially solve the above problems.

[0056] The embodiment of the present invention provides a method for identifying wellbore fractures based on dipole array waveform data, referring to Figure 1 As shown, the method may include the following steps:

[0057] Step S11, pre-processing 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, extracting the shear wave time difference curve and the shear wave arrival time curve from the sliding shear wave array waveform data.

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

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

[0061] Step S15: identifying the wellbore fractures of the target well based on the sliding back-propagation shear wave fracture imaging image.

[0062] The above method provided in the embodiment of the present invention has the overall inventive concept of performing sliding back propagation imaging on the sliding shear wave array waveform data in the dipole array waveform data, and then performing wellside fracture identification on the sliding back propagation shear wave fracture imaging image based on the sliding back propagation. The method can improve the resolution and recognition accuracy of fracture identification, and provide technical support for subsequent perforation and fracturing of target wells.

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

[0064] Step S20, obtaining conventional logging data and array acoustic logging data of the target well, using conventional logging data as an indicator; wherein conventional logging data may include: natural gamma logging curve (GR), caliper logging curve or density logging curve; array acoustic logging data may include: monopole array acoustic wave original waveform curve, monopole array acoustic wave gain curve, monopole array acoustic wave delay curve, and dipole array acoustic wave original waveform curve, dipole array acoustic wave gain curve, dipole array acoustic wave delay curve. Conventional logging curves in this embodiment, such as natural gamma logging curves, can be used to verify the results later, so as to verify the accuracy of the above method provided in the embodiment of the present invention.

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

[0066] Step S21 , preprocessing 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 process of preprocessing is as follows: first, the dipole array waveform data of the target well is gain restored based on the dipole array acoustic wave gain curve of the target well; then, the dipole array waveform data after gain recovery is delayed restored based on the dipole array acoustic wave delay curve of the target well, so that the dipole array waveform data is filled with zero before the zero moment; finally, the dipole array waveform data after delay recovery is bandpass filtered to obtain the sliding shear wave array waveform data of the target well.

[0068] In the specific implementation, the curve preprocessing of array acoustic logging data mainly includes gain recovery, delay recovery and bandpass filtering. First, gain recovery is performed, then delay recovery is performed, and finally bandpass filtering is performed. When collecting array acoustic waveforms underground, in order to ensure that the amplitude of the waveform is at the highest accuracy, automatic gain control is used during collection, so that the waveform amplitude is amplified from a decimal to an integer for easy storage and recording. Therefore, for the collected waveform, gain recovery must be performed first, and the gain parameter AGN is shown in formula (1) as follows:

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

[0070] When gain recovery is performed, the array acoustic waveform is divided by AGN. In the waveform acquisition process, in order to reduce the storage volume of recorded data, the waveform data for a period of time before the arrival of the first wave is often not collected. The purpose of delay recovery is to fill zero before the zero moment of the array acoustic waveform, so as to obtain waveform data at the accurate time. The purpose of bandpass filtering is to obtain the sliding shear wave data, and the frequency band range of the bandpass filter is selected to be 3000~5000Hz.

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

[0072] Step S22, extracting the shear wave time difference curve and the shear wave arrival time curve from the sliding shear wave array waveform data.

[0073] This step is based on the correlation coefficient method to extract the shear wave time difference curve and shear wave arrival time curve from the sliding shear wave array waveform data. Figure 6 As shown, Figure 6The first channel is the gamma channel, the second channel is the depth channel, the third channel is the sliding shear wave waveform, and the fourth channel is the shear wave time difference time-time difference correlation diagram obtained by the correlation method. The horizontal axis in the correlation diagram is the time difference value, with 0 on the left and 300μs / ft on the right. The lighter the color in the diagram, the lower the correlation coefficient, and the darker the color, the higher the correlation coefficient. The dark strip shown in the figure is the location of the shear wave time difference. The fifth channel is the shear wave time difference curve, and the sixth channel is the shear wave arrival time curve.

[0074] In the specific implementation, this step calculates the correlation function of the array waveform in two dimensions, time and time difference, and the time difference where the maximum function value is located is the shear wave time difference. Specifically, first open a time window for the multi-channel waveform, calculate the correlation function value of the multi-channel waveform data in the window, then gradually change the slope of the window within the slowness search range (DTMIN-DTMAX), and then calculate the correlation function value. After the entire time difference range is calculated, the time where the window is located is moved forward by one position, and then the same time difference search calculation is performed until the entire time range (CWBEGIN-CWEND) is calculated. The calculation formula of the correlation function Corr is shown in formula (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 time length of the time window, N represents the number of waveforms in the array waveform, and RRSP represents the distance between two adjacent receivers. Based on the shear wave time difference, the shear wave arrival time curve can be directly calculated based on the well diameter and the length of the logging tool.

[0077] Step S23: 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 sliding shear wave array waveform data.

[0078] Among them, the method of suppressing the shear wave direct wave may include shear wave variation and Radon domain filtering variation;

[0079] Step S24: Based on the shear wave time difference curve, continuous shock wave suppression is performed on the sliding shear wave array waveform data after the shear wave direct wave suppression.

[0080] Among them, the method of continuous shock wave suppression may include: two-dimensional filtering method, wavelet transform frequency division denoising method or SRME multiple wave attenuation.

[0081] Step S25: Based on the shear wave time difference curve, the gliding shear wave array waveform data after the continuous shock wave is suppressed is subjected to scattered noise suppression.

[0082] Among them, the method of suppressing scattered noise may include: inverse Q filtering and linear amplitude enhancement.

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

[0084] The above steps S23 to S26 are for suppressing interference signals on the sliding shear wave array waveform data to obtain the sliding back-transmitted shear wave signal of the target well. This method is actually a set of step-by-step suppression processing methods for various types of interference signals. The four interference signal suppression processing technologies of sliding shear wave direct wave suppression, continuous shock wave suppression, scattered noise suppression, and irrelevant noise suppression are installed and processed in sequence to achieve the purpose of interference signal suppression. Among them, the direct wave suppression of sliding shear waves can be achieved by signal processing methods such as shear wave transform and Radon domain filtering; continuous shock wave suppression can be achieved by signal processing methods such as two-dimensional filtering, wavelet transform frequency division denoising or SRME multiple wave attenuation; scattered noise suppression can be achieved by signal processing methods such as inverse Q filtering and linear amplitude enhancement; non-correlated noise suppression is usually achieved by signal processing methods such as common center point superposition and common reflection point superposition.

[0085] Reference Figure 3 As shown in the figure, the fourth track shows the sliding shear wave waveform (WVXX); the fifth track is the waveform after the sliding shear wave direct signal is filtered out by shear wave transform (SLXX), at this time, the sliding shear wave back transmission is highlighted; the sixth track is the waveform after the far-end amplitude of the sliding shear wave back transmission is compensated by the amplitude linear enhancement method (GNXX), from which it can be observed that the far-end waveform of the sliding shear wave back transmission, 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 the SRME multiple waves are used to attenuate the multiple oscillation signals of the sliding shear wave back transmission (PDXX). Compared with the sixth track, the oscillation period of the sliding shear wave back transmission in this track is reduced, which is manifested as the thinning of the waveform "strip" in the image; the eighth track is the waveform after the waveforms are superimposed by multiple receivers (CXX), at this time, the non-correlated noise is further suppressed, and the sliding shear wave back transmission signal is more clearly displayed.

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

[0087] In this step, the sliding back-propagation shear wave fracture imaging mainly uses the frequency-wave number domain separation method to process the sliding back-propagation shear wave, so as to achieve the purpose of separating the sliding back-propagation shear wave upward wave and downward wave. Finally, the sliding back-propagation shear wave fracture imaging image is obtained, with the sliding back-propagation shear wave upward reflection wave on the left, the wellbore in the middle, and the sliding back-propagation shear wave downward reflection wave on the right. Figure 3As shown in Figure 1, the ninth track is the imaging result of the sliding back-transmitted shear wave crack. This is an image obtained after the up-and-downward reflected wave separation process based on the sliding back-transmitted shear wave waveform in the eighth track. The left side of the image shows the up-and-downward wave of the sliding back-transmitted shear wave, and the right side shows the down-and-downward wave of the sliding back-transmitted shear wave.

[0088] Step S28: Identify the wellbore fractures of the target well based on the sliding back-propagation shear wave fracture imaging image.

[0089] This step is in the sliding back propagation shear wave crack imaging diagram, refer to Figure 3 and Figure 7 As shown in the figure, the inclined phase axis from the upper left to the lower right can be observed and marked with a straight line. The intersection of the straight line and the well axis is the position where the fracture passes through the wellbore, and the amplitude of the inclined phase axis represents the fracture parameters such as the fracture opening.

[0090] Also refer to Figure 7 As shown, Figure 3 The first track is the gamma curve track, the second track is the depth track, and the third track is the sliding back-transmitted shear wave through-well fracture imaging result. In the sliding back-transmitted shear wave through-well fracture imaging diagram, the left side is the sliding back-transmitted shear wave up-going wave, the middle vertical line represents the wellbore, and the right side is the sliding back-transmitted shear wave down-going wave. Multiple groups of inclined phase axes can be observed, all of which represent through-well fractures, and the depths of intersection with the wellbore are 5215m, 5233m, and 5243m respectively. The above three depths represent the depth of the fracture through the wellbore. The fracture depth position usually represents the reservoir layer, which will provide a reference for the subsequent perforation and fracturing scheme design of the well.

[0091] It should also be noted that the stepwise suppression processing method for various types of interference signals in steps S23 to S26 is the focus of interference signal processing in the embodiment of the present invention. The suppression of various types of interference signals and stepwise processing are two key points. In order to illustrate the advantages of this method, examples are given to illustrate the role of the above two key points.

[0092] Reference Figure 4 As shown in the figure, the scattered noise suppression processing step is missing, and the waveform SLXX (fifth track) is directly processed in the subsequent steps. In the final sliding back-transmission shear wave fracture imaging image IMAGEXX (eighth track), there is indeed a large amount of unfiltered scattered noise around the wellbore, which seriously affects the clear appearance of the effective back-transmission wave signal.

[0093] Reference Figure 5As shown in the figure, the interference signal suppression order does not follow the key point of the cascade processing of this method, but moves the suppression of multiple oscillation interference signals in the subsequent processing to the first step of suppression. Because the amplitude of the direct signal of the sliding shear wave is much higher than the amplitude of the sliding back-transmission wave at this time, the main purpose of this step becomes to suppress the multiple oscillation waves of the direct signal, and in this process, the sliding back-transmission wave is also suppressed as noise. This makes it difficult to identify the effective sliding shear wave back-transmission wave signal in the final sliding back-transmission shear wave crack imaging image IMAGEXX (track 9).

[0094] The method for identifying wellbore fractures based on dipole array waveform data provided in the embodiment of the present invention can effectively suppress multiple types of interference signals, thereby making the sliding shear wave back propagation signal more clearly visible, thereby accurately identifying wellbore fractures.

[0095] Based on the same inventive concept, an embodiment of the present invention also provides an oil and gas well fracturing process construction method, which may include: designing perforation positions and / or fracturing process plans based on the wellbore cracks of the target well; wherein the wellbore cracks of the target well are identified according to the above-mentioned method for identifying wellbore cracks based on dipole array waveform data.

[0096] Based on the same inventive concept, an embodiment of the present invention also provides a device for identifying wellbore 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] The extraction module 12 is used to extract the shear wave time difference curve and the shear wave arrival time curve from the sliding shear wave array waveform data;

[0099] The suppression module 13 is used to suppress interference signals on the sliding 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 sliding reverse transmission shear wave signal of the target well;

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

[0101] The identification module 15 is used to identify the wellbore fractures of the target well based on the sliding back-propagation shear wave fracture imaging image.

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

[0103] Based on the shear wave time difference curve and the shear wave arrival time curve, performing shear wave direct wave suppression on the sliding shear wave array waveform data;

[0104] Based on the shear wave time difference curve, continuous shock wave suppression is performed on the sliding shear wave array waveform data after the shear wave direct wave suppression;

[0105] Based on the shear wave time difference curve, the sliding shear wave array waveform data after continuous shock wave suppression is subjected to scattered noise suppression;

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

[0107] Among them, the methods of shear wave direct wave suppression include shear wave variation and Radon domain filter variation;

[0108] Methods for suppressing continuous shock 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, linear amplitude enhancement;

[0110] The methods of suppressing incoherent noise include: common center point superposition and common reflection point superposition.

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

[0112] Performing gain recovery on the dipole array waveform data of the target well based on the dipole array acoustic wave gain curve of the target well;

[0113] Performing delay recovery on the dipole array waveform data after gain recovery based on the dipole array acoustic wave delay curve of the target well, so as to fill zero before the zero time of the dipole array waveform data;

[0114] The dipole array waveform data after delayed recovery is subjected to bandpass filtering 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 sliding shear wave array waveform data based on the correlation coefficient method.

[0116] In another alternative embodiment, referring to Figure 8As shown, the device may also include: an acquisition module 10, the acquisition module 10 is used to acquire conventional logging data and array acoustic logging data of the target well, and use the conventional logging data as an indicator channel; wherein the conventional logging data may include: a natural gamma logging curve, a caliper logging curve or a density logging curve; the array acoustic logging data may include: a monopole array acoustic wave original waveform curve, a monopole array acoustic wave gain curve, a monopole array acoustic wave delay curve, and a dipole array acoustic wave original waveform curve, a dipole array acoustic wave gain curve, and a dipole array acoustic wave delay curve.

[0117] Based on the same inventive concept, an embodiment of the present invention further provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the method for identifying wellbore fractures based on dipole array waveform data is implemented.

[0118] Based on the same inventive concept, a computer device is also provided in an embodiment of the present invention, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the method for identifying wellbore fractures based on dipole array waveform data is implemented.

[0119] The principles of the problems solved by the above-mentioned device, medium, related equipment and oil and gas well fracturing process construction method in the embodiment of the present invention are similar to those of the above-mentioned method, so their implementation can refer to the implementation of the above-mentioned method, and the repeated parts will not be repeated.

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

[0121] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0122] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0123] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0124] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include these modifications and variations.

Claims

1. A method for identifying wellbore fractures based on dipole array waveform data, characterized in that: include: Preprocessing the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well; Extracting a shear wave time difference curve and a shear wave arrival time curve from the sliding shear wave array waveform data; Based on the shear wave time difference curve and / or the shear wave arrival time curve, interference signal suppression is performed on the sliding shear wave array waveform data to obtain the sliding reverse propagation shear wave signal of the target well; Based on the sliding back-propagation shear wave signal, separation imaging is performed in the frequency-wavenumber domain to obtain a sliding back-propagation shear wave fracture imaging map of the target well; The wellbore fractures of the target well are identified based on the sliding back-propagation shear wave fracture imaging image.

2. The method according to claim 1, characterized in that Based on the shear wave time difference curve and / or the shear wave arrival time curve, interference signal suppression is performed on the sliding shear wave array waveform data to obtain the sliding reverse propagation shear wave signal of the target well, including: Based on the shear wave time difference curve and the shear wave arrival time curve, performing shear wave direct wave suppression on the sliding shear wave array waveform data; Based on the shear wave time difference curve, continuous shock wave suppression is performed on the sliding shear wave array waveform data after the shear wave direct wave suppression; Based on the shear wave time difference curve, the sliding shear wave array waveform data after continuous shock wave suppression is subjected to scattered noise suppression; Based on the shear wave time difference curve, incoherent noise suppression is performed on the sliding shear wave array waveform data after the scattered noise is suppressed, so as to obtain the sliding reverse propagation shear wave signal of the target well.

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

4. 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: Performing gain recovery on the dipole array waveform data of the target well based on the dipole array acoustic wave gain curve of the target well; Performing delay recovery on the dipole array waveform data after gain recovery based on the dipole array acoustic wave delay curve of the target well, so as to fill zero before the zero time of the dipole array waveform data; The dipole array waveform data after delayed recovery is subjected to bandpass filtering to obtain the sliding shear wave array waveform data of the target well.

5. The method according to claim 1, characterized in that The step of extracting a shear wave time difference curve and a shear wave arrival time curve from the sliding shear wave array waveform data comprises: Based on the correlation coefficient method, a shear wave time difference curve and a shear wave arrival time curve are extracted from the sliding shear wave array waveform data.

6. The method according to any one of claims 1 to 5, characterized in that: Also includes: Acquiring conventional logging data and array acoustic logging data of the target well, and using the conventional logging data as an indicator channel; Wherein, the conventional logging data includes: natural gamma logging curve, caliper logging curve or density logging curve; The array acoustic logging data include: monopole array acoustic original waveform curve, monopole array acoustic gain curve, monopole array acoustic delay curve, as well as dipole array acoustic original waveform curve, dipole array acoustic gain curve, dipole array acoustic delay curve.

7. A method for fracturing a well, characterized in that: include: Design perforation locations and / or fracturing process plans based on wellside fractures of the target well; Wherein, 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 6.

8. A device for identifying wellbore fractures based on dipole array waveform data, characterized in that: include: A preprocessing module, used for preprocessing the dipole array waveform data of the target well to obtain the sliding shear wave array waveform data of the target well; An extraction module, used for extracting a shear wave time difference curve and a shear wave arrival time curve from the sliding shear wave array waveform data; A suppression module, for suppressing interference signals on the sliding 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 sliding reverse propagation shear wave signal of the target well; An imaging module, used for performing separation imaging in the frequency-wavenumber domain based on the sliding back-propagation shear wave signal, so as to obtain a sliding back-propagation shear wave fracture imaging map of the target well; An identification module is used to identify the wellbore fractures of the target well based on the sliding back-propagation shear wave fracture imaging image.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the method for identifying wellbore fractures based on dipole array waveform data according to any one of claims 1 to 6 is implemented.

10. 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, the method for identifying wellbore fractures based on dipole array waveform data according to any one of claims 1 to 6 is implemented.

Images