Transformed shear wave dynamic grid finite difference forward modeling method, device, equipment and medium

Abstract

The invention discloses a converted shear wave dynamic grid finite difference forward modeling method, device and equipment and a medium, and the method comprises the steps: setting forward modeling parameters, and inputting a geophysical parameter model; performing dynamic mesh generation processing on the input geophysical parameter model to obtain a geophysical parameter model after mesh optimization; performing converted shear wave finite difference numerical simulation based on a converted shear wave dynamic grid finite difference wave equation by adopting the grid optimized geophysical parameter model; and recording the amplitude of the converted shear wave at the detection point based on the finite difference wave equation of the dynamic grid of the converted shear wave according to preset observation system parameters to obtain a converted shear wave seismic record of forward modeling. Therefore, according to the method, the precision of the forward modeling wave field is hardly influenced, and the method can be widely applied to seismic data processing.

Description

Technical Field

[0001] The present invention relates to a converted shear wave dynamic grid finite-difference forward modeling method, device, equipment and medium, and relates to the field of geophysical exploration for oil and gas. Background Art

[0002] As a key link in seismic data processing, forward modeling cannot be ignored in the field of seismic exploration and development. The wave equation forward modeling method is an important basis for wave equation-based velocity modeling and migration imaging, and its efficiency has a significant impact on the cost of seismic data processing.

[0003] With the continuous deepening of the seismic exploration field, the data obtained from seismic data acquisition has gradually been upgraded from two-dimensional to three-dimensional and extended from single-component to multi-component, resulting in a significant increase in data volume. Therefore, how to improve the calculation efficiency, reduce the computer memory burden, and lower the cost of seismic data processing while ensuring the simulation accuracy has become a key problem to be solved in the field of oil and gas exploration and development.

[0004] Therefore, there is an urgent need for a new converted shear wave forward modeling method that can balance the forward modeling accuracy and calculation efficiency. Summary of the Invention

[0005] The present invention aims to solve at least one of the technical problems existing in the prior art. For this reason, in view of the above problems, the object of the present invention is to provide a converted shear wave dynamic grid finite-difference forward modeling method, device, equipment and medium, which can solve the problems of low calculation efficiency and serious numerical dispersion in the current conventional converted shear wave forward modeling technology, and efficiently and accurately simulate the converted shear wave seismic wave field.

[0006] To achieve the above object of the invention, the technical solution adopted by the present invention is as follows: In a first aspect, the present invention provides a converted shear wave dynamic grid finite-difference forward modeling method, including: Setting forward modeling parameters and inputting a geophysical parameter model; Performing dynamic grid meshing on the input geophysical parameter model to obtain a geophysical parameter model with optimized grid; Using the geophysical parameter model with optimized grid, performing converted shear wave finite-difference numerical simulation based on the converted shear wave dynamic grid finite-difference wave equation; According to the preset observation system parameters, recording the converted shear wave amplitude at the geophone based on the converted shear wave dynamic grid finite-difference wave equation to obtain the converted shear wave seismic record of the forward modeling.

[0007] Furthermore, the forward modeling parameters include the distribution position of geophones, the number of grid points in the horizontal and depth directions of the model and , grid spacing and , the dominant frequency of the source wavelet and the time sampling interval .

[0008] Furthermore, the geophysical parameter model includes a P-wave velocity model, an S-wave velocity model, and an S-wave velocity perturbation model.

[0009] Furthermore, perform dynamic grid meshing on the input geophysical parameter model to obtain a geophysical parameter model with optimized grid, including: Calculate the dynamic grid spacing of the velocity model in the depth direction according to the velocity model dynamic grid spacing formula ; Starting from , by gradually increasing the trial grid step size until , obtain a series of grid point depths , construct a new dynamic grid point depth sequence, and obtain a geophysical parameter model with optimized grid.

[0010] Furthermore, set the starting depth position of the dynamic grid meshing to be deeper than the deepest shot point and receiver point positions.

[0011] Furthermore, the velocity model dynamic grid spacing formula is: ; where is the minimum velocity value of the velocity model at a certain depth, is the longitudinal grid spacing of the dynamic grid velocity model, is the dominant frequency of the source wavelet used in the numerical simulation.

