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

等效水力隙宽和水力梯度是影响岩体裂隙网络渗流特性的重要因素。制作裂隙网络试验模型,建立高精度渗流试验系统;求解纳维-斯托克斯方程,模拟流体在裂隙网络内的流动状态,研究等效水力隙宽和水力梯度对非线性渗流特性的影响。结果表明,当水力梯度较小时,等效渗透系数保持恒定的常数,流体流动属于达西流动区域,流量与压力具有线性关系,可采用立方定律计算流体流动;当水力梯度较大时,等效渗透系数随着水力梯度的增加而急剧减少,流体流动进入强惯性效应流动区域,流量与压力具有强烈的非线性关系,可采用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电子地图相匹配,提取商服中心可达性、公共交通可达性、道路可达性、商服繁华度、设施便利性等解释变量,构建地价与房价结构方程模型,分析二者之间的结构关系。     

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

将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方程计算式建立目标函数反演获得的,然后通过测井交汇图技术把岩石孔隙结构参数与岩石基质压缩系数优选转换成纵横波速度与密度关系式,进而导出具体地区考虑碳酸盐岩孔隙结构孔隙度具体计算公式,测井与地震资料应用表明,基于碳酸盐岩孔隙结构的孔隙度预测方法的精度高于常规方法.     

TTI介质qP波逆时偏移中伪横波噪声压制方法

在对地下复杂构造介质,特别是盐丘侧翼及岩下区域进行成像时,相对于传统的各向同性逆时偏移和VTI逆时偏移,具有倾斜对称轴的TTI逆时偏移成像效果最优.不仅反射同相轴更加的连续,而且能量得到了更好的聚焦.传统的各向异性介质全弹性波RTM的计算量大且计算效率低.由于目前仍以纵波勘探为主,因此TTI逆时偏移qP波波动方程的选取显得尤为重要.为了提高计算效率,采用将沿着对称轴方向的横波速度设为零的方法,简化得到qP波波动方程.然而,这样会引入一种严重影响成像效果的低速度、低振幅的qSV波人为干扰.本文建立了qP波方程的完全匹配层控制方程,而后借助于辅助波场采用一种高效的压制伪横波噪声传播的方法,通过模型测试验证了该方法的有效性.     

地震波动方程的局部间断有限元方法数值模拟

近年来,随着地震波数值模拟对计算精度和效率的要求越来越高,间断有限元方法开始受到越来越多的关注.本文中,针对具有吸收边界条件的二维地震声波波动方程,作者提出了一种基于局部间断有限元方法的数值模拟算法.该算法在空间上使用局部间断有限元方法进行离散,在时间上采用了显式蛙跳格式.在这种时空离散的组合方式下,每个时间步上,此算法在空间剖分的每个单元上的求解计算是相互独立的,因而具有极高的并行性.通过数值算例,我们将该算法与连续有限元方法进行了比较.结果表明,本算法不仅具有对起伏构造的良好适应性,而且在计算效率和计算精度等方面,都具有优越性.     

基于逆散射级数法的鬼波压制方法

鬼波问题是影响海上地震资料分辨率和保真度提高的最重要因素之一.详细论述了逆散射级数理论和逆散射级数法鬼波压制原理,说明了逆散射级数方法进行鬼波压制理论的完善性和对鬼波描述的精确性.实现了基于逆散射级数理论的鬼波压制方法,方法以波动方程和Lippmann-Schwinger方程为基础,在频率-波数-波数域内构造与鬼波相关的压制算子,在不需要对地下介质作任何假设条件下实现地震数据驱动鬼波压制,并通过改善消除鬼波的压制算子,提高算法的稳定性.资料处理试验与处理结果分析表明,基于逆散射级数鬼波压制方法能在实现鬼波压制的同时较好保留有效反射波的信息,从而补偿地震资料低频损失和提高地震数据的保真度.数据处理试验还表明,研究方法能对低信噪比的地震资料进行有效的鬼波压制处理.建立了基于逆散射级数鬼波压制处理流程.