Elastic parameter mapping method, electronic device, storage medium and apparatus

Abstract

The application discloses an elastic parameter mapping method, an electronic device, a storage medium and an apparatus. The method comprises the following steps: establishing a seismic forward model, wherein the seismic forward model is a regular grid model; obtaining a spatial parameter of the seismic forward model; extracting a spatial boundary of the seismic forward model; performing irregular grid dissection on the seismic forward model based on the spatial boundary; extending an irregular grid boundary of the seismic forward model; obtaining a mapping relationship between regular grid nodes and irregular grid nodes of the seismic forward model based on the spatial parameter; and assigning elastic parameters on all the regular grid nodes to corresponding irregular grid nodes according to the mapping relationship, so as to complete the mapping of the elastic parameters from the regular grid to the irregular grid. The application solves the problems of unstable irregular grid dissection and low dissection grid quality caused by the dramatic change of the elastic parameters, and ensures the simulation accuracy of the seismic wave field forward simulation from the spatial grid aspect.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of application geophysical seismic wave field forward modeling, and more particularly, relates to an elastic parameter mapping method, an electronic device, a storage medium and an apparatus. BACKGROUND

[0002] Seismic wave field numerical forward modeling technology is one of the main technologies of seismic exploration, is based on wave equation, uses computer numerical calculation method to simulate the propagation process of seismic wave in underground medium (forward model), and is an important means to study the propagation law of seismic wave. With the continuous deepening of oil and gas exploration and development, the geological problems encountered are more and more complex. For the seismic wave field numerical forward modeling of complex geological forward model, a numerical method with high simulation efficiency and high simulation precision becomes particularly important. At present, the main seismic wave field numerical forward modeling methods include finite difference classification and finite element method, etc. Among them, the finite element method has higher simulation precision than the finite difference classification method, and is more suitable for seismic wave field numerical forward modeling work under complex geological conditions.

[0003] Generally, for the wave equation (hyperbolic equation), the following parts need to be completed to realize finite element numerical forward modeling calculation, including determining the solution region (i.e. seismic forward model), giving the initial value condition and boundary condition of the wave equation; constructing the weak form of the wave equation by using the virtual work principle; equation discretization and solution domain (i.e. seismic forward model) discretization; unit analysis, forming the unit stiffness matrix; synthesizing the overall stiffness matrix; processing boundary conditions; solving large linear equations. The solution of the large linear equations is completed to solve the wave equation, and then the propagation process of the simulated seismic wave field in the underground medium is realized. Among them, for the discretization of the seismic forward model, in order to accurately describe the complex (severe surface relief, severe lateral variation of underground medium, etc.) forward model, the finite element method generally selects the non-regular grid discretization method.

[0004] Before seismic wave field numerical forward modeling, the modeling of seismic forward elastic parameter model (referred to as seismic forward model) needs to be completed first. The established seismic forward model is usually stored in the form of regular grid. Therefore, before carrying out finite element seismic wave field numerical forward modeling, the regular grid format seismic forward model (i.e. elastic parameter model) needs to be mapped to the non-regular grid format seismic forward model. For the mapping of elastic parameters from regular grid format to non-regular grid format, the existing technology usually performs non-regular grid discretization on the forward model space according to the elastic parameter value of the seismic forward model as a constraint condition, and completes the mapping of the elastic parameters of the forward model.

[0005] Due to the fact that the elastic parameters of the seismic forward model under complex geological conditions change dramatically, the space grid of the forward model is divided based on the elastic parameter value as a constraint condition, in the area where the elastic parameters change dramatically, the scale of the irregular grid changes too much in the transition space, which easily causes the grid division to be too large, thus the grid division process is unstable and the quality of the divided grid is not high. The elastic parameters are stored by using the low-quality irregular grid, which easily causes the process of the finite element type seismic wave field forward simulation to be unstable and numerical dispersion to be generated, and to some extent, the numerical forward simulation precision of the finite element type method is affected.

[0006] The information disclosed in the Background section of the present application is only intended to deepen the understanding of the general background of the present application, and should not be regarded as acknowledging or implying in any form that the information constitutes prior art known to those skilled in the art. SUMMARY

