基于广义旋转法双九点格式标量波方程数值模拟

前人虽然基于传统旋转法提出了四阶精度最优17点差分格式,用于提高频率域地震波场数值模拟精度,但其仅适用于等网格间距的情形.这大大限制了该格式的使用范围.为了进一步提高17点有限差分格式的数值精度并将其推广到矩形网格,本文基于二阶精度有限差分算子利用广义旋转法提出了双九点格式,由于其差分格点与前人17点格式在分布上一样,所以也可称为二阶精度最优17点格式,但由差分格式构造原理上来讲,称其为双九点格式更妥.前人四阶精度17点格式仅为本文等网格间距情形时的特殊情况,本文方法单位波长网格点数仅需要2.2个即可.本文格式在继承传统旋转法良好几何旋转性质的同时,拥有平均导数方法适用于矩形网格的特点,和平均导数法所得的广义17点格式相比,本文格式数值精度更高,数值频散抑制性能和差分算子对称性更好.同时,本文双九点格式方法和思想对于后人借助传统九点格式的构造方法将其扩展到17点格式求解各类波动方程具有十分重要的意义.     

裂隙弱度参数的波场特征

本文将Schoenberg线性滑移裂隙介质理论应用到HTI(LSD)介质的弹性波波场数值模拟之中.主要针对描述裂隙状态的2个裂隙参数--法向弱度EN和切向弱度ET对波场的影响进行了有限元模拟.模拟结果揭示了EN主要影响快波而ET主要影响慢波的波场特征,通过EN和ET变化对应的裂隙含水程度说明了EN对裂隙中的流体敏感而E...     

起伏地表条件下各向异性地震波最短路径射线追踪

在地震波正反演研究中,考虑起伏地表和地震各向异性具有非常重要的理论意义和实际应用价值.本文在前人研究的基础上,将最短路径追踪算法引入到起伏地表各向异性介质模型的地震波走时计算中.模型剖分时,整体模型划分成正方形单元,起伏边界附近以不规则网格逼近,进而采用非规则节点布置实现非规则网格处的最短路径计算.追踪计算中采用Sena群速度近似公式,得到各向异性地震波的走时,实现了复杂地表情况下各向异性介质模型中地震波的射线追踪.理论模型计算结果显示,本文方法能够可靠地应用于复杂各向异性介质模型,具有较高的计算精度.     

东亚地区春季二氧化硫的输送与转化过程研究I.模式及其验证

结合新近评估的东亚地区污染源资料,作者利用一个耦合的区域化学输送模式系统以探讨东亚地区春季期间气象过程、气相与液相化学过程、非均相化学过程、气溶胶过程和干湿沉降过程对二氧化硫输送及转化过程的影响,并研究二氧化硫和硫酸盐气溶胶的空间分布及变化特征.模拟的二氧化硫和硫酸盐气溶胶的浓度值与2001年春季飞机和地面获取的观测值进行了比较.比较结果显示,模拟值与观测值具有很好的一致性,模式系统很好地反映了二氧化硫和硫酸盐气溶胶的分布特征和变化规律,再现了许多观测到的重要特征,为进一步分析模拟结果奠定了基础.     

关于“气煤弹性各向异性系数实验测试”一文的讨论

研究、确定气煤弹性各向异性的类型和大小,对煤田地震各向异性理论的研究具有重要意义."气煤弹性各向异性系数实验测试"一文通过对15块气煤岩样的超声测量,对气煤的弹性各向异性特征进行了分析,获得了对气煤地震勘探有指导意义的结论.但该文的一些理论概念和公式的使用及结论的获得有值得商榷的地方.本文针对此提出了一些修正意见,并利用该文的测试数据进行了补充分析,获得了气煤弹性各向异性相对完整的认识.     

最优化广义离散Shannon奇异核交错网格褶积微分算子地震波场模拟

笔者基于离散Shannon奇异核褶积微分算子(GDSCD)计算地震波速度应力方程的空间导数,推导了一阶GDSCD的具体形式,并提出了优化方法,即在频率域逼近平面波的真实导数,得到了不同半径和采样下限的最优窗函数系数,通过滤波响应分析算子精度,与多种数值方法对比以及模型测试表明,笔者构造的最优化GDSCD模拟地震波具有较高的计算效率和精度.     

大亚湾核电站厂址区域大气扩散研究

采用一个基于流体静力平衡的动力学诊断模式对大亚湾核电站厂址区域的气象场进行实际模拟,研究了复杂地形和下垫面对流场的动力和热力作用。数值实验表明：模式能很好地模拟中尺度局地环流的基本特征和变化规律。在模拟气象场的基础上,利用Monte Carlo多源模式计算了污染物的时空分布和干湿沉积量。模拟结果显示：模式很好地反映在海陆风和山谷风环流及排牙山影响下大气参数的非均匀性和非平稳性。     

地形性中尺度环流的数值模拟及验证

建立了一个高分辨三维边界层诊断模式以模拟下垫面热力和动力的 非均匀性对边界层结构的影响。模式模拟了广东大亚湾核电站厂址区域（12 km×12 km） 47个典型天气系统下的地形性中尺度环流,其大坑的模拟风向、风速和温度与观测值 相比,风向差小于17.5#+[o]的概率为87.5%,风速比值在0.7～1.3之间的概率为72.9%, 温度差小于0.4#+[o]C的概率为72.9%。     

Model analysis of seasonal variations in tropospheric ozone and carbon monoxide over East Asia

Temporal-spatial variations in tropospheric ozone concentrations over East Asia in the period from 1 January 2000 to 31 December 2004 were simulated by using the Models-3 Community Multi-scale Air Quality (CMAQ) modeling system with meteorological fields calculated by the Regional Atmospheric Modeling System (RAMS). The simulated concentrations of ozone and carbon monoxide were compared with ground level observations at two remote sites, Ryori (39.03°N, 141.82°E) and Yonagunijima (24.47°N, 123.02°E). The co...     

时间域全波场各向异性性参数反演

从各向异性弹性波的有限元正演方程出发，导出了反问题中时间域雅可比矩阵求解的计算公式。它具有与时间域有限元正演方程相同的表达形式，故可通过有限元正演计算来获得雅可比矩阵。研究了有限元正演算法的效率和精度、吸收边界条件等的问题，以提高反演系统的效率和精度。在此基础上，实现了叠前全波场各向异性弹性参数反演。计算表明，在初始模型偏离真实模型较大的情况下，层状模型和横向不均匀模型的反演结果均能准确地收敛到真实模型上。