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

Abstract

The invention discloses an elastic parameter mapping method and device, electronic equipment and a storage medium. The method comprises the steps that an earthquake forward modeling model is established, and the earthquake forward modeling model is a regular grid model; acquiring space parameters of the seismic forward modeling model; extracting a space boundary of the seismic forward modeling model; irregular mesh generation is carried out on the seismic forward modeling model based on the space boundary; expanding an irregular grid boundary of the seismic forward modeling model; obtaining a mapping relation between regular grid nodes and irregular grid nodes of the seismic forward modeling 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 relation, and completing the mapping of the elastic parameters from the regular grid to the irregular grid. According to the method, the problems of unstable irregular mesh generation, low mesh generation quality and the like caused by dramatic change of elastic parameters are solved, and the simulation precision of seismic wave field forward modeling simulation is ensured from the aspect of space meshes.

Description

Technical Field

[0001] The present invention belongs to the technical field of applied geophysical seismic wavefield forward simulation, and more specifically, relates to an elastic parameter mapping method, an electronic device, a storage medium, and a device. Background Art

[0002] The numerical forward simulation technology of seismic wavefield is one of the main technologies in seismic exploration. Based on the wave equation, it uses the method of computer numerical calculation to simulate the propagation process of seismic waves in underground media (forward model), and is an important means to study the propagation law of seismic waves. With the continuous deepening of oil and gas exploration and development, the geological problems encountered are becoming more and more complex. For the numerical forward simulation of seismic wavefield of complex geological forward models, numerical methods with high simulation efficiency and high simulation accuracy become particularly important. Currently, the main numerical forward simulation methods of seismic wavefield include various methods such as finite difference type and finite element type. Among them, the finite element type method has higher simulation accuracy than the finite difference type method and is more suitable for the numerical forward simulation of seismic wavefield under complex geological conditions.

[0003] Generally, for the finite element type numerical forward simulation calculation of the wave equation (hyperbolic equation), the following parts need to be completed, including determining the solution domain (i.e., the seismic forward model), giving the initial conditions and boundary conditions of the wave equation; constructing the weak form of the wave equation using the virtual work principle; discretizing the equation and the solution domain (i.e., the seismic forward model); element analysis to form the element stiffness matrix; synthesizing the global stiffness matrix; handling boundary conditions; and solving the large linear equations. By solving the large linear equations, the solution of the wave equation is completed, and then the propagation process of the simulated seismic wavefield in underground media is realized. Among them, for the discretization of the seismic forward model, in order to accurately describe complex (severe surface undulation, drastic lateral variation of underground media, etc.) forward models, the finite element method generally selects the non-regular grid discretization method.

[0004] Before the numerical forward simulation of seismic wavefield, it is first necessary to complete the modeling of the seismic forward elastic parameter model (abbreviated as the seismic forward model), and the established seismic forward model is usually stored in the form of a regular grid. Therefore, before the finite element type numerical forward simulation of seismic wavefield, it is necessary to map the regular grid format seismic forward model (i.e., the elastic parameter model) to the non-regular grid format seismic forward model. For the mapping of elastic parameters from the regular grid format to the non-regular grid format, the existing technology usually performs non-regular grid meshing on the forward model space based on the elastic parameter values of the seismic forward model as constraint conditions to complete the mapping of the elastic parameters of the forward model.

[0005] Due to the drastic changes in elastic parameters in the forward seismic model under complex geological conditions, when using the elastic parameter values as constraint conditions for spatial grid meshing of the forward model, in the area where the elastic parameters change drastically, the transition space of the irregular grid scale change is small, which easily causes excessive skewness of the grid meshing. Therefore, there are problems such as instability in the grid meshing process and low quality of the meshed grids. Storing elastic parameters using low-quality irregular grids easily leads to instability in the forward simulation process of the finite element type seismic wave field and easily generates numerical dispersion, which to a certain extent affects the numerical forward simulation accuracy of the finite element type method.

[0006] The information disclosed in the background art section of the present invention is only intended to deepen the understanding of the general background art of the present invention, and should not be regarded as an admission or any form of implication that this information constitutes the prior art known to those skilled in the art. Summary of the Invention

[0007] The object of the present invention is to propose an elastic parameter mapping method, an electronic device, a storage medium and a device, which realize mapping the elastic parameters of regular grids to irregular grids, solve problems such as unstable irregular grid meshing and low quality of meshed grids caused by drastic changes in elastic parameters, and ensure the simulation accuracy of the forward simulation of the seismic wave field from the aspect of spatial grids.

[0008] To achieve the above object, the present invention proposes an elastic parameter mapping method, an electronic device, a storage medium and a device.

[0009] According to the first aspect of the present invention, an elastic parameter mapping method is proposed, including:

[0010] Establish a forward seismic model, and the forward seismic model is a regular grid model;

[0011] Obtain the spatial parameters of the forward seismic model;

[0012] Extract the spatial boundary of the forward seismic model;

[0013] Perform irregular grid meshing on the forward seismic model based on the spatial boundary;

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

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

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

[0017] Optionally, directly extract the spatial parameters of the seismic forward model according to the header information of the seismic forward model;

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

[0019] Optionally, extract the spatial boundary of the seismic forward model by extracting the undulating surface coordinate points of the seismic forward model.

[0020] Optionally, the extraction of the undulating surface coordinate points of the seismic forward model includes:

[0021] Traverse all the horizontal coordinate points of the seismic forward model;

[0022] At each of the horizontal coordinate points, traverse all the depth coordinate points of the seismic forward model to obtain the depth coordinate points filled with elastic parameters;

[0023] All the depth coordinate points filled with elastic parameters and the corresponding horizontal coordinate points form the undulating surface coordinates.

[0024] Optionally, perform unstructured grid meshing on the seismic forward model by using a spatial boundary constraint method.

[0025] Optionally, the expansion of the unstructured grid boundary of the seismic forward model includes:

[0026] Traverse the unstructured grid of the seismic forward model to obtain the unstructured boundary grid of the seismic forward model;

[0027] Divide the unstructured boundary grid into six groups of boundary grid cells with different boundary positions of up, down, left, right, front, and back;

[0028] Define boundary grid nodes in the boundary grid cells, and expand the boundary grid cells outward according to the set expansion grid step length based on the boundary grid nodes until the expansion thickness of the boundary grid cells meets the requirements.

[0029] Optionally, the obtaining of the mapping relationship between the regular grid nodes and the unstructured grid nodes of the seismic forward model based on the spatial parameters includes:

[0030] Traverse all the unstructured grid nodes of the seismic forward model to obtain the horizontal coordinate values and depth coordinate values of all the unstructured grid nodes;

[0031] Calculate the multiples of the coordinate values of all the unstructured grid nodes and the spatial step length according to the spatial step length of the spatial parameters, the horizontal coordinate values, and the depth coordinate values;

[0032] Round each of the multiples, multiply the rounded multiples by the spatial step size to obtain the coordinate values of the regular grid points closest to each of the irregular grid points;

[0033] Store the coordinate values of each of the regular grid points into the sequence of the irregular grid points to complete the nearest mapping between the regular grid nodes and the irregular grid nodes.

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

[0035] A building module, configured to build a seismic forward model, where the seismic forward model is a regular grid model;

[0036] A first obtaining module, configured to obtain the spatial parameters of the seismic forward model;

[0037] An extraction module, configured to extract the spatial boundary of the seismic forward model;

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

[0039] An expansion module, configured to expand the irregular grid boundary of the seismic forward model;

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

[0041] A mapping module, configured to assign the elastic parameters on all the regular grid nodes to the 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.

[0042] According to a third aspect of the present invention, there is provided an electronic device, where the electronic device includes:

[0043] At least one processor; and,

[0044] A memory communicatively connected to the at least one processor; wherein,

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

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

[0047] The beneficial effects of the present invention are as follows: Aiming at the problems of unstable non - regular grid meshing process and low quality of meshed grids in the seismic forward modeling, the present invention conducts non - regular grid meshing on the seismic forward model, avoiding problems such as low grid meshing quality caused by changes in elastic parameters; based on spatial parameters, the mapping relationship between regular grid nodes and non - regular grid nodes of the seismic forward model is obtained, and based on this mapping relationship, the elastic parameter values at regular grid points are mapped to non - regular grid points, thereby completing the mapping of elastic parameters of the seismic forward model from regular grids to non - regular grids, solving problems such as unstable non - regular grid meshing and low quality of meshed grids caused by drastic changes in elastic parameters, and ensuring the simulation accuracy of seismic wave field forward simulation from the aspect of spatial grids.

