CN114491897A - Seismic wave numerical simulation method, device, medium and electronic equipment - Google Patents
Seismic wave numerical simulation method, device, medium and electronic equipment Download PDFInfo
- Publication number
- CN114491897A CN114491897A CN202011150064.1A CN202011150064A CN114491897A CN 114491897 A CN114491897 A CN 114491897A CN 202011150064 A CN202011150064 A CN 202011150064A CN 114491897 A CN114491897 A CN 114491897A
- Authority
- CN
- China
- Prior art keywords
- filling curve
- grid
- centroid
- unit
- hilbert
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
- G06T17/205—Re-meshing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a seismic wave numerical simulation method, a device, a medium and electronic equipment, wherein the method comprises the following steps: generating a full model unstructured grid; acquiring the centroid of the grid unit in the corresponding calculation range of the current simulated cannon from the full-model unstructured grid, and establishing a centroid unit; filling a space formed by grid units in a calculation range by using a Hilbert space filling curve, and establishing a corresponding relation between a centroid unit and a Hilbert space filling curve sequence; deleting the empty serial number in the Hilbert space filling curve sequence; dividing an area formed by grid units in a calculation range; and performing seismic wave numerical simulation calculation based on the divided areas. The invention divides the grids into different calculation areas by using the area division based on the space filling curve, keeps the adjacent characteristics between the grids, ensures that the grid units of each subarea are similar in number, improves the division efficiency of the non-structural grids, reduces the calculation consumption and has high seismic wave numerical simulation calculation efficiency.
Description
Technical Field
The invention belongs to the technical field of petroleum exploration, and particularly relates to a seismic wave numerical simulation method, a device, a medium and electronic equipment based on space filling curve region division.
Background
The seismic numerical simulation is a reliable method for researching the propagation characteristics of seismic waves and analyzing geological structures. The seismic numerical simulation method mainly comprises a finite difference method, a finite element method, a pseudo-spectrum method and the like. Both of these methods rely on mesh generation of the model. The finite element method has the characteristics of higher precision and stronger applicability. Finite element meshes in seismic wave numerical simulation are typically unstructured tetrahedral meshes.
Seismic numerical simulation is faced with very large numerical models. In the field of seismic exploration, the entire model typically ranges from hundreds of cubic kilometers, with the number of grids of finite elements on the scale of tens of millions to hundreds of millions. Therefore, a parallel computing method is generally used for simulation. The most common parallel computing strategy at present is to divide computing tasks equally into different computing units, such as different CPU cores of a common computer, or different network nodes of a super computer. Each computational unit will be assigned a fixed area of simulation tasks in the model. The different computing units can communicate the computing results with each other at certain time nodes. The parallel task is premised on the region division of the model.
There are currently two main implementations for region partitioning of models: firstly, the model is artificially divided into regions, and then grids are generated for each subregion; and secondly, mesh generation is carried out on the full model, and then region division is carried out on the meshes. The second method controls the size of the sub-regions divided by the region through an algorithm, so that the number of grids of each sub-region is equivalent, and parallel calculation is facilitated.
For the finite element method, because it adopts the non-structural grid, the grid ordering is generally not regular, and the distribution of the number of grid units in the model is directly related to the velocity size distribution and the geometric form of the model. Some open source packages are available to implement the area division function of the unstructured grid. For example, METIS uses a multi-level recursive bisection method, a multi-level K-way bisection method and a multi-constraint partitioning mechanism. However, when facing the simulation task of the seismic model, the software needs to consume huge computing resources, and the efficiency of region division is low.
In the field of seismic exploration, each simulation requires a portion of the model to be cut out for the purpose when the surface observation system is in a shot-rolling mode. Therefore, before the simulation algorithm is operated, the area division needs to be carried out on the currently cut grid every time, and when the three-dimensional multi-shot calculation is carried out, the process is more time-consuming and resource-consuming, and the calculation efficiency is lower.
Therefore, it is necessary to develop a seismic wave numerical simulation method with low computational consumption and high computational efficiency.
Disclosure of Invention
The invention aims to provide a seismic wave numerical simulation method which saves calculation consumption and has high calculation efficiency.
In a first aspect, the present invention provides a seismic wave numerical simulation method based on space filling curve region division, including: generating a full model unstructured grid; acquiring the centroid of the grid unit in the corresponding calculation range of the current simulated cannon from the full-model unstructured grid, and establishing a centroid unit; filling a space formed by grid units in the calculation range by using a Hilbert space filling curve, and establishing a corresponding relation between a centroid unit and a Hilbert space filling curve sequence; based on the corresponding relation between the centroid unit and the Hilbert space filling curve sequence, deleting the empty serial number in the Hilbert space filling curve sequence; dividing the region formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the empty serial numbers deleted; and performing seismic wave numerical simulation calculation based on the divided areas.
Optionally, obtaining a centroid of a grid unit in a calculation range corresponding to the current simulated cannon in the full-model unstructured grid, and establishing the centroid unit includes: acquiring grid units in a calculation range corresponding to the current simulated cannon from the full-model unstructured grid according to the simulation requirements of the observation system; and determining the mass center of each grid unit, numbering the grid units and establishing a mass center unit.
Optionally, the establishing a corresponding relationship between the centroid unit and the hilbert space filling curve sequence includes: establishing a corresponding relation between the centroid unit and the Hilbert space filling curve sequence based on the serial number of the centroid unit contained in each Hilbert unit; wherein each Hilbert cell of the Hilbert space-filling curve comprises at most one centroid cell.
Optionally, the dividing, based on the hilbert space filling curve sequence from which the empty sequence numbers are deleted, the region formed by the grid units in the calculation range includes: and according to the preset number of region division, dividing the Hilbert space filling curve sequence with the null sequence number deleted, so as to divide the region formed by the grid units in the corresponding calculation range.
Optionally, the full-model unstructured mesh is a triangular unit in a two-dimensional space and a tetrahedral unit in a three-dimensional space.
In a second aspect, the present invention also provides an electronic device, including: a memory storing executable instructions; and the processor runs the executable instructions in the memory to realize the seismic wave numerical simulation method based on space filling curve area division.
In a third aspect, the present invention further provides a computer-readable storage medium, where a computer program is stored, and when the computer program is executed by a processor, the method for seismic wave numerical simulation based on space filling curve area division is implemented.
In a fourth aspect, the present invention further provides a seismic wave numerical simulation apparatus based on space filling curve region division, including: the full model unstructured grid generation module generates a full model unstructured grid; the centroid unit establishing module is used for acquiring the centroid of the grid unit in the corresponding calculation range of the current simulated cannon in the full-model non-structural grid and establishing a centroid unit; the corresponding relation establishing module is used for filling a space formed by grid units in the calculation range by using a Hilbert space filling curve and establishing a corresponding relation between a centroid unit and a Hilbert space filling curve sequence; a null sequence number deleting module, which deletes the null sequence number in the Hilbert space filling curve sequence based on the corresponding relation between the centroid unit and the Hilbert space filling curve sequence; the dividing module is used for dividing the region formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the empty serial numbers deleted; and the calculation module is used for performing seismic wave numerical simulation calculation based on the divided areas.
Optionally, obtaining a centroid of a grid unit in a calculation range corresponding to the current simulated cannon in the full-model unstructured grid, and establishing the centroid unit includes: acquiring grid units in a calculation range corresponding to the current simulated cannon from the full-model unstructured grid according to the simulation requirements of the observation system; and determining and numbering each grid unit, and establishing a centroid unit.
Optionally, the establishing a corresponding relationship between the centroid unit and the hilbert space filling curve sequence includes: establishing a corresponding relation between the centroid unit and the Hilbert space filling curve sequence based on the serial number of the centroid unit contained in each Hilbert unit; wherein each Hilbert cell of the Hilbert space-filling curve comprises at most one centroid cell.
Optionally, the dividing, based on the hilbert space filling curve sequence from which the empty sequence number is deleted, the region formed by the grid cells within the calculation range includes: and according to the preset number of region division, dividing the Hilbert space filling curve sequence from which the empty sequence number is deleted, and dividing the region formed by the grid units in the corresponding calculation range.
Optionally, the full-model unstructured mesh is a triangular unit in a two-dimensional space and a tetrahedral unit in a three-dimensional space.
The invention has the beneficial effects that: the seismic wave numerical simulation method based on the space filling curve area division divides the grids into different calculation areas by using the area division based on the space filling curve, maintains the adjacent characteristics among the grids, ensures that the number of grid units of each sub-area is similar, improves the division efficiency of the non-structural grids, reduces the calculation consumption, and has high seismic wave numerical simulation calculation efficiency by carrying out numerical simulation calculation based on the divided areas.
The present invention has other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the invention.
Drawings
The above and other objects, features and advantages of the present invention will become more apparent by describing in more detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts throughout.
FIG. 1 shows a flow diagram of a seismic wave numerical simulation method based on space-filling curve region partitioning, according to an embodiment of the invention.
FIG. 2 is a diagram of a full-model unstructured grid of a seismic wave numerical simulation method based on space-filling curve region division according to an embodiment of the invention.
FIG. 3 is a schematic diagram of the corresponding calculation range of the current shot in the full-model unstructured grid based on the seismic wave numerical simulation method of space-filling curve region division according to an embodiment of the invention.
FIG. 4 is a diagram illustrating the generation of Hilbert space-filling curves of different orders and the resulting curves of a seismic wave numerical simulation method based on space-filling curve region partitioning according to an embodiment of the present invention.
Fig. 5 is a diagram illustrating a correspondence relationship between a centroid unit and a hilbert space-filling curve sequence in a seismic wave numerical simulation method based on space-filling curve region division according to an embodiment of the present invention.
FIG. 6 shows a two-dimensional and three-dimensional region-partitioned representation of a seismic wave numerical simulation method based on space-filling curve region partitioning according to an embodiment of the invention.
FIG. 7 shows a block diagram of a seismic wave numerical simulation apparatus based on space-filling curve region partitioning, according to an embodiment of the present invention.
Description of the reference numerals
102. A full model unstructured grid generation module; 104. a centroid unit establishing module; 106. a corresponding relation establishing module; 108. a module for deleting the empty serial number; 110. a dividing module; 112. and a calculation module.
Detailed Description
Preferred embodiments of the present invention will be described in more detail below. While the following describes preferred embodiments of the present invention, it should be understood that the present invention may be embodied in various forms and should not be limited by the embodiments set forth herein.
The invention provides a seismic wave numerical simulation method based on space filling curve area division, which comprises the following steps: generating a full model unstructured grid; acquiring the centroid of the grid unit in the corresponding calculation range of the current simulated cannon from the full-model unstructured grid, and establishing a centroid unit; filling a space formed by grid units in a calculation range by using a Hilbert space filling curve, and establishing a corresponding relation between a centroid unit and a Hilbert space filling curve sequence; deleting the empty serial numbers in the Hilbert space filling curve sequence based on the corresponding relation between the mass center unit and the Hilbert space filling curve sequence; dividing the region formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the empty serial numbers deleted; and performing seismic wave numerical simulation calculation based on the divided areas.
Specifically, a full-model unstructured grid is generated, the grid is a triangular unit in two dimensions and a tetrahedral unit in three dimensions, and the step is not specific to a grid generation method and only needs to generate a two-dimensional triangular or three-dimensional tetrahedral grid. Obtaining grid units and mass center coordinates in a corresponding calculation range of the current simulated cannon in a full-model non-structural grid to obtain mass center units, filling a space formed by the grid units in the calculation range by using a Hilbert space filling curve, and obtaining the corresponding relation between the mass center units and a Hilbert space filling curve sequence; since non-structural cells are not spatially regular, which may result in some cells of the Hilbert space-filling curve not containing centroid coordinates, these empty cells are culled when the neighboring centroids are closed. And cutting the Hilbert space filling curve sequence with the empty serial numbers deleted, so as to divide the region formed by the grid units in the corresponding calculation range, and performing seismic wave numerical simulation calculation based on the divided region.
In the shot-by-shot rolling observation mode of seismic exploration, the calculation range of each shot is different, region division needs to be carried out again according to the method, and seismic wave numerical simulation calculation is carried out based on the divided regions.
According to an exemplary embodiment, the seismic wave numerical simulation method based on the space filling curve region division divides grids into different calculation regions by using the region division based on the space filling curve, maintains the adjacent characteristics among the grids, ensures that the grid units of each sub-region are similar in number, improves the division efficiency of the non-structural grids, reduces the calculation consumption, performs numerical simulation calculation based on the divided regions, and has high seismic wave numerical simulation calculation efficiency.
As an alternative, obtaining the centroid of the grid unit in the corresponding calculation range of the current simulated cannon in the full-model unstructured grid, and establishing the centroid unit comprises the following steps: acquiring grid units in a calculation range corresponding to the current simulated cannon from the full-model unstructured grid according to the simulation requirements of the observation system; and determining the mass center of each grid unit, numbering the grid units and establishing a mass center unit.
Specifically, according to the simulation requirements of the observation system, the calculation range of the current simulated cannon and the grid units of the calculation range in the full-model unstructured grid are searched, and for the grid units in the current calculation range, the centroid coordinates of the grid units are used for replacing the grid units to obtain the centroid units. The purpose of this step is to reduce the computational consumption of seismic simulations, since a large portion of the area where the grid of the full model is simulated does not need to be simulated, resulting in a waste of resources. The main area is searched out based on the observation system, so that the calculation amount is reduced, and meanwhile, the simulation result is not influenced.
As an alternative, establishing the correspondence between the centroid unit and the hilbert space filling curve sequence includes: and establishing a corresponding relation between the centroid unit and the Hilbert space filling curve sequence based on the serial number of the centroid unit contained in each Hilbert unit.
Specifically, the number of a centroid unit included in each hilbert unit in the hilbert space filling curve is obtained, the hilbert units correspond to the included centroid unit numbers, and the corresponding relation between the centroid units and the hilbert space filling curve sequence is obtained.
Alternatively, each hilbert cell of the hilbert space-filling curve contains at most one centroid cell.
Specifically, the minimum order of the Hilbert space-filling curve is determined, so that each minimum unit of the Hilbert space-filling curve has only one centroid coordinate at most.
As an alternative, dividing the region formed by the grid cells within the calculation range based on the hilbert space filling curve sequence from which the empty sequence numbers are deleted includes: and according to the preset number of region division, dividing the Hilbert space filling curve sequence with the null sequence number deleted, so as to divide the region formed by the grid units in the corresponding calculation range.
Specifically, according to the set area division number, a one-dimensional Hilbert space filling curve is cut, namely a Hilbert space filling curve sequence with null sequence numbers deleted is cut, and meanwhile, areas formed by grid units represented by centroid units corresponding to the Hilbert space filling curve sequence are divided.
Alternatively, the full-model unstructured mesh is a triangular cell in two-dimensional space and a tetrahedral cell in three-dimensional space.
In a second aspect, the present invention also provides an electronic device, including: a memory storing executable instructions; and the processor runs executable instructions in the memory to realize the seismic wave numerical simulation method based on the space filling curve area division.
In a third aspect, the present invention further provides a computer-readable storage medium, where a computer program is stored, and when the computer program is executed by a processor, the method for seismic wave numerical simulation based on space-filling curve area division is implemented.
In a fourth aspect, the present invention further provides a seismic wave numerical simulation apparatus based on space filling curve region division, including: the full model unstructured grid generation module generates a full model unstructured grid; the centroid unit establishing module is used for acquiring the centroid of the grid unit in the corresponding calculation range of the current simulated cannon in the full-model non-structural grid and establishing a centroid unit; the corresponding relation establishing module is used for filling a space formed by grid units in the calculation range by using a Hilbert space filling curve and establishing the corresponding relation between a centroid unit and a Hilbert space filling curve sequence; the empty sequence number deleting module deletes the empty sequence number in the Hilbert space filling curve sequence based on the corresponding relation between the mass center unit and the Hilbert space filling curve sequence; the dividing module is used for dividing the region formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the null sequence numbers deleted; and the calculation module is used for performing seismic wave numerical simulation calculation based on the divided areas.
Specifically, a full-model unstructured grid is generated, the grid is a triangular unit in two dimensions and a tetrahedral unit in three dimensions, and the step is not specific to a grid generation method and only needs to generate a two-dimensional triangular or three-dimensional tetrahedral grid. Obtaining grid units and mass center coordinates in a corresponding calculation range of the current simulated cannon in a full-model non-structural grid to obtain mass center units, filling a space formed by the grid units in the calculation range by using a Hilbert space filling curve, and obtaining the corresponding relation between the mass center units and a Hilbert space filling curve sequence; since non-structural cells are not spatially regular, which may result in some cells of the Hilbert space-filling curve not containing centroid coordinates, these empty cells are culled when the neighboring centroids are closed. And cutting the Hilbert space filling curve sequence with the empty serial numbers deleted, so as to divide the region formed by the grid units in the corresponding calculation range, and performing seismic wave numerical simulation calculation based on the divided region.
In the shot-by-shot rolling observation mode of seismic exploration, the calculation range of each shot is different, region division needs to be carried out again according to the method, and seismic wave numerical simulation calculation is carried out based on the divided regions.
According to an exemplary embodiment, the seismic wave numerical simulation method based on the space filling curve region division divides grids into different calculation regions by using the region division based on the space filling curve, maintains the adjacent characteristics among the grids, ensures that the grid units of each sub-region are similar in number, improves the division efficiency of the non-structural grids, reduces the calculation consumption, performs numerical simulation calculation based on the divided regions, and has high seismic wave numerical simulation calculation efficiency.
As an alternative, obtaining the centroid of the grid unit in the corresponding calculation range of the current simulated cannon in the full-model unstructured grid, and establishing the centroid unit comprises the following steps: acquiring grid units in a calculation range corresponding to the current simulated cannon from the full-model unstructured grid according to the simulation requirements of the observation system; and determining the mass center of each grid unit, numbering the grid units and establishing a mass center unit.
Specifically, according to the simulation requirements of the observation system, the calculation range of the current simulated cannon and the grid units of the calculation range in the full-model unstructured grid are searched, and for the grid units in the current calculation range, the centroid coordinates of the grid units are used for replacing the grid units to obtain the centroid units. The purpose of this step is to reduce the computational consumption of seismic simulations, since a large portion of the area where the grid of the full model is simulated does not need to be simulated, resulting in a waste of resources. The main area is searched out based on the observation system, so that the calculation amount is reduced, and meanwhile, the simulation result is not influenced.
As an alternative, establishing the correspondence between the centroid unit and the hilbert space filling curve sequence includes: and establishing a corresponding relation between the centroid unit and the Hilbert space filling curve sequence based on the serial number of the centroid unit contained in each Hilbert unit.
Specifically, the number of a centroid unit included in each hilbert unit in the hilbert space filling curve is obtained, the hilbert units correspond to the included centroid unit numbers, and the corresponding relation between the centroid units and the hilbert space filling curve sequence is obtained.
Alternatively, each hilbert cell of the hilbert space-filling curve contains at most one centroid cell.
Specifically, the minimum order of the Hilbert space filling curve is determined, so that each minimum unit of the Hilbert space filling curve has only one centroid coordinate at most.
Alternatively, dividing the region formed by the grid cells in the calculation range based on the hilbert space filling curve sequence after deleting the empty sequence number comprises: and according to the preset number of region division, dividing the Hilbert space filling curve sequence with the null sequence number deleted, so as to divide the region formed by the grid units in the corresponding calculation range.
Specifically, according to the set area division number, a one-dimensional Hilbert space filling curve is cut, namely a Hilbert space filling curve sequence with null sequence numbers deleted is cut, and meanwhile, areas formed by grid units represented by centroid units corresponding to the Hilbert space filling curve sequence are divided.
Alternatively, the full-model unstructured mesh is a triangular cell in two-dimensional space and a tetrahedral cell in three-dimensional space.
Example one
FIG. 1 shows a flow diagram of a seismic wave numerical simulation method based on space-filling curve region partitioning, according to an embodiment of the invention. FIG. 2 is a diagram of a full-model unstructured grid of a seismic wave numerical simulation method based on space-filling curve region division according to an embodiment of the invention. FIG. 3 is a schematic diagram of the corresponding calculation range of the current shot in the full-model unstructured grid based on the seismic wave numerical simulation method of space-filling curve region division according to an embodiment of the invention. FIG. 4 is a diagram illustrating the generation of Hilbert space-filling curves of different orders and the resulting curves of a seismic wave numerical simulation method based on space-filling curve region partitioning according to an embodiment of the present invention. Fig. 5 is a diagram illustrating a correspondence relationship between a centroid unit and a hilbert space-filling curve sequence in a seismic wave numerical simulation method based on space-filling curve region division according to an embodiment of the present invention. FIG. 6 shows a two-dimensional and three-dimensional region-partitioned representation of a seismic wave numerical simulation method based on space-filling curve region partitioning, according to an embodiment of the invention.
Referring to fig. 1, fig. 2, fig. 3, fig. 4, fig. 5, and fig. 6, the method for numerically simulating seismic waves based on space-filling curve region division includes:
step 1: generating a full model unstructured grid;
the full-model unstructured grid is a triangular unit in a two-dimensional space and a tetrahedral unit in a three-dimensional space.
Step 2: acquiring the centroid of the grid unit in the corresponding calculation range of the current simulated cannon from the full-model unstructured grid, and establishing a centroid unit;
obtaining the centroid of the grid unit in the corresponding calculation range of the current simulated cannon in the full-model unstructured grid, wherein the establishing of the centroid unit comprises the following steps: acquiring grid units in a calculation range corresponding to the current simulated cannon from the full-model unstructured grid according to the simulation requirements of the observation system; and determining the mass center of each grid unit, numbering the grid units and establishing a mass center unit.
And step 3: filling a space formed by grid units in a calculation range by using a Hilbert space filling curve, and establishing a corresponding relation between a centroid unit and a Hilbert space filling curve sequence;
wherein, establishing the corresponding relation between the centroid unit and the Hilbert space filling curve sequence comprises: and establishing a corresponding relation between the centroid unit and the Hilbert space filling curve sequence based on the serial number of the centroid unit contained in each Hilbert unit.
Wherein, each Hilbert unit of the Hilbert space filling curve at most comprises a centroid unit.
And 4, step 4: deleting the empty serial numbers in the Hilbert space filling curve sequence based on the corresponding relation between the mass center unit and the Hilbert space filling curve sequence;
and 5: dividing the region formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the empty serial numbers deleted;
the method for dividing the region formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the null sequence numbers deleted comprises the following steps: and according to the preset number of region division, dividing the Hilbert space filling curve sequence with the null sequence number deleted, so as to divide the region formed by the grid units in the corresponding calculation range.
Step 6: and performing seismic wave numerical simulation calculation based on the divided areas.
Taking the full-model unstructured grid of a two-dimensional complex geological model as an example, the total grid number is 871921 cells. The model size is 42.62Km, Z10 Km. All grid cells are triangular cells, the cell arrangement has no specific direction, and the cell size has no certain rule. And searching the calculation range of the current simulated cannon and the corresponding grid unit in the full-model non-structure in the calculation range according to the simulation requirement of the observation system. As shown in FIG. 3, the model cutting range of the m-th shot and the k-th shot corresponding to the observation system with the observation range of 17km is a process of rolling along with the shots, the simulation of each shot needs to be performed again, the abscissa Position in FIG. 2 and FIG. 3 represents a point, and the Depth represents a Position. The number of units cut out from the simulation range of the k-cannon is 354924, and is only about 40% of the full model. Within the computational scope, the grid cells are replaced by the centroid coordinates of the cells. Fig. 4 shows Hilbert curves of different orders, and it can be seen that each order curve actually corresponds to an arrangement of small square cells, referred to herein as the cells of the Hilbert curve. And determining the minimum order of the Hilbert space filling curve, so that each minimum unit of the Hilbert space filling curve has only one centroid coordinate at most. In the step, the centroid of the triangular unit is taken as a slicing standard, and the order of the Hilbert curve is continuously increased, so that each Hilbert curve unit contains at most one centroid. And filling the area to be divided by using a Hilbert space filling curve, and establishing a corresponding relation between the centroid unit and the Hilbert curve sequence. As shown in fig. 5, after the regions to be divided are filled with Hilbert curves, the cell numbers of the triangular centroids contained in each Hilbert cell are shown, so that the corresponding relationship between the centroid cells and the Hilbert space filling curve sequence is established. And after the empty serial numbers on the Hilbert space filling curves are removed and the Hilbert units without mass center units in the graph 5 are removed, converting the two-dimensional units into one-dimensional sequences with space adjacent characteristics. And cutting the one-dimensional Hilbert space filling curve according to the set region division number, so as to divide the region formed by the centroid units in the corresponding calculation range. Fig. 6 shows two-dimensional and three-dimensional schematic diagrams of region partitioning. The seismic wave numerical simulation calculation is carried out according to the divided areas, in the figure 6, the abscissa Position represents a point, Depth represents a Position, X-Distance represents the Distance in the X direction, and Y-Distance represents the Distance in the Y direction.
Example two
FIG. 7 shows a block diagram of a seismic wave numerical simulation apparatus based on space-filling curve region partitioning, according to an embodiment of the present invention.
As shown in fig. 7, the seismic wave numerical simulation device based on space filling curve area division includes:
a full model unstructured grid generation module 102 for generating a full model unstructured grid;
a centroid unit establishing module 104, which obtains the centroid of the grid unit in the corresponding calculation range of the current simulated cannon in the full model non-structural grid, and establishes a centroid unit;
a corresponding relation establishing module 106, which uses a Hilbert space filling curve to fill the space formed by the grid units in the calculation range, and establishes a corresponding relation between the centroid unit and the Hilbert space filling curve sequence;
a null sequence number deleting module 108, configured to delete a null sequence number in the hilbert space-filling curve sequence based on a correspondence between the centroid unit and the hilbert space-filling curve sequence;
the dividing module 110 is configured to divide an area formed by the grid cells within the calculation range based on the hilbert space filling curve sequence from which the empty sequence number is deleted;
and the calculation module 112 is used for performing seismic wave numerical simulation calculation based on the divided areas.
Obtaining the centroid of the grid unit in the corresponding calculation range of the current simulated cannon in the full-model unstructured grid, wherein the establishing of the centroid unit comprises the following steps: acquiring grid units in a calculation range corresponding to the current simulated cannon from the full-model unstructured grid according to the simulation requirement of the observation system; and determining the mass center of each grid unit, numbering the grid units and establishing a mass center unit.
Wherein, establishing the corresponding relation between the centroid unit and the Hilbert space filling curve sequence comprises: and establishing a corresponding relation between the centroid unit and the Hilbert space filling curve sequence based on the serial number of the centroid unit contained in each Hilbert unit.
Wherein, each Hilbert unit of the Hilbert space filling curve at most comprises a centroid unit.
The method for dividing the area formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the null sequence numbers deleted comprises the following steps of: and according to the preset number of region division, dividing the Hilbert space filling curve sequence with the null sequence number deleted, so as to divide the region formed by the grid units in the corresponding calculation range.
The full-model unstructured grid is a triangular unit in a two-dimensional space and a tetrahedral unit in a three-dimensional space.
EXAMPLE III
The present disclosure provides an electronic device including: a memory storing executable instructions; and the processor runs executable instructions in the memory to realize the seismic wave numerical simulation based on the space filling curve area division.
An electronic device according to an embodiment of the present disclosure includes a memory and a processor.
The memory is to store non-transitory computer readable instructions. In particular, the memory may include one or more computer program products that 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), cache memory (cache), and/or the like. The non-volatile memory may include, for example, Read Only Memory (ROM), hard disk, flash memory, etc.
The processor may be a Central Processing Unit (CPU) or other form of processing unit having data processing capabilities and/or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the disclosure, the processor is configured to execute the computer readable instructions stored in the memory.
Those skilled in the art should understand that, in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus, an interface, and the like, and these well-known structures should also be included in the protection scope of the present disclosure.
For the detailed description of the present embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, which are not repeated herein.
Example four
The present disclosure provides a computer-readable storage medium storing a computer program which, when executed by a processor, implements the above-described seismic wave numerical simulation based on space-filling curve region division.
A computer-readable storage medium according to an embodiment of the present disclosure has non-transitory computer-readable instructions stored thereon. The non-transitory computer readable instructions, when executed by a processor, perform all or a portion of the steps of the methods of the embodiments of the disclosure previously described.
The computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROMs and DVDs), magneto-optical storage media (e.g., MOs), magnetic storage media (e.g., magnetic tapes or removable disks), media with built-in rewritable non-volatile memory (e.g., memory cards), and media with built-in ROMs (e.g., ROM cartridges).
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. 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 (9)
1. A seismic wave numerical simulation method based on space filling curve area division is characterized by comprising the following steps:
generating a full model unstructured grid;
acquiring the centroid of the grid unit in the corresponding calculation range of the current simulated cannon from the full-model unstructured grid, and establishing a centroid unit;
filling a space formed by grid units in the calculation range by using a Hilbert space filling curve, and establishing a corresponding relation between a centroid unit and a Hilbert space filling curve sequence;
based on the corresponding relation between the centroid unit and the Hilbert space filling curve sequence, deleting the empty serial number in the Hilbert space filling curve sequence;
dividing the region formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the empty serial numbers deleted;
and performing seismic wave numerical simulation calculation based on the divided areas.
2. The seismic wave numerical simulation method based on the space-filling curve area division according to claim 1, wherein the obtaining of the centroid of the grid cell within the calculation range corresponding to the current simulated shot in the full-model unstructured grid, and the establishing of the centroid cell comprises:
acquiring grid units in a calculation range corresponding to the current simulated cannon from the full-model unstructured grid according to the simulation requirements of the observation system;
and determining the mass center of each grid unit, numbering the grid units and establishing a mass center unit.
3. The method for seismic wave numerical simulation based on space-filling curve area division according to claim 1, wherein the establishing of the correspondence between the centroid unit and the hilbert space-filling curve sequence comprises:
establishing a corresponding relation between the centroid unit and the Hilbert space filling curve sequence based on the serial number of the centroid unit contained in each Hilbert unit;
wherein each Hilbert cell of the Hilbert space-filling curve comprises at most one centroid cell.
4. The seismic wave numerical simulation method based on the space-filling curve region division according to claim 3, wherein the deleting the empty sequence number in the hilbert space-filling curve sequence based on the correspondence between the centroid unit and the hilbert space-filling curve sequence comprises:
in the hilbert space filling curve sequence, the sequence numbers corresponding to the hilbert cells not including the centroid cell are deleted.
5. The seismic wave numerical simulation method based on the space-filling curve region division according to claim 1, wherein the dividing the region formed by the grid cells within the calculation range based on the hilbert space-filling curve sequence from which the empty sequence numbers are deleted comprises:
and according to the preset number of region division, dividing the Hilbert space filling curve sequence with the null sequence number deleted, so as to divide the region formed by the grid units in the corresponding calculation range.
6. The method of numerical seismic wave modeling based on the partitioning of the space-filling curve region of claim 1, wherein said full-model unstructured mesh is a triangular cell in two dimensions and a tetrahedral cell in three dimensions.
7. An electronic device, characterized in that the electronic device comprises:
a memory storing executable instructions;
a processor executing the executable instructions in the memory to implement the seismic wave numerical simulation method based on space-filling curve region partitioning according to any one of claims 1 to 6.
8. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program which, when executed by a processor, implements the seismic wave numerical simulation method based on space-filling curve region partitioning of any one of claims 1 to 6.
9. A seismic wave numerical simulation device based on space filling curve area division is characterized by comprising the following components:
the full model unstructured grid generation module generates a full model unstructured grid;
the centroid unit establishing module is used for acquiring the centroid of the grid unit in the corresponding calculation range of the current simulated cannon in the full-model non-structural grid and establishing a centroid unit;
the corresponding relation establishing module is used for filling the space formed by the grid units in the calculation range by using a Hilbert space filling curve and establishing the corresponding relation between the centroid unit and the Hilbert space filling curve sequence;
a null sequence number deleting module, which deletes the null sequence number in the Hilbert space filling curve sequence based on the corresponding relation between the centroid unit and the Hilbert space filling curve sequence;
the dividing module is used for dividing the region formed by the grid units in the calculation range based on the Hilbert space filling curve sequence with the empty serial numbers deleted;
and the calculation module is used for performing seismic wave numerical simulation calculation based on the divided areas.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011150064.1A CN114491897A (en) | 2020-10-23 | 2020-10-23 | Seismic wave numerical simulation method, device, medium and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011150064.1A CN114491897A (en) | 2020-10-23 | 2020-10-23 | Seismic wave numerical simulation method, device, medium and electronic equipment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114491897A true CN114491897A (en) | 2022-05-13 |
Family
ID=81470213
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011150064.1A Pending CN114491897A (en) | 2020-10-23 | 2020-10-23 | Seismic wave numerical simulation method, device, medium and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114491897A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116089076A (en) * | 2022-12-30 | 2023-05-09 | 清华大学 | Data processing method, device, computer equipment and storage medium |
-
2020
- 2020-10-23 CN CN202011150064.1A patent/CN114491897A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116089076A (en) * | 2022-12-30 | 2023-05-09 | 清华大学 | Data processing method, device, computer equipment and storage medium |
| CN116089076B (en) * | 2022-12-30 | 2023-11-28 | 清华大学 | Data processing method, device, computer equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| D’Ambrosio et al. | Simulating debris flows through a hexagonal cellular automata model: Sciddica s 3–hex | |
| CA2846327C (en) | Systems and methods for generating a large scale polygonal mesh | |
| EP2869096A2 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| CN106934826B (en) | Rock slope structure refined modeling and block identification method | |
| CN111768502A (en) | Non-structural grid two-dimensional flood simulation system based on GPU acceleration technology | |
| CN110992458A (en) | Large-scale oblique photography model organization and scheduling method | |
| Suchde et al. | Point cloud generation for meshfree methods: An overview | |
| CN116030218A (en) | Tetrahedral mesh dividing method, device, system and storage medium | |
| CN110221342A (en) | Seismic source location method, apparatus and storage medium based on three-dimensional velocity structure | |
| CN111695281A (en) | Particle rapid positioning method for tetrahedral mesh division finite element particle simulation | |
| CN114491897A (en) | Seismic wave numerical simulation method, device, medium and electronic equipment | |
| CN107886573B (en) | Slope three-dimensional finite element grid generation method under complex geological conditions | |
| US11301988B2 (en) | Reverse engineering data analysis system, and integrated circuit component data processing tool and method thereof | |
| CN109490948A (en) | Seismoacoustics wave equation vector parallel calculating method | |
| CN110968930B (en) | Geological variable attribute interpolation method and system | |
| CN109584369B (en) | Actual stratum full hexahedron grid generation method and device | |
| CN111859268B (en) | Magnetic tensor abnormal spatial domain fast forward modeling method based on grid point lattice | |
| Weier et al. | Generating and rendering large scale tiled plant populations | |
| Khamayseh et al. | Use of the spatial kD-tree in computational physics applications | |
| US8271518B2 (en) | Nearest neighbor search method | |
| CN110349265A (en) | A kind of tetrahedron topological net generation method and electronic equipment | |
| JP2008197921A (en) | Simulation device, simulation method and simulation program | |
| RU2611892C1 (en) | Method of three-dimensional simulation of specified hydrogeological feature implemented in computer system | |
| RU2815465C1 (en) | System for semi-deterministic modeling of fracturing based on random matrix | |
| EP3451191A1 (en) | Computer implemented method for manipulating a numerical model of a 3d domain |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |