基于水热耦合的青藏高原分布式水文模型 ——Ⅱ.考虑冰川和冻土的尼洋河流域水循环过程模拟

尼洋河流域是雅鲁藏布江第四大支流,受冰川、积雪和冻土影响,水循环关系极其复杂。为深入研究该区域内的水文循环过程,本文在寒区水循环模型（WEP-COR）的基础上,针对青藏高原气候和地质特点,构建了耦合“积雪-土壤-砂砾石层”连续体和“积雪-冰川”水热过程模拟的青藏高原分布式水循环模型（WEP-QTP）。在尼洋河流域通过对2013—2016年的流量过程模拟发现,工布江达和泥曲站的逐月流量Nash-Sutcliffe效率系数分别达到0.810和0.752,比改进前的0.430和0.095有明显提升;以2015年为例,对比WEP-COR和WEP-QTP模型发现,WEP-QTP模型在汛期特别是主汛前（冻土融化期）模拟的流量过程不会出现较大的波动,模拟得到的逐日流量Nash-Sutcliffe效率系数相比WEP-COR从-0.67提高到0.54。模型增强了地下水含水层的调节作用,使得流量过程更加平稳且接近实测,研究结果表明,WEP-QTP模型适用于青藏高原的水文模拟。     

利用陆态网络基准站资料研究山西断陷带近期地壳形变特征

利用陆态网络GPS基准站资料计算了山西断陷带内部以及相对于周边地块的基线时间序列,并基于均匀应变模型计算了山西北部至晋冀交界和山西南部至晋豫交界的应变时间序列,结果表明:(1)陆态网络山西基准站观测数据质量较高,能够在一定程度上反映山西断陷带的地壳微动态变化特征。(2)山西北部至晋冀交界区域第一剪应变和最小主应变自2012年10月起减小,山西大同至河北张家口和山西大同至北京延庆的基线同期出现缩短,可能表明该区域北东向的右旋剪切和压性应变在逐渐增强。(3)山西北部至晋冀交界区域张性剪切应变2014年5月起出现的短期增强与涿鹿4.3级地震震源机制较为一致,可能与此次地震有关,该短期变化在震后逐渐恢复。     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

Analysis on the noises from continuously monitoring GPS sites

In this paper,we analyze the time series of site coordinates of 27 continuously monitoring GPS sites covered bythe Crustal Movement Observation Network of China over the whole country.The data are obtained in the periodfrom the beginning of the observation to the November of 2005.On the basis of data processing,we analyze thepower spectrum density of coordinate component noise at each site and calculate the spectral indexes manifestingthe noise property of each component.The spectral indexes indicate that for most sites,the noise of time series ofeach coordinate component can be addressed by the model of white noise flicker noise;and for a small amountof sites,it can be described by the model of white noise flicker noise random walk noise.We also quantita-tively estimate each noise component in the model by using the criterion of maximum likelihood estimation.Theresult shows that the white noise in the time series of GPS site coordinates does not constitute the main part ofnoise.Therefore,the error estimation of site movement parameters is usually too small,or too optimistic if weconsider the white noise only.Correspondingly,if this factor is not fully considered in explaining these movementparameters,it might mislead the readers.     

S－闭空间上的函数

研究了当Ｘ为Ｓ－闭空间时，Ｃ（Ｘ）的子集Ｆ为紧的充要条件，从而扩充了著名的Ａｓｃｏｌｉ定理。     

北京市地理空间框架基准体系建设与应用

重点介绍了北京市地理空间框架基准体系的实现过程,包括空间坐标基准、高程基准、似大地水准面和全球卫星定位测绘服务系统建设。基准体系检测结果表明,北京市地理空间框架基准体系达到设计要求,在国内率先实现了测绘基准的现代化建设。     

基于向量式有限元的黏滞阻尼减震结构抗竖向连续倒塌动力反应分析

向量式有限元是以向量力学为基础的一种新的结构分析方法,在处理结构大变形等复杂行为时具有较大的优势。基于向量式有限元理论建立了黏滞阻尼单元,对附加黏滞阻尼器的平面钢框架结构进行了抗竖向连续倒塌动力分析,结合拆除构件法,采用MATLAB编制可以考虑初始变形的瞬时卸载法程序,实现结构在构件拆除前的静力分析和构件拆除后动力分析的全过程统一。研究了阻尼器布置位置和参数在结构抗竖向连续倒塌中的性能需求,以失效点竖向位移时程曲线、梁端转角、动力放大系数和结构塑性铰分布为参考指标,对比分析布置阻尼器前后钢框架结构的抗连续倒塌能力。结果表明向量式有限元是一种研究结构竖向连续倒塌动力响应的有效方法,合理布置阻尼器能够有效控制剩余结构的变形和振动,降低构件内力,减少塑性铰个数,较大地提高结构的抗竖向连续倒塌能力。     

