Seismic exploration imaging method, device and equipment based on anisotropic velocity model of underground layer medium, medium and product

Abstract

The invention discloses a seismic exploration imaging method, device, equipment, medium and product based on an anisotropic velocity model of a subsurface medium, and relates to the field of seismic exploration, and the method comprises the steps: extracting an isotropic velocity curve of a logging coordinate position, carrying out the function operation processing of the isotropic velocity curve and the velocity curve, and determining a velocity difference model; accumulating the velocity difference model to a known isotropic velocity model according to function operation, determining an anisotropic velocity model, and determining a subsurface medium variation coefficient and subsurface medium longitudinal wave anisotropy; carrying out VTI pre-stack depth migration by using the anisotropy speed, the underground layer medium variation coefficient and the underground layer longitudinal wave anisotropy, and determining an inclination angle parameter and an azimuth angle parameter of the underground layer medium; and TTI pre-stack depth migration imaging processing is carried out to complete high-precision imaging of the subsurface medium, so that the efficiency of solving the anisotropic speed is improved, and subsurface medium imaging better conforms to the geological evolution law, thereby carrying out the next exploration survey as soon as possible.

Description

Technical Field

[0001] The present application relates to the field of seismic surveying, and in particular to a seismic exploration imaging method, device, equipment, medium and product based on anisotropic velocity model of underground medium. Background Art

[0002] In seismic exploration, the characteristics of geological structural targets gradually change from large, shallow and simple to small, deep and complex. Therefore, the requirements for seismic data processing and imaging are getting higher and higher. However, a fact that cannot be ignored is that underground stratum media generally exhibit anisotropic characteristics. In order to achieve high-precision imaging of underground geology, anisotropic prestack depth migration processing technology has become an indispensable means. Velocity modeling technology is the core of prestack depth migration imaging, and accurate velocity models are the basis of advanced migration imaging algorithms. Seismic anisotropic media refers to media in which elastic waves propagate at different speeds in different directions. Compared with the isotropic assumption, the prestack depth migration imaging technology based on the TTI anisotropic assumption is closer to the actual underground layer medium, so it is widely used in industrial production. However, the development of anisotropic velocity modeling is relatively lagging.

[0003] Generally speaking, the study of anisotropy starts with VTI. In 1986, Thomsen proposed three parameters to characterize the elastic properties of P-wave VTI media: vertical propagation velocity Vp0, P-wave anisotropy ε, and coefficient of variation δ. For TTI anisotropic media, two more parameters of formation dip and azimuth θ were added based on the three parameters of VTI media. There are five parameters in total. Among them, the δ parameter is a very important parameter in TTI media. Although its influence is not as great as the Vp0 parameter, if the δ parameter is inaccurate, it will change the morphology of the structural boundary and the depth position, leading to geological interpretation errors and ultimately affecting the drilling results.

[0004] At present, the commonly used pre-stack depth migration imaging technology in seismic exploration is to extract the δ parameter using the thickness ratio method and then establish an anisotropic velocity model.

[0005] The calculation formula of δ parameter is:

[0006] Among them, δn is the coefficient of variation of the nth underground medium layer; △ZI is the thickness of the nth layer of the isotropic migration profile; △ZA is the thickness of the nth layer of the geological stratification of the logging data.

[0007] The thickness of the stratum imaged by prestack depth migration is proportional to the velocity, that is, the greater the velocity, the thicker the stratum medium to be processed and imaged. If the isotropic velocity V of the underground stratum medium is known, p Under the condition of The anisotropic velocity V of the underground medium can be obtainedp0 :

[0008] Among them, V p is the isotropic velocity; δ is the coefficient of variation of the underground medium.

[0009] By formula It can be seen that the anisotropic velocity V is obtained p0 First, the coefficient of variation δ is calculated based on the actual layer thickness of the well logging and the thickness interpreted by the seismic imaging. Then, the coefficient of variation δ and the isotropic velocity V p Find the anisotropic velocity V p0 In addition, it is well known that when performing prestack depth migration velocity inversion, the velocity of the underlying strata is affected by the velocity of the overlying strata, that is, the coefficient of variation δ of the first layer is calculated. 1 and the anisotropic velocity V p01 Then calculate the coefficient of variation δ of the second layer 2 and the anisotropic velocity V p02 , calculated layer by layer, and finally established the anisotropic velocity model V p0 The efficiency of obtaining the anisotropic velocity and the anisotropic velocity model by repeated iterations is low, and subsequent seismic exploration work cannot be carried out in a timely manner. Summary of the invention

