半拉格朗日平方守恒计算格式的研究

作者研究了一种新的显式平方守恒的半拉格朗日计算格式.该格式是欧拉空间平方守恒格式在半拉格朗日空间的推广和发展,并采用了守恒插值方法.此外,还给出这种格式在一维原始方程上的应用.     

黄土完全软化强度与残余强度的对比试验研究

为了解决复活蠕滑型黄土滑坡强度参数的取值问题,开展完全软化强度与残余强度的对比试验研究。以山西地区典型黄土为研究对象,采用预压固结法制备饱和重塑试样,并进行反复直剪强度试验获取完全软化强度和残余强度参数。试验结果表明：黄土的完全软化强度以黏聚力为零、颗粒未发生定向排列为主要特征。完全软化强度与二次固结应力和黏粒含量有关。试样在二次固结应力小于300kPa时的应力-应变曲线呈应变软化型,完全软化强度大于残余强度,黏粒含量高的试样应变软化更显著; 试样在二次固结应力大于等于300kPa时的应力应变曲线呈理想塑性型,完全软化强度近似等于残余强度。完全软化内摩擦角与残余内摩擦角的差值和预压固结应力、二次固结应力及黏粒含量有关。完全软化内摩擦角与残余内摩擦角的差值随二次固结应力的增大而减小,最终趋于0。当预压固结应力小于300kPa时,内摩擦角差值及黏粒含量对内摩擦角差值的影响随着预压固结应力的减小呈乘幂性增大; 当预压固结应力大于等于300kPa时,完全软化强度近似等于残余强度,可用完全软化强度近似代替残余强度。研究结论为复活蠕滑型黄土滑坡稳定性分析时强度参数的取值提供了一定的参考。     

黏性声波方程数值正演中的匹配Z变换完全匹配层吸收边界

本文将非分裂场形式的匹配Z变换PML引入黏声波方程数值正演中,用时域有限差分模拟检验了其在黏声波方程模拟中的效果.数值正演结果表明,在大角度入射时匹配Z变换完全匹配层比传统PML表现更佳,消除了大角度入射产生的低频虚假反射.长时间能量衰减计算证明匹配Z变换完全匹配层在黏声波方程模拟中具有105时间步的稳定性.     

基于完全竞争市场模式的中国页岩气资源管理体制分析与思考

页岩气是化石能源中相对清洁高效的一种天然气.美国的页岩气革命可以算是本世纪最具影响力的能源革命,使其实现了梦寐以求的“能源独立”.从美国的经验来看,页岩气革命可以认为是开放市场和解放思想的成果,是完全竞争市场模式理论在页岩气资源开发中的成功实践.中国页岩气技术可采储量为世界第一,但由于中国页岩气管理体制等原因的制约,目前中国只在个别区块取得了勘探突破,离页岩气的大规模商业开采还有比较长的路.通过对页岩气资源管理体制深入分析,总结和讨论传统油气管理体制的不足以及美国页岩气开发给我国页岩气资源开发的启示,运用西方经济学经典的完全竞争市场模式理论对页岩气资源管理体制改革及相关政策提出适当建议.     

不分裂卷积完全匹配层与旋转交错网格有限差分在孔隙弹性介质模拟中的应用

完全匹配层吸收边界在地震波模拟中已广泛使用,但常用的场分裂格式完全匹配层吸收边界(SPML)和传统的不分裂完全匹配层吸收边界(NPML)对极低频入射波或大角度入射波的边界吸收效果不好.一种无需分裂和显式卷积计算的完全匹配层吸收边界(CPML)不仅能够解决常规PML吸收边界的不足,而且具有存储量小、计算效率高、易于编程实现的特点.本文将这种完全匹配层(CPML)吸收边界引入到孔隙弹性介质速度-应力格式的旋转交错网格有限差分算法中,对完全匹配层吸收边界参数进行数值分析,得到一组优化的参数.孔隙弹性介质数值模拟结果表明这种不分裂卷积完全匹配层的吸收效果优于常规完全匹配层.     

高阶交错网格和PML吸收边界在横向各向同性介质地震波场模拟中的应用

