A fracture prediction method based on azimuthal ray poisson impedance body ellipse fitting

Abstract

The application relates to a fracture prediction method based on azimuthal ray Poisson impedance body ellipse fitting, and the method is characterized by comprising the following steps: constructing anisotropic ray Poisson impedance curves of each well in a fracture reservoir to be measured according to anisotropic three parameters and normal weak degree curves and tangent weak degree curves on the well; determining anisotropic ray Poisson impedance medium frequency bodies and anisotropic ray Poisson impedance low frequency bodies of each well in the fracture reservoir to be measured, and forming anisotropic ray Poisson impedance bodies of each well in the fracture reservoir to be measured; performing anisotropy prediction by using ellipse fitting based on the anisotropic ray Poisson impedance bodies, and forming a buried hill fracture reservoir plane distribution prediction result according to a target layer top and bottom time window; and the application can be widely applied to the technical field of oil and gas field development.

Description

Technical Field

[0001] This invention relates to the field of oil and gas field development technology, and in particular to a fracture prediction method based on azimuth ray Poisson impedance volume ellipse fitting. Background Technology

[0002] Utilizing the anisotropic properties of pre-stack wide-azimuth seismic data, such as amplitude, AVO gradient, and fast and slow shear waves, to predict the planar distribution of fracture zones is currently the mainstream approach for predicting fractured reservoirs in buried hills.

[0003] Poisson impedance, a function of Poisson ratio and density, reflects variations in subsurface lithology and hydrocarbon-bearing properties. It is widely used in pre-stack inversion based on narrow-azimuth seismic data, often guiding reservoir lithology and fluid prediction. Narrow-azimuth seismic data refers to post-stack seismic data volumes obtained by stacking narrow-azimuth pre-stack seismic gathers. It lacks anisotropic information such as amplitude, AVO gradient, and fast / slow shear waves, thus limiting it to traditional post-stack fracture prediction. The predicted fracture scale is coarse, failing to meet the accuracy requirements for fracture reservoir prediction during oil and gas field development. Furthermore, there are currently no research reports on fracture prediction using azimuth ray Poisson impedance ellipse fitting based on wide-azimuth seismic data. Summary of the Invention

[0004] To address the aforementioned problems, the purpose of this invention is to provide a fracture prediction method based on azimuth ray Poisson impedance ellipse fitting that can meet the accuracy requirements for fracture reservoir prediction during the oil and gas field development stage.

[0005] To achieve the above objectives, the present invention adopts the following technical solution: In a first aspect, a crack prediction method based on azimuth ray Poisson impedance volume ellipse fitting is provided, comprising:

[0006] Based on the anisotropic three parameters and the normal and tangential weakness curves on the well, anisotropic ray Poisson impedance curves of each well in the fractured reservoir to be tested are constructed.

[0007] The anisotropic ray Poisson impedance mid-frequency volume and anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested are determined, and the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested is formed.

[0008] Based on anisotropic ray Poisson impedance, anisotropy prediction is performed using elliptic fitting, and the planar distribution prediction results of buried hill fracture reservoirs are formed according to the top and bottom time windows of the target layer.

[0009] Furthermore, before constructing the anisotropic ray Poisson impedance curves for each well in the fractured reservoir to be tested, the following steps are also included:

[0010] Using anisotropic rock physics modeling formulas, and based on the well logging curves of the fractured reservoir to be tested, the anisotropic three parameters, as well as the normal and tangential weakness curves, are constructed in the well of the fractured reservoir to be tested.

[0011] Furthermore, the construction of anisotropic ray Poisson impedance curves for each well in the fractured reservoir under test, based on the anisotropic three parameters, normal weakness curve, and tangential weakness curve, includes:

[0012] Based on the incident angle range of the pre-stack seismic gather of the fractured reservoir to be tested, the incident angle θ is set to the average value of the upper 1 / 3, middle 1 / 3 and lower 1 / 3 of the incident angle range, and the azimuth angle φ is determined at the same time.

[0013] Based on the anisotropic ray Poisson impedance calculation formula, and according to the anisotropic three parameters, normal weakness curve and tangential weakness curve on the well, as well as the determined incident angle θ and azimuth angle φ, several anisotropic ray Poisson impedance curves of each well in the fractured reservoir to be tested are obtained.

[0014] Furthermore, the formula for calculating the anisotropic ray Poisson impedance is as follows:

[0015]

[0016] Where PI(θ,φ) is the Poisson impedance of the anisotropic ray, θ is the incident angle, and φ is the azimuth angle; Δ N Normal weakness; Δ T α represents the tangential weakness; α, β, and ρ represent the P-wave velocity, S-wave velocity, and density, respectively; A(θ) is the P-wave velocity dimension factor, B(θ) is the S-wave velocity dimension factor, and C(θ) is the density dimension factor; D(θ,φ) is the normal weakness proportionality coefficient; E(θ,φ) is the tangential weakness proportionality coefficient; and EI0 is the elastic impedance parameter.

[0017] Further, the determination of the anisotropic ray Poisson impedance mid-frequency volume and the anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested, and the formation of the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested, includes:

[0018] Based on the pre-stack azimuth seismic gather data of the fractured reservoir to be tested, the large-angle incident seismic data is partially superimposed by azimuth angle to form azimuth large-angle partially superimposed seismic data, which serves as the anisotropic ray Poisson impedance mid-frequency body for each well in the fractured reservoir to be tested.

[0019] Determine the anisotropic ray Poisson impedance low-frequency curves of each well in the fractured reservoir to be tested, and establish the anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested.

[0020] Based on the anisotropic ray Poisson impedance low-frequency volume and anisotropic ray Poisson impedance mid-frequency volume of each well in the fractured reservoir to be tested, anisotropic ray Poisson impedance volumes of each well in the fractured reservoir to be tested are formed.

[0021] Furthermore, the anisotropic ray Poisson impedance volume for each well in the fractured reservoir to be tested, based on the low-frequency and mid-frequency anisotropic ray Poisson impedance volumes of each well, comprises:

[0022] The Fourier formula is used to transform the low-frequency and mid-frequency anisotropic ray Poisson impedance data of each well in the fractured reservoir to the frequency domain.

[0023] Using a linear weighting method, the low-frequency anisotropic ray Poisson impedance data of each well in the fractured reservoir to be tested is weighted and fused with the mid-frequency anisotropic ray Poisson impedance data in the frequency domain to form the anisotropic ray Poisson impedance data of each well in the fractured reservoir to be tested.

[0024] Furthermore, the anisotropic prediction based on anisotropic ray Poisson impedance volume, using elliptic fitting, and forming a predicted result for the planar distribution of buried hill fracture reservoirs according to the top and bottom time windows of the target layer, includes:

[0025] Each anisotropic ray Poisson impedance volume at any sampling point in space corresponds to a value. Using the above values, a 6-point conventional ellipse fitting is performed on each sampling point to obtain the direction of the major axis, the length of the major axis, and the length of the minor axis. The direction of the major axis represents the fracture development direction of the sampling point, and the ratio of the length of the major axis to the length of the minor axis represents the fracture development intensity at that point. Thus, the fracture development intensity volume and fracture development direction volume of the fracture reservoir to be tested are obtained.

[0026] Using the top and bottom interpretation layers of the target layer as the top and bottom time windows for interlayer attribute extraction, the interlayer attributes of the fracture development intensity volume and fracture development direction volume of the fracture reservoir to be tested are extracted, and then the planar distribution prediction results of the buried hill fracture reservoir are obtained.

[0027] Secondly, a crack prediction system based on azimuth ray Poisson impedance ellipse fitting is provided, including:

[0028] An anisotropic ray Poisson impedance curve construction module is used to construct anisotropic ray Poisson impedance curves for each well in the fractured reservoir under test based on the anisotropic three parameters on the well, as well as the normal weakness curve and the tangential weakness curve.

[0029] An anisotropic ray Poisson impedance volume forming module is used to determine the anisotropic ray Poisson impedance mid-frequency volume and anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested, and to form the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested.

[0030] The anisotropy prediction module is used to perform anisotropy prediction based on anisotropic ray Poisson impedance volume and elliptic fitting, and to generate the planar distribution prediction results of buried hill fracture reservoirs according to the top and bottom time windows of the target layer.

[0031] Thirdly, a processing device is provided, including computer program instructions, wherein when the computer program instructions are executed by the processing device, they are used to implement the steps corresponding to the above-mentioned crack prediction method based on azimuth ray Poisson impedance volume ellipse fitting.

[0032] Fourthly, a computer-readable storage medium is provided, wherein computer program instructions are stored on the computer-readable storage medium, wherein the computer program instructions, when executed by a processor, are used to implement the steps corresponding to the above-mentioned crack prediction method based on azimuth ray Poisson impedance volume ellipse fitting.

[0033] The present invention has the following advantages due to the adoption of the above technical solutions:

[0034] 1. The well logging Poisson impedance curve in this invention has a clear formula, and the predicted Poisson impedance volume in different orientations can be calibrated with the well logging Poisson impedance curve.

[0035] 2. In this invention, the Poisson impedance volume varies significantly in different orientations, and the crack prediction based on the elliptical fitting of the orientation Poisson impedance volume has high accuracy.

[0036] In summary, this invention can be widely applied in the field of oil and gas field development technology. Attached Figure Description

[0037] Various other advantages and benefits will become apparent to those skilled in the art upon reading the following detailed description of preferred embodiments. The accompanying drawings are for illustrative purposes only and are not intended to limit the invention. Throughout the drawings, the same reference numerals denote the same parts. In the drawings:

[0038] Figure 1 This is a schematic diagram of a method flow provided in an embodiment of the present invention;

[0039] Figure 2 This is a schematic diagram of the anisotropic ray Poisson impedance curve of the buried hill fracture reservoir section of Well A in the Bohai M gas field, provided by an embodiment of the present invention.

[0040] Figure 3 This is a schematic diagram of an anisotropic ray Poisson impedance body well-connected profile constructed in the Bohai M gas field according to an embodiment of the present invention.

[0041] Figure 4This is a schematic diagram of the planar distribution of fractured reservoirs in the buried hill weathering zone, predicted by anisotropic ray Poisson impedance ellipse fitting in an embodiment of the present invention. Detailed Implementation

[0042] Exemplary embodiments of the invention will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the invention are shown in the drawings, it should be understood that the invention can be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided to enable a more thorough understanding of the invention and to fully convey the scope of the invention to those skilled in the art.

[0043] It should be understood that the terminology used herein is for the purpose of describing particular exemplary embodiments only and is not intended to be limiting. Unless the context clearly indicates otherwise, the singular forms “a,” “an,” and “described” as used herein may also include the plural forms. The terms “comprising,” “including,” “containing,” and “having” are inclusive and therefore indicate the presence of the stated features, steps, operations, elements, and / or components, but do not exclude the presence or addition of one or more other features, steps, operations, elements, components, and / or combinations thereof. The method steps, processes, and operations described herein are not construed as requiring them to be performed in a particular order described or illustrated unless the order of performance is explicitly indicated. It should also be understood that additional or alternative steps may be used.

[0044] Although terms such as first, second, third, etc., may be used in this document to describe multiple elements, components, regions, layers, and / or segments, these elements, components, regions, layers, and / or segments should not be limited by these terms. These terms may be used only to distinguish one element, component, region, layer, or segment from another. Unless the context clearly indicates otherwise, terms such as "first," "second," and other numerical terms used herein do not imply order or sequence. Therefore, the first element, component, region, layer, or segment discussed below may be referred to as the second element, component, region, layer, or segment without departing from the teachings of the exemplary embodiments.

[0045] Since narrow-azimuth seismic data can only be used for traditional post-stack fracture prediction, the predicted fracture scale is relatively coarse, which is difficult to meet the accuracy requirements of fracture reservoir prediction in the oil and gas field development stage. The fracture prediction method based on azimuth ray Poisson impedance volume ellipse fitting provided in this embodiment of the invention is based on wide-azimuth seismic data and uses azimuth ray Poisson impedance volume ellipse fitting to predict fractures. The predicted Poisson impedance volumes in different azimuths can be calibrated with well logging Poisson impedance curves, and the predicted fractures have high accuracy.

[0046] Example 1

[0047] like Figure 1 As shown, this embodiment provides a crack prediction method based on azimuth ray Poisson impedance volume ellipse fitting, including the following steps:

[0048] 1) Using the anisotropic rock physics modeling formula, based on the logging curves of the fractured reservoir to be tested, construct the three anisotropic parameters, normal weakness curve, and tangential weakness curve of the well in the fractured reservoir to be tested.

[0049] Specifically, the logging curves include P-wave velocity, S-wave velocity, density, and mineral composition logging curves.

[0050] Specifically, the three anisotropy parameters include ε, δ, and γ, where ε is the longitudinal wave anisotropy coefficient, used to measure the intensity of longitudinal wave anisotropy; δ is the longitudinal wave variation coefficient, used to measure the rate of change of longitudinal wave anisotropy in the vertical direction; and γ is the transverse wave anisotropy coefficient, used to measure the intensity of transverse wave anisotropy.

