An anisotropic elastic wave decoupling method, device and computer equipment

Abstract

The present specification relates to the technical field of exploration geophysics, and particularly relates to an anisotropic elastic wave decoupling method and device and computer equipment. The anisotropic elastic wave decoupling method comprises: determining a set of Thomson parameters included in an anisotropic model according to a received wave field decomposition request to be decomposed, the wave field decomposition request to be decomposed comprising a wave field to be decomposed; converting the set of Thomson parameters to obtain an initial set of elastic parameters; performing P-S wave velocity separation processing on the initial set of elastic parameters to obtain a target P-wave elastic parameter set and a target S-wave elastic parameter set; and substituting the target P-wave elastic parameter set and the target S-wave elastic parameter set into the anisotropic model to process the wave field to be decomposed to obtain a target P-wave matrix and a target S-wave matrix. The embodiments of the present specification simplify the P-S wave field decomposition process and reduce the calculation cost.

Description

TECHNICAL FIELD

[0001] The present specification relates to the technical field of exploration geophysics, and in particular to an anisotropic elastic wave decoupling method and device and computer equipment. BACKGROUND

[0002] With the development of computer technology and the increasingly complex oil exploration targets, the requirement for the accuracy of seismic data processing is higher and higher, thereby increasing the demand for multi-component elastic reverse time migration (ERTM) imaging technology, and the core technology of the ERTM is the wave field separation of elastic waves. In isotropic media, the wave field decomposition method of elastic waves includes Helmholtz decomposition, traveling wave separation method, and decoupling continuation equation, etc. However, in anisotropic media, the polarization direction of the wave field is not completely parallel or perpendicular to the propagation direction, and therefore, the wave field separation method of isotropy is not applicable.

[0003] In anisotropic media, the current mainstream wave field separation method is the wave number domain wave field separation method, which first needs to calculate the polarization direction of the wave number domain, and then projects the wave field in the wave number domain to the polarization direction of the wave field, and then obtains the separated wave field through inverse Fourier transform. Therefore, the wave number domain wave field separation method needs to perform a large number of Fourier transforms, resulting in a huge amount of calculation.

[0004] In the process of wave field separation in anisotropic media, how to avoid Fourier inverse transform, simplify the P-S wave field decomposition process, and reduce the calculation cost is a problem to be solved in the prior art. SUMMARY

[0005] To solve the problems in the prior art, the embodiments of the present specification provide an anisotropic elastic wave decoupling method, device, computer equipment and storage medium, which is applied to a wave field caused by an earthquake, converts a set of Thomson parameters included in an anisotropic model to determine a set of initial elastic parameters, and then performs P-S wave velocity separation processing on the initial elastic parameters to obtain a set of target P-wave elastic parameters and a set of target S-wave elastic parameters for the separation of P-S wave fields of a wave field to be decomposed. The P-S wave field decomposition process is simplified, and the calculation cost is reduced.

[0006] To solve the above technical problems, the specific technical solutions of the present specification are as follows:

[0007] On the one hand, the embodiments of the present specification provide an anisotropic elastic wave decoupling method applied to a wave field caused by an earthquake, which includes,

[0008] According to the received wave field decomposition request to be decomposed, a set of Thomsen parameters included in the anisotropy model is determined, and the wave field to be decomposed is included in the wave field decomposition request to be decomposed;

[0009] Conversion is performed on the set of Thomsen parameters to obtain an initial set of elastic parameters;

[0010] Horizontal and vertical wave velocity separation processing is performed on the initial set of elastic parameters to obtain a target set of longitudinal wave elastic parameters and a target set of transverse wave elastic parameters; and

[0011] The target set of longitudinal wave elastic parameters and the target set of transverse wave elastic parameters are substituted into the anisotropy model to process the wave field to be decomposed to obtain a target longitudinal wave matrix and a target transverse wave matrix.

[0012] Further, the initial set of elastic parameters further includes:

[0013]

[0014]

[0015]

[0016]

[0017]

[0018]

[0019] wherein the C 11 , the C 12 , the C 13 , the C 33 , the C 44 , the C 55 , and the C 66 are initial elastic parameters, the ε, the γ, and the δ are parameters of anisotropy, the v p0 is a longitudinal wave velocity of a medium along a symmetry axis direction, the v s0 is a transverse wave velocity of the medium along the symmetry axis direction, and the ρ is a medium density field.

[0020] Further, the horizontal and vertical wave velocity separation processing on the initial set of elastic parameters to obtain the target set of longitudinal wave elastic parameters and the target set of transverse wave elastic parameters further includes,

[0021] At least one elastic parameter to be separated is determined from the initial set of elastic parameters, and the elastic parameter to be separated includes an initial elastic parameter determined based on a transverse wave velocity and a longitudinal wave velocity;

[0022] determining a corresponding P-S wave velocity separation step for each of the at least one elastic parameter to be separated; and

[0023] performing P-S wave velocity separation processing on each of the at least one elastic parameter to be separated based on the P-S wave velocity separation step, to obtain the target set of P-wave elastic parameters and the target set of S-wave elastic parameters.

[0024] Further, in the case where the at least one elastic parameter to be separated is determined based on a square root operation, the P-S wave velocity separation step further comprises,

[0025] performing square root processing on the at least one elastic parameter to be separated, to obtain a separated elastic parameter; and

[0026] performing P-S wave velocity separation processing on the separated elastic parameter, to obtain a target P-wave elastic parameter and a target S-wave elastic parameter, for obtaining the target set of P-wave elastic parameters and the target set of S-wave elastic parameters.

[0027] Further, the target set of P-wave elastic parameters further comprises:

[0028]

[0029]

[0030]

[0031]

[0032]

[0033]

[0034] wherein the the the the the the and the are target P-wave elastic parameters, the ε and the δ are anisotropic parameters, the v p0 is a P-wave velocity of the medium along the direction of the symmetry axis, the v s0 is an S-wave velocity of the medium along the direction of the symmetry axis, and the ρ is a density field of the medium.

[0035] Further, the target set of S-wave elastic parameters further comprises:

[0036]

[0037]

[0038]

[0039]

[0040]

[0041]

[0042] Among them, the The The The The The and stated These are the target shear wave elastic parameters, where γ and δ are anisotropic parameters, and v is... p0 The longitudinal wave velocity of the medium along the axis of symmetry, v s0 Let ρ be the transverse wave velocity of the medium along the axis of symmetry, and let ρ be the density field of the medium.

