CN105223610B - 判断海底地震仪与航行舰船地震波信号耦合效果的方法 - Google Patents

判断海底地震仪与航行舰船地震波信号耦合效果的方法 Download PDF

Info

Publication number
CN105223610B
CN105223610B CN201510696626.5A CN201510696626A CN105223610B CN 105223610 B CN105223610 B CN 105223610B CN 201510696626 A CN201510696626 A CN 201510696626A CN 105223610 B CN105223610 B CN 105223610B
Authority
CN
China
Prior art keywords
submarine seismograph
coefficient
seismic
submarine
seismograph
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.)
Expired - Fee Related
Application number
CN201510696626.5A
Other languages
English (en)
Other versions
CN105223610A (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.)
760 RESEARCH INSTITUTE OF CSIC
Original Assignee
760 RESEARCH INSTITUTE OF CSIC
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 760 RESEARCH INSTITUTE OF CSIC filed Critical 760 RESEARCH INSTITUTE OF CSIC
Priority to CN201510696626.5A priority Critical patent/CN105223610B/zh
Publication of CN105223610A publication Critical patent/CN105223610A/zh
Application granted granted Critical
Publication of CN105223610B publication Critical patent/CN105223610B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

一种判断海底地震仪耦合效果的方法,属于海底地震仪检测航行舰船地震波信号能力的分析技术。其特征是:利用实验室实验和数据分析,确定海底地震仪与海底介质的振动系统耦合系数、振动系统阻尼系数、振动系统固有角频率和海底地震仪横切面比值;利用主观赋权法并结合航行舰船地震波信号特征确定耦合系数、阻尼系数、固有角频率和横切面比值的权重系数;基于航行舰船地震波信号特征,确定判定参数中4个判定值的取值,得到判定参数的计算方法;通过分析判定参数的计算结果,最终确定海底地震仪的耦合效果。

Description

判断海底地震仪与航行舰船地震波信号耦合效果的方法
技术领域
本发明属于航行舰船地震波探测领域,涉及到一种判断海底地震仪与航行舰船地震波信号耦合效果的方法,特别涉及一种利用耦合效果判定参数综合评估海底地震仪测试航行舰船地震波信号能力的方法。
背景技术
海底地震仪是测试航行舰船地震波信号的主要设备,传统的测试方式是将海底地震仪布放于海水与海底的分界面上,然而由于海底地震仪的设计不合理容易导致所测试的航行舰船地震波信号失真。目前国内针对海底地震仪的耦合效果判断仅参考耦合理论提出定性评价,缺乏基于判定参数的定量分析,从而导致试验所测数据严重失真,无法有效提取信号中的信息。
关于海底地震仪与海底介质的耦合理论已经比较成熟,国内外相关研究人员将海底地震仪假设为一刚体,放在均匀、各相同性的弹性半无限空间里。运动方式具有6个自由度,沿垂直轴与两个正交水平轴的运动和转动。对于每一个自由度都建立了特征频率和阻尼表示的振动模型,而且模型之间还存在相互作用。然而目前所建立的理论模型比较理想化,并且未结合实际的海底地震仪结构以及航行舰船地震波信号的特点,因此未能定量判定耦合效果。
发明内容
本发明的目的是提供一种判断海底地震仪与航行舰船地震波信号耦合效果的方法,解决海底地震仪检测航行舰船地震波信号能力的定量评价问题。
本发明的技术方案是:海底地震仪与航行舰船地震波信号的耦合效果判定参数由4个因素组成,4个因素分别为振动系统耦合系数、振动系统阻尼系数、振动系统固有角频率和海底地震仪横切面比值;振动系统耦合系数由海底地震仪排开水的质量、海底地震仪质量以及与海底地震仪一起运动的水的质量计算得到;振动系统阻尼系数和固有角频率在实验室测试得到,释放系留于海底地震仪上的浮球分别产生垂直向上和水平方向的力,通过分析实测数据计算得到阻尼系数和固有角频率;海底地震仪的横切面比值是通过测量海底地震仪的底宽和高度相比得到;依据每个因素在航行舰船地震波信号检测中的关键程度分别赋予一定的权重系数,最终计算出评价耦合效果的判定参数。
本发明的效果和益处是:
传统的判断耦合效果的方法主要是基于经验和理论分析,给出海底地震仪设计的参考意见,并不能准确的判定耦合效果;利用计算得到的判定参数,可以定量分析不同海底地震仪对航行舰船地震波信号的检测能力,从而提出结构设计的改进意见。
附图说明
图1海底地震仪阻尼系数和固有角频率测试示意图
图中:1实验水池;2海水;3定滑轮;4海底介质;5沉耦架;6海底地震仪;7浮球。
图2海底地震仪耦合系数测试示意图。
图中:1实验水池;2海底地震仪;3沉耦架;4海水;5海底介质。
具体实施方式
下面结合技术方案和附图详细叙述本发明的具体实施方式。
第一步:确定阻尼系数和固有角频率
阻尼系数和固有角频率测试实验如图1所示,通过实验室实验和测试数据的后处理确定海底地震仪的阻尼系数和固有角频率。实验水池底部面积大于海底地震仪底部面积的4倍,在海底地震仪将要布放的海域提取海底介质和海水样本,海底介质的体积大于实验水池底部面积与0.5米的乘积,海水介质的体积大于实验水池底部面积与海底地震仪1.5倍高度的乘积。在实验水池两个相互垂直的池壁上各安装1个定滑轮,定滑轮的高度为海底地震仪高度的1/2与0.5米之和,定滑轮距离所安装的池壁0.5米。在实验水池中铺上0.5米厚度的海底介质样本,设置海底地震仪开始采集数据,将海底地震仪放置与海底介质上。海底地震仪距离未安装定滑轮的两个池壁都为0.5米,海底地震仪中地震波传感器的两个水平方向分别垂直于安装定滑轮的池壁。在海底地震仪上系3条牵引绳,其中1条在海底地震仪正上方,另外2条为水平方向并通过定滑轮引出,2条水平牵引绳平行于水池池底。在实验水池中注入海水样本,水深为海底地震仪高度的1.5倍。在3条牵引绳末端系浮球,浮球直径250毫米,浮球完全浸入水中。定义地震波传感器两个水平方向为X方向和Y方向,定义地震波传感器垂直方向为Z方向,按传感器X、Y和Z方向的顺序释放浮球同时记录各自释放的时间。对海底地震仪的测试数据进行幅度谱分析,找到幅度谱峰值Amax所对应的频率f0,在幅度谱峰值两侧找到幅度对应的频率f1和f2,此时就可以计算得到海底地震仪的阻尼系数hc和固有角频率ω0和ω0=2πf0
第二步:计算耦合系数和横切面比值
耦合系数测试实验如图2所示,耦合系数由结构质量、结构排开水的质量及与海底地震仪一起运动的水的质量计算得到,横切面比值由海底地震仪的底宽和高度计算得到。将放置与实验水池中的海底地震仪提起,完全离开水面。测量海底地震仪的质量为Mo、底宽为Do、高度为Ho,测量海水的密度为ρw,测量实验水池的长度为L、宽度为D、水深为h0。将海底地震仪放置于实验水池的中心位置,再次测量实验水池水深为h1。计算得到结构排开水的质量Mw为ρw·L·D·(h1-h0),与海底地震仪一起运动的水的质量为0.5·ρw·L·D·(h1-h0),综合以上计算结果,可以得到耦合系数C为:
横切面比值Ro
第三步:确定权重系数
结合航行舰船地震波信号特征确定耦合系数、阻尼系数、固有角频率和横切面比值的权重系数,采用的赋权方法为主观赋权法。耦合系数表征海底地震仪与航行舰船地震波信号是否能够一起运动,直接决定海底地震仪与海底介质的耦合效果,权重系数设置为0.3;阻尼系数表征海底地震仪接收航行舰船地震波或其它干扰后是否能够快速回到稳定状态,决定海底地震仪的抗干扰能力,权重系数设置为0.15;固有角频率表征海底地震仪在哪些频点最敏感即引起的振动最剧烈,这是在检测航行舰船地震波信号时最应该引起重视的问题,它往往导致航行舰船地震波信号无法有效检测,权重系数设置为0.45;横切面比值表征海底地震仪遇到海水冲击时的摇摆程度,它可以引起海底地震仪中地震波传感器三个方向之间的互相耦合,设置权重系数为0.1。
第四步:确定判定参数
判定参数由耦合系数、阻尼系数、固有角频率和横切面比值共同确定,根据权重系数设置,得到判定参数Df的表达式为:
Df=0.3·C′+0.15·h′c+0.45·ω0′+0.1·R′o
式中,C′为耦合系数判定值,h′c为阻尼系数判定值,ω0′为固有角频率判定值,R′o为横切面比值判定值。
C为非负数,耦合系数判定值C′取值为:
当0≤C<0.1时,C′=1;当0.1≤C<0.5时,C′=0.5;当0.5≤C<5时,C′=0.1;当5≤C时,C′=0。
hc的最大值为1,阻尼系数判定值h′取值为:
当0.5≤hc≤1时,h′c=1;当0.1≤hc<0.5时,h′c=0.2;当hc<0.1时,h′c=0。
f0为非负数,固有角频率判定值ω0′取值为:
当f0>100Hz时,ω0′=1;当10Hz<f0≤100Hz时,ω0′=0.5;当5Hz<f0≤10Hz时,ω0′=0.2;当0Hz≤f0≤5Hz时,ω0′=0。
Ro为正数,横切面比值判定值R′o为:
当Ro>5时,R′o=1;当1<Ro≤5时,R′o=0.5;当0.5<Ro≤1时,R′o=0.2;当0<Ro≤0.5时,R′o=0。
如果判定参数中的4个判定值任意一个为0,则设定判定参数为0;如果4个判定值均不为0,则由判定值乘以各自的权重系数,得到判定参数。
第五步:判断耦合效果
判定参数介于0和1之间,值越大表示耦合效果越好。判定参数最小值为0表示所设计的海底地震仪耦合效果很差,不能用于信号检测。判定参数最大值为1表示所设计的海底地震仪耦合效果很好,可以检测真实信号。判定参数大于或等于0.6表示海底地震仪可以使用,但信号存在一定程度的失真;判定参数小于0.6表示海底地震仪不可使用,需要改进。

Claims (1)

1.一种判断海底地震仪与航行舰船地震波信号耦合效果的方法,其特征是:
一、确定阻尼系数和固有角频率:
通过实验室实验和测试数据的后处理确定海底地震仪的阻尼系数和固有角频率;实验水池底部面积大于海底地震仪底部面积的4倍,在海底地震仪将要布放的海域提取海底介质和海水样本,海底介质的体积大于实验水池底部面积与0.5米的乘积,海水介质的体积大于实验水池底部面积与海底地震仪1.5倍高度的乘积;在实验水池两个左右相互垂直的池壁上各安装1个定滑轮,即共安装2个定滑轮,定滑轮的高度为海底地震仪高度的1/2与0.5米之和,定滑轮距离所安装的池壁0.5米;在实验水池中铺上0.5米厚度的海底介质样本,设置海底地震仪开始采集数据,将海底地震仪放置在海底介质上;海底地震仪距离未安装定滑轮的两个池壁都为0.5米,海底地震仪中地震波传感器的两个水平方向分别垂直于安装定滑轮的池壁;在海底地震仪上系3条牵引绳,其中1条在海底地震仪正上方,另外2条为水平方向并通过定滑轮引出,2条水平牵引绳平行于实验水池池底;在实验水池中注入海水介质,水深为海底地震仪高度的1.5倍;在3条牵引绳末端系浮球,浮球直径250毫米,浮球完全浸入水中;定义地震波传感器两个水平方向为X方向和Y方向,定义地震波传感器垂直方向为Z方向,按地震波传感器X、Y和Z方向的顺序释放浮球同时记录各自释放的时间;对海底地震仪的测试数据进行幅度谱分析,找到幅度谱峰值Amax所对应的频率f0,在幅度谱峰值两侧找到幅度对应的频率f1和f2,此时就可以计算得到海底地震仪的阻尼系数hc和固有角频率ω0和ω0=2πf0
二、计算耦合系数和横切面比值:
耦合系数由结构质量、结构排开水的质量及与海底地震仪一起运动的水的质量计算得到,横切面比值由海底地震仪的底宽和高度计算得到;将放置在实验水池中的海底地震仪提起,完全离开水面;测量海底地震仪的质量为Mo、底宽为Do、高度为Ho,测量海水的密度为ρw,测量实验水池的长度为L、宽度为D、水深为h0;将海底地震仪放置于实验水池的中心位置,再次测量实验水池水深为h1;计算得到结构排开水的质量Mw为ρw·L·D·(h1-h0),与海底地震仪一起运动的水的质量M′w为0.5·ρw·L·D·(h1-h0),得到耦合系数C为:
横切面比值Ro
三、确定权重系数:
结合航行舰船地震波信号特征确定耦合系数、阻尼系数、固有角频率和横切面比值的权重系数,采用的赋权方法为主观赋权法;
(1)耦合系数表征海底地震仪与航行舰船地震波信号是否能够一起运动,直接决定海底地震仪与海底介质的耦合效果,权重系数设置为0.3;
(2)阻尼系数表征海底地震仪接收航行舰船地震波信号或其它干扰后是否能够快速回到稳定状态,决定海底地震仪的抗干扰能力,权重系数设置为0.15;
(3)固有角频率表征海底地震仪在哪些频点最敏感,即引起的振动最剧烈,这是在检测航行舰船地震波信号时最应该引起重视的问题,它往往导致航行舰船地震波信号无法有效检测,权重系数设置为0.45;
(4)横切面比值表征海底地震仪遇到海水冲击时的摇摆程度,它可以引起海底地震仪中地震波传感器三个方向之间的互相耦合,设置权重系数为0.1;
四、确定判定参数:
判定参数由耦合系数、阻尼系数、固有角频率和横切面比值共同确定,根据权重系数设置,得到判定参数Df的表达式为:
Df=0.3·C′+0.15·h′c+0.45·ω′0+0.1·R′o
式中,C′为耦合系数判定值,h′c为阻尼系数判定值,ω′0为固有角频率判定值,R′o为横切面比值判定值;
耦合系数C为非负数,耦合系数判定值C′取值为:
当0≤C<0.1时,C′=1;当0.1≤C<0.5时,C′=0.5;当0.5≤C<5时,C′=0.1;当5≤C时,C′=0;
阻尼系数hc的最大值为1,阻尼系数判定值h′取值为:
当0.5≤hc≤1时,h′c=1;当0.1≤hc<0.5时,h′c=0.2;当hc<0.1时,h′c=0;
固有角频率f0为非负数,固有角频率判定值ω′0取值为:
当f0>100Hz时,ω′0=1;当10Hz<f0≤100Hz时,ω′0=0.5;当5Hz<f0≤10Hz时,ω′0=0.2;当0Hz≤f0≤5Hz时,ω′0=0;
横切面比值Ro为正数,横切面比值判定值R′o为:
当Ro>5时,R′o=1;当1<Ro≤5时,R′o=0.5;当0.5<Ro≤1时,R′o=0.2;当0<Ro≤0.5时,R′o=0;
如果判定参数中的4个判定值任意一个为0,则设定判定参数为0;如果4个判定值均不为0,则由判定值乘以各自的权重系数,得到判定参数;
五、判断耦合效果:
(1)判定参数介于0和1之间,值越大表示耦合效果越好;
(2)判定参数最小值为0表示所设计的海底地震仪耦合效果很差;
(3)判定参数最大值为1表示所设计的海底地震仪耦合效果很好;
(4)判定参数大于或等于0.6表示海底地震仪可以使用,但信号存在一定程度的失真;
(5)判定参数小于0.6表示海底地震仪不可使用,需要改进。
CN201510696626.5A 2015-10-22 2015-10-22 判断海底地震仪与航行舰船地震波信号耦合效果的方法 Expired - Fee Related CN105223610B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510696626.5A CN105223610B (zh) 2015-10-22 2015-10-22 判断海底地震仪与航行舰船地震波信号耦合效果的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510696626.5A CN105223610B (zh) 2015-10-22 2015-10-22 判断海底地震仪与航行舰船地震波信号耦合效果的方法

Publications (2)

Publication Number Publication Date
CN105223610A CN105223610A (zh) 2016-01-06
CN105223610B true CN105223610B (zh) 2017-08-15

Family

ID=54992665

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510696626.5A Expired - Fee Related CN105223610B (zh) 2015-10-22 2015-10-22 判断海底地震仪与航行舰船地震波信号耦合效果的方法

Country Status (1)

Country Link
CN (1) CN105223610B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108873058B (zh) * 2018-04-27 2021-06-11 中国船舶工业系统工程研究院 一种海底地震仪样机及实测方法
CN108897040B (zh) * 2018-05-10 2020-05-08 中国石油天然气集团有限公司 一种确定陆检与海底之间的耦合性的方法及装置
CN108957522B (zh) * 2018-09-26 2019-03-19 中国海洋大学 一种基于克希霍夫积分的海底地震仪数据多次波预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4951264A (en) * 1986-05-16 1990-08-21 University Of Miami Method of measuring the shear modulus profile of a seabed
CN102759395A (zh) * 2012-07-16 2012-10-31 浙江大学 浮球式振动波传感器及其使用方法与用途
CN103257347A (zh) * 2012-02-17 2013-08-21 中国人民解放军海军装备研究院舰艇作战系统论证研究所 一种实际使用环境下的声纳作用距离指标的获取方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013003081A1 (en) * 2011-06-26 2013-01-03 The Regents Of The University Of California An inexpensive instrument for measuring wave exposure and water velocity

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4951264A (en) * 1986-05-16 1990-08-21 University Of Miami Method of measuring the shear modulus profile of a seabed
CN103257347A (zh) * 2012-02-17 2013-08-21 中国人民解放军海军装备研究院舰艇作战系统论证研究所 一种实际使用环境下的声纳作用距离指标的获取方法
CN102759395A (zh) * 2012-07-16 2012-10-31 浙江大学 浮球式振动波传感器及其使用方法与用途

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
舰船地震波环境噪声的统计特性;崔培 等;《舰船科学技术》;20110331;第33卷(第3期);第12-15页 *

