Anisotropy inversion initial model establishment method and device, electronic equipment and medium

Abstract

The invention discloses an anisotropic inversion initial model establishment method and device, electronic equipment and a medium. The method comprises the following steps: determining longitudinal wave impedance and a plurality of rock physical characteristics for anisotropic inversion; respectively establishing a plurality of expressions between the longitudinal wave impedance and the plurality of rock physical characteristics; obtaining a longitudinal wave impedance value according to the seismic data; respectively substituting the value of the longitudinal wave impedance into a plurality of expressions, and calculating to obtain corresponding rock physical characteristics; and establishing an anisotropic inversion initial model according to the longitudinal wave impedance and the plurality of rock physical characteristics. According to the method, starting from the initial model of the optimization theory, the rock physics constraint and the post-stack inversion technology are closely combined, the initial model suitable for the pre-stack anisotropy inversion is established, the rock physics relation between the initial models is guaranteed, and the problem that the pre-stack anisotropy inversion is caught in local optimum is solved to a certain extent.

Description

Technical Field

[0001] The present invention relates to the field of petroleum geophysical exploration, and more specifically, to a method, device, electronic equipment and medium for establishing an anisotropic inversion initial model. Background Art

[0002] Prestack seismic inversion uses the variation of amplitude with shot offset to obtain the lithology and fluid of the underground medium. These amplitude variations can be simulated using various approximations of the Zoeppritz equation or the Zoepptriz equation, each with specific assumptions and limitations. These geophysicists use different geophysical parameters to approximate the Zoeppritz equation. Typical prediction parameters are acoustic impedance and VP / VS ratio, acoustic impedance and elastic impedance, or Lamés parameters, the product of lambda times density and shear modulus times density. The above studies and applications only consider the isotropic case and ignore the effect of anisotropy on the AVO response. But in shales, which have intrinsic elastic anisotropy and are much stronger than other sedimentary rocks, the seismic velocity in this medium is angle-dependent, and the behavior of elastic waves is much more complicated than in the isotropic case. Wright (1987), Kim et al. (1993) and Thomsen (1993) analyzed the effect of anisotropy on seismic amplitude in VTI media. Thomsen (1993) showed that his weak anisotropy parameter can be used to add transverse isotropy terms to the Aki-Richards equation. Ruger (2002) gave a modified form of the P-wave reflection coefficient in VTI media. At present, reservoir anisotropy inversion modeling is usually obtained by interpolation of well logging curves. This technical solution does not fully utilize the rock physical constraints of the reservoir and aggravates the multi-solution nature of the inversion.

[0003] At present, there is still a need to develop a method for establishing anisotropic inversion initial model constrained by rock physics.

[0004] The information disclosed in the background technology section of the present invention is only intended to deepen the understanding of the general background technology of the present invention, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to those skilled in the art. Summary of the invention

[0005] The present invention proposes a method, device, electronic equipment and medium for establishing an anisotropic inversion initial model. Starting from the initial model of optimization theory, the rock physics constraints are closely combined with the post-stack inversion technology to establish an initial model suitable for pre-stack anisotropic inversion, ensure the rock physics relationship between the initial models, and overcome the problem of pre-stack anisotropic inversion falling into local optimum to a certain extent.

[0006] In a first aspect, an embodiment of the present disclosure provides a method for establishing an anisotropic inversion initial model, comprising:

[0007] Determine P-wave impedance and multiple rock physical characteristics for anisotropic inversion;

[0008] Establishing multiple expressions between the longitudinal wave impedance and multiple rock physical characteristics respectively;

[0009] Obtaining the value of the longitudinal wave impedance according to seismic data;

[0010] Substituting the value of the longitudinal wave impedance into a plurality of expressions respectively, and calculating the corresponding rock physical characteristics;

[0011] An anisotropic inversion initial model is established based on the longitudinal wave impedance and the multiple rock physical characteristics.

[0012] As a specific implementation of the embodiment of the present disclosure, the rock physical characteristics include shear wave impedance, density, and anisotropy parameters.

[0013] As a specific implementation of the embodiment of the present disclosure, the expressions are:

[0014] Zs=aZp+b

[0015] ρ=cZp+d

[0016] ε=eZp+f

[0017] δ=gZp+h