[0007] The purpose of the present application is to provide an elastic parameter mapping method, electronic equipment, storage medium and device, which realizes the mapping of the elastic parameters of the regular grid to the irregular grid, solves the problems of unstable irregular grid division and low quality of the divided grid caused by the dramatic change of the elastic parameters, and guarantees the simulation precision of the seismic wave field forward simulation from the aspect of the space grid.

[0008] To achieve the above-mentioned purpose, the present application provides an elastic parameter mapping method, electronic equipment, storage medium and device.

[0009] According to the first aspect of the present application, an elastic parameter mapping method is provided, comprising:

[0010] establishing a seismic forward model, wherein the seismic forward model is a regular grid model;

[0011] obtaining the space parameters of the seismic forward model;

[0012] extracting the space boundary of the seismic forward model;

[0013] dividing the seismic forward model into irregular grids based on the space boundary;

[0014] extending the irregular grid boundary of the seismic forward model;

[0015] obtaining the mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model based on the space parameters;

[0016] assigning the elastic parameters on all the regular grid nodes to the corresponding irregular grid nodes according to the mapping relationship, and completing the mapping of the elastic parameters from the regular grid to the irregular grid.

[0017] Optionally, the spatial parameters of the seismic forward model are directly extracted according to trace head information of the seismic forward model.

[0018] Alternatively, the spatial parameters of the seismic forward model are calculated according to spatial distance and spatial sampling points of the seismic forward model.

[0019] Optionally, the spatial boundary of the seismic forward model is extracted by extracting relief ground coordinate points of the seismic forward model.

[0020] Optionally, the extraction of the relief ground coordinate points of the seismic forward model comprises:

[0021] All horizontal coordinate points of the seismic forward model are traversed;

[0022] All depth coordinate points of the seismic forward model are traversed at each horizontal coordinate point to obtain depth coordinate points filled with elastic parameters;

[0023] All depth coordinate points filled with elastic parameters and corresponding horizontal coordinate points form the relief ground coordinate.

[0024] Optionally, the seismic forward model is subjected to irregular grid division by using a spatial boundary constraint method.

[0025] Optionally, the extension of the irregular grid boundary of the seismic forward model comprises:

[0026] The irregular grid of the seismic forward model is traversed to obtain an irregular boundary grid of the seismic forward model;

[0027] The irregular boundary grid is divided into six groups of boundary grid units in different boundary positions, i.e., upper, lower, left, right, front and back;

[0028] Boundary grid nodes are defined in the boundary grid units, and the boundary grid units are extended outward according to the boundary grid nodes at a set extension grid step until the extension thickness of the boundary grid units meets the requirements.

[0029] Optionally, the mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model is obtained based on the spatial parameters comprises:

[0030] All irregular grid nodes of the seismic forward model are traversed to obtain horizontal coordinate values and depth coordinate values of all the irregular grid nodes;

[0031] The coordinate values of all the irregular grid nodes and the spatial step are calculated according to the spatial step of the spatial parameters, the horizontal coordinate values and the depth coordinate values;

[0032] For each of the multiples, the multiple is multiplied by the space step to obtain a coordinate value of the regular grid point closest to each of the irregular grid points;

[0033] The coordinate value of each of the regular grid points is stored into the sequence of the irregular grid points, and the nearest mapping of the regular grid nodes to the irregular grid nodes is completed.

[0034] According to a second aspect of the present application, an elastic parameter mapping device is provided, comprising:

[0035] A establishing module is configured to establish a seismic forward model, the seismic forward model being a regular grid model;

[0036] A first obtaining module is configured to obtain a spatial parameter of the seismic forward model;

[0037] An extracting module is configured to extract a spatial boundary of the seismic forward model;

[0038] A grid division module is configured to perform irregular grid division on the seismic forward model based on the spatial boundary;

[0039] An extending module is configured to extend an irregular grid boundary of the seismic forward model;

[0040] A second obtaining module is configured to obtain a mapping relationship between regular grid nodes and irregular grid nodes of the seismic forward model based on the spatial parameter;

