简介:3D地震波走时计算是偏移、反偏移、层析等诸多地震勘探技术中的重要中间步骤。快速推进法计算3D地震波走时具有高效率、稳定性及适应能力强的特点,但快速推进法在震源附近区域的计算精度不高,降低了整个走时算法的计算精度。本文提出了一种联合3D走时计算方法来解决这一问题。该方法在震源附近小范围内使用计算精度较高的波前构建法计算走时,在剩余区域使用快速推进法计算走时,由于模型中绝大多数网格节点走时是通过快速推进法计算的,故新方法保留了快速推进法高效的特点,同时由于震源附近网格节点走时精度的提高,整个新算法的计算精度相对于快速推进法而言有了较大的改善。文中通过数值分析对上述结论进行了验证并使用三维岩丘模型验证了新方法的稳定性和适应能力。
简介:为实现对无人机遥感影像大规模光束法区域网平差的快速解算,采用逐点消元法并利用分块稀疏矩阵保存消元后的法方程,以减少内存的使用量;利用预处理共轭梯度算法实现快速解算,利用OpenMP技术实现预处理共轭梯度迭代计算多核并行处理。结果表明,将消元法和预处理共轭梯度算法用于无人机遥感影像大规模光束法区域网平差解算,既可节省内存,又可提高计算效率。
简介:利用深圳市2013年土地利用变更调查数据,结合景观格局指数与InVEST模型,构建适合快速城市化地区的湿地生态安全评价指标体系,综合评价深圳市湿地生态安全。研究结果显示,深圳市湿地破碎化严重;大鹏新区湿地的生态安全值为0.88,表明该区的湿地安全;盐田区湿地的生态安全值为0.62,表明该区的湿地较安全;坪山区、南山区、福田区湿地的生态安全值分别为0.48、0.44和0.43,表明该区的湿地处于低度预警状态;新宝安区、新龙岗区、罗湖区、光明新区、龙华区湿地的生态安全值分别为0.39、0.36、0.33、0.28和0.21,表明该区的湿地处于中度预警状态。深圳市湿地的生态安全受人为胁迫严重,亟待加强对深圳市湿地生态系统的保护与管理。
简介:利用自动气象站逐小时观测资料、实时探空观测资料及NCEP/NCAR(NationalCentersforEnvironmentalPrediction/NationalCenterforAtmosphericResearch)的1°×1°再分析资料,对2016年8月23—24日新疆维吾尔自治区巴州库尔勒地区一次罕见的短时强降水过程主要的环流系统、水汽输送、动力及湿位涡与垂直螺旋度的变化特征进行了诊断分析。结果表明:此次短时强降水天气过程发生在东部型南亚高压显著增强和西西伯利亚至巴尔喀什湖的长波槽缓慢东移的环流背景下,低层风向风速迅速辐合为此次短时强降水过程的发生提供了水汽和动力条件。此次强降水过程的水汽来源主要包括3个部分:乌拉尔山脊前偏北风引导冷空气南下与西风气流汇合于南疆西部地区的西路水汽输送、青藏高原西南侧低涡前部西南气流引导的西南路水汽输送及西太平洋副热带高压引导的偏南水汽输送,其中9426%的水汽来源于偏西与偏南气流。短时强降水过程发生前期,暴雨区上空左右两侧形成的中尺度环流圈是此次短时强降水过程发生的主要动力机制;垂直螺旋度的发展演变与强降水密切联系,当高层负的垂直螺旋度与低层正的垂直螺旋度配置耦合时,有利于短时强降水的发生。短时强降水过程发生在θse线密集且陡立的区域内,高层高值MPV1的下传触发了位势不稳定能量的释放,促进了强降水的产生。
简介:随着全张量重力梯度(FTG)测量技术的不断发展,重力梯度数据的三维反演技术在油气和矿产勘探中日益受到重视与关注。为了快速处理和解释大规模的高精度数据,图形处理器GPU(GraphicsProcessingUnit)和预处理分解技术(Preconditioningmethods)在地球物理反演中的使用变得十分重要。本文结合对称逐次超松弛(SSOR)技术与不完全乔列斯基分解共轭梯度算法(ICCG)提出改进的预处理共轭梯度法,并考虑到方法预处理分解占用额外的时间,开发该算法的GPU并行算法来提高加速效果。然后通过含噪的模型数据反演来证明改进的并行预处理方法在三维全张量重力梯度数据反演中的适应性。由此,基于NVIDIATeslaC2050GPU的并行SSOR-ICCG算法和在2.0GHzCPU上的串行程序比较,达到了大约25倍的加速比。最后,我们将该算法应用于美国路易斯安那州南方Vinton盐丘的实测航空重力梯度数据反演中,反演出良好的反演结果,验证了该方法在三维全张量重力梯度数据快速反演中的优势和可行性。
简介:利用NCEP再分析资料和常规观测资料,对2016年8月广西地区一次东风波暴雨过程成因进行分析,并基于多种物理量分析东风波及其诱生低涡的动力、热力特征。结果表明,东风波及其诱生低涡是此次暴雨过程的主要影响系统,对流层高层强辐散、低层强辐合的散度场垂直结构,低层强烈的水汽输送和积聚,以及中低层暖湿的大气状态是东风波及其诱生低涡西移加强,并造成相应区域强降水的重要原因。另外,利用重新定义的热力螺旋度物理量,分析其高值区与强降水落区在出现时间和空间上的对应关系,发现两者具有较好的协调性和一致性。因此,热力螺旋度的分布与演变情况对于暴雨预报和研究具有一定的指示意义。
简介:为深入了解青藏高原上中尺度对流复合体(MCC)及其向中尺度对流涡旋(MCV)转化过程中的动热力结构特征演变和发展机理,利用NCEPFNL再分析资料、FY-2E及TRMM卫星资料,对2013年7月22—23日青藏高原上的一次MCC转化MCV的过程进行数值模拟,对其发生的环境背景,以及过程中的涡度、温度和能量收支演变等进行诊断。结果表明:低层水平涡度向垂直涡度的转换,以及垂直方向上正涡度的输送,形成了垂直方向上有利于对流涡旋发展的正反气旋性环流配置。转化过程中大气中上部温度正异常主要来自于可分辨的凝结,温度升高使得高层等压面抬升;下方冷却异常主要来自于蒸发作用和垂直运动,低层温度降低引起等压面收缩下降,这样的配置有利于涡旋发展、对流上升运动加强以及降水发生。对流活动释放的潜热能是转化过程中的主要能量来源,高空急流入口区直接热力环流引发的有效位能向动能转化,也为MCC向MCV转化提供了能量。
简介:使用地面降水观测资料和高空探测资料,运用线性倾向估计法和相关性分析法,分析了2000—2012年武汉地区梅雨期(6—7月)暴雨的气候特征,并建立了暴雨天气预报方程。结果表明,6—7月武汉地区各旬暴雨日数呈先升后降趋势,其中最大值出现在6月下旬;暴雨集中期为6月下旬—7月中旬,占梅雨期总暴雨日数的61%;暴雨局地性强,雨强呈上升趋势。筛选出武汉站强天气威胁指数、百色站850hPa层的露点温度和怀化站850hPa层的风向、风速作为预报因子,建立武汉地区梅雨期暴雨预报方程,实际预报效果检验结果表明该方程的预报质量优于实际业务中常用的数值预报产品。
简介:基于SA多普勒天气雷达资料,使用滤波法、插值法和功率法对径向干扰回波进行剔除,并对3种方法的去除效果进行对比。结果表明:(1)滤波法适用于晴空径向干扰回波的去除,在有降水特别是降水回波与干扰回波混合的区域滤波法效果较差;(2)若只有单一方位存在径向干扰回波时使用插值法去除回波效果较好,而多个方位存在干扰回波时插值法效果较差;(3)功率法通过处理雷达接收功率,对连续多方位存在的干扰回波,以及径向干扰回波与降水回波重叠的情况,都有较好的识别及去除干扰回波的效果。在有降水时,降水回波越弱,对干扰的订正越强;降水回波越强,对干扰的订正越弱;晴空时,径向干扰回波基本被全部去除。
简介:在海洋磁力测量中,因受测量环境及各项改正的残余量等因素的系统性影响,不同测线间存在着不同程度的系统差,这种测线系统差会使测区磁场水平不统一,导致磁异常等值线图失真.为此,研究了基于虚拟检查线的海洋地磁测量数据调差方法,该方法尝试在不使用检查线数据情况下,沿垂直于测线方向选取一定数量的虚拟检查线,根据测线数据拟合出各条虚拟检查线上场的分布规律,并将测线方向拟合前后残差的平均值作为各条测线上的测线系统差,对各条测线数据进行改正.仿真计算与分析表明:对于磁异常变化较平缓的测区,该方法可以有效消除测线系统差;当磁异常变化剧烈时,该方法可以部分消除测线系统差.
简介:非均质高压气藏动态储量评价存在两大问题:(1)对于高压气藏而言,岩石弹性能量不可忽略,不考虑岩石弹性能量计算的动态储量相对于考虑岩石弹性能量所计算的动态储量,误差高达30%;(2)储层非均质性强,连通程度存在差异导致传统物质平衡法不适应,动态储量难以准确计算。针对以上问题,结合四川盆地安岳气田磨溪区块下寒武统龙王庙组气藏的地质与动态特征,应用最小二乘法确定并评价岩石有效压缩系数;划分连通单元,根据连通单元的特点有针对性地选取相应的动态储量计算方法,形成了以连通单元的划分为基础,实测压力梯度约束下的流压折算与全生命周期压力历史拟合相结合,全生命周期不稳定试井分析法、现代产量递减分析法与异常高压物质平衡法综合应用的非均质气藏动态储量评价方法,并对不同方法的适用性及影响结果准确性的因素进行评价,为非均质高压气藏动态储量评价提供了技术参考。