An anisotropic elastic wave field decomposition method, device and computer equipment

Abstract

This specification relates to the field of exploration geophysics, and more particularly to an anisotropic elastic wavefield decomposition method, apparatus, and computer equipment. The anisotropic elastic wavefield decomposition method includes processing the received source wavelet according to an anisotropic model to obtain an anisotropic elastic vector wavefield matrix; determining the scalar P-wave matrix and vector S-wave matrix based on the modified Helmholtz operator and the vector wavefield matrix; and performing anisotropic medium correction and vectorization processing on the scalar P-wave matrix and vector S-wave matrix to obtain the target vector P-wave matrix and target vector S-wave matrix corresponding to the anisotropic elastic wavefield. Using the embodiments in this specification, by processing the anisotropic elastic wavefield with the modified Helmholtz operator, and then performing anisotropic medium correction and vectorization processing, the calculation process is simplified and the computational cost is reduced.

Description

Technical Field

[0001] This specification relates to the field of exploration geophysics, and in particular to an anisotropic elastic wave field decomposition method, apparatus and computer equipment. Background Technology

[0002] Currently, during seismic exploration, P-wave and S-wave field separation is required. This separation operation plays a crucial role in seismic numerical simulation, multi-component seismic data processing, and interpretation.

[0003] In isotropic media, the P-wave and S-wave field matrices can be obtained by using Helmholtz decomposition (Helmholtz operator) and decoupling continuation equations. However, in anisotropic media, since the P-wave and S-wave are not completely parallel or perpendicular to the propagation direction of the wave vector, it is difficult to directly separate the P-wave and S-wave fields.

[0004] For P-wave and S-wave field separation in anisotropic media, methods include calculating the wave polarization direction in the wavenumber domain to achieve P-wave and S-wave field separation, and inversely transforming the wave polarization direction to the spatial domain and using spatial domain wavefield separation operators to achieve wavefield separation. However, both methods require multiple Fourier transforms, resulting in complex calculations and high computational costs.

[0005] How to reduce the computational cost of performing longitudinal and transverse wave field separation operations on anisotropic elastic wave fields in anisotropic media is a problem that urgently needs to be solved in the existing technology. Summary of the Invention

[0006] To address the problems in the prior art, embodiments of this specification provide an anisotropic elastic wave field decomposition method, apparatus, computer equipment, and storage medium. The modified Helmholtz operator is used to process the anisotropic elastic wave field, thereby performing anisotropic medium correction and vectorization processing, which simplifies the calculation process and reduces computational costs.

[0007] To solve the above-mentioned technical problems, the specific technical solution in this specification is as follows:

[0008] On the one hand, embodiments of this specification provide a method for decomposing anisotropic elastic wave fields, including,

[0009] Based on the anisotropic model, the received source wavelet is processed to obtain the anisotropic elastic vector wave field matrix;

[0010] Based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix, the scalar longitudinal wave matrix and the vector transverse wave matrix are determined; and

[0011] The scalar longitudinal wave matrix and the vector transverse wave matrix are corrected and vectorized for anisotropic media to obtain the target vector longitudinal wave matrix and the target vector transverse wave matrix corresponding to the anisotropic elastic vector wave field matrix.

[0012] Furthermore, determining the scalar longitudinal wave matrix and the vector transverse wave matrix based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix further includes:

[0013] A dot product is performed on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix to obtain the scalar longitudinal wave matrix; and

[0014] The vector transverse wave matrix is ​​obtained by performing a cross product between the modified Helmholtz operator and the anisotropic elastic vector wave field matrix.

[0015] Furthermore, the modified Helmholtz operator further includes,

[0016]

[0017] Among them, the For the modified Helmholtz operator, ε, γ, and δ are the parameters of anisotropy, and v p0 Let v be the longitudinal wave velocity of the medium along the axis of symmetry. s0 Let be the transverse wave velocity of the medium along the axis of symmetry.

[0018] Furthermore, the anisotropic medium correction and vectorization processing performed on the scalar longitudinal wave matrix and the vector transverse wave matrix to obtain the target vector longitudinal wave matrix and target vector transverse wave matrix corresponding to the anisotropic elastic vector wave field matrix further includes,

[0019] The scalar P-wave matrix and the vector S-wave matrix are corrected to obtain a corrected scalar P-wave matrix and a corrected vector S-wave matrix; and

[0020] The corrected scalar P-wave matrix and the corrected vector S-wave matrix are corrected and vectorized to obtain the target vector P-wave matrix and the target vector S-wave matrix.

[0021] Furthermore, the correction processing of the scalar P-wave matrix and the vector S-wave matrix to obtain the corrected scalar P-wave matrix and the corrected vector S-wave matrix further includes:

[0022] Based on the dispersion relation formula and the amplitude factor set, construct a set of transformation formulas;

[0023] Based on the set of transformation formulas, the scalar P-wave matrix and the vector S-wave matrix are transformed to obtain a wavenumber-domain scalar P-wave matrix and a wavenumber-domain vector S-wave matrix; and

[0024] Spatiotemporal transformation is performed on the wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix to obtain the corrected scalar P-wave matrix and the corrected vector S-wave matrix.

[0025] Furthermore, this set of conversion formulas further includes,

[0026]

[0027]

[0028]

[0029]

[0030] Among them, the The phase velocity of the longitudinal wave, the Let f(k) and g be the phase velocity of the transverse wave. p (k) and the g s (k) represent the amplitude factor, where k is the wave number, and v is the wave number. p0 The longitudinal wave velocity of the medium along the axis of symmetry, v s0 The transverse wave velocity of the medium along the axis of symmetry, the For the wavenumber domain correction scalar longitudinal wave matrix, the For the wavenumber domain correction vector shear wave matrix, the A p Let be the unit vector along the direction of the longitudinal wave particle vibration velocity, and the... is the anisotropic elastic vector wave field matrix in the wavenumber domain.

[0031] Furthermore, the correction and vectorization processing of the corrected scalar P-wave matrix and the corrected vector S-wave matrix to obtain the target vector P-wave matrix and the target vector S-wave matrix further includes,

[0032]

[0033]

[0034] Wherein, the U qPamp The target vector longitudinal wave matrix, the U qSamp The target vector shear wave matrix, the v p0 The longitudinal wave velocity of the medium along the axis of symmetry, v s0 The transverse wave velocity of the medium along the axis of symmetry, the For the modified Helmholtz operator, the qP amp To correct the scalar longitudinal wave matrix, and the qS amp To correct the vector shear wave matrix.

[0035] On the other hand, embodiments of this specification also provide an anisotropic elastic vector wave field decomposition apparatus, including,

[0036] The first processing unit is used to process the received source wavelet according to the anisotropic model to obtain the anisotropic elastic vector wave field matrix.

[0037] The determining unit is used to determine the scalar longitudinal wave matrix and the vector transverse wave matrix based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix; and

[0038] The second processing unit is used to perform anisotropic medium correction and vectorization processing on the scalar longitudinal wave matrix and the vector transverse wave matrix to obtain the target vector longitudinal wave matrix and the target vector transverse wave matrix corresponding to the anisotropic elastic vector wave field matrix.

[0039] On the other hand, embodiments of this specification also provide a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the above-described method.

[0040] On the other hand, embodiments of this specification also provide a computer-readable storage medium having computer instructions stored thereon, which, when executed by a processor, implement the above-described method.

[0041] Using the embodiments in this specification, by introducing a modified Helmholtz operator, dispersion relation formula, and amplitude factor set, the anisotropic elastic vector wave field is subjected to anisotropic medium P-wave separation, correction, and vectorization processing to obtain the target vector P-wave matrix and target vector S-wave matrix in the anisotropic medium. This simplifies the calculation process and reduces computational costs. Attached Figure Description

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

[0043] Figure 1 The figure shown is a schematic diagram of an implementation system for an anisotropic elastic wave field decomposition method according to an embodiment of this specification.

[0044] Figure 2 The diagram shown is a flowchart of an anisotropic elastic wave field decomposition method according to an embodiment of this specification.

[0045] Figure 3 The diagram shown is a flowchart of an anisotropic elastic wave field decomposition method according to another embodiment of this specification.

[0046] Figure 4 The diagram shown is a schematic diagram of an anisotropic elastic wave field decomposition method according to an embodiment of this specification.

[0047] Figure 5A The diagram shown is a schematic of another embodiment of an anisotropic elastic wave field decomposition method in this specification.

[0048] Figure 5B The diagram shown is a schematic diagram of a homogeneous anisotropic medium model according to an embodiment of this specification.

[0049] Figure 5C The diagram shown is a schematic representation of a high-order staggered mesh finite difference scheme according to an embodiment of this specification.

[0050] Figure 5D The diagram shown is a schematic diagram of the horizontal component of a particle vibration velocity according to an embodiment of this specification.

[0051] Figure 5E The diagram shown is a schematic diagram of the vertical component of a particle vibration velocity according to an embodiment of this specification.

[0052] Figure 5F The diagram shown is a schematic diagram of a calibrated scalar longitudinal wave according to an embodiment of this specification;

[0053] Figure 5G The diagram shown is a schematic diagram of a corrected vector shear wave according to an embodiment of this specification;

[0054] Figure 5H The diagram shown is a schematic diagram of the horizontal component of the longitudinal wave particle vibration velocity according to an embodiment of this specification.

[0055] Figure 5I The diagram shown is a schematic diagram of the vertical component of the longitudinal wave particle vibration velocity according to an embodiment of this specification.

[0056] Figure 5J The diagram shown is a schematic diagram of the horizontal component of the vibration velocity of a transverse wave particle according to an embodiment of this specification.

[0057] Figure 5K The diagram shown is a schematic diagram of the vertical component of the vibration velocity of a transverse wave particle according to an embodiment of this specification.

[0058] Figure 6 The diagram shown is a structural schematic of an anisotropic elastic wave field decomposition device according to an embodiment of this specification.

[0059] Figure 7 This is a schematic diagram of the structure of a computer device according to an embodiment of this specification.

[0060] [Explanation of Labels in the Attached Image]

[0061] 101. Data acquisition terminal;

[0062] 102. Server;

[0063] 103. User terminal;

[0064] 401. Anisotropic elastic vector wave field matrix;

[0065] 402. Set of vector formulas for the wavenumber domain;

[0066] 403. Wavenumber domain correction vectorization formula set;

[0067] 410. Dispersion Relationship Correction Algorithm;

[0068] 420. Spatiotemporal domain transformation algorithm;

[0069] 431. Target vector longitudinal wave matrix;

[0070] 432. Target vector shear wave matrix;

[0071] 501. Source wavelet;

[0072] 502. Anisotropic elastic vector wave field matrix;

[0073] 503. Correction scalar longitudinal wave matrix;

[0074] 504. Correction vector shear wave matrix;

[0075] 510. Anisotropic model;

[0076] 520. Correction algorithm;

[0077] 530. Correction vectorization algorithm;

[0078] 540. The modified Helmholtz operator;

[0079] 551. Target vector longitudinal wave matrix;

[0080] 552. Target vector shear wave matrix;

[0081] 610. First processing unit;

[0082] 620. Determine the unit;

[0083] 630. Second processing unit;

[0084] 702. Computer equipment;

[0085] 704. Processing equipment;

[0086] 706. Storage resources;

[0087] 708. Drive mechanism;

[0088] 710. Input / Output Module;

[0089] 712. Input devices;

[0090] 714. Output devices;

[0091] 716. Presentation equipment;

[0092] 718. Graphical User Interface;

[0093] 720. Network interface;

[0094] 722. Communication link;

[0095] 724. Communication bus. Detailed Implementation

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

[0097] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, apparatus, product, or device that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or devices.

[0098] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0099] Figure 1The diagram illustrates an implementation system for an anisotropic elastic wave field decomposition method according to an embodiment of this specification. The system may include a data acquisition terminal 101 and a server 102, which communicate via a network. This network may include a Local Area Network (LAN), a Wide Area Network (WAN), the Internet, or a combination thereof, and is connected to a website, user equipment (e.g., a computing device), and a backend system. The source wavelet can be input to the server 102 through the data acquisition terminal 101. The data acquisition terminal 101 may be, for example, an electronic device or a sensor. The source wavelet can be obtained from signals acquired by a sensor or input by a user through an electronic device. After receiving the source wavelet, the server 102 uses an anisotropic model to obtain the anisotropic elastic wave field and then uses a modified Helmholtz operator to transform, correct, and vectorize the anisotropic elastic wave field, obtaining the target vector P-wave matrix and the target vector S-wave matrix in the anisotropic medium.

[0100] If the acquisition terminal 101 is an electronic device, the target vector P-wave matrix and the target vector S-wave matrix can be sent to the electronic device. When the acquisition terminal 101 is a sensor, a schematic diagram of an implementation system for an anisotropic elastic wave field decomposition method may also include a user terminal 103, to which the target vector P-wave matrix and the target vector S-wave matrix are sent. The user terminal 103 and the server 102 communicate via a network, which may include a local area network (LAN), a wide area network (WAN), the Internet, or a combination thereof, and is connected to a website, user equipment (e.g., a computing device), and a backend system.

[0101] Alternatively, server 102 may be a node of a cloud computing system (not shown in the figure), or each server 102 may be a separate cloud computing system comprising multiple computers interconnected by a network and operating as a distributed processing system.

[0102] In an optional embodiment, the user terminal 103 may include an electronic device. The electronic device included in the user terminal 103 and the aforementioned data acquisition terminal 101 may be, but are not limited to, smartphones, data acquisition devices, desktop computers, tablets, laptops, smart speakers, digital assistants, augmented reality (AR) / virtual reality (VR) devices, smart wearable devices, and other similar electronic devices. Optionally, the operating system running on the electronic device may include, but is not limited to, Android, iOS, Linux, Windows, etc.

[0103] In addition, it should be noted that, Figure 1 The example shown is merely one application environment provided in this manual. In actual applications, it may also include multiple acquisition terminals 101 and multiple user terminals 103. This manual does not impose any restrictions.

[0104] like Figure 2 The diagram shows a flowchart of an anisotropic elastic wave field decomposition method according to an embodiment of this specification. While the anisotropic elastic wave field decomposition process is depicted in this figure, it can include more or fewer steps based on conventional or non-inventive methods. The order of steps listed in the embodiment is merely one possible execution order among many and does not represent the only possible order. In actual system or device products, the method can be executed sequentially or in parallel according to the embodiment or the accompanying drawings. Specifically, as shown... Figure 2 As shown, the method may include:

[0105] S210, based on the anisotropic model, the received source wavelet is processed to obtain the anisotropic elastic vector wave field matrix;

[0106] S220, based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix, determines the scalar longitudinal wave matrix and the vector transverse wave matrix;

[0107] S230 performs anisotropic medium correction and vectorization processing on the scalar longitudinal wave matrix and vector transverse wave matrix to obtain the target vector longitudinal wave matrix and target vector transverse wave matrix corresponding to the anisotropic elastic vector wave field.

[0108] Using the embodiments in this specification, an anisotropic model is employed to digitize the received source wavelet, yielding an anisotropic elastic vector wavefield matrix. A modified Helmholtz operator is then introduced to process this vector wave matrix, performing anisotropic medium correction and vectorization to obtain the target vector P-wave matrix and target vector S-wave matrix in the anisotropic medium. This simplifies the calculation process and reduces computational costs.

[0109] According to one embodiment of this specification, the anisotropic model includes anisotropic parameters, a first-order velocity-stress equation, and a staggered-grid finite-difference method. The staggered-grid finite-difference method is used to solve the first-order velocity-stress equation to obtain the corresponding matrix.

[0110] After receiving the source wavelet, an anisotropic model is determined for digitization in the anisotropic medium. Using the determined anisotropic model, the source wavelet is digitized to obtain the anisotropic elastic vector wave field matrix.

[0111] The original Helmholtz operator performs well in isotropic media for transverse and longitudinal wave field separation, but its performance is less than ideal in anisotropic media. Therefore, the original Helmholtz operator is modified to obtain a modified Helmholtz operator. This modification can be any treatment applied to the original Helmholtz operator to achieve ideal results in transverse and longitudinal wave field separation in anisotropic media. For example, it could involve normalization, weighting, or decomposition of the original Helmholtz operator.

[0112] The formula for separating transverse and longitudinal waves is obtained by processing the modified Helmholtz operator and the anisotropic elastic vector wave field matrix to obtain the scalar longitudinal wave matrix and the vector transverse wave matrix. This formula can be a pre-defined formula for separating transverse and longitudinal waves from the anisotropic elastic vector wave field matrix.

[0113] After determining the scalar P-wave matrix and the vector S-wave matrix, a pre-defined algorithm for correction and vectorization is obtained. Using this algorithm, the scalar P-wave matrix and the vector S-wave matrix are corrected and vectorized to obtain the target vector P-wave matrix and the target vector S-wave matrix corresponding to the anisotropic elastic vector wave field matrix in the anisotropic medium. This pre-defined algorithm for correction and vectorization may, for example, include a correction and vectorization algorithm based on dispersion relations.

[0114] According to another embodiment of this specification, in order to ensure that the modified Helmholtz operator can be completely degenerated into the original isotropic Helmholtz operator, the original Helmholtz operator is divided by the difference between the squares of the longitudinal and transverse wave velocities to obtain a modified Helmholtz operator for anisotropic media that eliminates amplitude distortion. The modified Helmholtz operator is shown in the following formula (1).

[0115]

[0116] in, For the modified Helmholtz operator, ε, γ, and δ are the parameters of anisotropy, v p0 Let v be the longitudinal wave velocity of the medium along the axis of symmetry. s0 Let be the transverse wave velocity of the medium along the axis of symmetry.

[0117] According to another embodiment of this specification, determining the scalar longitudinal wave matrix and the vector transverse wave matrix based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix includes: performing a dot product operation on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix to obtain the scalar longitudinal wave matrix; and performing a cross product operation on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix to obtain the vector transverse wave matrix.

[0118] In other words, the formulas for separating transverse and longitudinal waves can be represented by formulas (2) and (3) as follows.

[0119]