[0048] The system of the present invention has other characteristics and advantages, which will be obvious from the accompanying drawings incorporated herein and the subsequent detailed description, or will be described in detail in the accompanying drawings incorporated herein and the subsequent detailed description. These accompanying drawings and detailed description are used together to explain the specific principles of the present invention. Description of the Drawings

[0049] By describing the exemplary embodiments of the present invention in more detail in conjunction with the accompanying drawings, the above - mentioned and other objects, features, and advantages of the present invention will become more obvious. In the exemplary embodiments of the present invention, the same reference numerals generally represent the same components.

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

[0051] Figure 2 A flowchart showing the steps of an elastic parameter mapping method according to Embodiment 2 of the present invention is shown.

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

[0053] Figure 4 A schematic diagram showing the mapping of elastic parameters of the undulating ground surface from regular grids to non - regular grids according to Embodiment 2 of the present invention is shown.

[0054] Figure 5 A schematic diagram showing the mapping of elastic parameters of underground geological bodies from regular grids to non - regular grids according to Embodiment 2 of the present invention is shown.

[0055] Figure 6 A schematic diagram showing an elastic parameter mapping device according to Embodiment 3 of the present invention is shown. Detailed Description of the Invention

[0056] The present invention will be described in more detail below with reference to the accompanying drawings. Although the preferred 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 will be more thorough and complete, and will fully convey the scope of the present invention to those skilled in the art.

[0057] As Figure 1 shown, an elastic parameter mapping method according to the present invention includes:

[0058] Establishing a seismic forward model, where the seismic forward model is a regular grid model;

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

[0060] Extracting the spatial boundary of the seismic forward model;

[0061] Performing unstructured grid meshing on the seismic forward model based on the spatial boundary;

[0062] Expanding the unstructured grid boundary of the seismic forward model;

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

[0064] Assigning the elastic parameters on all regular grid nodes to the corresponding unstructured grid nodes according to the mapping relationship, thereby completing the mapping of the elastic parameters from the regular grid to the unstructured grid.

[0065] Specifically, the present invention first establishes a seismic forward modeling of a regular grid, and directly extracts spatial parameters such as the spatial step length of the forward model according to the header information of the seismic forward model, or calculates spatial parameters such as the spatial step length of the forward model based on the spatial distance and the number of spatial sampling points of the seismic forward model; then extracts the spatial boundary of the seismic forward model to depict the spatial contour of the model. Due to the particularity of the seismic forward model, except for the surface undulation of the model, the other several spatial planes are in a vertically orthogonal relationship. Therefore, by completing the extraction of the coordinate parameters of the undulating surface, the extraction of the spatial boundary of the seismic forward model is completed; then, with the help of the existing spatial boundary constraint method, the seismic forward model is meshed into an irregular grid. For example, the inelastic parameter constraint method is used to mesh the seismic forward model into an irregular grid. To meet the requirements of accurately depicting the undulating surface by the irregular grid and the variable-scale grid description of different depths of the model, the present invention adopts the method of extending downward along the surface. According to the specific model requirements, the purpose of changing the grid size with the increase of depth is achieved by setting the surface distance parameter; extend the irregular grid boundary of the seismic forward model. The numerical forward simulation of seismic waves generally simulates the wave field propagation process in a finite space. Therefore, a certain thickness of absorption boundary must be extended at the truncated boundary of the finite space model to eliminate the false reflection caused by model truncation; obtain the mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model based on the spatial parameters and the nearest mapping principle, and assign the elastic parameters on all regular grid nodes to the 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. The present invention adopts the method of replacing the elastic parameter constraint with the overall grid meshing of the model space for the irregular grid meshing of the forward model space, and maps the regular grid elastic parameters to the irregular grid according to the nearest mapping principle, solving the problems of unstable irregular grid meshing and low mesh quality caused by the drastic change of elastic parameters, and ensuring the simulation accuracy of the seismic wave field forward simulation from the aspect of the spatial grid.

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

[0067] Or, the spatial parameters of the seismic forward model are calculated based on the spatial distance and the number of spatial sampling points of the seismic forward model.

[0068] In one example, the spatial boundary of the seismic forward model is extracted by extracting the coordinate points of the undulating surface of the seismic forward model.

[0069] In one example, extracting the coordinate points of the undulating surface of the seismic forward model includes:

[0070] Traverse all the horizontal coordinate points of the seismic forward model;

[0071] Traverse all depth coordinate points of the seismic forward model at each horizontal coordinate point to obtain depth coordinate points filled with elastic parameters;

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

