本发明属于海洋物理学领域，具体涉及一种获取高真实感海浪泡沫光散射系数的方法。
背景技术：
随着多维光电探测手段的提高，实时获取与处理海面光电特性信息的需求不断增长。海浪泡沫作为一种常见的海面自然现象，由于海浪泡沫的存在会严重影响海面的光学特性，因此，探寻其光的散射特性是一项重要的基础性研究。在一定风速的影响下，破碎的海浪与空气进行混合，在海面上会产生一部分泡沫覆盖层。根据观测实验，当风速大于7米/秒时，海洋表面会产生白水，称为海浪泡沫，又称泡沫白帽，且海浪泡沫量会随着风速的增加而增加，当风速超过25米/秒时，约三分之一的海面将被泡沫覆盖。光与海浪泡沫作用时会发生多重反射与折射，且海浪泡沫的不同尺寸、浓度和盐度等属性对光的散射系数影响不同，从而会改变海面的光谱散射特性。
现阶段对海浪泡沫的研究主要集中在粗糙海面上电磁波的双站散射和热发射率的理论计算上，通常需要经过复杂的数值运算，如求解一定边界条件下的vrt(vectorradiativetransfer，矢量辐射传输)方程、波谱的傅里叶变换以及海面几何起伏的动态变换n-s(navier-stokes)方程等过程。
由于受海浪泡沫对光散射的影响，无论是主动式还是被动式的海上目标探测，海浪泡沫都会严重影响探测信号的能量分布。现有的海浪泡沫散射系数计算方法大都是将海浪泡沫假设成在海面上平铺的一层均匀电介质层，而海浪泡沫并不是均匀的，与实际情况不相符，这种海浪泡沫的假设对获得的海浪泡沫散射系数误差较大，可信度较低。对于涉及海浪监测的应用，会很大程度上减小光电目标检测系统的准确度。
技术实现要素：
为了解决现有技术中存在的上述问题，本发明提供了一种误差较小、可信度较高、监测准确较高的获取高真实感海浪泡沫光散射系数的方法。
为了实现上述发明目的，本发明采用的技术方案是：
一种获取高真实感海浪泡沫光散射系数的方法，包括：
获取海浪参数；
根据所述海浪参数，建立海浪运动模型，得到海浪高度场分布；
根据所述海浪高度场分布，叠加海浪泡沫的分布模型，得到海面泡沫的混合状态分布；
根据所述海面泡沫的混合状态分布，求解矢量辐射传输方程得到海浪泡沫的光散射系数。
进一步地，所述海面泡沫的混合状态分布符合高斯分布。
进一步地，所述海浪高度场分布为：
其中，n和m分别表示fft的网格大小，x为海浪空间指定网格点的位置矢量，h(x,t)为点x处在t时刻的波高，t是时间，k为波数二维矢量，
进一步地，所述海面泡沫的混合状态分布为：
ω1+ω2＝1，
其中，rbubble(x,y1)和rdrip(x,y2)分别表示海浪泡沫中的气泡和水滴在切面(x,y)处的尺寸大小，ω1和ω2分别表示气泡和水滴的混合比例，μ为海面波峰处的横坐标，h为海面的纵坐标，s为半波宽度，ηb和ηd分别为气泡和水滴的空间尺寸偏差常数，c1和c2分别表示气泡和水滴的体积浓度，t为温度常数。
进一步地，所述矢量辐射传输方程为：
其中，iλ为路径s上沿着方向
进一步地，所述海浪参数包括：海况参数、海浪运动参数、海浪泡沫参数和天气参数。
进一步地，所述海况参数包括海浪的道氏波级、近海面风速和海表温度廓线。
进一步地，所述海浪运动参数包括海浪的运动方向、风阻大小和风向。
进一步地，所述海浪泡沫参数包括泡沫的尺寸、浓度、混合比例和海水盐度。
进一步地，所述天气参数包括风速、压强、温度和湿度。
本发明从构建海浪高度场出发，考虑不同海况、天气条件、不同泡沫尺寸及混合比例对光在泡沫中传播的影响，建立气泡水滴混合状态的海浪泡沫的消溶模型，生成不同混合模式下海浪泡沫的物理模型，即海面泡沫的混合状态分布，继而构建矢量辐射传输方程，实现海浪泡沫光散射系数的实时计算，获得误差较小的海浪泡沫散射系数，可信度较高，提高光电目标检测系统的准确度，从而实现对海浪的准确监测。
附图说明
图1所示为本发明的获取高真实感海浪泡沫光散射系数的方法流程图。
图2所示为气泡浓度的运动范围划分。
图3所示为海浪泡沫分层示意图。
具体实施方式
下面结合具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例，凡基于本发明内容所实现的技术均属于本发明的范围。
实施例一
图1所示为本发明的获取高真实感海浪泡沫光散射系数的方法流程图，包括：
s1、获取海浪参数；
在一个具体实施方式中，所述海浪参数包括：海况参数、海浪运动参数、海浪泡沫参数和天气参数。所述海况参数包括海浪的道氏波级、近海面风速和海表温度廓线。所述海浪运动参数包括海浪的运动方向、风阻大小和风向。所述海浪泡沫参数包括泡沫的尺寸、浓度、混合比例和海水盐度。所述天气参数包括风速、压强、温度和湿度。
s2、根据所述海浪参数，建立海浪运动模型，得到海浪高度场分布；
本实施例使用fft(fastfouriertransformation快速傅里叶变换)方法进行海浪运动模型的建立，具体如下：
任意时刻的浪高h(x,t)对水平位置x＝(x,z)的所有正弦函数求和，可表示为：
式中：x为海浪空间指定网格点的位置矢量，h(x,t)为点x处在t时刻的波高，t是时间，k为波数二维矢量，k＝(kx,ky)；kx＝2πn/lx；kz＝2πm/lz，lx和lz是在x和z方向上的海平面范围，n和m是区间[-n/2,n/2]和[-m/2,m/2]内的整数点个数，以及n和m是fft的网格大小。
其中，
式中：
fft变化后，可以得到海平面上离散点x＝(nlx/n,mlz/m)的高度。
进一步，根据高斯随机数以空间分布的形式生成海浪高度场，波高可以在傅立叶域中快速实现：
式中：ξr和ξi是由高斯随机生成程序产生的平均值为0且均方误差为1的随机数。海浪的波高振幅大小遵循瑞利分布，由两个独立标准正态分布随机数构成。
波高的傅里叶振幅可以用下式表示：
式中：
将上式带入式：
对幅值通过ifft(逆变换快速傅里叶)算法计算，在一定范围内计算不同时间海拔顶点的高度谱，并获得初步的海浪高度场：
其中，n和m分别表示fft的网格大小，x为海浪空间指定网格点的位置矢量，h(x,t)为点x处在t时刻的波高，t是时间，k为波数二维矢量，
由于通过fft方法得到的海面高度场公式是使用的随机的方法生产一定区间内的高度场，对该公式进行修正如下：
采用p-m谱(p-mspectrumseawavemodel)进行采样，利用ittc(internationaltowingtankconference，国际拖曳水池会议)推荐的方向扩展函数构造方向谱，再将方向谱转化为波数谱，得到海浪高度场分布为：
其中，
s3、根据所述海浪高度场分布，叠加海浪泡沫的分布模型，得到海面泡沫的混合状态分布；
s2步骤得到了海浪高度场分布，此时需要在当前基础上叠加海浪泡沫的分布模型，波浪之上的泡沫是由破碎的海水和空气混合而成的白沫，白沫的密度远小于海水，因此白沫会漂浮于海浪之上。参看图2，图2示出了海浪泡沫的流动示意图，其中黑色实心圆形表示水滴，空心圆形表示气泡，在海浪的波峰产生水滴与气泡组成的泡沫，该区域为混合密集区；泡沫沿着海浪方向进行扩散，泡沫密度逐渐变小，该区域为扩散区，直到泡沫消失，即波谷处的平流消溶区。参看图3，海浪水体层3为海水层，泡沫覆盖层2和泡沫覆盖层1为形成海浪泡沫层，且泡沫覆盖层1的泡沫尺寸大于泡沫覆盖层2的泡沫尺寸，对海浪泡沫层进行分析后，最终得到海浪泡沫的分布情况如下。
本发明实施例中将水滴和气泡分别假设成实心与空心刚性球体，对海浪泡沫进行几何建模，得到海面泡沫的混合状态分布为：
ω1+ω2＝1，
其中，rbubble(x,y1)和rdrip(x,y2)分别表示海浪泡沫中的气泡和水滴在切面(x,y)处的尺寸大小，ω1和ω2分别表示气泡和水滴的混合比例，μ为海面波峰处的横坐标，h为海面的纵坐标，s为半波宽度，ηb和ηd分别为气泡和水滴的空间尺寸偏差常数，c1和c2分别表示气泡和水滴的体积浓度，t为温度常数。
s4、根据所述海面泡沫的混合状态分布，求解矢量辐射传输方程得到海浪泡沫的光散射系数；
本发明通过微分和积分的数值方法得到矢量辐射传输方程的各阶解，将各阶解代入海面泡沫的混合状态分布表达式中，得到海浪泡沫的光散射系数。
其中，所述矢量辐射传输方程为：
其中，iλ为路径s上沿着方向
本发明实施例从构建海浪高度场出发，考虑不同海况、天气条件、不同泡沫尺寸及混合比例对光在泡沫中传播的影响，建立气泡水滴混合状态的海浪泡沫的消溶模型，生成不同混合模式下海浪泡沫的物理模型，即海面泡沫的混合状态分布，继而构建矢量辐射传输方程，实现海浪泡沫光散射系数的实时计算，获得误差较小的海浪泡沫散射系数，可信度较高，提高光电目标检测系统的准确度，从而实现对海浪的准确监测。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。