[0120]

[0121] Where qP is the scalar longitudinal wave matrix. For the modified Helmholtz operator, U is the anisotropic elastic vector wave field matrix, and qS is the vector transverse wave matrix.

[0122] According to another embodiment of this specification, the anisotropic medium correction and vectorization processing of the scalar longitudinal wave matrix and the vector transverse wave matrix to obtain the target vector longitudinal wave matrix and the target vector transverse wave matrix corresponding to the anisotropic elastic vector wave field matrix includes: performing correction processing on the scalar longitudinal wave matrix and the vector transverse wave matrix to obtain the corrected scalar longitudinal wave matrix and the corrected vector transverse wave matrix; and performing correction and vectorization processing on the corrected scalar longitudinal wave matrix and the corrected vector transverse wave matrix to obtain the target vector longitudinal wave matrix and the target vector transverse wave matrix.

[0123] The correction and vectorization process for the determined scalar P-wave matrix and vector S-wave matrix may include, for example, first correcting the scalar P-wave matrix and vector S-wave matrix, and then correcting and vectorizing the corrected scalar P-wave matrix and corrected vector S-wave matrix to obtain the target vector P-wave matrix and target vector S-wave matrix. Specifically, the above-mentioned dispersion relation correction and vectorization algorithm may also include a dispersion relation-based correction algorithm and a dispersion relation-based correction vectorization algorithm. The dispersion relation-based correction algorithm is used to correct the scalar P-wave matrix and vector S-wave matrix to obtain the corrected scalar P-wave matrix and corrected vector S-wave matrix. The dispersion relation-based correction vectorization algorithm is used to correct and vectorize the corrected scalar P-wave matrix and corrected vector S-wave matrix to obtain the target vector P-wave matrix and target vector S-wave matrix.

[0124] Figure 3The diagram shows a flowchart of an anisotropic elastic vector wave field decomposition method according to another embodiment of this specification. While the process of anisotropic elastic vector wave field decomposition is described in this figure, it may include more or fewer operational steps based on conventional or non-inventive methods. Specifically, as shown... Figure 3 As shown, the method may include:

[0125] S321, Construct a set of transformation formulas based on the dispersion relation formula and the amplitude factor set;

[0126] S322, Based on the set of transformation formulas, the scalar longitudinal wave matrix and the vector transverse wave matrix are transformed to obtain the wavenumber domain scalar longitudinal wave matrix and the wavenumber domain vector transverse wave matrix.

[0127] S323 performs a spatiotemporal transformation on the wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix to obtain the corrected scalar P-wave matrix and the corrected vector S-wave matrix.

[0128] Using the embodiments in this specification, by introducing dispersion relation formulas and amplitude factor sets, the scalar P-wave matrix and vector S-wave matrix are corrected for anisotropic media, resulting in corrected scalar P-wave matrix and corrected vector S-wave matrix in anisotropic media. This avoids Fourier transform, simplifying the calculation process and reducing computational costs.

[0129] According to another embodiment of this specification, the dispersion relation is the relationship between the wavenumber modulus and the angular frequency and phase velocity, specifically as shown in the following formula (4).

[0130]

[0131] Where k is the wavenumber modulus, ω is the angular frequency, and v is the phase velocity.

[0132] The amplitude factor set includes multiple existing amplitude factors. It should be noted that amplitude factors are introduced here because the phase velocities of P-waves and S-waves are difficult to directly convert into the spatial domain. Therefore, an amplitude factor set is needed to process the phase velocities of the P-waves and S-waves so that they can be converted into the spatial domain for correction processing.

[0133] According to another embodiment of this specification, based on the dispersion relation formula and the amplitude factor set, the conversion formula set in the conversion formula set includes the following formulas (5)-(8).

[0134]

[0135]

[0136]

[0137]

[0138] in, Let be the phase velocity of the longitudinal wave. Let f(k) and g be the phase velocity of the transverse wave. p (k) and the g s (k) represents the amplitude factor, k is the wave number, and v p0 Let v be the longitudinal wave velocity of the medium along the axis of symmetry. s0 Let V be the transverse wave velocity of the medium along the axis of symmetry. For wavenumber domain correction of scalar P-wave matrix, For the wavenumber domain correction vector transverse wave matrix, A p Let be the unit vector in the direction of longitudinal wave propagation, and is the anisotropic elastic vector wave field matrix in the wavenumber domain.

[0139] Based on the set of conversion formulas, the scalar P-wave matrix and the vector S-wave matrix are converted to obtain the wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix. Specifically, the wavenumber domain conversion is performed on the scalar P-wave matrix and the vector S-wave matrix to obtain the wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix; the set of conversion formulas is then substituted into the wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix to obtain the wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix. The wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix are shown in formulas (9) and (10) below, respectively.

[0140]

[0141]

[0142] in, The wavenumber domain scalar longitudinal wave matrix, Let A be the wavenumber domain vector transverse wave matrix. p (k) is A p The modulus of (k), iA p (k) is the expression for the modified Helmholtz operator in the wavenumber domain, where k is the wavenumber and A is the wavenumber. p Let be the unit vector in the direction of longitudinal wave propagation, and is the anisotropic elastic vector wave field matrix in the wavenumber domain.

[0143] To obtain the corrected scalar P-wave matrix and the corrected vector S-wave matrix by performing a spatiotemporal transformation on the wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix, specifically, the spatiotemporal domain corrected scalar P-wave matrix formula and the spatiotemporal domain corrected vector S-wave matrix formula are obtained by performing a spatiotemporal domain transformation on the wavenumber domain scalar P-wave matrix and the wavenumber domain corrected vector S-wave matrix formula; the corrected scalar P-wave matrix and the corrected vector S-wave matrix formula are then solved to obtain the corrected scalar P-wave matrix and the corrected vector S-wave matrix. The spatiotemporal domain corrected scalar P-wave matrix formula and the spatiotemporal domain corrected vector S-wave matrix formula can be, for example, as shown in the following formulas (11) and (12). The spatiotemporal domain corrected scalar P-wave matrix formula and the spatiotemporal domain corrected vector S-wave matrix formula can be solved, for example, by using the staggered grid finite difference method.