[0043] Furthermore, the target P-wave elastic parameter set and the target S-wave elastic parameter set are substituted into the anisotropic model to process the wave field to be decomposed, and the resulting target P-wave matrix and target S-wave matrix further include,

[0044] Substituting the target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set into the anisotropic model, we obtain the target anisotropic equation set; and

[0045] The target anisotropy equations are solved using the staggered grid high-order finite difference method to obtain the target P-wave matrix and the target S-wave matrix.

[0046] On the other hand, embodiments of this specification also provide an anisotropic elastic wave decoupling device, applied to wave fields caused by earthquakes, including,

[0047] The first determining unit is used to determine the set of Thomson parameters included in the anisotropic model based on the received decomposition request of the wave field to be decomposed, wherein the decomposition request of the wave field to be decomposed includes the wave field to be decomposed.

[0048] A conversion unit is used to convert the Thomson parameter set to obtain an initial elasticity parameter set;

[0049] A separation unit is used to perform transverse and longitudinal wave velocity separation processing on the initial set of elastic parameters to obtain a target longitudinal wave elastic parameter set and a target transverse wave elastic parameter set; and

[0050] The processing unit is configured to substitute the target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set into the anisotropic model, and process the wave field to be decomposed to obtain a target longitudinal wave matrix and a target transverse wave matrix.

[0051] In another aspect, the embodiments of the present specification also provide a computer device, comprising a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and the processor implements the above method when executing the computer program.

[0052] In another aspect, the embodiments of the present specification also provide a computer readable storage medium having computer instructions stored thereon, and the computer instructions are executed by a processor to implement the above method.

[0053] According to the embodiments of the present specification, when receiving a wave field decomposition request to be decomposed, the set of Thomson parameters included in the anisotropic model is determined. Then, the initial elastic parameter set is determined by converting the set of Thomson parameters. Then, the target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set are obtained by performing the transverse wave and longitudinal wave velocity separation processing on the initial elastic parameter. Finally, the target longitudinal wave matrix and the target transverse wave matrix are determined by using the target longitudinal wave elastic parameter set, the target transverse wave elastic parameter set, and the anisotropic model. Thus, in the wave field separation process in the anisotropic medium, the Fourier inverse transform is avoided, the transverse wave and longitudinal wave field decomposition process is simplified, and the calculation cost is reduced. BRIEF DESCRIPTION OF DRAWINGS

[0054] In order to more clearly illustrate the technical solutions in the embodiments of the present specification or the prior art, the drawings needed in the embodiments or prior art description will be briefly introduced below. Obviously, the drawings in the following description are only some embodiments of the present specification, and other drawings can be obtained by those skilled in the art without creative labor.

[0055] Figure 1 An implementation system schematic diagram of an anisotropic elastic wave decoupling method according to an embodiment of the present specification is shown.

[0056] Figure 2 A flowchart of an anisotropic elastic wave decoupling method according to an embodiment of the present specification is shown.

[0057] Figure 3 A flowchart of an anisotropic elastic wave decoupling method according to another embodiment of the present specification is shown.

[0058] Figure 4 A principle diagram of an anisotropic elastic wave decoupling method according to an embodiment of the present specification is shown.

[0059] Figure 5AA schematic diagram of an staggered grid finite difference format of an embodiment of the present specification is shown.

[0060] Figure 5B A schematic diagram of a wave field snapshot corresponding to the horizontal component of the full wave field particle vibration velocity field of an embodiment of the present specification is shown.

[0061] Figure 5C A schematic diagram of a wave field snapshot corresponding to the horizontal component of the longitudinal wave particle vibration velocity field of an embodiment of the present specification is shown.

[0062] Figure 5D A schematic diagram of a wave field snapshot corresponding to the horizontal component of the transverse wave particle vibration velocity field of an embodiment of the present specification is shown.

[0063] Figure 5E A schematic diagram of a wave field snapshot corresponding to the vertical component of the full wave field particle vibration velocity field of an embodiment of the present specification is shown.

[0064] Figure 5F A schematic diagram of a wave field snapshot corresponding to the vertical component of the longitudinal wave particle vibration velocity field of an embodiment of the present specification is shown.

[0065] Figure 5G A schematic diagram of a wave field snapshot corresponding to the vertical component of the transverse wave particle vibration velocity field of an embodiment of the present specification is shown.

[0066] Figure 5H A schematic diagram of a single channel waveform curve corresponding to the horizontal component of the longitudinal wave particle vibration velocity field of an embodiment of the present specification is shown.

[0067] Figure 5I A schematic diagram of a single channel waveform curve corresponding to the horizontal component of the transverse wave particle vibration velocity field of an embodiment of the present specification is shown.

[0068] Figure 5J A schematic diagram of a single channel waveform curve corresponding to the vertical component of the longitudinal wave particle vibration velocity field of an embodiment of the present specification is shown.

[0069] Figure 5K A schematic diagram of a single channel waveform curve corresponding to the vertical component of the transverse wave particle vibration velocity field of an embodiment of the present specification is shown.

[0070] Figure 5L A schematic diagram of a model parameter of an embodiment of the present specification is shown.

[0071] Figure 5M A schematic diagram of a reverse time migration imaging result of an embodiment of the present specification is shown.

[0072] Figure 6 A schematic diagram of the structure of an anisotropic elastic wave decoupling device of an embodiment of the present specification is shown.

[0073] Figure 7 FIG. 1 is a structural schematic diagram of a computer device according to an embodiment of the present specification.

[0074]

Explanation of reference signs

[0075] 101, acquisition terminal;

[0076] 102, server;

[0077] 103, user terminal;

[0078] 401, wave field to be decomposed;

[0079] 410, anisotropy model;

[0080] 420, initial elastic parameter set;

[0081] 431, target longitudinal wave elastic parameter set;

[0082] 432, target transverse wave elastic parameter set;

[0083] 440, target anisotropy equation set;

[0084] 451, target longitudinal wave matrix;

[0085] 452, target transverse wave matrix;

[0086] 501, longitudinal wave velocity parameter field;

[0087] 502, transverse wave velocity parameter field;

[0088] 503, ε parameter field;

[0089] 504, δ parameter field;

[0090] 511, source longitudinal wave detection longitudinal wave image;

[0091] 512, source longitudinal wave detection transverse wave image;

[0092] 513, source transverse wave detection longitudinal wave image;

[0093] 514, source transverse wave detection transverse wave image;

[0094] 610, first determination unit;

[0095] 620, conversion unit;

[0096] 630, separation unit;

[0097] 640, processing unit;

[0098] 702, computer device;

[0099] 704, processing device;

[0100] 706, storage resource;

[0101] 708, drive mechanism;

[0102] 710, input / output module;

[0103] 712, input device;

[0104] 714, output device;

[0105] 716, presentation device;

[0106] 718, graphical user interface;

[0107] 720, network interface;

[0108] 722, communication link;

[0109] 724, communication bus. DETAILED DESCRIPTION

[0110] The technical solutions in the embodiments of the present specification will be described clearly and completely in combination with the drawings in the embodiments of the present specification. Obviously, the described embodiments are only part of the embodiments of the present specification, rather than all the embodiments. Based on the embodiments in the present specification, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present specification.

[0111] It should be noted that the terms "first", "second", and the like in the specification and claims of the present specification and the above-described drawings are used to distinguish similar objects, and do not necessarily have to be used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances, so that the embodiments of the present specification described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion, for example, a process, method, device, product, or apparatus including a series of steps or units does not have to be limited to those steps or units clearly listed, but can include other steps or units not clearly listed or inherent to these processes, methods, products, or apparatuses.

[0112] It should be noted that the steps shown in the flowchart of the 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 herein can be executed in an order different from that described herein.

[0113] Figure 1As shown in the embodiment of the present specification, a schematic diagram of an implementation system of an anisotropic elastic wave decoupling method can include a collection terminal 101 and a server 102. The collection terminal 101 and the server 102 communicate through a network. The network can include a local area network (LAN), a wide area network (WAN), the Internet, or a combination thereof, and is connected to a website, a user device (such as a computing device), and a backend system. The wave field to be decomposed can be input to the server 102 through the collection terminal 101. The collection terminal 101 can be an electronic device or a sensor. The wave field to be decomposed can be obtained from a signal collected by a sensor or input by a user through an electronic device. After receiving the wave field to be decomposed, the server 102 converts and separates the P and S wave velocities of the Thomson parameter set included in the anisotropic model to obtain a target P wave elastic parameter set and a target S wave elastic parameter set. Then, the target P wave elastic parameter set, the target S wave elastic parameter set, and the anisotropic model are used to determine a target P wave matrix and a target S wave matrix, which can be used to generate decomposed P and S wave images.

[0114] If the collection terminal 101 is an electronic device, the decomposed P and S wave images corresponding to the wave field to be decomposed can be sent to the electronic device. If the collection terminal 101 is a sensor, the schematic diagram of an implementation system of an anisotropic elastic wave decoupling method can further include a user terminal 103 to which the decomposed P and S wave images are sent. The user terminal 103 and the server 102 communicate through a network. The network can include a local area network (LAN), a wide area network (WAN), the Internet, or a combination thereof, and is connected to a website, a user device (such as a computing device), and a backend system.

[0115] Optionally, the server 102 can be a node of a cloud computing system (not shown in the figure), or each server 102 can be a separate cloud computing system, including multiple computers interconnected by a network and working as a distributed processing system.

[0116] In an optional embodiment, the user terminal 103 can include an electronic device. The electronic device included in the user terminal 103 and the electronic device that can be included in the acquisition terminal 101 are not limited to the type of electronic device such as a smartphone, an acquisition device, a desktop computer, a tablet computer, a notebook computer, a smart speaker, a digital assistant, an augmented reality (AR) / virtual reality (VR) device, a smart wearable device, and the like. Optionally, the operating system running on the electronic device can include, but is not limited to, an Android system, an IOS system, Linux, Windows, and the like.

[0117] In addition, it should be noted that Figure 1 The illustrated is only one application environment provided by the present specification, and in actual application, a plurality of acquisition terminals 101, a plurality of user terminals 103 can also be included, and the present specification is not limited.

[0118] As Figure 2 The flowchart shown is a flowchart of an anisotropic elastic wave decoupling method according to an embodiment of the present specification. In the present figure, the decoupling process performed on the wave field caused by the earthquake is described, but more or fewer operation steps can be included based on conventional or non-inventive labor. The order of steps listed in the embodiment is only one of the many execution orders, and does not represent the only execution order. In actual system or device product execution, the method order shown in the embodiment or the drawing can be executed in sequence or in parallel. Specifically as Figure 2 The method can include:

[0119] S210, according to the received wave field decomposition request to be decomposed, determining a set of Thomson parameters included in the anisotropic model, the wave field decomposition request to be decomposed including a wave field to be decomposed;

[0120] S220, converting the set of Thomson parameters to obtain an initial set of elastic parameters;

[0121] S230, performing P-S wave velocity separation processing on the initial set of elastic parameters to obtain a target P-wave elastic parameter set and a target S-wave elastic parameter set;

[0122] S240, substituting the target P-wave elastic parameter set and the target S-wave elastic parameter set into the anisotropic model to process the wave field to be decomposed to obtain a target P-wave matrix and a target S-wave matrix.

[0123] The anisotropic model is determined by determining a Thomsen parameter set included in the anisotropic model when receiving a wave field decomposition request to be decomposed. Then, conversion is performed on the Thomsen parameter set to determine an initial elastic parameter set. Then, the initial elastic parameter is subjected to P-S wave velocity separation processing to obtain a target P wave elastic parameter set and a target S wave elastic parameter set. The target P wave elastic parameter set, the target S wave elastic parameter set, and the anisotropic model are used to determine a target P wave matrix and a target S wave matrix. Thus, in the wave field separation process in the anisotropic medium, the Fourier inverse transform is avoided, the P-S wave field decomposition process is simplified, and the calculation cost is reduced.

[0124] According to one embodiment of the present specification, the anisotropic model represents any model for determining corresponding S wave and P wave matrices in an anisotropic medium by performing digital processing on a wave field to be decomposed. For example, the anisotropic model includes a first-order velocity stress equation and an equation solving algorithm. The equation solving algorithm is used to solve the first-order velocity stress equation based on the wave field to be decomposed, thereby obtaining the corresponding S wave and P wave matrices. The first-order velocity stress equation includes P-S wave velocities and a Thomsen parameter set (Thomsen parameter field).

[0125] The conversion of the Thomsen parameter set to obtain the initial elastic parameter set can be specifically based on the mapping relationship between each Thomsen parameter in the Thomsen parameter set and the initial elastic parameter shown in the following formula (1). The Thomsen parameter set includes a plurality of Thomsen parameters, and each Thomsen parameter is converted to determine the initial elastic parameter corresponding to each Thomsen parameter to obtain the initial elastic parameter set.

