苏北近岸海域表层沉积物粒度及其对环境动力的响应

对2007年9—11月在南黄海辐射沙脊区和江苏北部潮滩采集的表层底质样品进行粒度分析,获得了平均粒径、分选系数、偏态、峰态等粒度参数,并对其进行了聚类分析。结果表明,9 m以浅和20 m以深的样品可归为两种类型的沉积物,而水深9～20 m之间的沉积物两种类型都有,多样性较大。由Gao-Collins粒径趋势分析模型获得的沉积物净输运趋势显示,除外毛竹砂处沉积物有向外海输运的趋势之外,其他外围的断面均显示泥砂向辐射沙洲中心输运,这与历史资料中记录的结果相一致,表明粒径趋势分析适用于大型辐射沙洲海域。根据Stern-berg图分析,粒度参数的空间分布是由差异性输运造成,沉积物从射阳河口岸外向大丰潮间带上部输运的过程中,众数粒径几乎全部以悬移质的形式输运,而小庙洪向东灶潮滩上部输运的众数粒径以悬移方式输运的概率只有60%,从而使射阳河口岸外与大丰潮间带上部的沉积物粒度参数相同,而小庙洪与东灶潮间带上部的沉积物粒度参数有显著差异。     

利用地貌形态相似性评估外版海图水深质量

为检验评估外版海图水深表达海底地形的准确性,利用外版海图和参照数据的水深分别构建Delaunay三角网,对海底特征地貌形态进行提取,设计了位置、表面积、体积相似度指标,对外版海图和参照数据表达同一处特征地貌的相似程度进行了量化。实验结果表明:所提方法可以有效提取海底地形中的特征地貌,所提相似度指标可以较为合理地反映外版海图和参照数据表达同一处特征地貌的相似程度,可为把握外版海图水深表达海底地貌的准确性、进而定量评估其质量提供参考依据,克服了传统经验方法只能进行定性评估分析的不足。     

高分六号遥感影像水深反演能力评价

水深是重要的海洋要素,水深遥感反演是获取浅水水深的重要手段。当前水深遥感反 演应用以国外卫星数据为主,国产卫星数据的研究和应用较少。本文针对国产高分六号卫星 （GF-6） 数据,以三亚南山港为研究区域,分别建立单波段回归模型、双波段比值模型、多波段 回归模型,进行多光谱影像的水深反演能力研究,并与国外主流哨兵2 号卫星（Sentinel-2） 数 据进行实验比较。实验结果表明：GF-6 遥感影像具有较好的浅水水深反演能力和一定的反演精度,各波段水深探测能力依次为：绿波段跃蓝波段跃红波段跃近红外波段,反演方法效果依次为：多波段模型跃双波段模型跃单波段模型。相较于Sentinel-2 数据,GF-6 数据水深反演精度与其一致,这表明GF-6 影像具备替代国外遥感数据进行水深反演的能力和大规模应用的潜力。本文针对GF-6 影像水深反演能力的研究方法和分析,结果将为国产高分系列卫星数据的水深反演研究和应用提供有益的参考。     

基于梯度提升决策树算法的水深反演研究

传统的水深测量方法多通过舰载声纳实地探测的方法,灵活性较差且水深资料更新周期长,并且在某些海域,船只往往难以靠近从而无法完成测量。本文使用七连屿海域附近的WorldView-2多光谱遥感影像构建了基于梯度提升决策树(Gradient Boosting Decision Tree, GBDT)算法的水深反演模型,并利用单波束与人工测量相结合的水深数据,与传统的单波段模型、双波段模型以及BP神经网络水深反演模型的水深数据进行了水深反演精度对比。结果表明,在0～20 m深海域,GBDT模型反演精度高于其他模型,且更符合实际水深,其检验点的R²为0.9664, RMSE为0.94 m, MAE为0.75 m, RME为19%。     

一种局部海域扰动重力数据反演海底地形的非线性序列相关方法

