一种模拟非均匀介质中弹性波传播新的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,29(1):43-51

尼泊尔在地理位置上处于中国和印度包围之中，由于历史原因，其经济、政治、军事等长期受到制约。随着中国“一带一路”倡议的推进和中印近年来领土和国际政治问题的冲突与分歧，尼泊尔正在成为中印新的角力场。本文首先通过构建地缘影响力模型进行演变格局分析与趋势预测，并加入世界环境和区域环境两方面影响因素，运用定性与定量分析结合的方法，预测中、尼、印三国未来战略行动方向。综合以上分析发现：1）2000年—2016年，中国对尼泊尔地缘影响力较于印度对于尼泊尔地缘影响力上升趋势明显，预计在2022年至2023年中、印对尼泊尔地缘影响力达到平衡。2）硬实力是中国对尼泊尔地缘影响力上升的关键因素；软实力和交流互动力是未来中国增强在尼泊尔地缘影响力的发展因素；地理环境是中国对尼泊尔地缘影响力发展的制约因素。3）未来中国将继续实施积极性合作战略，印度形成对中国遏制性防范战略，尼泊尔维持被动性中印平衡战略。  相似文献   

“一带一路”沿线国家石油资源流动类型演变     

孙琪  赵媛  张新林  夏四友 《世界地理研究》2020,29(2):337-345

综合运用石油自给率、流动比率等方法，以“一带一路”沿线国家为研究对象，对1995—2014年各国石油资源流动类型时空演变进行研究。结果表明：1）1995—2014年石油净补给型国家数量较为稳定，无产型或低产型国家的石油贸易倾向变化相对较小；基本自给型国家数量占比最大，1995年占比1/2，其余四个年份占比稳定在1/3；净支出型国家数量显著增加，至2014年其中有5个国家保持不变，表明这些国家的石油产量具有较大的优势。2）1995年和2014年汇地国家均为39个国家，其中有34个国家一直保持汇地职能，空间分布上呈东多西少，南多北少，东西分化以西亚、中亚、东欧为界；有14个国家一直保持源地职能，其空间分布格局与汇地相应，西多东少，南多北少，主要分布于西亚、东欧及中亚；交流地数量较少，呈分散分布，且分布地不稳定。3）从数量变化来看，研究期内汇地略有增加，源地略有减少，交流地略有增加；“一带一路”源地总出口量与汇地总进口量的比呈逐年降低趋势。  相似文献   

基于GeoStudio的边坡稳定性分析及支护方案选择的理论探讨     

何保  宋帅 《地质与勘探》2019,55(5):1329-1335

本文选择大连市中山区长利巷地区的边坡,通过现场采取岩土样,进行室内试验得出各地层力学物理参数,对边坡进行稳定性分析。针对目前所拥有的模拟软件与方法,分析并筛选出适合该边坡的模拟软件。利用GeoStudio软件建模、模拟,使用Morgenstern-Price方法进行计算,通过分析得到计算结果,该边坡的安全安全系数为1. 127,认为其属于不稳定边坡。对于不稳定边坡,本文对主流的边坡支护方法进行分析筛选,拟采用扶壁式挡土墙支护方案对其进行预防。根据规范设计挡土墙尺寸,并对该挡土墙进行抗倾覆抗滑移验算,得出的验算安全系数满足规范要求,由此可得知该支护方案方法合理,技术可行。  相似文献   

考虑冲刷影响的筒型基础竖向极限承载力研究     

陆罗观  刘润  练继建  陈广思 《海洋工程》2019,37(3):69-77

大直径宽浅式筒型基础,阻水宽度大,在位工作期间受波浪海流作用,其周围土体易被冲刷。为研究单侧地基土体受冲刷后筒型基础的竖向极限承载力变化,通过引进冲刷率的概念,采用有限元方法研究了不同冲刷率下筒型基础的竖向极限承载力;并基于Meyerhof理论建立了计算不同冲刷率下筒型基础竖向极限承载力的极限平衡方法。研究结果表明,随着冲刷率增大,筒型基础的极限承载力出现不同程度的下降,当冲刷率为0.8时,即筒型基础单侧土体冲刷深度达6.4 m时,筒型基础的竖向极限承载力折减率为3.28%。建立的极限平衡算法可准确计算冲刷条件下筒型基础的竖向极限承载力。  相似文献   

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

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

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

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

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

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

基于GRACE RL06数据监测和分析南极冰盖27个流域质量变化   总被引：1，自引：0，他引：1       下载免费PDF全文

高春春  陆洋  史红岭  张子占  徐长仪  谭冰 《地球物理学报》2019,62(3):864-882

本文基于CSR最新公布的GRACE RL06版本数据,采用Slepian空域反演法估算了南极冰盖27个流域的质量变化.Slepian空域反演法结合了Slepian空间谱集中法和空域反演法的技术优势,能够有效降低GRACE在小区域反演时信号出现的严重泄漏和衰减,进而精确获得南极冰盖在每个流域的质量变化.相对于GRACE RL05版本数据,RL06在条带误差的控制上要更加优化,获得的南极冰盖质量变化时间序列也更加平滑,但在趋势估算上差别并不明显(小于10Gt/a).本文的估算结果显示:在2002年4月至2016年8月期间,整个南极冰盖质量变化速率为-118.6±16.3Gt/a,其中西南极为-142.4±10.5Gt/a,南极半岛为-29.2±2.1Gt/a,东南极则为52.9±8.6Gt/a.南极冰盖损失最大的区域集中在西南极Amundsen Sea Embayment(流域20-23),该地区质量变化速率为-203.5±4.1Gt/a,其次为南极半岛(流域24-27)以及东南极Victoria-Wilkes Land(流域13-15),质量变化速率分别为-29.2±2.1Gt/a和-19.0±4.7Gt/a,其中Amundsen Sea Embayment和南极半岛南部两个地区的冰排放呈现加速状态.南极冰盖质量显著增加的区域主要有西南极的Ellsworth Land(流域1)和Siple Coast(流域18和19)以及东南极的Coats-Queen Maud-Enderby Land(流域3-8),三个地区质量变化速率分别为17.2±2.4Gt/a、43.9±1.9Gt/a和62.7±3.8Gt/a,质量增加大多来自降雪累积,比如:Coats-Queen Maud-Enderby Land在2009年和2011年发生的大规模降雪事件,但也有来自冰川的增厚,如:Siple Coast地区Kamb冰流的持续加厚.此外,对GRACE估算的南极冰盖质量变化年际信号进行初步分析发现,GRACE年际信号与气候模型估算的冰盖表面质量平衡年际信号存在显著的线性相关关系,但与主要影响南极气候年际变化的气候事件之间却不存在线性相关关系,这说明南极冰盖质量变化的年际信号主要受冰盖表面质量平衡的支配,而气候事件对冰盖表面质量平衡的影响可能是复杂的非线性耦合过程.  相似文献   

不同非感应起电及感应起电参数化方案对青海东部一次雷暴云电荷结构影响的数值模拟研究          下载免费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样条插值拟合得到,因此在反演中它将会随着反射点位置的更新自然得到更新.与前人提出的保边缘层析算法或多层立体层析算法相比,本文提出的地层格架正则化无需引入混合正则化项或定义某种复杂的混合速度格式,更为直接也更容易实现.理论和实际数据算例证实了该正则化技巧的稳健性和可靠性,能够得到与实际地质构造特征更为一致的地质一致性反演结果.  相似文献