[0126]

[0127] wherein α0 and β0 are the vertical velocities of quasi P and quasi S waves along the vertical direction, v p0 = α0, v s0 = β0, ε, δ, and γ are parameters of anisotropy, which are all dimensionless Thomsen parameters, ε is a parameter describing the strength of P wave anisotropy, when ε is 0, the P wave has no anisotropy, δ is a parameter describing the relationship between the vertical velocity v p0 and the moveout velocity v pn , that is, δ is a physical quantity affecting the P wave velocity near the symmetry axis of the anisotropic medium, γ is a parameter describing the strength of S wave anisotropy or the degree of S wave splitting in the propagation process of the S wave, C 11 , C 12 , C 13 , C 33 , C 44 , C 55 , and C 66respectively.

[0128] According to the above formula (1), after the conversion of the plurality of Thomson parameters included in the Thomson parameter set, the initial elastic parameters obtained are functions of the shear wave velocity and the compressional wave velocity.

[0129] For the plurality of initial elastic parameters included in the initial elastic parameter set, the shear wave velocity and the compressional wave velocity included in the initial elastic parameter are separated to obtain the corresponding target compressional wave elastic parameter and the target shear wave elastic parameter. For example, for the initial elastic parameter being the sum of the shear wave velocity function and the compressional wave velocity function, the shear wave velocity and the compressional wave velocity are separated for the initial elastic parameter to obtain the target compressional wave elastic parameter being the compressional wave velocity function and the target shear wave elastic parameter being the shear wave velocity function. That is, the plurality of compressional wave elastic parameters included in the generated target compressional wave elastic parameter set are functions of the compressional wave velocity, and the plurality of shear wave elastic parameters included in the target shear wave elastic parameter set are functions of the shear wave velocity.

[0130] After the target compressional wave elastic parameter set and the target shear wave elastic parameter set are determined, the compressional wave elastic parameter set and the target shear wave elastic parameter set are substituted into the first-order velocity-stress equation of the anisotropic model to obtain a target anisotropic equation set, and then an equation solving algorithm is used to solve the target anisotropic equation set to obtain a target compressional wave matrix and a target shear wave matrix, and the data in the target compressional wave matrix and the target shear wave matrix are functions of time. Thus, the decomposed compressional wave and shear wave image can be obtained. It should be noted that the image is an image corresponding to the parameter field, also known as a wave field snapshot.

[0131] According to another embodiment of the present specification, according to the above formula (1), after the conversion of the plurality of Thomson parameters included in the Thomson parameter set, the initial elastic parameter set obtained includes, for example, the following formula (2) to formula (7).

[0132]

[0133]

[0134]

[0135]

[0136]

[0137]

[0138] wherein C 11 , C 12 , C 13 , C 33, C 44 , C 55 and C 66 are initial elastic parameters, and ε, γ and δ are anisotropic parameters, respectively, v p0 is the longitudinal wave velocity of the medium along the symmetry axis direction, v s0 is the transverse wave velocity of the medium along the symmetry axis direction, and p is the medium density field.

[0139] Figure 3 FIG. 3 shows a flowchart of an anisotropic elastic wave decoupling method according to another embodiment of the present specification. The anisotropic wave field decoupling method process is described in this figure, but more or fewer operation steps can be included based on conventional or non-inventive labor. Specifically, as shown in FIG. 3, the method can include: Figure 3

[0140] S331, determining at least one elastic parameter to be separated from the initial elastic parameter set, the elastic parameter to be separated including an initial elastic parameter determined based on the transverse wave velocity and the longitudinal wave velocity;

[0141] S332, determining a corresponding transverse-longitudinal wave velocity separation step for each of the at least one elastic parameter to be separated;

[0142] S333, performing transverse-longitudinal wave velocity separation processing on each of the elastic parameters to be separated based on the transverse-longitudinal wave velocity separation step, to obtain a target longitudinal wave elastic parameter set and a target transverse wave elastic parameter set.

[0143] With the embodiments of the present specification, a transverse-longitudinal wave velocity separation step corresponding to each of the at least one elastic parameter to be separated determined from the initial elastic parameter set is determined. Then, the corresponding elastic parameter to be separated is subjected to transverse-longitudinal wave velocity separation operation by using the transverse-longitudinal wave velocity separation step, and then a target longitudinal wave elastic parameter set and a target transverse wave elastic parameter set are obtained. The transverse-longitudinal wave velocity separation is performed for each of the initial elastic parameters included in the initial elastic parameter set, and then the transverse-longitudinal wave field separation is performed for the wave field to be decomposed. Compared with the separated transverse-longitudinal wave field obtained without transverse-longitudinal wave velocity separation, the transverse-longitudinal wave interference of the image corresponding to the constant longitudinal wave field obtained by the embodiments of the present specification is smaller.

[0144] According to another embodiment of the present specification, the elastic parameter to be separated is an initial elastic parameter corresponding to a function composed of the transverse wave velocity and the longitudinal wave velocity. For example, the elastic parameter to be separated can be the initial elastic parameter C 12 or C 13 included in formula (3) or formula (4), respectively.

[0145] ​For each initial elastic parameter in the initial set of elastic parameters, identification is performed respectively as to whether the shear wave velocity and the longitudinal wave velocity are included simultaneously, and in a case where it is determined that the target initial elastic parameter is a function constituted by the shear wave velocity and the longitudinal wave velocity, the target initial elastic parameter is taken as the elastic parameter to be separated.

[0146] For each elastic parameter to be separated, a corresponding shear wave and longitudinal wave velocity separation step is determined according to a multi-wave velocity constitution type corresponding to the elastic parameter to be separated. The multi-wave velocity constitution type represents a manner in which the shear wave velocity and the longitudinal wave velocity constitute a function, such as addition, multiplication, and division, and the like. For each type, a corresponding shear wave and longitudinal wave velocity separation step can be configured in advance. For example, for an addition type, as shown in the above formula (3), the shear wave and longitudinal wave velocity separation step can be that the function of the elastic parameter to be separated is arranged into a shear wave velocity function and a longitudinal wave velocity function directly with respect to the shear wave velocity and the longitudinal wave velocity. Further, the shear wave velocity function is taken as a target shear wave elastic parameter corresponding to the elastic parameter to be separated, and the longitudinal wave velocity is taken as a target longitudinal wave parameter corresponding to the elastic parameter to be separated.

