CN108107437B - 一种利用简正波耦合干涉的海洋环境监测方法 - Google Patents
一种利用简正波耦合干涉的海洋环境监测方法 Download PDFInfo
- Publication number
- CN108107437B CN108107437B CN201710894274.3A CN201710894274A CN108107437B CN 108107437 B CN108107437 B CN 108107437B CN 201710894274 A CN201710894274 A CN 201710894274A CN 108107437 B CN108107437 B CN 108107437B
- Authority
- CN
- China
- Prior art keywords
- marine environment
- interference
- order
- sound
- wave
- 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
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000008878 coupling Effects 0.000 title claims abstract description 20
- 238000010168 coupling process Methods 0.000 title claims abstract description 20
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 20
- 238000012544 monitoring process Methods 0.000 title claims abstract description 17
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 238000001914 filtration Methods 0.000 claims abstract description 5
- 230000005236 sound signal Effects 0.000 claims abstract description 3
- 230000021615 conjugation Effects 0.000 abstract description 4
- 230000002349 favourable effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 4
- 241001386813 Kraken Species 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
Abstract
本发明提供的是一种利用简正波耦合干涉的海洋环境监测方法。发射宽带声信号，得到海洋环境发生变化时，接收点处的声压场；利用模态滤波方法，将声压场进行模态分解，得到各阶简正波的声压分量；将各阶简正波声压分量和自身共轭相乘，得到各阶简正波和自身共轭的干涉图样；根据干涉图样，利用波导不变量法得到干涉图样中的干涉条纹斜率，并根据这一斜率和海洋环境变化的有关参数的关系，对海洋环境变化进行监测。本发明能够稳定、快速地计算实际中具有时空变化特性的海洋环境下的声场。本发明能够从复杂的声场干涉图样中得到清晰简单且物理意义明确的干涉图样更易于提取海洋环境变化的特征信息，有助于实现海洋环境变化的预警和监测。
Description
技术领域
本发明涉及的是一种海洋耦合简正波干涉特性以及海洋环境变化监测方法。
背景技术
声波在海洋中传播，产生了声场。实际海洋环境是时变和空变的，这会使得在其中传播的声信号的幅度和相位产生畸变，从而使接收到的声信号具有时频干涉和空频干涉特性，而这些干涉特性能够在一定程度上反映出海洋环境的变化情况及相关参数。然而，现有方法很难从声场中获取海洋环境变化的有关信息，因此，无法预报和测量引起这些海洋环境变化的有关因素。
为了从声场中提取海洋环境信息，一系列的声学模型和声场算法被相继提出，如声场的时频干涉、波导不变量分析法、波数积分方法、谱分解等。这些方法能够在一定程度上提取海洋环境变化的有关信息，但是往往依托于整个声场的信息，所得结果较为复杂，难以进行直观的分析。而各阶简正波之间的干涉能够得到清晰的干涉图样和明确的物理意义，便于分析。所以，为了得到清晰简单的干涉图样，以便对结果进行分析，对每一阶简正波的干涉结构的分析研究是很有必要的。
发明内容
本发明的目的在于提供一种操作性强，计算量适中，结果清晰、稳定且物理意义明确的利用简正波耦合干涉的海洋环境监测方法。
本发明的目的是这样实现的：
(1)发射宽带声信号，得到海洋环境发生变化时，接收点处的声压场p(R,z；f,r0)，其中R为接收器到声源的距离，z为接收点所在深度，f为声波频率，r0为孤立子内波到声源的距离；
(2)利用模态滤波方法，将步骤(1)得到的声压场p(R,z；f,r0)进行模态分解，得到各阶简正波的声压分量pn(R,z；f,r0)；
(3)将步骤(2)得到的各阶简正波声压分量pn(R,z；f,r0)和自身共轭相乘，得到各阶简正波和自身共轭的干涉图样
(4)根据步骤(3)得到的干涉图样，利用波导不变量法得到干涉图样中的干涉条纹斜率，并根据这一斜率和海洋环境变化的有关参数的关系，对海洋环境变化进行监测。
本发明提供了一种高效准确的海洋声场监测技术，能够从声场中充分提取海洋环境信息并由此监测海洋环境变化，如内波、锋面、涡旋、洋流等对声场带来的影响。现有的声场计算技术虽然能够进行声场计算和海洋环境信息提取，但往往依托于整个声场的信息，而不关注每一阶简正波的干涉结构，因此存在于各阶简正波之间的干涉中的重要信息往往会被忽略。本发明针对每一阶简正波的干涉结构，能够稳定、快速地从中提取海洋环境变化的有关信息，物理意义明确，图像清晰。相比其他声场计算技术，本发明能在复杂的声场干涉中，得到相对简单清晰且物理意义明确的简正波干涉图样，极大简化对声场干涉的分析以及海洋环境变化信息的提取。
本发明所采用的技术是，在海洋环境变化的情况下，将简正波耦合、模态滤波、简正波干涉和波导不变量法结合起来计算和分析声场，以达到快速简单提取海洋环境变化信息，并对海洋环境变化进行预警和监测的目的。
本发明基于的原理如下：
声场中的声压
其中，f为声波频率，r0为海洋环境变化量所在位置到声源的距离，zr为接收水听器所在深度，M为简正波的阶数，R为接收水听器到声源的距离，an(f；r0,zr,R)为第n阶简正波的幅度，φn(zr)为接收点处第n阶简正波的本征函数。
对海上试验测得的宽带声压场做模态分解，得到各阶简正波的声压分量
pn(f；r0,zr,R)＝an(f；r0,zr,R)φn(zr)(2)
由于海洋环境会随着时间和空间距离发生变化，这个变化会引起各阶简正波之间的耦合干涉，从而改变各阶简正波的幅度an(f；r0,zr,R)、各阶简正波的声压pn(f；r0,zr,R)以及声场的时频、空频干涉结构。为了得到声场的空频干涉特性，将各阶简正波的声压分量
pn(f；r0,zr,R)与其共轭
相乘，得到
海洋环境变化带来的简正波耦合干涉的空频干涉特性pnn(f；r0,zr,R)包含着海洋环境随着时间和空间变化的一些特征信息。因此，可以通过分析声场的耦合干涉特性pnn(f；r0,zr,R)，获得海洋环境的时变、空变特性，进而对海洋环境的变化进行监测。
接下来，以内波为例来说明，如何通过海洋环境变化引起简正波耦合的干涉特性，来进行海洋环境监测：
对于海洋环境中发射的宽带声信号，第n阶简正波的幅度
声压场
其中，f为声波频率，r0为孤立子内波相对于声源的水平距离，zr为接收水听器所在深度， M为简正波的阶数，R为接收水听器到声源的距离，kn为第n阶简正波的本征值，φn(zr)为第n阶简正波在接收水听器处的本征函数，An(R)为第n阶简正波在距离R处的幅度，满足
其中ln＝kn+iαn，为第n阶简正波的复本征值，αn为第n阶简正波的衰减系数，An(0)为声源处的第n阶简正波幅度，η0为孤立子内波的幅度，Znm为第n阶简正波和第m阶简正波之间的耦合矩阵，lnm＝ln-lm＝-lmn，knm＝kn-km＝-kmn，αnm＝αn-αm＝-αmn。
所以，计算结果保留到η0的一阶近似，第n阶简正波声压分量和自身共轭的乘积
其中
(3)式所表示的第n阶简正波和自身共轭乘积，由于内波的存在而导致各阶简正波之间发生了耦合干涉，在声波频率f-孤立子内波位置r0平面有明显的干涉条纹现象。
干涉图样的条纹斜率
经过近似处理，第一阶简正波的干涉条纹斜率
其中
Sp12(f)为声波频率为f时，第一阶简正波和第二阶简正波的相慢度差，即有
Sp12(f)＝Sp1(f)-Sp2(f) (10)
Sg12(f)为声波频率为f时，第一阶简正波和第二阶简正波的群慢度差，即有
Sg12(f)＝Sg1(f)-Sg2(f) (11)
根据(6)式，可以在干涉条纹斜率
Sp12(f)、Sg12(f)以及声波频率f已知的情况下，来估计孤立子内波的位置r0，有
根据
随着时间的变化，还可以估计孤立子内波沿着声源-接收点连线方向的运动速度。
因此，可以根据简正波的耦合干涉特性pnn(f；r0,zr,R)来监测海洋环境变化。
本发明有以下优势：
1、本发明能够稳定、快速地计算实际中具有时空变化特性的海洋环境下的声场。
2、本发明能够从复杂的声场干涉图样中得到清晰简单且物理意义明确的干涉图样更易于提取海洋环境变化的特征信息，有助于实现海洋环境变化的预警和监测。
附图说明
图1为发明内容程序流程图；
图2a-图2b为文中采用的海洋环境；图2a声速剖面，图2b孤立子内波形状；
图3a-图3b为声压场；图3a为KRAKEN声压场、图3b PE声压场；
图4a-图4b为第一阶简正波声压和自身共轭乘积；图4a为KRAKEN第一阶简正波声压、图 3b PE第一阶简正波声压；
图5为图4中的条纹斜率δr0/δf|1st。
具体实施方式
本发明提供的是一种利用简正波耦合干涉的海洋环境监测技术。所述的简正波耦合的干涉特性，由于受到实际海洋环境具有时变和空变特性的影响，在其中传播的声信号会产生简正波耦合干涉特性。
实际时变空变海洋环境下的声场中的声压
其中，f为声波频率，r0为海洋环境变化量所在位置到声源的距离，zr为接收水听器所在深度，M为简正波的阶数，R为接收水听器到声源的距离，an(f；r0,zr,R)为第n阶简正波的幅度，φn(zr)为接收点处第n阶简正波的本征函数。
对海上试验测得的宽带声压场做模态分解，得到各阶简正波的声压分量
pn(f；r0,zr,R)＝an(f；r0,zr,R)φn(zr) (2)
由于海洋环境会随着时间和空间距离发生变化，这个变化会引起各阶简正波之间的耦合干涉，从而改变各阶简正波的幅度an(f；r0,zr,R)、各阶简正波的声压pn(f；r0,zr,R)以及声场的时频、空频干涉结构。为了得到声场的空频干涉特性，我们用各阶简正波的声压分量 pn(f；r0,zr,R)和自身的共轭
相乘，得到
所述的海洋环境监测技术，利用海洋环境变化带来的简正波耦合干涉的空频干涉特性 pnn(f；r0,zr,R)，提取其中包含的海洋环境特征信息，获得海洋环境的时变、空变特性，进而对海洋环境的变化进行监测。
下面举例对本发明做更详细的描述。
(1)本发明利用编程实现功能的操作流程如图1所示。首先利用KRAKEN/PE方法计算特定海洋环境下的声场。本质上就是计算特定海洋环境下，信道的单位冲击相应。即根据图 2a-图2b所示的海洋环境，求出声压场p(f；r0,zr)(如图3a-图3b所示)，f是声波频率，r0是孤立子内波位置，zr为接收点所在深度。
(2)用模态滤波法，对(1)求出的声压场进行模态分解，得到各阶简正波的声压分量 pn(f；r0,zr)，其中n＝1,2,...,M，M为计算所考虑的简正波总阶数。然后计算pn(f；r0,zr)和自身共轭
的乘积
以第一阶简正波为例，所得结果如图4a-图4b所示。
(3)用波导不变量法，计算图4a-图4b中的干涉条纹斜率
以第一阶简正波为例，得到
根据这个关系，即可得到孤立子内波位置r0的估计值
(4)对于真实的海洋环境数据p(f；t,zr)，t为接收器从开始接收信号所经过的时间，可以根据上述方法得到各阶简正波的声压分量和自身的共轭乘积
然后求得干涉条纹斜率，并根据温度链测得孤立子内波的运动速度v，进而估计孤立子内波到声源的距离r0，从而对孤立子内波和其他的海洋环境变化进行预报和监测。
Claims (1)
1.一种利用简正波耦合干涉的海洋环境监测方法，其特征是：
(1)发射宽带声信号，得到海洋环境发生变化时，接收点处的声压场p(R,z；f,r0)，其中R为接收器到声源的距离，z为接收点所在深度，f为声波频率，r0为孤立子内波到声源的距离；
(2)利用模态滤波方法，将步骤(1)得到的声压场p(R,z；f,r0)进行模态分解，得到各阶简正波的声压分量pn(R,z；f,r0)；
(3)将步骤(2)得到的各阶简正波声压分量pn(R,z；f,r0)和自身共轭相乘，得到各阶简正波和自身共轭的干涉图样
(4)根据步骤(3)得到的干涉图样，利用波导不变量法得到干涉图样中的干涉条纹斜率，并根据这一斜率和海洋环境变化的有关参数的关系，对海洋环境变化进行监测。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710894274.3A CN108107437B (zh) | 2017-09-28 | 2017-09-28 | 一种利用简正波耦合干涉的海洋环境监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710894274.3A CN108107437B (zh) | 2017-09-28 | 2017-09-28 | 一种利用简正波耦合干涉的海洋环境监测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108107437A CN108107437A (zh) | 2018-06-01 |
| CN108107437B true CN108107437B (zh) | 2021-05-28 |
Family
ID=62206768
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710894274.3A Active CN108107437B (zh) | 2017-09-28 | 2017-09-28 | 一种利用简正波耦合干涉的海洋环境监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108107437B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108875234B (zh) * | 2018-06-28 | 2022-03-11 | 中国船舶科学研究中心（中国船舶重工集团公司第七0二研究所） | 应用于船舶三维声弹性分析的浅海声传播计算方法 |
| CN113191307B (zh) * | 2021-05-17 | 2023-09-29 | 江苏科技大学 | 一种基于贝叶斯估计的浅海波导简正波阶数判定方法 |
| CN113297801B (zh) * | 2021-06-15 | 2022-10-14 | 哈尔滨工程大学 | 一种基于steof-lstm的海洋环境要素预测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09274081A (ja) * | 1996-04-03 | 1997-10-21 | Tech Res & Dev Inst Of Japan Def Agency | 海底堆積層における音速測定方法 |
| JP2003161773A (ja) * | 2001-11-28 | 2003-06-06 | Mitsubishi Heavy Ind Ltd | 水中音響伝播特性模擬のための装置、方法及びプログラム |
| CN101604019A (zh) * | 2009-07-13 | 2009-12-16 | 中国船舶重工集团公司第七一五研究所 | 一种海洋环境与声场不确实性的表征和传递的快速计算方法 |
| CN102333052A (zh) * | 2011-05-26 | 2012-01-25 | 西北工业大学 | 一种适用于浅海低频条件的水声信号盲解卷方法 |
| CN106019288A (zh) * | 2016-06-17 | 2016-10-12 | 西北工业大学 | 基于简正波模态消频散变换的声源距离深度估计方法 |
-
2017
- 2017-09-28 CN CN201710894274.3A patent/CN108107437B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09274081A (ja) * | 1996-04-03 | 1997-10-21 | Tech Res & Dev Inst Of Japan Def Agency | 海底堆積層における音速測定方法 |
| JP2003161773A (ja) * | 2001-11-28 | 2003-06-06 | Mitsubishi Heavy Ind Ltd | 水中音響伝播特性模擬のための装置、方法及びプログラム |
| CN101604019A (zh) * | 2009-07-13 | 2009-12-16 | 中国船舶重工集团公司第七一五研究所 | 一种海洋环境与声场不确实性的表征和传递的快速计算方法 |
| CN102333052A (zh) * | 2011-05-26 | 2012-01-25 | 西北工业大学 | 一种适用于浅海低频条件的水声信号盲解卷方法 |
| CN106019288A (zh) * | 2016-06-17 | 2016-10-12 | 西北工业大学 | 基于简正波模态消频散变换的声源距离深度估计方法 |
Non-Patent Citations (2)
| Title |
|---|
| 一种基于模态匹配的浅海波导中宽带脉冲声源的被动测距方法;李晓曼 等;《物理学报》;20170418;第66卷(第9期);第1-16页 * |
| 浅海内波影响下的波导不变量变化特性分析;宋文华 等;《物理学报》;20141008;第63卷(第19期);第1-9页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108107437A (zh) | 2018-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110454687B (zh) | 一种基于改进vmd的管道多点泄漏定位方法 | |
| CN102262246B (zh) | 海上地震数据采集方法和装置 | |
| CN105308479B (zh) | 通过与偏移距相关的弹性fwi的多参数反演 | |
| CN108107437B (zh) | 一种利用简正波耦合干涉的海洋环境监测方法 | |
| CN103293515B (zh) | 一种舰船线谱噪声源纵向分布特性测量方法 | |
| JP6236282B2 (ja) | 異常検出装置、異常検出方法、及びコンピュータ読み取り可能な記憶媒体 | |
| US10877175B2 (en) | Seismic acquisition geometry full-waveform inversion | |
| CN104678384B (zh) | 一种波束域的声压差互相关谱分析水下目标速度估计方法 | |
| US10877174B2 (en) | Diving wave illumination using migration gathers | |
| Bonnel et al. | Waveguide mode amplitude estimation using warping and phase compensation | |
| Emmetière et al. | Understanding deep-water striation patterns and predicting the waveguide invariant as a distribution depending on range and depth | |
| CN109884592B (zh) | 一种面向低频高斯噪声源的声源定位仿真方法 | |
| US11455445B2 (en) | Method and system for non-intrusively inspecting a fluidic channel | |
| Le Courtois et al. | Autoregressive model for high-resolution wavenumber estimation in a shallow water environment using a broadband source | |
| Guarino et al. | Bottom attenuation coefficient inversion based on the modal phase difference between pressure and vertical velocity from a single vector sensor | |
| US8306762B2 (en) | Systems and methods for analysis of downhole data | |
| NL2022689B1 (en) | Method and system for non-intrusively determining cross-sectional variation for a fluidic channel | |
| Seo et al. | A novel sensing method of fault in moving machine | |
| KR101282692B1 (ko) | 충격음의 음장 표시 방법 | |
| Li et al. | Sparse Bayesian learning for horizontal wavenumber retrieval in underwater acoustical signal processing | |
| An et al. | Low resolution fourier synthesis modelling for underwater acoustic channel impulse response | |
| Lu et al. | Array shape estimation method based on extraction of single normal mode | |
| Wang et al. | An underwater acoustic channel simulator for UUV communication performance testing | |
| Zhao et al. | Single-hydrophone-based passive localization for vertically moving targets in shallow water | |
| Wu et al. | Warm-core eddy effects on sound propagation in the South China sea |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240111 Address after: Room 301, No. 727, Haijing Road, Lincheng Street, Zhoushan Marine Science City, Dinghai District, Zhoushan City, Zhejiang Province 316022 Patentee after: Haizhuoke (Zhejiang Free Trade Zone) Technology Co.,Ltd. Address before: 150001 Intellectual Property Office, Harbin Engineering University science and technology office, 145 Nantong Avenue, Nangang District, Harbin, Heilongjiang Patentee before: HARBIN ENGINEERING University |