[0018] Among them, Zp is the P-wave impedance, Zs is the S-wave impedance, ρ is the density, ε and δ are anisotropy parameters that characterize the degree of anisotropy of the underground medium; a and b are the rock physical characteristic parameters of S-wave impedance and P-wave impedance; c and d are the rock physical characteristic parameters of density and P-wave impedance; e and f are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance; g and h are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance.

[0019] As a specific implementation manner of the embodiment of the present disclosure, the longitudinal wave impedance is obtained through post-stack inversion based on the seismic data.

[0020] As a specific implementation of the embodiment of the present disclosure, it also includes:

[0021] After the corresponding rock physical characteristics are calculated, the rock physical characteristics are smoothed and filtered.

[0022] As a specific implementation of the embodiment of the present disclosure, the smoothing filter is a median filter.

[0023] As a specific implementation of the embodiment of the present disclosure, the time window of the median filter is between 100ms and 200ms.

[0024] In a second aspect, the present disclosure also provides an anisotropic inversion initial model establishment device, comprising:

[0025] Parameter preparation module, which determines the P-wave impedance and multiple rock physical characteristics for anisotropic inversion;

[0026] An expression establishment module, for establishing a plurality of expressions between the longitudinal wave impedance and a plurality of rock physical characteristics respectively;

[0027] A first calculation module obtains the value of the longitudinal wave impedance according to seismic data;

[0028] A second calculation module substitutes the value of the longitudinal wave impedance into a plurality of expressions to calculate the corresponding rock physical characteristics;

[0029] A modeling module is used to establish an anisotropic inversion initial model based on the longitudinal wave impedance and the multiple rock physical characteristics.

[0030] As a specific implementation of the embodiment of the present disclosure, the rock physical characteristics include shear wave impedance, density, and anisotropy parameters.

[0031] As a specific implementation of the embodiment of the present disclosure, the expressions are:

[0032] Zs=aZp+b

[0033] ρ=cZp+d

[0034] ε=eZp+f

[0035] δ=gZp+h

[0036] Among them, Zp is the P-wave impedance, Zs is the S-wave impedance, ρ is the density, ε and δ are anisotropy parameters that characterize the degree of anisotropy of the underground medium; a and b are the rock physical characteristic parameters of S-wave impedance and P-wave impedance; c and d are the rock physical characteristic parameters of density and P-wave impedance; e and f are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance; g and h are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance.

[0037] As a specific implementation manner of the embodiment of the present disclosure, the longitudinal wave impedance is obtained through post-stack inversion based on the seismic data.

[0038] As a specific implementation of the embodiment of the present disclosure, it also includes:

[0039] After the corresponding rock physical characteristics are calculated, the rock physical characteristics are smoothed and filtered.

[0040] As a specific implementation of the embodiment of the present disclosure, the smoothing filter is a median filter.

[0041] As a specific implementation of the embodiment of the present disclosure, the time window of the median filter is between 100ms and 200ms.

[0042] In a third aspect, an embodiment of the present disclosure further provides an electronic device, the electronic device comprising:

[0043] A memory storing executable instructions;

[0044] A processor runs the executable instructions in the memory to implement the anisotropic inversion initial model establishment method.

[0045] In a fourth aspect, an embodiment of the present disclosure further provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the anisotropic inversion initial model establishing method is implemented.

[0046] The methods and apparatus of the present invention have other features and advantages that will be apparent from, or will be described in detail in, the accompanying drawings and subsequent detailed descriptions incorporated herein, which together serve to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention in conjunction with the accompanying drawings, wherein like reference numerals generally represent like components throughout the exemplary embodiments of the present invention.

[0048] Figure 1 A flow chart showing the steps of a method for establishing an anisotropic inversion initial model according to an embodiment of the present invention.

[0049] Figure 2 A schematic diagram of an initial model of longitudinal wave impedance according to an embodiment of the present invention is shown.

[0050] Figure 3 A schematic diagram of a shear wave impedance initial model according to an embodiment of the present invention is shown.

[0051] Figure 4 A schematic diagram of a density initialization model according to an embodiment of the present invention is shown.

[0052] Figure 5 A schematic diagram of an initial model of anisotropy parameter ε according to an embodiment of the present invention is shown.