[0010] The purpose of this application is to provide a seismic exploration imaging method, device, equipment, medium and product based on the anisotropic velocity model of underground medium, which solves the problems of low precision in seismic exploration imaging processing and large error in well seismic depth.

[0011] To achieve the above objectives, this application provides the following solutions:

[0012] In a first aspect, the present application provides a seismic exploration imaging method based on anisotropic velocity model of underground medium, comprising:

[0013] Based on the isotropic velocity model, the isotropic velocity curve at the logging position is extracted, the logging velocity curve is edited, the isotropic velocity curve and the logging velocity curve are processed, and a velocity difference model is determined.

[0014] The velocity difference model is added to the isotropic velocity model according to a function operation, an anisotropic velocity model is determined and anisotropic velocity is obtained.

[0015] Based on the anisotropic velocity model, the coefficient of variation of the underground medium is determined according to the isotropic velocity model.

[0016] The anisotropy of the underground medium is determined based on the coefficient of variation of the underground medium.

[0017] VTI prestack depth migration is performed on the anisotropic velocity model, the coefficient of variation of the underground medium and the anisotropy of the underground layer to determine the dip parameters and azimuth parameters of the underground medium.

[0018] TTI pre-stack depth migration imaging processing is performed on the anisotropic velocity, the coefficient of variation of the underground medium, the anisotropy of the underground medium, the dip parameter and the azimuth parameter to complete high-precision imaging of the underground medium.

[0019] In a second aspect, the present application provides a seismic exploration imaging device based on anisotropic velocity model of underground medium, comprising:

[0020] The velocity difference model determination module is used to extract the isotropic velocity curve at the logging position based on the isotropic velocity model, edit the logging velocity curve, process the isotropic velocity curve and the logging velocity curve, and determine the velocity difference model.

[0021] The anisotropic velocity model determination module is used to accumulate the velocity difference model to the isotropic velocity model according to a function operation, determine the anisotropic velocity model and obtain the anisotropic velocity.

[0022] The underground layer medium variation coefficient determination module is used to determine the underground layer medium variation coefficient according to the isotropic velocity model.

[0023] The underground medium anisotropy determination module is used to determine the underground medium anisotropy based on the underground medium variation coefficient.

[0024] The module for determining the dip parameters and azimuth parameters of the underground medium is used to carry out VTI prestack depth migration on the anisotropic velocity model, the coefficient of variation of the underground medium and the anisotropy of the underground layer to determine the dip parameters and azimuth parameters of the underground medium.

[0025] The TTI pre-stack depth migration imaging processing module is used to carry out TTI pre-stack depth migration imaging processing on the anisotropic velocity, the coefficient of variation of the underground medium, the anisotropy of the underground medium, the dip parameter and the azimuth parameter to complete high-precision imaging of the underground medium.

[0026] In a third aspect, the present application provides a computer device comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of any of the above-mentioned seismic exploration imaging methods based on anisotropic velocity models of underground media.

[0027] In a fourth aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of any of the above-mentioned seismic exploration imaging methods based on anisotropic velocity models of underground media.

[0028] In a fifth aspect, the present application provides a computer program product, including a computer program, which, when executed by a processor, implements the steps of any of the above-mentioned seismic exploration imaging methods based on anisotropic velocity models of underground media.

[0029] According to the specific embodiments provided in this application, this application has the following technical effects:

[0030] The present application provides a seismic exploration imaging method, device, equipment, medium and product based on anisotropic velocity model of underground medium, which extracts isotropic velocity curve at the logging position from the known isotropic velocity model, edits V Sp Velocity curve, two velocity curves are processed by function operation to determine the velocity difference model; the velocity difference model is accumulated to the isotropic velocity model according to the function operation to determine the anisotropic velocity model; based on the anisotropic velocity model, the coefficient of variation of the underground medium is determined according to the isotropic velocity; the anisotropy of the underground medium is determined based on the coefficient of variation of the underground medium; the anisotropy of the underground medium is carried out by using the anisotropic velocity, the coefficient of variation of the underground medium and the anisotropy of the underground medium to carry out VTI prestack depth migration, determine the dip parameters and azimuth parameters of the underground medium, and use the seismic data for in-depth processing and analysis to fully mine useful information; the anisotropic velocity, the coefficient of variation of the underground medium, the anisotropy of the underground medium, the dip parameters and the azimuth parameters are carried out by TTI prestack depth migration imaging processing to complete the high-precision imaging of the underground medium, reduce the well seismic depth error, and meet the needs of subsequent seismic exploration. sp The velocity and the isotropic velocity at the logging position are used to obtain the velocity difference at the logging position. The velocity difference model is established by constraining the structure model layer. The velocity difference model is accumulated on the isotropic velocity model to obtain the anisotropic velocity model, so as to directly obtain the anisotropic velocity without repeatedly iterating to obtain the anisotropic velocity. This improves the efficiency of obtaining the anisotropic velocity and seismic exploration imaging processing, so that subsequent seismic exploration can be carried out as soon as possible. BRIEF DESCRIPTION OF THE DRAWINGS

[0031] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0032] Figure 1 A schematic diagram of a flow chart of a seismic exploration imaging method based on an anisotropic velocity model of underground medium in one embodiment of the present application;

[0033] Figure 2 A schematic diagram of extracting δ parameters using a thickness ratio method provided in an embodiment of the present application;

[0034] Figure 3 A schematic diagram of a flow chart for obtaining anisotropic velocity model Vp0 and coefficient of variation δ provided in one embodiment of the present application;

[0035] Figure 4 A schematic diagram of the comparison between the isotropic velocity and the Vsp velocity provided in one embodiment of the present application;

[0036] Figure 5 A schematic diagram of the comparison between anisotropic velocity and Vsp velocity provided in one embodiment of the present application;

[0037] FIG6 is a schematic diagram of a comparison diagram of isotropic velocity and anisotropic velocity provided in an embodiment of the present application, FIG6(a) is a schematic diagram of isotropic velocity; FIG6(b) is a schematic diagram of anisotropic velocity;

[0038] FIG7 is a schematic diagram comparing the imaging effects of isotropic prestack depth migration and anisotropic prestack depth migration provided in an embodiment of the present application, wherein FIG7(a) is a schematic diagram of the imaging effect of isotropic prestack depth migration, and FIG7(b) is a schematic diagram of the imaging effect of anisotropic prestack depth migration;

[0039] Figure 8 A schematic diagram of the structure of a computer device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0040] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0041] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0042] The velocity model obtained by conventional prestack depth migration processing is isotropic velocity, while the logging Vsp velocity is anisotropic velocity. sp velocity and the isotropic velocity V at the logging location p, calculate the speed difference V through the calculation function o , speed difference V o Establish velocity difference model V by constraining the structure model layer orr , speed difference model V orr Accumulated to the isotropic velocity model V through function operation p On the top, we get the anisotropic velocity model V p0 , so that the anisotropic velocity model can be quickly established.

[0043] Based on the above analysis, this application proposes a seismic exploration imaging method based on the anisotropic velocity model of underground media, which saves time, shortens the cycle, improves the accuracy of the anisotropic velocity model, and effectively improves the imaging quality in complex structural areas.

[0044] This application proposes a seismic exploration imaging method based on anisotropic velocity model of underground media, which is highly professional and rigorous. The construction of velocity model of underground media is the core of prestack depth migration imaging technology, and anisotropic velocity can more accurately reflect the actual situation of underground media. The specific steps of this method include the following aspects:

[0045] (1) First, based on the known isotropic velocity model, the present application extracts the isotropic velocity curve at the logging location and edits and smoothes the logging V sp Anisotropic velocity curve.