[0144]

[0145]

[0146] Among them, qP amp To correct the scalar longitudinal wave matrix, qS amp To correct the vector shear wave matrix, v p0 Let v be the longitudinal wave velocity of the medium along the axis of symmetry. s0 Let V be the transverse wave velocity of the medium along the axis of symmetry. U is the modified Helmholtz operator, and U is the anisotropic elastic vector wave field matrix.

[0147] Figure 4 The diagram shown is a schematic diagram of an anisotropic elastic wave field decomposition method according to an embodiment of this specification.

[0148] According to another embodiment of this specification, the correction and vectorization of the corrected scalar P-wave matrix and the corrected vector S-wave matrix to obtain the target vector P-wave matrix and the target vector S-wave matrix can be specifically achieved by: performing wavenumber domain transformation and decomposition on the anisotropic elastic vector wave field matrix to obtain the wavenumber domain vector P-wave field matrix and the wavenumber domain vector S-wave field matrix; correcting the wavenumber domain vector P-wave field matrix and the wavenumber domain vector S-wave field matrix using dispersion relations to obtain the dispersion-corrected P-wave matrix and the dispersion-corrected S-wave matrix; substituting the wavenumber domain target P-wave matrix and the wavenumber domain target S-wave field matrix obtained based on the amplitude factor, the wavenumber domain target P-wave matrix, and the wavenumber domain target S-wave field matrix into the dispersion-corrected P-wave matrix and the dispersion-corrected S-wave matrix and performing a spatiotemporal domain transformation to obtain the formulas for the target vector P-wave matrix and the target vector S-wave matrix; and then solving the formulas for the target vector P-wave matrix and the target vector S-wave matrix to obtain the target vector P-wave matrix and the target vector S-wave matrix. The formulas for the target vector P-wave matrix and the target vector S-wave matrix can be solved using the staggered grid finite difference method.

[0149] The wavenumber domain vector longitudinal wave field matrix and the wavenumber domain vector transverse wave field matrix can be represented, for example, by the following formulas (13) and (14).

[0150]

[0151]

[0152] in, The wavenumber domain vector longitudinal wave field matrix, Let A be the wavenumber domain vector transverse wave field matrix. p Let be the unit vector in the direction of longitudinal wave propagation, and is the anisotropic elastic vector wave field matrix in the wavenumber domain.

[0153] The dispersion-corrected longitudinal wave matrix and the dispersion-corrected transverse wave matrix can be represented by, for example, the following formula (15).

[0154]

[0155] in, The wavenumber domain vector longitudinal wave field matrix, Let f(k) and g be the wavenumber domain vector transverse wave field matrices. p (k) and the g s (k) represent the amplitude factors, v p0 Let v be the longitudinal wave velocity of the medium along the axis of symmetry. s0 Let qP be the transverse wave velocity of the medium along the axis of symmetry. ampTo correct the scalar longitudinal wave matrix, qS amp To correct the vector shear wave matrix, k is the wave number, A p Let A be the unit vector in the direction of longitudinal wave propagation, and let A be the vector in the direction of longitudinal wave propagation. p (k) represents the direction of longitudinal wave propagation in an anisotropic medium, where A p (k) is a function of k. Among them, C 11 C 44 C 33 and C 13 These are the stiffness matrix parameters in anisotropic media, k x and k z Let be the wavenumbers in the x and Z directions, respectively, and D = {(C 11 -C 44 )k x 2 -(C 33 -C 44 )k z 2 +4(C 13 +C 44 ) 2 k x 2 k z 2}

[0156] The target vector P-wave matrix and the target vector S-wave matrix in the wavenumber domain can be represented, for example, by the following formulas (16) and (17). That is, substitute formulas (16) and (17) into the above formula (15).

[0157]

[0158]

[0159] in, The target vector longitudinal wave matrix in the wavenumber domain. Let f(k) and g be the target vector transverse wave matrix in the wavenumber domain. p (k) and the g s (k) represent the amplitude factor, and The wavenumber domain vector longitudinal wave field matrix, is the wavenumber domain vector transverse wave field matrix.

[0160] The target P-wave matrix and target S-wave matrix in the wavenumber domain, obtained based on the amplitude factor, the target P-wave matrix in the wavenumber domain, and the target S-wave matrix in the wavenumber domain, are substituted into the dispersion-corrected P-wave matrix and dispersion-corrected S-wave matrix. For example, this could involve substituting the target P-wave matrix and target S-wave matrix in the wavenumber domain, obtained based on the amplitude factor, the target P-wave matrix in the wavenumber domain, and the target S-wave matrix in the wavenumber domain, into the dispersion-corrected P-wave matrix and dispersion-corrected S-wave matrix, and multiplying both sides of the equation by the imaginary number "i".

[0161] According to another embodiment of this specification, the correction scalar P-wave matrix and the correction vector S-wave matrix are corrected and vectorized to obtain the target vector P-wave matrix and the target vector S-wave matrix formulas, for example, the following formulas (18) and (19).

[0162]

[0163]

[0164] Among them, U qPamp U is the target vector longitudinal wave matrix. qSamp Let v be the target vector shear wave matrix. p0 Let v be the longitudinal wave velocity of the medium along the axis of symmetry. s0 Let V be the transverse wave velocity of the medium along the axis of symmetry. For the modified Helmholtz operator, qP amp To correct the scalar longitudinal wave matrix, and qS amp To correct the vector shear wave matrix.

[0165] like Figure 4 As shown, the dispersion relation correction vectorization algorithm may include, for example, performing wavenumber domain transformation and decomposition on the anisotropic elastic vector wavefield matrix 401 to obtain a set of wavenumber domain vector formulas 402, including the wavenumber domain vector longitudinal wavefield matrix and the wavenumber domain vector transverse wavefield matrix. This set of wavenumber domain vector formulas is then input into the dispersion relation correction algorithm 410 to obtain a set of wavenumber domain correction vectorization formulas 403. Furthermore, this set of wavenumber domain correction vectorization formulas 403 is input into the spatiotemporal domain transformation algorithm 420 to obtain the target vector longitudinal wave matrix 431 and the target vector transverse wave matrix 432.

