Method and device for weakening 90-degree ambiguity of fracture development orientation prediction, electronic equipment and storage medium

Abstract

The invention discloses a fracture development orientation prediction 90-degree ambiguity weakening method, electronic equipment, a storage medium and a fracture development orientation prediction 90-degree ambiguity weakening device. The method comprises the following steps: determining influence factors causing 90-degree fuzziness of fracture development orientation prediction based on an orientation elastic wave impedance equation; analyzing the influence factors based on a fracture orientation interpretation result of logging data or geological data and a fracture-induced orientation anisotropy model; and determining a strategy for weakening the 90-degree fuzziness of fracture development orientation prediction based on an analysis result. According to the method, starting from a root cause for generating fracture development orientation to predict 90-degree fuzziness, a mapping relation between a fracture orientation prediction result and fracture type reservoir geophysical parameters and combinations thereof is researched through model testing, and a strategy for weakening fracture development orientation to predict 90-degree fuzziness is determined through the mapping relation. Fracture orientation prediction with 90-degree ambiguity weakened is achieved, and reliability of fracture orientation seismic description is improved.

Description

Technical Field

[0001] This invention belongs to the field of geophysical exploration technology, and more specifically, relates to a method, electronic device, storage medium and apparatus for reducing 90-degree ambiguity in fracture development orientation prediction. Background Technology

[0002] To overcome the problem of inconsistent intercepts in all directions in conventional azimuth AVO gradient extraction methods, a stable azimuth AVO gradient unconstrained inversion method is proposed and has been successfully applied in actual work areas. To optimize traditional algorithms, previous researchers proposed an AVO inversion algorithm based on inverse operator estimation, which significantly improves the efficiency and accuracy of AVO inversion by directly inverting within a subdomain. AVAZ technology is an effective method based on seismic data analysis, mainly used for the prediction and monitoring of seismic fractures. It obtains AVAZ distribution maps in different directions by analyzing the velocity and amplitude variations of seismic waves, and then uses this data to predict the location and direction of seismic fractures. AVAZ technology has advantages such as high accuracy, high resolution, and non-invasiveness, and has been widely used in geological exploration, earthquake prediction, and earthquake engineering. While fracture prediction methods using AVAZ features are relatively simple, they require a high signal-to-noise ratio from seismic data, are significantly affected by noise, and suffer from ambiguity in predicting fracture azimuth. The amplitude of crack-induced reflection varies with orientation according to a cosine function. The crack orientation corresponds to the minimum or maximum point of the amplitude orientation variation, which is related to the medium parameters. Since the maximum and minimum points differ by 90 degrees, there is a 90-degree ambiguity in the AVAZ prediction of crack parameters.

[0003] The information disclosed in the background section of this invention is intended only to enhance the understanding of the general background of this invention, and should not be construed as an admission or in any way implying that such information constitutes prior art known to those skilled in the art. Summary of the Invention

[0004] The purpose of this invention is to propose a method, electronic device, storage medium, and apparatus for reducing 90-degree ambiguity in crack development azimuth prediction, thereby reducing 90-degree ambiguity in crack azimuth prediction and improving the reliability of crack azimuth seismic description.

[0005] To achieve the above objectives, the present invention proposes a method, electronic device, storage medium, and apparatus for reducing the 90-degree ambiguity in crack development orientation prediction.

[0006] According to a first aspect of the present invention, a method for weakening the 90-degree ambiguity in crack development orientation prediction is proposed, comprising:

[0007] Based on the azimuth elastic wave impedance equation, the influencing factors that lead to the 90-degree ambiguity in the prediction of crack development azimuth are determined.

[0008] The influencing factors are analyzed based on the fracture orientation interpretation results of well logging data or geological data and the fracture-induced orientation anisotropy model.

[0009] Based on the analysis results, a strategy was determined to reduce the 90-degree ambiguity in the prediction of crack development orientation.