[0073] Specifically, due to the particularity of the seismic forward model, except for the undulation of the model surface, the other several spatial planes are all in a vertically orthogonal relationship. Therefore, by extracting the undulating surface coordinates of the model, the spatial boundary of the seismic forward model is extracted. The present invention adopts a method of traversing the model space, traversing all horizontal coordinate points of the regular grid model, traversing all depth coordinate points at each horizontal coordinate point, and obtaining the depth coordinates of the undulating surface by judging the zero-point demarcation position of its elastic parameter value (whether there is elastic parameter filling in the model space), so as to complete the extraction of the coordinate parameters of the undulating surface, and thus complete the extraction of the spatial boundary of the seismic forward model.

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

[0075] Specifically, the spatial boundary constraint method includes an inelastic parameter constraint method.

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

[0077] Traverse the unstructured grid of the seismic forward model to obtain the unstructured boundary grid of the seismic forward model;

[0078] Divide the unstructured boundary grid into six groups of boundary grid cells with different boundary positions: upper, lower, left, right, front, and back;

[0079] Define boundary grid nodes in the boundary grid cells, and expand the boundary grid cells outward according to the defined grid boundary nodes at a set expansion grid step until the expansion thickness of the boundary grid cells meets the requirements.

[0080] Specifically, the expansion of the absorbing boundary of the unstructured grid. The numerical forward simulation of seismic waves generally simulates the wave field propagation process in a finite space. Therefore, at the truncated boundary of the finite space model, an absorbing boundary with a certain thickness needs to be expanded to eliminate the spurious reflections caused by model truncation. To expand the absorbing boundary, first traverse all unstructured grids of the model, define the boundary grid cells of the model, and divide the boundary cell grids into six groups of grid cells with different boundary positions: upper, lower, left, right, front, and back; secondly, define boundary grid nodes on the boundary grid cells; secondly, expand the cell grids outward according to the defined grid boundary nodes (the expansion grid step can be automatically calculated from the average step of the unstructured grid or can be input as a parameter); finally, loop the above three steps until the expansion boundary thickness meets the set requirements.

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

[0082] Traverse all the irregular grid nodes of the seismic forward modeling to obtain the horizontal coordinate values and depth coordinate values of all the irregular grid nodes;

[0083] Calculate the multiples of the coordinate values of all the irregular grid nodes and the spatial step according to the spatial step, horizontal coordinate value and depth coordinate value of the spatial parameters;

[0084] Round each multiple, and multiply the rounded multiple by the spatial step to obtain the coordinate values of the regular grid points closest to each irregular grid point;

[0085] Store the coordinate values of each regular grid point into the sequence of irregular grid points to complete the nearest mapping between the regular grid nodes and the irregular grid nodes.

[0086] Specifically, the present invention defines the mapping relationship between the regular grid coordinates and the irregular grid coordinates according to the "nearest" principle; traverses all the irregular grid node coordinates; calculates the multiple relationship between the irregular grid coordinate values and the regular grid spatial step according to the spatial step parameter of the regular grid; rounds the multiple relationship, and multiplies the rounded multiple relationship by the spatial step to obtain the regular grid point coordinates closest to the irregular grid point; stores the nearest grid point coordinate information into the sequence of irregular grid nodes, thereby completing the definition of the "nearest" mapping between the regular grid coordinates and the irregular grid coordinates.

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

[0088] Embodiment 1

[0089] This embodiment provides an elastic parameter mapping method, including:

[0090] Establish a forward seismic model, where the forward seismic model is a regular grid model; obtain the spatial parameters of the forward seismic model. In this embodiment, the spatial parameters such as the spatial step size of the forward model are directly extracted according to the header information of the forward seismic model, or the spatial step size of the forward model is calculated based on the spatial distance and the number of spatial sampling points of the forward seismic model; extract the spatial boundary of the forward seismic model. Due to the particularity of the forward seismic model, except for the undulating surface of the model, the other several spatial planes are in a vertically orthogonal relationship. Therefore, in this embodiment, the spatial boundary of the forward seismic model is extracted by extracting the coordinate points of the undulating surface of the forward seismic model. Traverse all the horizontal coordinate points of the forward seismic model, and at each horizontal coordinate point, traverse all the depth coordinate points of the forward seismic model 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 form the coordinate points of the undulating surface; use the spatial boundary constraint method to perform unstructured grid meshing on the forward seismic model based on the spatial boundary; expand the unstructured grid boundary of the forward seismic model. Generally, the numerical forward simulation of seismic waves simulates the wave field propagation process in a finite space. Therefore, at the truncated boundary of the finite space model, an absorbing boundary with a certain thickness must be expanded to eliminate the false reflection caused by model truncation. First, traverse all the unstructured grids of the model, traverse the unstructured grids of the forward seismic model to obtain the unstructured boundary grids of the forward seismic model; divide the unstructured boundary grids into six groups of boundary grid units with different boundary positions, namely up, down, left, right, front, and back; define boundary grid nodes in the boundary grid units, and expand the boundary grid units outward according to the set expansion grid step size based on the boundary grid nodes until the expansion thickness of the boundary grid units meets the requirements; obtain the mapping relationship between the regular grid nodes and the unstructured grid nodes of the forward seismic model based on the spatial parameters. Traverse all the unstructured grid nodes of the forward seismic model to obtain the horizontal coordinate values and depth coordinate values of all the unstructured grid nodes; calculate the multiple of the coordinate values of all the unstructured grid nodes and the spatial step size according to the spatial step size, horizontal coordinate values, and depth coordinate values of the spatial parameters; round each multiple, and multiply the rounded multiple by the spatial step size to obtain the coordinate values of the regular grid points closest to each unstructured grid point; store the coordinate values of each regular grid point in the sequence of unstructured grid points to complete the nearest mapping between the regular grid nodes and the unstructured grid nodes; assign the elastic parameters on all the regular grid nodes to the corresponding unstructured grid nodes according to the mapping relationship to complete the mapping of the elastic parameters from the regular grid to the unstructured grid.

[0091] Embodiment 2

[0092] As Figure 2 and Figure 3 shown, this embodiment provides an elastic parameter mapping method, including:

[0093] Step 1: Preparation of seismic forward model. The numerical forward simulation of seismic wave fields by the finite element method must be carried out based on the seismic forward model. Therefore, the first step in this embodiment is to prepare the data of the seismic forward model (elastic parameters), and directly extract the model space parameters such as the spatial step length of the forward model according to the header information of the forward model, or calculate the model space parameters such as the spatial step length of the forward model according to the spatial distance and the number of spatial sampling points of the seismic forward model.

[0094] Step 2: Extraction of the spatial boundary of the seismic forward model. Since the overall spatial dissection method of the seismic forward model is adopted in this embodiment to generate irregular spatial grids, after the data preparation is completed in the first step, the second step is to extract the spatial boundary of the seismic forward model and depict the spatial contour of the model. Due to the particularity of the seismic forward model, except for the surface undulation of the model, the other several spatial planes are in a vertically orthogonal relationship. Therefore, the technical point of this step is to extract the coordinates of the undulating surface of the model, and the extraction is carried out by traversing the model space. In the first layer, all horizontal coordinate points of the regular grid model are traversed. At each horizontal coordinate point, all depth coordinates are traversed in the second layer. By judging the zero-point demarcation position of its elastic parameter value (whether there is elastic parameter filling in the model space), the depth coordinates of the undulating surface are obtained, so as to complete the extraction of the coordinate parameters of the undulating surface, and thus complete the extraction of the spatial boundary of the seismic forward model.

[0095] Step 3: Overall grid dissection of the seismic forward model. Using the boundary of the seismic forward model extracted in the second step, the forward model is dissected into irregular grids by means of the existing spatial boundary constraint (non-elastic parameter constraint) method. To meet the requirements of accurately depicting the undulating surface by the irregular grid and the variable-scale grid description of different depths of the model, in this embodiment, the method of extending downward along the surface is adopted. According to the specific model requirements, the purpose of changing the grid size with the increase of depth is achieved by setting the surface distance parameter.

[0096] Step 4: Expansion of the absorbing boundary of the irregular grid. The numerical forward simulation of seismic waves generally simulates the wave field propagation process in a finite space. Therefore, a certain thickness of absorbing boundary must be extended at the truncated boundary of the finite space model to eliminate the false reflection caused by model truncation. To expand the absorbing boundary, first traverse all the irregular grids of the model, define the grid cells at the model boundary and divide the boundary cell grids into six groups of grid cells at different boundary positions: top, bottom, left, right, front, and back. Secondly, define the boundary grid nodes on the boundary grid cells. Secondly, expand the cell grids outward according to the defined grid boundary nodes (the expansion grid step length can be automatically calculated by the average step length of the irregular grid or can be input as a parameter). Finally, loop the above three steps until the expansion boundary thickness meets the set requirements.

