CN116164699B - 一种行星地震仪的主动源方位角自动识别方法 - Google Patents

一种行星地震仪的主动源方位角自动识别方法 Download PDF

Info

Publication number
CN116164699B
CN116164699B CN202310458339.5A CN202310458339A CN116164699B CN 116164699 B CN116164699 B CN 116164699B CN 202310458339 A CN202310458339 A CN 202310458339A CN 116164699 B CN116164699 B CN 116164699B
Authority
CN
China
Prior art keywords
planetary
seismograph
azimuth angle
lander
active source
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
Application number
CN202310458339.5A
Other languages
English (en)
Other versions
CN116164699A (zh
Inventor
张盛泉
游庆瑜
王娟
李少卿
张金海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN202310458339.5A priority Critical patent/CN116164699B/zh
Publication of CN116164699A publication Critical patent/CN116164699A/zh
Application granted granted Critical
Publication of CN116164699B publication Critical patent/CN116164699B/zh
Priority to US18/447,322 priority patent/US11899149B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C1/00—Measuring angles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/01—Measuring or predicting earthquakes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • 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
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Acoustics & Sound (AREA)
  • Manufacturing & Machinery (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种行星地震仪的主动源方位角自动识别方法，包括：依次截取行星地震仪接收到第一条着陆器支撑腿到第N条着陆器支撑腿的主动源初至信号前后各δt时长的三通道原始时间序列；基于行星地震仪布设到行星表面后的俯仰角θ1和翻滚角θ2将三通道原始时间序列转换到水平面；将水平面的时间序列转换到RTZ坐标上，并求震动信号在R方向和T方向的分量最大振幅；基于震动信号在R方向和T方向的分量最大振幅构造目标函数A，在0~360°范围内以预设角度为间隔扫描行星地震仪方位角φ，当目标函数A达到最小时，所对应的行星地震仪方位角φ即为最优估计值，无需人工干预即可实现行星地震仪方位角自动识别。

Description

一种行星地震仪的主动源方位角自动识别方法
技术领域
本发明涉及行星探测技术领域，更具体的说是涉及一种行星地震仪的主动源方位角自动识别方法。
背景技术
布设行星地震仪是探测行星内部结构和行星地下资源的最有效手段，行星地震仪的布设分为无人布设和有人布设两种。其中，无人布设将面临布设后其方位角无法自我感知的困难，必须设法获得仪器自身的方位角才能获得准确的行星地震信号定位和应用，因而如何实现行星地震仪的方位感知至关重要。此外，有人布设的行星地震仪可以通过宇航员用专业设备现场实测仪器的方位角。但无论哪种布设方式，都将面临仪器布设后受到振动、扰动或者不均匀沉降等带来的仪器姿态变化，而这种变化需要有方法能够及时感知，才能确保行星地震仪观测数据的精度和应用效果。
现有技术中，行星地震仪的方位感知尚处空白，且无自动识别算法。对地震仪的方位角识别方法有实地勘测、采用天然源P波极性分析法、波形模拟对比法、噪声互相关法等方式，但无地外行星地震仪主动源法的方位角自动识别方法。
因此，如何提供一种行星地震仪的主动源方位角自动识别方法是本领域技术人员亟需解决的问题。
发明内容
有鉴于此，本发明提供了一种行星地震仪的主动源方位角自动识别方法，无需人工干预即可实现行星地震仪方位角自动识别。
为了实现上述目的，本发明采用如下技术方案：
一种行星地震仪的主动源方位角自动识别方法，包括：
依次截取行星地震仪接收到第一条着陆器支撑腿到第N条着陆器支撑腿的主动源初至信号前后各δt时长的三通道原始时间序列；
基于行星地震仪布设到行星表面后的俯仰角θ1和翻滚角θ2将三通道原始时间序列转换到水平面；
将水平面的时间序列转换到RTZ坐标上，并求震动信号在R方向和T方向的分量最大振幅；
基于震动信号在R方向和T方向的分量最大振幅构造目标函数A，在0~360°范围内以预设角度为间隔扫描行星地震仪方位角φ，当目标函数A达到最小时，所对应的行星地震仪方位角φ即为最优估计值。
优选地，基于行星地震仪布设到行星表面后的俯仰角θ1和翻滚角θ2将三通道原始时间序列转换到水平面，具体计算公式为：
xhi= xi*cosθ1
yhi= yi*cosθ2
其中，xi表示行星地震仪的NS方向截取时间序列，yi表示行星地震仪的EW方向截取时间序列，i=1,2，...，N，N表示着陆器支撑腿的个数，xhi表示行星地震仪的NS方向截取时间序列xi在行星水平面上的分量，yhi表示行星地震仪的EW方向截取时间序列yi在行星水平面上的分量。
优选地，震动信号在R方向和T方向的分量最大振幅计算公式为：
Figure SMS_1
其中，θai为第i条着陆器支撑腿方向与行星地震仪N方向的夹角，Ari和Ati分别表示震动信号在R方向和T方向的分量最大振幅，i=1,2，...，N，N表示着陆器支撑腿的个数。
优选地，目标函数A计算公式为：
Figure SMS_2
其中，i=1,2，...，N，N表示着陆器支撑腿的个数；
行星地震仪方位角φ计算公式为：
φ = θci – θai
其中，i=1,2，...，N，N表示着陆器支撑腿的个数，震源入射方位角θai为第i条着陆器支撑腿方向与行星地震仪N方向的夹角，θci表示第i条着陆器支撑腿触地点的方位角。
优选地，在着陆器多个支撑腿的相同高度处安装固定震动激励装置，通过着陆器的控制器顺次开启和关闭各个震动激励装置，震动信号将由支撑腿传递至行星表面，再传递至行星地震仪内部，并被行星地震仪感知。
本发明具有以下效果：
1.通过主动源方式，可得到比噪声互相关法更高信噪比和更高灵敏度的探测数据，计算得到的方位角误差更小，结果更为可靠。
2.通过主动源方式，无需等待天然行星地震信号，即可自动、快速、高效获得行星地震仪的方位角，从而能够为后续天然行星地震信号的分析处理提供依据。
3.本发明为行星地震仪的月面或火星地面在轨定标提供了一种新途径，通过着陆器支撑腿定期产生的震动信号，为行星地震仪提供工作正常与否的判断依据，并可对其长期运行的性能参数进行定期主动标校。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的一种行星地震仪的主动源方位角自动识别方法流程图。
图2附图为本发明提供的布设于行星表面的行星地震仪与着陆器相对位置关系示意图。
图3附图为本发明提供的布设于行星表面的行星地震仪与着陆器的俯视图。
图4附图为本发明提供的震动激励装置发射波形图与行星地震仪接收波形图。
图5附图为本发明提供的主动源震动信号传播示意图。
其中，1、着陆器，2、支撑腿，3、震动激励装置，4、行星地震仪。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例公开了一种行星地震仪的主动源方位角自动识别方法，如图1所示，包括：
依次截取行星地震仪接收到第一条着陆器支撑腿到第N条着陆器支撑腿的主动源初至信号前后各δt时长的三通道原始时间序列；
基于行星地震仪布设到行星表面后的俯仰角θ1和翻滚角θ2将三通道原始时间序列转换到水平面；
将水平面的时间序列转换到RTZ坐标上，并求震动信号在R方向和T方向的分量最大振幅；
基于震动信号在R方向和T方向的分量最大振幅构造目标函数A，在0~360°范围内以预设角度为间隔扫描行星地震仪方位角φ，当目标函数A达到最小时，所对应的行星地震仪方位角φ即为最优估计值。
如图2所示，本发明在着陆器1多个支撑腿2的相同高度处安装固定震动激励装置3，通过着陆器1的控制器顺次开启和关闭各个震动激励装置3，震动信号将由支撑腿2传递至行星表面，再传递至行星地震仪4内部，并被行星地震仪4感知。该装置可根据实际需要定时开启或随机开启，从而实现行星地震仪方位的感知和校准。需要说明的是，着陆器支撑腿个数不限制，如图2所示，本实施例以四条支撑腿的着陆器为例。
如图3所示，当第一个震动激励装置3工作时，因为振动沿着着陆器1的腿部向地面传播并达到行星地震仪4的距离最短，因而会被行星地震仪4最先感知到。该震源虽然也会沿着着陆器1传播到其它着陆器腿并经地面被行星地震仪4观测到，但其振动会比直接安装震动激励的着陆器腿的到达时间更晚。因此，行星地震仪4记录到的初始到达振动的偏振方向主要与当前震动激励的着陆器支撑腿所在的方位有关，从而能够计算得到采集到的第一条着陆器腿震动信号。以此类推，行星地震仪4可以分别采集得到各个着陆腿的震动信号。
在本实施例中，依次截取行星地震仪接收到第一条着陆器支撑腿到第N条着陆器支撑腿的主动源初至信号前后各δt时长的三通道原始时间序列具体步骤为：
由发射震动信号和接收信号的频率和持续时间特征，截取行星地震仪接收到第一条着陆器支撑腿的主动源初至信号前后各δt时长的三通道原始时间序列，记为An1=｛x1，y1，z1｝，如图4所示，其中，行星地震仪的NS方向截取时间序列为x1，EW方向截取时间序列为y1，垂直方向截取时间序列为z1。假设着陆器有N条支撑腿，则第二、第三到第N条支撑腿主动源震动信号的截取记录数据记为An2、 An3、…、 AnN，第i条支撑腿主动源震动信号的截取记录数据记为Ani=｛xi，yi，zi｝，其中i=1,2，...，N，N表示着陆器支撑腿的个数。
在本实施例中，基于行星地震仪布设到行星表面后的俯仰角θ1和翻滚角θ2将三通道原始时间序列转换到水平面，具体包括：
由行星地震仪内部姿态传感器得到其布设行星表面后的俯仰角θ1和翻滚角θ2，对于截取记录数据Ani=｛xi，yi，zi｝其中i=1,2，...，N，将三通道原始时间序列xi，yi，zi按下列公式转换到水平面，记为xhi，yhi，zhi
Figure SMS_3
通过上述步骤，将第一条着陆器支撑腿到第N条着陆器支撑腿对应的三通道原始时间序列分别转换到水平面。
在本实施例中，将水平面的时间序列转换到RTZ坐标上，并求震动信号在R方向和T方向的分量最大振幅具体包括：
将水平面的时间序列xhi和yhi，其中i=1,2，...，N，以震源入射方位角θai，θai定义为第i条着陆器支撑腿方向（震源方向）与行星地震仪N方向的夹角，按下述公式转换到RTZ坐标上，并求取其最大振幅。
RTZ坐标系定义如下：R方向定义为从产生激励信号的支撑腿触地点沿着行星大圆路径指向行星地震仪的方向，Z方向为垂直水平面向上，T方向垂直于R和Z方向所构成的平面。震动信号在R方向的分量最大振幅记为Ari, 震动信号在T方向的分量最大振幅记为Ati：
Figure SMS_4
通过上述步骤得到从第一条着陆器支撑腿到第N条着陆器支撑腿对应的震动信号在R方向的分量最大振幅和在T方向的分量最大振幅。
在本实施例中，目标函数A为：
Figure SMS_5
如图5所示，行星地震仪方位角可表示为φ = θci – θai，其中i=1,2，...，N，震源入射方位角θai为第i条着陆器支撑腿方向与行星地震仪N方向的夹角。而 θci为已知第i个支撑腿触地点的方位角，可由着陆器方位角和其支撑腿尺寸分布及行星地震仪布设位置计算获得。
在0~360°范围内以1°（或者更小）为间隔扫描φ ，当A值达到最小时，所对应的行星地震仪绝对方位角φ即为最优估计值，最终实现行星地震仪方位角的自动感知。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (5)

1.一种行星地震仪的主动源方位角自动识别方法，其特征在于，包括：
依次截取行星地震仪接收到第一条着陆器支撑腿到第N条着陆器支撑腿的主动源初至信号前后各δt时长的三通道原始时间序列；
基于行星地震仪布设到行星表面后的俯仰角θ1和翻滚角θ2，将三通道原始时间序列转换到水平面；
将水平面的时间序列转换到RTZ坐标上，并求震动信号在R方向和T方向的分量最大振幅；
基于震动信号在R方向和T方向的分量最大振幅构造目标函数A，在0~360°范围内以预设角度为间隔扫描行星地震仪方位角φ，当目标函数A达到最小时，所对应的行星地震仪方位角φ即为最优估计值。
2.根据权利要求1所述的一种行星地震仪的主动源方位角自动识别方法，其特征在于，基于行星地震仪布设到行星表面后的俯仰角θ1和翻滚角θ2将三通道原始时间序列转换到水平面，具体计算公式为：
xhi= xi*cosθ1
yhi= yi*cosθ2
其中，xi表示行星地震仪的NS方向截取时间序列，yi表示行星地震仪的EW方向截取时间序列，i=1,2，...，N，N表示着陆器支撑腿的个数，xhi表示行星地震仪的NS方向截取时间序列xi在行星水平面上的分量，yhi表示行星地震仪的EW方向截取时间序列yi在行星水平面上的分量。
3.根据权利要求2所述的一种行星地震仪的主动源方位角自动识别方法，其特征在于，震动信号在R方向和T方向的分量最大振幅计算公式为：
Figure QLYQS_1
；
其中，θai为第i条着陆器支撑腿方向与行星地震仪N方向的夹角，Ari和Ati分别表示震动信号在R方向和T方向的分量最大振幅，i=1,2，...，N，N表示着陆器支撑腿的个数。
4.根据权利要求3所述的一种行星地震仪的主动源方位角自动识别方法，其特征在于，目标函数A计算公式为：
Figure QLYQS_2
；
其中，i=1,2，...，N，N表示着陆器支撑腿的个数；
行星地震仪方位角φ计算公式为：
φ = θci – θai
其中，i=1,2，...，N，N表示着陆器支撑腿的个数，震源入射方位角θai为第i条着陆器支撑腿方向与行星地震仪N方向的夹角，θci表示第i条着陆器支撑腿触地点的方位角。
5.根据权利要求1所述的一种行星地震仪的主动源方位角自动识别方法，其特征在于，在着陆器多个支撑腿的相同高度处安装固定震动激励装置，通过着陆器的控制器顺次开启和关闭各个震动激励装置，震动信号将由支撑腿传递至行星表面，再传递至行星地震仪内部，并被行星地震仪感知。
CN202310458339.5A 2023-04-26 2023-04-26 一种行星地震仪的主动源方位角自动识别方法 Active CN116164699B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202310458339.5A CN116164699B (zh) 2023-04-26 2023-04-26 一种行星地震仪的主动源方位角自动识别方法
US18/447,322 US11899149B1 (en) 2023-04-26 2023-08-10 Method for automatically identifying active source azimuth of planetary seismometer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310458339.5A CN116164699B (zh) 2023-04-26 2023-04-26 一种行星地震仪的主动源方位角自动识别方法

Publications (2)

Publication Number Publication Date
CN116164699A CN116164699A (zh) 2023-05-26
CN116164699B true CN116164699B (zh) 2023-06-27

Family

ID=86413545

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310458339.5A Active CN116164699B (zh) 2023-04-26 2023-04-26 一种行星地震仪的主动源方位角自动识别方法

Country Status (2)

Country Link
US (1) US11899149B1 (zh)
CN (1) CN116164699B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116481487B (zh) * 2023-06-20 2023-09-12 中国科学院地质与地球物理研究所 基于行星日晷的行星地震仪方位角自动识别方法与装置
CN117949394B (zh) * 2024-03-21 2024-06-04 中国科学院地质与地球物理研究所 一种木星星震干涉偏振成像系统及成像方法
CN117930333B (zh) * 2024-03-22 2024-06-14 中国科学院地质与地球物理研究所 基于负反馈控制的数字调零方法及应用其的宽频带行星地震仪

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103645490A (zh) * 2013-12-28 2014-03-19 淮海工学院 智能化卫星连续运行参考站姿态监测方法
JP2015118065A (ja) * 2013-12-20 2015-06-25 ホーチキ株式会社 地震計
CN105388515A (zh) * 2015-10-28 2016-03-09 中国石油天然气股份有限公司 一种检波器的质控方法及装置
CN108594291A (zh) * 2018-06-22 2018-09-28 南方科技大学 一种海底地震仪
CN212364597U (zh) * 2020-06-05 2021-01-15 武汉大学 用于雷达卫星和gnss卫星的高精度定标定位装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115453149A (zh) * 2022-09-14 2022-12-09 天津云圣智能科技有限责任公司 无人机中加速度计的校准方法、装置及电子设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015118065A (ja) * 2013-12-20 2015-06-25 ホーチキ株式会社 地震計
CN103645490A (zh) * 2013-12-28 2014-03-19 淮海工学院 智能化卫星连续运行参考站姿态监测方法
CN105388515A (zh) * 2015-10-28 2016-03-09 中国石油天然气股份有限公司 一种检波器的质控方法及装置
CN108594291A (zh) * 2018-06-22 2018-09-28 南方科技大学 一种海底地震仪
CN212364597U (zh) * 2020-06-05 2021-01-15 武汉大学 用于雷达卫星和gnss卫星的高精度定标定位装置

Also Published As

Publication number Publication date
CN116164699A (zh) 2023-05-26
US11899149B1 (en) 2024-02-13

Similar Documents

Publication Publication Date Title
CN116164699B (zh) 一种行星地震仪的主动源方位角自动识别方法
JP6907448B2 (ja) タンクおよび容器のための音響検量アレイ
US9291733B2 (en) Device and method for determining S-wave attenuation in near-surface condition
Iezzi et al. Three‐dimensional acoustic multipole waveform inversion at Yasur Volcano, Vanuatu
McKee et al. Infrasound signal detection and back azimuth estimation using ground‐coupled airwaves on a seismo‐acoustic sensor pair
Chen et al. Improving GNSS-acoustic positioning by optimizing the ship’s track lines and observation combinations
US10401518B2 (en) Estimation of direct arrival signals based on predicted direct arrival signals and measurements
US11016213B2 (en) Gradient-based 4D seabed acquisition positioning
US20100142326A1 (en) Method For Identifying Gas Leaks Using A Stationary Seabed Placed Steered Beamformed Acoustic Antenna With Active Narrow Beam Transmitter Interrogation Capability
CN103969686A (zh) 用于海洋地震振动源的远场特征波形的确定的装置和方法
US9285495B2 (en) Methods and apparatus for node positioning during seismic survey
Mordret et al. Seismic stereometry reveals preparatory behavior and source kinematics of intermediate‐size earthquakes
AU2009333192B2 (en) System and method for reducing signature variation of seismic sources
US10132946B2 (en) Methods and systems that combine wavefields associated with generalized source activation times and near-continuously recorded seismic data
CA2855301A1 (en) System and method for determining the position of a sensor in seismic exploration
US20200333485A1 (en) Randomizing sweeps in a marine survey
US20130286781A1 (en) Method and Apparatus for Selective Seismic Detection of Elongated Targets
KR101519088B1 (ko) 해양 환경에서의 3차원 탄성파 탐사 방법 및 시스템
RU2457514C1 (ru) Способ определения предвестника цунами
US10228479B2 (en) Dynamic gain adjustments in seismic surveys
US8634270B2 (en) Determining sea conditions in marine seismic spreads
Maher et al. Evidence for near-source nonlinear propagation of volcano infrasound from Strombolian explosions at Yasur Volcano, Vanuatu
US20160131782A1 (en) Parameter variation improvement for seismic data using sensitivity kernels
Chen et al. Six‐Component Earthquake Synchronous Observations Across Taiwan Strait: Phase Velocity and Source Location
RU2484499C1 (ru) Способ определения глубин акватории гидролокатором бокового обзора и гидролокатор бокового обзора для его осуществления

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