[0010] Optionally, the calculation expression for the azimuth elastic wave impedance equation is as follows:

[0011]

[0012] in,

[0013]

[0014] EI0=V p0 ×ρ0;

[0015]

[0016]

[0017]

[0018] φ is the observation azimuth angle, θ is the incident angle of the seismic wave, φ sym Let be the crack dip azimuth angle, ρ be the density of the isotropic background medium, ρ0 be the average density of the isotropic background medium, and V be the density of the crack. p0 V is the average longitudinal wave velocity of the isotropic background medium. S0 V is the average transverse wave velocity of the isotropic background medium. p V is the longitudinal wave velocity of the isotropic background medium. S Let EI0 be the transverse wave velocity of the isotropic background medium, EI0 be the average wave impedance of the isotropic background medium, and δ be the wave velocity of the isotropic background medium. N Let Δδ be the crack normal weakness parameter of the isotropic background medium. N The difference in the normal weakness parameter of the cracks in two adjacent isotropic background media is represented by δ. V Let Δδ be the vertical weakness parameter of the crack in the isotropic background medium. V The difference in vertical weakness parameters of cracks in two adjacent isotropic background media, δ H The horizontal weakness parameter of the crack in the isotropic background medium. The coefficients are the second-order Fourier coefficients, φ well For interpretation results of fracture orientation in well logging or geological data, sign[·] represents the sign. For weakness parameters, Let g be the average of the squares of the transverse and longitudinal wave velocities of two adjacent isotropic background media, and g = (VS / V P ) 2 .

[0019] Optionally, the influencing factors are:

[0020] The uncertainty in the sign of the coefficients of the second-order term in the Fourier series expansion of elastic impedance, i.e. middle The uncertainty of the sign.

[0021] Optionally, the analysis results are:

[0022] and The symbols are consistent.

[0023] Optionally, the strategy includes:

[0024] pass extract The symbol is used to determine The symbol is used to reduce the 90-degree ambiguity in predicting the orientation of crack development.

[0025] According to a second aspect of the present invention, a device for weakening the 90-degree ambiguity in crack development orientation prediction is provided, comprising:

[0026] The first determination module is used to determine the influencing factors that cause the 90-degree ambiguity in the prediction of crack development azimuth based on the azimuth elastic wave impedance equation.

[0027] The analysis module is used to analyze the influencing factors based on the fracture orientation interpretation results of well logging data or geological data and the fracture-induced orientation anisotropy model.

[0028] The second determination module is used to determine a strategy to reduce the 90-degree ambiguity in the prediction of crack development orientation based on the analysis results.

[0029] Optionally, the calculation expression for the azimuth elastic wave impedance equation is as follows:

[0030]

[0031] in,

[0032]

[0033] EI0=V p0 ×ρ0;

[0034]

[0035]

[0036] φ is the observation azimuth angle, θ is the incident angle of the seismic wave, φ sym Let be the crack dip azimuth angle, ρ be the density of the isotropic background medium, ρ0 be the average density of the isotropic background medium, and V be the density of the crack. p0 V is the average longitudinal wave velocity of the isotropic background medium. S0 V is the average transverse wave velocity of the isotropic background medium. p V is the longitudinal wave velocity of the isotropic background medium. S Let EI0 be the transverse wave velocity of the isotropic background medium, EI0 be the average wave impedance of the isotropic background medium, and δ be the wave velocity of the isotropic background medium. N Let Δδ be the crack normal weakness parameter of the isotropic background medium. N The difference in the normal weakness parameter of the cracks in two adjacent isotropic background media is represented by δ. V Let Δδ be the vertical weakness parameter of the crack in the isotropic background medium. V The difference in vertical weakness parameters of cracks in two adjacent isotropic background media, δ H The horizontal weakness parameter of the crack in the isotropic background medium. The coefficients are the second-order Fourier coefficients, φ well For interpretation results of fracture orientation in well logging or geological data, sign[·] represents the sign. For weakness parameters, Let g be the average of the squares of the transverse and longitudinal wave velocities of two adjacent isotropic background media, and g = (V S / V P ) 2 .

