一个三维斜压水动力模型的建立及其在日本琵琶湖中的应用

水温的空间分布对很多湖泊水文现象关系甚大,本文建立了日本琵琶湖的三维斜压水动力数值模式。模拟计算结果表明：１）在不同风速驱动下,湖的上、下两层可形成不同的环流型;２）除形成环流处外,其它部分流速、流向在上、下层表现出明显的垂直切变;３）闭合式温度场储存位能,即使在无风和初始静止流场情况下,也可驱动湖流并形成若干环流特征;４）风对水的驱动作用是形成北湖北部稳定环流的决定性因子,并且本身可构造出温度场     

崇文门站施工允许地表沉降及三维有限元模拟分析

建设中的北京地铁五号线崇文门站下穿既有地铁一号线区间隧道,为保证既有地下铁道正常运营和地下结构的安全,需严格控制新建车站施工引起的地层位移。结合新建车站暗挖施工及既有线地铁的实际情况,分析和计算穿越既有隧道施工地表允许沉降值。为了控制地表沉降,确保施工安全,拟采用小刚度管棚或大刚度顶管作超前预支护,用3D-S igm a三维有限元软件进行施工效应的计算模拟,掌握大管棚和顶管预支护洞室的力学效应,预测车站施工引起既有隧道的沉降量。最后分析比较两种预支护的支护效果,为选择合理的预支护方案提供参考。     

华北克拉通地壳三维密度结构与地质含义

利用重力异常揭示地壳三维密度结构是重力资料地质解释的目标和任务,密度反演的好坏至关重要.本文对华北克拉通进行了重力异常多尺度密度反演研究,首先利用小波变换对重力异常进行多尺度分解,接着利用功率谱分析方法估算各层场源的平均深度,然后利用广义密度反演方法进行各层密度反演,对华北地区进行了多层密度反演,得到其密度结构,并进行地质解释和油气藏分析.结果显示了重力场多尺度密度反演方法的有效性,对华北地区密度结构的地质含义进行了初步分析,位于造山带中的低密度异常主要反映沿造山带展布的花岗岩体.对比华北东部区浅层密度扰动与油气坳陷位置,发现油气坳陷都表现为低幅度的低密度异常区.说明利用小波多尺度反演提取的密度信息对油气勘探的部署有一定指导意义.     

新疆布腊图勘查区地震磁法勘探煤层火烧边界

煤层在燃烧过程中,其围岩受到高温烘烤,致使煤层上、下部围岩裂隙发育,失去了形成层状地震反射的条件,同时在磁性、导电性等方面物理性质也发生显著变化,为地震勘探及磁法勘探的实施提供了理论依据。该文通过实例介绍了利用三维地震勘探和磁法勘探相结合准确圈定各主采煤层火烧区边界的方法及经验。     

斜入射平面SH波在层状半空间中沉积谷地周围的三维散射

采用直接刚度法计算自由场动力响应,以层状半空间中移动斜线均布荷载动力格林函数模拟散射波场,采用间接边界元方法求解了层状半空间中沉积谷地对斜入射平面SH波的三维散射问题.由于文中采用的层状场地三维动力刚度矩阵是精确的,且用于模拟散射波场的均布移动斜线荷载可以直接施加在沉积交界面处而不存在奇异性,所以本文方法具有很高精度.文中以均匀半空间和基岩上单一土层中沉积谷地对入射SH波的散射为例进行了数值计算,研究表明,沉积谷地对地震波的三维散射与二维散射之间存在本质差别;层状半空间中沉积谷地与均匀半空间中沉积谷地附近地表位移存在显著差异.     

基于非结构网格的电阻率三维带地形反演

地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探.     

黄河三角洲海岸冲淤及泥沙输运模型数值模拟分析