[0012] Furthermore, convert the S-wave dynamic grid finite-difference wave equation to: ; where is the P-wave background velocity, is the S-wave background velocity, is the S-wave velocity perturbation, is the incident P-wave wavefield, is the source time function, is the reflected converted S-wave wavefield, is the S-wave reflection source at the conversion point, is the derivative of the dynamic grid depth coordinate with respect to the original grid depth coordinate .

[0013] In the second aspect, the present invention also provides a converted S-wave dynamic grid finite-difference forward simulation device, including: An input unit, configured to set forward modeling parameters and input a geophysical parameter model; A grid optimization unit, configured to perform dynamic grid meshing processing on the input geophysical parameter model to obtain a geophysical parameter model with optimized grid; A numerical simulation unit, configured to perform converted shear wave finite-difference numerical simulation based on the converted shear wave dynamic grid finite-difference wave equation by using the geophysical parameter model with optimized grid; A forward modeling unit, configured to record the converted shear wave amplitude at the geophone points based on the converted shear wave dynamic grid finite-difference wave equation according to the preset acquisition system parameters, and obtain the converted shear wave seismic record of the forward modeling.

[0014] In a third aspect, the present invention further provides an electronic device, including: at least one processor; and a memory communicatively connected to the processor; wherein, the memory stores instructions executable by the processor, and when the instructions are executed by the processor, the processor is enabled to execute any one of the methods described above.

[0015] In a fourth aspect, the present invention further provides a computer-readable storage medium storing one or more programs, where the one or more programs include computer instructions for causing a computer to execute any one of the methods described above.

[0016] Due to the above technical solutions adopted by the present invention, it has the following characteristics: Compared with the prior art, the present invention has the dual advantages of high simulation accuracy and high calculation efficiency. When performing forward modeling, the calculation amount and the computer memory occupation are significantly reduced, and a faster calculation speed is achieved. At the same time, the method of the present invention can keep the wave field accuracy of the forward modeling almost unaffected, and can be widely applied to seismic data processing. Description of the Drawings

[0017] By reading the detailed description of the preferred embodiments below, various other advantages and benefits will become clear to those of ordinary skill in the art. The drawings are only for the purpose of showing the preferred embodiments and are not considered to be a limitation of the present invention. Throughout the drawings, the same reference numerals are used to represent the same components. In the drawings: Figure 1 is a flowchart of an embodiment of the present invention; Figure 2 is a diagram of the input longitudinal wave velocity model of an embodiment of the present invention; Figure 3 is a diagram of the input shear wave velocity model in an embodiment of the present invention; Figure 4 is a diagram of the input shear wave velocity perturbation in an embodiment of the present invention; Figure 5It is a diagram of the dynamic grid longitudinal wave velocity model in the embodiments of the present invention; Figure 6 It is a diagram of the dynamic grid shear wave velocity model in the embodiments of the present invention; Figure 7 It is a diagram of the dynamic grid shear wave velocity perturbation in the embodiments of the present invention; Figure 8 It is a simulated seismic record of the dynamic grid converted shear wave in the embodiments of the present invention; Figure 9 It is a structural diagram of the electronic device in the embodiments of the present invention. Detailed implementation manners

[0018] It should be understood that the terms used herein are for the purpose of describing specific example embodiments only and are not intended to be limiting. Unless the context clearly indicates otherwise, the singular forms "a", "an", and "the" as used herein may also include the plural forms. The terms "comprising", "including", "containing", and "having" are inclusive and thus specify the presence of the stated features, steps, operations, elements, and / or components, but do not preclude the presence or addition of one or more other features, steps, operations, elements, components, and / or combinations thereof. The method steps, processes, and operations described herein are not to be construed as necessarily requiring them to be performed in the particular order described or illustrated, unless the order of performance is explicitly stated. It should also be understood that additional or alternative steps may be used.

[0019] Although the terms first, second, third, etc. may be used in this document to describe multiple elements, components, regions, layers, and / or sections, these elements, components, regions, layers, and / or sections should not be limited by these terms. These terms may be used only to distinguish one element, component, region, layer, or section from another. Unless the context clearly indicates otherwise, terms such as "first", "second", and other numerical terms do not imply an order or sequence when used in this document. Thus, the first element, component, region, layer, or section discussed below may be referred to as the second element, component, region, layer, or section without departing from the teachings of the example embodiments.

[0020] For ease of description, spatial relative relationship terms may be used in this document to describe the relationship of one element or feature shown in the figure with respect to another element or feature, such as "inside", "outside", "inner side", "outer side", "below", "above", etc. Such spatial relative relationship terms are intended to include different orientations of the device in use or operation in addition to the orientations depicted in the figure.

[0021] Due to the serious dispersion phenomenon in the conventional coarse-grid converted S-wave finite-difference forward simulation and the huge computational cost of using the conventional fine grid. The present invention provides a converted S-wave dynamic grid finite-difference forward simulation method, device, equipment and medium that take into account both computational efficiency and simulation accuracy, including: setting forward simulation parameters and inputting a geophysical parameter model; performing dynamic grid meshing on the input geophysical parameter model to obtain an optimized geophysical parameter model; using the optimized geophysical parameter model to perform converted S-wave finite-difference numerical simulation based on the converted S-wave dynamic grid finite-difference wave equation; and recording the converted S-wave amplitude at the geophone point based on the converted S-wave dynamic grid finite-difference wave equation according to the pre-set observation system parameters to obtain the converted S-wave seismic record of the forward simulation. Therefore, the present invention can effectively improve the computational efficiency of the converted S-wave finite-difference forward simulation while maintaining the numerical simulation accuracy.

[0022] Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although the exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided so that the present invention can be more thoroughly understood and the scope of the present invention can be fully conveyed to those skilled in the art.

[0023] Example 1: As Figure 1 shown, the converted S-wave dynamic grid finite-difference forward simulation method provided in this embodiment includes: S1. Set forward simulation parameters such as the observation system and seismic wavelet, and input a geophysical parameter model, where the geophysical parameter model includes a P-wave velocity model, an S-wave velocity model, and an S-wave velocity perturbation model of the original conventional grid, and the sizes of all rectangular grids in the model are the same.

[0024] In this embodiment, the forward simulation parameters include, but are not limited to, the distribution positions of shot points and geophone points, the number of grid points in the horizontal and depth directions of the model and , the grid spacing and , the main frequency of the source wavelet and the time sampling interval , etc.

[0025] S2. Perform dynamic grid meshing on the input geophysical parameter model to obtain an optimized geophysical parameter model.

[0026] In this embodiment, resample the P-wave velocity model, S-wave velocity model, and S-wave velocity perturbation model of the original conventional grid along the depth direction, and perform dynamic grid meshing. The formula for calculating the dynamic grid spacing of the velocity model is: (1); Wherein, is the minimum velocity value of the longitudinal wave velocity model (shear wave velocity model or shear wave velocity perturbation model) at a certain depth, is the longitudinal grid spacing of the dynamic grid longitudinal wave velocity model (shear wave velocity model or shear wave velocity perturbation model), is the main frequency of the seismic source wavelet adopted in the numerical simulation.

[0027] In this embodiment, the specific process of performing dynamic grid meshing processing on all input geophysical parameter models based on the above formula (1) to obtain the geophysical parameter model with optimized grid is as follows: First, calculate the dynamic grid spacing of the velocity model in the depth direction according to formula (1) ; Then, starting from , by gradually increasing the trial grid step size until , a series of grid point depths are obtained, and a new dynamic grid point depth sequence is constructed. Therefore, the conventional grid coordinates are converted into grid coordinates with a dynamically varying grid spacing in the depth direction. At the grid point depth , the interpolated new velocity model is the dynamic grid model.

[0028] Furthermore, the starting depth position of the dynamic grid meshing should be deeper than the deepest shot point and geophone point positions.

[0029] S3. Using the optimized longitudinal and shear wave velocity models and shear wave velocity perturbation model of the grid, perform converted shear wave finite difference numerical simulation based on the converted shear wave dynamic grid finite difference wave equation.

[0030] In this embodiment, based on the parameter model of the dynamic grid meshing, the converted shear wave is simulated, and the converted shear wave seismic record at time is calculated. The conversion of the conventional grid coordinates to the dynamic grid coordinates can be understood as a coordinate mapping, expressed as: (2); In the formula, , represents the dynamic grid position coordinates obtained after coordinate conversion, is the mapping from the conventional grid depth to the dynamic grid depth coordinates, ([[]]END]] x , y ) is the original grid position coordinates.

[0031] The derivative of the dynamic grid depth coordinate with respect to the original grid depth coordinate is expressed as: (3); According to the derivative rule of the inverse function, derive the first-order and second-order derivatives of with respect to (4); (5); Derive that the second-order derivative of the wave field with respect to depth in the dynamic grid coordinate system is: (6); Under the conventional grid meshing, the simulated converted shear wave wave equation is expressed as: (7); Substitute formula (6) into formula (7), and the expression form of the dynamic grid converted shear wave wave equation can be obtained, that is, the wave equation used for forward modeling in the present invention, namely the converted shear wave dynamic grid finite-difference wave equation is: (8); where is the longitudinal wave background velocity, is the shear wave background velocity, is the shear wave velocity perturbation, is the incident longitudinal wave field, is the source time function, is the reflected converted shear wave field, is the shear wave reflection source at the conversion point, and the above formula can be solved by the high-order finite-difference method.

[0032] S4. According to the preset observation system parameters, record the converted shear wave amplitude at the geophone based on the above-mentioned converted shear wave dynamic grid finite-difference wave equation, and obtain the converted shear wave seismic record of the forward modeling.

[0033] The following specifically illustrates the specific application of the converted shear wave dynamic grid finite-difference forward modeling method of the present invention through specific embodiments.

[0034] As Figure 1 shows the method flow of the present invention. The input longitudinal wave velocity model is as Figure 2 shown, the input shear wave velocity model is as Figure 3 shown, and the input shear wave velocity perturbation model is as Figure 4As shown in the figure. The forward simulation parameters are as follows: the number of horizontal and vertical grid points in the velocity model are 1534 and 552 respectively, the horizontal and vertical grid spacings are both 6 m, the time sampling interval is 0.5 ms, and the source uses a Ricker wavelet with a dominant frequency of 15 Hz. The pre-set acquisition system parameters are as follows: the horizontal coordinate of the shot point is 4596 m, the depth is 24 m, the number of geophone points is 767, the geophone interval is 12 m, and they are evenly distributed at a depth of 24 m near the surface. Figure 5 is the longitudinal wave velocity model of the dynamic grid, Figure 6 is the transverse wave velocity model of the dynamic grid, Figure 7 is the transverse wave velocity perturbation model of the dynamic grid. The number of longitudinal grid points in the dynamic grid model is reduced from 552 to 296, and the number of grid points is reduced by 46.3%. The percentage reduction in the number of grid points is the memory size that can be reduced by the present invention, which is 46.3%. Applying the converted shear wave dynamic grid finite-difference forward simulation technology to conduct the converted shear wave forward simulation, a converted shear wave seismic record is obtained, as Figure 8 shown in the figure. In summary, the dynamic grid method adopted by the present invention can effectively improve the calculation efficiency and reduce the memory occupied by the calculation without reducing the forward simulation accuracy.

[0035] Embodiment 2: The above Embodiment 1 provides a converted shear wave dynamic grid finite-difference forward simulation method. Correspondingly, this embodiment provides a converted shear wave dynamic grid finite-difference forward simulation device. The device provided in this embodiment can implement the converted shear wave dynamic grid finite-difference forward simulation method of Embodiment 1, and this device can be implemented by software, hardware, or a combination of software and hardware. For the convenience of description, when describing this embodiment, various units are described separately according to their functions. Of course, in implementation, the functions of each unit can be implemented in the same or multiple software and / or hardware. For example, the device may include integrated or separate functional modules or functional units to execute the corresponding steps in the methods of Embodiment 1. Since the device in this embodiment is basically similar to the method embodiment, the description process of this embodiment is relatively simple, and the relevant parts can refer to the partial description of Embodiment 1. The embodiment of the converted shear wave dynamic grid finite-difference forward simulation device provided by the present invention is only illustrative.