[0037] Optionally, the influencing factors are:

[0038] The uncertainty in the sign of the coefficients of the second-order term in the Fourier series expansion of elastic impedance, i.e. middle The uncertainty of the sign.

[0039] According to a third aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0040] At least one processor; and,

[0041] A memory communicatively connected to the at least one processor; wherein,

[0042] The memory stores instructions executable by the at least one processor, which, when executed by the at least one processor, enables the at least one processor to perform the method for reducing 90-degree ambiguity in crack development orientation prediction as described in any of the first aspects.

[0043] According to a fourth aspect of the invention, a non-transitory computer-readable storage medium is provided, which stores computer instructions for causing a computer to execute the method for reducing 90-degree ambiguity in crack development orientation prediction as described in any of the first aspects.

[0044] The beneficial effects of this invention are as follows: This invention determines the influencing factors that lead to 90-degree ambiguity in fracture development azimuth prediction through the azimuth elastic wave impedance equation, and then analyzes the influencing factors based on the fracture azimuth interpretation results of well logging data or geological data and the fracture-induced azimuth anisotropy model, determines the strategy to weaken the 90-degree ambiguity in fracture development azimuth prediction, realizes the weakening of 90-degree ambiguity in fracture azimuth prediction, and improves the reliability of fracture azimuth seismic description.

[0045] The system of the present invention has other features and advantages that will be apparent from or will be set forth in detail in the accompanying drawings and following detailed description, which together serve to explain the particular principles of the invention. Attached Figure Description

[0046] The above and other objects, features and advantages of the present invention will become more apparent from the accompanying drawings, in which like reference numerals generally denote like parts.

[0047] Figure 1 A flowchart illustrating the steps of a method for reducing 90-degree ambiguity in crack development orientation prediction according to the present invention is shown.

[0048] Figure 2 a and Figure 2 b shows the second-order coefficients of the Fourier series expansion of elastic impedance according to Embodiment 2 of the present invention. A diagram illustrating the relationship.

[0049] Figure 3 A rose plot of the well perimeter fracture azimuth seismic prediction results according to Embodiment 2 of the present invention is shown.

[0050] Figure 4 A rose plot of wellbore interpretation results for the azimuth of fractures according to Embodiment 2 of the present invention is shown. Detailed Implementation

[0051] The invention will now be described in more detail with reference to the accompanying drawings. While preferred 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 so that the invention will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.

[0052] like Figure 1 As shown, a method for reducing 90-degree ambiguity in crack development orientation prediction according to the present invention includes:

[0053] Based on the azimuth elastic wave impedance equation, the influencing factors that lead to the 90-degree ambiguity in the prediction of crack development azimuth are determined.

[0054] The influencing factors are analyzed based on the fracture orientation interpretation results of well logging data or geological data and the fracture-induced orientation anisotropy model.

[0055] Based on the analysis results, a strategy was determined to reduce the 90-degree ambiguity in the prediction of crack development orientation.

[0056] Specifically, this invention, based on the azimuth elastic wave impedance equation, fixes the incident angle and analyzes... As the observation azimuth changes, it can be observed that under a fixed incident angle θ, if At the observation azimuth angle equal to the crack orientation φ sym hour Reaching the maximum value; if At the observation azimuth angle equal to the crack orientation φ sym hour Reaching the minimum value; In practical applications, only the absolute values ​​of the coefficients of the second-order term of the Fourier series expansion of elastic impedance can be obtained, and the sign cannot be directly determined. Therefore, the uncertainty of the sign of the coefficients of the second-order term of the Fourier series expansion of elastic impedance leads to a 90° ambiguity in the crack direction estimation. That is, the influencing factor is the uncertainty of the sign of the coefficients of the second-order term of the Fourier series expansion of elastic impedance; based on the second term of the Fourier series reflection coefficient To investigate the azimuth variation, a crack-induced azimuth anisotropy model was designed. The absolute value of the second-order Fourier coefficients of this model is 10. 5 The fracture azimuth is 30°, and the interpretation of the fracture azimuth based on well logging data or geological data is as follows: Analyzing the symbols, we obtain and Based on the consistent sign analysis results, a strategy was determined to weaken the 90-degree ambiguity in the prediction of crack development orientation, namely, extracting... The symbol can be used to determine The true positive and negative signs are used to achieve crack location prediction with reduced 90-degree ambiguity, thereby improving the reliability of crack location seismic description.

[0057] In one example, the calculation expression for the azimuth elastic wave impedance equation is as follows:

[0058]

[0059] in,

[0060]

[0061] EI0=V p0 ×ρ0;

[0062]

[0063]

[0064] φ is the observation azimuth angle, θ is the incident angle of the seismic wave, φ sym Let be the crack dip azimuth angle, ρ be the density of the isotropic background medium, ρ0 be the average density of the isotropic background medium, and V be the density of the crack. p0 V is the average longitudinal wave velocity of the isotropic background medium. S0 V is the average transverse wave velocity of the isotropic background medium. P V is the longitudinal wave velocity of the isotropic background medium. S Let EI0 be the transverse wave velocity of the isotropic background medium, EI0 be the average wave impedance of the isotropic background medium, and δ be the wave velocity of the isotropic background medium. N Let Δδ be the crack normal weakness parameter of the isotropic background medium. N The difference in the normal weakness parameter of the cracks in two adjacent isotropic background media is represented by δ. V Let Δδ be the vertical weakness parameter of the crack in the isotropic background medium. V The difference in vertical weakness parameters of cracks in two adjacent isotropic background media, δ H The horizontal weakness parameter of the crack in the isotropic background medium. The coefficients are the second-order Fourier coefficients, φ well For interpretation results of fracture orientation in well logging or geological data, sign[·] represents the sign. For weakness parameters, Let g be the average of the squares of the transverse and longitudinal wave velocities of two adjacent isotropic background media, and g = (V S / V P ) 2 .

[0065] In one example, the influencing factors are:

[0066] The uncertainty in the sign of the coefficients of the second-order term in the Fourier series expansion of elastic impedance, i.e. middle The uncertainty of the sign.

[0067] In one example, the analysis result is:

[0068] and The symbols are consistent.

[0069] In one example, the strategy includes:

[0070] pass extract The symbol is used to determine The symbol is used to reduce the 90-degree ambiguity in predicting the orientation of crack development.

[0071] The present invention will be further described below with reference to the accompanying drawings and specific embodiments, but this is not intended to limit the invention. It should be noted that, unless otherwise specified, the embodiments and features described in the embodiments of the present invention can be combined with each other.

[0072] Example 1

[0073] This embodiment provides a method for reducing the 90-degree ambiguity in crack development orientation prediction, including:

[0074] The factors influencing the 90-degree ambiguity in crack development azimuth prediction are determined based on the azimuth elastic wave impedance equation; the calculation expression of the azimuth elastic wave impedance equation is as follows:

[0075]

[0076] in,

[0077]

[0078] EI0=V p0 ×ρ0;

[0079]

[0080] φ is the observation azimuth angle, θ is the incident angle of the seismic wave, φ sym Let be the crack dip azimuth angle, ρ be the density of the isotropic background medium, ρ0 be the average density of the isotropic background medium, and V be the density of the crack. p0 V is the average longitudinal wave velocity of the isotropic background medium. S0 V is the average transverse wave velocity of the isotropic background medium. p V is the longitudinal wave velocity of the isotropic background medium. SLet EI0 be the transverse wave velocity of the isotropic background medium, EI0 be the average wave impedance of the isotropic background medium, and δ be the wave velocity of the isotropic background medium. N Let Δδ be the crack normal weakness parameter of the isotropic background medium. N The difference in the normal weakness parameter of the cracks in two adjacent isotropic background media is represented by δ. V Let Δδ be the vertical weakness parameter of the crack in the isotropic background medium. V The difference in vertical weakness parameters of cracks in two adjacent isotropic background media, δ H The horizontal weakness parameter of the crack in the isotropic background medium. The coefficients are the second-order Fourier coefficients, φ well For interpretation results of fracture orientation in well logging or geological data, sign[·] represents the sign. For weakness parameters, Let g be the average of the squares of the transverse and longitudinal wave velocities of two adjacent isotropic background media, and g = (V S / V P ) 2 .