完全匹配层吸收边界条件通常可以很好地吸收模型边界的地震反射波,但对于横向各向同性介质的模拟效果欠佳,且尚在发展阶段。为此,文中推导了横向各向同性介质中弹性动力学波动方程,给出了施加完全匹配层(PML)吸收边界条件下时间域二阶、空间域十阶精度的高阶交错网格的有限差分形式,并分别建立了均匀的垂直向对称轴的横向各向同性介质(VTI介质)和倾斜向对称轴的横向各向同性介质(TTI介质)模型。计算结果表明,对于对称轴为任意角度的横向各向同性介质,当PML边界层厚度达到一定的数值时,可以很好地抑制人工边界所产生的地震波反射效应,且PML的吸收效果不会被入射角与入射波频率影响。     

探地雷达FDTD数值模拟中不分裂卷积完全匹配层对倏逝波的吸收效果研究

介绍了CPML边界条件的原理,推导了CPML的GPR正演FDTD差分公式,对比分析了Berenger PML、UPML、CPML三种PML对倏逝波的吸收性能.开展了PML边界中关键参数κ和α的选取实验,确定了参数的取值范围与选取原则.然后,以二维TM波为例,研究了倏逝波产生的机理,分析了决定逝波性吸收性能的影响因素.均匀介质的波场快照、检测点的反射误差及全局反射误差对比,说明了3种边界条件对传输波都具有较好的吸收能力,而对低频倏逝波的吸收表现迥异,其中CPML因为引入了参数α,对倏逝波的吸收效果最佳,但离散化造成的全域误差也最大.最后,应用加载UPML和CPML边界条件的FDTD程序,开展了GPR二维剖面法、宽角法矩状地电模型及三维复杂模型的正演,展示了倏逝波反射对雷达正演剖面及波场快照的影响.进一步对比了UPML与CPML对倏逝波的吸收表现优劣,结果显示,CPML可有效减少边界反射误差,并能取得满意的精度,综合考虑对倏逝波的吸收、全域误差、编程难易程度等因素,在GPR正演中推荐使用CPML.     

瑞利面波数值模拟中的PML吸收边界条件

建立了弹性介质情况下完全匹配层(PML)吸收边界的2×12阶速度-应力交错网格有限差分算法,讨论了PML吸收边界条件的构建及其有限差分算法实现.通过与未加吸收边界及加常规指数衰减吸收边界3种情况下比较的波场模拟计算表明, PML吸收边界具有吸收更干净且能够吸收各种角度的边界反射等优点,其吸收率(吸收能量与未吸收能量之比)达到99.99%,很好地消除了周期折叠效应,使得所要计算的波场特征变得非常清晰,瑞利面波清楚地显示在波形记录上.     

面向拓扑分析的海洋流场临界点提取算法研究

临界点是海洋流场拓扑结构中的重要构成要素,基于临界点的特征提取对于揭示海洋流场拓扑特征、开展海洋流场拓扑分析具有重要意义。本文基于临界点理论和Sperner引理,综合改进后的双线性插值算法和Sperner完全标号法,对海洋流场数据进行了临界点特征提取。首先,在双线性插值算法中添加滑动窗口处理,筛选临界点的候选网格单元,并采用聚合思想通过降低网格分辨率解决了网格插值中的二义性问题,同时考虑了0值网格存在的9种情形,通过迭代聚合思想滑动筛选候选网格单元,解决了插值网格均为0的情况。其次,提出了基于Sperner完全标号的最小值法临界点提取规则,将速度向量模最小的网格中心作为临界点,解决了实际流场物理场景中非0值的临界点提取。对两次提取结果进行合并、去重等处理,可以得到较为全面的临界点提取与分类结果。最后,通过对多个海域、不同深度流场数据的实验结果分析,证明了综合后的临界点提取方法的有效性及可行性。     

不完全阶次球谐位系数模型的研究

利用局部重力资料确定不完全阶次位系数模型的方法进行研究．首先根据Molodensky边值问题的数学理论和函数展为级数式的一般原理,论述了方法在理论上的正确性;然后根据模型结构,提出了位系数补偿因子,并给出其确定的原则及计算公式;接着构制一个适合局部区域的不完全阶次的高阶位系数模型的方法进行探讨;最后对其精度进行检验．结果显示局部区域2160阶不完全阶次模型对实测垂线偏差二分量的恢复精度分别为±2．56s和±2．61s．