[0166] Specifically, the dispersion relation correction algorithm 410 may include, for example, the above formulas (15)-(17). That is, the target vector P-wave matrix and target vector S-wave matrix in the wavenumber domain, obtained based on the amplitude factor, the target vector P-wave matrix in the wavenumber domain, and the target vector S-wave matrix in the wavenumber domain, are substituted into the dispersion correction P-wave matrix and dispersion correction S-wave matrix [substituting formulas (16) and (17) into the above formula (15)] to obtain the wavenumber domain correction vectorization formula set 403.

[0167] The spatiotemporal domain transformation algorithm 420 is any algorithm that can convert the wavenumber domain formula set into a spatiotemporal domain formula set. After performing spatiotemporal domain transformation on the wavenumber domain correction vectorization formula set 403 using the spatiotemporal domain transformation algorithm 420, the target vector P-wave matrix formula and the target vector S-wave matrix formula are obtained; then, the target vector P-wave matrix formula and the target vector S-wave matrix formula are solved to obtain the target vector P-wave matrix 431 and the target vector S-wave matrix 432.

[0168] Figure 5A The diagram shown is a schematic of another embodiment of an anisotropic elastic wave field decomposition method in this specification. Figure 5B The diagram shown is a schematic diagram of a homogeneous anisotropic medium model according to an embodiment of this specification. Figure 5C The diagram shown is a schematic representation of a high-order staggered mesh finite difference scheme according to an embodiment of this specification. Figure 5D The diagram shown is a schematic representation of the horizontal component of a particle vibration velocity according to an embodiment of this specification. Figure 5D In the figure, the horizontal axis “vx” represents the horizontal component of the particle vibration velocity; Figure 5E The diagram shown is a schematic representation of the vertical component of a particle vibration velocity according to an embodiment of this specification. Figure 5E In the figure, the horizontal axis “vz” represents the vertical component of the particle vibration velocity; Figure 5F The diagram shown is a schematic representation of a calibrated scalar longitudinal wave according to an embodiment of this specification. Figure 5F In the diagram, the horizontal axis “P” represents the correction scalar longitudinal wave; Figure 5G The diagram shown is a schematic representation of a corrected vector shear wave according to an embodiment of this specification. Figure 5F In the diagram, the horizontal axis “S” represents the correction vector shear wave; Figure 5H The diagram shown is a schematic representation of the horizontal component of the longitudinal wave particle vibration velocity according to an embodiment of this specification. Figure 5H In the figure, the horizontal axis “PVX” represents the horizontal component of the vibration velocity of the longitudinal wave particles; Figure 5I The diagram shown is a schematic representation of the vertical component of the longitudinal wave particle vibration velocity according to an embodiment of this specification. Figure 5I In the figure, the horizontal axis “PVZ” represents the vertical component of the longitudinal wave particle vibration velocity; Figure 5J The diagram shown is a schematic representation of the horizontal component of the vibration velocity of a transverse wave particle according to an embodiment of this specification. Figure 5JIn the figure, the horizontal axis “SVX” represents the horizontal component of the vibration velocity of the transverse wave particles; Figure 5K The diagram shown is a schematic representation of the vertical component of the vibration velocity of a transverse wave particle according to an embodiment of this specification. Figure 5K In the figure, the horizontal axis “SVZ” represents the vertical component of the vibration velocity of the transverse wave particles.

[0169] According to another embodiment of this specification, after receiving the source wavelet 501, the source wavelet 501 is processed using an anisotropic model 510 to obtain an anisotropic elastic vector wave field matrix 502. Based on the anisotropic elastic vector wave field matrix 502 and the modified Helmholtz operator 540, a scalar P-wave matrix and a vector S-wave matrix are obtained. The scalar P-wave matrix and the vector S-wave matrix are input into a correction algorithm 520 to obtain a corrected scalar P-wave matrix 503 and a corrected vector S-wave matrix 504. Then, the corrected scalar P-wave matrix 503 and the corrected vector S-wave matrix 504 are input into a correction vectorization algorithm 530 to obtain a target vector P-wave matrix 551 and a target vector S-wave matrix 552. The correction algorithm 520 can be, for example, the above formulas (11) and (12), and the correction vectorization algorithm 530 can be, for example, the above formulas (18) and (19). In addition, the correction algorithm 520 can include, for example, a correction algorithm based on dispersion relations. The correction vectorization algorithm 530 may include, for example, a correction vectorization algorithm based on dispersion relation.

[0170] For example, the P-wave and S-wave field separation performance described in this specification was tested using an anisotropic model (homogeneous two-dimensional VTI medium model). The example was a test with a horizontal and vertical range of 4 km, and P-wave and S-wave velocities and densities of 3000 m / s, 2000 m / s, and 1000 kg / m², respectively. 3 The anisotropy parameters ε and δ are 0.15 and 0.1, respectively. First, based on the elastic wave equation obtained from a source wavelet with a dominant frequency of 30Hz, a staggered-grid finite-difference scheme is used for numerical simulation to obtain the vibration velocity vector wave field of the particles. The spatial sampling interval is set to 10m, the time sampling interval is set to 1ms, and the source location is (2000m, 2000m). Figure 5B As shown, construct as follows Figure 5C The model of a homogeneous anisotropic medium is shown. This leads to the following: Figure 5D The diagram shows the horizontal component of the total particle vibration velocity at 0.55s. Figure 5E The diagram shows the vertical component of the total particle vibration velocity at 0.55s.

[0171] Furthermore, substituting the scalar P-wave matrix and the vector S-wave matrix into the above formulas (11) and (12), and using the spatial staggered grid difference scheme defined by the staggered grid finite difference method, the formulas obtained by substituting the scalar P-wave matrix and the vector S-wave matrix into the above formulas (11) and (12) are discretized to obtain the corrected scalar P-wave matrix and the corrected vector S-wave matrix. Specifically, Figure 5F The diagram shown is a schematic of the corrected scalar P-wave corresponding to the corrected scalar P-wave matrix. Figure 5G This is a schematic diagram of the correction vector shear wave corresponding to the correction vector shear wave matrix.