[0147] According to another embodiment of the present specification, in a case where the elastic parameter to be separated is determined based on a square root operation, the shear wave and longitudinal wave velocity separation step may, for example, include square root processing with respect to the elastic parameter to be separated to obtain a separated elastic parameter, and shear wave and longitudinal wave velocity separation processing with respect to the separated elastic parameter to obtain a target longitudinal wave elastic parameter and a target shear wave elastic parameter, for obtaining a target set of longitudinal wave elastic parameters and a target set of shear wave elastic parameters.

[0148] The elastic parameter to be separated determined based on the square root operation is, for example, as shown in the above formula (4).

[0149] According to another embodiment of the present specification, based on the shear wave and longitudinal wave velocity separation step, the shear wave and longitudinal wave velocity separation processing is performed with respect to each elastic parameter to be separated, and the determined target set of longitudinal wave elastic parameters can include the following formulas (8) to (13).

[0150]

[0151]

[0152]

[0153]

[0154]

[0155]

[0156] wherein, and are target longitudinal wave elastic parameters, and ε and δ are anisotropic parameters, respectively, v p0 is the longitudinal wave velocity of the medium along the symmetry axis direction, v s0 is the transverse wave velocity of the medium along the symmetry axis direction, and p is the medium density field.

[0157] According to another embodiment of the present specification, based on the transverse-longitudinal wave velocity separation step, the transverse-longitudinal wave velocity separation processing is performed on each elastic parameter to be separated, and the target transverse wave elastic parameter set can include the following formulas (14)-(18).

[0158]

[0159]

[0160]

[0161]

[0162]

[0163]

[0164] wherein, and are target transverse wave elastic parameters, and γ and δ are anisotropic parameters, respectively, v p0 is the longitudinal wave velocity of the medium along the symmetry axis direction, v s0 is the transverse wave velocity of the medium along the symmetry axis direction, and p is the medium density field.

[0165] As can be seen from the above formulas, for the elastic parameter C 12 to be separated, which is of the addition type, the transverse wave velocity and the longitudinal wave velocity are arranged, and the function of the elastic parameter to be separated is arranged into the target longitudinal wave elastic parameter and the target transverse wave elastic parameter Similarly, for the elastic parameter C 13 to be separated, which is of the square root type, the square root processing is performed on the elastic parameter to be separated to obtain the separated elastic parameter. The transverse-longitudinal wave velocity separation processing is performed on the separated elastic parameter to obtain the target longitudinal wave elastic parameter and the target transverse wave elastic parameter

[0166] Figure 4 The principle diagram of an anisotropic elastic wave decoupling method according to an embodiment of the present specification is shown.

[0167] According to another embodiment of the present specification, the target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set are substituted into the anisotropic model to process the wave field to be decomposed to obtain the target longitudinal wave matrix and the target transverse wave matrix, which includes: substituting the target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set into the anisotropic model to obtain a target anisotropic equation group; and using an interlaced grid high-order finite difference method to solve the target anisotropic equation group to obtain the target longitudinal wave matrix and the target transverse wave matrix.

[0168] The anisotropic model may be, for example, a first-order velocity-stress equation and an equation solving algorithm, which may be, for example, an interlaced grid high-order finite difference method.

[0169] It should be noted that the first-order velocity-stress equation of anisotropy is shown in the following formula (20).

[0170]

[0171] wherein v=(v x , v y , v z ) T is the particle vibration velocity component, τ=(τ xy , τ yy , τ zz , τ yz , τ xz , τ xy ) is the stress component, and T is the matrix transposition operation, is the particle vibration velocity v=(v x , v y , v z ) T the first-order derivative with respect to time, is the first-order derivative with respect to time of the stress τ=(τ xy , τ yy , τ zz , τ yz , τ xz , τ xy ), v x , v y , v z are the components of the velocity in the x, y, z directions, τ xy , τ yy , τ zz are three normal stress components, τ yz , τ xz , τ xy are three shear stress components, ρ is the medium density field, L is the partial derivative operator matrix, C is the elastic parameter matrix, and the partial derivative operator matrix and the elastic parameter matrix are respectively:

[0172]

[0173] wherein, are the derivatives in x, y, z directions, respectively, C 11 , C 12 , C 13 , C 33 , C 44 , C 55 , and C 66 are initial elastic parameters, respectively.

[0174] The obtained target longitudinal wave elastic parameter set only related to longitudinal wave velocity and target transverse wave elastic parameter set only related to transverse wave velocity are substituted into the first-order velocity stress equation to obtain a target anisotropic equation set. The target anisotropic equation set includes a longitudinal wave wave field equation composed of longitudinal wave stress and longitudinal wave particle vibration velocity and a transverse wave wave field equation composed of transverse wave stress and transverse wave particle vibration velocity.

[0175] Specifically, the longitudinal wave wave field equation can be shown as the following formula (21) for example, and the transverse wave wave field equation can be shown as the following formula (22) for example.

[0176]

[0177]

[0178] wherein, is a longitudinal wave stress component, is a transverse wave stress component, and are a horizontal component and a vertical component of longitudinal wave particle vibration velocity; and are a horizontal component and a vertical component of transverse wave particle vibration velocity, ε, γ, and δ are anisotropic parameters, respectively, v p0 is a longitudinal wave velocity of the medium along the symmetry axis direction, v s0 is a transverse wave velocity of the medium along the symmetry axis direction, and ρ is a medium density field, x, y, and z are three corresponding directions, respectively.

[0179] It is to be noted that the total wave field is obtained by adding the longitudinal wave wave field and the transverse wave wave field, and the longitudinal wave field equation is obtained by the following formula (23).

[0180]

