一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

基于点云数据的不规则文物三维重建——以经幢为例     

喜文飞  李灵索  杨勇  李国柱  李春梅  李金平 《测绘通报》2020,(12):87-89

海量点云数据可以较好地构建文物的三维建模，针对文物复杂的表面结构，传统的算法很难达到精细建模，为了解决复杂模型表面的空间三角网构网问题，本文结合空间三角形的法向量构建三角网，以空间三角网表面的曲率为约束条件，构建物体表面的空间三角网，然后对构建的空间三角网进行精细建模。以昆明市文物——经幢为例进行三维重建试验，结果表明，新方法适合室内文物的三维建模，可以达到较好的效果。  相似文献   

融合多源数据的三维建模方法及精度分析     

祖立辉 《北京测绘》2020,(3):347-351

以地面相机影像、低空无人机影像和地面三维激光点云数据为融合对象,研究空地影像融合、低空影像融合地面三维激光点云和低空影像融合三维白模进行三维建模的方法,并对其三维建模精度和建模效果进行综合分析。这3种融合方法进行三维建模丰富了三维模型的色彩纹理信息,提高了三维模型的几何结构精度,能满足较大场景的三维实景模型的建模要求,为数字化古建筑(古文物)保护提供具有较高价值的参考数据。  相似文献   

航空重力垂直梯度探测潜艇方法研究     

张志强  郑晗  崔银锋 《海洋测绘》2019,39(4):6-9

航空重力梯度测量属于被动探测,抗干扰能力强,如果能和其他探潜手段相配合将极大地提高航空搜潜的效率。针对航空重力梯度测量是否能够用于探测潜艇的问题,依据俄亥俄级弹道导弹核潜艇的结构特点,建立了适用于重力梯度计算的潜艇模型,分别给出了潜艇外壳、内部质量亏损产生的重力垂直梯度的计算方法,并对不同精度重力梯度仪可探测的潜艇重力垂直梯度值进行了计算,从航空反潜的角度对探测潜艇效果进行了分析。经计算表明,如航空重力梯度仪精度达到10~(-2) E,将具备一定的实际探潜效能;如精度达到10~(-4) E,反潜机搜索宽度可与现有航空磁性探潜相当。  相似文献   

基于低空摄影测量的机载LiDAR点云修复方法     

李冬森  王永志  陈建忠 《测绘与空间地理信息》2019,42(4)

针对机载LiDAR点云存在数据缺失造成的空洞问题,研究利用低空摄影测量技术,探索机载LiDAR点云空洞修复的方法。利用低空摄影测量手段获取的遥感影像可以生成高精度的修复点云,并通过将修复点云融合到原始LiDAR点云中,实现对机载LiDAR点云空洞的修复。该方法操作简单、效率高,适用于大面积机载LiDAR点云数据的批量修复,能够为城市三维精细化建模提供重要的数据支撑。  相似文献   

非连续变形分析中块体大转动问题研究     

巩师林  凌道盛  胡成宝  钮家军 《岩土力学》2020,41(11):3810-3822

原有非连续变形分析（DDA）采用一阶近似后的位移增量表达式更新块体构形，推导相关子矩阵，且对不同时步计算出的应变增量直接叠加，当模拟的块体发生大转动时往往会产生较大误差。为考虑块体转动与变形的耦合作用，引入先变形、后转动的块体位移增量表达式。重新推导了惯性力子矩阵，将块体转动时的离心力与科氏力加到荷载矩阵中。计算时对应变分量及其相关变量进行坐标变换与修正，并采用新引入的位移增量表达式计算块体顶点位移，进行后接触修正与更新块体构形。数值算例表明，改进后的程序能够消除转动带来的误差，自动考虑了块体转动时离心力和科氏力引起的变形，应变计算精度更高。改进方法克服了块体体积自由膨胀、应变场畸变等问题，给出了合理的块体应变。  相似文献   

GPM和TRMM卫星日降水数据在黄河源区的适用性分析     

刘兆晨  杨梅学  王学佳  程立真 《冰川冻土》2020,42(2):575-586

通过评估GPM计划三种日降水产品（IMERG-E、 IMERG-L和IMERG-F）和TRMM卫星、 两种日降水产品（TMPA 3B42和TMPA 3B42RT）在黄河源及其周边区域38个台站的适用性， 探究了五种产品探测精度和海拔高度及雨强的相关关系， 结果表明： 在与实测资料的一致性和偏差方面， GPM卫星产品要全面优于TMPA产品。在TRMM卫星产品中， 3B42产品明显优于3B42RT。五种产品的相关系数均表现出明显的从东南到西北递减的趋势， 均方根误差北部普遍低于南部。IMERG产品的探测率（POD）和探测成功率（CSI）都要普遍高于TMPA产品， 而误报率（FAR）则是TMPA 产品更低， 表现更好。五种产品均在个别台站出现了严重误报的情况， 这些台站主要分布在研究区的西北部。IMERG三种产品对于海拔高度的依赖程度具有很强的一致性， 而3B42RT产品对海拔高度几乎没有依赖。除3B42RT产品外， 其余四种产品的偏差均随雨强的增加而增大。在探测率方面， IMERG产品对小雨、 中雨和大雨的探测能力均优于TMPA产品。  相似文献   

