CN112833863A - 一种海洋孤立内波引起海表面高度变化的测量方法 - Google Patents
一种海洋孤立内波引起海表面高度变化的测量方法 Download PDFInfo
- Publication number
- CN112833863A CN112833863A CN202011644526.5A CN202011644526A CN112833863A CN 112833863 A CN112833863 A CN 112833863A CN 202011644526 A CN202011644526 A CN 202011644526A CN 112833863 A CN112833863 A CN 112833863A
- Authority
- CN
- China
- Prior art keywords
- flow velocity
- relative flow
- ocean current
- observation platform
- sea surface
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
-
- 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/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Hydrology & Water Resources (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开一种海洋孤立内波引起海表面高度变化的测量方法,包括如下步骤:步骤1:获得浮标观测平台所在位置点的不同时间、不同层数的海流相对浮标观测平台的相对流速,对相对流速依次进行质量控制和去除背景海流信息,得到相对流速;步骤2:获得浮标观测平台的垂直加速度与水平加速度,根据水平加速度获得移动速度,将移动速度和相对流速进行叠加,得到真实的海流流速;步骤3:根据真实的海流流速判断是否出现内波,若出现内波,则继续执行,否则结束处理或跳转至步骤1;步骤4:对垂直加速度连续进行两次时间积分,得到海洋孤立内波发生时间段内引起的海表面高度变化。本发明可以较高精度的测量出孤立内波引起的海表面高度变化的信息。
Description
技术领域
本发明涉及海洋监测技术领域,具体涉及一种海洋孤立内波引起海表面高度变化的测量方法。
背景技术
海洋孤立内波(也可以简称为孤立内波)属于海洋内波其中一类,海洋孤立内波是一种发生在海水密度层结的海洋中的波动,其最大振幅出现在海水内部,波长在几百米至几十千米之间,属于海洋小尺度运动。海洋孤立内波引起的瞬时强烈海流在全水深均有发生,特别是海洋孤立内波从开始发生到结束这段时间可能会引起剧烈的海表面高度变化,此种剧烈的海表面高度变化可能会多作业设备以及作业人员带来安全隐患,是海洋工程作业安全的潜在威胁。例如,搭载监测设备或作业人员的小型载体容易在剧烈变化的海面上发送倾覆,从而带来安全隐患。因此,有必要对海洋孤立内波进行预警,以提高作业安全性。
目前,就海洋孤立内波的观测有相关的预警方案,如申请号为201210578674.0、201110147411.X、201210551965.0的中国发明申请专利,以及授权号为ZL201710370753.5的中国发明专利。此外,也还有对海洋孤立内波进行测量相关的方法,如申请号为201210578674.0、201110312123.5的中国发明申请专利。但是,关于测量孤立内波从开始出现到消失的这段时间所引起海表面高度变化的方法比较缺乏,到目前为止,一直还没有一个直接的观测方法,仅仅对海洋内波监测是不够的,还需要一个能够比较直接的对海表面高度变化进行测量的方法。
发明内容
针对现有技术的不足,本发明的目的在于提供一种海洋孤立内波引起海表面高度变化的测量方法,其能够解决测量因海洋孤立内波而引起海洋表面高度变化的问题;
实现本发明的目的的技术方案为:一种海洋孤立内波引起海表面高度变化的测量方法,包括如下步骤:
步骤1:获得浮标观测平台所在位置点的不同时间t、不同层数q的海流相对浮标观测平台的相对流速Vt,q,对相对流速Vt,q依次进行质量控制和去除背景海流信息,经过质量控制后得到相对流速V′t,q,经过去除背景海流后得到相对流速V″t,q;
步骤2:获取所述浮标观测平台不同时间t自身的垂直加速度Gzt和自身水平加速度Ght,对水平加速度Ght进行时间积分后得到浮标观测平台的水平移动速度Vt,将相对流速V′t,q和水平移动速度Vt叠加,得到真实的海流流速Ft,q;
步骤3:判断当前的海流流速Ft,q是否符合内波出现的条件,若符合,则判断出现海洋孤立内波并进入步骤4,否则,判断未出现海洋孤立内波并结束处理或跳转至步骤1;
步骤4:对垂直加速度Gzt连续进行两次时间积分,得到浮标观测平台垂向移动距离,将浮标观测平台垂向移动距离作为海洋孤立内波发生时间段内引起的海表面高度变化。
进一步地,所述质量控制包括以下步骤:
步骤111:取不同时间t、不同层数q的相对流速Vt,q;
步骤112:若相对流速Vt,q满足Percent Good<阈值m或水平流速Ot,l>阈值n,则删除该数据,若不满足,则进行步骤113;
步骤113:对相对流速Vt,q按观测点进行三点滑动平均滤波后,进行步骤114;
步骤114:若水平流速Ot,l相邻上下两层的速度差HVG>阈值R,则删除该数据,否则进行步骤115;
步骤115:若相对流速Vt,q缺值时,则输出V′t,q=(Vt+1,q+Vt-1,q)/2,若相对流速Vt,q不缺值时,则V′t,q=Vt,q。
进一步地,所述去除背景海流信息的具体实现包括以下步骤:
步骤121:当ADCP采集数据超过预设时间Ty后,取经过质量控制后的不同时间、相同层数的相对流速V′t,q;
步骤123:除去背景流速之后的相对流速V″t,q,V″t,q=V′t,q-Vt背景。
进一步地,所述Ty=360分钟。
进一步地,所述将相对流速V′t,q和水平移动速度Vt叠加,得到真实的海流流速Ft,q,其具体实现过程包括如下步骤:
步骤211:取不同时间t的浮标观测平台的水平加速度Ght;
步骤212:对水平加速度Ght进行中值滤波,然后进行步骤213;
步骤213:对经过中值滤波后的Ght对时间t积分,获得浮标观测平台的水平移动速度Vt;
步骤214:对V′t,q和Vt进行叠加,得到真实的海流流速Ft,q。
进一步地,所述步骤4中,在进行第一次时间积分之前,还包括对垂直加速度Gzt进行滤波,然后再对滤波后的垂直加速度Gzt连续进行两次时间积分。
进一步地,所述滤波为中值滤波。
本发明的有益效果为:本发明可以较高精度的测量出孤立内波引起的海表面高度变化的信息,同时有效提高孤立内波监测数据的精度。
附图说明
图1为本发明的流程示意图。
具体实施方式
下面,结合附图以及具体实施方案,对本发明做进一步描述。
如图1所示,一种海洋孤立内波引起海表面高度变化的测量方法,包括如下步骤:
步骤1:获得浮标观测平台所在位置点的不同时间t、不同层数q的海流相对浮标观测平台的相对流速Vt,q,对相对流速Vt,q依次进行质量控制和去除背景海流信息,经过质量控制后得到相对流速V′t,q,经过去除背景海流信息后得到相对流速V″t,q。其中,相对流速Vt,q可通过海流流速仪进行测量获得,海流流速仪可以采用声学多普勒流速剖面仪(ADCP)、声学多普勒流速剖面仪(ADV)、倒置回声测深仪(IES)等设备。本实施例中,在浮标观测平台上安装采样率为300kHz的ADCP和三维重力传感器,浮标观测平台设置在水深2500米位置并且放置与孤立内波的来波方向,采样间隔设置为3分钟,每个采样样本为30个ping,ping时间间隔为1秒,层数为30层,层厚4m,ADCP输出数据包括25层海流数据(即相对流速)和海表温度数据。
采样间隔的具体数值可以根据实际情况进行调整,目前,3分钟的采样间隔是可以观测到海洋内波的最小采样时间间隔。至于ping的数量可以按实际情况进行增减,其他也是类似。
将ADCP安装在浮标上面,ADCP朝下(即朝向海洋)观测,通常可以观测1-30层的海洋内波。
本步骤中,所述质量控制目的在于,通过筛选相关数据得到具有更高数据质量的相对流速V′t,q,实质上相当于是一个去噪过程,其具体实现包括以下步骤:
步骤111:取不同时间t、不同层数q的相对流速Vt,q。
步骤112:若相对流速Vt,q满足Percent Good<阈值m或水平流速Ot,l>阈值n,则删除该数据,若不满足,则进行步骤113。Percent Good是海流流速仪自动输出的数据,是良好波束测量的百分比,用来衡量所测数据质量的好坏,也即是数据良好率。
步骤113:对相对流速Vt,q按观测点进行三点滑动平均滤波后,进行步骤114。三点滑动滤波也即是每三个观测点的相对流速Vt,q进行平均,也即是对相同层数、不同时间时间点的相对流速Vt,q进行三点中值滤波。
步骤114:若水平流速Ot,l相邻上下两层的速度差HVG>阈值R,则删除该数据,否则进行步骤115。
步骤115:若相对流速Vt,q缺值时,则输出V′t,q=(Vt+1,q+Vt-1,q)/2,若相对流速Vt,q不缺值时,则V′t,q=Vt,q。也即,若缺值,则采用下一时刻和上一时刻两者的平均值作为当前时刻的相对流速V′t,q(也即作为质量控制后的输出值),若不缺值,则采用当前时刻的相对流速作为输出值。
所述去除背景海流信息的具体实现包括以下步骤:
步骤121:当ADCP采集数据超过预设时间Ty后(通常为360分钟),取经过质量控制后的不同时间、相同层数的相对流速V′t,q。
步骤123:除去背景流速之后的相对流速V″t,q=V′t,q-Vt背景。
在步骤121-步骤123中,设置具体时间进行相对应的处理,是基于这样的考虑:
ADCP观测获得的海流信息中包括了孤立内波信号和非内波信号(也即是实施例中提到的背景流速等背景信号),非海洋内波信号在海流信息中会影响和混淆海洋内波信号的识别,因此需要提前去除,以便于后面海洋内波的识别。预设时间Ty选取360分钟即6个小时,主要原因为:海洋内波引起的海流流速变化持续时间一般小于60分钟,背景信号的变化周期一般是超过360分钟的,那么获取的背景流信息主要通过时间平均来获得,360分钟时间足够长,如果这段时间有海洋内波信号,360分钟的平均也足够把它平滑掉,同时360分钟的长度平均又可以基本描述背景信号的变化。因此,通过预设时间Ty可以有效将包括背景流进行平滑掉,同时又能基本描述出背景信号变化。
步骤2:获取所述浮标观测平台不同时间t自身的垂直加速度Gzt和自身水平加速度Ght,对水平加速度Ght进行时间积分后得到浮标观测平台的水平移动速度Vt,叠加步骤1的相对流速V′t,q和水平移动速度Vt,得到经过校正后真实的海流流速Ft,q。其中,垂直加速度Gzt也可以称之为垂向加速度。
叠加步骤1的相对流速V′t,q和水平移动速度Vt,得到经过校正后真实的海流流速Ft,q,其具体实现过程包括如下步骤:
步骤211:取不同时间t的浮标观测平台的水平加速度Ght。
步骤212:对水平加速度Ght进行中值滤波,然后进行步骤213。
步骤214:对V′t,q(也即是观测得到的海流流速)和Vt进行叠加,得到真实的海流流速Ft,q,即有Ft,q=V′t,q+Vt。
通过本步骤处理,可以有效消除由于仪器(例如ADCP)偏移引起的观测误差,而获得真实的海流流速,并且为后续基于真实的海流流速进行海洋孤立内波预警强度打下坚实的数据基础。
步骤3:根据海流流速Ft,q判断当前的海流流速Ft,q是否符合内波出现的条件,若符合,则判断出现海洋孤立内波并进入步骤4,否则,判断未出现海洋孤立内波并结束处理或跳转至步骤1重新处理。
其中,根据海流流速Ft,q判断当前的海流流速Ft,q是否符合内波出现的条件,可以采用申请号为2017103707535的中国发明申请专利中公开的技术方案,具体可以采用该发明申请专利中提到的水平流速判别法或垂直流速判别法,由于属于现有的公开技术内容,在此不进行赘述。
步骤4:对垂直加速度Gzt连续进行两次时间积分,得到浮标观测平台垂向移动距离,将浮标观测平台垂向移动距离作为海洋孤立内波发生时间段内引起的海表面高度变化,从而完成海表面高度变化的测量。也即先将垂直加速度Gzt进行第一次时间积分,得到浮标观测平台在垂直方向的移动速度Vzt,然后对Vzt再次进行第二次积分,从而得到海洋孤立内波从开始发生到结束时间段内引起的海表面高度变化Hzt,也即是孤立内波引起的海表面高度变化,
优选地,在本步骤中,在进行第一次时间积分之前,还包括对垂直加速度Gzt进行滤波,然后再对滤波后的垂直加速度Gzt连续进行两次时间积分,其中滤波采用中值滤波。
本发明可以较高精度的测量出孤立内波引起的海表面高度变化的信息,同时有效提高孤立内波监测数据的精度。
本说明书所公开的实施例只是对本发明单方面特征的一个例证,本发明的保护范围不限于此实施例,其他任何功能等效的实施例均落入本发明的保护范围内。对于本领域的技术人员来说,可根据以上描述的技术方案以及构思,做出其它各种相应的改变以及变形,而所有的这些改变以及变形都应该属于本发明权利要求的保护范围之内。
Claims (7)
1.一种海洋孤立内波引起海表面高度变化的测量方法,其特征在于,包括如下步骤:
步骤1:获得浮标观测平台所在位置点的不同时间t、不同层数q的海流相对浮标观测平台的相对流速Vt,q,对相对流速Vt,q依次进行质量控制和去除背景海流信息,经过质量控制后得到相对流速V′t,q,经过去除背景海流后得到相对流速V″t,q;
步骤2:获取所述浮标观测平台不同时间t自身的垂直加速度Gzt和自身水平加速度Ght,对水平加速度Ght进行时间积分后得到浮标观测平台的水平移动速度Vt,将相对流速V′t,q和水平移动速度Vt叠加,得到真实的海流流速Ft,q;
步骤3:判断当前的海流流速Ft,q是否符合内波出现的条件,若符合,则判断出现海洋孤立内波并进入步骤4,否则,判断未出现海洋孤立内波并结束处理或跳转至步骤1;
步骤4:对垂直加速度Gzt连续进行两次时间积分,得到浮标观测平台垂向移动距离,将浮标观测平台垂向移动距离作为海洋孤立内波发生时间段内引起的海表面高度变化。
2.根据权利要求1所述的海洋孤立内波引起海表面高度变化的测量方法,其特征在于,所述质量控制包括以下步骤:
步骤111:取不同时间t、不同层数q的相对流速Vt,q;
步骤112:若相对流速Vt,q满足Percent Good<阈值m或水平流速Ot,l>阈值n,则删除该数据,若不满足,则进行步骤113;
步骤113:对相对流速Vt,q按观测点进行三点滑动平均滤波后,进行步骤114;
步骤114:若水平流速Ot,l相邻上下两层的速度差HVG>阈值R,则删除该数据,否则进行步骤115;
步骤115:若相对流速Vt,q缺值时,则输出V′t,q=(Vt+1,q+Vt-1,q)/2,若相对流速Vt,q不缺值时,则V′t,q=Vt,q。
4.根据权利要求1所述的海洋孤立内波引起海表面高度变化的测量方法,其特征在于,所述Ty=360分钟。
5.根据权利要求1所述的海洋孤立内波引起海表面高度变化的测量方法,其特征在于,所述将相对流速V′t,q和水平移动速度Vt叠加,得到真实的海流流速Ft,q,其具体实现过程包括如下步骤:
步骤211:取不同时间t的浮标观测平台的水平加速度Ght;
步骤212:对水平加速度Ght进行中值滤波,然后进行步骤213;
步骤213:对经过中值滤波后的Ght对时间t积分,获得浮标观测平台的水平移动速度Vt;
步骤214:对V′t,q和Vt进行叠加,得到真实的海流流速Ft,q。
6.根据权利要求1所述的海洋孤立内波引起海表面高度变化的测量方法,其特征在于,所述步骤4中,在进行第一次时间积分之前,还包括对垂直加速度Gzt进行滤波,然后再对滤波后的垂直加速度Gzt连续进行两次时间积分。
7.根据权利要求6所述的海洋孤立内波引起海表面高度变化的测量方法,其特征在于,所述滤波为中值滤波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011644526.5A CN112833863B (zh) | 2020-12-31 | 2020-12-31 | 一种海洋孤立内波引起海表面高度变化的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011644526.5A CN112833863B (zh) | 2020-12-31 | 2020-12-31 | 一种海洋孤立内波引起海表面高度变化的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112833863A true CN112833863A (zh) | 2021-05-25 |
CN112833863B CN112833863B (zh) | 2022-08-16 |
Family
ID=75927166
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011644526.5A Active CN112833863B (zh) | 2020-12-31 | 2020-12-31 | 一种海洋孤立内波引起海表面高度变化的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112833863B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114646304A (zh) * | 2022-03-11 | 2022-06-21 | 广州海洋地质调查局 | 一种基于多波束数据的海洋内波识别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110060525A1 (en) * | 2009-09-06 | 2011-03-10 | Chung-Chu Teng | System for monitoring, determining, and reporting directional spectra of ocean surface waves in near real-time from a moored buoy |
CN107218931A (zh) * | 2017-05-23 | 2017-09-29 | 广东贝达海洋科学有限公司 | 一种基于在线监测数据的海洋内波预警***及方法 |
CN209938881U (zh) * | 2019-04-11 | 2020-01-14 | 南京信息工程大学 | 基于六维加速度传感器的海洋波浪浮标 |
CN110986897A (zh) * | 2019-12-27 | 2020-04-10 | 武汉迈普时空导航科技有限公司 | 一种基于组合导航技术的涌浪测量方法及*** |
CN111611720A (zh) * | 2020-05-28 | 2020-09-01 | 自然资源部第一海洋研究所 | 一种适用于真实海洋的浪场非高斯态程度评估方法 |
-
2020
- 2020-12-31 CN CN202011644526.5A patent/CN112833863B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110060525A1 (en) * | 2009-09-06 | 2011-03-10 | Chung-Chu Teng | System for monitoring, determining, and reporting directional spectra of ocean surface waves in near real-time from a moored buoy |
CN107218931A (zh) * | 2017-05-23 | 2017-09-29 | 广东贝达海洋科学有限公司 | 一种基于在线监测数据的海洋内波预警***及方法 |
CN209938881U (zh) * | 2019-04-11 | 2020-01-14 | 南京信息工程大学 | 基于六维加速度传感器的海洋波浪浮标 |
CN110986897A (zh) * | 2019-12-27 | 2020-04-10 | 武汉迈普时空导航科技有限公司 | 一种基于组合导航技术的涌浪测量方法及*** |
CN111611720A (zh) * | 2020-05-28 | 2020-09-01 | 自然资源部第一海洋研究所 | 一种适用于真实海洋的浪场非高斯态程度评估方法 |
Non-Patent Citations (3)
Title |
---|
MAHMOUD M. SADEK,ET AL.: "Two-dimensional instability of the bottom boundary layer under a solitary wave", 《PHYSICS OF FLUIDS》 * |
刘国栋: "波浪浮标数据处理方法研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑(2012年)》 * |
陈敏: "内孤立波与深海浮式平台相互作用特性研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑(2020年)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114646304A (zh) * | 2022-03-11 | 2022-06-21 | 广州海洋地质调查局 | 一种基于多波束数据的海洋内波识别方法 |
CN114646304B (zh) * | 2022-03-11 | 2022-11-08 | 广州海洋地质调查局 | 一种基于多波束数据的海洋内波识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112833863B (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107218931B (zh) | 一种基于在线监测数据的海洋内波预警***及方法 | |
WO2006104087A1 (ja) | 所定水温域の深度分布予測方法、回遊性魚類の漁場予測方法及び回遊性魚類の漁場予測情報配信システム | |
CN111024207B (zh) | 一种矢量水听器线谱自动检测与判决方法 | |
CN108562891B (zh) | 一种深海低信噪比条件下声源深度自主实时跟踪方法 | |
CN112833863B (zh) | 一种海洋孤立内波引起海表面高度变化的测量方法 | |
CN107064929B (zh) | 一种利用s波段多普勒雷达探测海面浪高的方法 | |
KR20110031353A (ko) | 타깃 파라미터의 수동적 결정 방법 및 장치 | |
CN112683245B (zh) | 一种海洋孤立内波预警强度校正方法 | |
JP2009075017A (ja) | 表層流速推定方法、装置並びにプログラム | |
Wadhams et al. | Sidescan sonar imagery of the winter marginal ice zone obtained from an AUV | |
JP6567665B2 (ja) | 船舶の航路(パス)の全地点における各々のドリフト(浮動)ベクトルを推定する方法 | |
Andres et al. | Icebergs and sea ice detected with inverted echo sounders | |
CN116642468A (zh) | 基于无人机航空摄影和无人船水上水下一体化扫描方法 | |
JP2009229432A (ja) | Gps受信体を利用した津波観測システム | |
WO2022034748A1 (ja) | 水中騒音監視装置、水中騒音監視方法及び記録媒体 | |
KR101293840B1 (ko) | 파랑정보의 실시간 공유를 통한 선박 안전운항방법 | |
JP3726112B2 (ja) | Gps式波向計測方法 | |
Peterson et al. | Sea-ice thickness measurement: recent experiments using helicopter-borne EM-induction sensors | |
JP6983960B2 (ja) | 信号処理装置及び信号処理方法 | |
KR101641129B1 (ko) | x 밴드 레이더를 이용한 해류속도 및 연안 수심 측정 시스템 및 그 시스템을 이용한 해류속도 및 연안 수심 측정방법 | |
CN110618452A (zh) | 一种基于小波技术的自适应富钴结壳厚度提取方法 | |
JP2007003396A (ja) | レーダ信号処理装置 | |
Gubnitsky et al. | Inter-pulse estimation for sperm whale click detection | |
CN114721042B (zh) | 一种去除浅剖多次波的方法 | |
CN116697945B (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 |