黄河三角洲是陆海交界地带,陆地海洋相互作用显著,泥沙所引起的岸滩演变、港口航道的淤积、水动力环境的改变等问题比较复杂。该文采用三维HEM-3D数值模型对黄河三角洲海域流场变化、盐度、悬浮泥沙浓度及海岸冲淤分布进行了潮流周期内的数值模拟分析。结果表明,黄河三角洲海域的悬浮泥沙浓度分布与潮流场变化和河口泥沙输入有密切的关系,在三角洲北部受五号桩外强潮流区的影响,近岸海底的泥沙发生明显的再悬浮,并在涨潮流向南输送,含沙量达1.5 g/L左右。在三角洲南部(现行河口区域和莱州湾区域),受现行河口入海泥沙扩散的影响显著。海域年冲淤分布,在北部废弃三角洲区域,由于海洋动力作用强烈,浅水冲刷,海底侵蚀显著,形成了明显呈沿岸展布的侵蚀中心,在侵蚀中心以外,侵蚀快速减弱。在现行河口区域,以淤积为主,在羽状流扩散的控制下,泥沙入海后向南输运,至莱州湾区域逐渐减弱。数值模拟的结果与卫星遥感解译的岸线变化基本一致,效果良好。     

基于三台测风激光雷达的大气湍流和三维风场研究

为了获取大气湍流和空间三维风场结构,利用3台同型号的测风激光雷达开展协同观测试验。（1）利用虚拟铁塔协同观测技术开展大气湍流探测,与香河102 m铁塔安装的三维超声风速仪观测结果做对比,32 m处高频（10 Hz）风速的相关系数高达0.92,平均误差为0.77 m/s,均方根误差为0.41 m/s;大气湍流强度（TKE）的相关系数高达0.99,平均误差为?0.02 m2/s2,均方根误差为0.08 m2/s2,并且协同观测的高频风速与三维超声风速仪的观测结果具有相同的频谱结构。（2）利用扫描协同观测技术开展三维风场探测,与铁塔上的常规测风设备相比,其90 m高度处的水平风速和风向的相关系数分别为0.92和0.93,平均误差为?0.41 m/s和0°,均方根误差为0.73 m/s和34°。相比于单台测风激光雷达,基于3台测风激光雷达协同观测技术具有一定的优势:不需要风场水平均匀的假设、探测精度更高等。但其对观测环境的要求较高:观测路径上不能有遮挡、观测必须协同等。在科研业务应用中,需要根据实际的观测需求合理制定观测方案。      

瞬变电磁Crank-Nicolson FDTD三维正演

时域有限差分(FDTD)方法使用Yee网格剖分电磁场的空间采样,通过时间步迭代实现电磁场数值模拟,具有内存消耗低、计算简单等特点,常用于瞬变电磁三维正演.然而,常规FDTD方法的时间迭代步长Δt受Courant-Friedrich-Lewy(CFL)条件严格限制,过多的迭代次数以及过密的采样往往导致计算速度慢、累积误差不断增大.本文提出一种不受CFL条件约束的无条件稳定隐式差分算法Crank-Nicolson FDTD(CN-FDTD)用于瞬变电磁三维正演.基于Crank-Nicolson差分方法对Maxwell方程组重新离散,空间网格仍然采用Yee元胞,时间步进采用在整时间步电场、磁场同时采样的策略,建立无条件稳定FDTD格式,突破CFL条件限制.与常规FDTD交替采样相比,CN-FDTD电场、磁场同时采样的策略构成的隐式差分格式,需要求解大型稀疏矩阵方程组.通常,瞬变电磁三维正演模型中产生的矩阵阶数往往较大,需要占用大量内存和求解时间.为解决上述问题,采用Crank-Nicolson-cycle-sweep-uniform(CNCSU-FDTD)方法近似求解CN-FDTD方程,在保...     

基于雷声到达时间差的单站闪电通道三维定位系统

利用麦克风阵列采集雷声信号,设计了一套由麦克风阵列和便携式数据采集存储设备组成的单站闪电通道三维定位系统,采用CCF相关函数法计算了雷声脉冲信号到达不同麦克风的时间差,通过最小二乘法获得声源的方向角和仰角信息,并结合声源与麦克风阵列的距离定位声源的三维位置。对一次包括8次回击的人工触发闪电过程进行了观测,得到多个声源点的三维定位信息。系统的定位结果与高速摄像的观测结果具有较好的一致性,定位方法合理、可行,为进一步研发便携式闪电通道三维定位系统奠定了基础。