中国东部层积云发展过程中云微物理特征的演变   总被引：1，自引：0，他引：1       下载免费PDF全文

李思聪  李昀英  孙国荣  宋文婷 《地球物理学报》2019,62(12):4513-4526

基于2007—2010年的CloudSat卫星观测数据,以云层液态水路径为指标将层积云的发展过程划分为五个阶段,对比研究了中国东部降水与非降水层积云发展过程中云微物理特征和云微物理机制的演变,并分析了其海陆差异.研究表明:非降水层积云中,云滴增长主要通过凝结过程完成,但云滴的凝结增长有限,难以形成降水,在非降水层积云发展的旺盛阶段,云层中上部云滴发生较弱的碰并过程.降水层积云中云滴碰并增长活跃,当云层液态水路径小于500 g·m~(-2)时,云滴在从云顶下落至云底的过程中持续碰并,并在云底附近出现云水向雨水的转化;当降水层积云液态水路径超过500 g·m~(-2)时,云滴碰并增长主要发生在云层上部,在云层中部,云液态水含量、液态粒子数浓度和液态粒子有效半径达到最大,云水向雨水的转化最为活跃.层积云微物理特征的海陆差异主要是由海陆上空气溶胶浓度和云中上升气流强度不同导致的.在非降水层积云中下部,陆地丰富的气溶胶为云滴凝结增长提供了充足的云凝结核,因而云微物理量的量值在陆地上空更大,而在云层中上部,云滴凝结增长达到极限,海洋充足的水汽输送使云微物理量的量值在海洋上空更大.当降水层积云液态水路径大于500 g·m~(-2)时,陆地层积云中更强的上升气流使大量云滴在云层中上部累积滞留,云滴碰并增长活跃,云层中上部云微物理量的量值在陆地上空更大.  相似文献   

不同非感应起电及感应起电参数化方案对青海东部一次雷暴云电荷结构影响的数值模拟研究          下载免费PDF全文

李江林  余晔  李万莉  李亚珺 《地球物理学报》2019,62(7):2366-2381

本研究利用加入起电、放电参数化方案的数值模式（Weather Research and Forecasting Model（Version 3.7.1），WRF3.7.1_ELEC），通过设计五组不同非感应起电及感应起电参数化方案敏感性试验，对发生在青藏高原东北部青海大通地区的一次雷暴过程进行模拟研究，对比分析了不同非感应起电机制及感应起电机制对雷暴云电荷结构的影响.结果表明：在雷暴云发展旺盛阶段，Saunders（S91）、Riming Rate（RR）、和Saunders和Peck（SP98）三种非感应起电方案模拟的雷暴云最低层均为负电荷区，而混合方案（Brooks and SP98，BSP）模拟的雷暴云最低层为正电荷区，主电荷区自下而上为"+-+-"排列的四层电荷结构.与甚高频辐射源定位法推算的结果对比，BSP方案模拟的本次高原雷暴云电荷结构更接近实际情况；几种不同非感应起电方案模拟的主电荷区外围与主电荷区电荷结构不同，说明在雷暴发展的不同阶段雷暴云的电荷结构是不同的；几种非感应起电方案模拟的电荷结构不尽相同，主要是由于霰、冰和雪粒子在不同高度所带电荷的极性及电量的大小不同，霰粒子的电荷密度对低层的影响较大，冰粒子和雪粒子的电荷密度对中上层的影响较大；加入感应起电机制后，雷暴云电荷结构分布几乎没有变化，但能使雷暴云发展旺盛阶段低层和中层的正负电荷区电荷密度有所加强.  相似文献   

地层格架正则化约束下的二维立体层析反演          下载免费PDF全文

张力起  杨锴  邢逢源  李振伟  倪瑶 《地球物理学报》2019,62(2):634-647

通过把地层格架信息作用于立体层析Fréchet导数矩阵,使得更新后的速度模型呈现出符合地质规律的块状特征.地层格架信息基于立体层析反演中得到的反射点位置进行非规则B样条插值拟合得到,因此在反演中它将会随着反射点位置的更新自然得到更新.与前人提出的保边缘层析算法或多层立体层析算法相比,本文提出的地层格架正则化无需引入混合正则化项或定义某种复杂的混合速度格式,更为直接也更容易实现.理论和实际数据算例证实了该正则化技巧的稳健性和可靠性,能够得到与实际地质构造特征更为一致的地质一致性反演结果.  相似文献