CN108680915A - 雷达波束下含破碎浪及泡沫海面散射分区并行计算方法 - Google Patents

雷达波束下含破碎浪及泡沫海面散射分区并行计算方法 Download PDF

Info

Publication number
CN108680915A
CN108680915A CN201810116607.4A CN201810116607A CN108680915A CN 108680915 A CN108680915 A CN 108680915A CN 201810116607 A CN201810116607 A CN 201810116607A CN 108680915 A CN108680915 A CN 108680915A
Authority
CN
China
Prior art keywords
foam
wave
sea
broken
scattering
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
Application number
CN201810116607.4A
Other languages
English (en)
Other versions
CN108680915B (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN201810116607.4A priority Critical patent/CN108680915B/zh
Publication of CN108680915A publication Critical patent/CN108680915A/zh
Application granted granted Critical
Publication of CN108680915B publication Critical patent/CN108680915B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/006—Theoretical aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

一种雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，包括步骤：建立三维含破碎浪及泡沫电大尺寸非线性海面复合模型；对雷达主波束进行一级子波束分解分区，对子波束进行角谱展开；对不含泡沫海面二级分区，建立确定性面散射系数模型；对含静态泡沫海面进行二级分区，建立确定性体面散射系数模型；对含动态泡沫的破碎浪进行二级分区，建立确定性破碎浪及泡沫体面散射系数模型；利用CUDA对两重分区分别并行计算，并利用散射系数的非相干叠加得到总散射系数。优点是，雷达波束照射下利用波束分解对含破碎浪及泡沫海面进行分区，利用角谱展开结合相应的高频算法提供一种准确有效的计算含破碎浪及泡沫的电大尺寸非线性海面散射特性的方法。

Description

雷达波束下含破碎浪及泡沫海面散射分区并行计算方法
技术领域
本发明涉及一种雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，利用波束分解和角谱展开可针对含有浪涌、泡沫、破碎浪等复杂海态、海况下的海面结构进行分区，并结合不同的算法并行处理，有效提高计算效率。
背景技术
小擦地角下海杂波特性研究一直都是军事探测和民用遥感领域的热点问题，而雷达波束以小擦地角入射时会在海面形成巨大的照射区域，使得物理分析和数学求解变得异常复杂。巨大照射面积下往往会包含多种海态和海况，可能既有普通的浪涌又有卷曲破碎浪和泡沫，不同海表面对雷达后向回波的影响不同，很难用一种方法同时高效地求解各种海表面形态。
为了充分描述海面不同区域尤其是破碎区域的散射特征，利用波束分解结合角谱展开方法模拟真实的雷达照射波束，波束分解将大尺度波束分解成不同幅度的小尺度波束，各小尺度波束通过相干叠加得到和大尺多波束相同的效果。用角谱展开方法将不同照射区域的子波束展开为不同平面波的叠加。对不含破碎浪的普通海面利用面元双尺度法进行求解；利用IPO-EEC方法求解破碎浪多次散射和绕射的影响；针对海面上方出现泡沫覆盖区域，利用矢量辐射传输理论，求解泡沫层的影响。
目前还没有一种能够针对不同区域采用不同方法，并考虑不同因素的影响，利用并行技术加速进行计算的方法，现存方法存在计算效率低，计算结果不准确的问题。
发明内容
本发明所要解决的技术问题是：提供一种雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，建立含破碎浪及泡沫的非线性海面复合模型，并利用波束分解对含破碎浪及泡沫海面进行分区，利用角谱展开结合不同算法分别对不同子区域下的面元并行计算，利用散射系数分非相干叠加得到含破碎浪及泡沫海面总散射系数。
为解决上述技术问题，本发明的技术方案是：雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，包括以下步骤：
(1)建立三维含破碎浪及泡沫电大尺寸非线性海面复合模型。
(2)对雷达主波束进行一级子波束分解分区，对子波束进行角谱展开。
(3)对不含泡沫海面二级分区，建立确定性面散射系数模型；对含静态泡沫海面进行二级分区，建立确定性体面散射系数模型。
(4)对含动态泡沫的破碎浪进行二级分区，建立确定性破碎浪及泡沫体面散射系数模型。
(5)利用CUDA(统一计算设备架构)对两重分区分别并行计算，并利用散射系数的非相干叠加得到总散射系数。
所述的步骤1)中建立三维含破碎浪及泡沫电大尺寸非线性海面复合模型，利用动态泡沫和静态泡沫的覆盖率结合斜率判据分别得到动态泡沫及破碎浪的分布位置以及静态泡沫在海面的分布位置；
所述的步骤2)对雷达主波束进行波束分解，使之分解为一系列具有不同振幅及照射区域的子波束，根据子波束照射区域对含破碎浪及泡沫电大尺寸非线性海面复合模型进行一级分区，分为含动态泡沫的破碎浪部分，含静态泡沫的海面部分以及不含泡沫的普通海面部分，基于GPU及CUDA，将三部分分别放置不同的grid中，使得不同区域可利用不同算法并行计算。同时利用傍轴近似，将子波束展开为一系列不同振幅的平面波，作为以下高频算法的入射激励。
所述的步骤3)中，对不含泡沫的电大尺寸海面进行二级分区，利用面元双尺度法求解确定性海面面元散射系数；对含静态泡沫的电大尺寸海面进行二级分区，利用矢量辐射传输方法结合面元双尺度法，求解确定性体面散射系数。基于CUDA，将每一面元分别放置不同的thread中并行运算。
所述的步骤4)中，对含泡沫破碎浪进行二级分区，对平面波照射区域利用基于面元的迭代物理光学法(IPO)及等效边缘电磁流法(EEC)求解破碎浪散射系数，其中利用IPO方法考虑破碎浪与海面耦合作用，利用EEC考虑破碎浪劈绕射作用；结合利用矢量辐射传输理论，求解确定性破碎浪上动态泡沫层的影响。基于CUDA，将每一面元分别放置不同的thread中并行运算。
所述的步骤5)中，基于步骤(1)的含破碎浪及泡沫复合海面模拟及步骤(2)(3)、(4)分别利用IPO-EEC法计算破碎浪散射系数和双尺度面元法计算得到的普通海面散射系数并考虑动态泡沫和静态泡沫对其的影响进行非相干叠加，利用CUDA对两级分区分别并行计算，可得到含破碎浪及泡沫海面散射系数。
采用了上述技术方案，本发明的有益效果在于，在雷达波束照射下，基于波束分解对含破碎浪海面进行分区，利用角谱展开结合高频算法含破碎浪及泡沫非线性电大尺寸海面分区并行计算提供一种准确有效的计算含破碎浪及泡沫电大尺寸复合海面散射系数并行计算方法。
附图说明
图1是本发明含破碎浪电大尺寸海面并行计算流程图；
图2是本发明中风速为20m/s时，破碎浪及泡沫在电大尺寸海面上的分布图。
具体实施方式
下面结合实例对本发明进一步说明。
参照图1和图2，本发明的具体实现步骤如下：
步骤1，建立多尺度非线性海面模型,面元个数为M×N，得到海面高度起伏沿风向的斜率，从而得到斜率沿风向概率密度分布函数。根据下式得到动态泡沫及静态泡沫的覆盖率
其中，对于动态泡沫有a＝5，对于静态泡沫有a＝0.8，Λ(c,u10)为单位面积单位速度间隔内破碎波的平均长度，g为重力加速度，u10为海面上方10m处的风速，c表示波速。
根据下式，可以得到不同风速下的动态泡沫及静态泡沫的覆盖厚度：
其中，
对二维破碎浪x轴沿风向进行旋转，并沿风向垂直方向进行三维扩展，得到三维含风向破碎浪。根据动态泡沫覆盖率，求解沿风向斜率最大的面元个数为NC＝FC(U10)×M×N,并求出这NC面元中的最小斜率，标记为c1，若满足则用破碎浪替换该海面面元，从而得到含破碎浪海面结构建模，并在破碎浪处覆盖满足动态泡沫厚度的泡沫层。根据静态泡沫层覆盖率，求解沿风向斜率次大的面元个数为NS＝FS(U10)×M×N，并求出这NS个面元中的最小斜率，标记为c2，若满足则将该面元覆盖满足静态泡沫厚度的泡沫层，如图2所示为风速20m/s时的海面破碎浪及泡沫分布，海面尺寸为192m×192m，灰度图中黑色背景为普通海面，白色区域为破碎浪及动态泡沫分布位置，灰色区域为静态泡沫分布位置。
步骤2，将主波束用子波束展开为如下形式：
其中为主波束电场，表示子波束电场，Wmn为子波束权重因子。ΔX与ΔY分别为子波束沿x轴和y轴的偏移量。将整个海面分为含动态泡沫的破碎浪部分，含静态泡沫的海面部分以及不含泡沫的普通海面部分，结合子波束的照射区域分别计算。对子波束进一步利用角谱展开为平面波形式：
其中wi(kyi,kxi)为展开的各个平面波复振幅。
步骤3，利用将面元双尺度法计算一级分区确定性面元散射场
其中，为海面重力波的镜像散射场，为海面张力波的Bragg散射场，下标mn表示第mn个面元，从而得到确定性面元的散射系数对含静态泡沫区域利用矢量辐射传输理论可以得到单个面元含静态泡沫覆盖的破碎浪的散射系数：
其中
其中，θi为单个面元的局部入射角，ke与ks分别为单个粒子的消光系数及散射系数，Rmn-hh与Rmn-vv分别为经过倾斜面元调制的水平和垂直极化菲涅尔反射系数，-2kedssecθi为经过厚度为ds的静态泡沫粒子层的散射和吸收的衰减因子。
步骤4，利用迭代物理光学法求解破碎浪劈尖处组成的腔体散射及与之相邻海面部分的多次散射。利用阻抗劈绕射系数得到阻抗劈等效边缘电磁流，用电磁流辐射积分可得劈尖绕射场。利用场的矢量叠加性可得到单个破碎浪的散射场：
其中，为不同破碎浪破碎波腔体多次散射以及与之相邻海面耦合散射的散射场，为不同区域破碎浪劈尖绕射场。从而得到确定性破碎浪的散射系数利用矢量辐射传输理论，可以得到单个破碎浪含动态泡沫覆盖的破碎浪的散射系数：
其中
其中，σbw为单个破碎波的散射系数，-2kedcsecθi为经过厚度为dc的泡沫粒子层的散射和吸收的衰减因子。
步骤5，利用IPO-EEC法计算破碎浪散射系数和双尺度面元法计算得到的普通海面散射系数并考虑动态泡沫和静态泡沫对其的影响进行非相干叠加如下：
其中，不含破碎浪及泡沫的海面散射系数，为考虑了动态泡沫影响的破碎波散射系数，为考虑了静态泡沫影响的海面散射系数，A为雷达照射区域面积。利用CUDA对两级分区分别并行计算，可得到含破碎浪及泡沫海面散射系数。
采用了上述技术方案，本发明的有益效果在于，基于波束分解对含破碎浪海面进行分区，利用角谱展开结合高频算法含破碎浪电大尺寸海面分区并行计算提供一种准确有效的计算含破碎浪及泡沫电大尺寸复合海面散射系数并行计算方法。本发明不局限于上述具体的实施方式，本领域的普通技术人员从上述构思出发，不经过创造性的劳动，所作出的种种变换，均落在本发明的保护范围之内。