[0181] wherein, τ = (τ xx , τ yy , τ zz , τ yz , τ xz , τ xy are stress components, v = (vx , v y , v z ) T is the particle vibration velocity component of the longitudinal wave wave field, and is the particle vibration velocity component of the transverse wave wave field, and and is the particle vibration velocity component of the longitudinal wave wave field, and is the particle vibration velocity component of the transverse wave wave field, and and is the stress component of the longitudinal wave wave field, and is the stress component of the transverse wave wave field.

[0182] Further, the staggered grid high-order finite difference method is used to solve the above formula (21) and formula (22), and the target longitudinal wave matrix and the target transverse wave matrix are obtained.

[0183] As shown in Figure 4 , when the wave field to be decomposed 401 is received, the corresponding anisotropic model 410 is determined, and the conversion as shown in formula (1) is performed on the Thomson parameter set included in the anisotropic model 410, to obtain the initial elastic parameter set 420 as shown in formula (2)-formula (7).

[0184] The initial elastic parameter set 420 is subjected to transverse-longitudinal wave velocity separation processing, to obtain the target longitudinal wave elastic parameter set 431 as shown in formula (8)-formula (13), and the target transverse wave elastic parameter set 432 as shown in formula (14)-formula (19).

[0185] The target longitudinal wave elastic parameter set 431 and the target transverse wave elastic parameter set 432 are substituted into the anisotropic model 410 for processing the wave field to be decomposed 401, to obtain the target anisotropic equation group 440 as shown in formula (21) and formula (22).

[0186] The staggered grid high-order finite difference method is used to solve the target anisotropic equation group 440, to obtain the target longitudinal wave matrix 451 and the target transverse wave matrix 452, which can be used for imaging.

[0187] Figure 5A Fig. 1 shows a staggered grid finite difference format diagram according to an embodiment of the present specification; Figure 5B Fig. 2 shows a wave field snapshot diagram corresponding to the horizontal component of the full wave field particle vibration velocity field according to an embodiment of the present specification; Figure 5C Fig. 3 shows a wave field snapshot diagram corresponding to the horizontal component of the longitudinal wave particle vibration velocity field according to an embodiment of the present specification; Figure 5D Fig. 4 shows a wave field snapshot diagram corresponding to the horizontal component of the transverse wave particle vibration velocity field according to an embodiment of the present specification; Figure 5E Fig. 5 shows a wave field snapshot diagram corresponding to the vertical component of the full wave field particle vibration velocity field according to an embodiment of the present specification; and Figure 5FFig. 1 is a schematic diagram of a wave field snapshot corresponding to a vertical component of a particle vibration velocity field of a longitudinal wave according to an embodiment of the present specification; Figure 5G Fig. 2 is a schematic diagram of a wave field snapshot corresponding to a vertical component of a particle vibration velocity field of a transverse wave according to an embodiment of the present specification; Figure 5H Fig. 3 is a schematic diagram of a single-channel waveform curve corresponding to a horizontal component of a particle vibration velocity field of a longitudinal wave according to an embodiment of the present specification; Figure 5I Fig. 4 is a schematic diagram of a single-channel waveform curve corresponding to a horizontal component of a particle vibration velocity field of a transverse wave according to an embodiment of the present specification; Figure 5J Fig. 5 is a schematic diagram of a single-channel waveform curve corresponding to a vertical component of a particle vibration velocity field of a longitudinal wave according to an embodiment of the present specification; Figure 5K Fig. 6 is a schematic diagram of a single-channel waveform curve corresponding to a vertical component of a particle vibration velocity field of a transverse wave according to an embodiment of the present specification; Figure 5L Fig. 7 is a schematic diagram of a model parameter according to an embodiment of the present specification; Figure 5M Fig. 8 is a schematic diagram of a reverse-time migration imaging result according to an embodiment of the present specification.

[0188] Upon receiving a wave field to be decomposed, an anisotropic model is constructed based on the wave field to be decomposed. The model includes a staggered grid finite difference format as shown in Figure 5A The anisotropic model has a size of 4000m x 4000m, and the velocities of the longitudinal and transverse waves are 3500m / s and 2000m / s respectively, and the density is 1g / cm 3 The anisotropic parameters ε = 0.15 and δ = 0.1. Numerical simulation is performed by using a staggered grid spatial tenth-order difference and a time second-order difference format, the spatial sampling interval is 10m in both the horizontal and vertical directions, the time sampling interval is 1ms, the source is an explosion source, and the source wavelet is a 20Hz Ricker wavelet. The source is located at the center of the model.

[0189] It should be noted that the staggered grid finite difference format as shown in Figure 5A is defined according to the following formula (24).

[0190]

[0191] where τ xx , τ yy , τ zz , τ yz , τ xz , τ xy are stress components, v x , v y , v z are velocity components, i, j, k are spatial integer grid points, and h is a spatial half grid point.

[0192] Based on the above, the Figure 5AThe defined difference format obtains the staggered grid high-order finite difference form as shown in the following formula (25).

[0193]

[0194] Wherein, τ xx , τ yy , τ zz , τ xy , τ yz , τ xz are stress components, v x , v y , v z are velocity components, Δt is a time sampling interval, n is an integer time point, and are half time points, Δx, Δy, Δz are sampling intervals along x, y, z directions respectively, ρ is medium density, C 11 , C 12 , C 13 , C 33 , C 44 , C 66 are initial elastic parameters, F x , F y , F z are forward staggered grid difference formats in x, y, z directions, B x , B y , B z are backward staggered grid difference coefficients in x, y, z directions, and the specific forms are as follows:

[0195]

[0196]

[0197]

[0198]

[0199]

[0200]

[0201] Wherein, is a 2M-order staggered grid difference coefficient, U i,j,k is a wave field value at a grid point (i, j, k), U i+m,j,k , U i,j+m,k , U i,j,k+m are wave field values at grid points (i+m, j, k), (i, j+m, k), (i, j, k+m) respectively, and Δx, Δy, Δz are sampling intervals in x, y, z directions.

[0202] For the to-be-decomposed wave field obtained based on the source wavelet, the method of the embodiments of the present specification is used for transverse and longitudinal wave field separation to obtain a wave field snapshot corresponding to a 0.4s full wave field particle vibration velocity field. Specifically, the wave field snapshot corresponding to the 0.4s full wave field particle vibration velocity field includes a wave field snapshot corresponding to a horizontal component (VX) of the full wave field particle vibration velocity field as shown in FIG. 8A; a wave field snapshot corresponding to a horizontal component (PVX) of the longitudinal wave particle vibration velocity field as shown in FIG. 8B; a wave field snapshot corresponding to a horizontal component (SVX) of the transverse wave particle vibration velocity field as shown in FIG. 8C; a wave field snapshot corresponding to a vertical component (VZ) of the full wave field particle vibration velocity field as shown in FIG. 8D; a wave field snapshot corresponding to a vertical component (PVZ) of the longitudinal wave particle vibration velocity field as shown in FIG. 8E; and a wave field snapshot corresponding to a vertical component (SVZ) of the transverse wave particle vibration velocity field as shown in FIG. 8F. The wave field snapshots shown in FIGS. 8A-8F can show that the results obtained by using the embodiments of the present specification for numerical simulation have a certain crosstalk, but the residual amplitude is relatively weak, and the overall transverse and longitudinal wave field separation effect is good. Figure 5B Figure 5C Figure 5D Figure 5E Figure 5F Figure 5G As can be seen from the wave field snapshots shown in FIGS. 8A-8F, although the results obtained by using the embodiments of the present specification for numerical simulation have a certain crosstalk, the residual amplitude is relatively weak, and the overall transverse and longitudinal wave field separation effect is good. Figure 5B - Figure 5G

[0203] In order to further verify the effectiveness of the transverse and longitudinal wave field separation of the method of the embodiments of the present specification, the decomposed longitudinal and transverse single-channel waveforms at the grid point (x=2400m, y=2400m) are extracted. Specifically, a single-channel waveform curve corresponding to the horizontal component of the longitudinal wave particle vibration velocity field as shown in FIG. 9A; a single-channel waveform curve corresponding to the horizontal component of the transverse wave particle vibration velocity field as shown in FIG. 9B; a single-channel waveform curve corresponding to the vertical component of the longitudinal wave particle vibration velocity field as shown in FIG. 9C; and a single-channel waveform curve corresponding to the vertical component of the transverse wave particle vibration velocity field as shown in FIG. 9D. As can be seen from the single-channel waveform curves shown in FIGS. 9A-9D, the residual amplitude relative relationship (the residual part is the part in the circle) can be directly observed. Among them, for the x component, the maximum relative amplitude of the residual transverse wave in the longitudinal wave is 4.2%, and the maximum relative amplitude of the residual longitudinal wave in the transverse wave is 1.7%, and for the z component, they are 3.2% and 5.1% respectively. It can be seen that the residual amplitude is generally small, and since the stacking is performed in the reverse time migration imaging process, the interference of random noise can be effectively reduced, so that the method of the embodiments of the present specification applied to the reverse time migration imaging will not produce obvious imaging artifacts. Figure 5H Figure 5I Figure 5J Figure 5K Figure 5H - Figure 5K

[0204] ​​​​​​​​​​​Further, an anisotropy model (VTI-BGP model) provided by BGP is used for imaging to verify the correctness of the embodiments of the present specification, the model has a horizontal width of 40 km and a vertical depth of 6.7 km, the grid spacing is 25 m horizontally and 10 m vertically, 800 receivers and shots are evenly distributed on the surface of the model with an interval of 50 m, and a Ricker wavelet (Ricker wavelet) with a main frequency of 20 Hz is selected as the source wavelet. Figure 5L The P and S wave velocities and anisotropy parameters of the VTI-BGP model are shown, including a P wave velocity parameter field 501, an S wave velocity parameter field 502, an ε parameter field 503, and a δ parameter field 504.

[0205] Based on the target P wave elastic parameter set and the target S wave elastic parameter set obtained based on the VTI-BGP model based on the embodiments of the present specification, P and S wave field separation is performed to obtain four reverse time migration imaging results as shown in Figure 5M The four imaging results include a source P wave receiver P wave image 511 (PP imaging result), a source P wave receiver S wave image 512 (PS imaging result), a source S wave receiver P wave image 513 (SP imaging result), and a source S wave receiver S wave image 514 (SS imaging result). As can be seen from the imaging results shown in Figure 5M As can be seen from the imaging results shown in

[0206] Figure 6 Fig. 1 shows a structural schematic diagram of an anisotropic elastic wave decoupling device according to an embodiment of the present specification. As shown in Figure 6 Fig. 1 shows a structural schematic diagram of an anisotropic elastic wave decoupling device according to an embodiment of the present specification. As shown in

[0207] The first determination unit 610 is configured to determine a Thomson parameter set included in the anisotropy model according to the received wave field decomposition request to be decomposed, and the wave field decomposition request to be decomposed includes a wave field to be decomposed.

[0208] The conversion unit 620 is configured to convert the Thomson parameter set to obtain an initial elastic parameter set.

[0209] The separation unit 630 is configured to perform P and S wave velocity separation processing on the initial elastic parameter set to obtain a target P wave elastic parameter set and a target S wave elastic parameter set.

[0210] The processing unit 640 is configured to substitute the target P wave elastic parameter set and the target S wave elastic parameter set into the anisotropy model, and process the wave field to be decomposed to obtain a target P wave matrix and a target S wave matrix.

[0211] Since the principle of solving the problem of the above device is similar to that of the above method, the implementation of the above device can be referred to the implementation of the above method, and the repeated parts will not be described again.

[0212] As Figure 7 FIG. 7 shows a diagram of a computer device 702 in accordance with an embodiment of the present specification. The apparatus in the present specification can be the computer device in the present embodiment, which executes the method of the present specification. The computer device 702 can 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 can also include any storage resources 706 for storing any kind of information, such as code, settings, data, etc. Without limitation, for example, the storage resources 706 can include any one or combination of the following: any type of RAM, any type of ROM, flash memory devices, hard disks, optical disks, etc. More generally, any storage resource can store information using any technology. Further, any storage resource can provide volatile or non-volatile retention of information. Further, any storage resource can represent a fixed or removable component of the computer device 702. In one case, the computer device 702 can perform any operation of the associated instructions when the processing device 704 executes the associated instructions stored in any storage resource or combination of storage resources. The computer device 702 also includes one or more drive mechanisms 708 for interacting with any storage resources, such as a hard drive mechanism, an optical drive mechanism, etc.

[0213] The computer device 702 can also include an input / output module 710 (I / O) for receiving various inputs (via input devices 712) and for providing various outputs (via output devices 714). One particular output mechanism can include a presentation device 716 and an associated graphical user interface (GUI) 718. In other embodiments, the input / output module 710 (I / O), the input devices 712, and the output devices 714 can also not be included, just as a computer device in a network. The computer device 702 can 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 above-described components together.

[0214] The communication links 722 can be implemented in any manner, for example, through a local area network, a wide area network (e.g., the Internet), a point-to-point connection, etc., or any combination thereof. The communication links 722 can include any combination of hardwired links, wireless links, routers, gateway functionality, name servers, etc., governed by any protocol or combination of protocols.

[0215] The present specification also provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the above method.

[0216] The embodiment of the present specification further provides a computer program product, the computer program product comprising a computer program, the computer program being executed by a processor to implement the method described above.

[0217] Those skilled in the art will understand that the embodiment of the present specification can be provided as a method, a system, or a computer program product. Therefore, the present specification can take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present specification can take the form of a computer program product implemented on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROMs, optical storage devices, etc.) containing computer-usable program code.