[0046] (2) Secondly, this application uses the smoothed V sp The velocity is compared point by point with the isotropic velocity extracted in step (1), and the velocity difference between them is obtained through function operation. The velocity difference value is obtained by subtracting the velocities of two points at the same depth.

[0047] (3) Based on the consideration of drilling stratification information and the interpretation of geological strata, this application establishes a structural model and constrains the velocity difference in space to generate a velocity difference model.

[0048] (4) Finally, the present application adds the velocity difference model to the isotropic velocity model through function calculation, thereby generating the final anisotropic velocity model.

[0049] like Figure 1 As shown, the embodiment of the present application provides a seismic exploration imaging method based on anisotropic velocity model of underground medium, which specifically includes:

[0050] S1: Based on the isotropic velocity model, extract the isotropic velocity curve at the logging position, edit the logging velocity curve, process the isotropic velocity curve and the logging velocity curve, and determine the velocity difference model.

[0051] S2: adding the velocity difference model to the isotropic velocity model according to function operation, determining the anisotropic velocity model and obtaining the anisotropic velocity.

[0052] S3: Based on the anisotropic velocity model, determine the coefficient of variation of the underground medium according to the isotropic velocity model.

[0053] S4: Determine the anisotropy of the underground medium based on the coefficient of variation of the underground medium.

[0054] S5: Perform VTI prestack depth migration on the anisotropic velocity model, the coefficient of variation of the underground medium and the anisotropy of the underground medium to determine the dip parameters and azimuth parameters of the underground medium.

[0055] S6: Perform TTI pre-stack depth migration imaging processing on the anisotropic velocity, the coefficient of variation of the underground medium, the anisotropy of the underground medium, the dip parameter and the azimuth parameter to complete high-precision imaging of the underground medium.

[0056] Further, in an exemplary embodiment, step S1 may be replaced by the following steps.

[0057] S101: Based on a known isotropic velocity model, an isotropic velocity curve at a logging coordinate position is extracted, and the logging velocity curve is edited.

[0058] S102: Compare the logging velocity curve and the isotropic velocity curve to determine the velocity difference.

[0059] The Vsp velocity of the well logging is compared with the isotropic velocity function, and the velocity difference between them is obtained through function operation.

[0060] S103: Based on the velocity difference, a velocity difference model is constructed according to the logging geological stratification and geological layer interpretation.

[0061] Combined with logging geological stratification, the seismic profile layers are interpreted, and based on this as a constraint, the velocity difference is used to establish a velocity difference model.

[0062] S104: Establish a structural model based on the drilling stratification information and interpretation of the underground layer medium.

[0063] Based on drilling stratification information, the geological strata, the contact relationship between different geological strata, the fracture location and the thickness of the strata are tracked and interpreted, and a spatial structural model is established using computer technology.

[0064] S105: Determine a velocity difference model based on the constructed model and constraining the velocity difference in space.

[0065] The velocity difference model is added to the isotropic velocity model through function operation to obtain the anisotropic velocity model.

[0066] Further, in an exemplary embodiment, step S3 may be replaced by the following steps.

[0067] S301: Utilize formula Determine the coefficient of variation of underground medium; where V p is the isotropic velocity; V p0 is the anisotropic velocity; δ is the coefficient of variation of the underground medium.

[0068] Further, in an exemplary embodiment, step S4 may be replaced by the following steps.

[0069] S401: Determine the anisotropy of the underground medium according to a preset multiple of the coefficient of variation of the underground medium.

[0070] S402: Perform VTI prestack depth migration based on the anisotropy of the underground medium to determine prestack depth migration gathers.

[0071] S403: Determine the anisotropy of the underground medium according to the VTI prestack depth migration gathers; the accurate value of the anisotropy of the underground medium is determined by the leveling condition of the VTI prestack depth migration gathers.

[0072] Further, in an exemplary embodiment, step S5 may be replaced by the following steps.

[0073] S501: Directly extracting the dip parameters and azimuth parameters of the underground medium according to the seismic imaging data volume of the anisotropic velocity shift.

