CN108152800B - 一种适用于包含多杂质的河冰的河冰rcs的确定方法 - Google Patents
一种适用于包含多杂质的河冰的河冰rcs的确定方法 Download PDFInfo
- Publication number
- CN108152800B CN108152800B CN201711405750.7A CN201711405750A CN108152800B CN 108152800 B CN108152800 B CN 108152800B CN 201711405750 A CN201711405750 A CN 201711405750A CN 108152800 B CN108152800 B CN 108152800B
- Authority
- CN
- China
- Prior art keywords
- scattering
- rcs
- ice
- determining
- obtaining
- 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
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种适用于包含多杂质的河冰的河冰RCS的确定方法,包括:获取具有随机分布的球形泥沙散射粒子的VRT参数;获取具有随机分布的椭球形气泡散射粒子的VRT参数;通过算术平均方法获得不同河冰介质层的介电常数;确定不同介质交界面的透射率矩阵;将球形泥沙散射粒子的VRT参数、椭球形气泡散射粒子的VRT参数以及透射率矩阵和预设的辐射传输方程,获得不同散射情形下电磁波的散射强度;根据所述不同散射情形下电磁波的散射强度得到不同散射情形下河冰的RCS;根据所述不同散射情形下河冰的RCS确定总的RCS。本发明能够计算河冰上下层表面RCS、河冰内部不同介质交界面的RCS、不同介质层中散射粒子贡献的RCS以及总的RCS。本发明具有测量精度高的特点。
Description
技术领域
本发明涉及一种利用微波散射测量包含多杂质的河冰的河冰RCS的确定方法,属于微波测量领域。
背景技术
在黄河内蒙古河段由于纬度及气温的影响,河冰每年都会封冻,在流凌期可能还会发生流凌灾害。黄河的组成成分与世界其他河流还有很大不同,黄河的输沙量位居世界第一,在这种情况下,黄河冰的组成成分也会有一定的泥沙。在目前的研究中,多国的学者对河冰中气泡的电磁散射进行大量研究,认为气泡是造成河冰体散射的重要原因。但是对于黄河冰这种除了包含气泡还包含较多泥沙杂质的研究还是非常少。
发明内容
有鉴于上述技术问题,本发明提供了一种适用于包含多杂质的河冰的河冰RCS的确定方法,其中RCS(英文全称:Radar Cross Section,雷达散射截面积),其中河冰是指河水因热量变化产生的结冰现象。本发明实施例通过求解辐射传输方程得到电磁波的散射强度,进一步确定河冰的RCS。具体的实施过程包括:
S2、获取具有随机分布的球形泥沙散射粒子的VRT参数,具体包括散射振幅函数、散射吸收参数、相矩阵;
S3、获取具有随机分布的椭球形气泡散射粒子的VRT参数,具体包括散射振幅函数、散射吸收参数、相矩阵;
S4、通过算术平均方法获得不同河冰介质层的介电常数;
S5、确定不同介质交界面的透射率矩阵;
S8、将不同散射情形计算出的RCS相加算出总的RCS;
进一步地,步骤S2中,求解球形泥沙粒子VRT参数的具体方法为:
S21、根据瑞利近似ka<<1,a为散射体的半径,不考虑介质层中球形散射粒子的内部散射,以此求出散射电场,再根据散射电场与入射电场的关系求出散射振幅函数fpq;
S23、计算球形散射体散射粒子散射系数κs2、吸收系数κa1以及背景介质的吸收系数κa2;
进一步地,步骤S3中,求解椭球形气泡粒子VRT参数的具体方法为:
S31、计算过程与球形粒子相似,只是由于椭球形粒子不是非常对称计算复杂度要大一些,根据瑞利近似kc<<1,c为散射体的最大半长轴,不考虑介质层中球形散射粒子的内部散射,以此求出散射电场,再根据散射电场与入射电场的关系求出散射振幅函数fpq;
S33、计算椭球形散射体散射粒子散射系数κs1和消光系数κe1;
进一步地,步骤S5中,求解不同介质交界面的透射率矩阵具体步骤为:
S51、求解空气/纯冰+气泡表面的透射矩阵,假设空气中的介电常数与真空中的相同,纯冰+气泡的介电常数由S4给出;
S52、求解纯冰+气泡/纯冰+泥沙表面的透射矩阵,S4给出两介质层的介电常数;
S53、求解纯冰+泥沙/纯冰+气泡表面的透射矩阵,S4给出两介质层的介电常数;
S54、求解纯冰+气泡/空气表面的透射矩阵,纯冰+气泡的介电常数由S4给出,空气中的介电常数与真空中的相同;
进一步地,步骤S7中,求解不同散射情形下河冰的RCS具体方法为:
S71、利用积分方程模型计算低粗糙度(ks<3)空气/冰表面的RCS:σ1;
S72、计算纯冰+空气泡层中气泡贡献的RCS:σ2;
S73、计算具有低粗糙度的纯冰+空气泡/纯冰+泥沙两层介质分界面的RCS:σ3;
S74、计算纯冰+泥沙层中泥沙贡献的RCS:σ4;
S75、计算具有低粗糙度的纯冰+泥沙/水表面的RCS:σ5。
本发明能够计算河冰上下层表面RCS、河冰内部不同介质交界面的RCS、不同介质层中散射粒子贡献的RCS以及总的RCS。
本发明实施例利用微波散热测量的方式确定包含多杂质的河冰的河冰RCS,其具有测量精度高,且适用于多杂质的河冰。
附图说明
图1为河冰的分层模型图;
图2为本发明实施例中的河冰RCS的确定方法的原理流程图。
具体实施方式
为使本领域技术人员更好的理解本发明的技术方案,下面结合附图和具体实施方式对本发明作详细说明。
下面,结合图1和图2对本发明的具体实施过程进行更加详细地说明:本发明实施例中的河冰RCS的确定方法可以包括:
辐射传输理论主要研究电磁场在随机介质中多次散射、吸收和传播的过程。矢量辐射传输理论可以表述为一个强度为的电磁波在入射方向上传播,在这个方向上消耗掉的总能量等于沿着这个方向上散射粒子的散射、吸收和背景介质的吸收以及来自散射方向上的散射之和。矢量辐射传输方程借助斯托克斯矢量考虑电磁场强度叠加忽略相位的影响。在主动遥感中矢量辐射传输方程可以写成如下形式:
κe为消光矩阵,表示电磁波在传播过程中散射粒子散射的能量κs、散射粒子吸收能量κas以及背景介质吸收能量κba之和,即κe=κs+κas+κba。为相矩阵,本身与相位无关,它表示入射辐射强度和散射辐射强度的耦合,即多次散射。
空间中的电磁波可以用如下形式来表示
角括号表示综系平均;v和h表示电磁波的垂直极化分量和水平极化分量;在计算相矩阵时,需要用到散射振幅函数fpq,p和q表示v或者h。
S2、获取具有随机分布的球形泥沙散射粒子的VRT参数,具体包括散射振幅函数、散射吸收参数、相矩阵;
S21、散射振幅函数描述了散射电场与入射电场的关系。对于不同形状的散射粒子,散射振幅函数是不同的。对于黄河冰来说其主要的散射粒子为气泡和泥沙,气泡可近似为随机分布的椭球状粒子,而泥沙可近似为随机分布的球形粒子。泥沙的尺寸相对较小,在L波段和C波段下满足瑞利近似。瑞利近似下球形散射粒子的散射场为
根据散射场与入射场的对应关系可以求出瑞利近似下的球形粒子散射振幅函数表示为
a为散射粒子的半径,s和i分别表示反射和入射。
S22、球形散射粒子的相矩阵为
相矩阵中的元素可以按照定义将散射振幅函数带入到相矩阵中求出。
S23、由于瑞利近似是一种低频近似,它的散射振幅函数达不到足够的精确,因此不能用光学定理来计算消光κe。散射粒子散射系数κs2、吸收系数κa1以及背景介质的吸收系数κa2可通过如下公式计算
κa2=2k”(1-fs)
式中k”为波数的虚部,fs为散射粒子的占空比,k'为波数的实部,n0为单位体积内散射粒子的个数,εs为散射粒子的介电常数。散射粒子的吸收系数与背景介质的吸收系数之和为总的吸收系数κa,κa=κa1+κa2。
S3、获取具有随机分布的椭球形气泡散射粒子的VRT参数,具体包括散射振幅函数、散射吸收参数、相矩阵;
S31、河冰中的气泡可以近似为椭球状散射粒子,在L波段和C波段下,气泡的几何尺寸相比波长是比较小的,所以仍然可以采取瑞利近似。瑞利近似下椭球粒子的散射场形式为
在瑞利近似下,椭球粒子的散射振幅函数为:
其中
式中Rs=[(s+a2)(s+b2)(s+c2)]0.5,a、b、c分别为椭球气泡的三个半长轴,可以近似气泡的两个半长轴a和b的长度是相等的且半长轴c的长度大于a和b,这样得到
g1=0.5(1-g3)
式中的方向矢量内乘可表示为
式中的β、γ为椭球粒子的坐标系与主坐标系之间的欧拉旋转角。
S32、椭球状粒子的相矩阵为
S33、椭球粒子的消光矩阵可以由光学定理给出
椭球粒子的散射系数为
下标中的h和v分别表示水平极化下和垂直极化下散射粒子的散射系数。
S4、通过算术平均方法确定不同河冰介质层的介电常数;
对于纯冰+气泡和纯冰+泥沙介质层,由于散射粒子的加入使得冰的介电常数发生变化。由于气泡和泥沙的介电常数与纯冰的介电常数相差较小,只相当于在纯冰中加入了一个小的扰动,所以可用简单的算术平均方法计算介质层的介电常数。
S5、获得不同介质交界面的透射率矩阵;
电磁波在传播的过程中会遇到不同介质的表面,这时部分电磁波能量被分界面反射,形成反射波;另外一部分将透过分界面继续传播,形成透射波。在VRT中反射与透射能量的多少可以由反射率矩阵和透射率矩阵表达。透射系数的形式为:
式中m、n代表不同介质,Tmn(θm)表示从介质m到介质n的透射率矩阵,介电常数的上标'表示介电常数的实部。θm和θn可以通过Snell定理求出。
S6、将S2与S3获得的VRT参数以及S5中的透射率矩阵带入到辐射传输方程中,确定不同散射情形下电磁波的散射强度Iq s;
将辐射传输方程应用到如图1所示河冰模型当中,在此我们只考虑单次散射情况,即求出辐射传输方程的一阶解。河冰的全部散射强度可以近似为以下几个部分:(1)空气/冰表面的散射强度Is1;(2)纯冰+空气泡层中气泡的体散射强度Is2;(3)纯冰+空气泡/纯冰+泥沙两层介质分界面的散射强度Is3;(4)纯冰+泥沙层中泥沙的体散射强度Is4;(5)纯冰+泥沙/水表面的面散射强度Is5。总的散射强度I为以上五种散射强度之和。
S71、利用积分方程模型确定低粗糙度(ks<3)空气/冰表面的RCS:σ1;
对于粗糙度不太强的随机粗糙面(一般为ks<3),可以使用积分方程模型计算单次散射的RCS,计算公式为:
其中
fpp=2Rp/cosθ0p为h或v
Rh=(η0cosθ0-η1cosθ1)/(η0cosθ0+η1cosθ1)
Rv=(η1cosθ0-η0cosθ1)/(η1cosθ0+η0cosθ1)
式中εr为空气介质与纯冰+气泡层的介电常数比。
Rv与Rh为菲涅儿反射系数。
S72、确定纯冰+空气泡层中气泡贡献的RCS:σ2;
电磁波从空气层进入冰层发生折射,在冰层中发生衰减,当电磁波遇到散射体时能量散射到各个方向,电磁波的衰减用消光矩阵表示,则由气泡贡献的RCS可以表示为
其中
式中L01和L10为介质0和介质1由于粗糙度引起的损耗,P1为纯冰+气泡层的相矩阵。
S73、确定具有低粗糙度的纯冰+空气泡/纯冰+泥沙两层介质分界面的RCS:σ3;
电磁波在经过空气/冰表面的面散射损耗以及纯冰+气泡内部的吸收损耗,电磁波照射到纯冰+气泡/纯冰+泥沙表面后经过散射。此时的RCS可以表示为
在计算σ1时要将里面的参数换成纯冰+气泡层中的参数。
S74、确定纯冰+泥沙层中泥沙贡献的RCS:σ4;
电磁波经过空气层、纯冰+气泡层,在纯冰+泥沙层中入射到泥沙上造成散射。泥沙粒子的RCS为
下标2代表第二层介质(纯冰+泥沙)的参数,下标s代表散射,式中
S75、确定具有低粗糙度的纯冰+泥沙/水表面的RCS:σ5;
纯冰+泥沙/水表面可以近似为具有较小的粗糙度,电磁波在纯冰+泥沙/水表面贡献的RCS为
在计算σ1时要将里面的参数换成纯冰+泥沙层中的参数。
S8、将不同散射情形计算出的RCS相加算出总的RCS;
经过以上计算可以得到总的RCS,即
σ=σ1+σ2+σ3+σ4+σ5
综上所述,本发明能够确定河冰上下层表面RCS、河冰内部不同介质交界面的RCS、不同介质层中散射粒子贡献的RCS以及总的RCS。
另外,本发明实施例中的河冰RCS确定方法可以应用在RCS确定装置中,并通过其中的数据处理模块执行上述河冰RCS确定方法的操作,从而执行包括多杂质的河冰的河冰RCS的确定过程,具有测量精度高的特点。
以上实施例仅为本发明的示例性实施例,不用于限制本发明,本发明的保护范围由权利要求书限定。本领域技术人员可以在本发明的实质和保护范围内,对本发明做出各种修改或等同替换,这种修改或等同替换也应视为落在本发明的保护范围内。
Claims (8)
1.一种适用于包含多杂质的河冰的河冰RCS的确定方法,其特征在于,包括:
步骤1,获取具有随机分布的球形泥沙散射粒子的VRT参数;
步骤2,获取具有随机分布的椭球形气泡散射粒子的VRT参数;
步骤3,通过算术平均方法获得不同河冰介质层的介电常数;
步骤4,确定不同介质交界面的透射率矩阵;
步骤5,基于球形泥沙散射粒子的VRT参数、椭球形气泡散射粒子的VRT参数以及透射率矩阵和预设的辐射传输方程,获得不同散射情形下电磁波的散射强度;
步骤6,根据所述不同散射情形下电磁波的散射强度得到不同散射情形下河冰的RCS;
步骤7,根据所述不同散射情形下河冰的RCS确定总的RCS。
2.根据权利要求1所述的方法,其特征在于,步骤1包括:
步骤11,根据散射电场与入射电场的关系求出瑞利近似下的散射振幅函数;
步骤12,根据入射方向的电磁波强度与来自散射方向的电磁波强度的耦合关系求出相矩阵,并根据散射振幅函数求出相矩阵中的元素;
步骤13,获得球形散射体散射粒子散射系数、吸收系数以及背景介质的吸收系数。
3.根据权利要求1所述的方法,其特征在于,步骤2包括:
步骤21,根据散射电场与入射电场的关系求出瑞利近似下的散射振幅函数;
步骤22,根据入射方向的电磁波强度与来自散射方向的电磁波强度的耦合关系求出相矩阵,并根据散射振幅函数求出相矩阵中的元素;
步骤23,获得椭球形散射体散射粒子散射系数和消光系数。
4.根据权利要求1所述的方法,其特征在于,步骤4包括:
步骤41,确定空气/纯冰+气泡表面的透射矩阵;
步骤42,确定纯冰+气泡/纯冰+泥沙表面的透射矩阵;
步骤43,确定纯冰+泥沙/纯冰+气泡表面的透射矩阵;
步骤44,确定纯冰+气泡/空气表面的透射矩阵;
其中,所述步骤41-44中使用到的各种介电常数由步骤3确定。
5.根据权利要求1所述的方法,其特征在于,步骤6包括:
步骤61,利用积分方程模型获得低粗糙度空气/冰表面的RCS;
步骤62,确定纯冰+空气泡层中气泡贡献的RCS;
步骤63,确定具有低粗糙度的纯冰+空气泡/纯冰+泥沙两层介质分界面的RCS;
步骤64,确定纯冰+泥沙层中泥沙贡献的RCS;
步骤65,确定具有低粗糙度的纯冰+泥沙/水表面的RCS。
6.根据权利要求5所述的方法,其特征在于,步骤7包括:将步骤61-65中得到的各RCS求和得到总的RCS。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711405750.7A CN108152800B (zh) | 2017-12-22 | 2017-12-22 | 一种适用于包含多杂质的河冰的河冰rcs的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711405750.7A CN108152800B (zh) | 2017-12-22 | 2017-12-22 | 一种适用于包含多杂质的河冰的河冰rcs的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108152800A CN108152800A (zh) | 2018-06-12 |
CN108152800B true CN108152800B (zh) | 2021-11-12 |
Family
ID=62465286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711405750.7A Active CN108152800B (zh) | 2017-12-22 | 2017-12-22 | 一种适用于包含多杂质的河冰的河冰rcs的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108152800B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112307679B (zh) * | 2020-11-23 | 2022-11-04 | 内蒙古工业大学 | 一种构建河冰厚度反演微波散射模型的方法及装置 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3798247B2 (ja) * | 1998-08-07 | 2006-07-19 | ザ・ジョーンズ・ホプキンス・ユニバーシティ | 並列ドップラ処理を行うレーダ式の氷測深器 |
CN101493522B (zh) * | 2009-02-27 | 2012-02-01 | 武汉理工大学 | 海面反常波检测方法及装置 |
US8947295B2 (en) * | 2012-06-27 | 2015-02-03 | The United States Of America As Represented By The Secretary Of The Air Force | Low clutter method for bistatic RCS measurements |
KR101460591B1 (ko) * | 2013-05-03 | 2014-11-12 | 한남대학교 산학협력단 | 특성 벡터 융합을 이용한 레이더 표적 인식 방법 및 장치 |
CN104569926B (zh) * | 2015-01-06 | 2017-07-07 | 北京环境特性研究所 | 一种地物波谱rcs测量的控制方法及系统 |
CN105486664B (zh) * | 2015-12-31 | 2018-04-20 | 浙江大学 | 探测海洋浮游植物生物量和poc的激光雷达装置及方法 |
CN105928905B (zh) * | 2016-06-14 | 2018-11-09 | 国家海洋局第二海洋研究所 | 颗粒物后向散射系数偏振敏感性测量装置 |
CN106326579B (zh) * | 2016-08-27 | 2019-03-15 | 安徽大学 | 一种基于电磁散射技术的地海环境可视化生成方法 |
CN106569191A (zh) * | 2016-10-17 | 2017-04-19 | 西北工业大学 | 一种利用高分辨率成像获取目标rcs的方法 |
-
2017
- 2017-12-22 CN CN201711405750.7A patent/CN108152800B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108152800A (zh) | 2018-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kunkel et al. | Monte Carlo analysis of multiply scattered lidar returns | |
CN110095784B (zh) | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 | |
Guo et al. | Applications of dense media radiative transfer theory for passive microwave remote sensing of foam covered ocean | |
CN102460113A (zh) | 监测介质的颗粒尺寸分布的超声法方法 | |
CN110568449B (zh) | 一种风生粗糙海面激光反射、透射矩阵计算方法 | |
Niu et al. | Analytic solution of long wave propagation over a submerged hump | |
CN108152800B (zh) | 一种适用于包含多杂质的河冰的河冰rcs的确定方法 | |
Sundarabalan et al. | Modeling the underwater light field fluctuations in coastal oceanic waters: Validation with experimental data | |
CN110376572A (zh) | 一种机载激光雷达水下探测深度模拟方法和装置 | |
Bao et al. | Study of a laser echo in an inhomogeneous dust environment with a continuous field Monte Carlo radiative transfer model | |
Rosenkranz | Radiative transfer solution using initial values in a scattering and absorbing atmosphere with surface reflection | |
Zang et al. | Analysis of computation model of particle deposition on transmittance for photovoltaic panels | |
Kajino et al. | NHM-Chem, the Japan MeteorologicalAgency's regional meteorology–chemistry model (v1. 0): model description and aerosol representations | |
Tsukui et al. | Effective roughness on the sea surface for determining variability characteristics of reflected sound waves | |
Kuo | Acoustic wave scattering from two solid boundaries at the ocean bottom: reflection loss | |
Mon | Backward and forward scattering of microwaves by ice particles: A review | |
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 | |
Indira Devi et al. | A study of VLF wave propagation characteristics in the earth-ionosphere waveguide | |
Benjamin et al. | Extinction of Millimeter wave on Two Dimensional Slices of Foam-Covered Sea-surface | |
Kim | Influence of surf noise on oceanic ambient noise in coastal waters with sandy seabed | |
Zhou et al. | Polarimetric passive microwave remote sensing of wind vectors with foam-covered rough ocean surfaces | |
CN114791992B (zh) | 一种基于射线理论的深海目标前向散射声场计算方法 | |
Fois et al. | Future ocean scatterometry at very strong winds | |
Li | Research on the electromagnetic scattering from foam sea based on small slope approximation | |
Shaker et al. | Scattering of infrared by prismatic cloud ice crystals |
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 |