[0041] A mapping module is configured to assign elastic parameters on all the regular grid nodes to corresponding irregular grid nodes according to the mapping relationship, and complete mapping of the elastic parameters from the regular grid to the irregular grid.

[0042] According to a third aspect of the present application, an electronic device is provided, comprising:

[0043] At least one processor; and,

[0044] A memory connected with the at least one processor in communication; wherein,

[0045] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the elastic parameter mapping method according to any one of the first aspect.

[0046] According to a fourth aspect of the present application, a non-transitory computer readable storage medium is provided, the non-transitory computer readable storage medium storing computer instructions for causing a computer to execute the elastic parameter mapping method according to any one of the first aspect.

[0047] The present application has the beneficial effects that: the present application is aimed at the problems of unstable irregular grid division process and low quality of divided grid, performs irregular grid division on the seismic forward model, avoids the problems of low quality of grid division caused by elastic parameter change, obtains the mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model based on the space parameter, maps the elastic parameter value on the regular grid node to the irregular grid node based on the mapping relationship, thereby completing the mapping of the elastic parameter of the seismic forward model from the regular grid to the irregular grid, and solves the problems of unstable irregular grid division and low quality of divided grid caused by the dramatic change of the elastic parameter, thereby ensuring the simulation accuracy of the seismic wave field forward simulation from the space grid aspect.

[0048] The system of the present application has other characteristics and advantages, which will be apparent from or set forth in the accompanying drawings and the detailed description that follows, which together serve to explain certain principles of the present application. BRIEF DESCRIPTION OF DRAWINGS

[0049] The above and other objects, features and advantages of the present application will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings in which like reference characters refer to like parts throughout the several views.

[0050] Figure 1 A flow chart showing the steps of an elastic parameter mapping method according to the present application is shown.

[0051] Figure 2 A flow chart showing the steps of an elastic parameter mapping method according to embodiment 2 of the present application is shown.

[0052] Figure 3 A schematic diagram showing an elastic parameter mapping method according to embodiment 2 of the present application is shown.

[0053] Figure 4 A schematic diagram showing the mapping of the elastic parameter of the relief surface from the regular grid to the irregular grid according to embodiment 2 of the present application is shown.

[0054] Figure 5 A schematic diagram showing the mapping of the elastic parameter of the underground geological body from the regular grid to the irregular grid according to embodiment 2 of the present application is shown.

[0055] Figure 6 A schematic diagram showing an elastic parameter mapping device according to embodiment 3 of the present application is shown. DETAILED DESCRIPTION

[0056] The application will be described in more detail with reference to the drawings. Although the preferred embodiments of the application are shown in the drawings, it is understood that the application can be carried out in various forms and should not be limited by the embodiments set forth herein. On the contrary, these embodiments are provided so that this application will be thorough and complete, and will fully convey the scope of the application to those skilled in the art.

[0057] As shown in the drawings, according to an elastic parameter mapping method of the application, comprising: Figure 1

[0058] establishing a seismic forward model, the seismic forward model being a regular grid model;

[0059] acquiring spatial parameters of the seismic forward model;

[0060] extracting spatial boundaries of the seismic forward model;

[0061] performing irregular grid dissection on the seismic forward model based on the spatial boundaries;

[0062] extending irregular grid boundaries of the seismic forward model;

[0063] acquiring a mapping relationship between regular grid nodes and irregular grid nodes of the seismic forward model based on the spatial parameters;

[0064] assigning elastic parameters on all regular grid nodes to corresponding irregular grid nodes according to the mapping relationship, to complete the mapping of the elastic parameters from the regular grid to the irregular grid.

