基于Duffing振子混沌系统的地震速度分析方法

强随机噪声干扰是导致地震勘探资料低信噪比的主要原因,如何在强随机噪声干扰下获取有效的信息是值得关注的问题.Duffing振子混沌系统是一个非线性的动力学系统,其对强随机噪声具有免疫能力,而对特定的周期性信号具有敏感性.本文提出一种基于Duffing振子混沌系统的速度分析方法.对CMP道集按照时距曲线关系进行移动窗口截取,将所截取的信号构建为待测信号加入Duffing振子混沌系统,通过相图网格分割方法(GPM)判断系统状态的改变,从而在强随机噪声背景下获得高分辨率的速度谱.理论模型和实际资料的处理结果表明,与传统的水平叠加速度分析方法相比,本方法能够在强随机噪声背景下获得更准确的速度分析结果.     

InSAR约束下的2008年汶川地震同震和震后形变分析

2008年5月12日,青藏高原东缘的龙门山断裂带上发生了M_w7.9级汶川地震.本文通过分析覆盖汶川地震震中区域的ALOS/PALSAR像对的方位向偏移量来选择无明显电离层扰动影响的像对进行干涉处理,获取了高精度、连续的InSAR地表形变场.在此基础上,结合高精度GPS同震形变数据,采用同震、黏弹性松弛震后形变联合反演模型同时确定了汶川地震的同震滑动分布和龙门山地区的流变结构参数.研究结果表明,汶川地震是一个断层破裂非常复杂的地震事件,其中,北川段、岳家山段、虹口段和汉旺段的滑动以逆冲为主,而青川段以右旋走滑为主.滑动主要发生在10 km深度以上的区域,最大滑动量位于虹口段的东北端,达10.7 m.地震释放的总能量为9.28×10²⁰ N·m（M_w7.91）,与地震学的结果一致.联合反演模型确定的龙门山地区中下地壳的黏性系数为2×10¹⁸ Pa·s,为青藏高原东部地区的黏性系数提供了一个可靠的下限值.如果有更长时间的震后形变观测时间序列,将为该区域提供更为可靠的流变结构.     

一阶弹性波方程数值模拟中的混合吸收边界条件

Liu和Sen（2010和2012）在地震波场数值模拟中提出一种混合吸收边界条件. 该方法具有计算量小、容易实现及吸收效果好等优点. 但现有的混合吸收边界条件是针对二阶位移-应力方程设计的,存在稳定性问题. 本文首先推导了两种速度-应力单程波方程：二阶Higdon单程波方程和一阶Higdon单程波方程. 进而提出基于一阶弹性波方程的混合吸收边界条件方法. 在内部区域和边界之间引入一个过渡区域,通过单程波与双程波方程平滑过渡来消除人工边界反射. 为了改善混合吸收边界条件的吸收效果和稳定性,我们采用了能同时吸收纵、横波反射的一阶单程波方程和与变量位置有关的加权系数. 为了验证混合吸收边界条件的有效性,将其与常规分裂完全匹配层（PML）方法进行了比较. 数值模拟结果表明,与PML边界条件相比,混合吸收边界条件在耗用更小计算时间和存储量的前提下,可以获得更好的吸收效果. 另外,本文提出的两种混合吸收边界条件中,混合一阶Higdon吸收边界条件具有更好的稳定性.     

基于优化时空域高阶有限差分方法的粘滞声波叠前逆时偏移(英文)

叠前逆时偏移是对地下介质进行精确成像的方法之一。由于实际地下介质具有粘滞性,研究粘滞声波叠前逆时偏移具有一定的现实意义。逆时偏移的步骤之一是求解波动方程,对地震波场进行正向和反向外推,因此,精确、高效地求解波动方程对逆时偏移的成像效果和计算效率具有重要影响。本文中,我们利用基于优化时空域频散关系的高阶有限差分方法求解粘滞声波方程。频散分析和数值模拟的结果证明了优化时空域有限差分方法具有较高的精度,可以很好地压制数值频散。利用混合吸收边界条件处理边界反射,然后利用震源归一化互相关成像条件进行成像,并利用拉普拉斯滤波方法去除低频噪音。数值模型的测试结果显示,在考虑地下介质的粘滞性时,粘滞声波方程逆时偏移比声波方程逆时偏移具有更高的成像分辨率。另外,在进行波场外推的时候,采用自适应变长度的有限差分算子计算空间导数,在不影响求解精度的情况下,有效地提高了计算效率。     

多块体形状堆石体碾压颗粒破碎数值模拟

通过6种典型堆石块体的形状近似,分别采用两种接触本构模型建立了多块体形状堆石体离散元数值模型,研究堆石体在碾压荷载作用下的颗粒破碎过程,建立颗粒破碎的量化计算方法,分析碾压前、后堆石级配曲线的变化,讨论接触本构模型和颗粒形状对块体破碎的影响。模拟结果显示,碾压荷载下堆石颗粒以张拉破碎为主,随着碾压遍数的增加,局部开始出现剪切破碎;提出的颗粒破碎量化计算方法,在大粒径范围对粒径变化幅度预测偏大,但级配曲线整体趋势与实测结果比较吻合;相比较于接触连接模型,平行连接模型与现场碾压试验结果更接近。6种块体形状的数值结果显示,随着形状系数的增加,在其他条件不变的情况下颗粒破碎率逐渐降低,其中类长方形块体在碾压荷载作用下颗粒破碎最明显。相比较于纯圆颗粒或者单一非圆颗粒,采用的6种颗粒形状建立的堆石体振动碾压离散元模型,更接近现场实际情况。     

林下植被遥感反演研究进展

林下植被在森林生态系统碳、水和营养元素的累积和循环方面有着重要作用与科学意义。多角度、高光谱和激光雷达遥感系统凭借对森林分层结构的敏感性,成为量化林下植被的重要手段。本文综述了森林林下植被的遥感反演研究进展：首先讨论了林下植被的定义,其次对当前林下植被的遥感反演现状作了深入分析,总结了多角度、高光谱和激光雷达遥感观测的森林背景的反射率、林下植被的叶面积指数、高度和覆盖度的遥感反演原理和方法。基于不同卫星观测角度下森林冠层和森林背景对总反射的贡献差异,林下反射率可通过多个角度的观测数据进行反演。此外,借助激光雷达穿透冠层直接观测林下植被的优势,总结了激光点云数据和回波波形信息反演林下植被的覆盖度和高度的方法,以及今后使用遥感技术反演的难点和获取林下植被信息的主要发展方向。     

国外空间规划研究与实践的新动向及对我国的启示

空间规划是可持续发展必不可少的公共管理工具。20世纪90年代以来,人们认识到空间规划应该具有协调和整合空间发展的功能,而不仅仅是部门性的土地利用管理体系;应该关注地域的可持续发展,而不是一味地追求增长;应该关注不同利益人之间共识的形成,而不是简单地分析一定地域空间开发的适宜性和先后顺序。20世纪90年代以来国外空间规划实践的一些新动向： a) 重新分配空间规划的权利和职能; b)注重全纳性的参与; c) 重视空间规划的协调与整合; d)强化空间规划的环境维度; e)增强空间规划的灵活性。基于此对我国的空间规划实践得到几点启示： a)空间发展的战略框架和行动计划相结合; b)形成适合于我国的规划参与机制; c)在增强规划约束性的基础上增加灵活性; d)强化区域层次的空间规划; e)突出空间规划的公共管理职能。     

加拿大北方森林火烧迹地遥感分析

森林火灾是加拿大北方森林地区最主要的扰动因素,对北方生态系统起着主导作用。基于MODIS数据,采用全球扰动指数算法(MGDI),对加拿大萨斯喀彻温省和亚伯达省2004-2011年的森林火烧迹地进行检测和分析。通过与MODIS标准火烧迹地产品以及加拿大林业局数据进行比较,扰动指数算法检测的火烧迹地面积比MODIS标准产品更接近于林业局的统计数据。分析表明,在2004-2011年间,由于火灾原因,整个研究区森林面积平均每年减少76.36万hm2,占该区域森林总面积的3.36%。萨斯喀彻温省平均每年燃烧的森林面积为46.83万hm2,亚伯达省为29.53万hm2。其中,2006、2008、2010和2011年是火灾的高峰年份。火烧迹地主要集中在生态交错带的北方保护区、针叶林保护区、针叶林平原区,以及北方平原东北部的伍德布法罗国家森林保护区。     

叠前三参数非高斯反演方法研究

针对地球物理反演中广泛采用的"噪声高斯分布假设",本文研究了叠前地震资料中噪声的非高斯分布特征,提出了针对非高斯噪声的地震叠前非高斯反演概念和思想,构造了能同时压制高斯和非高斯噪声的混合范数作为反演目标函数,采用改进的Powell算法进行求解,有效地抑制了叠前地震资料中的高斯和非高斯混合噪声.模型试算和实际地震数据的反演结果验证了方法的正确性和算法的可靠性.     

黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.