CN116738182B - 海底沉积物声衰减范围值的预测方法、装置和存储介质 - Google Patents

海底沉积物声衰减范围值的预测方法、装置和存储介质 Download PDF

Info

Publication number
CN116738182B
CN116738182B CN202310811107.3A CN202310811107A CN116738182B CN 116738182 B CN116738182 B CN 116738182B CN 202310811107 A CN202310811107 A CN 202310811107A CN 116738182 B CN116738182 B CN 116738182B
Authority
CN
China
Prior art keywords
sediment
data
range value
sound attenuation
attenuation coefficient
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
CN202310811107.3A
Other languages
English (en)
Other versions
CN116738182A (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.)
First Institute of Oceanography MNR
Sun Yat Sen University
Original Assignee
First Institute of Oceanography MNR
Sun Yat Sen University
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 First Institute of Oceanography MNR, Sun Yat Sen University filed Critical First Institute of Oceanography MNR
Priority to CN202310811107.3A priority Critical patent/CN116738182B/zh
Publication of CN116738182A publication Critical patent/CN116738182A/zh
Application granted granted Critical
Publication of CN116738182B publication Critical patent/CN116738182B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/02—Investigating particle size or size distribution
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
    • G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/11—Analysing solids by measuring attenuation of acoustic waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/27—Regression, e.g. linear or logistic regression
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Computation (AREA)
  • Pathology (AREA)
  • Artificial Intelligence (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Dispersion Chemistry (AREA)
  • Mathematical Physics (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Software Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种海底沉积物声衰减范围值的预测方法、计算机装置及存储介质，海底沉积物声衰减范围值的预测方法包括海底原位声学测量获取沉积物声衰减系数数据，沉积物物理参数测量，运用支持向量回归算法建立沉积物声衰减系数高频范围值预测模型，以及声衰减范围值预测等步骤。本发明相对于使用经验公式进行计算的相关技术，预测得到的声衰减系数值更加接近于真实值，预测结果误差更低，精度更高，而且声衰减系数上限预测模型和声衰减系数下限预测模型的数据库可扩展，允许后续添加实测数据，进一步对声衰减系数上限预测模型和声衰减系数下限预测模型进行优化升级，从而提高预测精度。本发明广泛应用于海底沉积物声学技术领域。

Description

海底沉积物声衰减范围值的预测方法、装置和存储介质
技术领域
本发明涉及海底沉积物声学技术领域，尤其是一种海底沉积物声衰减范围值的预测方法、计算机装置及存储介质。
背景技术
声衰减系数是海底沉积物最重要的声学性质之一，是声传播理论的基础参数，声衰减系数的准确预报决定着地声模型的精度，是海洋声场研究不可缺少的研究内容。目前海底沉积物的声衰减系数预测主要是根据声衰减系数经验方程进行预测，但是由于经验方程都是根据不同区域的声学数据建立的方程，具有一定的局域性，根据不同海域建立的声衰减系数经验方程并不能通用。而且经验方程大多数是单参数或者双参数方程，而单参数或者双参数并不能完全代表海底沉积物的性质，其预测结果与实际测量结果相差较大。其次，经验方程大多是基于取样测量构建的，但是由于取样测量会对沉积物产生扰动，且改变了其原始沉积环境，取样测量有一定的误差。
发明内容
针对目前的海底沉积物声衰减系数测量技术存在的误差大等技术问题，本发明的目的在于提供一种海底沉积物声衰减范围值的预测方法、计算机装置及存储介质。
一方面，本发明实施例包括一种海底沉积物声衰减范围值的预测方法，包括：
S1、海底原位声学测量获取沉积物声衰减系数数据的步骤：
S101.首先在科考船甲板上对原位测量设备设定发射周期、采样间隔、采样时间和发射
功
率参数；
S102.利用动力定位系统将科考船停泊在工作站位，然后利用A架和船载地质绞车设备，将原位测量设备吊放到海底，确认原位测量设备触底后，原位测量设备开始工作，原位测量设备的声学发射电路激发发射换能器发射声波，测量频率频点为n点；声波信号穿过海底沉积物后，被原位测量设备的接收换能器接收，并存储在存储单元中；
S103.利用A架和船载地质绞车设备将原位测量设备从海底提升到船甲板上，完成一次原位测量；
S104.利用电脑连接原位测量设备的存储单元，将原位系统测量接收的信号导出，后续在电脑中对信号进行分析；
S105.基于公式(1)计算声衰减系数αs：
式中，αs为海底沉积物的声衰减系数，其单位为dB/m，As为采集的沉积物与海水中信号的振幅峰值，Aw为同一通道采集的海水中信号的振幅峰值；采用n次发射和接收，测量获取该站位沉积物声衰减系数的范围值；
S106.在原位测量时，根据GPS卫星信号采集作业站位的经纬度坐标，根据经纬度坐标和Generic Mapping Tools软件计算得到作业站位离最近海岸的距离，根据船载多波束系统获取站位的水深信息，综合得到站位信息；
S2.沉积物物理参数测量的步骤：
S201.在原位测量完成后，在同一站位，利用重力取样器、箱式取样器或抓斗采集相同站位的海底沉积物柱状样品；
S202.对沉积物样品进行沉积物物理参数的测量，具体测量海底沉积物的粒组含量、粒度系数、孔隙度、密度参数，其中粒组含量包括砾石、砂、粉砂和黏土的含量，粒度系数包括平均粒径、中值粒径、偏态值、峰态值、分选系数、曲率系数和不均匀系数；
S203.采用环刀法密度测定，选取代表性试样，利用电热烘箱在105-110℃恒温下进行烘干，烘干时间8小时以上；对含有机质超过干土质量5％的土，将温度控制在65～70℃的恒温下烘至恒重，得到容重；采用比重瓶法测定土粒比重，根据所做的土粒比重试验结果，结合试样的容重来计算孔隙度；
S204.沉积物粒径测量时，采用筛析法和Malvern Mastersizer激光粒度仪测试样品的全样粒度，获得粒组含量、粒度系数数据；
S3.运用支持向量回归算法建立沉积物声衰减系数高频范围值预测模型的步骤：
S301.对于多类沉积环境，根据每一类沉积环境获取的原位测量数据，将物理参数、站位信息、粒组含量和粒度系数整合为一个特征矩阵T，以特征矩阵T作为输入特征数据，将步骤S1中原位测量获取的声衰减系数αs的最大值整理为数据矩阵Lmax，将声衰减系数αs的最小值整理为数据矩阵Lmin，将数据矩阵Lmax和数据矩阵Lmin作为训练标签；
S302.在模型训练时，首先选择沉积环境，将该沉积环境下获取的数据特征矩阵T和数据矩阵Lmax作为训练数据，输入到支持向量回归算法中，分别对各类沉积环境的预测模型进行训练，并根据数据特征矩阵T和数据矩阵Lmax之间的代价函数cost function对各类沉积的预测环境模型进行调参和优化，分别得到各类沉积环境的声衰减预测的max模型；
S303.将数据特征矩阵T和数据矩阵Lmin作为训练数据，输入到支持向量回归算法中，分别对各类沉积环境的预测模型进行训练，并根据数据特征矩阵模型T和数据矩阵Lmin之间的代价函数cost function对各类沉积环境的预测模型进行调参和优化，分别得到各类沉积环境的声衰减预测的min模型；
S304.将max模型和min模型根据公式f(声衰减范围值)＝{f(min模型)，f(max模型)}整合为三类沉积环境的声衰减系数范围值的预测模型；
S4.声衰减范围值预测的步骤：
根据需要预测区域进行沉积环境选择，确定好沉积环境后，将站位信息、物理参数、粒组含量和粒度系数输入到步骤S3构建的声衰减系数范围值预测模型中，得到沉积物的声衰减系数的范围值。
进一步地，原位测量的频率范围为10kHz-120kHz。
进一步地，n为不小于3的整数。
进一步地，步骤S202在陆地实验室中进行。
进一步地，环刀法中所使用的环刀内径为61.8mm，高为20mm。
进一步地，沉积环境包括大陆架、大陆坡和深海盆地。
另一方面，本发明实施例还包括一种计算机装置，包括存储器和处理器，所述存储器用于存储至少一个程序，所述处理器用于加载所述至少一个程序以执行实施例中的一种海底沉积物声衰减范围值的预测方法。
另一方面，本发明实施例还包括一种存储介质，其中存储有处理器可执行的程序，所述处理器可执行的程序在由处理器执行时用于执行实施例中的一种海底沉积物声衰减范围值的预测方法。
本发明的有益效果是：实施例中的海底沉积物声衰减范围值的预测方法，相对于使用经验公式进行计算的相关技术，预测得到的声衰减系数值更加接近于真实值，预测结果误差更低，精度更高，而且声衰减系数上限预测模型和声衰减系数下限预测模型的数据库可扩展，允许后续添加实测数据，进一步对声衰减系数上限预测模型和声衰减系数下限预测模型进行优化升级，从而提高预测精度。
附图说明
图1为实施例中海底沉积物声衰减范围值的预测方法的步骤图；
图2为实施例中海底沉积物声衰减范围值的预测方法的原理示意图。
具体实施方式
本实施例中，参照图1，海底沉积物声衰减范围值的预测方法包括以下步骤：
S1.海底原位声学测量获取沉积物声衰减系数数据；
S2.沉积物物理参数测量；
S3.运用支持向量回归算法建立沉积物声衰减系数高频范围值预测模型；
S4.声衰减范围值预测。
本实施例中，步骤S1-S4的原理如图2所示。
在执行步骤S1，也就是海底原位声学测量获取沉积物声衰减系数数据这一步骤时，具体可以执行以下步骤：
S101.首先在科考船甲板上对原位测量设备设定发射周期、采样间隔、采样时间和发射功率参数；
S102.利用动力定位系统将科考船停泊在工作站位，然后利用A架和船载地质绞车设备，将原位测量设备吊放到海底，确认原位测量设备触底后，原位测量设备开始工作，原位测量设备的声学发射电路激发发射换能器发射声波，测量频率频点为n点；声波信号穿过海底沉积物后，被原位测量设备的接收换能器接收，并存储在存储单元中；
S103.利用A架和船载地质绞车设备将原位测量设备从海底提升到船甲板上，完成一次原位测量；
S104.利用电脑连接原位测量设备的存储单元，将原位系统测量接收的信号导出，后续在电脑中对信号进行分析；
S105.基于公式(1)计算声衰减系数αs：
式中，αs为海底沉积物的声衰减系数，其单位为dB/m，As为采集的沉积物与海水中信号的振幅峰值，Aw为同一通道采集的海水中信号的振幅峰值；采用n次发射和接收，测量获取该站位沉积物声衰减系数的范围值；
S106.在原位测量时，根据GPS卫星信号采集作业站位的经纬度坐标，根据经纬度坐标和Generic Mapping Tools软件计算得到作业站位离最近海岸的距离，根据船载多波束系统获取站位的水深信息，综合得到站位信息。
在执行步骤S2，也就是沉积物物理参数测量这一步骤时，具体可以执行以下步骤：
S201.在原位测量完成后，在同一站位，利用重力取样器、箱式取样器或抓斗采集相同
站
位的海底沉积物柱状样品；
S202.对沉积物样品进行沉积物物理参数的测量，具体测量海底沉积物的粒组含量、粒
度
系数、孔隙度、密度参数，其中粒组含量包括砾石、砂、粉砂和黏土的含量，粒度系数包括平均粒径、中值粒径、偏态值、峰态值、分选系数、曲率系数和不均匀系数；
S203.采用环刀法密度测定，选取代表性试样，利用电热烘箱在105-110℃恒温下进行烘干，烘干时间8小时以上；对含有机质超过干土质量5％的土，将温度控制在65～70℃的恒温下烘至恒重，得到容重；采用比重瓶法测定土粒比重，根据所做的土粒比重试验结果，结合试样的容重来计算孔隙度；
S204.沉积物粒径测量时，采用筛析法和Malvern Mastersizer激光粒度仪测试样品的全样粒度，获得粒组含量、粒度系数数据；
在执行步骤S3，也就是运用支持向量回归算法建立沉积物声衰减系数高频范围值预测模型这一步骤时，具体可以执行以下步骤：
S301.对于多类沉积环境，根据每一类沉积环境获取的原位测量数据，将物理参数、站位信息、粒组含量和粒度系数整合为一个特征矩阵T，以特征矩阵T作为输入特征数据，将步骤S1中原位测量获取的声衰减系数αs的最大值整理为数据矩阵Lmax，将声衰减系数αs的最小值整理为数据矩阵Lmin，将数据矩阵Lmax和数据矩阵Lmin作为训练标签；
S302.在模型训练时，首先选择沉积环境，将该沉积环境下获取的数据特征矩阵T和数据矩阵Lmax作为训练数据，输入到支持向量回归算法中，分别对各类沉积环境的预测模型进行训练，并根据数据特征矩阵T和数据矩阵Lmax之间的代价函数cost function对各类沉积的预测环境模型进行调参和优化，分别得到各类沉积环境的声衰减预测的max模型；
S303.将数据特征矩阵T和数据矩阵Lmin作为训练数据，输入到支持向量回归算法中，分别对各类沉积环境的预测模型进行训练，并根据数据特征矩阵模型T和数据矩阵Lmin之间的代价函数cost function对各类沉积环境的预测模型进行调参和优化，分别得到各类沉积环境的声衰减预测的min模型；
S304.将max模型和min模型根据公式f(声衰减范围值)＝{f(min模型)，f(max模型)}整合为三类沉积环境的声衰减系数范围值的预测模型；
在执行步骤S4，也就是声衰减范围值预测这一步骤时，具体可以执行以下步骤：
根据需要预测区域进行沉积环境选择，确定好沉积环境后，将站位信息、物理参数、粒组含量和粒度系数输入到步骤S3构建的声衰减系数范围值预测模型中，得到沉积物的声衰减系数的范围值。
本实施例中的海底沉积物声衰减范围值的预测方法的原理如图2所示。参照图2所示的原理，本实施例中的海底沉积物声衰减范围值的预测方法具有以下优点：
1、预测精度高
相对于传统的经验公式，如Hamilton经验公式，本实施例中的海底沉积物声衰减范围值的预测方法使用经过训练的声衰减系数上限预测模型和声衰减系数下限预测模型进行预测，预测得到的声衰减系数值更加接近于真实值，预测结果误差更低，精度更高。具体地，在训练得到声衰减系数上限预测模型和声衰减系数下限预测模型时，使用了原位测量(也就是在海底原位环境下直接测量)所得到的数据，因此所获得的声衰减系数上限预测模型和声衰减系数下限预测模型具有更高的精度。
2、预测模型可以进一步升级
在以后的工作中，可以逐渐地将实测数据添加到声衰减系数上限预测模型和声衰减系数下限预测模型的数据库中，从而进一步对声衰减系数上限预测模型和声衰减系数下限预测模型进行优化升级，从而提高预测精度。
3、便于操作
对于操作人员来说，本实施例中的海底沉积物声衰减范围值的预测方法操作简便，只需要输入相关参数就可以得到底质沉积物的声衰减系数，不需要进行其他操作，也没有繁杂的计算公式。
可以通过编写执行本实施例中的一种海底沉积物声衰减范围值的预测方法的计算机程序，将该计算机程序写入至计算机装置或者存储介质中，当计算机程序被读取出来运行时，执行本实施例中的一种海底沉积物声衰减范围值的预测方法，从而实现与实施例中的一种海底沉积物声衰减范围值的预测方法相同的技术效果。
需要说明的是，如无特殊说明，当某一特征被称为“固定”、“连接”在另一个特征，它可以直接固定、连接在另一个特征上，也可以间接地固定、连接在另一个特征上。此外，本公开中所使用的上、下、左、右等描述仅仅是相对于附图中本公开各组成部分的相互位置关系来说的。在本公开中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。此外，除非另有定义，本实施例所使用的所有的技术和科学术语与本技术领域的技术人员通常理解的含义相同。本实施例说明书中所使用的术语只是为了描述具体的实施例，而不是为了限制本发明。本实施例所使用的术语“和/或”包括一个或多个相关的所列项目的任意的组合。
应当理解，尽管在本公开可能采用术语第一、第二、第三等来描述各种元件，但这些元件不应限于这些术语。这些术语仅用来将同一类型的元件彼此区分开。例如，在不脱离本公开范围的情况下，第一元件也可以被称为第二元件，类似地，第二元件也可以被称为第一元件。本实施例所提供的任何以及所有实例或示例性语言(“例如”、“如”等)的使用仅意图更好地说明本发明的实施例，并且除非另外要求，否则不会对本发明的范围施加限制。
应当认识到，本发明的实施例可以由计算机硬件、硬件和软件的组合、或者通过存储在非暂时性计算机可读存储器中的计算机指令来实现或实施。所述方法可以使用标准编程技术-包括配置有计算机程序的非暂时性计算机可读存储介质在计算机程序中实现，其中如此配置的存储介质使得计算机以特定和预定义的方式操作——根据在具体实施例中描述的方法和附图。每个程序可以以高级过程或面向对象的编程语言来实现以与计算机系统通信。然而，若需要，该程序可以以汇编或机器语言实现。在任何情况下，该语言可以是编译或解释的语言。此外，为此目的该程序能够在编程的专用集成电路上运行。
此外，可按任何合适的顺序来执行本实施例描述的过程的操作，除非本实施例另外指示或以其他方式明显地与上下文矛盾。本实施例描述的过程(或变型和/或其组合)可在配置有可执行指令的一个或多个计算机系统的控制下执行，并且可作为共同地在一个或多个处理器上执行的代码(例如，可执行指令、一个或多个计算机程序或一个或多个应用)、由硬件或其组合来实现。所述计算机程序包括可由一个或多个处理器执行的多个指令。
进一步，所述方法可以在可操作地连接至合适的任何类型的计算平台中实现，包括但不限于个人电脑、迷你计算机、主框架、工作站、网络或分布式计算环境、单独的或集成的计算机平台、或者与带电粒子工具或其它成像装置通信等等。本发明的各方面可以以存储在非暂时性存储介质或设备上的机器可读代码来实现，无论是可移动的还是集成至计算平台，如硬盘、光学读取和/或写入存储介质、RAM、ROM等，使得其可由可编程计算机读取，当存储介质或设备由计算机读取时可用于配置和操作计算机以执行在此所描述的过程。此外，机器可读代码，或其部分可以通过有线或无线网络传输。当此类媒体包括结合微处理器或其他数据处理器实现上文所述步骤的指令或程序时，本实施例所述的发明包括这些和其他不同类型的非暂时性计算机可读存储介质。当根据本发明所述的方法和技术编程时，本发明还包括计算机本身。
计算机程序能够应用于输入数据以执行本实施例所述的功能，从而转换输入数据以生成存储至非易失性存储器的输出数据。输出信息还可以应用于一个或多个输出设备如显示器。在本发明优选的实施例中，转换的数据表示物理和有形的对象，包括显示器上产生的物理和有形对象的特定视觉描绘。
以上所述，只是本发明的较佳实施例而已，本发明并不局限于上述实施方式，只要其以相同的手段达到本发明的技术效果，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明保护的范围之内。在本发明的保护范围内其技术方案和/或实施方式可以有各种不同的修改和变化。

Claims (8)

1.一种海底沉积物声衰减范围值的预测方法，其特征在于，所述海底沉积物声衰减范围值的预测方法包括：
S1、海底原位声学测量获取沉积物声衰减系数数据的步骤：
S101.首先在科考船甲板上对原位测量设备设定发射周期、采样间隔、采样时间和发射功率参数；
S102.利用动力定位系统将科考船停泊在工作站位，然后利用A架和船载地质绞车设备，将原位测量设备吊放到海底，确认原位测量设备触底后，原位测量设备开始工作，原位测量设备的声学发射电路激发发射换能器发射声波，测量频率频点为n点；声波信号穿过海底沉积物后，被原位测量设备的接收换能器接收，并存储在存储单元中；
S103.利用A架和船载地质绞车设备将原位测量设备从海底提升到船甲板上，完成一次原位测量；
S104.利用电脑连接原位测量设备的存储单元，将原位系统测量接收的信号导出，后续在电脑中对信号进行分析；
S105.基于公式(1)计算声衰减系数αs：
式中，αs为海底沉积物的声衰减系数，其单位为dB/m，As为采集的沉积物与海水中信号的振幅峰值，Aw为同一通道采集的海水中信号的振幅峰值；采用n次发射和接收，测量获取该站位沉积物声衰减系数的范围值；
S106.在原位测量时，根据GPS卫星信号采集作业站位的经纬度坐标，根据经纬度坐标和Generic Mapping Tools软件计算得到作业站位离最近海岸的距离，根据船载多波束系统获取站位的水深信息，综合得到站位信息；
S2.沉积物物理参数测量的步骤：
S201.在原位测量完成后，在同一站位，利用重力取样器、箱式取样器或抓斗采集相同站位的海底沉积物柱状样品；
S202.对沉积物样品进行沉积物物理参数的测量，具体测量海底沉积物的粒组含量、粒度系数、孔隙度、密度参数，其中粒组含量包括砾石、砂、粉砂和黏土的含量，粒度系数包括平均粒径、中值粒径、偏态值、峰态值、分选系数、曲率系数和不均匀系数；
S203.采用环刀法密度测定，选取代表性试样，利用电热烘箱在105-110℃恒温下进行烘干，烘干时间8小时以上；对含有机质超过干土质量5％的土，将温度控制在65～70℃的恒温下烘至恒重，得到容重；采用比重瓶法测定土粒比重，根据所做的土粒比重试验结果，结合试样的容重来计算孔隙度；
S204.沉积物粒径测量时，采用筛析法和Malvern Mastersizer激光粒度仪测试样品的全样粒度，获得粒组含量、粒度系数数据；
S3.运用支持向量回归算法建立沉积物声衰减系数高频范围值预测模型的步骤：
S301.对于多类沉积环境，根据每一类沉积环境获取的原位测量数据，将物理参数、站位信息、粒组含量和粒度系数整合为一个特征矩阵T，以特征矩阵T作为输入特征数据，将步骤S1中原位测量获取的声衰减系数αs的最大值整理为数据矩阵Lmax，将声衰减系数αs的最小值整理为数据矩阵Lmin，将数据矩阵Lmax和数据矩阵Lmin作为训练标签；
S302.在模型训练时，首先选择沉积环境，将该沉积环境下获取的数据特征矩阵T和数据矩阵Lmax作为训练数据，输入到支持向量回归算法中，分别对各类沉积环境的预测模型进行训练，并根据数据特征矩阵T和数据矩阵Lmax之间的代价函数cost function对各类沉积的预测环境模型进行调参和优化，分别得到各类沉积环境的声衰减预测的max模型；
S303.将数据特征矩阵T和数据矩阵Lmin作为训练数据，输入到支持向量回归算法中，分别对各类沉积环境的预测模型进行训练，并根据数据特征矩阵模型T和数据矩阵Lmin之间的代价函数cost function对各类沉积环境的预测模型进行调参和优化，分别得到各类沉积环境的声衰减预测的min模型；
S304.将max模型和min模型根据公式f(声衰减范围值)＝{f(min模型)，f(max模型)}整合为三类沉积环境的声衰减系数范围值的预测模型；
S4.声衰减范围值预测的步骤：
根据需要预测区域进行沉积环境选择，确定好沉积环境后，将站位信息、物理参数、粒组含量和粒度系数输入到步骤S3构建的声衰减系数范围值预测模型中，得到沉积物的声衰减系数的范围值。
2.根据权利要求1所述的海底沉积物声衰减范围值的预测方法，其特征在于，原位测量的频率范围为10kHz-120kHz。
3.根据权利要求1所述的海底沉积物声衰减范围值的预测方法，其特征在于，n为不小于3的整数。
4.根据权利要求1所述的海底沉积物声衰减范围值的预测方法，其特征在于，步骤S202在陆地实验室中进行。
5.根据权利要求1所述的海底沉积物声衰减范围值的预测方法，其特征在于，环刀法中所使用的环刀内径为61.8mm，高为20mm。
6.根据权利要求1所述的海底沉积物声衰减范围值的预测方法，其特征在于，沉积环境包括大陆架、大陆坡和深海盆地。
7.一种计算机装置，其特征在于，包括存储器和处理器，所述存储器用于存储至少一个程序，所述处理器用于加载所述至少一个程序以执行权利要求1-6任一项所述的一种海底沉积物声衰减范围值的预测方法。
8.一种计算机可读存储介质，其中存储有处理器可执行的程序，其特征在于，所述处理器可执行的程序在由处理器执行时用于执行权利要求1-6任一项所述的一种海底沉积物声衰减范围值的预测方法。
CN202310811107.3A 2023-07-03 2023-07-03 海底沉积物声衰减范围值的预测方法、装置和存储介质 Active CN116738182B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310811107.3A CN116738182B (zh) 2023-07-03 2023-07-03 海底沉积物声衰减范围值的预测方法、装置和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310811107.3A CN116738182B (zh) 2023-07-03 2023-07-03 海底沉积物声衰减范围值的预测方法、装置和存储介质

Publications (2)

Publication Number Publication Date
CN116738182A CN116738182A (zh) 2023-09-12
CN116738182B true CN116738182B (zh) 2024-02-27

Family

ID=87908049

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310811107.3A Active CN116738182B (zh) 2023-07-03 2023-07-03 海底沉积物声衰减范围值的预测方法、装置和存储介质

Country Status (1)

Country Link
CN (1) CN116738182B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117150183B (zh) * 2023-10-27 2024-02-13 西北工业大学青岛研究院 海底沉积物孔隙度计算方法、装置、电子设备及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107247291A (zh) * 2017-06-21 2017-10-13 上海海洋大学 海底浅地层声能衰减模型构建及其两个重要声能衰减特征参数提取方法
CN208520817U (zh) * 2018-07-27 2019-02-19 国家海洋局第一海洋研究所 沉积物声速和声衰减测量装置
CN110210100A (zh) * 2019-05-27 2019-09-06 中国科学院南海海洋研究所 一种高精度海底底质沉积物声速预测方法
WO2021227605A1 (zh) * 2020-05-13 2021-11-18 中国科学院海洋研究所 基于rov的深海沉积物声学参数原位探测系统及方法
WO2022159698A1 (en) * 2021-01-22 2022-07-28 Aramco Services Company Method and system for image-based reservoir property estimation using machine learning

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101535841A (zh) * 2006-02-24 2009-09-16 高尔韦爱尔兰国立大学 单通道地震反射轮廓数据在定量确定海床和分层底基的沉积物物理性质中的应用
CN208887630U (zh) * 2018-06-14 2019-05-21 中国海洋大学 一种深海工程地质环境原位长期观测装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107247291A (zh) * 2017-06-21 2017-10-13 上海海洋大学 海底浅地层声能衰减模型构建及其两个重要声能衰减特征参数提取方法
CN208520817U (zh) * 2018-07-27 2019-02-19 国家海洋局第一海洋研究所 沉积物声速和声衰减测量装置
CN110210100A (zh) * 2019-05-27 2019-09-06 中国科学院南海海洋研究所 一种高精度海底底质沉积物声速预测方法
WO2021227605A1 (zh) * 2020-05-13 2021-11-18 中国科学院海洋研究所 基于rov的深海沉积物声学参数原位探测系统及方法
WO2022159698A1 (en) * 2021-01-22 2022-07-28 Aramco Services Company Method and system for image-based reservoir property estimation using machine learning

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
High-Frequency Dependence of Acoustic Properties of Three Typical Sediments in the South China Sea;Jingqiang Wang 等;Marine Science and Engineering;1-12 *
基于遗传BP神经网络的海底沉积物声速预报;陈文景 等;海洋学报;第38卷(第1期);116-123 *

Also Published As

Publication number Publication date
CN116738182A (zh) 2023-09-12

Similar Documents

Publication Publication Date Title
CN116910473B (zh) 海底沉积物宽频声速预测方法、计算机装置和存储介质
Williams et al. Acoustic backscattering experiments in a well characterized sand sediment: Data/model comparisons using sediment fluid and Biot models
CN110210100B (zh) 一种高精度海底底质沉积物声速预测方法
CN116520431B (zh) 浅表海底沉积物宽频分层声速结构构建方法、装置和介质
KR101435617B1 (ko) 해저 지질 조사용 현장 음파 전달속도 측정 시스템
CN116738182B (zh) 海底沉积物声衰减范围值的预测方法、装置和存储介质
Wang et al. Review of underwater acoustic propagation models.
Kint et al. Uncertainty assessment applied to marine subsurface datasets
CN116644295B (zh) 海底沉积物声衰减宽频预测方法、计算机装置和存储介质
Lucieer et al. Seafloor mapping field manual for multibeam sonar
Muhamad et al. Marine Habitat Mapping using Multibeam Echosounder Survey and Underwater Video Observations: A Case Study from Tioman Marine Park
GB2600431A (en) Method of estimating polymetallic nodule abundance
CN116519799B (zh) 宽频域海底沉积物取样测量声速值校正方法、装置和介质
Ton et al. Beach and nearshore monitoring techniques
Bassett et al. Calibration and Processing of Nortek Signature1000 Echosounders
Osler et al. A sediment probe for the rapid assessment of seabed characteristics
JP2010122081A (ja) 水測能力予測システム、水測能力予測方法、水測能力予測プログラム及び水測機器搭載移動体
Margolina et al. BRS Sound Exposure Modeling Tool: a system for planning, visualization and analysis
Cordero et al. Characterizing free and open-source tools for ocean-mapping
McClatchie et al. Spatial and temporal variability of Antarctic krill: implications for stock assessment
Kwon et al. Spatial coherence analysis of underwater ambient noise measured at the Yellow Sea
Ocalan Seabed sediment volume computations using single beam echo sounder with dual frequencies based on different interpolation methods
Sun et al. Characterisation and prediction of dredge-generated sediment plume dynamics and fate
CN117150183B (zh) 海底沉积物孔隙度计算方法、装置、电子设备及介质
Ardani et al. Optimization of bathymetry estimates for nearshore hydrodynamic models using bayesian methods

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