山西奇村地热田资源评价

山西奇村地热田开采已有30年历史，从1988年大量开采地热水以来，使浅层热水生产井水温持续下降。本文建立了奇村地热田的概念模型，利用LUMPFIT软件，采用热田多孔抽水试验资料，对地热田在不同热水开采量下，地下水位降深及温度变化进行了预测，提出了合理的开采量，对地热田开采有重要指导意义。     

不同电阻率测深方法对旁侧不均匀体的反映

利用三维有限元数值模拟,重点研究了旁侧三维不均匀体所产生的视电阻率的影响问题,展示了4种常规电阻率测深法在旁侧不均匀体上的异常分布特点;探讨了测深点正下方目标体异常与旁侧不均匀体异常的区别问题,指出旁侧不均匀体的影响是产生电阻率法假异常的主要原因之一.     

锚杆动测问题的数值解

建立了锚杆动测问题的有限元模型，并对锚杆中的缺陷进行了数值模拟，获得了位移场的数值解。对该解在时域和频域范围内的讨论表明，锚杆中的缺陷将引起动测曲线一系列的变化，为现场动测曲线的识别提供了依据。     

深厚软弱地基上相邻桩箱基础高层建筑地震反应的数值模拟

采用二维有限元整体分析法对深厚软弱地基上相邻的双幢桩箱基础高层建筑的地震反应进行了数值模拟，探讨了SSI效应对双幢高层建筑地震反应的影响规律；并分析了深厚软弱地基上双幢高层建筑和单幢高层建筑地震反应之间的差异。一般地，SSI效应使深厚软弱地基上相邻的双幢高层建筑各楼层的绝对加速度反应要比单幢高层建筑的反应略小一些，但也可能使双幢高层建筑中部楼层的绝对加速度反应大于单幢高层建筑的反应。与此相反，SSI效应一般使相邻的双幢高层建筑各楼层的相对位移反应有所增大；但也可能使双幢高层建筑各楼层的相对位移反应有所减小，两者的差异与场地条件、输入地震动特性有关。此外，还探讨了双幢高层建筑的间距对其地震反应的影响。随着间距的增大，双幢高层建筑的地震反应趋近于单幢高层建筑的反应。输入加速度峰值越大，场地土越软，相邻的双幢高层建筑的间距对其地震反应的影响越大。     

1995年永登地震山体变形计算机模拟与机理分析

1995年甘肃永登5．8级地震引起了黄土山梁的酥裂震陷，说明黄土在较小震级下也会有失稳现象．利用三维土体的有限元计算方法，对永登地震山体变形进行了计算机模拟，并定量地分析了变形的机理。     

长江三角洲稻田生态系统综合增温潜势源汇交替的数值分析

将生物地球化学C, N循环模式与植被生态模式耦合, 建立了生物地球化学-植被生态耦合动力模式. 验证表明模型可以准确模拟稻田生态系统中作物生长及CH₄, N₂O排放过程. 利用模式分析了不同施肥量对CO₂吸收及CH₄, N₂O排放的影响, 以及3种温室气体综合增温潜势的演变规律, 发现了生长期内稻田综合增温潜势源汇交替现象, 提出最大汇施肥量以及零排放施肥量等概念. 利用模式研制了长江三角洲地区不同施肥方案, 为该地区可持续发展提供了科学依据.     

末次盛冰期东亚气候的数值模拟

