求解弹性波方程的辛RKN格式

将弹性波方程变换至Hamilton体系,构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nystrm(RKN)格式,运用根数理论得到此格式的阶条件方程组.通过给定系数的限定条件,得到方程的对称解.为了使时间离散误差达到极小,提出数值频率与真实频率比较,通过Taylor展开,得到关于辛系数的限定方程,求解方程组得到最小频散辛RKN格式.对比分析时间演进方程的稳定性,得到使库朗数达到极大值的限定方程,求解方程组得到最稳定辛RKN格式.发现此两种格式为同一格式.新得到的辛RKN格式不依赖于空间离散方法,为了对比的需要,选取有限差分法进行空间离散.在频散、稳定性分析中,与常见辛格式对比,从理论上分析了本文提出的格式在数值频散压制、稳定性提升等方面的优势,数值实验进一步证实了理论分析的正确性.     

VTI介质纯P波混合法正演模拟及稳定性分析

各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性.     

基于广州实证的后改革时代中国城市住房权问题

社会公正是当前地理学的一个热点也是传统研究课题。住房问题往往是社会公正研究和实践的重点案例。中国住房制度改革之后,中国住房拥有率及人均住房面积均大幅度提高。然而在住房分配的过程中,不平等和不公正现象时有发生。如何区别平均和公正,如何看待中国住房改革中的社会公平与不公共存的现象,是了解中国新城市现象的一个重要途径。本文从近期西方政治学和地理学关于住房权（housing rights）、社会公正以及城市权利等文献的讨论出发,通过在广州的实证研究,运用结构方程模型和Pratt指数对涉及住房权主要方面的住房分异的发生途径及其公正性进行综合分析。实证研究说明在住房分配领域,存在着制度因素延续下来所导致的分异。一些制度因素已经转化为市场因素,直接或间接影响个人居住条件和住房所有权的获得,并通过后者进一步影响住房分异。与此同时,纯粹市场因素的影响则相对薄弱。上述结论有助于了解当前中国城市住房权的核心问题,认识社会公正的发展前景,以及探寻减少社会不公正的制度利器。     

基于Gassmann方程的流体替换方法在珊溪水库地震研究中的应用

引入基于Gassmann方程的流体替换方法,在分析地震波P波速度、波速比与岩石孔隙度和饱和度关系的基础上,应用于珊溪水库地震波速比和P波速度变化特征研究,得到:(1)珊溪水库震中区岩石始终处于接近水饱和的饱水状态,波速比和P波速度"下降-回升"的变化实质上反映了震中区岩石"孔隙度增大(饱和度减小)-饱和度增大"的变化,每一丛地震的波速比由极小值逐渐增大为极大值是由于岩石从不饱和状态变化到饱和状态;(2)根据每一丛地震波速比的变化,计算得到珊溪水库流体扩散率αs=1.06×104 cm2 s-1,该数值与美国南卡罗莱纳水库、巴西Acu水库、广东新丰江水库的流体扩散率基本一致;(3)震源区岩石孔隙度上限值为8.7%~2.0%,该数值与华东勘测设计研究院通过室内岩石物理力学性质试验测定的珊溪水库坝址区新鲜流纹斑岩的孔隙度平均值一致。     

页岩中气体的超临界等温吸附研究

页岩气等温吸附实验多为临界温度以上的吸附实验,其得到的吸附量为过剩吸附量。为了研究页岩气超临界等温吸附机理,运用重力法,在临界温度以上,分别进行了甲烷和二氧化碳在页岩中的高压等温吸附实验。在分析经典型吸附和超临界吸附区别的基础上,通过修改的超临界等温吸附模型(Langmuir方程和微孔充填(Dubinin Radushkevich,D-R))对实验数据进行了拟合。结果表明:简单的Langmuir方程可近似拟合甲烷吸附实验数据,但精度不高,且无法拟合二氧化碳的吸附数据;将吸附相密度作为可优化参数,修改的微孔充填模型和Langmuir模型能很好地拟合甲烷和二氧化碳的吸附数据,其中修改的微孔充填模型拟合效果最好,且回归得到的超临界甲烷吸附相密度同文献报道的一致,表明吸附气可能以微孔充填的形式存在。      