2008年汶川8.0级地震孕震机理研究

用三维流变非连续变形与有限元相结合（DDA+FEM）的方法,在青藏高原及邻近地区三维构造块体相互制约的大背景中,考虑了龙门山断裂带东西两侧地势、地壳厚度和分层的明显变化,及断裂带东侧四川盆地及鄂尔多斯块体坚硬地壳阻挡的影响,通过用GPS资料做位移速率边界约束和震源机制约束,计算得到研究区的速度场和应力场与该地区GPS测量结果和震源机制分布结果基本一致.在此基础上,模拟计算现今构造块体边界断层上表征剪应力及法向应力等综合影响的危险度分布.结果表明,上、中地壳层危险度分布中危险度较高的地段多数与近几十年来发生的七级以上大震区域基本一致.包括2008年汶川8.0级等大震的发震断层.通过分别对龙门山断裂带东西两侧的两种不同构造格局进行试算表明,龙门山断裂带东西两侧地势、地壳厚度、分层与物性明显变化对汶川大震的孕育发生均起了关键性作用.计算得到的应变率强度分布图可见,高原东部整个边缘地带均接近应变率强度的陡变带.其中以龙门山断裂带上的陡变最为明显,西侧应变率强度是东侧的近4倍,而且断裂带东侧应变率强度等值线衰减比西侧快.反映了汶川大震逆冲型发震断层地区独特的特征.此外,由计算得到的应变能密度分布图可见,龙门山断裂带在上、中地壳层中均位于宽度相同、其走向与龙门山断裂带走向一致的高应变能密度带中,在上地壳层这个带的东西两侧则是应变能密度较低的地区,而在中地壳层,其强度在断裂带东侧逐渐向东衰减,西侧应变能密度高,而东侧应变能密度较低.表明在印度板块强烈推挤作用和高原各构造块体相互制约及龙门山断裂带东西两侧特殊构造环境中,高原地壳物质向东水平运动,受到龙门山断裂带东侧介质刚性强度较大的四川盆地阻挡,使得汶川大震发震断层在大震前已积累了相当水平的应变能,并同时处于力学上的不稳定状态.     

连续弯道水流模拟中二次流修正效果评价

为了选择适用于复杂连续弯道的二次流修正方法,选取了线性方法中两种典型的计算模式,通过在正交曲线坐标系二维浅水动量方程中增加扩散应力项开发了考虑二次流影响的平面二维水流模型。基于不同复杂度的4个连续弯道试验的水位和流速分布资料,通过理论分析和数值试验对比定量评估了传统二维模型和两种二次流修正方法在连续弯道水流模拟中的效果。测试结果表明,Delft3D模型的二次流修正方法自由度较高,适用于不同复杂度的连续弯道水流模拟,而Lien模型二次流修正方法适用于微弯或中弯,应用于中弯时需要慎重选用,不适用于急弯连续弯道模拟。两种方法比较,建议连续弯道水流模拟中优选Delft3D模型的二次流修正方法。     

不连续面在双重介质热-水-力三维耦合分析中的有限元数值实现

在建立双重介质热-水-力耦合微分控制方程的基础上,提出了裂隙岩体热-水-力耦合的三维力学模型,对不同介质分别建立以节点位移、水压力和温度为求解量的三维有限元格式,开发了双重介质热-水-力耦合分析的的三维有限元计算程序,在有限元数值分析中不连续面应力计算采用等厚度空间8节点节理单元进行离散,而不连续面渗流和热能计算时采用平面4节点等参单元进行离散,这样保证了不同介质之间的水量、热量交换和两类模型接触处节点水头、温度和位移相等。通过高温岩体地热开发算例,揭示了在热-水-力耦合作用下不连续面处于低应力区,其张开度随运行时间的延长呈非线性增加,非稳定渗流阶段不连续面显著地控制着渗流场的整体分布,它的水头远高于拟连续岩体介质的水头,而进入稳定渗流阶段不连续面的控渗作用不明显,由于高温岩体地热开发系统中存在大规模的热量补给,不连续面对岩体温度场分布的影响并不显著。