用美国国家大气研究中心(NCAR)的CCM3全球气候模式研究了末次盛冰期(LGM)和现代情景下的东亚季风和地面水分特征以及青藏高原冰川扩张. 结果表明: 在LGM时, 我国北方和西太平洋地区冬季风显著加强, 南方地区冬季风变化不大; 而对LGM时期的夏季风, 我国南方和南海地区显著减弱, 北方变化不显著; LGM时期季风的这种变化, 使我国东北、华北大部分地区、黄土高原和青藏高原东部年降水量比现代显著减少, 造成这些地区当时地面净失去更多水分, 使当地变干燥, 其中青藏高原东部、黄土高原西部地面变干燥最显著; 而在LGM时期青藏高原中部一些地区由于蒸发减少使地面变湿润, 有利于当时这些地区的湖面上升; 此外, LGM时期冬季青藏高原 绝大部分地区积雪明显比现代厚, 通过分析模拟资料计算的冰川平衡线高度发现: 尽管我们模拟出LGM时期较小的降温幅度, 但是通过模式中大气物理过程青藏高原降水和气温之间保持平衡, LGM时期当地冰川平衡线高度与现代相比降低了300~900 m, 即从现代的5400 m以上降为4600~5200 m, 指示着LGM时期青藏高原冰川的大规模扩张.     

具有脆性形变和蠕动变形的复合单断层系模型中的应力、滑动与地震

对弹性半空间中的几个含有不连续断层系的三维弹性连续介质模型，在沿单元垂直走滑断层发生的滑动过程进行了数值模拟，模型的几何形状和所使用的边界条件大致与圣安德烈斯断层中部的情形相符合，数值模拟综合考虑了脆性形变和蠕动变形，断层带的总形变率为蠕变和摩擦滑动率之和。脆性断层的性质由破裂段（数值单元）的不同地震应力降的分布来给定。所假设的分布描述了两种理想化的情形，它们分别对应于不同断层滑动的极限状态：（1）以大尺度范围为特征的强无序状态，表示不成熟的断层带及扩展空间域。（2）相对有序的状态，以小尺度范围为特征，表示成熟的完全滑动断层。假定断层的蠕变特性在所有的情形中都相同。并且是用蠕滑速率与应力的幂次律关系来刻画的。利用这种脆性蠕变过程和模型参数可以得到一个在12．5km深度含有“脆性－延性”过渡带的应力－深度剖面，以及一个在1857年地震破裂西北65km处存在有“脆性－蠕变”过渡带的沿走向的应力剖面。模拟的震源空间分布在统计意义上与观测数据相类似，结果表明，断层强烈不均匀性的特征尺度范围会在断层系统的地球地震响应中显著表现出来，若尺度范围小，则会使震级－频率分布接近于特征地震分布，在时间上会层致如地震空区说所预言的大地震的准周期分布，而大尺度范围的地震活动性符合古登堡－里克特的震级－频度关系以及大地震的随机或是成簇的时间分布特征。数值模拟的结果表明，各种形式的震级－频度关系以及地震时间分布的统计特征都可以通过表征断层带不均匀性的尺度范围统一起来，这可以通过借助一个给定的断层带或是岩石层的结构特征给出一个清晰的物理解释，并且得到地震及断层观测资料的支持。在一些模拟中，小地震的震级－频度关系曲线明显低于自的古登堡－里克特关系曲线。这说明在构造加载过程中，小地震通过平滑应力的长波成分，为大地震的发生创造了条件。这一过程是通过与小地震的大量断裂有关的短波应力的消失过程来完成的，大地震的大尺度破裂的发生正好导致与上述长，短波应力的平滑与消失变化相反的过程。     

地震学中的统计学原理（二）

3.4 模型检验模型检验即检查选用模型的特征与实际样本特征是否一致。选择模型与检验模型采用不相交的数据集较为理想。在做不到这一点的地方,检验需要考虑到模型已被调整到尽可能与数据密切拟合的事实,至少在某种意义上如此。     

哈尔滨市磨盘山水利枢纽区三维渗流数值模拟模型研究

哈尔滨市磨盘山水库防渗型式拟采用垂直砼防渗墙及帷幕灌浆方案。为了分析左岸不同长度防渗对水库渗漏量及坝基渗透稳定性的影响。对磨盘山水利枢纽区渗流场进行了三维数值模拟研究。结果表明，随着左岸防渗长度的增大。渗漏量及下游渗流出口渗透比降基本呈线性趋势减小。当左岸防渗长度为200m时。渗漏量及各渗透比降均小于允许值，为推荐的渗控方案。