[0081] The influencing factor is the uncertainty in the sign of the coefficients of the second-order term in the Fourier series expansion of elastic impedance, i.e. middle The uncertainty of the sign.

[0082] The influencing factors were analyzed based on the fracture orientation interpretation results from well logging or geological data and the fracture-induced orientation anisotropy model; the absolute value of the second-order Fourier coefficient of the fracture-induced orientation anisotropy model was 10. 5 The crack is located at an angle of 30°.

[0083] The analysis results are as follows: and The symbols are consistent.

[0084] Based on the analysis results, a strategy was determined to reduce the 90-degree ambiguity in the prediction of crack development orientation.

[0085] pass extract The symbol is used to determine The symbol is used to reduce the 90-degree ambiguity in predicting the orientation of crack development.

[0086] Example 2

[0087] This embodiment provides a method for reducing the 90-degree ambiguity in crack development orientation prediction, including:

[0088] Based on the azimuth elastic wave impedance equation, at a fixed incident angle, through analysis... As the observation azimuth changes, it can be observed that under a fixed incident angle θ, if At the observation azimuth angle equal to the crack orientation φ sym hour Reaching the maximum value; if At the observation azimuth angle equal to the crack orientation φ sym hour The minimum value is reached; in practical applications, only the absolute value of the coefficients of the second-order term of the Fourier series expansion of elastic impedance can be obtained, and the sign cannot be directly determined. Therefore, the uncertainty of the sign of the coefficients of the second-order term of the Fourier series expansion of elastic impedance leads to the 90° ambiguity in the crack direction estimation. That is, the influencing factor is the uncertainty of the sign of the coefficients of the second-order term of the Fourier series expansion of elastic impedance.

[0089] The calculation expression for the azimuth elastic wave impedance equation is as follows:

[0090]

[0091] in,

[0092]

[0093] EI0=V p0 ×ρ0;

[0094]

[0095] φ is the observation azimuth angle, θ is the incident angle of the seismic wave, φ sym Let be the crack dip azimuth angle, ρ be the density of the isotropic background medium, ρ0 be the average density of the isotropic background medium, and V be the density of the crack. p0 V is the average longitudinal wave velocity of the isotropic background medium. S0 V is the average transverse wave velocity of the isotropic background medium. p V is the longitudinal wave velocity of the isotropic background medium. S Let EI0 be the transverse wave velocity of the isotropic background medium, EI0 be the average wave impedance of the isotropic background medium, and δ be the wave velocity of the isotropic background medium. N Let Δδ be the crack normal weakness parameter of the isotropic background medium. N The difference in the normal weakness parameter of the cracks in two adjacent isotropic background media is represented by δ. V Let Δδ be the vertical weakness parameter of the crack in the isotropic background medium. V The difference in vertical weakness parameters of cracks in two adjacent isotropic background media, δ H The horizontal weakness parameter of the crack in the isotropic background medium. The coefficients are the second-order Fourier coefficients, φ wellFor interpretation results of fracture orientation in well logging or geological data, sign[·] represents the sign. For weakness parameters, Let g be the average of the squares of the transverse and longitudinal wave velocities of two adjacent isotropic background media, and g = (V S / V P ) 2 .