[0074] S502: Alternatively, the isotropic velocity pre-stack depth migration imaging data is applied to the time domain by applying the isotropic velocity ratio to obtain imaging data in the time domain; the time domain imaging data is converted to the depth domain by applying the anisotropic velocity to determine the inclination parameters and azimuth parameters of the underground medium.

[0075] A seismic exploration imaging method based on anisotropic velocity model of underground medium, the specific process is as follows Figure 3 As shown, it is different from the current popular method of calculating the coefficient of variation δ first and then calculating the anisotropic velocity V p0 , a method of iteratively inverting from shallow to deep layers layer by layer. The specific technical details of this application are as follows:

[0076] Extract the isotropic velocity function at the well coordinate position, logging V spThe velocity function is compared with the extracted velocity function, and the velocity difference between the two is obtained by functional operation; referring to the drilling data stratification, the seismic data imaging layer is interpreted, the structural model is established, the velocity difference function is constrained along the interpreted layer in space, and the velocity difference model is established; the velocity difference model is accumulated to the original isotropic velocity model through functional operation to complete the establishment of the anisotropic velocity model.

[0077] In Linux environment, the operation code for obtaining the speed difference function is:

[0078] awk'{intA,B;scanf("%d%d\n",&A,&B);printf("%d\n",AB);return;}.

[0079] In Linux environment, the calculation code for establishing anisotropic velocity model is as follows:

[0080] awk'{int C, D; scanf ("%d%d\n", &C, &D); printf ("%d\n", C+D); return;}.

[0081] Among them, A is the Vsp velocity function; B is the isotropic velocity function of the Vsp coordinate position; C is the velocity difference model; and D is the isotropic velocity model.

[0082] Taking the seismic geological survey project in a certain area as an example, TTI anisotropic pre-stack depth migration imaging processing was carried out. It is known that at the logging coordinate position, according to the isotropic velocity model that was previously corrected and optimized, the isotropic velocity function was extracted from it. By comparison, it was found that there was an obvious difference between the Vsp velocity and the isotropic velocity, such as Figure 4 As shown, the green line is the isotropic velocity function (curve) extracted at the logging location; the blue line is the Vsp velocity function (curve). In order to solve this problem, this application performs a function operation on the logging Vsp velocity and the extracted isotropic velocity to obtain the velocity difference between the two. Then, based on the drilling stratification data and the seismic profile layer interpretation, a relevant structural model is established. This application will use the structural model to spatially constrain the velocity difference at the logging location and establish a corresponding velocity difference model. This velocity difference model will be used for function operations and accumulated with the isotropic velocity model to obtain the establishment of an anisotropic velocity model. As shown Figure 5As shown in Figure 6, the comparison between the anisotropic velocity and the Vsp velocity at the logging position shows that there is a good correlation between the two. The blue line is the Vsp velocity function (curve); the velocity difference model is accumulated to the isotropic velocity model through function operation to obtain the anisotropic velocity model; the green line is the anisotropic velocity function (curve) obtained by extracting the anisotropic velocity model at the logging position. As shown in Figure 6, the difference between the isotropic velocity model and the anisotropic velocity model is compared. Figure 6 (a) is a schematic diagram of the isotropic velocity; Figure 6 (b) is a schematic diagram of the anisotropic velocity. It is found that the anisotropic velocity model can not only accurately present the details of the velocity changes within the same set of formations, but also has a gradual change trend between different formations. By obtaining the value of the coefficient of variation δ body, and taking 1.5 times of it as the value of the anisotropic parameter ε. In order to perform TTI prestack depth migration imaging processing, the isotropic prestack depth migration imaging data needs to be converted to the time domain using the isotropic velocity described above, and then converted to the depth domain according to the corresponding anisotropic velocity Vp0 ratio. Then, the present application extracts attribute parameters such as dip angle θ and azimuth angle φ of the formation based on the depth domain seismic data after proportional conversion. After the extraction of anisotropic parameters Vp0, δ, ε, θ, and φ is completed, TTI prestack depth migration imaging processing is performed. As shown in Figure 7, the isotropic prestack depth migration imaging section and the TTI anisotropic prestack depth migration imaging section are compared. Figure 7 (a) is a schematic diagram of the isotropic prestack depth migration imaging effect, and Figure 7 (b) is a schematic diagram of the anisotropic prestack depth migration imaging effect. The results show that the anisotropic prestack depth migration imaging wave group has clear features, significant fractures, and a small error with the well stratification depth, and the imaging quality is high. The above is the steps and results of TTI anisotropic prestack depth migration imaging processing based on an earthquake geological survey project in a certain area.

