陇东黄土高原下垫面不均匀性指标的建立及其对大孔径闪烁仪(LAS)观测感热通量的影响

目前定量研究下垫面不均匀性对大孔径闪烁仪(LAS)、涡动相关仪(ECS)观测感热通量差异的影响还比较少.本文利用黄土高原庆阳观测站2012年6、7月典型晴天两主风向范围E-SE和SW-W的陆面过程数据,在建立了下垫面不均匀性定量化指标的基础上,分析了下垫面不均匀性与地表温度变率的关系及其对LAS、ECS观测感热通量差异的影响.结果表明:庆阳站下垫面不均匀性大小η和地表温度变率г具有很好的相关性,相关系数达到0.566以上,证明了新建立指标的合理性.下垫面不均匀性大小η和LAS、ECS观测的感热通量差异HLAS-HECS具有很好的一致性,相关系数达到0.634.下垫面不均匀性对LAS和ECS的感热通量差异影响显著,下垫面越不均匀LAS和ECS测量的感热通量差异越大.对两主风向分别进行分析,在E-SE风向范围η和HLAS-HECS的相关系数为0.430,HLAS和HEC拟合的线性趋势系数为1.279,在SW-W风向范围η和HLAS-HECS的相关系数为0.680,HLAS和HEC拟合的线性趋势系数为1.297.下垫面不均匀性的影响程度越大,LAS和ECS观测的感热通量差异越大.     

基于新的差分结构的时-空域高阶有限差分波动方程数值模拟方法

传统的高阶有限差分波动方程数值模拟方法采用高阶差分算子近似空间偏导数,能有效抑制空间频散.然而,传统的有限差分法仅采用二阶差分算子近似时间偏导数,这使得地震波场沿时间外推的精度较低.当采用较大的时间采样间隔,传统的有限差分法模拟波场会出现明显的时间频散,甚至不稳定.本文基于新的差分结构和中心网格剖分,发展了一种空间任意偶数阶精度、时间四阶和六阶精度的时空域有限差分方法.基于对离散后的频散关系进行泰勒展开,本文推导了时空域高阶有限差分算子的差分系数.相速度分析表明时间四阶、六阶精度的差分方法能显著地减小传统时间二阶精度差分方法的时间频散.在相同的精度下与传统差分法比较,本文发展的时间四阶、六阶有限差分方法的计算效率比传统方法高.均匀和非匀均介质中的波场数值模拟实验进一步证实本文研究的时空高阶有限差分方法的优越性.     

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

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

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

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

上地壳纵横波速度结构相关反演成像方法

基于纵横波初至走时数据的层析成像方法越来越广泛地被应用于揭示不同构造域壳幔速度结构特征.我们从同一地质体的纵横波速度属性相关这一基本思想出发,提出一种相关反演成像的方法:纵横波速度反演交替进行,在迭代反演过程中每通过一次反演获得相应的纵波速度(或横波速度)结构后,更新相应的纵横波速度比模型以及相应的横波(或纵波)速度反演的初始模型,然后继续开展后续横波(或纵波)速度反演工作.在反演过程中依据纵横波速度的相关性信息和射线路径长度将走时残差以不同权重分配到射线路径经过的单元,依据网格节点周围平均的慢度扰动更新速度模型.正反演过程分别基于有限差分走时计算方法和反投影成像方法.两种典型模型试验表明,该技术应用于上地壳速度结构反演成像过程,可有效提高反演结果的可靠性,在很大程度上避免了常规单独反演纵波和横波速度过程容易带来的畸变和失真.该方法应用于重建青藏高原西部札达—泉水沟深地震测深(DSS)剖面下方的上地壳速度结构,揭示出与青藏高原西缘板块碰撞相关的上地壳速度结构特征.     

基于Love波相速度反演南北地震带地壳上地幔结构

收集了南北地震带区域地震台网中292个地震台站2008年1月至2011年3月期间的地震波形数据,由频时分析方法提取了Love波相速度频散曲线,经过反演得到了研究区内的Love波相速度分布.根据Love波纯路径频散,采用线性反演方法对0.25°×0.25°的网格点进行了一维S波速度结构反演,利用线性插值获取了南北地震带地区的三维S波速度结构.结果显示了松潘—甘孜地体和川滇菱形块体地区的下地壳具有明显的S波低速层分布,该异常分布特征支持解释青藏高原隆升及其地壳物质运移的下地壳流模型.在100至120km深度上,川滇菱形块体西北部呈现较强的S波高速异常,这可能是印度岩石圈板块沿喜马拉雅东构造结下插至该区域所致,该区域下地壳的低速软弱物质与上地幔的高速强硬物质形成了鲜明对比,暗示了地壳和上地幔可能具有不同的构造运动和变形方式,这为该区域的壳幔动力学解耦提供了条件.     

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

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

西北内陆河流域典型生态系统通量数据空间代表性研究

涡度通量数据的空间代表性问题是影响其数据质量不确定性的重要原因。选取了我国西北地区具有代表性的四个通量观测站生长季数据,采用FSAM模型对各站点通量贡献源区范围以及通量贡献最大值点进行了计算,给出了在不稳定与稳定两种大气层结条件下各站点源区的主要分布范围(P=80%):盈科站源区范围相对最小(90~200 m);阿柔站源区范围约为(120~140 m);关滩站源区范围约为(500~600 m);敦煌站源区范围约为(140~280 m)。同时结合各站点仪器架设高度以及不同的下垫面生态环境等因素,对通量源区的变化进行了初步分析,仪器的架设高度直接影响涡度源区范围;不同的下垫面生态环境形成不同的近地面湍流状态进而对通量源区的分布产生间接影响。结果表明各站点观测的通量信息基本均来源于所感兴趣的研究区域。     

基于区域航空运输可达性的机场空间布局

构建表现机场腹地与其航空运输网络间相互关系的结构方程模型,明确“二三产业总产值/总人口数/社会消费品零售总额”与“通航机场数/航班密度/通航里程”2个变量组之间的数量关系。基于腹地陆路交通网络和机场航空运输网络,列出机场陆域可达性和航线可达性的数学表达式,并基于乘客的主观感受组合2类可达性指标得到区域航空运输可达性指标,将该指标作为比选机场空间布局方案的标准。以烟台/威海地区为例对不同机场布局方案的场景进行分析,结果表明：在目前的需求形态下烟台/威海地区应设置1个机场,以使旅客享有较好的航空运输服务,使整个区域的航空运输可达性最优。     

相邻历元误差相关的抗差卡尔曼滤波算法分析

根据卡尔曼滤波理论以及抗差理论,本文推导出了相邻历元误差相关的抗差卡尔曼滤波模型,其对观测值中含有粗差有良好的抗差性。通过对含有粗差的变形监测数据分析,与相邻历元误差相关的卡尔曼滤波模型进行比较,采用本文构造的抗差卡尔曼滤波模型处理数据,无论是否有粗差存在观测值里,变形计算结果与实际结果大体一致,粗差对计算结果的影响不敏感。在对变形监测数据分析时,可得出卡尔曼滤波方法估计的状态向量,没有寄存大量的以往观测数据,而是使用最近的观测数据,经过不断的预测和改正,把新的状态展示在系统中。