CN114545505A - Dynamic ocean wave field information acquisition method based on particle grid thought - Google Patents
Dynamic ocean wave field information acquisition method based on particle grid thought Download PDFInfo
- Publication number
- CN114545505A CN114545505A CN202210032055.5A CN202210032055A CN114545505A CN 114545505 A CN114545505 A CN 114545505A CN 202210032055 A CN202210032055 A CN 202210032055A CN 114545505 A CN114545505 A CN 114545505A
- Authority
- CN
- China
- Prior art keywords
- grid
- wave field
- detector
- wave
- field information
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 239000002245 particle Substances 0.000 title claims abstract description 23
- 238000001228 spectrum Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000004088 simulation Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, e.g. survey design
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Oceanography (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a dynamic ocean wave field information collection method based on a particle grid thought, and relates to the technical field of seismic wave field simulation. The method comprises the following steps: acquiring a wave field grid at each moment; generating a detector position according to the combination of the Gerstner wave and the sea wave spectrum; searching a corresponding wave field grid according to the position of the detector, and calculating wave field information at the position of the detector according to the wave field information of the wave field grid; collecting wave field values near the detector according to the PIC idea; and integrating the wave field information values to obtain the seismic record. The invention can simulate the fluctuation state of the geophone along with the sea surface, and the wave field information of the corresponding position is obtained in an interpolation mode in the state.
Description
Technical Field
The invention relates to the technical field of seismic wave field simulation, in particular to a dynamic ocean wave field information acquisition method based on a particle grid thought.
Background
The Gerstner model was introduced for the first time in 1986 by Fournier et al in the field of computer image processing, and describes the motion state of individual particles on the sea surface mainly from a kinetic perspective, and is generally expressed using parametric equations:
where (x, y) represents the current position of the calculated particle, (x)0，y0) Representing the initial position of the calculated particle, r the radius of motion (which can also be understood as the amount of energy) and w the frequency of motion (angular velocity). In analytic geometry, the parametric equation describes the motion trajectory of a circle.
In conjunction with the wave spectrum, take the PM spectrum as an example:
where α and β are constants, and are generally taken as a ═ 0.0081, β ═ 0.74,
representing the gravitational acceleration, w is the frequency of the waves, U represents the wind speed 19.5 meters above the sea surface.
The position of each particle on the sea surface, i.e. the position of the geophone described in the present invention, can be determined by substituting r into the above equation.
The Particle grid (PIC) method is a method proposed in the last 70 th century to solve the problems of fluid simulation, etc., and the PIC has been applied in many fields after decades of development. The core of the PIC method is the algorithm of particle-to-grid (P2G) and grid-to-particle (G2P).
As shown in fig. 6, the PIC method divides the object solution into two parts, one part is a particle part and the other part is a grid part. The specific operation method is that the initial conditions on the particle part are distributed to the grid nodes in a weighting mode, calculation is carried out on the grid nodes, then the calculation result is returned to the particle, the particle is updated, and the required result can be obtained by circulating the steps.
In typical seismic numerical modeling calculations, geophones are often placed on grid points to receive wavefield information (e.g., pressure, velocity components, etc.) at the grid points. However, the method has great limitation, and in a simulated marine environment, when the geophone vibrates in a small range in a grid, the geophone is often fixed on the grid point during numerical calculation and cannot fluctuate along with the sea surface, so that more specific wave field information cannot be acquired.
Disclosure of Invention
In view of the above, the invention provides a dynamic marine wave field information acquisition method based on a particle grid thought, which can simulate the fluctuation state of a detector along with the sea surface, and obtain wave field information of corresponding positions in an interpolation mode in the state.
In order to achieve the purpose, the invention adopts the following technical scheme:
a dynamic ocean wave field information collection method based on a particle grid idea comprises the following steps:
acquiring a wave field grid at each moment;
generating a detector position according to the combination of the Gerstner wave and the sea wave spectrum;
searching a corresponding wave field grid according to the position of the wave detector, and calculating wave field information of the grid where the wave detector is located according to the wave field information of the wave field grid;
collecting wave field values near the detector in the grid according to the PIC thought and the calculated value;
and (4) integrating wave field information values, namely recording the detector value corresponding to each moment as a row of a result matrix, and recording all moments in the result matrix to finally obtain the seismic record.
Optionally, the wavefield information includes pressure, velocity component, grid spacing, and grid size for each differential grid point.
Optionally, the wavefield grid acquiring process includes: obtaining wave field grid of each time by solving wave equation by finite difference method
Optionally, the wavefield grid is a weight grid.
Optionally, when the PIC idea is used, the kernel function is selected according to actual conditions.
According to the technical scheme, compared with the prior art, the dynamic ocean wave field information acquisition method based on the particle grid idea is disclosed, the initial detector distribution generated by using the Gerstner method accords with common knowledge, then the PIC method is used for interpolation calculation, the influence of the action range is considered, the wave field information obtained after interpolation has smaller error, the seismic records under the condition of the dynamic detectors can be displayed, and data support basically consistent with the reality is provided for research in aspects of geophysical field simulation and the like.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
FIG. 1 is a flow chart of an implementation of the present invention;
FIG. 2 is a seismic record (homogeneous medium) of the embodiment at a distance of 19.5m from the sea surface at 6 m/s;
FIG. 3 is a seismic record (homogeneous medium) of the embodiment at a distance of 19.5m from the sea surface at 8 m/s;
FIG. 4 is a seismic record (homogeneous medium) of the embodiment at a distance of 19.5m from the sea surface at 12 m/s;
FIG. 5 is a seismic record (double layer medium) for an example at 8m/s at 19.5m from the sea surface;
FIG. 6 is a diagram of an implementation concept of a particle grid method according to an embodiment;
FIG. 7 is a schematic diagram of grid points and detector locations according to an embodiment;
FIG. 8 is a schematic diagram of a two-layer model of an embodiment;
FIG. 9 is a schematic diagram of a homogeneous model according to an embodiment.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the invention mainly utilizes the idea of grid-to-particle conversion. The detector is first assumed to be a particle that can move at will and the wavefield recording is assumed to be a grid. As shown in fig. 7, where the black nodes are grid nodes and the gray nodes are detector locations. The wavefield information on the grid is then distributed onto the detectors using a kernel function (weight function). The form of the kernel function (taking cubic spline kernel as an example) is:
wherein w (r, h) is weight, r is distance between grid node and detector, and h is kernel function radius.
According to the weight, the wave field value on the detector can be obtained:
rec＝∑wp (5)
where rec is the value collected by the detector, w is the weight, and p is the wave field value at the grid point.
Specifically, the present embodiment provides a dynamic marine wave field information collection method based on a particle grid concept, which can show seismic records under a dynamic geophone condition, and provide data support basically consistent with the actual data for research in geophysical field simulation and other aspects, please refer to fig. 1 to 8, and includes:
the method comprises the following steps: and solving the wave equation by using a finite difference method to obtain wave field information at each moment. Extracting wave field information of each moment, including pressure, velocity component and the like of each differential grid point, and simultaneously recording information such as grid distance, grid size and the like;
step two: setting the position of the detector according to requirements, and assuming that the detector moves along with the sea surface continuously, the position of the detector is coincided with the sea surface at any time, so that the position of the detector is set by combining Gerstner waves and a sea wave spectrum;
step three: calculating each detector, firstly judging the grid position of the detector, searching a weight grid according to the grid position, and calculating wave field information at the position of the detector according to the wave field information on the grid;
step four: and integrating the wave field values to finally obtain corresponding seismic records.
To better illustrate the effects of the above embodiments, a specific example is given below:
example (c): assuming that the wind speed at 19.5m from the sea surface in a certain sea area is 6m/s, 8m/s or 12m/s, the calculation is performed according to the method, the first model is a uniform medium model (fig. 9), and the finite difference grid size is firstly set as follows: setting the seismic source position as follows: the boundary condition is a PML absorption boundary. The detectors are set to 1000 evenly distributed at 400-direction grid node positions. The simulation flow is shown in figure 1. The final result is: the seismic record is shown in figure 3 under the condition of 6m/s wind speed; the seismic record is shown in figure 4 under the condition of 8m/s wind speed; the seismic record is shown in figure 4 at a wind speed of 12 m/s. The second model is a two-layer model (fig. 8), and as the above conditions, there is a speed change interface at the 200-direction node, and assuming that the wind speed is 8m/s, the result can be shown in fig. 5.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (5)
1. A dynamic ocean wave field information collection method based on a particle grid idea is characterized by comprising the following steps:
acquiring a wave field grid at each moment;
generating a detector position according to the combination of the Gerstner wave and the sea wave spectrum;
searching a corresponding wave field grid according to the position of the detector, and calculating wave field information of the position of the detector in a 3 x 3 grid according to the wave field information of the wave field grid;
collecting wave field values near the detector in a 3 multiplied by 3 grid according to PIC thought;
and (4) integrating wave field information values, namely recording the detector value corresponding to each moment as a row of a result matrix, and recording all moments in the result matrix to finally obtain the seismic record.
2. The method of claim 1, wherein the wavefield information comprises pressure, velocity component, grid spacing, and grid size for each differential grid point.
3. The method for collecting dynamic ocean wave field information based on particle grid thought as claimed in claim 1, wherein the wave field grid obtaining process is: and solving the wave equation by using a finite difference method to obtain the wave field grid at each moment.
4. The method of claim 1, wherein the wavefield grid is a 3 x 3 weight grid.
5. The method of claim 1, further comprising selecting kernel functions based on the fact when using PIC concept.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210032055.5A CN114545505A (en) | 2022-01-12 | 2022-01-12 | Dynamic ocean wave field information acquisition method based on particle grid thought |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210032055.5A CN114545505A (en) | 2022-01-12 | 2022-01-12 | Dynamic ocean wave field information acquisition method based on particle grid thought |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114545505A true CN114545505A (en) | 2022-05-27 |
Family
ID=81670776
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210032055.5A Pending CN114545505A (en) | 2022-01-12 | 2022-01-12 | Dynamic ocean wave field information acquisition method based on particle grid thought |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114545505A (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5999488A (en) * | 1998-04-27 | 1999-12-07 | Phillips Petroleum Company | Method and apparatus for migration by finite differences |
| US6125330A (en) * | 1997-09-05 | 2000-09-26 | Schlumberger Technology Corporation | Method of determining the response caused by model alterations in seismic simulations |
| WO2007147150A2 (en) * | 2006-06-15 | 2007-12-21 | Lightspeed Logic, Inc. | Simultaneous dynamical integration applied to detailed placement |
| WO2007149717A2 (en) * | 2006-06-08 | 2007-12-27 | Lightspeed Logic, Inc. | Morphing for global placement using integer linear programming |
| CN108181652A (en) * | 2017-12-29 | 2018-06-19 | 中国海洋大学 | A kind of subsea node seismic data uplink and downlink wave field numerical method |
| KR102006325B1 (en) * | 2018-12-21 | 2019-08-01 | 한국해양과학기술원 | Method for simulation of sea state changes using base spectrum |
| WO2020160749A1 (en) * | 2019-02-04 | 2020-08-13 | Siemens Industry Software Nv | Optimising ship noise radiation using digital twins and controls |
| CN111797552A (en) * | 2020-06-12 | 2020-10-20 | 中国海洋大学 | Numerical data simulation method for undulating sea surface seismic wave field based on sea wave spectrum |
| CN111948708A (en) * | 2019-10-18 | 2020-11-17 | 中国石油大学(北京) | Seismic wave field forward modeling method for dipping in undulating surface of boundary |
| CN113536638A (en) * | 2021-07-22 | 2021-10-22 | 北京大学 | High-precision seismic wave field simulation method based on discontinuous finite elements and staggered grids |
-
2022
- 2022-01-12 CN CN202210032055.5A patent/CN114545505A/en active Pending
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6125330A (en) * | 1997-09-05 | 2000-09-26 | Schlumberger Technology Corporation | Method of determining the response caused by model alterations in seismic simulations |
| US5999488A (en) * | 1998-04-27 | 1999-12-07 | Phillips Petroleum Company | Method and apparatus for migration by finite differences |
| WO2007149717A2 (en) * | 2006-06-08 | 2007-12-27 | Lightspeed Logic, Inc. | Morphing for global placement using integer linear programming |
| WO2007147150A2 (en) * | 2006-06-15 | 2007-12-21 | Lightspeed Logic, Inc. | Simultaneous dynamical integration applied to detailed placement |
| CN108181652A (en) * | 2017-12-29 | 2018-06-19 | 中国海洋大学 | A kind of subsea node seismic data uplink and downlink wave field numerical method |
| KR102006325B1 (en) * | 2018-12-21 | 2019-08-01 | 한국해양과학기술원 | Method for simulation of sea state changes using base spectrum |
| WO2020160749A1 (en) * | 2019-02-04 | 2020-08-13 | Siemens Industry Software Nv | Optimising ship noise radiation using digital twins and controls |
| CN111948708A (en) * | 2019-10-18 | 2020-11-17 | 中国石油大学(北京) | Seismic wave field forward modeling method for dipping in undulating surface of boundary |
| CN111797552A (en) * | 2020-06-12 | 2020-10-20 | 中国海洋大学 | Numerical data simulation method for undulating sea surface seismic wave field based on sea wave spectrum |
| CN113536638A (en) * | 2021-07-22 | 2021-10-22 | 北京大学 | High-precision seismic wave field simulation method based on discontinuous finite elements and staggered grids |
Non-Patent Citations (2)
| Title |
|---|
| 刘煜: "质点—网格海冰模式在渤海的数值预报试验", 《海洋预报》, vol. 22, 31 May 2005 (2005-05-31), pages 36 * |
| 杨志钢: "基于Gerstner函数和JONSWAP海浪谱的海浪仿真", 《哈尔滨商业大学学报（自然科学版）》, vol. 37, no. 6, 30 November 2021 (2021-11-30), pages 677 - 678 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111239802B (en) | Deep learning speed modeling method based on seismic reflection waveform and velocity spectrum | |
| CN106772583B (en) | A kind of earthquake diffracted wave separation method and device | |
| CN107908913B (en) | Earth power digital-analog method based on parallel computer | |
| EP3736776B1 (en) | Apparatus, system and method for the generation of polygonal meshes | |
| CN111833450B (en) | Ultrasonic three-dimensional rapid reconstruction and analysis method fused with finite element analysis method | |
| CN115906256A (en) | Reservoir landslide surge numerical simulation method and system | |
| CN106959467A (en) | Seimic wave velocity inversion method and device | |
| CN114839673A (en) | Separation method, separation system and computer equipment for multi-seismic-source efficient acquisition wave field | |
| CN114330163B (en) | Typhoon-ionosphere disturbance dynamics model modeling method for high-frequency ground wave beyond-visual-range radar | |
| CN115540828A (en) | Internal wave forecasting method based on wall sensor | |
| Tarbell et al. | An example of the initialization of the divergent wind component in a mesoscale numerical weather prediction model | |
| CN117113725B (en) | Energy obtaining analysis method of parameter pendulum type wave energy device | |
| CN112101440B (en) | Coronage substance ejection detection method based on improved U-Net network | |
| CN114545505A (en) | Dynamic ocean wave field information acquisition method based on particle grid thought | |
| CN108983290B (en) | It is a kind of three-dimensional Method in Transverse Isotropic Medium in travel when determine method and system | |
| CN116503716A (en) | Radar image derivatization and database capacity expansion method | |
| CN115619643A (en) | Construction method for super-resolution network model of rotator image | |
| Kirwan Jr et al. | Ring evolution in general circulation models from path analysis | |
| CN111257833A (en) | Sound source identification method based on Laplace norm for fast iterative shrinkage threshold | |
| CN110988991A (en) | Elastic parameter inversion method, device and system | |
| CN110231152A (en) | Derrick detection method and apparatus based on unmanned plane | |
| CN116660982B (en) | Full waveform inversion method based on attention convolution neural network | |
| CN115685363A (en) | Multi-parameter gravity field data fusion method | |
| Sokołowski et al. | Wind Turbines Modeling as the Tool for Developing Algorithms of Processing their Video Recordings | |
| CN116993752B (en) | Semantic segmentation method, medium and system for live-action three-dimensional Mesh model |
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 |