[0053] Figure 6A schematic diagram of an initial model of anisotropy parameter δ according to an embodiment of the present invention is shown.

[0054] Figure 7 A block diagram of an anisotropic inversion initial model building device according to an embodiment of the present invention is shown.

[0055] Description of reference numerals:

[0056] 201, parameter preparation module; 202, expression establishment module; 203, first calculation module; 204, second calculation module; 205, modeling module. DETAILED DESCRIPTION

[0057] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.

[0058] To facilitate understanding of the solutions and effects of the embodiments of the present invention, four specific application examples are given below. Those skilled in the art should understand that the examples are only for facilitating understanding of the present invention, and any specific details thereof are not intended to limit the present invention in any way.

[0059] Example 1

[0060] Figure 1 A flow chart showing the steps of a method for establishing an anisotropic inversion initial model according to an embodiment of the present invention.

[0061] like Figure 1 As shown, the method for establishing the initial model of anisotropic inversion includes: step 101, determining the P-wave impedance and multiple rock physical characteristics used for anisotropic inversion; step 102, respectively establishing multiple expressions between the P-wave impedance and the multiple rock physical characteristics; step 103, obtaining the value of the P-wave impedance according to seismic data; step 104, substituting the value of the P-wave impedance into multiple expressions respectively, and calculating the corresponding rock physical characteristics; step 105, establishing the initial model of anisotropic inversion according to the P-wave impedance and the multiple rock physical characteristics.

[0062] In one example, the rock physical characteristics include shear impedance, density, and anisotropy parameters.

[0063] In an example, the expressions are:

[0064] Zs=aZp+b

[0065] ρ=cZp+d

[0066] ε=eZp+f

[0067] δ=gZp+h

[0068] Among them, Zp is the P-wave impedance, Zs is the S-wave impedance, ρ is the density, ε and δ are anisotropy parameters that characterize the degree of anisotropy of the underground medium; a and b are the rock physical characteristic parameters of S-wave impedance and P-wave impedance; c and d are the rock physical characteristic parameters of density and P-wave impedance; e and f are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance; g and h are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance.

[0069] In one example, P-wave impedance is obtained from seismic data through post-stack inversion.

[0070] In one example, it also includes:

[0071] After the corresponding rock physical characteristics are calculated, the rock physical characteristics are smoothed and filtered.

[0072] In one example, the smoothing filter is a median filter.

[0073] In one example, the time window of the median filter is between 100 ms and 200 ms.

[0074] Specifically, the rock physical parameters of complex reservoirs are used to establish the rock physical characteristics of the inversion variables used for anisotropic inversion. The rock physical characteristics are expressed by the following formula:

[0075] Zs=aZp+b

[0076] ρ=cZp+d

[0077] ε=eZp+f

[0078] δ=gZp+h

[0079] Among them, the inversion variables used for anisotropic inversion include Zp, Zs, ρ, ε, δ, Zp is the P-wave impedance, Zs is the S-wave impedance, ρ is the density, ε and δ are the anisotropic parameters characterizing the anisotropy degree of the underground medium; a and b are the rock physical characteristic parameters of the S-wave impedance Zs and the P-wave impedance Zp; c and d are the rock physical characteristic parameters of the density ρ and the P-wave impedance Zp; e and f are the rock physical characteristic parameters of the anisotropic parameter ε and the P-wave impedance Zp; g and h are the rock physical characteristic parameters of the anisotropic parameter δ and the P-wave impedance Zp; the value range of the rock physical characteristic parameters is:

[0080] a∈[0.18,0.85],b∈[-0.71,0.42];

[0081] c∈[0.01,0.03],d∈[2.36,2.46];

[0082] e∈[-0.14, 0.51], f∈[-4.8, 1.7];

[0083] g∈[-0.16,0.42],h∈[-4.01,1.6];

[0084] The P-wave impedance Zp0 is obtained by post-stack inversion of seismic data. The inversion variables Zs0, ρ0, ε0, and δ0 used for anisotropic inversion are calculated based on the rock physical characteristics and rock physical characteristic parameters. The specific calculation method is as follows:

[0085] Zs0=aZp0+b

[0086] ρ0=cZp0+d

[0087] ε0=eZp0+f

[0088] δ0=gZp0+h