[0096] (1) The coefficients of the second-order terms based on the Fourier series expansion of elastic impedance The inability to directly determine the sign of the crack leads to ambiguity in crack orientation estimation. Therefore, a crack-induced azimuth anisotropy model is designed, with the absolute value of its second-order Fourier coefficients being 10. 5 The fracture azimuth is 30°. Based on the azimuth variation of the second term of the Fourier series reflection coefficient, and the interpretation of the fracture azimuth based on well logging or geological data, [the following is applied]. Analyze the symbols, such as Figure 2 As shown in (a), when the signs of the second-order Fourier coefficients are positive, The value is greater than zero. For example, Figure 2 As shown in (b), when the signs of the second-order Fourier coefficients are negative, The value is less than zero. Therefore, The sign of the coefficient is consistent with that of the second-order Fourier coefficient. Based on this analysis, a strategy to weaken the 90-degree ambiguity in the prediction of crack development orientation was determined, namely, extracting... The symbol can be used to determine The true positive and negative signs are used to achieve crack location prediction with reduced 90-degree ambiguity, thereby improving the reliability of crack location seismic description.

[0097] In the prediction of a fractured reservoir in a certain basin, a method for reducing 90-degree ambiguity in fracture development azimuth prediction, as described in this embodiment, was applied. This fractured reservoir mainly develops two phases of faults, primarily third- and fourth-order faults. The first phase of faults mainly consists of mountain peaks and derived faults, all trending northeast-trending reverse faults. The second phase of faults mainly consists of high-steep fold belt-controlled faults and derived small faults, with smaller scale and displacement, all trending northwest-trending reverse faults. Imaging logging shows that horizontal fractures are developed at the bottom of well D1, with few high-angle fractures, which is consistent with the prediction results. Figure 3 This is a rose plot of the wellbore fracture azimuth seismic prediction results after using the weakening method for 90-degree ambiguity prediction of fracture development azimuth in this embodiment. Figure 4 Rose plot of well logging interpretation results for the azimuth of fractures around the well. Figure 3The results indicate that the main azimuths of the fracture azimuth prediction results include 345 degrees and 60 degrees, and the number of fractures with a main azimuth of 345 degrees should be much greater than the number of fractures with a main azimuth of 60 degrees. Well logging interpretation results show that the main azimuths of the perimeter fractures include 345 degrees and 60 degrees, which is consistent with the earthquake prediction results. Although the difference between the number of fractures with a main azimuth of 345 degrees and the number of fractures with a main azimuth of 60 degrees is not very large, the reason for this difference may be that the fractures with a main azimuth of 60 degrees are low-angle fractures. While earthquakes can predict the azimuth of low-angle fractures, because low-angle fractures contribute less to the change in earthquake amplitude azimuth, the predicted number of low-angle fractures is relatively smaller when the number is the same as that of high-angle fractures.

[0098] Example 3

[0099] This embodiment provides a device for reducing the 90-degree ambiguity in crack development orientation prediction, including:

[0100] The first determination module is used to determine the influencing factors that cause the 90-degree ambiguity in the prediction of crack development azimuth based on the azimuth elastic wave impedance equation.

[0101] The analysis module is used to analyze influencing factors based on the fracture orientation interpretation results and fracture-induced orientation anisotropy model based on well logging data or geological data.

[0102] The second determination module is used to determine a strategy to reduce the 90-degree ambiguity in the prediction of crack development orientation based on the analysis results.

[0103] The calculation expression for the azimuth elastic wave impedance equation is as follows:

[0104]

[0105] in,

[0106]

[0107] EI0=V p0 ×ρ0;

[0108]

[0109]