[0218] The present specification is described with reference to flowcharts and / or block diagrams of the method, device (system), and computer program product according to the embodiment of the present specification. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of flows and / or blocks in the flowcharts 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, a special-purpose computer, an embedded processor, or other programmable data processing apparatus to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing apparatus produce a device implemented in the flowcharts and / or block diagrams. Figure 1 one or more flows and / or blocks Figure 1 an apparatus that performs the functions specified in the flow(s) or block(s).

[0219] These computer program instructions can also be stored in a computer-readable memory that can direct the computer or other programmable data processing apparatus to work in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including instruction apparatus, which implements the flowcharts and / or block diagrams. Figure 1 one or more flows and / or blocks Figure 1 an apparatus that performs the functions specified in the flow(s) or block(s).

[0220] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus, so that a series of operation steps are performed on the computer or other programmable data processing apparatus to produce a computer-implemented process, so that the instructions executed on the computer or other programmable data processing apparatus provide a process for implementing the flowcharts and / or block diagrams. Figure 1 one or more flows and / or blocks Figure 1 an apparatus that performs the functions specified in the flow(s) or block(s).

[0221] The above specific embodiments further illustrate the purpose, technical solutions and beneficial effects of the present specification, and it should be understood that the above are only specific embodiments of the present specification and are not used to limit the protection scope of the present specification, and any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present specification shall be included in the protection scope of the present specification.