[0089] Post-stack inversion is the wave impedance inversion in the existing technology;

[0090] The inversion variables Zp0, Zs0, ρ0, ε0, and δ0 are smoothed and filtered. The smoothing filter is a median filter in the prior art, and the filtering time window is selected between 100ms and 200ms.

[0091] Example 2

[0092] The present invention also provides an anisotropic inversion initial model establishment device, comprising:

[0093] Parameter preparation module, which determines the P-wave impedance and multiple rock physical characteristics for anisotropic inversion;

[0094] An expression building module, for building multiple expressions between the longitudinal wave impedance and multiple rock physical characteristics;

[0095] The first calculation module obtains the value of longitudinal wave impedance according to the seismic data;

[0096] The second calculation module substitutes the value of the longitudinal wave impedance into multiple expressions to calculate the corresponding rock physical characteristics;

[0097] The modeling module establishes the initial anisotropic inversion model based on the longitudinal wave impedance and multiple rock physical characteristics.

[0098] In one example, the rock physical characteristics include shear impedance, density, and anisotropy parameters.

[0099] In an example, the expressions are:

[0100] Zs=aZp+b

[0101] ρ=cZp+d

[0102] ε=eZp+f

[0103] δ=gZp+h

[0104] Among them, Zp is the P-wave impedance, Zs is the S-wave impedance, ρ is the density, ε and δ are anisotropy parameters that characterize the degree of anisotropy of the underground medium; a and b are the rock physical characteristic parameters of S-wave impedance and P-wave impedance; c and d are the rock physical characteristic parameters of density and P-wave impedance; e and f are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance; g and h are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance.

[0105] In one example, P-wave impedance is obtained from seismic data through post-stack inversion.

[0106] In one example, it also includes:

[0107] After the corresponding rock physical characteristics are calculated, the rock physical characteristics are smoothed and filtered.

[0108] In one example, the smoothing filter is a median filter.

[0109] In one example, the time window of the median filter is between 100 ms and 200 ms.

[0110] Specifically, the rock physical parameters of complex reservoirs are used to establish the rock physical characteristics of the inversion variables used for anisotropic inversion. The rock physical characteristics are expressed by the following formula:

[0111] Zs=aZp+b

[0112] ρ=cZp+d

[0113] ε=eZp+f

[0114] δ=gZp+h

[0115] Among them, the inversion variables used for anisotropic inversion include Zp, Zs, ρ, ε, δ, Zp is the P-wave impedance, Zs is the S-wave impedance, ρ is the density, ε and δ are the anisotropic parameters characterizing the anisotropy degree of the underground medium; a and b are the rock physical characteristic parameters of the S-wave impedance Zs and the P-wave impedance Zp; c and d are the rock physical characteristic parameters of the density ρ and the P-wave impedance Zp; e and f are the rock physical characteristic parameters of the anisotropic parameter ε and the P-wave impedance Zp; g and h are the rock physical characteristic parameters of the anisotropic parameter δ and the P-wave impedance Zp; the value range of the rock physical characteristic parameters is:

[0116] a∈[0.18,0.85],b∈[-0.71,0.42];

[0117] c∈[0.01,0.03],d∈[2.36,2.46];

[0118] e∈[-0.14, 0.51], f∈[-4.8, 1.7];

[0119] g∈[-0.16,0.42],h∈[-4.01,1.6];

[0120] The P-wave impedance Zp0 is obtained by post-stack inversion of seismic data. The inversion variables Zs0, ρ0, ε0, and δ0 used for anisotropic inversion are calculated based on the rock physical characteristics and rock physical characteristic parameters. The specific calculation method is as follows:

[0121] Zs0=aZp0+b

[0122] ρ0=cZp0+d

[0123] ε0=eZp0+f

[0124] δ0=gZp0+h

[0125] Post-stack inversion is the wave impedance inversion in the existing technology;

[0126] The inversion variables Zp0, Zs0, ρ0, ε0, and δ0 are smoothed and filtered. The smoothing filter is a median filter in the prior art, and the filtering time window is selected between 100ms and 200ms.

[0127] Example 3

[0128] Aiming at the problem that pre-stack anisotropic inversion is prone to fall into local optimum, the present invention starts from the initial model of optimization theory, closely combines rock physics constraints with post-stack inversion technology, establishes an initial model suitable for pre-stack anisotropic inversion, and ensures the rock physics relationship between the initial models, in order to overcome the problem of pre-stack anisotropic inversion falling into local optimum to a certain extent. First, a linear relationship model is established using the rock physics relationship of complex reservoirs, and then inversion constraints are established using post-stack inversion, and the inversion constraints are filtered to obtain the initial model of anisotropic inversion.

[0129] Taking actual seismic data of a certain region as an example, an initial model suitable for pre-stack anisotropic inversion is established to further illustrate the effect of the present invention.

[0130] Figure 2 A schematic diagram of an initial model of longitudinal wave impedance according to an embodiment of the present invention is shown, which is obtained by post-stack wave impedance inversion in the prior art.

[0131] Figure 3A schematic diagram of the initial shear wave impedance model according to an embodiment of the present invention is shown, using the rock physics relationship Zs0 = aZp0 + b. Figure 2 The longitudinal wave impedance is transformed to obtain a=0.53, b=-0.004.

[0132] Figure 4 A schematic diagram of the density initialization model according to an embodiment of the present invention is shown, using the rock physics relationship ρ0=cZp0+d. Figure 2 The longitudinal wave impedance is transformed to obtain c=0.001, d=2.65.

[0133] Figure 5 A schematic diagram of an initial model of anisotropic parameter ε according to an embodiment of the present invention is shown. Using the rock physics relationship ε0=eZp0+f, Figure 2 The longitudinal wave impedance is transformed to obtain e=-0.027, f=0.756.

[0134] Figure 6 A schematic diagram of an initial model of anisotropic parameter δ according to an embodiment of the present invention is shown, using the rock physics relationship δ0 = gZp0 + h, Figure 2 The longitudinal wave impedance is transformed to obtain g=-0.014, h=0.345.

[0135] 1. Given angle gathers and logging curves, the logging curves include P-wave impedance, S-wave impedance, density, anisotropy parameters ε and δ;

[0136] 2. Perform rock physics analysis on the logging curves to establish rock physics characteristics and rock physics characteristic parameters;

[0137] 3. Perform post-stack inversion to obtain the longitudinal wave impedance Zp0;

[0138] 4. Use the rock physical characteristic parameters in step 2 to convert the longitudinal wave impedance Zp0 in step 3 to obtain Zs0, ρ0, ε0, δ0;

[0139] 5. Perform median filtering on Zp0, Zs0, ρ0, ε0, δ0 in the above step 4 to obtain the smoothed Zp1, Zs1, ρ1, ε1, δ1, thereby obtaining the initial model for anisotropic inversion.

[0140] Example 4

[0141] Figure 7 A block diagram of an anisotropic inversion initial model building device according to an embodiment of the present invention is shown.

[0142] like Figure 7As shown, the anisotropic inversion initial model building device includes:

[0143] Parameter preparation module 201, determining the P-wave impedance and multiple rock physical characteristics for anisotropic inversion;

[0144] An expression building module 202 builds multiple expressions between the longitudinal wave impedance and multiple rock physical characteristics respectively;

[0145] A first calculation module 203 obtains the value of longitudinal wave impedance according to seismic data;

[0146] The second calculation module 204 substitutes the value of the longitudinal wave impedance into a plurality of expressions to calculate the corresponding rock physical characteristics;

[0147] The modeling module 205 establishes an initial anisotropic inversion model based on the longitudinal wave impedance and multiple rock physical characteristics.

[0148] As an option, rock physical characteristics include shear wave impedance, density, and anisotropy parameters.

[0149] As an alternative, the expressions are:

[0150] Zs=aZp+b

[0151] ρ=cZp+d

[0152] ε=eZp+f

[0153] δ=gZp+h

[0154] Among them, Zp is the P-wave impedance, Zs is the S-wave impedance, ρ is the density, ε and δ are anisotropy parameters that characterize the degree of anisotropy of the underground medium; a and b are the rock physical characteristic parameters of S-wave impedance and P-wave impedance; c and d are the rock physical characteristic parameters of density and P-wave impedance; e and f are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance; g and h are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance.

[0155] As an alternative, P-wave impedance can be obtained from seismic data by post-stack inversion.

