2014年云南景谷6.6级地震前后活动能量场分析

本文将地震释放能量作为地震活动场的基本要素, 采用自然正交函数展开方法, 系统研究2014年景谷6.6级地震震中区域地震活动能量场的时空特征。 结果表明, 景谷6.6级地震前能量场的前4个主要特征值对应的典型能量场时间因子都出现明显上升变或下降异常, 空间等值线形成局部能量高值异常危险区, 可能是景谷 6.6级地震空间异常的中短期特征。 我们还进一步研究了景谷6.6级地震前后空间异常随时间的演变, 给出了异常发展变化模式, 最后应用研究地震资料的累积次数从物理角度解释异常形成的机制。     

2015年尼泊尔8.1级地震前后拉萨地电场观测数据变化分析*

2015年4月25日尼泊尔8.1级强烈地震前后,拉萨地磁台(以下简称“拉萨台”)地电场观测数据出现了长时间、大幅度异常变化。在对拉萨台地电场观测条件、观测系统以及观测结果进行初步讨论的基础上,利用“能量合成累加”和“MSA功率谱分析”等数据处理方法,从“时间域”和“频率域”两个方面,对该异常变化的主要特征和演化过程进行了分析。结果显示,尼泊尔8.1级地震前后,拉萨台的地电场观测数据异常变化,经历了“趋势变化-扰动变化-发震期-恢复期”等发展阶段,以及“低频变化-高频变化-平稳变化-高频变化”等演化过程。综合分析认为,拉萨台地电场观测异常变化特征和演化过程,与相关现象和机理研究结果基本吻合,能够完整反映强烈地震孕育、发生和发展全过程,是具有一定客观性、典型性和代表性的珍贵资料。     

等效水力隙宽和水力梯度对岩体裂隙网络非线性渗流特性的影响

等效水力隙宽和水力梯度是影响岩体裂隙网络渗流特性的重要因素。制作裂隙网络试验模型,建立高精度渗流试验系统;求解纳维-斯托克斯方程,模拟流体在裂隙网络内的流动状态,研究等效水力隙宽和水力梯度对非线性渗流特性的影响。结果表明,当水力梯度较小时,等效渗透系数保持恒定的常数,流体流动属于达西流动区域,流量与压力具有线性关系,可采用立方定律计算流体流动;当水力梯度较大时,等效渗透系数随着水力梯度的增加而急剧减少,流体流动进入强惯性效应流动区域,流量与压力具有强烈的非线性关系,可采用Forchheimer方程计算流体流动。随着等效水力隙宽的增加,区别线性和非线性流动区域的临界水力梯度呈幂函数关系递减。当水力梯度小于临界水力梯度时,控制方程可选立方定律;当水力梯度大于临界水力梯度时,控制方程可选Forchheimer方程,其参数A和B可根据经验公式计算得到。其研究结果可为临界水力梯度的确定及流体流动控制方程的选取提供借鉴意义。     

基于逐减随机震源采样法的频率域二维黏滞声波方程全波形反演

全波形反演方法利用叠前地震波场的运动学和动力学信息重建地下速度结构,具有揭示复杂地质背景下构造与岩性细节信息的潜力。然而,巨大的计算量是阻碍其发展的一个瓶颈问题。为此,研究者们提出了震源编码技术来减少计算量,但是此方法在模型更新过程中会引进随机串扰噪声,降低反演结果准确性。所以,在保证计算精度的情况下,本文提出了采用逐减随机震源采样的方法来高效计算全波形反演问题。笔者将此方法应用于频率域二维黏滞声波波动方程全波形反演,开始了在频率域进行随机震源采样类方法的研究,计算过程中共使用了依次增大的8个频率段;并应用Overthrust模型来验证此类随机震源采样法的正确性。实验结果表明：基于逐减随机震源采样法的反演结果与实际Overthrust模型的拟合误差为0.065 65,而应用基于全部震源的全波形反演方法得到的反演结果与实际Overthrust模型的拟合误差为0.064 64,两者差别不大;但计算用时由740 min减少到291.2 min,即计算效率提高了2.54倍。为了更好地确定方法的有效性,将其应用于Marmousi模型进行试算。模型试算结果表明：基于逐减随机震源和基于全部震源得到的反演结果与实际Marmousi模型的拟合误差分别为0.080 12和0.078 97,相差不大;但计算用时由1 218.9 min减少到274.4 min,计算效率提高了4.44倍。综上,在保证反演精度的情况下,基于逐减随机震源采样法的频率域全波形反演方法大大减少了计算量,具有不可替代的计算优势,并且没有引进随机串扰噪声。     

北京市房价与地价的动态关系——基于结构方程模型的实证分析

城市房价与地价之间的关系错综复杂,不仅受多种因素的交织影响,相互之间也存在动态关系。研究房价与地价关系的传统方法,如Granger因果检验和回归分析等,无法刻画房价与地价之间多维的网络状关系,相比之下,结构方程模型能同时处理多个内生潜变量,且不受观测指标共线性的影响,为刻画地价与房价的交互作用提供了新的工具。从住房与土地市场的供需传导机制出发,推导出房价与地价的结构模型,以北京市为例,运用2003-2013年居住用地价格和2014年在售楼盘价格,与北京市GIS电子地图相匹配,提取商服中心可达性、公共交通可达性、道路可达性、商服繁华度、设施便利性等解释变量,构建地价与房价结构方程模型,分析二者之间的结构关系。     

