基于起伏地层的曲面重力场快速高精度正演     

郑翾宇  柳建新  陈波  陈龙伟  郭荣文 《物探化探计算技术》2019,(2)

Parker公式在起伏地层重力场正演方面以其简洁、快速等优点被广泛应用于地形改正以及界面反演等,但传统的Parker公式正演方法存在一定的缺陷,在起伏地层模型重力场的正演计算中,由于正演式中e的指数项数值稳定性较差,导致正演结果精度不高。在前人研究的基础上,对Parker公式进行了改进,通过增加上、下界面平均值的方法,提高正演的数值稳定性,改进后的正演精度得到大幅提升。在此基础上,提出了一种曲面观测算法,将观测面由传统的水平面推广到任意起伏面,实现了上下界面起伏地层在起伏观测面上的快速高精度重力正演计算和带地形的密度界面起伏模型的快速高精度重力正演模拟。为了保障数值精度,采用了Gauss-FFT算法,模型实验均取得很好效果。  相似文献   

直流电阻率2.5维有限元数值模拟     

彭艳华  柳建新  刘海飞  孙丽影 《地球物理学进展》2014,(4)

从稳定电流场满足的基本微分方程出发,推导了地面点源电场电位边值问题对应的等价变分问题.用网格离散积分区域时,以四边形为最小研究单元.针对电场在电源附近衰减快,变化梯度大,数值模拟结果在电源附近位置有较大误差问题,采用边界/奇异校正算法改善模拟结果在电源附近的奇异性问题,同时减少边界条件的影响,提高模拟结果精度.从波数域电位求取空间域电位,利用最优化离散波数进行反付氏变换.模型算例结果验证了该算法的可行性与有效性.  相似文献   

四川木里县项脚沟“7·5”特大型泥石流特征及发展趋势分析     

杨华铨  柳金峰  孙昊  赵万玉  张文涛 《中国地质灾害与防治学报》2024,(1):100-107

2021年7月5日,凉山州木里县项脚沟暴发特大型泥石流灾害。文章通过野外调查和特征参数计算,研究了本次泥石流的形成条件、暴发过程和暴发特征。根据实地考察判断,泥石流主要在森林火灾、短时强降雨和沟道地形条件等因素共同作用下暴发,为暴雨径流冲刷引发的火后泥石流。暴雨形成的洪水冲刷坡面、侵蚀沟道,导致沟道两岸坍塌滑坡,堵溃效应明显,泥石流规模扩大。项脚沟过火面积达74.61%,高烈度火烧区面积达57.98%,泥石流临界降雨强度为77.84 mm/h,累计降雨达141.60 mm。根据泥石流特征参数计算结果,此次泥石流密度范围为1.83～1.93 g/cm³,属黏性泥石流,主沟下游出口流速为7.22 m/s,峰值流量759.08 m³/s。结合雨洪法和形态调查法结果分析,此次泥石流重现周期为百年一遇。结合泥石流发育趋势,文章认为流域仍有可能暴发大型泥石流,提出了上游固坡、中游调控、下游排导的防治建议。  相似文献   

地形畸变对偶极-偶极电阻率测深拟断面图的影响   总被引：4，自引：3，他引：1       下载免费PDF全文

吴小平  柳建新  段无悔  蓝燕 《地球物理学进展》2010,25(2):635-641

偶极－偶极剖面方法是一种非常流行的物探方法,但就地形对该方法的实质性影响及修正研究得比较少.我们就几种典型的地形畸变情况进行了正演计算,并绘制了拟断面图.通过对这些图形分析,可以看到在地形畸变处的拟断面图不能反映地下地质体的情况,且有效拟断面图是不连续的,需要分段处理,分级解释.  相似文献   

基于FPGA的CT重建加速技术综述   总被引：1，自引：1，他引：0  

邓靖飞  李建新  李磊  闫镔 《CT理论与应用研究》2010,19(2):25-33

图像重建速度是X-CT成像系统中的一个重要指标。随着X-CT系统的最新发展,对重建物体规模以及重建速度要求的增加,CT重建加速的研究已成为CT研究的热点之一。本文通过软件加速与硬件加速两方面概述CT重建加速的几种主要方法,并结合FPGA的结构特点论述对CT重建加速的适用性以及相对于其他方法的优势。介绍了用FPGA实现重建加速的基本机制以及国内外的研究情况,并对数据结构、数据访问与存储的效率、高级语言等关键问题进行了分析。最后对用FPGA实现CT重建加速的前景进行了展望。  相似文献   

改进的小波分频重构算法在石油地震勘探中的应用     

柳建新  李杰  杨俊 《地球物理学进展》2010,25(6)

在地震资料处理过程中,为了提高信号的信噪比以及保留数据的可靠频段,通常采用频率域滤波的方法.但是,传统的频率域滤波方法是基于傅里叶变换,在窗口与窗口之间的接合部位会出现假异常.且在滤除噪声的同时,很可能同时也将有用信息一起滤除掉. 本文针对上述问题,将小波域分频处理与重构的方法应用到石油地震勘探数据处理中.在提高信号信噪比的同时又可兼顾其它频段信息,使处理出的地震剖面信噪比提高,频带有所拓宽,从而达到最佳重构.文中先后给出了小波变换及其逆变换的数值计算方法以及最佳小波和最佳重构系数的求取方法.最后通过实例分析,其结果表明,该方法在理论上是正确的,实际应用效果也是成功的.  相似文献   

海拉尔盆地乌尔逊-贝尔凹陷挤压构造的发现及其地质意义          下载免费PDF全文

刘志宏柳行军  王芃吴相梅  朱德丰万传彪 《地学前缘》2009,16(4):138-146

海拉尔盆地是叠置于内蒙古-大兴安岭古生代碰撞造山带之上的中-新生代盆地,总体上呈北东向展布,由三坳、两隆五个一级构造单元组成,自西向东依次为扎贲诺尔坳陷、嵯岗隆起、贝尔湖坳陷、巴彦山隆起和呼和湖坳陷, 其中乌尔逊凹陷、贝尔凹陷是位于贝尔湖坳陷凹陷南部的两个二级构造单元。研究表明,乌尔逊-贝尔凹陷在早白垩世大磨拐河期-伊敏期早期、伊敏期末期分别经历了北西 南东向、近东西向两次较长时间的挤压作用,使伸展断陷中的部分控陷正断层发生反转作用,在研究区形成了新的断层传播褶皱、断层转折褶皱以及叠瓦构造和双重构造、突发构造、构造三角带等。由于挤压构造的几何学、运动学特征的不同及隆升速率与同构造生长地层沉积速率的比值在空间上的差异,同挤压地层层序在不同构造部位也表现出不同的沉积特点,其构造特征与造山带十分相似。乌尔逊-贝尔凹陷中两次挤压构造的发现,说明中国东北地区的地壳厚度在早白垩世并非一直处于伸展减薄状态,而是在上述两次挤压时期地壳处于叠置增厚的特点,在空间上形成挤压坳陷型盆地和陆内(板内)造山带相间分布的构造格局。  相似文献   

基于VC++大地电磁数据处理软件开发   总被引：4，自引：3，他引：1       下载免费PDF全文

谢维  柳建新  张东风  王贵财  童孝忠 《地球物理学进展》2009,24(6):2242-2249

研究了GDP大地电磁系统AVG的数据格式,根据实际需要、结合先进的计算机技术和当今比较流行的Visual C++6.0开发工具开发了相应的软件—GDP大地电磁数据处理软件,其中主要功能有：(1)多类型参数数据的曲线显示、参数值的实时显示以及读取存储功能;(2)可视化数据编辑功能;(3)电阻率深度提取功能.笔者在软件系统的设计中,嵌入了模块化思想,为今后完善和移植GDP大地电磁系统数据处理软件奠定了基础.最后通过GDP大地电磁系统数据处理软件对青海某矿区的实测CSAMT资料进行校正处理,效果明显,说明GDP大地电磁系统数据处理软件有较强的实用性.  相似文献   

复杂二维/三维大地电磁的有限单元法正演模拟策略   总被引：1，自引：0，他引：1  

童孝忠  柳建新  郭荣文 《CT理论与应用研究》2009,18(1):47-54

复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性。对于复杂的二维和三维大地电磁正演问题,我们采用有限单元法进行求解。有限单元法最后形成一个线性方程组,系数矩阵是大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点;为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分;在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维和三维模型电磁响应的计算,验证了算法的正确性。  相似文献   

基于多种探测资料对武汉一次短时暴雪天气的监测分析   总被引：2，自引：0，他引：2  

李德俊  唐仁茂  向玉春  柳草  陈英英  徐桂荣 《高原气象》2012,31(5):1386-1392

利用Thies Clima激光雨滴谱仪、MP3000A微波辐射仪、多普勒雷达和人工加密观测资料,分析了2011年2月12日武汉一次暴雪天气过程的演变特征。结果表明:(1)在这次短时暴雪过程中,先后出现了降雨、雨夹雪和纯降雪3个阶段,激光雨滴谱仪监测到该降雨阶段的雨滴谱较宽,最大直径达4.5mm,数浓度比较小,为2～5 429个.m-3.mm-1,雨滴谱型呈双峰型分布特点;(2)微波辐射仪对水汽相态监测显示,08:12-09:34(世界时)期间在1.5～3.5km过冷层比湿维持14～16g.kg-1大值区,之后迅速由雨夹雪过程(仅16min)过渡到纯降雪阶段,在65min内武汉站降雪量达到5.035mm,占整个降雪量的80%以上;(3)雷达反射率因子与粒子谱宽和数浓度密切相关,若后两者越大则对应的雷达回波强度也越大;同时还发现,与微波辐射仪反演的水汽密度、激光雨滴谱仪测得的降水强度也有较好的对应关系,若将这些结合起来对降水的相态变化、强度、量级和持续时间有很好的监测能力。  相似文献