CN116519799B - 宽频域海底沉积物取样测量声速值校正方法、装置和介质 - Google Patents
宽频域海底沉积物取样测量声速值校正方法、装置和介质 Download PDFInfo
- Publication number
- CN116519799B CN116519799B CN202310797298.2A CN202310797298A CN116519799B CN 116519799 B CN116519799 B CN 116519799B CN 202310797298 A CN202310797298 A CN 202310797298A CN 116519799 B CN116519799 B CN 116519799B
- Authority
- CN
- China
- Prior art keywords
- sound velocity
- sediment
- sampling
- measurement
- situ
- 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
- 239000013049 sediment Substances 0.000 title claims abstract description 101
- 238000005070 sampling Methods 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000012937 correction Methods 0.000 claims abstract description 70
- 238000005259 measurement Methods 0.000 claims abstract description 55
- 238000012625 in-situ measurement Methods 0.000 claims abstract description 44
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 239000013535 sea water Substances 0.000 claims description 48
- 238000011065 in-situ storage Methods 0.000 claims description 27
- 238000012360 testing method Methods 0.000 claims description 18
- 230000005484 gravity Effects 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 13
- 238000011835 investigation Methods 0.000 claims description 7
- 239000004927 clay Substances 0.000 claims description 6
- 238000011545 laboratory measurement Methods 0.000 claims description 4
- 238000000611 regression analysis Methods 0.000 claims description 3
- 150000003839 salts Chemical class 0.000 claims description 3
- 239000000758 substrate Substances 0.000 abstract description 3
- 238000004590 computer program Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- 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/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H5/00—Measuring propagation velocity of ultrasonic, sonic or infrasonic waves, e.g. of pressure 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
- G01N29/4463—Signal correction, e.g. distance amplitude correction [DAC], distance gain size [DGS], noise filtering
-
- 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
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种宽频域海底沉积物取样测量声速值校正方法、装置和介质,本发明根据大量数据进行统计回归,构建出变声速比情况下的声速比校正函数,并充分考虑频率和底质类型等因素,校正结果具有更高的精度,可以针对不同类型的海底沉积物,分别构建宽频域的声速校正函数,校正函数的适用对象范围广,操作简便,只需要输入相关参数就可以将取样测量结果校正至原位测量结果,不需要进行其他操作,也没有繁杂的计算公式。本发明广泛应用于海底沉积物声学技术领域。
Description
技术领域
本发明涉及海底沉积物声学技术领域,尤其是一种宽频域海底沉积物取样测量声速值校正方法、计算机装置及存储介质。
背景技术
海底沉积物的声学特性参数一直是水声通讯、潜艇、水雷探测及海洋工程勘查等领域研究的重要内容,目前的主要研究方式为海底原位测量以及海底取样后实验室测量。由于受到温度压力环境因素的改变,取样测量结果通常被认为失真,需要通过某种方法校正至原位状态下的真实值,如何校正到海底原位测量值一直是沉积声学研究的重点和难点。
发明内容
针对目前需要将海底沉积物的声学测量结果校正至海底原位测量值的技术问题,本发明的目的在于提供一种宽频域海底沉积物取样测量声速值校正方法、计算机装置及存储介质。
一方面,本发明实施例包括一种宽频域海底沉积物取样测量声速值校正方法,包括:
S1.海底沉积物原位声学测量的步骤
S101.首先在科考船甲板上对原位测量设备设定发射周期、采样间隔、采样时间和发射功率参数;
S102.利用动力定位***将科考船停泊在工作站位,然后利用A架和船载地质绞车设备,将原位测量设备吊放到海底,确认原位测量设备触底后,原位测量设备开始工作,原位测量设备的声学发射电路激发发射换能器发射声波,发射声波中心测量频率分别为8kHz、10kHz、16 kHz、25kHz、31.5 kHz、50 kHz、80 kHz和100kHz;声波信号穿过海底沉积物后,被原位测量设备的接收换能器接收,并存储在存储单元中;
S103.利用A架和船载地质绞车设备将原位测量设备从海底提升到船甲板上,完成一次原位测量;
S104.利用电脑连接原位测量设备的存储单元,将原位***测量接收的信号导出,后续在电脑中对信号进行分析;
S105.基于公式(1)声波传播声压和距离计算沉积物的声速:
(1)
式中,V 原位为沉积物原位测量声速值,L1和L2为原位测量***发射换能器和两个接收换能器的距离,∆L是声波在海底沉积物中传播的距离差,t1和t2为声波穿过海底沉积物到达两个接收换能器的时间,∆t为声波到达两个接收换能器的时间差;
S106.原位测量***在海底测试时,利用温盐深仪CTD或者声速剖面仪SVP设备同步采集近海底海水的温度、声速或深度信息;
S2.海底沉积物声学取样测量的步骤
S201.在相同站位进行海底沉积物取样工作,采样采用1套重力式取样器进行作业,当重力式取样器到达海底时,靠重力使采样筒***海底沉积物中,然后将取样器提升到船甲板上,从而获取沉积物柱状样品。
S202.利用声波信号发射采集***在实验室中对获取的沉积物柱状样品进行声学特性测试,测量中心测量频率分别为25kHz,31.5 kHz,50 kHz,80 kHz,及100kHz。声波信号发射采集***主要包括声波发射和采集电路、功率放大器、发射和接收换能器、柱状样品声学测试平台。在实验室采用轴向测试方法,首先将柱状样品进行分割成约30 cm长度的样品段L,放置在声学测试平台上进行测试,发射换能器T和接收换能器R分别放置在柱状样品两端,然后发射声波,完成对L段沉积物的声学测试,沉积物声速计算公式为:
(2)
式中:V 取样为沉积物取样测量声速值,单位为m/s;t为声信号初至时间,单位为μs;t0为换能器零声时修正值,单位为μs;L为沉积物的长度,单位为mm;实验室测量环境为标准大气压,温度为23℃;
S3.沉积物声速比计算的步骤
S301.根据Mackenzie公式计算海水的声速值:
V(T, S, z) = a1 + a2T + a3T2 + a4T3 + a5(S - 35) + a6z + a7z2 + a8T(S- 35) + a9Tz3
这里T表示温度,单位为摄氏度;S表示盐度,单位为千分比;z表示深度,单位为米;
S302.实验室取样测量的温度是T=23℃,盐度采用和海底海水一样的盐度S=S海底,深度z=0,代入到Mackenzie公式中,得到实验室海水的声速值V实验室海水;
S303.将步骤S1采集到的海底海水的温度T海底,盐度S海底,和深度z数据代入到Mackenzie公式中,可以得到海底海水的声速值V海底海水;
S304.分别计算实验室中的速度比和海底原位的速度比:
(3)
(4)
S4.设定校正方程和校正因子的步骤
S401.执行步骤S1获得多种不同的沉积物类型各自的n个原位测量声速值,执行步骤S2获得多种不同的沉积物类型各自的n个海底海水速度值,执行步骤S3获得多种不同的沉积物类型各自的n个取样测量声速值和n个实验室海水速度值,执行步骤S3计算得到n个V实验室声速比和V海底声速比数据;
S402.选取步骤S401中的多个沉积物类型中的一个沉积物类型,将计算得到的这一类沉积物的n个的V实验室声速比和V海底声速比数据进行回归分析,构建以下校正方程:
(5)
其中为校正因子,根据一类沉积物的n个的V实验室声速比和V海底声速比数据可以计算得到这类沉积物的校正因子/>数值;
S403.重复S401和S402从而得到多个沉积物类型各自的校正方程和校正因子;
S5.获取宽频域的校正因子曲线G(f)的步骤
S501.将不同测量频率25kHz、31.5 kHz、50 kHz,80 kHz及100kHz的原位测量值和取样测量值,分别重复执行步骤S1、S2和S3,获得不同频率的校正因子值;
S502.将不同频率的校正因子值与频率f=25kHz,31.5 kHz、50 kHz、80 kHz及100kHz进行拟合,得到一条/>和频率的回归曲线函数/>=G(f),回归曲线函数的频率范围是25kHz-100kHz;
S6.将取样测量声速值校正到原位测量声速值的步骤
S601.首先利用取样样品在实验室中测量,获取得到多种沉积物类型中一类沉积物的取样测量声速值及实验室海水声速值,并计算得到实验室速度比,在取样测量时,同步测量海底海水的温度、盐度和深度信息;
S602.利用Mackenzie公式计算得到海底海水的声速值V海底海水;
S603.根据实验室取样测量频率f,结合回归曲线=G(f),确定校正因子/>的值,然后将海底海水速度值、取样测量声速值和实验室海水速度值带入到校正方程中,计算得到对应频率的海底原位声速值:
(6)
(7)
S604.将宽频域取样测量声速值校正到原位测量声速值,宽频域的频率范围是25kHz-100kHz。
进一步地,重力式取样器中的采样管的长度为6m。
进一步地,重力式取样器中的压载铅块的重量为0.6t以上。
进一步地,重力式取样器中的PVC衬管的内径为75-110mm。
进一步地,常数a1, a2, ..., a9的取值为:
a1 = 1448.96, a2 = 4.591, a3 = -5.304×10-2, a4 = 2.374×10-4, a5 =1.340, a6 = 1.630×10-2, a7 = 1.675×10-7, a8 = -1.025×10-2, a9 = -7.139×10-13。
进一步地,n是不小于100的正整数。
进一步地,海底沉积物类型包括黏土质粉砂、粉砂质黏土、砂质粉砂和粉砂质砂。
进一步地,V实验室声速比和V海底声速比的小数点位数保留6位以上。
另一方面,本发明实施例还包括一种计算机装置,包括存储器和处理器,所述存储器用于存储至少一个程序,所述处理器用于加载所述至少一个程序以执行实施例中的一种宽频域海底沉积物取样测量声速值校正方法。
另一方面,本发明实施例还包括一种存储介质,其中存储有处理器可执行的程序,所述处理器可执行的程序在由处理器执行时用于执行实施例中的一种宽频域海底沉积物取样测量声速值校正方法。
本发明的有益效果是:实施例中的宽频域海底沉积物取样测量声速值校正方法,提出通过大量数据统计回归,构建出变声速比情况下的声速比校正函数,并充分考虑频率和底质类型等因素,校正结果具有更高的精度,可以针对不同类型的海底沉积物,分别构建宽频域的声速校正函数,校正函数的适用对象范围广,操作简便,只需要输入相关参数就可以将取样测量结果校正至原位测量结果,不需要进行其他操作,也没有繁杂的计算公式。
附图说明
图1为实施例中宽频域海底沉积物取样测量声速值校正方法的步骤图;
图2为实施例中宽频域海底沉积物取样测量声速值校正方法的原理示意图。
具体实施方式
本实施例中,参照图1,宽频域海底沉积物取样测量声速值校正方法包括以下步骤:
S1.海底沉积物原位声学测量;
S2.海底沉积物声学取样测量;
S3.沉积物声速比计算;
S4.设定校正方程和校正因子;
S5.获取宽频域的校正因子曲线G(f);
S6.将取样测量声速值校正到原位测量声速值。
本实施例中,步骤S1-S6的原理如图2所示。
在执行步骤S1,也就是海底沉积物原位声学测量这一步骤时,具体可以执行以下步骤:
S101.首先在科考船甲板上对原位测量设备设定发射周期、采样间隔、采样时间和发射功率参数;
S102.利用动力定位***将科考船停泊在工作站位,然后利用A架和船载地质绞车设备,将原位测量设备吊放到海底,确认原位测量设备触底后,原位测量设备开始工作,原位测量设备的声学发射电路激发发射换能器发射声波,发射声波中心测量频率分别为8kHz、10kHz、16 kHz、25kHz、31.5 kHz、50 kHz、80 kHz和100kHz;声波信号穿过海底沉积物后,被原位测量设备的接收换能器接收,并存储在存储单元中;
S103.利用A架和船载地质绞车设备将原位测量设备从海底提升到船甲板上,完成一次原位测量;
S104.利用电脑连接原位测量设备的存储单元,将原位***测量接收的信号导出,后续在电脑中对信号进行分析;
S105.基于公式(1)声波传播声压和距离计算沉积物的声速:
(1)
式中,V 原位为沉积物原位测量声速值,L1和L2为原位测量***发射换能器和两个接收换能器的距离,∆L是声波在海底沉积物中传播的距离差,t1和t2为声波穿过海底沉积物到达两个接收换能器的时间,∆t为声波到达两个接收换能器的时间差;
S106.原位测量***在海底测试时,利用温盐深仪CTD或者声速剖面仪SVP设备同步采集近海底海水的温度、声速或深度信息。
在执行步骤S2,也就是海底沉积物声学取样测量这一步骤时,具体可以执行以下步骤:
S201.在相同站位进行海底沉积物取样工作,采样采用1套重力式取样器进行作业。其中,
重力式取样器中的采样管的长度为6m,压载铅块的重量为0.6t以上,PVC衬管的内径为75-110mm。当重力式取样器到达海底时,靠重力使采样筒***海底沉积物中,然后将取样器提升到船甲板上,从而获取沉积物柱状样品。
S202.利用声波信号发射采集***在实验室中对获取的沉积物柱状样品进行声学特性测试,测量中心测量频率分别为25kHz,31.5 kHz,50 kHz,80 kHz,及100kHz。声波信号发射采集***主要包括声波发射和采集电路、功率放大器、发射和接收换能器、柱状样品声学测试平台。在实验室采用轴向测试方法,首先将柱状样品进行分割成约30 cm长度的样品段L,放置在声学测试平台上进行测试,发射换能器T和接收换能器R分别放置在柱状样品两端,然后发射声波,完成对L段沉积物的声学测试,沉积物声速计算公式为:
(2)
式中:V 取样为沉积物取样测量声速值,单位为m/s;t为声信号初至时间,单位为μs;t0为换能器零声时修正值,单位为μs;L为沉积物的长度,单位为mm;实验室测量环境为标准大气压,温度为23℃。
在执行步骤S3,也就是沉积物声速比计算这一步骤时,具体可以执行以下步骤:
S301.根据Mackenzie公式计算海水的声速值:
V(T, S, z) = a1 + a2T + a3T2 + a4T3 + a5(S - 35) + a6z + a7z2 + a8T(S- 35) + a9Tz3
这里T表示温度,单位为摄氏度;S表示盐度,单位为千分比;z表示深度,单位为米。
本实施例中,常数a1, a2, ..., a9的取值为:
a1 = 1448.96, a2 = 4.591, a3 = -5.304×10-2, a4 = 2.374×10-4, a5 =1.340, a6 = 1.630×10-2, a7 = 1.675×10-7, a8 = -1.025×10-2, a9 = -7.139×10-13。
S302.实验室取样测量的温度是T=23℃,盐度采用和海底海水一样的盐度S=S海底,深度z=0,代入到Mackenzie公式中,得到实验室海水的声速值V实验室海水;
S303.将步骤S1采集到的海底海水的温度T海底,盐度S海底,和深度z数据代入到Mackenzie公式中,可以得到海底海水的声速值V海底海水;
S304.分别计算实验室中的速度比和海底原位的速度比:
(3)
(4)
在执行步骤S4,也就是设定校正方程和校正因子这一步骤时,具体可以执行以下步骤:
S401.执行步骤S1获得多种不同的沉积物类型各自的n个原位测量声速值,执行步骤S2获得多种不同的沉积物类型各自的n个海底海水速度值,执行步骤S3获得多种不同的沉积物类型各自的n个取样测量声速值和n个实验室海水速度值,执行步骤S3计算得到n个V实验室声速比和V海底声速比数据;本实施例中,n是不小于100的正整数;多个海底沉积物类型包括黏土质粉砂、粉砂质黏土、砂质粉砂和粉砂质砂等中国近海常见的海底沉积物类型;
S402.选取步骤S401中的多个沉积物类型中的一个沉积物类型,将计算得到的这一类沉积物的n个的V实验室声速比和V海底声速比数据进行回归分析,构建以下校正方程:
(5)
其中为校正因子,根据一类沉积物的n个的V实验室声速比和V海底声速比数据可以计算得到这类沉积物的校正因子/>数值;本实施例中,V实验室声速比和V海底声速比的小数点位数保留6位以上,以获得较高精度;
S403.重复S401和S402从而得到多个沉积物类型各自的校正方程和校正因子。
在执行步骤S5,也就是获取宽频域的校正因子曲线G(f)这一步骤时,具体可以执行以下步骤:
S501.将不同测量频率25kHz、31.5 kHz、50 kHz,80 kHz及100kHz的原位测量值和取样测量值,分别重复执行步骤S1、S2和S3,获得不同频率的校正因子值;
S502.将不同频率的校正因子值与频率f=25kHz,31.5 kHz、50 kHz、80 kHz及100kHz进行拟合,得到一条/>和频率的回归曲线函数/>=G(f),回归曲线函数的频率范围是25kHz-100kHz。
在执行步骤S6,也就是将取样测量声速值校正到原位测量声速值这一步骤时,具体可以执行以下步骤:
S601.首先利用取样样品在实验室中测量,获取得到多种沉积物类型中一类沉积物的取样测量声速值及实验室海水声速值,并计算得到实验室速度比,在取样测量时,同步测量海底海水的温度、盐度和深度信息;
S602.利用Mackenzie公式计算得到海底海水的声速值V海底海水;
S603.根据实验室取样测量频率f,结合回归曲线=G(f),确定校正因子/>的值,然后将海底海水速度值、取样测量声速值和实验室海水速度值带入到校正方程中,计算得到对应频率的海底原位声速值:
(6)
(7)
S604.将宽频域取样测量声速值校正到原位测量声速值,宽频域的频率范围是25kHz-100kHz。
本实施例中的宽频域海底沉积物取样测量声速值校正方法具有以下优点:
1、预测精度高
传统提出的预测方法认为声速比是保持不变的,与实际情况不符。实施例中的宽频域海底沉积物取样测量声速值校正方法提出通过大量数据统计回归,构建出变声速比情况下的声速比校正函数,并充分考虑频率和底质类型等因素,校正结果具有更高的精度。
2、校正方法的可扩展性
实施例中的宽频域海底沉积物取样测量声速值校正方法可以针对不同类型的海底沉积物,分别构建宽频域的声速校正函数,校正函数的适用对象范围广。
3、便于操作
对于操作人员来说,实施例中的宽频域海底沉积物取样测量声速值校正方法操作简便,只需要输入相关参数就可以将取样测量结果校正至原位测量结果,不需要进行其他操作,也没有繁杂的计算公式。
可以通过编写执行本实施例中的一种宽频域海底沉积物取样测量声速值校正方法的计算机程序,将该计算机程序写入至计算机装置或者存储介质中,当计算机程序被读取出来运行时,执行本实施例中的一种宽频域海底沉积物取样测量声速值校正方法,从而实现与实施例中的一种宽频域海底沉积物取样测量声速值校正方法相同的技术效果。
需要说明的是,如无特殊说明,当某一特征被称为“固定”、“连接”在另一个特征,它可以直接固定、连接在另一个特征上,也可以间接地固定、连接在另一个特征上。此外,本公开中所使用的上、下、左、右等描述仅仅是相对于附图中本公开各组成部分的相互位置关系来说的。在本公开中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。此外,除非另有定义,本实施例所使用的所有的技术和科学术语与本技术领域的技术人员通常理解的含义相同。本实施例说明书中所使用的术语只是为了描述具体的实施例,而不是为了限制本发明。本实施例所使用的术语“和/或”包括一个或多个相关的所列项目的任意的组合。
应当理解,尽管在本公开可能采用术语第一、第二、第三等来描述各种元件,但这些元件不应限于这些术语。这些术语仅用来将同一类型的元件彼此区分开。例如,在不脱离本公开范围的情况下,第一元件也可以被称为第二元件,类似地,第二元件也可以被称为第一元件。本实施例所提供的任何以及所有实例或示例性语言(“例如”、“如”等)的使用仅意图更好地说明本发明的实施例,并且除非另外要求,否则不会对本发明的范围施加限制。
应当认识到,本发明的实施例可以由计算机硬件、硬件和软件的组合、或者通过存储在非暂时性计算机可读存储器中的计算机指令来实现或实施。所述方法可以使用标准编程技术-包括配置有计算机程序的非暂时性计算机可读存储介质在计算机程序中实现,其中如此配置的存储介质使得计算机以特定和预定义的方式操作——根据在具体实施例中描述的方法和附图。每个程序可以以高级过程或面向对象的编程语言来实现以与计算机***通信。然而,若需要,该程序可以以汇编或机器语言实现。在任何情况下,该语言可以是编译或解释的语言。此外,为此目的该程序能够在编程的专用集成电路上运行。
此外,可按任何合适的顺序来执行本实施例描述的过程的操作,除非本实施例另外指示或以其他方式明显地与上下文矛盾。本实施例描述的过程(或变型和/或其组合)可在配置有可执行指令的一个或多个计算机***的控制下执行,并且可作为共同地在一个或多个处理器上执行的代码(例如,可执行指令、一个或多个计算机程序或一个或多个应用)、由硬件或其组合来实现。所述计算机程序包括可由一个或多个处理器执行的多个指令。
进一步,所述方法可以在可操作地连接至合适的任何类型的计算平台中实现,包括但不限于个人电脑、迷你计算机、主框架、工作站、网络或分布式计算环境、单独的或集成的计算机平台、或者与带电粒子工具或其它成像装置通信等等。本发明的各方面可以以存储在非暂时性存储介质或设备上的机器可读代码来实现,无论是可移动的还是集成至计算平台,如硬盘、光学读取和/或写入存储介质、RAM、ROM等,使得其可由可编程计算机读取,当存储介质或设备由计算机读取时可用于配置和操作计算机以执行在此所描述的过程。此外,机器可读代码,或其部分可以通过有线或无线网络传输。当此类媒体包括结合微处理器或其他数据处理器实现上文所述步骤的指令或程序时,本实施例所述的发明包括这些和其他不同类型的非暂时性计算机可读存储介质。当根据本发明所述的方法和技术编程时,本发明还包括计算机本身。
计算机程序能够应用于输入数据以执行本实施例所述的功能,从而转换输入数据以生成存储至非易失性存储器的输出数据。输出信息还可以应用于一个或多个输出设备如显示器。在本发明优选的实施例中,转换的数据表示物理和有形的对象,包括显示器上产生的物理和有形对象的特定视觉描绘。
以上所述,只是本发明的较佳实施例而已,本发明并不局限于上述实施方式,只要其以相同的手段达到本发明的技术效果,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。在本发明的保护范围内其技术方案和/或实施方式可以有各种不同的修改和变化。
Claims (4)
1.一种宽频域海底沉积物取样测量声速值校正方法,其特征在于,所述海底沉积物类型包括黏土质粉砂、粉砂质黏土、砂质粉砂和粉砂质砂,所述宽频域海底沉积物取样测量声速值校正方法包括:
S1.海底沉积物原位声学测量的步骤;
S2.海底沉积物声学取样测量的步骤;
S3.沉积物声速比计算的步骤;
S4.设定校正方程和校正因子的步骤;
S5.获取宽频域的校正因子曲线G(f)的步骤;
S6.将取样测量声速值校正到原位测量声速值的步骤;
所述海底沉积物原位声学测量的步骤,包括:
S101.首先在科考船甲板上对原位测量设备设定发射周期、采样间隔、采样时间和发射功率参数;
S102.利用动力定位***将科考船停泊在工作站位,然后利用A架和船载地质绞车设备,将原位测量设备吊放到海底,确认原位测量设备触底后,原位测量设备开始工作,原位测量设备的声学发射电路激发发射换能器发射声波,发射声波中心测量频率分别为8kHz、10kHz、16 kHz、25kHz、31.5 kHz、50 kHz、80 kHz和100kHz;声波信号穿过海底沉积物后,被原位测量设备的接收换能器接收,并存储在存储单元中;
S103.利用A架和船载地质绞车设备将原位测量设备从海底提升到船甲板上,完成一次原位测量;
S104.利用电脑连接原位测量设备的存储单元,将原位***测量接收的信号导出,后续在电脑中对信号进行分析;
S105.基于公式(1)声波传播声压和距离计算沉积物的声速:
(1)
式中,V 原位为沉积物原位测量声速值,L1和L2为原位测量***发射换能器和两个接收换能器的距离,∆L是声波在海底沉积物中传播的距离差,t1和t2为声波穿过海底沉积物到达两个接收换能器的时间,∆t为声波到达两个接收换能器的时间差;
S106.原位测量***在海底测试时,利用温盐深仪CTD或者声速剖面仪SVP设备同步采集近海底海水的温度、声速或深度信息;
所述海底沉积物声学取样测量的步骤,包括:
S201.在相同站位进行海底沉积物取样工作,采样采用1套重力式取样器进行作业,当重力式取样器到达海底时,靠重力使采样筒***海底沉积物中,然后将取样器提升到船甲板上,从而获取沉积物柱状样品;
S202.利用声波信号发射采集***在实验室中对获取的沉积物柱状样品进行声学特性测试,测量中心测量频率分别为25kHz,31.5 kHz,50 kHz,80 kHz,及100kHz;声波信号发射采集***主要包括声波发射和采集电路、功率放大器、发射和接收换能器、柱状样品声学测试平台;在实验室采用轴向测试方法,首先将柱状样品进行分割成样品段,放置在声学测试平台上进行测试,发射换能器T和接收换能器R分别放置在柱状样品两端,然后发射声波,完成对L段沉积物样品段的声学测试,沉积物声速计算公式为:
(2)
式中:V 取样为沉积物取样测量声速值,单位为m/s;t为声信号初至时间,单位为μs;t0为换能器零声时修正值,单位为μs;L为沉积物的长度,单位为mm;实验室测量环境为标准大气压,温度为23℃;
所述沉积物声速比计算的步骤,包括:
S301.根据Mackenzie公式计算海水的声速值:
V(T, S, z) = a1 + a2T + a3T2 + a4T3 + a5(S - 35) + a6z + a7z2 + a8T(S -35) + a9Tz3
这里T表示温度,单位为摄氏度;S表示盐度,单位为千分比;z表示深度,单位为米;
S302.实验室取样测量的温度是T=23℃,盐度采用和海底海水一样的盐度S=S海底,深度z=0,代入到Mackenzie公式中,得到实验室海水的声速值V实验室海水;
S303.将步骤S1采集到的海底海水的温度T海底,盐度S海底,和深度z数据代入到Mackenzie公式中,可以得到海底海水的声速值V海底海水;
S304.分别计算实验室中的速度比和海底原位的速度比:
(3)
(4);
所述设定校正方程和校正因子的步骤,包括:
S401.执行步骤S1获得多种不同的沉积物类型各自的n个原位测量声速值,执行步骤S2获得多种不同的沉积物类型各自的n个海底海水速度值,执行步骤S3获得多种不同的沉积物类型各自的n个取样测量声速值和n个实验室海水速度值,执行步骤S3计算得到n个V实验室声速比和V海底声速比数据;
S402.选取步骤S401中的多个沉积物类型中的一个沉积物类型,将计算得到的这一类沉积物的n个的V实验室声速比和V海底声速比数据进行回归分析,构建以下校正方程:
(5)
其中为校正因子,根据一类沉积物的n个的V实验室声速比和V海底声速比数据可以计算得到这类沉积物的校正因子/>数值;
S403.重复S401和S402从而得到多个沉积物类型各自的校正方程和校正因子;
所述获取宽频域的校正因子曲线G(f)的步骤,包括:
S501.将不同测量频率25kHz、31.5 kHz、50 kHz,80 kHz及100kHz的原位测量值和取样测量值,分别重复执行步骤S1、S2和S3,获得不同频率的校正因子值;
S502.将不同频率的校正因子值与频率f=25kHz,31.5 kHz、50 kHz、80 kHz及100kHz进行拟合,得到一条/>和频率的回归曲线函数/>=G(f),回归曲线函数的频率范围是25kHz-100kHz;
所述将取样测量声速值校正到原位测量声速值的步骤,包括:
S601.首先利用取样样品在实验室中测量,获取得到多种沉积物类型中一类沉积物的取样测量声速值及实验室海水声速值,并计算得到实验室速度比,在取样测量时,同步测量海底海水的温度、盐度和深度信息;
S602.利用Mackenzie公式计算得到海底海水的声速值V海底海水;
S603.根据实验室取样测量频率f,结合回归曲线=G(f),确定校正因子/>的值,然后将海底海水速度值、取样测量声速值和实验室海水速度值带入到校正方程中,计算得到对应频率的海底原位声速值:
(6)
(7)
S604.将宽频域取样测量声速值校正到原位测量声速值,宽频域的频率范围是25kHz-100kHz。
2.根据权利要求1所述的宽频域海底沉积物取样测量声速值校正方法,其特征在于:
重力式取样器中的采样管的长度为6m;
重力式取样器中的压载铅块的重量为0.6t以上;
重力式取样器中的PVC衬管的内径为75-110mm;
常数a1, a2, ..., a9的取值为:
a1 = 1448.96, a2 = 4.591, a3 = -5.304×10-2, a4 = 2.374×10-4, a5 = 1.340,a6 = 1.630×10-2, a7 = 1.675×10-7, a8 = -1.025×10-2, a9 = -7.139×10-13;
n是不小于100的正整数;
V实验室声速比和V海底声速比的小数点位数保留6位以上。
3.一种计算机装置,其特征在于,包括存储器和处理器,所述存储器用于存储至少一个程序,所述处理器用于加载所述至少一个程序以执行权利要求1或2所述的一种宽频域海底沉积物取样测量声速值校正方法。
4.一种计算机可读存储介质,其中存储有处理器可执行的程序,其特征在于,所述处理器可执行的程序在由处理器执行时用于执行权利要求1或2所述的一种宽频域海底沉积物取样测量声速值校正方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310797298.2A CN116519799B (zh) | 2023-07-03 | 2023-07-03 | 宽频域海底沉积物取样测量声速值校正方法、装置和介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310797298.2A CN116519799B (zh) | 2023-07-03 | 2023-07-03 | 宽频域海底沉积物取样测量声速值校正方法、装置和介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116519799A CN116519799A (zh) | 2023-08-01 |
CN116519799B true CN116519799B (zh) | 2024-01-09 |
Family
ID=87403292
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310797298.2A Active CN116519799B (zh) | 2023-07-03 | 2023-07-03 | 宽频域海底沉积物取样测量声速值校正方法、装置和介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116519799B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09274081A (ja) * | 1996-04-03 | 1997-10-21 | Tech Res & Dev Inst Of Japan Def Agency | 海底堆積層における音速測定方法 |
CN103196998A (zh) * | 2013-04-18 | 2013-07-10 | 国家***第二海洋研究所 | 一种柱状样沉积物声速测量方法 |
KR101653650B1 (ko) * | 2016-03-08 | 2016-09-05 | 한국해양과학기술원 | 해저 퇴적층 음향특성 측정 장치 |
CN110210100A (zh) * | 2019-05-27 | 2019-09-06 | 中国科学院南海海洋研究所 | 一种高精度海底底质沉积物声速预测方法 |
CN112881525A (zh) * | 2021-01-18 | 2021-06-01 | 湖南国天电子科技有限公司 | 一种海底沉积物声学特性原位测量装置、***和方法 |
CN113191087A (zh) * | 2021-02-02 | 2021-07-30 | 中国人民解放军海军大连舰艇学院 | 联合历史剖面声速与实测表层声速的航渡式测深数据剖面声速改正方法 |
WO2021227605A1 (zh) * | 2020-05-13 | 2021-11-18 | 中国科学院海洋研究所 | 基于rov的深海沉积物声学参数原位探测***及方法 |
CN115639584A (zh) * | 2022-08-30 | 2023-01-24 | 中国石油大学(华东) | 一种海底基准站高精度定位方法、定位***及设备 |
CN115839998A (zh) * | 2023-01-06 | 2023-03-24 | 自然资源部第一海洋研究所 | 基于宽频换能器的海底沉积物样品声速测试***及装置 |
CN116186492A (zh) * | 2023-03-27 | 2023-05-30 | 中国人民解放军国防科技大学 | 一种深海水文环境数据拟合方法、装置及计算机介质 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
IL138073A0 (en) * | 2000-08-24 | 2001-10-31 | Glucon Inc | Photoacoustic assay and imaging system |
US6983208B2 (en) * | 2003-11-24 | 2006-01-03 | Mgd Technologies, Inc. | Method and apparatus for combined measurements of concentration, distribution and flow velocity of suspended solids |
US8170812B2 (en) * | 2007-10-16 | 2012-05-01 | Daniel Measurement And Control, Inc. | Method and system for detecting deposit buildup within an ultrasonic flow meter |
US9829595B2 (en) * | 2009-02-06 | 2017-11-28 | Westerngeco L.L.C. | Particle motion sensor-based streamer positioning system |
CN104411250B (zh) * | 2013-01-23 | 2017-06-23 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的动作方法 |
US11327201B2 (en) * | 2016-10-31 | 2022-05-10 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Porosity prediction based on effective stress |
US11661844B2 (en) * | 2020-10-07 | 2023-05-30 | Saudi Arabian Oil Company | Method and apparatus for fluid characterization and holdup estimation using acoustic waves |
-
2023
- 2023-07-03 CN CN202310797298.2A patent/CN116519799B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09274081A (ja) * | 1996-04-03 | 1997-10-21 | Tech Res & Dev Inst Of Japan Def Agency | 海底堆積層における音速測定方法 |
CN103196998A (zh) * | 2013-04-18 | 2013-07-10 | 国家***第二海洋研究所 | 一种柱状样沉积物声速测量方法 |
KR101653650B1 (ko) * | 2016-03-08 | 2016-09-05 | 한국해양과학기술원 | 해저 퇴적층 음향특성 측정 장치 |
CN110210100A (zh) * | 2019-05-27 | 2019-09-06 | 中国科学院南海海洋研究所 | 一种高精度海底底质沉积物声速预测方法 |
WO2021227605A1 (zh) * | 2020-05-13 | 2021-11-18 | 中国科学院海洋研究所 | 基于rov的深海沉积物声学参数原位探测***及方法 |
CN112881525A (zh) * | 2021-01-18 | 2021-06-01 | 湖南国天电子科技有限公司 | 一种海底沉积物声学特性原位测量装置、***和方法 |
CN113191087A (zh) * | 2021-02-02 | 2021-07-30 | 中国人民解放军海军大连舰艇学院 | 联合历史剖面声速与实测表层声速的航渡式测深数据剖面声速改正方法 |
CN115639584A (zh) * | 2022-08-30 | 2023-01-24 | 中国石油大学(华东) | 一种海底基准站高精度定位方法、定位***及设备 |
CN115839998A (zh) * | 2023-01-06 | 2023-03-24 | 自然资源部第一海洋研究所 | 基于宽频换能器的海底沉积物样品声速测试***及装置 |
CN116186492A (zh) * | 2023-03-27 | 2023-05-30 | 中国人民解放军国防科技大学 | 一种深海水文环境数据拟合方法、装置及计算机介质 |
Non-Patent Citations (3)
Title |
---|
"海底表层沉积物声速的环境因素影响特性";邹大鹏 等;《声学学报》;第46卷(第2期);227-236 * |
"深海海底沉积物声学特性原位测量试验研究";王景强 等;《地球物理学报》;第63卷(第12期);4463-4472 * |
海底沉积物声学性质原位测量***海上试验研究;李艳华;郭常升;李会银;谭宝海;;海洋科学;34(09);55-60 * |
Also Published As
Publication number | Publication date |
---|---|
CN116519799A (zh) | 2023-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116520431B (zh) | 浅表海底沉积物宽频分层声速结构构建方法、装置和介质 | |
KR101435617B1 (ko) | 해저 지질 조사용 현장 음파 전달속도 측정 시스템 | |
US10101304B2 (en) | Phased array system capable of computing gains for non-measured calibration points | |
CN101473196A (zh) | 测量装置及其使用方法 | |
EA018239B1 (ru) | Акустическое измерение толщины с использованием газа в качестве связующей среды | |
Wang et al. | Frequency dependence of sound speed and attenuation in fine-grained sediments from 25 to 250 kHz based on a probe method | |
CN109991590B (zh) | 一种在有限空间压力罐内测试换能器低频发射特性的***与方法 | |
CN116910473B (zh) | 海底沉积物宽频声速预测方法、计算机装置和存储介质 | |
CN116738182B (zh) | 海底沉积物声衰减范围值的预测方法、装置和存储介质 | |
CN106383173B (zh) | 一种水泥声阻抗计算方法和装置 | |
CN116519799B (zh) | 宽频域海底沉积物取样测量声速值校正方法、装置和介质 | |
JPWO2007132509A1 (ja) | 共振現象を利用した超音波探査方法 | |
CN116644295B (zh) | 海底沉积物声衰减宽频预测方法、计算机装置和存储介质 | |
Lecours et al. | Assessing the spatial data quality paradox in the deep-sea | |
CN111665296A (zh) | 基于emat测量超声换能器三维辐射声场的方法及装置 | |
CN102901472B (zh) | 单波束测深仪检测方法 | |
CN103261525B (zh) | 用于确定装料区内混合物中的沉积层体积的方法和装置 | |
CN109297530A (zh) | 一种深海全水柱多环境要素信息融合方法及处理终端 | |
CN104790939A (zh) | 一种获取胶结率的方法和装置 | |
CN113091878A (zh) | 用于水下声速测量的装置及其使用方法 | |
RU2563603C1 (ru) | Способ определения чувствительности по полю гидроакустического приемника | |
JP4500973B2 (ja) | コンクリート構造物の圧縮強度測定方法及び測定装置 | |
JP2007205848A (ja) | 超音波法によるコンクリート構造物のひび割れ深さ探査方法及びそのひび割れ深さ探査装置 | |
CN115951423B (zh) | 声波测井接收换能器一致性检测方法、装置、设备及介质 | |
Chiu et al. | Error analysis on normal incidence reflectivity measurement and geoacoustic inversion of ocean surficial sediment properties |
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 |