[0172] Specifically, by substituting the scalar longitudinal wave matrix and the vector transverse wave matrix into the above formulas (11) and (12), the formula obtained is shown in the following formula (20).

[0173]

[0174] Among them, qP amp To correct the scalar longitudinal wave matrix, and These are the correction vector transverse wave matrices of the particle vibration velocity vector field along the x, y, and z directions, respectively. x v y v z R1, R2, and R3 represent the components of the particle vibration velocity vector field along the x, y, and z directions, respectively, and R1, R2, and R3 are the stiffness matrix parameters of the anisotropic medium, respectively.

[0175] Discretize the formulas obtained by substituting the scalar longitudinal wave matrix and the vector transverse wave matrix into the above formulas (11) and (12) to obtain the formulas corresponding to the corrected scalar longitudinal wave matrix and the corrected vector transverse wave matrix as shown in the following formula (21).

[0176]

[0177] in, These represent the forward and backward staggered grid difference schemes along the x-direction, respectively. These represent the forward and backward staggered grid difference schemes along the y-direction, respectively. Let qP represent the forward and backward staggered mesh difference schemes along the z-direction, respectively. amp To correct the scalar longitudinal wave matrix, and These are the correction vector transverse wave matrices of the particle vibration velocity vector field along the x, y, and z directions, respectively.

[0178] Furthermore, substituting the corrected scalar P-wave matrix and the corrected vector S-wave matrix into the above formulas (18) and (19), and using the spatial staggered grid difference scheme defined by the staggered grid finite difference method, the formulas obtained by substituting the corrected scalar P-wave matrix and the corrected vector S-wave matrix into the above formulas (18) and (19) are discretized to obtain the target vector P-wave matrix and the target vector S-wave matrix. Specifically, Figure 5H The diagram shows the horizontal components of the target vector P-wave corresponding to the horizontal components in the target vector P-wave matrix. Figure 5I This is a schematic diagram of the vertical component of the target vector longitudinal wave, corresponding to the vertical component in the target vector longitudinal wave matrix. Figure 5J The diagram shown is a schematic representation of the horizontal components of the target vector shear wave, corresponding to the horizontal components in the target vector shear wave matrix. Figure 5K This is a schematic diagram of the vertical component of the target vector shear wave, corresponding to the vertical component in the target vector shear wave matrix. From the above... Figures 5H-5K The obtained analysis results show that this specification can obtain the phase-preserving target vector P-wave matrix and target vector S-wave matrix, and the P-wave and S-wave have almost no crosstalk.

[0179] Substituting the corrected scalar P-wave matrix and the corrected vector S-wave matrix into the above formulas (18) and (19), the resulting formula is shown in formula (22) below. The formula corresponding to the corrected scalar P-wave matrix and the corrected vector S-wave matrix is ​​shown in formula (22) below.

[0180]

[0181] in, and These represent the vibration velocities of the longitudinal wave particles, v and v', respectively. qPamp and the vibration velocity v of the transverse wave particles qSamp x, y, z components, qP amp To correct the scalar longitudinal wave matrix, and These are the correction vector transverse wave matrices of the particle vibration velocity vector field along the x, y, and z directions, respectively. p0 Let v be the longitudinal wave velocity of the medium along the axis of symmetry. s0 Let r1 be the transverse wave velocity of the medium along the axis of symmetry, and r2 and r3 be the stiffness matrix parameters of the anisotropic medium, respectively.

[0182] The formulas obtained by substituting the correction scalar P-wave matrix and the correction vector S-wave matrix into the above formulas (18) and (19) are discretized to obtain the formulas corresponding to the target vector P-wave matrix and the target vector S-wave matrix as shown in the following formula (23).

[0183]

[0184] in, and These represent the vibration velocities of the longitudinal wave particles, v and v', respectively. qPamp and the vibration velocity v of the transverse wave particles qSamp x, y, z components, qP amp To correct the scalar longitudinal wave matrix, and These are the correction vector transverse wave matrices of the particle vibration velocity vector field along the x, y, and z directions, respectively. p0 Let v be the longitudinal wave velocity of the medium along the axis of symmetry. s0 Let r1 be the transverse wave velocity of the medium along the axis of symmetry, and r2 and r3 be the stiffness matrix parameters of the anisotropic medium, respectively.

[0185] Figure 6 The diagram shown is a structural schematic of an anisotropic elastic wave field decomposition device according to an embodiment of this specification. Figure 6 As shown, including,

[0186] The first processing unit 610 is used to process the received source wavelet according to the anisotropic model to obtain the anisotropic elastic vector wave field matrix.

[0187] Determining unit 620, used to determine the scalar P-wave matrix and vector S-wave matrix based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix; and

[0188] The second processing unit 630 is used to perform anisotropic medium correction and vectorization processing on the scalar longitudinal wave matrix and the vector transverse wave matrix to obtain the target vector longitudinal wave matrix and the target vector transverse wave matrix corresponding to the anisotropic elastic wave field.

[0189] Since the principle of the above-mentioned device in solving the problem is similar to that of the above-mentioned method, the implementation of the above-mentioned device can refer to the implementation of the above-mentioned method, and the repeated parts will not be described again.

[0190] like Figure 7The diagram illustrates the structure of a computer device according to an embodiment of this specification. The apparatus described in this specification can be the computer device in this embodiment, performing the methods described above. The computer device 702 may include one or more processing devices 704, such as one or more central processing units (CPUs), each of which can implement one or more hardware threads. The computer device 702 may also include any storage resource 706 for storing information of any kind, such as code, settings, data, etc. Without limitation, for example, the storage resource 706 may include any one or more combinations of the following: any type of RAM, any type of ROM, flash memory, hard disk, optical disk, etc. More generally, any storage resource can use any technology to store information. Furthermore, any storage resource can provide volatile or non-volatile retention of information. Further, any storage resource may represent a fixed or removable component of the computer device 702. In one case, when the processing device 704 executes associated instructions stored in any storage resource or combination of storage resources, the computer device 702 can perform any operation of the associated instructions. The computer device 702 also includes one or more drive mechanisms 708 for interacting with any storage resource, such as a hard disk drive mechanism, an optical disk drive mechanism, etc.

