基于Casiofx-9860GSD计算器的公路线路计算软件设计与开发

基于Casio fx-9860 G SD型计算器的特性,笔者利用VS2008开发出能够读取、处理维地道路系统生成的Excel直曲表,并能生成计算器能读取的线元串列数据库文件的软件。然后利用生成的串列文件在计算器中进行线路任意里程正反算。实现了公路平面设计参数的自动转换、输入。本文着重介绍软件开发、工程实例中遇到的问题和解决方法。通过多个工程实践证明,该系统可以很好地帮助测量技术人员快速、准确地进行横断面测量、线路中边桩放样、线路地形图等工作,提高了线路测量的效率和规范性。     

考虑关断时间的回线源激发TEM三维时域有限差分正演

从麦克斯韦旋度方程出发可以直接导出瞬变电磁场扩散方程,然而扩散方程不含电场对时间的一阶导数,不能构成显式的时域有限差分方程,借鉴du Fort-Frankel有限差分离散方法引入虚拟位移电流项构建显式时域有限差分方程.对Wang和Hohmann的经典时域算法进行了两点改进:第一,通过将矩形回线源电流密度加入麦克斯韦方程组的安培环路定理方程,实现回线源瞬变电磁激发源加入;第二,在计算中考虑关断时间.第一点改进使时域有限差分方程考虑了一次场的计算,并且源的计算不再依赖均匀半空间模型响应作为初始条件,使算法能够适应表层电阻率不均匀时的三维复杂模型.由于实际观测中不可能出现阶跃电流的关断形式,第二点改进可以方便设置发射电流下降沿.采用改进的三维时域有限差分正演算法对均匀半空间模型、四类三层模型、均匀半空间中含有低阻块体模型进行了计算并分别与解析解、线性数字滤波解、积分方程解和Wang的三维时域有限差分解进行了对比验证.以H模型为例,采用建立的三维时域有限差分正演算法计算了不同关断时间的斜阶跃脉冲回线源瞬变电磁中心点感应电动势衰减曲线.以实际地质资料为基础,构建包含两层采空区的三维复杂模型,以1 μs的极短关断时间进行了复杂模型定回线源瞬变电磁响应计算,并计算了该复杂模型的视电阻率曲线.     

基于辛算法模拟探地雷达在复杂地电模型中的传播

近年来,探地雷达(GPR)凭借其快速、高效、无破损等特点,已经广泛应用于浅地层目标探测中.数值模拟是研究探地雷达电磁波在地下结构中传播规律的有效手段.辛算法是一种保持Hamilton系统总能量不变的时域数值计算方法.本文提出了基于一阶显式辛分块龙格库塔方法的探地雷达数值模拟方法.通过对比本文算法与时域有限差分方法计算结果可知,在同等计算精度下,本文算法可以节省25%的计算时间.并基于本文算法对两个复杂GPR模型进行正演模拟,得到模拟GPR探测wiggle图,这有助于更好的理解和分析实测雷达数据.     

基于球体层状介质模型的大地电磁正演

目前大地电磁测深法的一维正演理论,是基于平面波垂直入射水平层状介质的假设模型。但由于地球是一个球体,因此有必要研究基于球状介质模型的大地电磁正演理论。这里详细推导了基于球体层状介质模型的大地电磁正演公式,计算了若干理论模型。通过同基于水平层状介质模型的大地电磁正演结果对比,验证了正演公式的正确性。同时,指出当探测周期增加到上万秒时,阻抗相位会增大;而当探测周期增加到数十万秒时,视电阻率会减小。     

基于解释量版的渤海油田巨厚砂砾岩定量表征

渤海油田古近系三角洲沉积发育巨厚（大于100 m）的砂岩和砾岩储层,储层横向变化快,内部非均质性强,地震预测难题大。文章针对此类问题,提出一种基于厚度解释量版的砂砾岩定量表征方法。根据已有钻井资料建立井点处地质模型,结合已有地质资料分析储层可能的变化形式,形成地质模型库,对地质模型库进行正演获得相应的正演地震资料库,然后提取正演地震资料的反射振幅属性,建立不同砂砾岩变化形式下的厚度解释量版,最后提取实际地震资料的反射振幅属性,将其投射到厚度解释量版上,实现对巨厚砂砾岩储层的定量表征。该方法克服了常规地震反演、地震属性等方法难以有效刻画巨厚砂砾岩储层的难题,实现了对具有干涉效应的巨厚砂砾岩储层的厚度及其内部变化的精确定量表征。     

基于模型正演的折射静校正方法研究与应用

随着地震勘探区域地质条件复杂化,静校正在资料处理中的作用日显重要.本文模拟研究了折射渡在不同地质条件中的特性,完成了正演模型下的折射波静校正量计算,并将静校正后记录与理论模型记录对比,分析该方法的准确性,总结折射波静校正方法的技术难点,提出一种更准确优越的解决方法.最后,运用该方法选取最合理的参数,并配合先进的处理软件...     

过套管电阻率测井响应数值模拟研究进展

套管井电阻率测井仪的研制成功是套管井测井技术的重要进展,在老井地层评价中起到了越来越重要的作用.对过套管电阻率测井响应数值模拟,本文综述了近十几年来的研究现状和进展,分析了过套管电阻率测井的测量方法,比较了过套管电阻率测井响应的计算方法,讨论了多层状地层、有限长套管测井响应计算的结果,对我们最近计算的套管接箍测井响进行了分析.针对过套管电阻率测井响应计算问题,指出了在实际应用中尚需解决的实际问题和解决的基本思路,并力图展望今后的研究发展趋势和解决问题的方法.     

低阻覆盖层对高频大地电磁测深的影响

低阻覆盖层对电磁波具有屏蔽效应,会降低高频大地电磁测深的探测能力。通过一维正演模拟,对无低阻覆盖层时探测能力和有低阻覆盖层时探测能力的变化进行对比研究,表明有低阻覆盖时,探测能力普遍下降;在相同探测能力条件下,高频大地电磁测深可探测异常体上覆围岩厚度与异常体厚度的比值,和围岩与异常体电阻率比值呈线性关系。通过綦江县某铁...     

气煤孔隙度与Thomsen各向异性参数高阶交错网格数值模拟

对气煤煤样孔隙率测试的基础上,讨论了各向异性系数与孔隙率的关系,表明各向异性的大小与孔隙率之间有一定的关系。应用高阶交错网格差分算法,对一阶速度-应力弹性波方程进行模拟,分析了其稳定性和边界条件。通过对Thomsen各向异性参数的模拟表明,主要控制的是P波的各向异性强弱,δ主要控制水平垂直正交方向上的各向异性强弱。对煤田各向异性研究具有重要意义。     

时间域拉东变换的实现及其改进

介绍了时间域拉东正反变换的基本原理及该项技术在地震数据处理中的应用。以Ng和Perz的算法为基础,以开发地震勘探处理与解释系统为实例,解决了实验中遇到的拉东正变换时收敛精度、反变换数据无损还原和截断信号引起的能量扩散、低信噪比情况下的滤波等问题,进一步研究了改进算法在其他地形下应用问题。