[0097] The fifth step is to define the mapping relationship between regular grid coordinates and irregular grid coordinates according to the "nearest" principle. Traverse all irregular grid node coordinates; calculate the multiple relationship between irregular grid coordinate values ​​and regular grid spatial step length according to the regular grid spatial step length parameter obtained in the first step; round the multiple relationship, and obtain the regular grid point coordinate closest to the irregular grid point by multiplying the rounded multiple relationship with the spatial step length; store the nearest grid point coordinate information in the irregular grid node sequence, thereby completing the "nearest" mapping definition between regular grid coordinates and irregular grid coordinates.

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

[0099] The earthquake forward model adopted in this embodiment belongs to the geological conditions of the piedmont belt, and its geological conditions are composed of a series of roughly parallel imbricate thrust belts, with the front mountain fault and the central fault as the main sliding surface, forming a large-scale, multi-level superimposed thrust-nappe structural belt, which has typical thrust-nappe structural characteristics. On top of the thrust-nappe structural characteristics, two structural styles such as typical imbricate structures and Feilaifeng structures are also developed, showing a forward-spreading thrust-nappe structural belt. Its complex geological characteristics make the lateral changes of the elastic parameters of the corresponding earthquake forward model particularly violent. At the same time, the terrain gradient of the piedmont belt has a steep change feature, and the terrain gradient changes sharply from 500m to 5000m within the range of 30-50km, and its surface shows a violent undulating feature.

[0100] In the earthquake forward model with piedmont characteristics, due to the drastic changes in the lateral elastic parameters of the underground and the drastic surface undulations, the elastic parameter constraint grid generation method is prone to produce low-quality grid units with large skewness. The overall grid generation method can easily form a high-quality grid, and the elastic parameters are mapped to the irregular grid using the "nearest" principle, thereby avoiding the problem of low-quality irregular grids. Figure 4 As shown in the figure, the use of the overall meshing method to divide the model space does not affect the accurate description of the irregular grid on the undulating surface. At the same time, the irregular grid is not constrained by elastic parameters and is a high-quality grid that is close to orthogonal. Figure 5 As shown, the elastic parameters of the geological body maintain 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] Establishing a module for establishing a seismic forward model, which is a regular grid model;

[0104] A first acquisition module, configured to acquire spatial parameters of a seismic forward modeling model;

[0105] An extraction module, configured to extract spatial boundaries of a seismic forward modeling model;

[0106] A grid meshing module, configured to perform unstructured grid meshing on a seismic forward modeling model based on the spatial boundaries;

[0107] An extension module, configured to extend unstructured grid boundaries of a seismic forward modeling model;

[0108] A second acquisition module, configured to acquire a mapping relationship between regular grid nodes and unstructured grid nodes of a seismic forward modeling model based on the spatial parameters;

[0109] A mapping module, configured to assign elastic parameters on all regular grid nodes to corresponding unstructured grid nodes according to the mapping relationship, so as to complete the mapping of elastic parameters from regular grids to unstructured grids.

[0110] Embodiment 4

[0111] This embodiment provides an electronic device, which includes:

[0112] At least one processor; and,

[0113] A memory communicatively connected to the at least one processor; wherein,

[0114] The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor, so that the at least one processor can execute the elastic parameter mapping method in Embodiment 1.

[0115] The electronic device according to an embodiment of the present disclosure includes a memory and a processor, and the memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, and the computer program products may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, random access memory (RAM) and / or cache memory, etc. The non-volatile memory may include, for example, read-only memory (ROM), hard disk, flash memory, etc.

[0116] The processor may be a central processing unit (CPU) or other forms of processing units with data processing capabilities and / or instruction execution capabilities, and may 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 should understand that, in order to solve the technical problem of how to obtain good user experience effects, the present embodiment may also include well-known structures such as communication buses and interfaces, and these well-known structures should also be included in the protection scope of the present disclosure.

[0118] For the detailed description of this embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, and details will not be repeated here.

[0119] Embodiment 5

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

[0121] According to the computer-readable storage medium of an embodiment of the present disclosure, non-temporary computer-readable instructions are stored thereon. When the non-temporary computer-readable instructions are run by a processor, all or part of the steps of the methods of the foregoing embodiments of the present disclosure are executed.