[0036] Specifically, the present invention provides a converted shear wave dynamic grid finite-difference forward simulation device, including: An input unit, configured to set forward simulation parameters and input a geophysical parameter model; A grid optimization unit, configured to perform dynamic grid meshing processing on the input geophysical parameter model to obtain a geophysical parameter model with optimized grid; A numerical simulation unit, configured to perform a converted shear wave finite-difference numerical simulation based on the converted shear wave dynamic grid finite-difference wave equation by using a geophysical parameter model with optimized grid. A forward simulation unit, configured to record the converted shear wave amplitude at a geophone based on the converted shear wave dynamic grid finite-difference wave equation according to preset observation system parameters, and obtain a converted shear wave seismic record of the forward simulation.

[0037] Embodiment 3: This embodiment provides an electronic device corresponding to the converted shear wave dynamic grid finite-difference forward simulation method provided in Embodiment 1. The electronic device can be an electronic device for a client, such as a mobile phone, a laptop computer, a tablet computer, a desktop computer, etc., to execute the method of Embodiment 1.

[0038] As Figure 9 shown, the electronic device includes a processor, a memory, a communication interface, and a bus. The processor, the memory, and the communication interface are connected through the bus to complete communication with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, an Extended Industry Standard Component (EISA) bus, etc. The memory stores a computer program that can run on the processor. When the processor runs the computer program, it executes the method of Embodiment 1. The implementation principle and technical effects are similar to those of Embodiment 1 and will not be elaborated here. Those skilled in the art can understand that Figure 9 the structure shown in

[0039] In a preferred embodiment, when the logical instructions in the above-mentioned memory can be implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present application, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present application. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), and optical discs that can store program codes.

[0040] In a preferred embodiment, the processor can be various types of general-purpose processors such as a central processing unit (CPU) and a digital signal processor (DSP), which are not limited herein.

[0041] Embodiment 4: This embodiment provides a computer-readable storage medium storing one or more programs. The one or more programs include computer instructions. When the computer instructions are executed by a computer, the computer is caused to execute the method provided in Embodiment 1 above.

[0042] Embodiment 5: This embodiment provides a computer program product. The computer program product may be a computer program stored on a computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer is capable of executing the method provided in Embodiment 1 above. Its implementation principle and technical effects are similar to those of Embodiment 1 and will not be elaborated herein.

[0043] In a preferred embodiment, the computer-readable storage medium can be a tangible device that holds and stores instructions used by an instruction execution device. For example, it can be, but is not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any combination of the above. The computer-readable storage medium stores computer program instructions that cause the computer to execute the method provided in Embodiment 1 above.

[0044] This application is described with reference to the flowcharts and / or block diagrams of methods, apparatuses (devices), and computer program products according to embodiments of the present application. It should be understood that each flow and / or block in the flowchart and / or block diagram, as well as the combination of flows and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, such that the instructions executed by the processor of the computer or other programmable data processing devices generate a device for implementing the functions specified in one flow Figure 1 or more flows and / or blocks Figure 1 or a device for implementing the functions specified in one block or more blocks.

[0045] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing devices to work in a specific manner, such that the instructions stored in the computer-readable memory generate a manufactured article including an instruction device that implements the functions specified in one flow Figure 1 or more flows and / or blocks Figure 1 or a device for implementing the functions specified in one block or more blocks.

[0046] These computer program instructions can also be loaded onto a computer or other programmable data processing devices, such that a series of operation steps are executed on the computer or other programmable devices to generate a computer-implemented process, and thus the instructions executed on the computer or other programmable devices provide steps for implementing the functions specified in one flow Figure 1 or more flows and / or blocks Figure 1 or a device for implementing the functions specified in one block or more blocks.