Claims (6)

1.一种雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，其特征在于，包括以下步骤：
(1)建立三维含破碎浪及泡沫电大尺寸非线性海面复合模型；
(2)对雷达主波束进行一级子波束分解分区，对子波束进行角谱展开；
(3)对不含泡沫海面二级分区，建立确定性面散射系数模型；对含静态泡沫海面进行二级分区，建立确定性体面散射系数模型；
(4)对含动态泡沫的破碎浪进行二级分区，建立确定性破碎浪及泡沫体面散射系数模型；
(5)利用CUDA对两重分区分别并行计算，并利用散射系数的非相干叠加得到总散射系数。
2.根据权利要求1所述的雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，其特征在于，步骤1)中的含破碎浪及泡沫非线性海面几何结构建模，利用动态泡沫和静态泡沫的覆盖率结合斜率判据分别得到动态泡沫及破碎浪的分布位置以及静态泡沫在海面的分布位置。
3.根据权利要求1所述的雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，其特征在于，步骤(2)中将雷达主波束进行波束分解,使之分解为一系列具有不同振幅及照射区域的子波束，根据子波束照射区域对含破碎浪及泡沫电大尺寸非线性海面复合模型进行一级分区，分为含动态泡沫的破碎浪部分，含静态泡沫的海面部分以及不含泡沫的普通海面部分，基于GPU(图形处理单元)及CUDA，将三部分分别放置不同的grid(线程网)中，使得不同区域可利用不同算法并行计算；同时利用傍轴近似，将子波束展开为一系列不同振幅的平面波，作为以下高频算法的入射激励。
4.根据权利要求1所述的雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，其特征在于，步骤(3)中对不含泡沫的电大尺寸海面进行二级分区，利用面元双尺度法求解确定性海面面元散射系数；对含静态泡沫的电大尺寸海面进行二级分区，利用矢量辐射传输方法结合面元双尺度法，求解确定性体面散射系数；基于CUDA，将每一面元分别放置不同的thread(线程)中并行运算。
5.根据权利要求1所述的雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，其特征在于，步骤(4)中对含泡沫破碎浪进行二级分区，对平面波照射区域利用基于面元的迭代物理光学法(IPO)及等效边缘电磁流法(EEC)求解破碎浪散射系数，其中利用IPO方法考虑破碎浪与海面耦合作用，利用EEC考虑破碎浪劈绕射作用；结合利用矢量辐射传输理论，求解确定性破碎浪上动态泡沫层的影响；基于CUDA，将每一面元分别放置不同的thread中并行运算。
6.根据权利要求1所述的雷达波束下含破碎浪及泡沫海面散射分区并行计算方法，其特征在于，步骤(2)(3)、(4)、(5)中利用CUDA对两级分区分别并行计算，利用散射系数的非相干叠加得到各区域散射系数分布以及总散射系数，提高计算效率。
CN201810116607.4A 2018-02-06 2018-02-06 雷达波束下含破碎浪及泡沫海面散射分区并行计算方法 Active CN108680915B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810116607.4A CN108680915B (zh) 2018-02-06 2018-02-06 雷达波束下含破碎浪及泡沫海面散射分区并行计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810116607.4A CN108680915B (zh) 2018-02-06 2018-02-06 雷达波束下含破碎浪及泡沫海面散射分区并行计算方法