在不考虑地壳均衡条件下,从重力场基本理论出发推导了重力异常/扰动与海深信息之间的近似解析关系,其中扰动重力反演海底地形的模型较之重力异常更加简洁.通过对反演模型级数展开式的深入分析表明,在重力数据分辨率数值小于海深分辨率时并不明显收敛,因此在用重力数据反演高分辨率海底地形时应顾及展开级数的高阶项并采取相应处理方法.在大量统计分析基础上,构建了局部海域扰动重力数据反演高分辨率海底地形的非线性序列相关方法,该方法利用少量船测重力/水深测线数据获取非线性相关参数,而后利用卫星测高重力数据按非线性相关函数模型反演海底地形.试验分析表明,2次非线性相关函数反演效果表现较优,对于海底地形平缓海域,1°×1°区域内单条船载测线数据获得的相关参数即可达到2%相对精度.对于海底地形复杂海域,论文方法适用范围应该尽量缩小,30′×30′区域范围基于单条船载测线数据反演的相对精度可优于7%.从论文试验结果分析,海底地形与重力数据之间存在2次非线性相关函数特征,且这种特征在高分辨率情况下仍然适用,利用非线性函数模型并结合稀疏船载测量数据、密集卫星测高重力数据可为局部海域高分辨率海底地形反演提供一种可行的解算方...     

乌梁素海沉水植物群落光谱特征及冠层水深影响分析

沉水植物对于改善富营养化水体和重建水生生态系统起着至关重要的作用.应用遥感技术可以实时、大面积监测沉水植物的分布和生长情况,而冠层水深直接影响沉水植物在湖泊、河流中的准确遥感解译.本研究基于实测光谱数据,分析了乌梁素海沉水植物光谱特征,并研究了冠层水深对乌梁素海沉水植物反射光谱的影响,建立了乌梁素海沉水植物冠层水深反演模型.结果表明:1)挺水植物在短波红外1662 nm和2223 nm附近分别有一个反射峰,这是挺水植物区别于沉水植物和漂浮藻类的重要波段; 0深度沉水植物(WDC=0)与漂浮藻类的光谱反射率非常接近,但是在绿波段(550～690 nm)有明显差异,因此,可以利用绿波段和短波红外波段的光谱特征来区分挺水植物、沉水植物和漂浮藻类.2)沉水植物群落的光谱反射率随冠层水深的增加而降低,在700～900 nm波段范围内变化最为明显,且在700～735 nm波段附近,沉水植物群落光谱反射率与冠层水深呈显著负相关.3)在建立的单波段/波段比沉水植物冠层水深反演模型中,波段比反演模型要优于单波段反演模型,波段比反演模型的决定系数R2 0.70,均方根误差13.70 cm,平均相对误差28%,反演精度较好,适用于10～60 cm沉水植物冠层水深的反演.4)利用波段响应函数,将实测光谱反射率积分到Landsat-8 OLI波段上,建立OLI了冠层水深反演模型,其中,波段比幂函数模型反演效果最好,R2为0.49,均方根误差为18.17 cm,平均相对误差40.05%.可用于精确大气校正后乌梁素海沉水植物冠层水深的反演.     

基于共轭梯度算法的重力梯度数据三维聚焦反演研究

重力梯度数据相对于传统重力数据,能够更细致、准确地描述地球浅部构造和研究矿产资源分布等信息.本文采用共轭梯度算法,在加权密度域求解重力梯度数据三维聚焦反演最优化问题,以恢复地下三维密度分布,目标函数包括数据不拟合函数和最小支撑稳定函数.首先,在推导目标函数对加权密度的一阶导数时,为了得到更合理的计算公式,我们考虑变加权函数中含有密度变量;此外,本文通过密度上下限约束,改善了传统聚焦反演中聚焦因子选取困难的问题.新算法获得的反演结果,对聚焦因子的选择约束较少,相比传统聚焦算法,能够更容易的获得理想结果.将方法应用于理论模型验证其有效性和正确性,并应用本文方法处理文顿盐丘地区的航空全张量重力梯度数据,得到了与已知地质信息匹配的密度分布,表明本文方法具有处理实际数据的能力.     

基于印模法的地面磁电阻率数据三维非线性共轭梯度反演

本文提出了能提高异常体分辨能力,同时得到绝对电导率的地面磁电阻率数据三维反演方法.磁电阻率响应用准直流的低频磁场代替;数值模拟由频率域电场满足的Helmholtz方程出发,采用三维交错网格有限差分法;长直导线源作为发射源,其中源的计算包含在背景场中;结合地面磁电阻率数据各分量的特点,选择y分量进行反演研究;反演采用三维非线性共轭梯度反演技术,为了提高异常体的深度分辨能力,进行迭代重构反演;用印模法对初始模型进行重构,采用的是辅模型在浅部,元模型在深部的组合方式.从合成数据和实际数据的反演结果可以得到以下的认识:(1)由频率域麦克斯韦方程组出发,低频磁场数据反演可以直接得到电导率,而不是相对电导率之比;(2)采用印模法组合初始模型,进行迭代重构反演,可以提高地面磁电阻率数据反演对异常体的分辨能力,确定埋深位置,同时不会丧失对于浅部异常体的分辨能力;(3)在结合印模法的地面磁电阻率数据三维反演中,深部异常体的分辨能力受地表不均匀导电体影响较小;(4)确定印模深度可以采用上一次重构反演结束时的模型变化量,通过相邻两次重构反演结束时的模型变化量之差来确定迭代重构是否终止.因为静磁场与重力场在数学上的相似性,本文的反演方法可以被运用到重力场等位场的地面数据的反演中.     

地面矩形大定源电磁法频率域三维非线性共轭梯度反演

本文实现了地面矩形大定源三维频率域反演.矩形大定源三维模型响应计算采用交错网格有限差分技术.正演的微分方程为异常电场满足的非齐次Helmholtz方程,方程右手边源项中的大定源产生的背景格林函数由虚界面法结合虚框法计算.频率域三维反演采用非线性共轭梯度反演技术.反演的数据类型为垂直磁场的频率域响应Hz的实部和虚部分量.数值结果表明,(1)三维模型正演模拟数值结果与前人一致,为三维反演奠定基础;(2)针对两个三维导电模型,分别进行了三维反演数值试算.反演结果可以清晰恢复出异常体的电阻率和位置信息,表明地面矩形大定源三维频率域非线性共轭梯度反演具有可行性.本文研究的意义在于,在电磁响应时频转换技术的基础上,如果将野外实测的瞬变电磁数据变换为对应的频率响应,则结合本文提出的三维反演技术,可以为矩形大定源瞬变电磁数据的三维解释提供一个新的思路.     

基于GPU并行的时间域全波形优化共轭梯度法快速GPR双参数反演

探地雷达（GPR）时间域全波形反演计算量巨大,内存要求高,在微机上计算难度大.本文中作者基于GPU并行加速的维度提升反演策略,采用优化的共轭梯度法,避免了Hessian矩阵的计算,在普通微机上实现了时间域全波形二维GPR双参数（介电常数和电导率）快速反演.论文首先推导了二维TM波的时域有限差分法（FDTD）的交错网格离散差分格式及波场更新策略.然后,基于Lagrange乘数法,将约束问题转化为无约束最小问题,构建了共轭梯度法反演目标函数,采用Fletcher-Reeves公式与非精确线搜索Wolfe准则,确保了梯度方向修正因子及迭代步长选取的合理性.而GPU并行计算及维度提升反演策略的应用,数倍地提升了反演速度.最后,开展了3个模型的合成数据的反演实验,分别从观测方式、梯度优化及天线频率等方面,分析了这些因素对雷达全波形反演的影响,说明双参数的反演较单一的介电常数反演,能提供更丰富的信息约束,有效提高模型重建的精度.