[0083] The application of this application in the field of geological survey has been proven to be practical and feasible. In actual projects, this application can bring significant results to geological surveys, further verifying its effectiveness and reliability in practice. Case practice in a certain area shows that this application can establish anisotropic velocity models with high precision, improve the quality of underground space imaging in practical applications, and reduce the errors of anisotropic pre-stack depth migration imaging depth and drilling stratification depth. Therefore, this application is of great significance to geological survey work.

[0084] The embodiment of the present application provides a seismic exploration device based on anisotropic velocity model of underground medium, specifically comprising:

[0085] The velocity difference model determination module is used to extract the isotropic velocity curve at the logging position based on the isotropic velocity model, edit the logging velocity curve, process the isotropic velocity curve and the logging velocity curve, and determine the velocity difference model.

[0086] The anisotropic velocity model determination module is used to accumulate the velocity difference model to the isotropic velocity model according to a function operation, determine the anisotropic velocity model and obtain the anisotropic velocity.

[0087] The underground layer medium variation coefficient determination module is used to determine the underground layer medium variation coefficient according to the isotropic velocity model.

[0088] The underground medium anisotropy determination module is used to determine the underground medium anisotropy based on the underground medium variation coefficient.

[0089] The module for determining the dip parameters and azimuth parameters of the underground medium is used to carry out VTI prestack depth migration on the anisotropic velocity model, the coefficient of variation of the underground medium and the anisotropy of the underground layer to determine the dip parameters and azimuth parameters of the underground medium.

[0090] The TTI pre-stack depth migration imaging processing module is used to carry out TTI pre-stack depth migration imaging processing on the anisotropic velocity, the coefficient of variation of the underground medium, the anisotropy of the underground medium, the dip parameter and the azimuth parameter to complete high-precision imaging of the underground medium.

[0091] In an exemplary embodiment, a computer device is provided. The computer device may be a server or a terminal. The internal structure diagram thereof may be as follows: Figure 8 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. The processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store video tag processing data. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a video tag processing method is implemented.

[0092] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0093] In an exemplary embodiment, a computer program product is provided, including a computer program, and when the computer program is executed by a processor, the steps in the above method embodiments are implemented.

[0094] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0095] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A seismic exploration imaging method based on anisotropic velocity model of underground medium, characterized in that: The seismic exploration imaging method based on the anisotropic velocity model of underground medium includes: Based on the isotropic velocity model, extracting the isotropic velocity curve at the logging position, editing the logging velocity curve, processing the isotropic velocity curve and the logging velocity curve, and determining the velocity difference model; Accumulating the velocity difference model to the isotropic velocity model according to a function operation, determining an anisotropic velocity model and obtaining anisotropic velocity; Based on the anisotropic velocity model, determining the coefficient of variation of the underground medium according to the isotropic velocity model; Determining the anisotropy of the underground medium based on the coefficient of variation of the underground medium; Perform VTI prestack depth migration on the anisotropic velocity model, the coefficient of variation of the underground medium and the anisotropy of the underground layer to determine the dip parameters and azimuth parameters of the underground medium; TTI pre-stack depth migration imaging processing is performed on the anisotropic velocity, the coefficient of variation of the underground medium, the anisotropy of the underground medium, the dip parameter and the azimuth parameter to complete high-precision imaging of the underground medium.

2. The seismic exploration imaging method based on the anisotropic velocity model of underground medium according to claim 1, characterized in that: Based on the isotropic velocity model, an isotropic velocity curve at a logging position is extracted, the logging velocity curve is edited, the isotropic velocity curve and the logging velocity curve are processed, and a velocity difference model is determined, which specifically includes: Based on the known isotropic velocity model, the isotropic velocity curve at the logging coordinate position is extracted and the logging velocity curve is edited; comparing the logging velocity curve with the isotropic velocity curve to determine a velocity difference; Based on the velocity difference, a velocity difference model is constructed according to the logging geological stratification and geological horizon interpretation.