Publications (2)

Publication Number Publication Date
CN108680915A true CN108680915A (zh) 2018-10-19
CN108680915B CN108680915B (zh) 2021-08-10

Family

ID=63800132

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810116607.4A Active CN108680915B (zh) 2018-02-06 2018-02-06 雷达波束下含破碎浪及泡沫海面散射分区并行计算方法

Country Status (1)

Country Link
CN (1) CN108680915B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112130122A (zh) * 2020-09-01 2020-12-25 武汉大学 一种空基高频雷达海面散射系数的估计方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040143395A1 (en) * 2001-11-07 2004-07-22 Takeshi Matsuoka Method of observing sea ice
CN101493522A (zh) * 2009-02-27 2009-07-29 武汉理工大学 海面反常波检测方法及装置
CN102306217A (zh) * 2011-08-12 2012-01-04 北京理工大学 基于非线性一维海面分形模型的电磁散射系数估计方法
CN104992064A (zh) * 2015-07-09 2015-10-21 西安电子科技大学 多尺度电大区域粗糙海面后向雷达回波分区并行计算方法
CN105701284A (zh) * 2016-01-11 2016-06-22 西安电子科技大学 时变多尺度电大区域海面电磁散射矢量场的并行计算方法
CN106225693A (zh) * 2016-08-29 2016-12-14 中国科学院遥感与数字地球研究所 一种细粒子气溶胶光学厚度和气溶胶类型同时反演方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040143395A1 (en) * 2001-11-07 2004-07-22 Takeshi Matsuoka Method of observing sea ice
CN101493522A (zh) * 2009-02-27 2009-07-29 武汉理工大学 海面反常波检测方法及装置
CN102306217A (zh) * 2011-08-12 2012-01-04 北京理工大学 基于非线性一维海面分形模型的电磁散射系数估计方法
CN104992064A (zh) * 2015-07-09 2015-10-21 西安电子科技大学 多尺度电大区域粗糙海面后向雷达回波分区并行计算方法
CN105701284A (zh) * 2016-01-11 2016-06-22 西安电子科技大学 时变多尺度电大区域海面电磁散射矢量场的并行计算方法
CN106225693A (zh) * 2016-08-29 2016-12-14 中国科学院遥感与数字地球研究所 一种细粒子气溶胶光学厚度和气溶胶类型同时反演方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YA-QIU JIN ET AL.: "Observation of scattering and emission from fluctuating, rough surface of water basin by using a system of scatterometer-radiometer at X band", 《IEEE 1999 INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》 *
李清亮 等: "《雷达地海杂波测量与建模》", 31 December 2017 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112130122A (zh) * 2020-09-01 2020-12-25 武汉大学 一种空基高频雷达海面散射系数的估计方法
CN112130122B (zh) * 2020-09-01 2022-09-13 武汉大学 一种空基高频雷达海面散射系数的估计方法