[0051] 2) Based on the anisotropic ray Poisson impedance calculation formula, and according to the anisotropic three parameters and the normal and tangential weakness curves of the well, the anisotropic ray Poisson impedance curves of each well in the fractured reservoir to be tested are constructed, specifically as follows:

[0052] 2.1) Based on the incident angle range of the pre-stack seismic gather of the fractured reservoir to be tested, let the incident angle θ be equal to the average value of the upper 1 / 3 (representing large incident angle), middle 1 / 3 (representing medium incident angle), and lower 1 / 3 (representing small incident angle) of the incident angle range, respectively. At the same time, let the azimuth φ be equal to 30° (representing 0-60° azimuth), 90° (representing 60-120° azimuth), 150° (representing 120-180° azimuth), 210° (representing 180-240° azimuth), 270° (representing 240-300° azimuth) and 330° (representing 300-360° azimuth), respectively.

[0053] 2.2) Based on the anisotropic ray Poisson impedance calculation formula, according to the anisotropic three parameters, normal weakness curve and tangential weakness curve on the well, as well as the determined incident angle θ and azimuth angle φ, 18 anisotropic ray Poisson impedance curves of each well in the fractured reservoir to be tested are obtained.

[0054] Specifically, the formula for calculating the Poisson impedance of anisotropic rays is:

[0055]

[0056] Where PI(θ,φ) is the Poisson impedance of the anisotropic ray, θ is the incident angle, and φ is the azimuth angle; Δ N Normal weakness; Δ TLet θ be the tangential weakness; α, β, and ρ be the P-wave velocity, S-wave velocity, and density, respectively; A(θ) be the P-wave velocity dimension factor, B(θ) be the S-wave velocity dimension factor, and C(θ) be the density dimension factor; D(θ,φ) be the normal weakness proportionality coefficient; E(θ,φ) be the tangential weakness proportionality coefficient; and EI0 be the elastic impedance parameter, where:

[0057] A(θ) = sec 2 θ (2)

[0058] B(θ) = -8gsin 2 θ (3)

[0059] C(θ) = 1 - 4gsin 2 θ (4)

[0060]

[0061] E(θ,φ)=2gcos 2 φsin 2 θ(1-sin 2 φsin 2 θ)tan 2 θ(1-g)] (6)

[0062] Where g is a function of the three parameters of anisotropy.

[0063] Specifically, the anisotropic ray Poisson impedance curves include ray Poisson impedance curves for small-angle, medium-angle, and large-angle incident rays in azimuths of 0-60°, 60-120°, 120-180°, 180-240°, 240-300°, and 300-360°.

[0064] 3) Based on the pre-stack azimuth seismic gather data of the fractured reservoir to be tested, the large-angle incident seismic data are partially superimposed by azimuth angle to form azimuth-based large-angle partially superimposed seismic data, which serves as the anisotropic ray Poisson impedance mid-frequency body for each well in the fractured reservoir to be tested, specifically:

[0065] 3.1) The pre-stack azimuth seismic gather data of the fractured reservoir to be tested were divided into azimuth gathers of 6 sectors according to azimuth angles of 0-60°, 60-120°, 120-180°, 180-240°, 240-300° and 300-360°.

[0066] 3.2) Extract the pre-stack sub-azimuth seismic gather data for each azimuth gather in the range of 10-30° incident angle according to the incident angle, forming 6 sub-azimuth sub-incident angle gathers.

[0067] 3.3) The six azimuth and incident angle gathers are stacked to obtain six azimuth large angle partial stacked seismic data, which are used as the anisotropic ray Poisson impedance mid-frequency body of each well in the fractured reservoir to be tested.

[0068] 4) Determine the anisotropic ray Poisson impedance low-frequency curves of each well in the fractured reservoir to be tested, and establish the anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested using spatial three-dimensional interpolation. Specifically:

[0069] 4.1) Design a 0-10Hz low-pass filter to perform low-pass filtering on the anisotropic ray Poisson impedance curves of each well in the fractured reservoir to be tested, and obtain the low-frequency anisotropic ray Poisson impedance curves of each well.

[0070] 4.2) Using conventional interpolation methods for layer-constrained wells, the anisotropic ray Poisson impedance low-frequency curves of each well are interpolated to obtain the anisotropic ray Poisson impedance low-frequency data volume, which is the anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested.

[0071] 5) The anisotropic ray Poisson impedance low-frequency volume and the anisotropic ray Poisson impedance mid-frequency volume of each well in the fractured reservoir to be tested are weighted and fused in the frequency domain to form the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested, specifically:

[0072] 5.1) Using the Fourier formula, the anisotropic ray Poisson impedance low-frequency volume and anisotropic ray Poisson impedance mid-frequency volume of each well in the fractured reservoir to be tested are transformed to the frequency domain.

[0073] 5.2) Using a linear weighting method, the low-frequency and mid-frequency anisotropic ray Poisson impedance data of each well in the fractured reservoir under test are weighted and fused in the frequency domain to form the anisotropic ray Poisson impedance data of each well in the fractured reservoir under test:

[0074] 5.2.1) Average the anisotropic ray Poisson impedance low-frequency volumes of each well in the fractured reservoir under test at each azimuth angle. The data volume of this average volume is on the order of a.

[0075] 5.2.2) Average the anisotropic ray Poisson impedance mid-frequency volume of each well in the fractured reservoir under test at each azimuth angle. The data volume of this average volume is on the order of b.

[0076] 5.2.3) Calculate the fusion coefficient s:

[0077] s = a / b (7)

[0078] 5.2.4) In the frequency domain, the anisotropic ray Poisson impedance low-frequency volume and anisotropic ray Poisson impedance mid-frequency volume of each well in the fractured reservoir to be tested, in the 0-60°, 60-120°, 120-180°, 180-240°, 240-300° and 300-360° azimuths, are weighted and fused according to a weight of 1:s to form the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested, thereby ensuring the energy balance of the data volume in each azimuth before and after fusion.

[0079] 6) Based on anisotropic ray Poisson impedance, anisotropy prediction is performed using elliptic fitting, and the planar distribution prediction results of buried hill fracture reservoirs are formed according to the top and bottom time windows of the target layer, specifically:

[0080] 6.1) Any sampling point in space corresponds to a value in one of the six anisotropic ray Poisson impedance volumes in the 0-60° azimuth, 60-120° azimuth, 120-180° azimuth, 180-240° azimuth, 240-300° azimuth, and 300-360° azimuth. Using the above six values, a six-point conventional ellipse fitting is performed on each sampling point to obtain the major axis direction, major axis length, and minor axis length of the ellipse. The major axis direction represents the fracture development direction of the sampling point, and the ratio of the major axis length to the minor axis length represents the fracture development intensity of the point. Thus, the fracture development intensity volume and fracture development direction volume of the fracture reservoir to be tested are obtained.

[0081] 6.2) Using the top and bottom interpretation layers of the target layer as the top and bottom time windows for interlayer attribute extraction, the interlayer attributes (average values) of the fracture development intensity volume and fracture development direction volume of the fracture reservoir to be tested are extracted, and then the plane distribution prediction results of the buried hill fracture reservoir (fracture development intensity attribute, fracture development direction attribute) are obtained.

[0082] The following is a detailed explanation of how the fracture prediction method based on azimuth ray Poisson impedance volume ellipse fitting of the present invention is used to predict the Archean buried hill fracture reservoir in the Bohai M gas field, taking it as a specific example:

[0083] Using anisotropic rock physics modeling formulas, based on the P-wave velocity, S-wave velocity, density, and mineral composition logging curves of Well A in the Archean buried hill fracture reservoir of the M gas field, the anisotropic three parameters ε, δ, and γ of Well A, as well as the normal weakness Δ, are constructed. N Curve and tangential weakness Δ T curve.

[0084] 2) Based on the anisotropic ray Poisson impedance calculation formula, the anisotropic ray Poisson impedance curve of well A is constructed according to the three anisotropic parameters of well A and the normal and tangential weakness curves.

[0085] Specifically, 18 anisotropic ray Poisson impedance curves were constructed for Well A, representing small-angle, medium-angle, and large-angle incident ray Poisson impedance curves at azimuths of 0-60°, 60-120°, 120-180°, 180-240°, 240-300°, and 300-360°. Figure 2 As shown. The anisotropic ray Poisson impedance curves for the other seven wells in the M gas field were calculated using the same method.

[0086] 3) Based on the pre-stack azimuth seismic gather data of the M gas field, the seismic data stacks with an incident angle of 10-30° are selected to form 10-30° partial incident angle stacked seismic data volumes in azimuths of 0-60°, 60-120°, 120-180°, 180-240°, 240-300° and 300-360° (i.e., data in the 10-30° incident angle range are stacked in 6 azimuths to obtain 10-30° partial incident angle stacked seismic data volumes in 6 azimuths), and these 6 partial incident angle stacked seismic data volumes are used as anisotropic ray Poisson impedance mid-frequency volumes.

[0087] 4) Extract the anisotropic ray Poisson impedance curves of each well at large angles of incidence in the 0-60°, 60-120°, 120-180°, 180-240°, 240-300°, and 300-360° azimuths. Six curves can be obtained for each well. Perform 0-15Hz low-pass filtering on the above six curves of each well, and perform three-dimensional spatial interpolation on the six low-pass filtered curves to obtain the anisotropic ray Poisson impedance low-frequency volumes in the 0-60°, 60-120°, 120-180°, 180-240°, 240-300°, and 300-360° azimuths.

[0088] 5) The anisotropic ray Poisson impedance low-frequency and mid-frequency volumes in the above six orientations are transformed to the frequency domain using the Fourier transform formula, and then weighted and fused using a linear weighting method to form anisotropic ray Poisson impedance volumes in the six orientations, such as... Figure 3 As shown, this is a cross-section of anisotropic ray Poisson impedance volume in six directions.

[0089] 6) Based on anisotropic ray Poisson impedance volumes in six azimuths, anisotropy prediction is performed using ellipse fitting. Each sampling point has six azimuths: 0-60°, 60-120°, 120-180°, 180-240°, 240-300°, and 300-360°, resulting in six data samples. Ellipse fitting is performed on the six data samples at each sampling point. The ratio of the major axis to the minor axis of the fitted ellipse represents the anisotropy intensity at that sampling point, and the direction of the major axis represents the anisotropy direction at that sampling point. Based on the planar properties of the ellipse fitting results extracted from the top and bottom layers of the buried hill weathering zone in the M gas field, planar prediction results for the fractured reservoir in the buried hill weathering zone of the M gas field are formed, such as... Figure 4 As shown.

[0090] Example 2

[0091] This embodiment provides a crack prediction system based on azimuth ray Poisson impedance volume ellipse fitting, including:

[0092] The anisotropic ray Poisson impedance curve construction module is used to construct anisotropic ray Poisson impedance curves for each well in the fractured reservoir under test, based on the anisotropic three parameters on the well, as well as the normal and tangential weakness curves.

[0093] An anisotropic ray Poisson impedance volume forming module is used to determine the anisotropic ray Poisson impedance mid-frequency volume and anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested, and to form the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested.

[0094] The anisotropy prediction module is used to perform anisotropy prediction based on anisotropic ray Poisson impedance volume and elliptic fitting, and to generate the planar distribution prediction results of buried hill fracture reservoirs according to the top and bottom time windows of the target layer.

[0095] In a preferred embodiment, it further includes:

[0096] The anisotropic rock physics modeling module is used to construct the three anisotropic parameters, normal weakness curve, and tangential weakness curve of the well in the fractured reservoir to be tested, based on the well logging curve of the fractured reservoir to be tested, using anisotropic rock physics modeling formulas.

[0097] In a preferred embodiment, the anisotropic ray Poisson impedance forming module includes:

[0098] The intermediate frequency volume building unit is used to perform partial superposition of large-angle incident seismic data by azimuth angle based on the pre-stack azimuth seismic gather data of the fractured reservoir to be tested, forming partial superposition of large-angle incident seismic data by azimuth angle, which serves as the anisotropic ray Poisson impedance intermediate frequency volume of each well in the fractured reservoir to be tested.

[0099] The low-frequency volume establishment unit is used to extract the anisotropic ray Poisson impedance low-frequency curves of each well in the fractured reservoir to be tested, and to establish the anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested using a spatial three-dimensional interpolation method.

[0100] The weighted fusion unit is used to perform weighted fusion of the low-frequency anisotropic ray Poisson impedance data of each well in the fractured reservoir under test with the mid-frequency anisotropic ray Poisson impedance data in the frequency domain, so as to form the anisotropic ray Poisson impedance data of each well in the fractured reservoir under test.

[0101] The system provided in this embodiment is used to execute the above-described method embodiments. For specific processes and details, please refer to the above embodiments, which will not be repeated here.

[0102] Example 3

[0103] This embodiment provides a processing device corresponding to the crack prediction method based on azimuth ray Poisson impedance ellipse fitting provided in Embodiment 1. The processing device can be applied to client processing devices, such as mobile phones, laptops, tablets, desktop computers, etc., to execute the method of Embodiment 1.

[0104] The processing device includes a processor, a memory, a communication interface, and a bus. The processor, memory, and communication interface are connected via the bus to enable communication between them. The memory stores a computer program that can run on the processing device. When the processing device runs the computer program, it executes the crack prediction method based on azimuth ray Poisson impedance ellipse fitting provided in Embodiment 1.

[0105] In some implementations, the memory may be high-speed random access memory (RAM), and may also include non-volatile memory, such as at least one disk storage device.

[0106] In other implementations, the processor can be any type of general-purpose processor, such as a central processing unit (CPU) or a digital signal processor (DSP), and there is no limitation here.

[0107] Furthermore, the logical instructions in the aforementioned memory can be implemented as software functional units and sold or used as independent products, and can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0108] Those skilled in the art will understand that the structure of the above-described computing device is only a partial structure related to the solution of this application and does not constitute a limitation on the computing device to which the solution of this application is applied. A specific computing device may include more or fewer components, or combine certain components, or have different component arrangements.

[0109] Example 4

[0110] This embodiment provides a computer program product corresponding to the crack prediction method based on azimuth ray Poisson impedance ellipse fitting provided in Embodiment 1. The computer program product may include a computer-readable storage medium on which computer-readable program instructions for executing the crack prediction method based on azimuth ray Poisson impedance ellipse fitting described in Embodiment 1 are loaded.

[0111] A computer-readable storage medium can be a tangible device that holds and stores instructions for use by an instruction execution device. A computer-readable storage medium can be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any combination thereof.

[0112] The computer-readable storage medium provided in the above embodiments has a similar implementation principle and technical effect to the above method embodiments, and will not be described again here.

[0113] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. 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... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0114] 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.

[0115] 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.

[0116] The above embodiments are only used to illustrate the present invention. The structure, connection method and manufacturing process of each component can be varied. All equivalent transformations and improvements made on the basis of the technical solution of the present invention should not be excluded from the protection scope of the present invention.

Claims

1. A crack prediction method based on azimuth ray Poisson impedance ellipse fitting, characterized in that, include: Based on the anisotropic three parameters and the normal and tangential weakness curves on the well, anisotropic ray Poisson impedance curves of each well in the fractured reservoir to be tested are constructed. The anisotropic ray Poisson impedance mid-frequency volume and anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested are determined, and the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested is formed. Based on anisotropic ray Poisson impedance, anisotropy prediction is performed using elliptic fitting, and the planar distribution prediction results of buried hill fracture reservoirs are formed according to the top and bottom time windows of the target layer.

2. The crack prediction method based on azimuth ray Poisson impedance ellipse fitting as described in claim 1, characterized in that, Before constructing the anisotropic ray Poisson impedance curves for each well in the fractured reservoir to be tested, the following steps are also included: Using anisotropic rock physics modeling formulas, and based on the well logging curves of the fractured reservoir to be tested, the anisotropic three parameters, as well as the normal and tangential weakness curves, are constructed in the well of the fractured reservoir to be tested.

3. The crack prediction method based on azimuth ray Poisson impedance ellipse fitting as described in claim 1, characterized in that, The process involves constructing anisotropic ray Poisson impedance curves for each well in the fractured reservoir under test, based on the anisotropic three parameters, normal weakness curve, and tangential weakness curve. This includes: Based on the range of incident angles of the pre-stack seismic gathers of the fractured reservoir to be tested, let the incident angles be... It equals the average of the upper 1 / 3, middle 1 / 3, and lower 1 / 3 of the range of incident angle values, and simultaneously determines the azimuth angle. ; Based on the Poisson impedance calculation formula for anisotropic rays, and according to the three anisotropic parameters above the well, the normal weakness curve and the tangential weakness curve, and the determined incident angle... and azimuth This yielded several anisotropic ray Poisson impedance curves for each well in the fractured reservoir to be tested.

4. The crack prediction method based on azimuth ray Poisson impedance ellipse fitting as described in claim 3, characterized in that, The formula for calculating the anisotropic ray Poisson impedance is as follows: in, For anisotropic rays, Poisson impedance, Angle of incidence It is the azimuth angle; Normal weakness; Tangential weakness; , , These are longitudinal wave velocity, transverse wave velocity, and density, respectively. The dimension factor for P-wave velocity. The transverse wave velocity dimension factor, Density dimension factor; This is the normal weakness proportionality coefficient; This is the tangential weakness proportionality coefficient; This is the elastic impedance parameter.