[0110] φ is the observation azimuth angle, θ is the incident angle of the seismic wave, φ sym Let be the crack dip azimuth angle, ρ be the density of the isotropic background medium, ρ0 be the average density of the isotropic background medium, and V be the density of the crack. p0 V is the average longitudinal wave velocity of the isotropic background medium. S0 V is the average transverse wave velocity of the isotropic background medium. p V is the longitudinal wave velocity of the isotropic background medium. SLet EI0 be the transverse wave velocity of the isotropic background medium, EI0 be the average wave impedance of the isotropic background medium, and δ be the wave velocity of the isotropic background medium. N Let Δδ be the crack normal weakness parameter of the isotropic background medium. N The difference in the normal weakness parameter of the cracks in two adjacent isotropic background media is represented by δ. V Let Δδ be the vertical weakness parameter of the crack in the isotropic background medium. V The difference in vertical weakness parameters of cracks in two adjacent isotropic background media, δ H The horizontal weakness parameter of the crack in the isotropic background medium. The coefficients are the second-order Fourier coefficients, φ well For interpretation results of fracture orientation in well logging or geological data, sign[·] represents the sign. For weakness parameters, Let g be the average of the squares of the transverse and longitudinal wave velocities of two adjacent isotropic background media, and g = (V S / V P ) 2 .

[0111] The influencing factors are:

[0112] The uncertainty in the sign of the coefficients of the second-order term in the Fourier series expansion of elastic impedance, i.e. middle The uncertainty of the sign.

[0113] Example 4

[0114] This disclosure also provides an electronic device, which includes:

[0115] At least one processor; and,

[0116] A memory communicatively connected to the at least one processor; wherein,

[0117] The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the method for reducing the 90-degree ambiguity in crack development orientation prediction in Embodiment 1.

[0118] An electronic device according to embodiments of the present disclosure includes a memory and a processor. The memory is used to store non-transitory computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

[0119] The processor may be a central processing unit (CPU) or other form of processing unit with data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of this disclosure, the processor is used to execute computer-readable instructions stored in the memory.

[0120] Those skilled in the art will understand that, in order to solve the technical problem of how to achieve a good user experience, this embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included within the protection scope of this disclosure.

[0121] For a detailed description of this embodiment, please refer to the corresponding descriptions in the foregoing embodiments, which will not be repeated here.

[0122] Example 5

[0123] This disclosure provides a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the method for reducing 90-degree ambiguity in crack development orientation prediction in Embodiment 1.

[0124] A computer-readable storage medium according to embodiments of the present disclosure stores non-transitory computer-readable instructions. When these non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the methods described in the foregoing embodiments of the present disclosure are performed.

[0125] The aforementioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or portable hard drive), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).

[0126] The various embodiments of the present invention have been described above. These descriptions are exemplary and not exhaustive, nor are they limited to the disclosed embodiments. Many modifications and variations will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for weakening the 90-degree ambiguity in predicting the orientation of crack development, characterized in that, include: Based on the azimuth elastic wave impedance equation, the influencing factors that lead to the 90-degree ambiguity in the prediction of crack development azimuth are determined. The influencing factors are analyzed based on the fracture orientation interpretation results of well logging data or geological data and the fracture-induced orientation anisotropy model. Based on the analysis results, a strategy was determined to reduce the 90-degree ambiguity in the prediction of crack development orientation.

2. The method for weakening the 90-degree ambiguity in crack development orientation prediction according to claim 1, characterized in that, The calculation expression for the azimuth elastic wave impedance equation is as follows: in, EI0=V p0 ×ρ0; φ is the observation azimuth angle, θ is the incident angle of the seismic wave, φ sym Let be the crack dip azimuth angle, ρ be the density of the isotropic background medium, ρ0 be the average density of the isotropic background medium, and V be the density of the crack. p0 V is the average longitudinal wave velocity of the isotropic background medium. S0 V is the average transverse wave velocity of the isotropic background medium. p V is the longitudinal wave velocity of the isotropic background medium. S Let EI0 be the transverse wave velocity of the isotropic background medium, EI0 be the average wave impedance of the isotropic background medium, and δ be the wave velocity of the isotropic background medium. N Let Δδ be the crack normal weakness parameter of the isotropic background medium. N The difference in the normal weakness parameter of the cracks in two adjacent isotropic background media is represented by δ. V Let Δδ be the vertical weakness parameter of the crack in the isotropic background medium. V The difference in vertical weakness parameters of cracks in two adjacent isotropic background media, δ H The horizontal weakness parameter of the crack in the isotropic background medium. The coefficients are the second-order Fourier coefficients, φ well For interpretation results of fracture orientation in well logging or geological data, sign[·] represents the sign. For weakness parameters, Let g be the average of the squares of the transverse and longitudinal wave velocities of two adjacent isotropic background media, and g = (V S / V P ) 2 .