基于单程波偏移算子的地表相关多次波成像

在常规地震资料处理中,多次反射波被视为噪声并从地震数据中去除,以免在之后的地震资料解释中造成误解.而事实上,多次波也是地震信号,是照明波场的一部分,能够对地下构造成像的精度做出贡献.本文分析了多次波在传统单程波叠前深度偏移中产生构造假象的机制和表现,为实现基于单程波偏移算子的多次波成像,修改了单程波叠前深度偏移的边界条件,即将输入的震源波场用包含多次波的记录来替代,输入的记录波场用预测出的表层相关多次波来替代,实现了基于单程波偏移算子的地表相关多次波成像,并从理论上给出了其成像依据.通过基于二范式最小能量差原则求取的匹配因子,将多次波成像结果与一次波成像结果进行匹配叠加,应用多次波成像来弥补一次波成像的不足.简单模型验证了基于单程波偏移算子的多次波成像方法的有效性,最后对Sigsbee2B模型进行了一次波与多次波联合成像试算,盐边界高陡构造成像质量得到了明显改善.     

First-principles derivation of reactive transport modeling parameters for particle tracking and PDE approaches

Both Eulerian and Lagrangian reactive transport simulations in natural media require selection of a parameter that controls the “promiscuity” of the reacting particles. In Eulerian models, measurement of this parameter may be difficult because its value will generally differ between natural (diffusion-limited) systems and batch experiments, even though both are modeled by reaction terms of the same form. And in Lagrangian models, there previously has been no a priori way to compute this parameter. In both cases, then, selection is typically done by calibration, or ad hoc. This paper addresses the parameter selection problem for Fickian transport by deriving, from first principles and D (the diffusion constant) the reaction-rate-controlling parameters for particle tracking (PT) codes and for the diffusion–reaction equation (DRE). Using continuous time random walk analysis, exact reaction probabilities are derived for pairs of potentially reactive particles based on D and their probability of reaction provided that they collocate. Simultaneously, a second PT scheme directly employing collocation probabilities is derived. One-to-one correspondence between each of D, the reaction radius specified for a PT scheme, and the DRE decay constant are then developed. These results serve to ground reactive transport simulations in their underlying thermodynamics, and are confirmed by simulations.     

固体介质中孤立波的传播及演化特征

由于地球介质中广泛存在断裂、微裂缝等地质现象,实际地震资料中会出现形似孤立波的非线性地震现象。因此,对固体介质中孤立波的研究有利于解释这些非线性地震现象的形成机制。本文基于KdV方程,以雷克子波作为初始条件,采用伸缩子机理构建体力模型,利用有限差分的方法模拟了孤立波的演化过程。理论结果表明,非线性地震纵波可以从雷克子波逐渐演化成孤立波,而且地震波的初始振幅和频散系数对模拟结果也有重要影响。通过与实际资料的对比也能说明这种演化的可能性。同时根据方程系数矩阵中元素带状分布的特征,采用稀疏矩阵的存取方法,可以减小计算内存,提高计算效率。     

基于声学全波形反演的油气藏地震成像方法

岩性油气藏在我国天然气勘探开发中占有非常重要的位置,其分布区域的成像是合理布设井位,提高钻井成功率的关键之一.本文首先基于地下介质的声学近似和波场回传理论,利用频率域单程声波方程延拓计算地震波场,进行全波形反演,获得地层密度和体积模量的定量成像,并依据油气藏物性特征和流体饱和多孔介质岩石物理模型,简要讨论了孔隙度和饱和度与密度及体积模量的关系,明确了地震油气藏成像新概念.在此基础上,定义了基于流体体积模量和孔隙度的成像函数,进行油气藏成像.理论模型计算表明该方法是可行的.通过对西部地区某气田二维地震数据处理,实现了致密砂岩气藏成像,钻井结果证实了气藏区域成像位置的准确性和方法的有效性.