复杂二维/三维大地电磁的有限单元法正演模拟策略

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

弹性波正演模拟中改进的非分裂式PML实现方法(英文)

在弹性波有限差分正演模拟中,吸收边界条件常用来吸收截断边界处引入的不期望边界反射,其中完全匹配层(PML)吸收边界条件被认为是目前最理想的吸收边界条件。但是PML吸收边界条件的传统实现却存在着很大不足:全局分裂式PML吸收边界条件实现简单但是需要占用太多内存;局部分裂式PML吸收边界条件需要考虑多个边界和角点区域,编程实现非常复杂;非分裂式PML吸收边界条件由于涉及卷积运算,计算量很大。本文基于非分裂式PML吸收边界条件,结合复频移伸展函数,提出了一种新的数值实现方法,其计算方程简单、占用内存小、编程实现容易,是对PML介质理论数值实现的改进和完善。     

弹性波正演模拟中PML吸收边界条件的改进

在弹性波有限差分正演模拟中, 完全匹配层(PML) 吸收边界条件是使用广泛、吸收效果最好的吸收边界条件.在目前的两种PML实现方法中, 分裂形式的完全匹配层(SPML) 方法计算存储量大、编程实现复杂; 非分裂形式的完全匹配层(NPML) 方法计算效率低、计算过程复杂.针对传统PML吸收边界条件在实现过程中存在的问题, 推导出了一种简洁有效的非卷积实现的NPML吸收边界条件, 既不需要对场分量进行分裂, 也不需要做复杂的卷积运算.分析结果表明, 本文实现的NPML吸收边界条件不仅具有良好的吸收衰减性能, 而且计算方程简单, 编程实现容易, 占有内存更小.      

天山南坡科其卡尔巴契冰川度日因子变化特征研究

度日模型是估算冰川消融的一种简单而有效的方法.根据科其卡尔巴契冰川2003年的观测资料,分析了该冰川度日因子的空间变化规律及其影响因素.研究表明:各高度上的度日因子,介于2.0～9.7mm·℃^-1·d^-1之间变化,平均值为5.7mm·℃^-1·d^-1,与青藏高原各冰川及其它地区冰川相比较小;随着海拔的增高,度日因子随之递增;随平均气温的升高而随之递减.由于冰面状况复杂,度日因子变化幅度较大,裸冰区的度日因子明显大于表碛覆盖区.人为测量误差、反照率、地形等对度日因子的影响也不容忽视.     

植物中单体烷烃碳同位素组成与其生长环境的关系

用色谱一同位素质谱议对青藏高原东北部高寒草甸生态系统中的两种高等模特一火绒草（Ｌｅｏｎｔｏｐｏｄｉｕｍ ｌｅｏｎｔｏｐｏｄｉｏｉｄｉｅｓ）和蒲公英（Ｔａｒａｘａｃｕｍ ｍｏｎｇｏｌｉｕｍ）的烷烃进行分析,讨论环境因素对植物单体烷烃碳同位素组成（δ＾１３Ｃｉａ）的影响。研究发现,两种植物奇碳数正构烷烃的δ＾１３Ｃｉａ值都比较小,且有随海拔高度的增加而增大的趋势,而偶碳数正构烷烃的δ＾１３Ｃｉａ值可能     

粗差验后方差的无偏估计与最优稳健估计

在正态粗差假设下导出了粗差验后方差的无偏估计,对误差工膨胀模型和误差均值移动模型,两者的无偏估计公式是相同的。这证明了李德仁验后方差的朱建军方差不是无偏的。由于偏方定义的彭方法是正态粗差假设下的最优稳健估计。     

地质雷达正演中的频散压制和吸收边界改进方法

从麦克斯韦方程组出发,建立了地质雷达的时域有限差分法（ＦＤＴＤ）数学模型,导出了理想 散关系和超级吸收边界条件。理想频散关系 考虑了ＦＤＴＤ法的收敛性和稳定性,也考虑了高频电磁波在Ｙｅｅ氏网格中的传播特点;超吸收边界条件则用磁场分量来提高电场分量精度。数值试验表明,理想频散关系能真实地反映雷达波在地下介质中的传播规律,超吸收边界条件能有效减小截断边界的伪反射,提高正演精度。将之应用于实际计算,取得     

石油烃降解细菌对正烷烃的降解作用研究

报道使用气相色谱法测定了石油烃降解 细菌对柴油的正烷烃的降解作用。结果表明,石油烃降解细菌对正烷烃有明显的降解作用, 混合菌株的降解率明显高于单菌株的降解率;在20℃的条件下,经过21d后,绝大部分的正 烷烃被降解。总的降解率为94.93%, 其中细菌的降解率为75.67%, 理化因子降解率为19.26% ;温度对正烷烃的降解率有明显的影响,在10℃时,其降解速度较慢,20℃降解加快,35℃ 温度使正烷烃的降解速度最快。     

弱各向异性介质弹性波的准各向同性近似正演模拟

准各向同性(QI)近似可用于弱各向异性介质的正演模拟.本文通过运用QI方法的零阶和一阶近似,计算了VTI介质模型的地震记录.得出的地震记录与标准各向同性射线理论(IRT)和基于伪谱法的三维地震正演模拟得出的地震记录作了比较,可以认为是精确的合成地震记录.