3. The seismic exploration imaging method based on anisotropic velocity model of underground medium according to claim 2, characterized in that: Based on the velocity difference, according to the logging geological stratification and geological horizon interpretation, a velocity difference model is constructed, which specifically includes: Establish structural model based on drilling stratification information and interpretation of underground strata media; Based on the construction model and spatially constraining the velocity difference value, a velocity difference model is determined.

4. The seismic exploration imaging method based on anisotropic velocity model of underground medium according to claim 1, characterized in that: Based on the anisotropic velocity model, determining the coefficient of variation of underground medium according to the isotropic velocity model specifically includes: Using the formula Determine the coefficient of variation of underground medium; where V p is the isotropic velocity; V p0 is the anisotropic velocity; δ is the coefficient of variation of the underground medium.

5. The seismic exploration imaging method based on anisotropic velocity model of underground medium according to claim 1, characterized in that: Determining the anisotropy of the underground medium based on the underground medium variation coefficient specifically includes: Determining the anisotropy of the underground medium according to a preset multiple of the coefficient of variation of the underground medium; The anisotropy of the underground medium is used to perform VTI prestack depth migration to determine prestack depth migration gathers; The anisotropy of the underground medium is determined according to the VTI pre-stack depth migration gathers; the accurate value of the anisotropy of the underground medium is determined by the leveling condition of the VTI pre-stack depth migration gathers.

6. The seismic exploration imaging method based on anisotropic velocity model of underground medium according to claim 1, characterized in that: The TTI pre-stack depth migration imaging process is performed on the anisotropic velocity, the coefficient of variation of the underground medium, the anisotropy of the underground medium, the dip parameter and the azimuth parameter to complete high-precision imaging of the underground medium, specifically including: Directly extracting the dip parameters and azimuth parameters of the underground layer medium according to the seismic imaging data volume of the anisotropic velocity shift; Alternatively, the pre-stack depth migration imaging data of the isotropic velocity is applied to the time domain by applying the isotropic velocity ratio to obtain imaging data in the time domain; the imaging data in the time domain is converted to the depth domain by applying the anisotropic velocity to determine the inclination parameters and azimuth parameters of the underground medium.

7. A seismic exploration imaging device based on anisotropic velocity model of underground medium, characterized in that: The seismic exploration imaging device based on the anisotropic velocity model of underground medium specifically comprises: A velocity difference model determination module is used to extract the isotropic velocity curve at the logging position based on the isotropic velocity model, edit the logging velocity curve, process the isotropic velocity curve and the logging velocity curve, and determine the velocity difference model; An anisotropic velocity model determination module is used to accumulate the velocity difference model to the isotropic velocity model according to a function operation, determine the anisotropic velocity model and obtain the anisotropic velocity; A module for determining the coefficient of variation of underground medium, used for determining the coefficient of variation of underground medium according to the isotropic velocity model; A module for determining anisotropy of underground medium, used for determining anisotropy of underground medium based on the coefficient of variation of underground medium; A module for determining the dip parameters and azimuth parameters of the underground medium, for carrying out VTI prestack depth migration on the anisotropic velocity model, the coefficient of variation of the underground medium and the anisotropy of the underground layer, to determine the dip parameters and azimuth parameters of the underground medium; The TTI pre-stack depth migration imaging processing module is used to carry out TTI pre-stack depth migration imaging processing on the anisotropic velocity, the coefficient of variation of the underground medium, the anisotropy of the underground medium, the dip parameter and the azimuth parameter to complete high-precision imaging of the underground medium.

8. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the seismic exploration imaging method based on the anisotropic velocity model of the underground medium as described in any one of claims 1-6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the seismic exploration imaging method based on the anisotropic velocity model of the underground medium described in any one of claims 1 to 6 is implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the seismic exploration imaging method based on the anisotropic velocity model of the underground medium described in any one of claims 1 to 6 is implemented.

Images