Also Published As

Publication number Publication date
CN108680915B (zh) 2021-08-10

Similar Documents

Publication Publication Date Title
CN104992064B (zh) 多尺度电大区域粗糙海面后向雷达回波分区并行计算方法
CN110095784A (zh) 一种复杂环境影响下的海洋-低层大气激光传输建模方法
CN104656078A (zh) 一种基于射线追踪的海面舰船的方法
Liang et al. The fast EPILE combined with FBM for electromagnetic scattering from dielectric targets above and below the dielectric rough surface
CN115267720A (zh) 海上舰船目标复合电磁散射rcs计算方法
CN108680915A (zh) 雷达波束下含破碎浪及泡沫海面散射分区并行计算方法
Lu et al. Estimating ice particle scattering properties using a modified Rayleigh‐Gans approximation
Zhao et al. Electromagnetic scattering of coated objects over sea surface based on SBR-SDFSM
Isleifson et al. A Monte Carlo method for simulating scattering from sea ice using FVTD
Bian et al. Efficient computation of broadband noise propagation using Gaussian beam tracing method
Keiffer et al. A time domain rough surface scattering model based on wedge diffraction: Application to low-frequency backscattering from two-dimensional sea surfaces
Zhang et al. Simulation of full-polarization electromagnetic backscattering characteristics of large number of high-density chaff clouds
CN104765028B (zh) 一种高斯随机起伏海底界面混响信号仿真方法
US6175815B1 (en) Storage reduction method for fast multipole field calculations
Benjamin et al. Millimeter-wave propagation and attenuation in closed packed sea-foam layer and complex dielectric constant of sea-foam using split-step fourier transform
Fuerkaiti et al. Aircraft community noise prediction in 3D environment using Gaussian beam tracing
Burrage et al. Predicting radiometric effects of a Rough Sea surface, whitecaps, foam, and spray using SURFER 2D
Daniel et al. Target observability improvement in multi-static maritime forward scatter radar
Bennani et al. Bistatic Radar Cross Section of an complex target on sea surface
Zhang et al. A volume-surface composite scattering model for nonlinear ocean surface with breaking waves and foam layers under high wind conditions
CN113126053B (zh) 用于评估舰船雷达散射截面的方法及装置
Xu et al. Multi-beam Line Measurement Problem Based on Geometric Modeling and Genetic Algorithm
RU80020U1 (ru) Установка для моделирования отражения радиоволн от морской поверхности
Zurk Electromagnetic wave propagation and scattering in dense, discrete random media with application to remote sensing of snow
Rino et al. Low‐frequency acoustic scatter from subsurface bubble clouds

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
CB03 Change of inventor or designer information

Inventor after: Zhang Xiaoxiao

Inventor after: Wu Zhensen

Inventor after: Su Xiang

Inventor before: Zhang Xiaoxiao

CB03 Change of inventor or designer information