缔合勒让德函数的标准向前按列递推公式

针对推导扇谐项（次等于阶的项）和准扇谐项（次等于阶减1的项）存在不严谨性的问题,给出基于fnALFs的2个导数定积分按列递推公式,结果与已有的研究结论完全一致,但推导过程更简明,且该定积分公式可从第2阶开始递推。另外,本文还给出另外2个fnALFs定积分计算公式。     

基于PML边界下的弹性及黏弹性TTI介质波场模拟

基于一阶速度-应力波动方程,采用高阶交错网格有限差分数值模拟方法,对弹性及黏弹性TTI介质进行正演数值模拟。模拟时采用完全匹配层吸收边界条件(PML)消除边界反射。同时设计了层状介质模型、断层模型,通过模型的正演计算,得到了不同时刻的地震波波场快照及合成地震记录,分析其波场运动学及动力学特征。模拟结果表明,交错网格有限差分法可以很好地完成对复杂介质的波场模拟,具有较高的精度和可靠性。     

GEM10B地球模型全球自由空气重力异常场及勒让德函数计算的研究

本文给出了GEM 10 B地球重力模型2-36阶、2-12阶以及13-36阶等三幅全球自由空气重力异常图.在2-12阶和2-36阶重力图中,全球重力异常呈现出双层环状结构特征.内环是环绕太平洋地区的正异常带,外环是紧密包围着环太平洋正异常带的巨大负异常带.在13-36阶重力图中,全球大洋海岭地区几乎全部呈低幅值（0-10mGal）正异常. 本文分析了利用递推公式计算勒让德函数的稳定性问题,给出一组适合计算机运用的完全正规化缔合勒让德函数计算公式.采用该组公式计算（双精度运算）的完全正规化缔合勒让德函数值在θ≥1°,L≤300阶时,其ε_L误差小于10^-13. 为了解决对高阶勒让德函数值计算精度的估计问题,本文提出以完全正规化缔合勒让德函数的σ_l=1性质为基础的σ_L函数精度判定法.该方法简便易行.     

高阶交错网格非稳态相移叠后深度偏移

为更好地研究复杂构造和速度分布条件下地震资料的精确偏移.在波场模拟方面,本文采用了高阶交错网格有限差分[1]推导了二维各向同性介质声波方程的数值模拟公式.并给出了声波方程完全匹配层吸收边界[2-5](PML).在偏移方面,本文介绍了一种能适应介质速度横向变化的非稳态相移算子(NSPS)[6]及其叠后深度偏移方法.并对简单的四层模型和Marmousi模型进行试算.偏移结果表明,该方法可以有效地处理复杂条件下的各向同性介质声波偏移问题.     

用于地震波场模拟的PML边界衰减因子研究(英文)

在地震波场数值模拟中,为了消除由人为边界产生的边界反射,需要引进边界吸收条件。本文从声波方程完全匹配层吸收边界的经典方法出发,基于高斯函数任意阶光滑的特点,提出了一种高斯型衰减因子,分析比较该因子与一般衰减因子的性质,并基于均匀与层状速度模型分别进行了数值模拟计算。数值结果显示,当选择相同的PML边界吸收层层数时,高斯型衰减因子的吸收效果明显优于一般的衰减因子,边界反射更少,信噪比更高;对比最近提出的正弦型衰减因子,在信噪比接近的情况下,高斯型衰减因子所需的PML吸收层层数更少。     

高阶旋转交错网格有限差分方法模拟TTI介质中横波分裂

笔者给出了一种能够模拟弹性波在任意各向异性介质中传播的二维三分量高阶有限差分算法.相对于常规交错网格有限差分方法,旋转交错网格有限差分方法在介质具有强差异性时能更精确地模拟地震波的传播,避免常规交错网格中因对弹性系数进行插值而带来的误差.采用高阶旋转交错网格有限差分方法模拟并分析了零偏移距横波分裂现象随裂缝介质方位角和倾角变化的响应特征.结果表明:结合完全匹配层(PML)吸收边界条件的高阶旋转交错网格有限差分方法能获得高精度的地震波场模拟数据,并且在边界具有良好的吸收效果;横波分裂现象主要受裂缝走向与波的极化方向之间的夹角影响,受裂缝倾角影响较小,且快慢横波的能量也跟裂缝走向与波极化方向间的夹角有关.具有倾斜对称轴的横向各向同性(TTI)介质倾角的变化可能会导致记录中波到达时的变化,影响快慢横波的时差.利用横波分裂的能量分布和方位各向异性特征,可以帮助检测裂缝的方位角和倾角.横波在多层TTI介质中传播时会发生多次分裂的现象.     

