CN112505774A - 一种地震声波数值模拟中的组合边界方法 - Google Patents
一种地震声波数值模拟中的组合边界方法 Download PDFInfo
- Publication number
- CN112505774A CN112505774A CN202011477440.8A CN202011477440A CN112505774A CN 112505774 A CN112505774 A CN 112505774A CN 202011477440 A CN202011477440 A CN 202011477440A CN 112505774 A CN112505774 A CN 112505774A
- Authority
- CN
- China
- Prior art keywords
- boundary
- pml
- attenuation
- setting
- wave
- 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.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000010521 absorption reaction Methods 0.000 claims abstract description 28
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000000694 effects Effects 0.000 abstract description 24
- 230000002238 attenuated effect Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005284 excitation Effects 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种地震声波数值模拟中的组合边界方法,包括给定数值模拟速度模型,设置参数,读取速度模型,确定震源位置同时加入震源项;在有限差分数值模拟中加入PML边界条件,设置PML衰减因子在计算区域为0,衰减区域为一个数值,使用基于PML边界的波动方程进行波场模拟;在衰减区域的外边界处设置二阶CE边界条件,吸收经过衰减带后而未被完全吸收的入射波,经过二次吸收后,产生的边界反射更弱;输出得到不同时刻的波场快照。该方法对二阶CE边界条件进行改进,与PML边界条件组合,相对于单纯使用PML边界,这种组合方法对外边界处的入射波能量强弱要求不高,可在减少衰减带厚度的同时提高计算效率的效果。
Description
技术领域
本发明属于地震声波数值模拟技术领域,具体涉及一种地震声波数值模拟过程中使用的一种组合边界方法,该方法将CE边界条件和PML边界条件组合为一种新的组合边界条件。
背景技术
在地震声波数值模拟计算过程中选取的有限计算区域会产生边界反射,干扰正常波场计算,因此需要引入人工边界条件降低边界反射的影响。在数值模拟中,人工边界条件按方法原理可以分为吸收类和衰减类边界条件两大类。其中,吸收类边界条件通过在边界处使用单程波方程模拟地震波的传播,使入射波只向计算区域外传播而不会产生边界反射。衰减类边界条件的方法原理是在计算区域外引入一层衰减区域,入射波在衰减层内传播时,入射波能量逐渐衰减从而不会产生边界反射。
基于傍轴近似原理的Clayton-Engquist吸收边界条件(以下简称CE边界条件)是吸收类边界条件中最经典的一种边界条件。CE边界条件在边界处使用不同阶数(通常为一阶或二阶)的单程波方程吸收边界反射波的能量,在计算区域使用正常差分方程计算波场,利用单程波方程吸收边界反射,占用内存小,方法简便且在程序上易于实现。但是,CE边界条件有两点不足之处:1、低阶CE边界条件精度低,吸收效果不好,高阶CE边界条件不稳定且计算难度增大;2、吸收类边界条件受入射角的限制较大研究表明,CE边界条件对入射角小于45度的入射波具有较好的吸收效果。
在衰减类边界条件的方法中,完全匹配层法(PML)应用效果最优,近年来得到广泛的应用。与其他边界条件相比,完全匹配层法具有最佳的吸收效果,但由于在计算区域外构建了一层衰减区域,若要实现较好的吸收效果对于衰减层的厚度也有一定的要求,衰减带厚度的增加意味着计算量的增大以及计算效率的降低。
发明内容
本发明的目的在于提供一种地震声波数值模拟中的组合边界方法,以克服在地震波场数值模拟中,衰减层厚度会影响完全匹配层法吸收效果的问题。本发明对二阶CE边界条件进行改进,与PML边界条件组合,先使用PML边界条件衰减入射波,再使用二阶CE边界吸收未衰减完的入射波,得到一种具有较好吸收效果的新组合边界。
本发明的目的是通过以下技术方案实现的:
一种地震声波数值模拟中的组合边界方法,包括以下步骤:
A、给定数值模拟速度模型,设置x、z方向网格数、有限差分时间步长以及空间步长等参数,读取速度模型,确定震源位置同时加入震源项;
B、在有限差分数值模拟中加入PML边界条件,在计算区域四周构建一个衰减带,设置PML衰减因子在计算区域为0,衰减区域为一个数值,使用基于PML边界的波动方程进行波场模拟;
C、在衰减区域的外边界处设置二阶CE边界条件,吸收经过衰减带后而未被完全吸收的入射波,经过二次吸收后,产生的边界反射更弱;
D、经过上述步骤运算,可以输出得到不同时刻的波场快照。
与现有技术相比,本发明的有益效果在于:
本发明地震声波数值模拟中的组合边界方法,结合完全匹配层法实现效果好以及二阶CE边界条件易实现,占用内存小的优点,先使用PML边界条件衰减入射波,再在衰减区域外边界处使用二阶CE边界条件吸收未衰减完全的入射波;相对于单纯使用PML边界,这种组合方法对到达外边界处的入射波能量强弱要求不高,因此可以使用较窄的衰减带,从而达到在减少衰减带厚度的同时提高计算效率的效果。
附图说明
图1为本发明地震声波数值模拟中的组合边界方法的模拟流程;
图2为均匀介质模型,震源位置(900m,900m);
图3a-图3c为均匀介质模型下仅使用PML边界条件与使用组合边界条件吸收效果对比;
图4为三层均匀介质模型,震源位置(900m,240m);
图5a-图5c为三层均匀介质模型下仅使用PML边界条件与使用组合边界条件吸收效果对比。
具体实施方式
下面结合实施例对本发明作进一步说明:
研究发现,CE边界条件对入射角小于45度的入射波具有较好的吸收效果。因此,本发明对二阶CE边界条件进行改进,与PML边界条件组合,先使用PML边界条件衰减入射波,再使用二阶CE边界吸收未衰减完的入射波,得到一种具有较好吸收效果的新组合边界。
A、进行地震声波波场有限差分数值模拟,给定数值模拟速度模型,输入x、z方向网格数,时间步长以及空间步长等参数,读取速度模型,确定震源位置同时加入震源项;
B、在波场数值模拟过程中加入PML边界条件,首先在计算区域四周构建一个衰减带,设置PML衰减因子,使用基于PML边界的波动方程进行波场模拟,基于PML边界的波动方程即相对于常规的波动方程加入了一个衰减项;在计算区域衰减因子设置为0,衰减区域衰减因子设置为一个数值,这样地震波在计算区域内传播不会有影响,而当地震波传播到衰减带内时会逐步衰减。
C、在衰减区域的外边界处即衰减带的最外层设置二阶CE边界条件,即用单向傍轴近似波动方程模拟波能量由内向外的单向移动,吸收经过衰减带后而未被完全吸收的入射波,经过二次吸收后,产生的边界反射更弱。因为,对入射波有二次吸收,所以对到达衰减带边界处的波能量强弱的要求不大,可以在一定程度上减少衰减带厚度。
D、经过上述步骤运算,可以输出得到不同时刻的波场快照,具体描述声波到达边界处的边界反射吸收效果。
实施例
在相同的衰减带厚度下,分析仅使用PML边界条件与使用组合边界条件对边界反射的吸收效果,以及增加衰减带厚度后,PML边界条件的吸收效果。设置了三种速度模型,模型网格大小设置为300*300,网格间距为6m*6m,时间步长为0.8ms,选用主频30赫兹的雷克子波作为激发震源,采用时间2阶空间12阶的高阶有限差分进行声波波场数值模拟。
均匀介质模型下,取衰减带厚度为120m,采用组合边界条件和PML边界条件下的吸收效果如图3a,3b,采用衰减带厚度为210m的PML边界条件的吸收效果如图3c所示。
三层均匀介质模型下,取衰减带厚度为120m,采用组合边界条件和PML边界条件下的吸收效果如图5a,5b,采用衰减带厚度为210m的PML边界条件的吸收效果如图5c所示。
通过图3a与图3b、图5a与图5b的对比可以看出,在采用相同衰减带厚度时,组合边界条件具有较好的吸收效果。通过图3a与图3c、图5a与图5c的对比可以看出,在采用相对较少的衰减厚度时也能达到很好的吸收效果,说明组合边界的应用在保证吸收效果的同时可以有效的提高计算效率。验证了本发明提出的组合边界的有效性。
Claims (1)
1.一种地震声波数值模拟中的组合边界方法,包括以下步骤:
A、给定数值模拟速度模型,设置x、z方向网格数、有限差分时间步长以及空间步长参数,读取速度模型,确定震源位置同时加入震源项;
B、在有限差分数值模拟中加入PML边界条件,在计算区域四周构建一个衰减带,设置PML衰减因子在计算区域为0,衰减区域为一个数值,使用基于PML边界的波动方程进行波场模拟;
C、在衰减区域的外边界处设置二阶CE边界条件,吸收经过衰减带后而未被完全吸收的入射波,经过二次吸收后,产生的边界反射更弱;
D、经过上述步骤运算,可以输出得到不同时刻的波场快照。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011477440.8A CN112505774B (zh) | 2020-12-15 | 2020-12-15 | 一种地震声波数值模拟中的组合边界方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011477440.8A CN112505774B (zh) | 2020-12-15 | 2020-12-15 | 一种地震声波数值模拟中的组合边界方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112505774A true CN112505774A (zh) | 2021-03-16 |
CN112505774B CN112505774B (zh) | 2022-12-27 |
Family
ID=74973577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011477440.8A Active CN112505774B (zh) | 2020-12-15 | 2020-12-15 | 一种地震声波数值模拟中的组合边界方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112505774B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6687659B1 (en) * | 2000-03-24 | 2004-02-03 | Conocophillips Company | Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications |
CN103698814A (zh) * | 2013-12-31 | 2014-04-02 | 中国海洋石油总公司 | 一种用于变密度声波方程的混合吸收边界条件的实现方法 |
CN103823239A (zh) * | 2013-10-13 | 2014-05-28 | 中国石油集团西北地质研究所 | 频率域优化混合交错网格有限差分正演模拟方法 |
CN104237944A (zh) * | 2014-10-09 | 2014-12-24 | 王兵 | 一种适用于交错网格有限差分的全吸收pml方法 |
CN105447225A (zh) * | 2015-11-06 | 2016-03-30 | 中国海洋大学 | 一种应用于声波有限差分数值模拟的组合吸收边界条件 |
US20160238723A1 (en) * | 2015-02-13 | 2016-08-18 | Valeriy Brytik | Efficient and Stable Absorbing Boundary Condition In Finite-Difference Calculations |
CN108073732A (zh) * | 2016-11-10 | 2018-05-25 | 中国石油化工股份有限公司 | 获得稳定的近最佳匹配层吸收边界条件的方法 |
CN109188517A (zh) * | 2018-09-03 | 2019-01-11 | 中国海洋大学 | 基于Higdon余弦型加权的混合吸收边界条件方法 |
US20190115002A1 (en) * | 2017-10-16 | 2019-04-18 | The Hong Kong University Of Science And Technology | Sound absorber with stair-stepping structure |
CN109725351A (zh) * | 2018-12-18 | 2019-05-07 | 中国石油天然气集团有限公司 | 一种3d弹性波混合吸收边界条件的确定方法、装置及系统 |
CN111538083A (zh) * | 2020-06-08 | 2020-08-14 | 吉林大学 | 基于速度梯度的崎岖海底界面的光滑处理方法 |
-
2020
- 2020-12-15 CN CN202011477440.8A patent/CN112505774B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6687659B1 (en) * | 2000-03-24 | 2004-02-03 | Conocophillips Company | Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications |
CN103823239A (zh) * | 2013-10-13 | 2014-05-28 | 中国石油集团西北地质研究所 | 频率域优化混合交错网格有限差分正演模拟方法 |
CN103698814A (zh) * | 2013-12-31 | 2014-04-02 | 中国海洋石油总公司 | 一种用于变密度声波方程的混合吸收边界条件的实现方法 |
CN104237944A (zh) * | 2014-10-09 | 2014-12-24 | 王兵 | 一种适用于交错网格有限差分的全吸收pml方法 |
US20160238723A1 (en) * | 2015-02-13 | 2016-08-18 | Valeriy Brytik | Efficient and Stable Absorbing Boundary Condition In Finite-Difference Calculations |
CN105447225A (zh) * | 2015-11-06 | 2016-03-30 | 中国海洋大学 | 一种应用于声波有限差分数值模拟的组合吸收边界条件 |
CN108073732A (zh) * | 2016-11-10 | 2018-05-25 | 中国石油化工股份有限公司 | 获得稳定的近最佳匹配层吸收边界条件的方法 |
US20190115002A1 (en) * | 2017-10-16 | 2019-04-18 | The Hong Kong University Of Science And Technology | Sound absorber with stair-stepping structure |
CN109188517A (zh) * | 2018-09-03 | 2019-01-11 | 中国海洋大学 | 基于Higdon余弦型加权的混合吸收边界条件方法 |
CN109725351A (zh) * | 2018-12-18 | 2019-05-07 | 中国石油天然气集团有限公司 | 一种3d弹性波混合吸收边界条件的确定方法、装置及系统 |
CN111538083A (zh) * | 2020-06-08 | 2020-08-14 | 吉林大学 | 基于速度梯度的崎岖海底界面的光滑处理方法 |
Non-Patent Citations (2)
Title |
---|
靳双双 等: "基于改进混合吸收边界条件的弹性波方程正演模拟", 《2018年中国地球科学联合学术年会论文集(三十)——专题61:应用地球物理学前沿、专题62:地球重力场及其地学应用》 * |
高刚 等: "完全匹配层人工边界条件中的衰减因子分析", 《石油物探》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112505774B (zh) | 2022-12-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Van Renterghem et al. | Efficient FDTD-PE model for sound propagation in situations with complex obstacles and wind profiles | |
Van Renterghem | Efficient outdoor sound propagation modeling with the finite-difference time-domain (FDTD) method: a review | |
Redondo et al. | Time domain simulation of sound diffusers using finite-difference schemes | |
CN104749628A (zh) | 一种基于弥散黏滞性波动方程的吸收边界反射方法 | |
Yevick et al. | Nonlocal boundary conditions for finite-difference parabolic equation solvers | |
CN112505774B (zh) | 一种地震声波数值模拟中的组合边界方法 | |
Habault et al. | Active control in an anechoic room: Theory and first simulations | |
CN104237944B (zh) | 一种适用于交错网格有限差分的全吸收pml方法 | |
CN105447225A (zh) | 一种应用于声波有限差分数值模拟的组合吸收边界条件 | |
Fang et al. | An unsplit complex frequency-shifted perfectly matched layer for second-order acoustic wave equations | |
CN208473013U (zh) | 一种新型吸声尖劈 | |
Larner et al. | The prediction of the diffuse field sound absorption of perforated panel systems | |
Song et al. | An optimization method for material sound absorption performance based on surrogate model | |
CN108505641A (zh) | 一种新型吸声尖劈 | |
Liu et al. | Study on energy propagation and noise radiation in plates containing the array of acoustic black holes | |
Gao et al. | Application of the Reflectionless Discrete Perfectly Matched Layer for Acoustic Wave Simulation | |
Bonomo et al. | Modeling the acoustic scattering from axially symmetric fluid, elastic, and poroelastic objects due to nonsymmetric forcing using COMSOL Multiphysics | |
Zhu et al. | Waveform Prediction of Underwater Acoustic Propagation Using Finite-difference Time-domain Method | |
CN106094038A (zh) | 适用于tti介质的频率域有限元全吸收pml方法 | |
Kim et al. | Long range propagation modeling of offshore wind turbine construction noise using Finite Element and Parabolic Equation models | |
Min et al. | A coherent image source method for sound prediction in long spaces with a sound absorbent ceiling | |
Zhao et al. | Studies on effective and stable absorbing boundary conditions in ultrasonic wave modeling | |
Jin et al. | Effects of the Rayleigh Secular Function on Time-Harmonic Asymptotic Solutions Due to Horizontal Vibration Sources | |
Wang et al. | Theoretical and quantitative evaluation of hybrid PML-ABCs for seismic wave simulation | |
Yang et al. | An improved convolution perfectly matched layer for elastic second-order wave equation |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |