基于偏振植被指数的植被-土壤混合像元室内实验研究

混合像元是遥感影像中普遍存在的一种现象,对其组成和各成分比例的反演一直是遥感研究中的重难点,而国内利用偏振植被指数对混合像元的研究几乎没有涉及。本次研究通过对不同面积比例的植被-土壤混合像元偏振反射高光谱特征进行分析,讨论不同条件下的植被-土壤混合像元偏振高光谱特性,并利用偏振反射比计算了12种0°偏振态下的偏振植被指数,分别构建了植被指数与植被面积比例以及光谱特征参数与植被面积比例的数学模型。结果表明,混合像元中植被面积比例和偏振角均对其偏振高光谱有一定的影响;865 nm的偏振反射比与植被占像元面积比例的相关性最好,采用多项式进行拟合时,其决定系数达到0.99,适合进行植被占像元面积比例的反演;偏振光谱“红边”处的一阶微分值与植被像元比例存在良好的线性相关,R²=0.974;植被面积比例与植被指数和光谱特征参数呈现良好的相关关系,其中P-DVI和光谱吸收指数（SAI）与植被面积比例的拟合效果最好,决定系数分别为0.99以及0.94,适合进行植被-土壤混合像元中植被面积比例的反演。     

基于交错网格Fourier伪谱微分矩阵算子的地震波场模拟GPU加速方案

作为高精度波形反演或逆时偏移的重要组成部分,地震波数值模拟对计算速度和效率提出了更高要求.GPU通用计算技术的产生及其内在数据并行性,为高效地震波数值模拟应用和研究得以有效开展奠定了基础.本文借助交错网格的Fourier伪谱微分矩阵算子和GPU上高效矩阵乘法,实现了复杂介质中地震波模拟的高效算法.数值试验表明,优化后的GPU计算相比CPU单核计算在大规模二维地震波场计算中获得至少100x以上的加速比.这对我们快速分析目标反射层在地震剖面中同相轴位置,制定优化采集方案具有重要意义.     

双相介质地震波场数值模拟的迭积微分算子及其PML边界条件

将基于计算数学中Forsyte 广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度--应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演模拟方法,来解决复杂非均匀介质模型中的波场传播问题.为了大幅衰减人工边界引起的反射,本文将完全匹配层(Perfectly Matched Layer,PML)吸收边界条件引入到所构建的方法中,以解决迭积微分算子法的边界问题.以二维波动方程为例,用迭积微分算子法实现了双相介质的地震波场正演模拟,模拟结果表明,双相介质模型较好地解释了含流体孔隙特性.同时也表明迭积微分算子法是一种非常实用、有效的数值模拟方法.     

Hi-net速度微分与KiK-net 加速度对比研究

首先介绍了日本防灾科学技术研究所Hi-net和KiK-net台网布设情况和数据格式, 然后利用Hi-net和KiK-net共同记录到的282次M_j≥5.0、 震源深度≤100 km的地震, 分析了Hi-net速度微分与KiK-net加速度对数峰值比、 对数反应谱比的分布特征及其随震级、 距离的变化关系, 讨论了速度微分结果在幅值和频谱两方面的不确定性. 结果表明, 数字型速度记录微分与观测加速度无论是峰值还是频谱均存在一定差异, 不能将其直接用作加速度. 考虑到对数峰值比和对数反应谱比在不同震级档、 距离档存在一定规律, 因此对速度微分结果的峰值和特定周期的反应谱除以一定的系数, 可以降低数字型速度记录数值微分造成的不确定性. 另外, 也可采用概率的形式描述这种不确定性.      

核脉冲信号数字高斯成形模型的建立与仿真

详细阐述了一种以Sallen—Key滤波器为基础的新型数字高斯成形的方法。运用二阶微分方程建立高斯成形数学模型,在模型解算过程中利用数值微分算法得到一个通用的输入信号与输出信号的数值递推解,从而实现了核脉冲信号的数字高斯成形。通过标准指数信号的模拟以及实际测量信号的模拟,确定了数字高斯成形的参数。     

约束优化的微分进化算法在波动方程反问题中的应用

本文提出了一种约束优化的微分进化算法,该算法使得种群在演化过程中能保持较好的多样性,且参数设置简单,不容易陷入局部最优,并能在较短时间内找到约束优化问题的最优解.结合有限元法对波动方程反问题进行研究,在对一维及二维波动方程反问题的数值模拟中都得到了较好的结果,表明了该算法具有较好的稳定性和较强的抗噪能力.     

地震波数值模拟中的最优Shannon 奇异核褶积微分算子

如何有效协调地震波数值模拟过程中的计算精度和计算效率,一直是许多地球物理工作者密切关注且努力寻求解决的问题.这一问题求解的好坏,将直接影响到人们正确认识地下构造的地球物理响应.本文从离散Shannon奇异核理论出发,推导了基于Shannon奇异核的交错网格褶积微分算子,分析并求出了影响算子精度窗函数的最优化参数.通过将...     

基于辛格式奇异核褶积微分算子的地震标量波场模拟

本文在对地震波场进行模拟时,采用辛差分格式对波动方程进行时间离散,采用奇异核褶积微分算子对波动方程进行空间离散.该方法尽管增加了一些计算量,但提高了计算精度和稳定性;相对于其他非辛算法,它是全局保结构的,并且具有较强的长时间跟踪能力.该方法为解决大尺度、长时程地震波场的高精度模拟问题提供了一种新的、有效的选择.     

地震波传播的褶积微分算子法数值模拟

为了解决地震波场数值模拟中的速度与精度的匹配及局部信息与全局信息的耦合等问题, 该文借助新推出的基于Forsyte广义正交多项式的褶积微分算子, 将计算数学中的Forsyte多项式褶积微分算子应用到地震波传播的数值模拟中.复杂非均匀介质模型数值模拟结果说明了该方法的可行性和优越性.该方法同时具有广义正交多项式褶积微分算子的高精度和有限差分短算子算法的高速度.通过对算子长度的调节及算子系数的优化, 可同时兼顾波场解的全局信息与局部信息.      

欧拉梁动力反应初边值问题的微分求积解

动荷载作用下欧拉梁动响应的计算是一个初边值问题,通常很难得到解析解,传统数值方法一般是把空间和时间分别离散进行求解,计算相对复杂,效率也不高.针对分布动荷载作用下欧拉梁的振动偏微分方程,采用传统微分求积法,在空间和时间上同时进行离散;对于所有非0阶的初/边值条件,采用嵌入法在权系数计算中予以考虑.算例的数值结果与精确解的对比证明采用传统微分求积法处理此问题是可行的,而且是高效的.对于实际工程中的其他类似问题,该方法同样适用.