[0122] The above-mentioned computer-readable storage medium includes but is not limited to: optical storage media (such as CD-ROMs and DVDs), magneto-optical storage media (such as MOs), magnetic storage media (such as magnetic tapes or external hard drives), media with built-in rewritable non-volatile memories (such as memory cards), and media with built-in ROMs (such as ROM cartridges).

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

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

Claims

1. An elastic parameter mapping method, characterized in that, it includes: Establish a forward seismic model, and the forward seismic model is a regular grid model; Obtain the spatial parameters of the forward seismic model; Extract the spatial boundary of the forward seismic model; Perform unstructured grid meshing on the forward seismic model based on the spatial boundary; Expand the unstructured grid boundary of the forward seismic model; Obtain the mapping relationship between the regular grid nodes and the unstructured grid nodes of the forward seismic model based on the spatial parameters; Assign the elastic parameters on all the regular grid nodes to the corresponding unstructured grid nodes according to the mapping relationship, and complete the mapping of the elastic parameters from the regular grid to the unstructured grid.

2. The elastic parameter mapping method according to claim 1, characterized in that, directly extract the spatial parameters of the forward seismic model according to the header information of the forward seismic model; or calculate the spatial parameters of the forward seismic model according to the spatial distance and the number of spatial sampling points of the forward seismic model.

3. The elastic parameter mapping method according to claim 1, characterized in that, extract the spatial boundary of the forward seismic model by extracting the undulating surface coordinate points of the forward seismic model.

4. The elastic parameter mapping method according to claim 4, characterized in that, the extraction of the undulating surface coordinate points of the forward seismic model includes: traverse all the horizontal coordinate points of the forward seismic model; traverse all the depth coordinate points of the forward seismic model at each horizontal coordinate point, and obtain the depth coordinate points filled with elastic parameters; all the depth coordinate points filled with elastic parameters and the corresponding horizontal coordinate points form the undulating surface coordinates.

5. The elastic parameter mapping method according to claim 1, characterized in that, perform unstructured grid meshing on the forward seismic model by using a spatial boundary constraint method.

6. The elastic parameter mapping method according to claim 1, characterized in that, the expansion of the unstructured grid boundary of the forward seismic model includes: traverse the unstructured grid of the forward seismic model, and obtain the unstructured boundary grid of the forward seismic model; divide the unstructured boundary grid into six groups of boundary grid cells with different boundary positions of up, down, left, right, front and back; define boundary grid nodes in the boundary grid cells, and expand the boundary grid cells outward according to the set expansion grid step length according to the boundary grid nodes until the expansion thickness of the boundary grid cells meets the requirements.

7. The elastic parameter mapping method according to claim 1, characterized in that, the obtaining of the mapping relationship between the regular grid nodes and the unstructured grid nodes of the forward seismic model based on the spatial parameters includes: traverse all the unstructured grid nodes of the forward seismic model, and obtain the horizontal coordinate values and depth coordinate values of all the unstructured grid nodes; calculate the multiples of the coordinate values of all the unstructured grid nodes and the spatial step length according to the spatial step length of the spatial parameters, the horizontal coordinate values and the depth coordinate values; Round each of the multiples, multiply the rounded multiple by the spatial step size to obtain the coordinate value of the regular grid point closest to each of the irregular grid points; Store the coordinate value of each regular grid point into the sequence of the irregular grid points to complete the nearest mapping between the regular grid nodes and the irregular grid nodes.

8. An electronic device, characterized in that, the electronic device includes: at least one processor; and, a memory communicatively connected to 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 execute the elastic parameter mapping method according to any one of claims 1-7.

9. A non-transitory computer-readable storage medium, characterized in that, the non-transitory computer-readable storage medium stores computer instructions for causing a computer to execute the elastic parameter mapping method according to any one of claims 1-7.

10. An elastic parameter mapping device, characterized in that, comprising: a building module for building a seismic forward model, the seismic forward model being a regular grid model; a first obtaining module for obtaining the spatial parameters of the seismic forward model; an extracting module for extracting the spatial boundary of the seismic forward model; a grid meshing module for performing irregular grid meshing on the seismic forward model based on the spatial boundary; an expanding module for expanding the irregular grid boundary of the seismic forward model; a second obtaining module for obtaining the mapping relationship between the regular grid nodes and the irregular grid nodes of the seismic forward model based on the spatial parameters; a mapping module for assigning the elastic parameters on all the regular grid nodes to the 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.

Images