5. The crack prediction method based on azimuth ray Poisson impedance ellipse fitting as described in claim 1, characterized in that, The process of determining the anisotropic ray Poisson impedance mid-frequency volume and anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested, and forming the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested, includes: Based on the pre-stack azimuth seismic gather data of the fractured reservoir to be tested, the large-angle incident seismic data is partially superimposed by azimuth angle to form azimuth large-angle partially superimposed seismic data, which serves as the anisotropic ray Poisson impedance mid-frequency body for each well in the fractured reservoir to be tested. Determine the anisotropic ray Poisson impedance low-frequency curves of each well in the fractured reservoir to be tested, and establish the anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested. Based on the anisotropic ray Poisson impedance low-frequency volume and anisotropic ray Poisson impedance mid-frequency volume of each well in the fractured reservoir to be tested, anisotropic ray Poisson impedance volumes of each well in the fractured reservoir to be tested are formed.

6. The crack prediction method based on azimuth ray Poisson impedance ellipse fitting as described in claim 5, characterized in that, The anisotropic ray Poisson impedance volume for each well in the fractured reservoir under test is formed based on the low-frequency and mid-frequency anisotropic ray Poisson impedance volumes of each well in the fractured reservoir under test, including: The Fourier formula is used to transform the low-frequency and mid-frequency anisotropic ray Poisson impedance data of each well in the fractured reservoir to the frequency domain. Using a linear weighting method, the low-frequency anisotropic ray Poisson impedance data of each well in the fractured reservoir to be tested is weighted and fused with the mid-frequency anisotropic ray Poisson impedance data in the frequency domain to form the anisotropic ray Poisson impedance data of each well in the fractured reservoir to be tested.

7. The crack prediction method based on azimuth ray Poisson impedance ellipse fitting as described in claim 1, characterized in that, The method, based on anisotropic ray Poisson impedance, uses elliptic fitting for anisotropic prediction and, according to the top and bottom time windows of the target layer, forms a predicted result for the planar distribution of buried hill fracture reservoirs, including: Each anisotropic ray Poisson impedance volume at any sampling point in space corresponds to a value. Using the above values, a 6-point conventional ellipse fitting is performed on each sampling point to obtain the direction of the major axis, the length of the major axis, and the length of the minor axis. The direction of the major axis represents the fracture development direction of the sampling point, and the ratio of the length of the major axis to the length of the minor axis represents the fracture development intensity at that point. Thus, the fracture development intensity volume and fracture development direction volume of the fracture reservoir to be tested are obtained. Using the top and bottom interpretation layers of the target layer as the top and bottom time windows for interlayer attribute extraction, the interlayer attributes of the fracture development intensity volume and fracture development direction volume of the fracture reservoir to be tested are extracted, and then the planar distribution prediction results of the buried hill fracture reservoir are obtained.

8. A crack prediction system based on azimuth ray Poisson impedance ellipse fitting, characterized in that, include: An anisotropic ray Poisson impedance curve construction module is used to construct anisotropic ray Poisson impedance curves for each well in the fractured reservoir under test based on the anisotropic three parameters on the well, as well as the normal weakness curve and the tangential weakness curve. An anisotropic ray Poisson impedance volume forming module is used to determine the anisotropic ray Poisson impedance mid-frequency volume and anisotropic ray Poisson impedance low-frequency volume of each well in the fractured reservoir to be tested, and to form the anisotropic ray Poisson impedance volume of each well in the fractured reservoir to be tested. The anisotropy prediction module is used to perform anisotropy prediction based on anisotropic ray Poisson impedance volume and elliptic fitting, and to generate the planar distribution prediction results of buried hill fracture reservoirs according to the top and bottom time windows of the target layer.

9. A processing device, characterized in that, It includes computer program instructions, wherein when executed by a processing device, the computer program instructions are used to implement the steps corresponding to the crack prediction method based on azimuth ray Poisson impedance volume ellipse fitting as described in any one of claims 1-7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer program instructions, wherein when executed by a processor, the computer program instructions are used to implement the steps corresponding to the crack prediction method based on azimuth ray Poisson impedance volume ellipse fitting as described in any one of claims 1-7.

Images