[0156] As an option, it also includes:

[0157] After the corresponding rock physical characteristics are calculated, the rock physical characteristics are smoothed and filtered.

[0158] As an optional solution, the smoothing filter is a median filter.

[0159] As an optional solution, the time window of the median filter is between 100ms and 200ms.

[0160] Example 5

[0161] This embodiment provides an electronic device, which includes: a memory storing executable instructions; and a processor, which runs the executable instructions in the memory to implement the above-mentioned anisotropic inversion initial model establishment method.

[0162] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.

[0163] The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

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

[0165] Those skilled in the art should be able to understand that in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus and an interface, and these well-known structures should also be included in the protection scope of the present disclosure.

[0166] For detailed description of this embodiment, reference may be made to the corresponding descriptions in the aforementioned embodiments, which will not be repeated here.

[0167] Example 6

[0168] This embodiment provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the anisotropic inversion initial model establishment method is implemented.

[0169] According to the computer-readable storage medium of the embodiment of the present disclosure, non-transitory computer-readable instructions are stored thereon. When the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the above-mentioned methods of each embodiment of the present disclosure are executed.

[0170] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (e.g., memory card) and media with built-in ROM (e.g., ROM box).

[0171] Those skilled in the art should understand that the purpose of the above description of the embodiments of the present invention is only to exemplarily illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any given examples.

[0172] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. A method for establishing an anisotropic inversion initial model, It is characterized in that include: Determine P-wave impedance and multiple rock physical characteristics for anisotropic inversion; Establishing multiple expressions between the longitudinal wave impedance and multiple rock physical characteristics respectively; Obtaining the value of the longitudinal wave impedance according to seismic data; Substituting the value of the longitudinal wave impedance into a plurality of expressions respectively, and calculating the corresponding rock physical characteristics; An anisotropic inversion initial model is established based on the longitudinal wave impedance and the multiple rock physical characteristics.

2. The method for establishing anisotropic inversion initial model according to claim 1, in, The rock physical characteristics include shear wave impedance, density, and anisotropy parameters.

3. The anisotropic inversion initial model establishment method according to claim 2, in, The expressions are: Zs=aZp+b ρ=cZp+d ε=eZp+f δ=gZp+h Among them, Zp is the P-wave impedance, Zs is the S-wave impedance, ρ is the density, ε and δ are anisotropy parameters that characterize the degree of anisotropy of the underground medium; a and b are the rock physical characteristic parameters of S-wave impedance and P-wave impedance; c and d are the rock physical characteristic parameters of density and P-wave impedance; e and f are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance; g and h are the rock physical characteristic parameters of anisotropy parameters and P-wave impedance.

4. The method for establishing anisotropic inversion initial model according to claim 1, in, The longitudinal wave impedance is obtained through post-stack inversion based on the seismic data.

5. The method for establishing anisotropic inversion initial model according to claim 1, in, Also includes: After the corresponding rock physical characteristics are calculated, the rock physical characteristics are smoothed and filtered.

6. The method for establishing anisotropic inversion initial model according to claim 5, in, The smoothing filter is a median filter.

7. The method for establishing anisotropic inversion initial model according to claim 6, in, The time window of the median filter is between 100ms and 200ms.

8. An anisotropic inversion initial model building device, It is characterized in that include: Parameter preparation module, which determines the P-wave impedance and multiple rock physical characteristics for anisotropic inversion; An expression establishment module, for establishing a plurality of expressions between the longitudinal wave impedance and a plurality of rock physical characteristics respectively; A first calculation module obtains the value of the longitudinal wave impedance according to seismic data; A second calculation module substitutes the value of the longitudinal wave impedance into a plurality of expressions to calculate the corresponding rock physical characteristics; A modeling module is used to establish an anisotropic inversion initial model based on the longitudinal wave impedance and the multiple rock physical characteristics.

9. An electronic device, It is characterized in that The electronic device comprises: A memory storing executable instructions; A processor, wherein the processor runs the executable instructions in the memory to implement the anisotropic inversion initial model establishment method according to any one of claims 1 to 7.

10. A computer-readable storage medium, It is characterized in that The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the anisotropic inversion initial model establishment method according to any one of claims 1 to 7 is implemented.

Images