准P波方程紧致交错网格井间地震波场模拟及边界条件

研究井间地震波场的形成过程以及波场的传播机理、规律,对于指导实际井间地震勘探有着重要的意义.基于具有垂直对称轴的横向各向同性(VTI)介质中的一阶准P波方程,应用具有无条件稳定性质的紧致交错网格隐式差分格式求解该方程.重点研究了紧致交错网格求解该方程的完全匹配层(PML)吸收边界条件,在此基础上实现了VTI介质中一阶准P波方程的井间地震波场的正演模拟.数值算例表明:紧致交错网格能精准模拟VTI介质中准P波的传播过程,得到高精度的正演结果.一阶准P波方程能以足够的精度描述VTI介质中准P波特征.完全匹配层吸收边界能有效地解决人工边界问题,是一种高效的边界吸收算法.     

完全规格化缔合勒让德函数常用递推算法的适用性

完全规格化缔合勒让德函数递推算法的适用性是衡量算法优劣的重要标志。从第一相对数值精度、第二相对数值精度和计算速度等方面对4种常用的递推算法--标准向前列递推算法、标准向前行递推算法、跨阶次递推算法和Belikov递推算法的适用性进行分析。结果表明,标准向前行递推算法适用范围最小;对于cosθ∈[－1,1],在1 900阶内,标准向前列递推算法、跨阶次递推算法和Belikov递推算法均适用,且第1种算法速度最快;在3 000阶内,跨阶次递推算法和Belikov递推算法适用,且后者更优。     

测量上海产业网络的点入度和点出度——超越后工业化社会的迷思

利用上海市2007年144个部门投入产出表,计算基于投入产出完全消耗系数的点入度和点出度,将产业划分为Ⅰ型（高点入度、高点出度）,Ⅱ型（低点入度、高点出度）,Ⅲ型（低点入度、低点出度）,Ⅳ型（高点入度、低点出度）。研究发现：① Ⅰ型、Ⅱ型、Ⅲ型和Ⅳ型制造业点入度均值分别比各类型服务业高0.52、0.45、0.04、0.27,表明制造业相对服务业具有更强的后向拉动效应。② Ⅱ型和Ⅳ型服务业点出度均值分别高于Ⅱ型和Ⅳ型制造业0.29、0.16,而Ⅰ型和Ⅲ型服务业点出度均值分别小于Ⅰ型和Ⅲ型制造业0.32、0.17,表明服务业的前向拉动效应存在着分异。③ Ⅰ型、Ⅳ型制造业及服务业的点出度对GDP均具有正向的相关性,提高对这些产业的使用有利于提升GDP。Ⅱ型服务业的点入度对GDP具有负向的相关性,点入度较低的房地产开发经营业、银行业、证券业GDP分别为该类服务业前3名。Ⅱ型服务业点出度、Ⅲ型制造业和服务业点入度和点出度对GDP的回归系数均未通过显著性检验。尽管制造业GDP占比下降,制造业仍在上海产业网络中发挥着重要的作用。上海应超越后工业化的迷思,降低对房地产业的依赖,加强制造业与服务业的融合发展,特别是加强金融对制造业等实体经济的支持,提升产业网络关联多样性和非关联多样性,促进航空、航天、机器人等新兴产业的发展。     

基于城镇化视角的绿洲城市用水变化驱动效应分析

城镇化是驱动区域用水变化的双刃剑,既可以通过刺激生产来促进用水总量增长,也可以通过优化用水结构和提高用水效率来抑制用水总量增长。以黑河流域张掖市为例,通过构建城镇化驱动用水变化的完全分解模型,定量测度了2000-2011年张掖市人口和经济城镇化过程对用水总量和用水效率变化的驱动效应。结果表明：张掖市人口和经济城镇化过程对用水总量均具有明显的减量效应,对用水效率均具有明显的增量效应,而且与人口和经济城镇化速度高度相关。因此建议干旱区绿洲城市在水资源约束下,应将城镇化作为优化人口结构、经济结构和用水结构,提高经济效率以及用水效率,并最终实现水-生态-经济良性循环的抓手。该结论对于科学认识干旱区绿洲城市城镇化与水资源利用的关系具有一定的参考价值。