3. The method for weakening the 90-degree ambiguity in crack development orientation prediction according to claim 2, characterized in that, The influencing factors are: The uncertainty in the sign of the coefficients of the second-order term in the Fourier series expansion of elastic impedance, i.e. middle The uncertainty of the sign.

4. The method for weakening the 90-degree ambiguity in crack development orientation prediction according to claim 3, characterized in that, The analysis results are as follows: and The symbols are consistent.

5. The method for weakening the 90-degree ambiguity in crack development orientation prediction according to claim 4, characterized in that, The strategy includes: pass extract The symbol is used to determine The symbol is used to reduce the 90-degree ambiguity in predicting the orientation of crack development.

6. An electronic device, characterized in that, The electronic device includes: At least one processor; and, A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor to enable the at least one processor to perform the method for reducing 90-degree ambiguity in crack development orientation prediction as described in any one of claims 1-5.

7. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer instructions for causing a computer to perform the method for reducing 90-degree ambiguity in crack development orientation prediction as described in any one of claims 1-5.

8. A device for reducing 90-degree ambiguity in predicting the orientation of crack development, characterized in that, include: The first determination module is used to determine the influencing factors that cause the 90-degree ambiguity in the prediction of crack development azimuth based on the azimuth elastic wave impedance equation. The analysis module is used to analyze the influencing factors based on the fracture orientation interpretation results of well logging data or geological data and the fracture-induced orientation anisotropy model. The second determination module is used to determine a strategy to reduce the 90-degree ambiguity in the prediction of crack development orientation based on the analysis results.

9. The device for weakening the 90-degree ambiguity in crack development orientation prediction according to claim 1, characterized in that, The calculation expression for the azimuth elastic wave impedance equation is as follows: in, EI0=V p0 ×ρ0; φ is the observation azimuth angle, θ is the incident angle of the seismic wave, φ sym Let be the crack dip azimuth angle, ρ be the density of the isotropic background medium, ρ0 be the average density of the isotropic background medium, and V be the density of the crack. p0 V is the average longitudinal wave velocity of the isotropic background medium. S0 V is the average transverse wave velocity of the isotropic background medium. p V is the longitudinal wave velocity of the isotropic background medium. S Let EI0 be the transverse wave velocity of the isotropic background medium, EI0 be the average wave impedance of the isotropic background medium, and δ be the wave velocity of the isotropic background medium. N Let Δδ be the crack normal weakness parameter of the isotropic background medium. N The difference in the normal weakness parameter of the cracks in two adjacent isotropic background media is represented by δ. V Let Δδ be the vertical weakness parameter of the crack in the isotropic background medium. V The difference in vertical weakness parameters of cracks in two adjacent isotropic background media, δ H The horizontal weakness parameter of the crack in the isotropic background medium. The coefficients are the second-order Fourier coefficients, φ well For interpretation results of fracture orientation in well logging or geological data, sign[·] represents the sign. For weakness parameters, Let g be the average of the squares of the transverse and longitudinal wave velocities of two adjacent isotropic background media, and g = (V S / V P ) 2 .

10. The device for reducing 90-degree ambiguity in crack development orientation prediction according to claim 1, characterized in that, The influencing factors are: The uncertainty in the sign of the coefficients of the second-order term in the Fourier series expansion of elastic impedance, i.e. middle The uncertainty of the sign.

Images