Claims

1. An anisotropic elastic wave decoupling method applied to a wavefield caused by an earthquake, characterized in that, Comprising: According to the received wave field decomposition request to be decomposed, determine the thomson parameter set included in the anisotropy model, the wave field decomposition request to be decomposed includes the wave field to be decomposed; For the thomson parameter set, conversion is carried out to obtain an initial elastic parameter set; For the initial elastic parameter set, transverse wave and longitudinal wave velocity separation processing is carried out to obtain a target longitudinal wave elastic parameter set and a target transverse wave elastic parameter set; And The target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set are substituted into the anisotropy model, and the wave field to be decomposed is processed to obtain a target longitudinal wave matrix and a target transverse wave matrix; The target longitudinal wave elastic parameter set comprises: Wherein, the , the , the , the , the , the , and the are target longitudinal wave elastic parameters, the and the are anisotropic parameters, the is a longitudinal wave velocity of the medium along the symmetry axis direction, the is a transverse wave velocity of the medium along the symmetry axis direction, and the is a medium density field; The target transverse wave elastic parameter set comprises: wherein the , the , the , the , the , the , and the are target shear wave elastic parameters, the , and the are anisotropic parameters, the is a P-wave velocity of the medium along the symmetry axis direction, the is a S-wave velocity of the medium along the symmetry axis direction, and the is a medium density field.

2. The method of claim 1, wherein, The initial elastic parameter set comprises: wherein the , the , the , the , the , the , and the are initial elastic parameters, the , the , and the are anisotropic parameters, the is the longitudinal wave velocity of the medium along the symmetry axis direction, the is the transverse wave velocity of the medium along the symmetry axis direction, and the is the medium density field.

3. The method of claim 1, wherein, The initial elastic parameter set, transverse wave and longitudinal wave velocity separation processing is carried out to obtain a target longitudinal wave elastic parameter set and a target transverse wave elastic parameter set, comprising: From the initial elastic parameter set, at least one elastic parameter to be separated is determined, the elastic parameter to be separated includes an initial elastic parameter determined based on transverse wave velocity and longitudinal wave velocity; For each elastic parameter to be separated in the at least one elastic parameter to be separated, a corresponding transverse wave and longitudinal wave velocity separation step is determined; and Based on the transverse wave and longitudinal wave velocity separation step, the transverse wave and longitudinal wave velocity separation processing is carried out on each elastic parameter to be separated to obtain the target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set.

4. The method of claim 3, wherein, In the case of the elastic parameter to be separated being determined based on square root operation, the transverse wave and longitudinal wave velocity separation step comprises: Square root processing is carried out on the elastic parameter to be separated to obtain a separated elastic parameter; and The target longitudinal wave elastic parameter and the target transverse wave elastic parameter are obtained by carrying out transverse wave and longitudinal wave velocity separation processing on the separated elastic parameter, which are used to obtain the target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set.

5. The method of claim 1, wherein, The target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set are substituted into the anisotropy model, and the wave field to be decomposed is processed to obtain a target longitudinal wave matrix and a target transverse wave matrix, comprising: The target longitudinal wave elastic parameter set and the target transverse wave elastic parameter set are substituted into the anisotropy model to obtain a target anisotropy equation set; and The target anisotropy equation set is solved by using staggered grid high-order finite difference method to obtain the target longitudinal wave matrix and the target transverse wave matrix.

6. An anisotropic elastic wave decoupling device applied in a wave field caused by an earthquake, characterized in that, Comprising: A first determination unit is configured to determine a thomson parameter set included in an anisotropy model according to a received wave field decomposition request to be decomposed, the wave field decomposition request to be decomposed including a wave field to be decomposed; A conversion unit is configured to convert the thomson parameter set to obtain an initial elastic parameter set; A separation unit is configured to carry out transverse wave and longitudinal wave velocity separation processing on the initial elastic parameter set to obtain a target longitudinal wave elastic parameter set and a target transverse wave elastic parameter set; And a processing unit, configured to substitute the target set of P-wave elastic parameters and the target set of S-wave elastic parameters into the anisotropic model, and process the wave field to be decomposed to obtain a target P-wave matrix and a target S-wave matrix; the target set of P-wave elastic parameters comprises: wherein the , the , the , the , the , the , and the are target P-wave elastic parameters, the and the are anisotropic parameters, the is the P-wave velocity of the medium along the symmetry axis direction, the is the S-wave velocity of the medium along the symmetry axis direction, and the is the medium density field; the target set of S-wave elastic parameters comprises: wherein the , the , the , the , the , the , and the are target transverse wave elastic parameters, the and the are anisotropic parameters, the is a longitudinal wave velocity of the medium along the symmetry axis direction, the is a transverse wave velocity of the medium along the symmetry axis direction, and the is a medium density field.

7. A computer device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, the processor executes the computer program to implement 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, and the computer program is run by the processor to execute the method of any one of claims 1-5.

Images