Also Published As

Publication number Publication date
CN105223610A (zh) 2016-01-06

Similar Documents

Publication Publication Date Title
CN104360046B (zh) 一种码头混凝土结构内部隐患的综合物探联合诊断方法
Lanckriet et al. Comprehensive field study of swash-zone processes. II: Sheet flow sediment concentrations during quasi-steady backwash
CN105223610B (zh) 判断海底地震仪与航行舰船地震波信号耦合效果的方法
CN103630931B (zh) 从近场测量和建模假想特征计算假想源特征的方法和系统
Faraci et al. Wave-and current-dominated combined orthogonal flows over fixed rough beds
CN105486487A (zh) 一种波浪检测系统
An et al. Detecting local scour using contact image sensors
CN111381292B (zh) 一种预测砂岩含烃储层的测井解释方法与装置
Yu et al. Laboratory evaluation of time‐domain reflectometry for bridge scour measurement: comparison with the ultrasonic method
CN112835124A (zh) 基于成像测井和阵列声波测井资料的裂缝有效性评价方法
CN108508093A (zh) 一种工件缺陷高度的检测方法及系统
US20230251221A1 (en) In-situ evaluation method and system for loess collapsibility based on non-destructive time-domain reflection technology
CN215180930U (zh) 海底底质小尺度原位声学成像系统
Abelev et al. Numerical investigation of dynamic free-fall penetrometers in soft cohesive marine sediments using a finite difference approach
CN115980758A (zh) 一种集成于海底原位装备的泥线位置测试仪及识别方法
CN205920228U (zh) 一种走航式深海天然气水合物探测取样装置
Truong et al. Geophysical monitoring techniques for underwater landslide in 1g models
CN112904409B (zh) 一种基于分布式声波传感技术的水下沉积物沉积特性测量装置与方法
Ertel et al. Advances in pile integrity testing
Mukhlas et al. Bag and floater motions of a fabric membrane cage
Perillo et al. Acoustic imaging of experimental subaqueous sediment-laden flows and their deposits
Zhu et al. Damage detection of pile foundations in marine engineering based on a multidimensional dynamic signature
CN219162374U (zh) 一种集成于海底原位装备的泥线位置测试仪
CN205352655U (zh) 一种入、反射波时域分离系统
Zhang et al. Comparison of turbulence schemes for prediction of wave-induced near-bed sediment suspension above a plane bed

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Cui Pei

Inventor after: Shao Jun

Inventor after: Bai Chunzhi

Inventor after: Jiang Kaina

Inventor after: Fu Yibo

Inventor after: Wu Yunchao

Inventor after: Wang Xu

Inventor before: Cui Pei

Inventor before: Bai Chunzhi

Inventor before: Wu Yunchao

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
CB03 Change of inventor or designer information

Inventor after: Cui Pei

Inventor after: Shao Jun

Inventor after: Bai Chunzhi

Inventor after: Jiang Kaina

Inventor after: Fu Yibo

Inventor after: Wu Yunchao

Inventor after: Wang Xu

Inventor after: Li Biao

Inventor before: Cui Pei

Inventor before: Shao Jun

Inventor before: Bai Chunzhi

Inventor before: Jiang Kaina

Inventor before: Fu Yibo

Inventor before: Wu Yunchao

Inventor before: Wang Xu

CB03 Change of inventor or designer information
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170815

Termination date: 20171022

CF01 Termination of patent right due to non-payment of annual fee