[0065] ​Specifically, the present application firstly establishes a regular grid seismic forward model, and directly extracts spatial parameters such as spatial step length of the forward model according to trace information of the seismic forward model, or calculates the spatial parameters such as the spatial step length of the forward model according to spatial distance and spatial sampling point number of the seismic forward model; then extracts spatial boundaries of the seismic forward model to depict spatial profile of the model. Due to the particularity of the seismic forward model, in addition to the model surface relief, other several spatial surfaces are in vertical orthogonal relationship, so the extraction of the spatial boundaries of the seismic forward model is completed by completing the extraction of the coordinate parameters of the relief surface; then the seismic forward model is divided into irregular grids by means of existing spatial boundary constraint methods, for example, the seismic forward model is divided into irregular grids by using a non-elastic parameter constraint method. In order to meet the requirement of accurate depiction of the relief surface by the irregular grid and the variable-scale grid description of the model at different depths, the present application adopts the downward continuation along the surface, and realizes the purpose of changing the grid size with the increase of the depth by setting the surface distance parameter according to the specific model requirement; the irregular grid boundary of the seismic forward model is expanded. The seismic wave numerical forward simulation generally simulates the wave field propagation process in a limited space, so the absorption boundary of a certain thickness must be expanded at the truncated boundary of the limited space model to eliminate the false reflection caused by the model truncation; the mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model is obtained based on the spatial parameters and the nearest mapping principle, and the elastic parameters on all the regular grid nodes are assigned to the corresponding irregular grid nodes according to the mapping relationship, so as to complete the mapping of the elastic parameters from the regular grid to the irregular grid. The present application adopts the method of model space overall grid division to replace the elastic parameter constraint to divide the seismic forward model into irregular grids, and maps the regular grid elastic parameters into the irregular grid according to the nearest mapping principle, so as to solve the problems of unstable irregular grid division and low grid quality caused by the dramatic change of the elastic parameters, and ensure the simulation accuracy of the seismic wave field forward simulation from the spatial grid aspect.

[0066] In one example, the spatial parameters of the seismic forward model are directly extracted according to trace information of the seismic forward model.

[0067] Alternatively, the spatial parameters of the seismic forward model are calculated according to spatial distance and spatial sampling point number of the seismic forward model.

[0068] In one example, the spatial boundaries of the seismic forward model are extracted by extracting relief surface coordinate points of the seismic forward model.

[0069] In one example, the relief surface coordinate points of the seismic forward model include:

[0070] All horizontal coordinate points of the seismic forward model are traversed;

[0071] traversing all depth coordinate points of the seismic forward model at each horizontal coordinate point, obtaining the depth coordinate points filled with the elastic parameters;

[0072] all depth coordinate points filled with the elastic parameters and corresponding horizontal coordinate points form the relief surface coordinates.

[0073] Specifically, due to the particularity of the seismic forward model, in addition to the model surface relief, the other several spatial surfaces are vertically orthogonal, so the extraction of the spatial boundary of the seismic forward model is completed by extracting the model relief surface coordinates. The present application traverses all horizontal coordinate points of the regular grid model in a traversal model space manner, traverses all depth coordinate points at each horizontal coordinate point, obtains the depth coordinates of the relief surface by judging the zero point boundary position of the elastic parameter value (whether the model space is filled with the elastic parameters), and thus completes the coordinate parameter extraction of the relief surface, thereby completing the extraction of the spatial boundary of the seismic forward model.

[0074] In one example, the spatial boundary constraint method is used to perform irregular grid dissection on the seismic forward model.

[0075] Specifically, the spatial boundary constraint method includes a non-elastic parameter constraint method.

[0076] In one example, the irregular grid boundary of the extended seismic forward model includes:

[0077] traversing the irregular grid of the seismic forward model to obtain the irregular boundary grid of the seismic forward model;

[0078] dividing the irregular boundary grid into six groups of boundary grid units with different boundary positions of upper, lower, left, right, front and back;

[0079] defining boundary grid nodes in the boundary grid unit, and expanding the boundary grid unit outward according to the boundary grid nodes at a set expansion grid step length until the expansion thickness of the boundary grid unit meets the requirements.