塔克拉玛干沙漠腹地地表辐射与能量平衡

利用2013年塔克拉玛干沙漠塔中流动沙面地表辐射、土壤热通量、土壤温湿度和湍流通量观测资料,分析了沙漠腹地地表辐射和能量收支特征及闭合状况。结果表明：除潜热通量外,其余地表辐射各分量和能量平衡分量的月平均日变化结果整体均表现为标准的单峰型日循环形态,其中R_s↓与R_s↑变化同步,R_l↑、R_l↓滞后R_s↓0.5 ~ 1 h。各分量均表现出夏季高、春秋季次之、冬季低的季节波动性。干旱和极低的植被覆盖造成沙漠腹地全年潜热通量始终较为微弱,约占净辐射的2.8%,感热通量成为能量的主要消耗形式,约占净辐射的49%。偶尔的降水会刺激潜热通量突然增加。地表反照率相对较高且稳定,日变化呈早晚大、正午小的“U”型趋势,并具有明显的冬季高、夏季低的季节波动性,年均值0.28,月均值0.25~0.32。能量残差各月的日变化也均呈单峰曲线,日出后和日落前能量闭合程度最佳,并出现过闭合现象,全年夏季小,春秋季次之,冬季较大,月平均日峰值5.1~99.9 W·m^-2。土壤表层热储存是影响该地区能量平衡的重要因子之一,考虑表层土壤热存储后,地表能量闭合率达75.3%,能量闭合率夏季 > 春季 > 秋季 > 冬季,白天相比夜间有大幅提升。     

三维复电阻率电磁场正演研究

将Cole—Cole模型表示的地下异常体各剖分单元复电阻率,引入到体积分方程求解均匀大地三维电磁场计算中,实现电偶源地面激发、接收复电阻率三维电磁场正演模拟。     

Soil erosion and its response to the changes of precipitation and vegetation cover on the Loess Plateau

Soil erosion is a major threat to our terrestrial ecosystems and an important global environmental problem. The Loess Plateau in China is one of the regions that suffered more severe soil erosion and undergoing climate warming and drying in the past decades. The vegetation restoration named Grain-to-Green Program has now been operating for more than 10 years. It is necessary to assess the variation of soil erosion and the response of precipita- tion and vegetation restoration to soil erosion on the Loess Plateau. In the study, the Revised Universal Soil Loss Equation （RUSLE） was applied to evaluate annual soil loss caused by water erosion. The results showed as follows. The soil erosion on the Loess Plateau between 2000 and 2010 averaged for 15.2 t hm-2 a 1 and was characterized as light for the value less than 25 t hm-2 a-1. The severe soil erosion higher than 25 t hm-2 a-~ was mainly distributed in the gully and hilly regions in the central, southwestern, and some scattered areas of earth-rocky mountainous areas on the Loess Plateau. The soil erosion on the Loess Plateau showed a deceasing trend in recent decade and reduced more at rates more than 1 t hm 2 a 1 in the areas suffering severe soil loss. Benefited from the improved vegetation cover and ecological construction, the soil erosion on the Loess Plateau was significantly declined, es- pecially in the east of Yulin, most parts of Yah＇an prefectures in Shaanxi Province, and the west of Luliang and Linfen prefectures in Shanxi Province in the hilly and gully regions. The variation of vegetation cover responding to soil erosion in these areas showed the relatively higher contribution than the precipitation. However, most areas in Qingyang and Dingxi pre- fectures in Gansu Province and Guyuan in Ningxia Hui Autonomous Region were predomi- nantly related to precipitation.     

缓变地形下β效应的Rossby代数孤立波

正压流体中,从准地转位涡方程出发采用摄动方法和时空伸长变换推导了在缓变地形下β效应的Rossby代数孤立波方程,得到Rossby波振幅满足带有缓变地形非齐次Benjamin-Davis-Ono(BDO)方程的结论.通过分析孤立Rossby波振幅的演变,指出了β效应、地形效应是诱导Rossby孤立波产生的重要因素,说明了在缓变地形强迫效应和非线性作用相平衡的假定下,Rossby孤立波振幅的演变满足非齐次BDO方程,给出在切变基本气流下缓变地形和正压流体中Rossby波的相互作用.     

基于碳酸盐岩孔隙结构预测孔隙度方法研究

许多研究都已经证实在碳酸盐岩储层中孔隙结构对声波速度影响很大,因此在孔隙度反演时必须考虑孔隙结构的影响.本文通过对Gassmann方程的合理简化并引入Eshelby-Walsh干燥岩石椭球包体近似公式,推导出包含岩石孔隙结构参数、饱和岩石压缩系数与岩石基质压缩系数三个参数的计算孔隙度的新公式,该式岩石基质压缩系数是通过Russell流体因子和Gassmann-Boit-Geer tsma方程计算式建立目标函数反演获得的,然后通过测井交汇图技术把岩石孔隙结构参数与岩石基质压缩系数优选转换成纵横波速度与密度关系式,进而导出具体地区考虑碳酸盐岩孔隙结构孔隙度具体计算公式,测井与地震资料应用表明,基于碳酸盐岩孔隙结构的孔隙度预测方法的精度高于常规方法.