CN116879952B - 点源弹性波海底反射系数的计算方法、存储介质和设备 - Google Patents
点源弹性波海底反射系数的计算方法、存储介质和设备 Download PDFInfo
- Publication number
- CN116879952B CN116879952B CN202311139493.2A CN202311139493A CN116879952B CN 116879952 B CN116879952 B CN 116879952B CN 202311139493 A CN202311139493 A CN 202311139493A CN 116879952 B CN116879952 B CN 116879952B
- Authority
- CN
- China
- Prior art keywords
- reflection coefficient
- point source
- submarine
- elastic wave
- calculation
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000004590 computer program Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000010587 phase diagram Methods 0.000 description 4
- 238000004378 air conditioning Methods 0.000 description 3
- HOWHQWFXSLOJEF-MGZLOUMQSA-N systemin Chemical compound NCCCC[C@H](N)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC(O)=O)C(=O)OC(=O)[C@@H]1CCCN1C(=O)[C@H]1N(C(=O)[C@H](CC(O)=O)NC(=O)[C@H](CCCN=C(N)N)NC(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H]2N(CCC2)C(=O)[C@H]2N(CCC2)C(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H](CCC(N)=O)NC(=O)[C@@H](NC(=O)[C@H](C)N)C(C)C)CCC1 HOWHQWFXSLOJEF-MGZLOUMQSA-N 0.000 description 3
- 108010050014 systemin Proteins 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Remote Sensing (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Geology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种点源弹性波海底反射系数的计算方法、存储介质和设备，属于地球物理探测、声学探测领域，所述方法包括初始化海底反射系数的计算精度与计算范围；对参数空间进行离散化并获取点源弹性波海底反射系数；联立等效方程与传统点源弹性波海底反射系数计算方程；对等效方程待定系数进行求解；得到第一步精度下的点源海底反射系数的简洁表达式；利用得到的表达式，在给定的计算精度范围内进行点源弹性波海底反射系数计算。本发明方法在保证计算精度的前提下，规避了传统点源海底弹性波反射系数计算过程复杂度较高、效率低下的问题，实现了点源反射系数的简洁、高效计算，推动了点源反射系数的实际应用。
Description
技术领域
本发明属于地球物理探测、声学探测领域。具体地涉及一种点源弹性波海底反射系数的计算方法、存储介质和设备。
背景技术
海底弹性波反射系数蕴含了丰富的地震波场（或声场）与海底沉积物弹性参数信息，海底反射系数的计算是海底反射波场分析与海底弹性参数反演的基础。目前，海底反射系数的计算方法大多以平面波假设为基础，而震源（声源）大多以等效点源为主；这使得以平面波为基础的计算方法只能在远场、高频、小入射角时对海底反射系数进行较为准确地描述。平面波假设无法准确描述海底反射系数的频变与相位特征，同时也无法恰当给出大入射角情况下的反射系数振幅信息。平面波反射系数虽计算复杂度低，但平面波假设的局限性限制了其在上述场景下的适用性。点源反射系数计算方法虽能够规避平面波假设的局限性，但其计算过程复杂度较高、效率低下，这阻碍了点源反射系数的实际应用。因此，需针对点源海底反射系数计算复杂且效率低下的问题，对原有计算方法进行改进。在保证计算精度的前提下给出简洁、高效的点源弹性波海底反射系数计算方法，以满足海底反射波场分析与海底弹性参数反演的需求。
发明内容
为解决上述点源海底弹性波反射系数计算过程复杂度较高、效率低下的问题，本发明提出了一种点源弹性波海底反射系数的计算方法、存储介质和设备，所述方法将点源海底反射系数的求解转化为等效方程待定系数的求解，通过预先对等效方程待定系数求解获得简洁的点源海底反射系数表达式，进而利用该表达式进行高效的点源弹性波海底反射系数计算。
本发明是通过如下技术方案实现的：
一种点源弹性波海底反射系数的计算方法，所述方法具体包含以下步骤：
第一步、初始化海底反射系数的计算精度与计算范围；
第二步、对参数空间进行离散化并获取点源弹性波海底反射系数；
第三步、联立等效方程与传统点源弹性波海底反射系数计算方程；
第四步、对第三步中等效方程待定系数进行求解；
第五步、得到第一步精度下的点源海底反射系数的简洁表达式；
第六步、利用第五步得到的表达式，在第一步给定的计算精度范围内进行点源弹性波海底反射系数计算。
进一步地，第二步中对参数空间进行离散化得到，对应的
点源弹性波海底反射系数记为；为各参数的离散点，以第一步设定
的精度为间隔均匀分布于初始化的参数空间计算范围内；为由传统点源弹性波海底反射
系数计算方法得到的海底反射系数。传统点源弹性波海底反射系数计算方法包括反射球面
波积分法、波动方程法、反射率计算法。
进一步地，第三步中等效方程与传统点源弹性波海底反射系数计算方程通过下式进行联立：
；
式（1）中表示等效方程表达式，为等效方程待定系数，为海底反射
系数， 表示对变量进行赋值，为由传统点源弹性波海底反射系数计算方法得到的海
底反射系数；在参数空间范围内进行联立，得到：
；
R为点源弹性波海底反射系数记。
进一步地，第四步中对等效方程待定系数进行求解，通过以下算式完成：
；
式中，为由高斯核函数构成的对称方阵，其中，R为点源弹性波海底反射系数记。
进一步地，第五步中得到点源海底反射系数的简洁表达式如下：
；
式中，，为等效方程待定系数。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序适用于由处理器加载并执行所述点源弹性波海底反射系数的计算方法。
一种计算机设备，所述设备包括存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，运行所述点源弹性波海底反射系数的计算方法。
本发明与现有技术相比本发明的有益效果：
本发明提出一种点源弹性波海底反射系数的计算方法，将点源海底反射系数的求解转化为等效方程待定系数的求解，通过预先对等效方程待定系数求解获得简洁的点源海底反射系数表达式，本发明计算方法所需计算时间大幅缩短，极大提高了计算效率，在保证计算精度的前提下规避了点源海底弹性波反射系数计算过程复杂度较高、效率低下的问题，推动了点源反射系数的实际应用。
附图说明
图1为本发明计算方法流程图；
图2为离散化参数图；
图3为海底反射系数的振幅图；
图4为海底反射系数的相位图；
图5为等效方程待定系数的振幅图；
图6为等效方程待定系数的相位图；
图7为不同纵波速度情况下的反射系数图，其中，A为振幅图，B为相位图；
图8为不同横波速度情况下的反射系数图，其中，A为振幅图，B为相位图；
图9为不同密度情况下的反射系数图，其中，A为振幅图，B为相位图；
图10为不同频率情况下的反射系数，其中，A为振幅图，B为相位图；
图11 为不同传播距离情况下的反射系数，其中，A为振幅图，B为相位图。
具体实施方式
下面对本发明作进一步详细描述，所提出的实施例仅是本申请的一部分方案而不是全部实施例。基于本发明中的实施例，相关技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例都属于本发明的保护范围。
一种点源弹性波海底反射系数的计算方法，如图1所示，主要包括以下步骤：
第一步、初始化海底反射系数的计算精度与计算范围。
本实施例中通过对输入参数（包括：海底弹性参数、频率、传播距离）进行初始化来
约束海底反射系数的计算精度与计算范围。参数空间计算范围用表示，，计算精
度用表示，。其中，分别
表示海底纵波速度、海底横波速度、海底密度、频率和传播距离，各变量均采用国际单位制。
第二步、对参数空间进行离散化并获取点源弹性波海底反射系数。
按照步骤1）的初始化参数对输入参数进行离散化，得到，
显然，本实施例中离散样点数，参数空间包含海底纵波速度、海底横波速度、海底密
度、频率和传播距离。本实施例离散参数如图2所示，图中横坐标范
围为1~675表示离散样点编号，与右侧坐标轴1450~1550的数值范围对应的粗黑实线表示海
底纵波速度，与右侧坐标轴1500~1600的数值范围对应的虚线表示海底密度，与左侧坐
标轴100~200的数值范围对应的点划线表示海底横波速度，与左侧坐标轴90~110的数值
范围对应的细实线表示传播距离，与左侧坐标轴40~60的数值范围对应的虚线表示频率。
利用上述离散化参数，由传统点源弹性波海底反射系数计算方法可得到的海底反
射系数，传统点源弹性波海底反射系数计算方法包括反射球面波积分法、波动方程法、反射
率计算法，各种方法得到的计算结果是相互吻合的。本实施例在0°~70°的入射角范围内，得
到的与上述离散化参数对应的海底反射系数见图3、图4；图3表示海底
反射系数的模，图4表示海底反射系数的相位。
第三步、联立等效方程与传统点源弹性波海底反射系数计算方程；
通过变量赋值将传统点源弹性波海底反射系数计算方程与等效方程进行联立：
；
式（5）中表示等效方程表达式，为等效方程待定系数，为海底反射
系数， 表示对变量进行赋值。在参数空间范围内进行联立，得到：
；
本实施例将式写为，其中为由高斯核函数构成的对称方阵，其中。
第四步、对步骤3）中等效方程待定系数进行求解，得到等效方程待定系数，
本实施例计算得到的等效方程待定系数如图5、图6所示，其中图5表示等效方程待
定系数的模，图6表示等效方程待定系数的相位。
第五步、得到第一步精度下的点源海底反射系数的简洁表达式：
；
式中，，为高斯函数。
第六步、利用第五步得到的简洁表达式，在第一步给定的计算精度范围内进行点源弹性波海底反射系数计算。
如表1所示，本实施例对计算范围中心值（数值1）及中心偏离值（数值2、数值3）进行点源弹性波海底反射系数计算，计算结果见图7~图11，图中虚线表示平面波计算结果，实线表示传统点源弹性波海底反射系数计算方法的计算结果，“+”表示本发明计算方法的计算结果。图7为不同纵波速度情况下的反射系数，即纵波速度分别为“数值1”、“数值2”、“数值3”，其它参数均为“数值1”时的计算结果；图8为不同横波速度情况下的反射系数，即横波速度分别为“数值1”、“数值2”、“数值3”，其它参数均为“数值1”时的计算结果；图9~图11依次类推。显然，本发明计算结果与传统点源弹性波海底反射系数计算方法的计算结果相吻合；同时平面波计算结果具有较大误差，无法准确描述海底反射系数的频变与相位特征，同时也无法恰当给出大入射角情况下的反射系数振幅信息。平面波假设对比式（7）和反射球面波积分法计算时间发现，本发明计算方法所需计算时间大幅缩短，极大提高了计算效率。
从上述结果可以看出，本发明所提出的点源弹性波海底反射系数的高效计算方法不仅具有更低的计算复杂度，更高的计算效率，还可达到与现有计算方法等同的精度，避免了平面波计算误差带来的局限性，推动了点源反射系数的实际应用。
表1 输入参数表
Claims (4)
1.一种点源弹性波海底反射系数的计算方法，其特征在于，所述方法具体包含以下步骤：
第一步、初始化海底反射系数的计算精度与计算范围；
第二步、对参数空间进行离散化并获取点源弹性波海底反射系数；
第三步、联立等效方程与传统点源弹性波海底反射系数计算方程；所述的等效方程与传统点源弹性波海底反射系数计算方程通过下式进行联立：
f(W,m)＝r←r′ (1)
式(1)中f(w,m)＝r表示等效方程表达式，W为等效方程待定系数，r为海底反射系数，←表示对变量进行赋值，r′为由传统点源弹性波海底反射系数计算方法得到的海底反射系数；在参数空间范围内进行联立，得到：
f(W,M)＝R (2)
式中，R为点源弹性波海底反射系数；
第四步、对第三步中等效方程待定系数进行求解，通过以下算式完成：
W＝G-1R (3)
式中，G为由高斯核函数构成的对称方阵其中i,j∈{k|k＞0&k∈N}，R为点源弹性波海底反射系数；
第五步、得到第一步精度下的点源海底反射系数的简洁表达式，如下：
r＝g(m)·W (4)
式中，W为等效方程待定系数；
第六步、利用第五步得到的表达式，在第一步给定的计算精度范围内进行点源弹性波海底反射系数计算。
2.根据权利要求1所述的一种点源弹性波海底反射系数的计算方法，其特征在于，第二步中对参数空间进行离散化得到M＝[m1 m2…mn]T，对应的点源弹性波海底反射系数记为R＝[r′1 r′2…r′n]T；m为各参数的离散点，以第一步设定的精度为间隔均匀分布于初始化的参数空间计算范围内；r′为由传统点源弹性波海底反射系数计算方法得到的海底反射系数。
3.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序适用于由处理器加载并执行权利要求1-2任何一项所述点源弹性波海底反射系数的计算方法。
4.一种计算机设备，所述设备包括存储器和处理器，所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，运行权利要求1-2任何一项所述点源弹性波海底反射系数的计算方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311139493.2A CN116879952B (zh) | 2023-09-06 | 2023-09-06 | 点源弹性波海底反射系数的计算方法、存储介质和设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311139493.2A CN116879952B (zh) | 2023-09-06 | 2023-09-06 | 点源弹性波海底反射系数的计算方法、存储介质和设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116879952A CN116879952A (zh) | 2023-10-13 |
| CN116879952B true CN116879952B (zh) | 2023-11-21 |
Family
ID=88260870
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311139493.2A Active CN116879952B (zh) | 2023-09-06 | 2023-09-06 | 点源弹性波海底反射系数的计算方法、存储介质和设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116879952B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003087878A8 (en) * | 2002-04-16 | 2004-11-11 | Shell Int Research | Marine time-lapse seismic surveying |
| WO2013169117A1 (en) * | 2012-05-07 | 2013-11-14 | Fugro Norway As | Method and apparatus for sea current aided, enhanced azimuth seismic data acquisition |
| CN104570072A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种粘弹性介质中的球面pp波反射系数建模方法 |
| CN105629301A (zh) * | 2016-03-29 | 2016-06-01 | 中国地质大学（北京） | 薄层弹性波反射系数快速求解方法 |
| CN107576388A (zh) * | 2017-08-22 | 2018-01-12 | 哈尔滨工程大学 | 一种浅海信道下三维结构声源辐射声场预报方法 |
| CN108594302A (zh) * | 2018-07-26 | 2018-09-28 | 广州海洋地质调查局 | 一种地震子波的提取方法及处理终端 |
| CN116520431A (zh) * | 2023-07-03 | 2023-08-01 | 自然资源部第一海洋研究所 | 浅表海底沉积物宽频分层声速结构构建方法、装置和介质 |
| CN116660996A (zh) * | 2023-07-24 | 2023-08-29 | 山东科技大学 | 一种基于深度学习的漂流式浅海局部地声参数预测方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6549853B2 (en) * | 2001-03-22 | 2003-04-15 | Council Of Scientific & Industrial Research | Method for determining seafloor roughness using multibeam echosounder |
| US6833698B2 (en) * | 2002-05-15 | 2004-12-21 | Chevrontexaco U.S.A. Inc. | Methods of decoupling diffusion effects from relaxation times to determine properties of porous media containing fluids |
-
2023
- 2023-09-06 CN CN202311139493.2A patent/CN116879952B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003087878A8 (en) * | 2002-04-16 | 2004-11-11 | Shell Int Research | Marine time-lapse seismic surveying |
| WO2013169117A1 (en) * | 2012-05-07 | 2013-11-14 | Fugro Norway As | Method and apparatus for sea current aided, enhanced azimuth seismic data acquisition |
| CN104570072A (zh) * | 2013-10-16 | 2015-04-29 | 中国石油化工股份有限公司 | 一种粘弹性介质中的球面pp波反射系数建模方法 |
| CN105629301A (zh) * | 2016-03-29 | 2016-06-01 | 中国地质大学（北京） | 薄层弹性波反射系数快速求解方法 |
| CN107576388A (zh) * | 2017-08-22 | 2018-01-12 | 哈尔滨工程大学 | 一种浅海信道下三维结构声源辐射声场预报方法 |
| CN108594302A (zh) * | 2018-07-26 | 2018-09-28 | 广州海洋地质调查局 | 一种地震子波的提取方法及处理终端 |
| CN116520431A (zh) * | 2023-07-03 | 2023-08-01 | 自然资源部第一海洋研究所 | 浅表海底沉积物宽频分层声速结构构建方法、装置和介质 |
| CN116660996A (zh) * | 2023-07-24 | 2023-08-29 | 山东科技大学 | 一种基于深度学习的漂流式浅海局部地声参数预测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 频率-波数域递推计算海底多分量地震记录中的反射系数;顾汉明, 王家映, 朱光明;地球物理学报(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116879952A (zh) | 2023-10-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109858523B (zh) | 一种神经网络和射线理论的浅海声速剖面反演方法 | |
| CN110750934A (zh) | 深海弹性结构与环境耦合声辐射预报方法 | |
| CN111830493B (zh) | 水中目标收发分置中高频声目标强度预报系统及方法 | |
| CN113687307B (zh) | 低信噪比及混响环境下的自适应波束成形方法 | |
| CN116520431B (zh) | 浅表海底沉积物宽频分层声速结构构建方法、装置和介质 | |
| Wei et al. | Simulation of ultrasound beam formation of baiji (Lipotes vexillifer) with a finite element model | |
| WO2022226856A1 (zh) | 浅海多层海底地声参数反演方法、装置、计算机设备及存储介质 | |
| JP6718098B2 (ja) | 位置推定装置及び方法 | |
| CN107016159A (zh) | 本征值确定方法及装置 | |
| CN111401180B (zh) | 神经网络识别模型训练方法、装置、服务器及存储介质 | |
| CN115630446A (zh) | 用于潜航器结构低频辐射噪声实时仿真的快速推演方法 | |
| CN116879952B (zh) | 点源弹性波海底反射系数的计算方法、存储介质和设备 | |
| CN116660996B (zh) | 一种基于深度学习的漂流式浅海局部地声参数预测方法 | |
| Sessarego et al. | Acoustic scattering by an elastic spherical shell near the seabed | |
| CN116359990A (zh) | 一种基于变密度声波方程的深度域反演方法 | |
| CN112729532B (zh) | 浅海运动声源近场传递函数获取方法、装置、设备及介质 | |
| CN112434467B (zh) | 一种水下旋转体目标低目标强度外形设计方法 | |
| CN115469314A (zh) | 一种均匀圆环阵稳健水下目标方位跟踪方法及系统 | |
| Martinelli | An application of the level set method to underwater acoustic propagation | |
| CN114974298A (zh) | 一种基于深度学习的海洋声源定位及海洋参数反演方法 | |
| CN115292886A (zh) | 一种基于声场干涉现象的弹性海底声源运动参数估计系统及估计方法 | |
| CN114841892A (zh) | 一种基于全连接网络的稀疏导波数据恢复方法 | |
| CN113721245A (zh) | 海底水平阵阵形校正方法及处理器 | |
| CN116859460B (zh) | 一种适用于弹性波近场反射的海底物性参数反演方法 | |
| CN114444353A (zh) | 基于有限元方法的二维声场模型的声场获取方法 |
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 |