[0080] Specifically, the irregular grid absorbs the boundary expansion. Seismic wave numerical forward modeling generally simulates the wave field propagation process in a limited space, so the truncated boundary of the limited space model must be expanded to a certain thickness to absorb the boundary, so as to eliminate the false reflection caused by model truncation. The expansion of the absorbing boundary first traverses all irregular grids of the model, defines the boundary grid unit of the model, and divides the boundary unit grid into six groups of grid units with different boundary positions of upper, lower, left, right, front and back; secondly, the boundary grid nodes are defined on the boundary grid unit; thirdly, the unit grid is expanded outward according to the defined grid boundary nodes (the expansion grid step length can be automatically calculated from the average step length of the irregular grid, or can be input as a parameter); and finally, the above three steps are repeated until the expansion boundary thickness meets the set requirements.

[0081] In one example, the mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model based on the spatial parameters comprises:

[0082] All irregular grid nodes of the seismic forward model are traversed to obtain horizontal coordinate values and depth coordinate values of all irregular grid nodes;

[0083] The multiple of the coordinate values of all irregular grid nodes and the spatial step length is calculated according to the spatial step length, the horizontal coordinate values and the depth coordinate values of the spatial parameters;

[0084] The integer of each multiple is obtained, and the coordinate values of the nearest regular grid node to each irregular grid node are obtained by multiplying the integer and the spatial step length;

[0085] The coordinate values of each regular grid node are stored in the sequence of irregular grid nodes, and the nearest mapping of the regular grid nodes and the irregular grid nodes is completed.

[0086] Specifically, the mapping relationship between the regular grid coordinates and the irregular grid coordinates is defined according to the 'nearest' principle; all irregular grid node coordinates are traversed; the multiple relationship of the irregular grid coordinate values and the regular grid spatial step length is calculated according to the spatial step length parameter of the regular grid; the multiple relationship is rounded, and the nearest regular grid node coordinate to the irregular grid node is obtained by multiplying the rounded multiple relationship and the spatial step length; the nearest grid node coordinate information is stored in the irregular grid node sequence, and the 'nearest' mapping definition of the regular grid coordinates and the irregular grid coordinates is completed.

[0087] The application will be further described below in conjunction with the drawings and specific embodiments, but not as a limitation of the application. It should be noted that the embodiments in the application and the features in the embodiments can be combined with each other without conflict.

[0088] Embodiment 1

[0089] The embodiment provides an elastic parameter mapping method, comprising:

[0090] The seismic forward model is a regular grid model. Spatial parameters of the seismic forward model are obtained. In this embodiment, the spatial step length of the forward model and other spatial parameters are directly extracted according to the trace information of the seismic forward model, or the spatial step length of the forward model and other spatial parameters are calculated according to the spatial distance and the number of spatial sampling points of the seismic forward model. The spatial boundary of the seismic forward model is extracted. Due to the particularity of the seismic forward model, in addition to the model surface relief, the other several spatial surfaces are vertically orthogonal, so the spatial boundary of the seismic forward model is extracted by extracting the relief surface coordinate points of the seismic forward model. All horizontal coordinate points of the seismic forward model are traversed, and all depth coordinate points of the seismic forward model are traversed at each horizontal coordinate point to obtain depth coordinate points filled with elastic parameters. All depth coordinate points filled with elastic parameters and corresponding horizontal coordinate points form the relief surface coordinates. The spatial boundary constraint method is used to perform irregular grid division on the seismic forward model based on the spatial boundary. The irregular grid boundary of the seismic forward model is expanded. Seismic wave numerical forward modeling generally simulates the wave field propagation process in a limited space, so the absorption boundary of a certain thickness must be expanded at the truncated boundary of the limited space model to eliminate the false reflection caused by the model truncation. The irregular boundary grid of the seismic forward model is obtained by traversing the irregular grid of the model. The irregular boundary grid is divided into six groups of boundary grid units in different boundary positions of upper, lower, left, right, front and back. The boundary grid nodes are defined in the boundary grid units, and the boundary grid units are expanded outward according to the boundary grid nodes and the set expansion grid step length until the expansion thickness of the boundary grid units meets the requirements. The mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model is obtained based on the spatial parameters. The horizontal coordinate values and the depth coordinate values of all irregular grid nodes of the seismic forward model are obtained by traversing all irregular grid nodes of the seismic forward model. The coordinate values of all irregular grid nodes and the multiples of the spatial step length are calculated according to the spatial step length, the horizontal coordinate values and the depth coordinate values of the spatial parameters. The coordinate values of the nearest regular grid points to each irregular grid point are obtained by taking the integer of each multiple and multiplying the integer by the spatial step length. The coordinate values of each regular grid point are stored in the sequence of the irregular grid points to complete the nearest mapping of the regular grid nodes and the irregular grid nodes. The elastic parameters on all regular grid nodes are assigned to the corresponding irregular grid nodes according to the mapping relationship, and the mapping of the elastic parameters from the regular grid to the irregular grid is completed.

[0091] Embodiment 2

[0092] As shown in Figure 2 and Figure 3 , the embodiment provides an elastic parameter mapping method, which comprises:

[0093] The first step is to prepare a seismic forward model. The finite element method seismic wave field numerical forward modeling is based on a seismic forward model, so the first step of the embodiment is to prepare seismic forward model (elastic parameter) data, and to directly extract the forward model space step and other model space parameters according to the forward model trace information, or to calculate the forward model space step and other model space parameters according to the seismic forward model space distance and the space sampling point number.

[0094] The second step is to extract the spatial boundary of the seismic forward model. Since the embodiment adopts the method of spatially dividing the whole seismic forward model to generate a spatial irregular grid, after the first step of data preparation, the second step is to extract the spatial boundary of the seismic forward model to depict the spatial profile of the model. Due to the particularity of the seismic forward model, in addition to the model surface relief, the other several spatial surfaces are vertically orthogonal, so the technical point of this step is to extract the relief surface coordinates of the model. The first layer traverses all the horizontal coordinate points of the regular grid model, and the second layer traverses all the depth coordinates at each horizontal coordinate point. The depth coordinates of the relief surface are obtained by judging the zero point of the elastic parameter value (whether the model space is filled with elastic parameters), so as to complete the extraction of the relief surface coordinates and the extraction of the spatial boundary of the seismic forward model.

[0095] The third step is to divide the whole grid of the seismic forward model. The seismic forward model boundary extracted in the second step is used to constrain the non-regular grid division of the forward model by using the existing spatial boundary constraint (non-elastic parameter constraint) method. In order to meet the requirement of accurately depicting the relief surface by the non-regular grid and describing the variable scale grid at different depths of the model, the embodiment adopts the method of downward continuation along the surface, and sets the surface distance parameter to realize the purpose of changing the grid size with the increase of depth according to the specific model requirement.

[0096] The fourth step is to expand the absorbing boundary of the non-regular grid. Seismic wave numerical forward modeling generally simulates the wave field propagation process in a limited space, so the absorbing boundary of a certain thickness must be expanded at the truncated boundary of the limited space model to eliminate the false reflection caused by the model truncation. The absorbing boundary is expanded by first traversing all the non-regular grids of the model, defining the boundary grid elements of the model, and dividing the boundary element grid into six groups of grid elements with different boundary positions, i.e. upper, lower, left, right, front and back. Then, the boundary grid nodes are defined on the boundary grid elements. Then, the unit grid is expanded outward according to the defined grid boundary nodes (the expansion grid step can be automatically calculated from the average step of the non-regular grid, or can be input as a parameter). Finally, the above three steps are repeated until the expansion boundary thickness meets the set requirement.

[0097] Fifth, based on the "nearest" principle, define the mapping relationship between regular and irregular grid coordinates. Iterate through all irregular grid node coordinates; calculate the ratio between the irregular grid coordinate values ​​and the regular grid space step size obtained in the first step; round the ratio, and multiply the rounded ratio by the space step size to obtain the coordinates of the nearest regular grid point; store the nearest grid point coordinate information in the irregular grid node sequence, thus completing the "nearest" mapping definition between regular and irregular grid coordinates.

[0098] Step 6: Elastic parameter mapping. Based on the mesh mapping relationship obtained in Step 5, the elastic parameter at each regular mesh coordinate point is assigned to the corresponding irregular mesh point, completing the mapping of the elastic parameter from the regular mesh to the irregular mesh.

[0099] The seismic forward model used in this embodiment belongs to the geological conditions of the piedmont zone, which consists of a series of roughly parallel imbricate thrust zones. The foremountain fault and central fault form a large-scale, multi-level superimposed thrust-nappe tectonic zone with typical thrust structure characteristics. Above the thrust structure, two other structural styles have developed: typical imbricate structures and flyover structures, exhibiting a forward-extending reverse thrust-nappe tectonic zone. Its complex geological features result in particularly dramatic lateral variations in the elastic parameters of the corresponding seismic forward model. Simultaneously, the topographic gradient of the piedmont zone exhibits a steep change, ranging from 500m to 5000m within a 30-50km range, resulting in a highly undulating surface.

[0100] For seismic forward modeling of piedmont zones, due to the dramatic variations in lateral elastic parameters and the severe surface undulations, mesh generation using elastic parameter-constrained meshing methods easily produces low-quality mesh cells with large skewness. However, global mesh generation methods can easily generate high-quality meshes. Furthermore, by utilizing the "nearest" principle to map elastic parameters onto irregular meshes, the problem of low-quality irregular meshes is avoided. For example... Figure 4 As shown, using the global meshing method to partition the model space does not affect the accurate description of undulating surfaces by irregular meshes. Furthermore, irregular meshes are not constrained by elastic parameters and all exhibit near-orthogonal, high-quality meshes. For example... Figure 5 As shown, the elastic parameters of the geological body retain their variation characteristics when mapped from a regular grid to an irregular grid.

[0101] Example 3

[0102] like Figure 6 As shown, this embodiment provides an elastic parameter mapping device, including:

[0103] A module is created to build a seismic forward model, which is a regular grid model.

[0104] a first obtaining module, configured to obtain a spatial parameter of a seismic forward model;

[0105] an extracting module, configured to extract a spatial boundary of the seismic forward model;

[0106] a meshing module, configured to perform irregular meshing on the seismic forward model based on the spatial boundary;

[0107] an extending module, configured to extend an irregular mesh boundary of the seismic forward model;

[0108] a second obtaining module, configured to obtain a mapping relationship between regular mesh nodes and irregular mesh nodes of the seismic forward model based on the spatial parameter;

[0109] a mapping module, configured to assign elastic parameters on all regular mesh nodes to corresponding irregular mesh nodes according to the mapping relationship, so as to complete mapping of the elastic parameters from the regular mesh to the irregular mesh.

[0110] Embodiment 4

[0111] The embodiment provides an electronic device, which comprises:

[0112] at least one processor; and

[0113] a memory connected with the at least one processor in communication; wherein

[0114] the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to perform the elastic parameter mapping method in the embodiment 1.

[0115] The electronic device according to the embodiment of the present disclosure comprises a memory and a processor, the memory is used to store non-transitory computer readable instructions. Specifically, the memory can comprise one or more computer program products, which can comprise 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, etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.

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

[0117] Those skilled in the art shall understand that, in order to solve the technical problem of how to obtain a good user experience effect, the embodiment can also include well-known structures such as a communication bus, an interface, and the like, which shall also be included in the protection scope of the present disclosure.

[0118] The detailed description of the present embodiment can refer to the corresponding description in the foregoing embodiments, which will not be repeated here.

[0119] Embodiment 5

[0120] The present embodiment provides a non-transitory computer readable storage medium storing computer instructions for causing a computer to execute the elastic parameter mapping method in the embodiment 1.

[0121] The computer readable storage medium according to the embodiments of the present disclosure has non-transitory computer readable instructions stored thereon. When the non-transitory computer readable instructions are run by a processor, all or part of the steps of the method of the embodiments of the present disclosure are executed.

[0122] The computer readable storage medium described above includes, but is not limited to, an optical storage medium (for example, CD-ROM and DVD), a magneto-optical storage medium (for example, MO), a magnetic storage medium (for example, magnetic tape or a mobile hard disk), a medium with a built-in rewritable non-volatile memory (for example, a memory card), and a medium with a built-in ROM (for example, a ROM cartridge).

[0123] The above has described the embodiments of the present application, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes are obvious to those skilled in the art without departing from the scope and spirit of the described embodiments.