[0191] Computer device 702 may also include an input / output module 710 (I / O) for receiving various inputs (via input device 712) and providing various outputs (via output device 714). A specific output mechanism may include a presentation device 716 and an associated graphical user interface (GUI) 718. In other embodiments, the input / output module 710 (I / O), input device 712, and output device 714 may be omitted, and the device may function solely as a computer device within a network. Computer device 702 may also include one or more network interfaces 720 for exchanging data with other devices via one or more communication links 722. One or more communication buses 724 couple the components described above together.

[0192] Communication link 722 can be implemented in any way, such as via a local area network, a wide area network (e.g., the Internet), a point-to-point connection, or any combination thereof. Communication link 722 may include any combination of hardwired links, wireless links, routers, gateway functions, name servers, etc., governed by any protocol or combination of protocols.

[0193] This specification also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method.

[0194] This specification also provides a computer program product, which includes a computer program that, when executed by a processor, implements the above-described method.

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

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

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

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

[0199] The above specific embodiments further illustrate the purpose, technical solutions, and beneficial effects of this specification. It should be understood that the above are merely specific embodiments of this specification and are not intended to limit the scope of protection of this specification. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this specification should be included within the scope of protection of this specification.

Claims

1. A method for decomposing anisotropic elastic wave fields, characterized in that, include: Based on the anisotropic model, the received source wavelet is processed to obtain the anisotropic elastic vector wave field matrix; Based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix, the scalar longitudinal wave matrix and the vector transverse wave matrix are determined. as well as The scalar longitudinal wave matrix and the vector transverse wave matrix are corrected and vectorized for anisotropic media to obtain the target vector longitudinal wave matrix and the target vector transverse wave matrix corresponding to the anisotropic elastic vector wave field matrix. The determination of the scalar longitudinal wave matrix and the vector transverse wave matrix based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix includes: A dot product is performed on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix to obtain the scalar longitudinal wave matrix; and The vector transverse wave matrix is ​​obtained by performing a cross product operation on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix. The modified Helmholtz operator includes: Among them, the For the modified Helmholtz operator, the stated The above and stated These are the parameters for anisotropy, respectively. The longitudinal wave velocity of the medium along the axis of symmetry, and the... Let be the transverse wave velocity of the medium along the axis of symmetry.

2. The method according to claim 1, characterized in that, The step of performing anisotropic medium correction and vectorization processing on the scalar longitudinal wave matrix and the vector transverse wave matrix to obtain the target vector longitudinal wave matrix and the target vector transverse wave matrix corresponding to the anisotropic elastic vector wave field matrix includes: The scalar P-wave matrix and the vector S-wave matrix are corrected to obtain a corrected scalar P-wave matrix and a corrected vector S-wave matrix; and The corrected scalar P-wave matrix and the corrected vector S-wave matrix are corrected and vectorized to obtain the target vector P-wave matrix and the target vector S-wave matrix.

3. The method according to claim 2, characterized in that, The step of performing correction processing on the scalar P-wave matrix and the vector S-wave matrix to obtain the corrected scalar P-wave matrix and the corrected vector S-wave matrix includes: Based on the dispersion relation formula and the amplitude factor set, construct a set of transformation formulas; Based on the set of transformation formulas, the scalar P-wave matrix and the vector S-wave matrix are transformed to obtain a wavenumber-domain scalar P-wave matrix and a wavenumber-domain vector S-wave matrix; and Spatiotemporal transformation is performed on the wavenumber domain scalar P-wave matrix and the wavenumber domain vector S-wave matrix to obtain the corrected scalar P-wave matrix and the corrected vector S-wave matrix.

4. The method according to claim 3, characterized in that, The set of conversion formulas includes: Among them, the The phase velocity of the longitudinal wave, the The phase velocity of the transverse wave is... , and stated These are the amplitude factors, respectively. For wave number, the The longitudinal wave velocity of the medium along the axis of symmetry is... The transverse wave velocity of the medium along the axis of symmetry, the For the wavenumber domain correction scalar longitudinal wave matrix, the For the wavenumber domain correction vector shear wave matrix, the Let be the unit vector along the direction of the longitudinal wave particle vibration velocity, and the... is the anisotropic elastic vector wave field matrix in the wavenumber domain.

5. The method according to claim 2, characterized in that, The step of correcting and vectorizing the corrected scalar P-wave matrix and the corrected vector S-wave matrix to obtain the target vector P-wave matrix and the target vector S-wave matrix includes: Among them, the The target vector longitudinal wave matrix, the The target vector shear wave matrix, the The longitudinal wave velocity of the medium along the axis of symmetry is... The transverse wave velocity of the medium along the axis of symmetry, the For the modified Helmholtz operator, the stated To correct the scalar longitudinal wave matrix, and the aforementioned To correct the vector shear wave matrix.

6. An anisotropic elastic wave field decomposition device, characterized in that, include: The first processing unit is used to process the received source wavelet according to the anisotropic model to obtain the anisotropic elastic vector wave field matrix. The determining unit is used to determine the scalar longitudinal wave matrix and the vector transverse wave matrix based on the modified Helmholtz operator and the vector wave field matrix; as well as The second processing unit is used to perform anisotropic medium correction and vectorization processing on the scalar longitudinal wave matrix and the vector transverse wave matrix to obtain the target vector longitudinal wave matrix and the target vector transverse wave matrix corresponding to the anisotropic elastic vector wave field matrix. The determination of the scalar longitudinal wave matrix and the vector transverse wave matrix based on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix includes: A dot product is performed on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix to obtain the scalar longitudinal wave matrix; and The vector transverse wave matrix is ​​obtained by performing a cross product operation on the modified Helmholtz operator and the anisotropic elastic vector wave field matrix. The modified Helmholtz operator includes: Among them, the For the modified Helmholtz operator, the stated The above and stated These are the parameters for anisotropy, respectively. The longitudinal wave velocity of the medium along the axis of symmetry, and the... Let be the transverse wave velocity of the medium along the axis of symmetry.

7. 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 computer program, it implements the method of any one of claims 1-5.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, performs the method of any one of claims 1-5.

Images