CN116859460A - 一种适用于弹性波近场反射的海底物性参数反演方法 - Google Patents
一种适用于弹性波近场反射的海底物性参数反演方法 Download PDFInfo
- Publication number
- CN116859460A CN116859460A CN202311126579.1A CN202311126579A CN116859460A CN 116859460 A CN116859460 A CN 116859460A CN 202311126579 A CN202311126579 A CN 202311126579A CN 116859460 A CN116859460 A CN 116859460A
- Authority
- CN
- China
- Prior art keywords
- inversion
- plane wave
- result
- submarine
- physical property
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000000704 physical effect Effects 0.000 title claims abstract description 28
- 230000009897 systematic effect Effects 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 5
- 239000003086 colorant Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002310 reflectometry Methods 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种适用于弹性波近场反射的海底物性参数反演方法，属于地球物理探测、声学探测领域，所述方法为初始化待反演参数空间的范围及离散化间隔；参数空间离散化并计算各离散点近场反射系数；获取各离散点的平面波反演误差；求解任意平面波反演结果及其系统误差之间的数学关系；输入待反演的观测数据并对待反演的数据进行平面波反演；利用所述的数学关系计算所述反演误差；对反演结果进行校正；输出校正后的反演结果。本发明方法能够很好地与海底物性参数真实值吻合，可正确指示海底物性参数的分布情况；最终得到弹性波近场反射情况下的可靠反演结果。
Description
技术领域
本发明属于地球物理探测、声学探测领域，具体地涉及一种适用于弹性波近场反射的海底物性参数反演方法。
背景技术
海底物性参数对海洋资源探测、地下储气监测、海洋工程建设、海底环境监测、滑坡灾害评估等具有重要意义。海底物性参数决定了海底地震波（或声波）反射系数随入射角的变化关系，因此利用海底反射系数随入射角的变化关系可以开展海底沉积物物性参数的反演。目前，现有反演方法主要以平面波理论为基础，而在获取不同入射角反射系数的过程中，所采用的场源大多以等效点源为主。虽然平面波理论计算简洁高效，但平面波理论只能在远场情况下时对点源海底反射系数进行较为准确地描述；这使得现有反演方法在利用点源近场反射数据进行反演时会产生较大的系统误差。
发明内容
为解决现有反演方法在利用点源近场反射数据进行反演时会产生较大系统误差的问题，本发明提出了一种适用于弹性波近场反射的海底物性参数反演方法。通过对平面波反演结果及其系统误差的理论关系进行求解，进而实现平面波反演误差的校正，破平面波假设的限制，最终得到弹性波近场反射情况下的可靠反演结果。
本发明是通过如下技术方案实现的：
一种适用于弹性波近场反射的海底物性参数反演方法，所述方法具体包含以下步骤：
第一步、初始化待反演参数空间的范围及离散化间隔；
第二步、参数空间离散化并计算各离散点近场反射系数；
第三步、获取各离散点的平面波反演误差；
第四步、求解任意平面波反演结果及其系统误差之间的数学关系；
第五步、输入待反演的观测数据并对待反演的数据进行平面波反演；
第六步、利用第四步的数学关系计算第五步的反演误差；
第七步、利用第六步的计算结果对第五步反演结果进行校正；
第八步、输出第七步校正后的反演结果即为最终反演结果；
进一步地，第三步对第二步中各离散点的平面波反演误差进行计算，通过下式实现：
；
式中，表示各离散点平面波反演误差；为第二步中对参数空间进行离散化得
到的离散化参数，为第二步得到的近场海底反射系数， 中的各个分量由中的相应分量在第二步中通过点源海底反射
系数计算得到：
；
式 中、分别表示平面波反演、点源海底反射系数正演，。
进一步地，第四步求取任意平面波反演结果及其系统误差之间的数学关系，
通过下式完成：
；
式中为式（1）的计算结果，，表示矩阵的逆矩阵，的表达式如下：
。
进一步地，第六步中利用第四步的数学关系计算第五步的反演误差，计算公式如下：
；
式中为第五步中输入的待反演的观测数据，表示第五步中平面波反演结
果的系统误差。
进一步地，第七步中利用第六步的计算结果对第五步反演结果进行校正，通过下式实现：
；
式中为校正后的反演结果即近场反射情况下的海底物性参数反演结果，表示第五步中平面波反演结果的系统误差，表示平面波反演，为第五步中输入
的待反演的观测数据。
本发明与现有技术相比的有益效果：本发明方法通过对平面波反演结果及其系统误差的理论关系进行求解，进而实现平面波反演误差的校正，破平面波假设的限制，能够很好地与海底物性参数真实值吻合，可正确指示海底物性参数的分布情况；最终得到弹性波近场反射情况下的可靠反演结果。可满足近场反射情况下海底物性参数反演的需求。
附图说明
图1为本发明反演方法流程图；
图2为离散化参数图；
图3为海底反射系数图；
图4为各离散点平面波反演误差图；
图5为待反演观测数据图；
图6为平面波反演结果图；
图7 平面波反演结果的系统误差图；
图8 近场反射情况下的海底物性参数反演结果。
具体实施方式
下面对本发明作进一步详细描述，所提出的实施例仅是本申请的一部分方案而不是全部实施例。基于本发明中的实施例，相关技术人员在没有做出创造性劳动的前提下所获得的所有其它实施例都属于本发明的保护范围。
一种适用于弹性波近场反射的海底物性参数反演方法，如图1所示，主要包括以下步骤：
第一步、初始化待反演参数空间的范围及离散化间隔；
本实施例中待反演参数包括海底纵波速度、横波速度、密度，待反演参数空
间用表示为，，离散化间隔
用表示，。其中，各变量均采用国际单位制。
第二部、参数空间离散化并计算各离散点近场反射系数；
按照步骤1）的初始化参数对输入参数空间进行离散化，得到，其中，，。本例中离散参数如图2所示，
图中横坐标范围为1~2535表示离散样点号，纵坐标从上到下依次为纵波速度、横波速度、密度。
利用上述离散化参数，可通过反射球面波积分法、波动方程法或反射率计算法计
算海底近场反射系数，本实施例采用球面波积分法在1°~65°的入射角范围内，得到的与上
述离散化参数对应的海底近场反射系数，见图3，横坐标范围为1~2535
表示离散样点号，纵坐标范围为1°~65°表示反射系数入射角，不同灰度颜色表示反射系数
值得大小。
第三步、获取各离散点的平面波反演误差；
对各离散点的平面波反演误差进行计算，通过下式实现：
；
式中中表示各离散点平面波反演误差，本例计算结果见图4，图中横坐标表示
离散样点号，纵坐标从上到下依次为纵波速度误差、横波速度误差、密度误差；为步骤2）
中对参数空间进行离散化得到的离散化参数，为步骤2）得到的
近场海底反射系数， 中的各个分量由中的相应分量在步骤2）中
通过点源海底反射系数计算得到：
；
式 中、分别表示平面波反演、点源海底反射系数正演，。
第四步、求解任意平面波反演结果及其系统误差之间的数学关系；
求取任意平面波反演结果及其系统误差之间的数学关系，通过下式完成：
；
式中为式计算结果，，为高斯
函数；表示矩阵的逆矩阵，的表达式如下：
；
至此，我们完成了本实施例任意平面波反演结果及其系统误差之间的数学关系求解。
第五步、输入待反演的观测数据并对待反演的数据进行平面波反演；
本实施例中输入待反演观测数据如图5所示，图中横坐标编号1~600表示600组待反演输入数据，横坐标1°~65°表示输入数据的入射角范围，灰度颜色表示反射系数值。
由待反演数据进行平面波反演得到的反演结果如图6所示，图中横坐标编号1~600与图5中横坐标一一对应，纵坐标从上到下依次为纵波速度、横波速度、密度；实线表示平面波反演结果，虚线表示观测数据所对应的海底物性参数真实值。显然，平面波反演结果与海底物性参数真实值之间存在较大误差，无法正确指示海底物性参数的分布情况。
第六步、利用第四步的数学关系计算第五步的反演误差；
利用第四步的数学关系计算第五步的反演误差，计算公式如下：
；
式中为第五步中输入的待反演的观测数据，表示第五步中平面波反演
结果的系统误差。本例中平面波反演结果的系统误差计算结果如图7所示，图中横坐标
意义与图6一致，纵坐标从上到下分别表示纵波速度误差、横波速度误差、密度误差。
第七步、利用第六步的计算结果对第五步反演结果进行校正；
利用第六步的计算结果对第五步反演结果进行校正，通过下式实现：
；
式中为校正后的反演结果即近场反射情况下的海底物性参数反演结果，其
它各符号与前述步骤的相应符号含义相同。本实施例中校正后的反演结果如图8所示，图中
纵横坐标意义同图6；实线表示本发明反演结果，虚线表示观测数据所对应的海底物性参数
真实值。显然，本发明结果能够很好地与海底物性参数真实值吻合，可正确指示海底物性参
数的分布情况；也就是说本发明克服了现有技术的不足、可满足近场反射情况下海底物性
参数反演的需求。
第八步输出第七步校正后的反演结果即为最终反演结果；
输出，即为最终反演结果。
从上述结果可以看出，以平面波理论为基础的反演方法在利用海底近场反射数据进行反演时会产生较大的系统误差的问题。本发明所提出的反演方法能够有效对平面波反演误差进行校正，突破了现有平面波反演方法的近场反演局限性，得到可靠的近场反射情况下的反演结果。
Claims (5)
1.一种适用于弹性波近场反射的海底物性参数反演方法，其特征在于，所述方法具体包含以下步骤：
第一步、初始化待反演参数空间的范围及离散化间隔；
第二步、参数空间离散化并计算各离散点近场反射系数；
第三步、获取各离散点的平面波反演误差；
第四步、求解任意平面波反演结果及其系统误差之间的数学关系；
第五步、输入待反演的观测数据并对待反演的数据进行平面波反演；
第六步、利用第四步的数学关系计算第五步的反演误差；
第七步、利用第六步的计算结果对第五步反演结果进行校正；
第八步、输出第七步校正后的反演结果即为最终反演结果。
2.根据权利要求1所述的一种适用于弹性波近场反射的海底物性参数反演方法，其特征在于，第三步对第二步中各离散点的平面波反演误差进行计算，通过下式实现：
；
式中，表示各离散点平面波反演误差；/>为第二步中对参数空间进行离散化得到的离散化参数/>，/>为第二步得到的近场海底反射系数，/>中的各个分量由/>中的相应分量在第二步中通过点源海底反射系数计算得到：
；
式中、/>分别表示平面波反演、点源海底反射系数正演，/>。
3.根据权利要求2所述的一种适用于弹性波近场反射的海底物性参数反演方法，其特征在于，第四步求取任意平面波反演结果及其系统误差/>之间的数学关系，通过下式完成：
；
式中为式（1）的计算结果，/>，/>表示矩阵/>的逆矩阵，/>的表达式如下：
。
4.根据权利要求3所述的一种适用于弹性波近场反射的海底物性参数反演方法，其特征在于，第六步中利用第四步的数学关系计算第五步的反演误差，计算公式如下：
；
式中为第五步中输入的待反演的观测数据，/>表示第五步中平面波反演结果的系统误差。
5.根据权利要求1所述的一种适用于弹性波近场反射的海底物性参数反演方法，其特征在于，第七步中利用第六步的计算结果对第五步反演结果进行校正，通过下式实现：
；
式中为校正后的反演结果即近场反射情况下的海底物性参数反演结果，/>表示第五步中平面波反演结果的系统误差，/>表示平面波反演，/>为第五步中输入的待反演的观测数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311126579.1A CN116859460B (zh) | 2023-09-04 | 2023-09-04 | 一种适用于弹性波近场反射的海底物性参数反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311126579.1A CN116859460B (zh) | 2023-09-04 | 2023-09-04 | 一种适用于弹性波近场反射的海底物性参数反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116859460A true CN116859460A (zh) | 2023-10-10 |
| CN116859460B CN116859460B (zh) | 2023-11-21 |
Family
ID=88219408
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311126579.1A Active CN116859460B (zh) | 2023-09-04 | 2023-09-04 | 一种适用于弹性波近场反射的海底物性参数反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116859460B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11966001B1 (en) * | 2023-09-06 | 2024-04-23 | First Institute Of Oceanography, Mnr | Calculation method, storage medium and device for seabed reflection coefficient of point source elastic wave |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140039799A1 (en) * | 2012-08-06 | 2014-02-06 | Christine E. Krohn | Seismic Inversion for Formation Properties and Attentuation Effects |
| CN104597490A (zh) * | 2015-01-28 | 2015-05-06 | 中国石油大学(北京) | 基于精确Zoeppritz方程的多波AVO储层弹性参数反演方法 |
| CN104820218A (zh) * | 2015-05-26 | 2015-08-05 | 中国科学院声学研究所东海研究站 | 一种基于频域自相关的浅海海底单参数反演方法 |
| CN106154276A (zh) * | 2016-06-27 | 2016-11-23 | 西北工业大学 | 基于海底混响和传播损失的深海海底参数反演方法 |
| US20170269246A1 (en) * | 2016-03-16 | 2017-09-21 | John E. Anderson | Method to Estimate and Remove Direct Arrivals From Arrayed Marine Sources |
| CN109521469A (zh) * | 2018-11-16 | 2019-03-26 | 国家海洋局第海洋研究所 | 一种海底沉积物弹性参数的正则化反演方法 |
-
2023
- 2023-09-04 CN CN202311126579.1A patent/CN116859460B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140039799A1 (en) * | 2012-08-06 | 2014-02-06 | Christine E. Krohn | Seismic Inversion for Formation Properties and Attentuation Effects |
| CN104597490A (zh) * | 2015-01-28 | 2015-05-06 | 中国石油大学(北京) | 基于精确Zoeppritz方程的多波AVO储层弹性参数反演方法 |
| CN104820218A (zh) * | 2015-05-26 | 2015-08-05 | 中国科学院声学研究所东海研究站 | 一种基于频域自相关的浅海海底单参数反演方法 |
| US20170269246A1 (en) * | 2016-03-16 | 2017-09-21 | John E. Anderson | Method to Estimate and Remove Direct Arrivals From Arrayed Marine Sources |
| CN106154276A (zh) * | 2016-06-27 | 2016-11-23 | 西北工业大学 | 基于海底混响和传播损失的深海海底参数反演方法 |
| CN109521469A (zh) * | 2018-11-16 | 2019-03-26 | 国家海洋局第海洋研究所 | 一种海底沉积物弹性参数的正则化反演方法 |
Non-Patent Citations (1)
| Title |
|---|
| YANGTING LIU ET AL.: "Seafoor elastic‑parameter estimation based on a regularized AVO inversion method", MARINE GEOPHYSICAL RESEARCH * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11966001B1 (en) * | 2023-09-06 | 2024-04-23 | First Institute Of Oceanography, Mnr | Calculation method, storage medium and device for seabed reflection coefficient of point source elastic wave |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116859460B (zh) | 2023-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116859460B (zh) | 一种适用于弹性波近场反射的海底物性参数反演方法 | |
| CN109858523B (zh) | 一种神经网络和射线理论的浅海声速剖面反演方法 | |
| CN114397643B (zh) | 一种基于超短基线水声定位系统的声线修正方法 | |
| CN110398744A (zh) | 一种基于声信号的海洋温跃层特征参数寻优反演方法 | |
| Guérin et al. | An improved Lagrangian model for the time evolution of nonlinear surface waves | |
| CN109521469B (zh) | 一种海底沉积物弹性参数的正则化反演方法 | |
| CN107016159A (zh) | 本征值确定方法及装置 | |
| Rogers et al. | A study of dissipation of wind-waves by mud at Cassino Beach, Brazil: Prediction and inversion | |
| CN111239819A (zh) | 一种基于地震道属性分析的带极性直接包络反演方法 | |
| CN115951361A (zh) | 基于大数据平台的海上风电桩的智能地形扫测方法 | |
| CN115841076A (zh) | 基于bp神经网络模型的浅海分层海底地声参数反演方法 | |
| CN117198330A (zh) | 一种声源识别方法、系统及电子设备 | |
| CN116879952B (zh) | 点源弹性波海底反射系数的计算方法、存储介质和设备 | |
| CN111273349B (zh) | 一种用于海底浅部沉积层的横波速度提取方法及处理终端 | |
| CN116663428B (zh) | 一种基于神经网络的面波成像联合反演方法 | |
| CN116977842A (zh) | 一种面向浅海分区环境的sar图像水深估计方法、系统、设备及介质 | |
| CN114974298A (zh) | 一种基于深度学习的海洋声源定位及海洋参数反演方法 | |
| CN113221651B (zh) | 利用声传播数据和无监督机器学习的海底沉积分类方法 | |
| Zalesny et al. | Numerical algorithm of data assimilation based on splitting and adjoint equation methods | |
| TANG et al. | Application of LVQ neural network combined with the genetic algorithm in acoustic seafloor classification | |
| CN113534255B (zh) | 一种任意形态不连续面自适应表达的方法 | |
| Martinelli | An application of the level set method to underwater acoustic propagation | |
| Duan et al. | PINNs for Sound Propagation and Sound Speed Field Estimation Simultaneously | |
| Steen et al. | Uncertainty analysis in ship-model resistance test | |
| CN116011356B (zh) | 一种基于ConvLSTM与OpenFOAM数值计算耦合的海底管缆冲刷预测方法 |
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 |