[0124] The above has described the embodiments of the present application, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes are obvious to those skilled in the art without departing from the scope and spirit of the described embodiments.

Claims

1. An elastic parameter mapping method, characterized by, The method comprises the following steps: establishing a seismic forward model, which is a regular grid model; obtaining spatial parameters of the seismic forward model; extracting the spatial boundary of the seismic forward model; performing irregular grid division on the seismic forward model based on the spatial boundary; extending the irregular grid boundary of the seismic forward model; obtaining a mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model based on the spatial parameters; assigning the elastic parameters on all the regular grid nodes to the corresponding irregular grid nodes according to the mapping relationship, and completing the mapping of the elastic parameters from the regular grid to the irregular grid.

2. The elastic parameter mapping method of claim 1, wherein, directly extracting the spatial parameters of the seismic forward model according to the trace information of the seismic forward model; alternatively, calculating the spatial parameters of the seismic forward model according to the spatial distance and the spatial sampling points of the seismic forward model.

3. The method of claim 1, wherein extracting the spatial boundary of the seismic forward model by extracting the undulating ground surface coordinate points of the seismic forward model.

4. The elastic parameter mapping method of claim 3, wherein, The method of extracting the undulating ground surface coordinate points of the seismic forward model comprises the following steps: traversing all the horizontal coordinate points of the seismic forward model; traversing all the depth coordinate points of the seismic forward model at each horizontal coordinate point to obtain the depth coordinate points filled with elastic parameters; all the depth coordinate points filled with elastic parameters and the corresponding horizontal coordinate points constitute the undulating ground surface coordinates.

5. The method of claim 1, wherein performing irregular grid division on the seismic forward model by using a spatial boundary constraint method.

6. The method of claim 1, wherein The method of extending the irregular grid boundary of the seismic forward model comprises the following steps: traversing the irregular grid of the seismic forward model to obtain the irregular boundary grid of the seismic forward model; dividing the irregular boundary grid into six groups of boundary grid units with different boundary positions, i.e., upper, lower, left, right, front and back; defining a boundary grid node in the boundary grid unit, and expanding the boundary grid unit outward according to the boundary grid node at a set expansion grid step until the expansion thickness of the boundary grid unit meets the requirements.

7. The method of claim 1, wherein The method of obtaining a mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model based on the spatial parameters comprises the following steps: traversing all the irregular grid nodes of the seismic forward model to obtain the horizontal coordinate values and the depth coordinate values of all the irregular grid nodes; calculating the coordinate values of all the irregular grid nodes and the multiples of the spatial step according to the spatial step, the horizontal coordinate values and the depth coordinate values of the spatial parameters; taking the integer of each multiple, multiplying the integer multiple by the spatial step to obtain the coordinate values of the nearest regular grid node to each irregular grid node; storing the coordinate values of each regular grid node into the sequence of the irregular grid node to complete the nearest mapping between the regular grid nodes and the irregular grid nodes.

8. An electronic device, comprising: The electronic device comprises: at least one processor; and a memory connected in communication with the at least one processor; wherein The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to perform the elastic parameter mapping method of any one of claims 1-7.

9. A non-transitory computer-readable storage medium, comprising: The non-transitory computer readable storage medium stores computer instructions for causing a computer to perform the elastic parameter mapping method of any one of claims 1-7.

10. An elastic parameter mapping apparatus, characterized by Comprise: The establishing module is used for establishing a seismic forward model, and the seismic forward model is a regular grid model; The first obtaining module is used for obtaining a spatial parameter of the seismic forward model; The extracting module is used for extracting a spatial boundary of the seismic forward model; The grid division module is used for performing non-regular grid division on the seismic forward model based on the spatial boundary; The extending module is used for extending a non-regular grid boundary of the seismic forward model; The second obtaining module is used for obtaining a mapping relationship between a regular grid node and a non-regular grid node of the seismic forward model based on the spatial parameter; The mapping module is used for assigning elastic parameters on all the regular grid nodes to corresponding non-regular grid nodes according to the mapping relationship, and completing the mapping of the elastic parameters from the regular grid to the non-regular grid.

Images