[0047] Each embodiment in this specification is described in a progressive manner. The same or similar parts among the embodiments can be referred to each other, and the differences between each embodiment and other embodiments are emphasized. In the description of this specification, the descriptions with reference to terms such as "a preferred embodiment", "furthermore", "specifically", "in this embodiment", etc. mean that the specific features, structures, materials, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the embodiments of this specification. In this specification, the schematic expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described can be combined in a suitable manner in any one or more embodiments or examples. In addition, without contradiction, those skilled in the art can combine and combine the different embodiments or examples described in this specification and the features of different embodiments or examples.

[0048] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some of the technical features. These modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for converting shear wave dynamic grid finite difference forward modeling, characterized in that: include: Set forward modeling parameters and input geophysical parameter model; Perform dynamic meshing processing on the input geophysical parameter model to obtain the geophysical parameter model after mesh optimization; The geophysical parameter model after grid optimization is used to perform the converted shear wave finite difference numerical simulation based on the converted shear wave dynamic grid finite difference wave equation. According to the preset observation system parameters, the converted shear wave amplitude is recorded at the detection point based on the converted shear wave dynamic grid finite difference wave equation to obtain the converted shear wave seismic record of the forward simulation.

2. The converted shear wave dynamic grid finite difference forward modeling method according to claim 1 is characterized in that: Forward modeling parameters include the distribution of receiver points, the number of grid points in the horizontal and depth directions of the model, and the number of grid points in the horizontal and depth directions of the model. and , Grid Spacing and , source wavelet dominant frequency and time sampling interval .

3. The converted shear wave dynamic grid finite difference forward modeling method according to claim 1 is characterized in that: The geophysical parameter model includes the P-wave velocity model, S-wave velocity model and S-wave velocity disturbance model.

4. The converted shear wave dynamic grid finite difference forward modeling method according to claim 3 is characterized in that: Perform dynamic meshing processing on the input geophysical parameter model to obtain a geophysical parameter model after mesh optimization, including: Calculate the dynamic grid spacing of the velocity model in the depth direction according to the dynamic grid spacing formula of the velocity model ; from Initially, by gradually increasing the trial grid step size until , get a series of grid point depths , construct a new dynamic grid point depth sequence and obtain the geophysical parameter model after grid optimization.

5. The converted shear wave dynamic grid finite difference forward simulation method according to claim 4 is characterized in that: Set the starting depth of dynamic meshing to a position deeper than the deepest shot and receiver points.

6. The converted shear wave dynamic grid finite difference forward modeling method according to claim 4 is characterized in that: The dynamic grid spacing formula of the velocity model is: ； in, is the minimum velocity value of the velocity model at a certain depth, is the longitudinal grid spacing of the dynamic grid velocity model, is the dominant frequency of the source wavelet used in the numerical simulation.

7. The converted shear wave dynamic grid finite difference forward modeling method according to claim 6 is characterized in that: Convert the shear wave dynamic grid finite difference wave equation to: ； in, is the longitudinal wave background velocity, is the shear wave background velocity, is the shear wave velocity disturbance, is the incident longitudinal wave field, is the source time function, is the reflected converted shear wave field, is the shear wave reflection source at the transition point, The dynamic mesh depth coordinate is relative to the original mesh depth coordinate The derivative of .

8. A finite difference forward modeling device for converting shear waves, characterized in that: include: An input unit is configured to set forward modeling parameters and input geophysical parameter models; A grid optimization unit is configured to perform dynamic grid generation processing on the input geophysical parameter model to obtain a geophysical parameter model after grid optimization; The numerical simulation unit is configured to use the geophysical parameter model after grid optimization to perform a converted shear wave finite difference numerical simulation based on a converted shear wave dynamic grid finite difference wave equation; The forward modeling unit is configured to record the converted shear wave amplitude at the detection point based on the converted shear wave dynamic grid finite difference wave equation according to the preset observation system parameters, and obtain the converted shear wave seismic record of the forward modeling.

9. An electronic device, characterized in that: include: at least one processor; and a memory communicatively connected to the processor; wherein the memory stores instructions executable by the processor, and the instructions are executed by the processor so that the processor can execute the method according to any one of claims 1-7.

10. A computer-readable storage medium storing one or more programs, characterized in that: The one or more programs include computer instructions for causing a computer to execute the method according to any one of claims 1-7.

Images