网格剖分及其精度和计算量分析

网格剖分程度直接影响着地震波正演数值模拟的计算精度及其计算量.以均匀介质模型为例,分析不同网格大小对波场模拟精度和计算量的影响,得出精细化网格剖分是实现高精度地震波正演模拟的有效方法,然而其计算量较大.以均匀倾斜介质模型为例,探讨倾斜地层网格剖分问题,数值实例分析不同震源频率对不同网格剖分方案引起的波场传播精度的影响,...     

各向同性介质中射线理论与有限差分地震波场模拟方法比较研究

地震波场数值模拟不仅是研究复杂地区地震资料采集、处理和解释的有效辅助手段,而且是研究地球深部精细构造和地球深部探测的有效工具.射线理论和波动方程理论是地震波场数值模拟的理论基础.射线理论主要刻画地震波在介质中传播的走时场、地震射线等运动学属性;波动方程理论通过求解波动方程来描述地震波在介质中传播的弹性动力学响应(能量衰减、相位特征、偏振属性、以及全波形等).基于波动方程理论的波场数值模拟由于能够引入丰富的波场信息,使得人们对不同介质中地震波的传播过程有了较全面的了解.本文以二维层状均匀介质模型为例,通过射线追踪法和交错网格有限差分法模拟得到的波场快照图、单炮地震记录剖面、合成理论地震图的分析比较,不但对地震波在各向同性层状均匀介质中的传播规律和特点有了深刻的理解和认识,同时又可以相互验证两种不同方法的正确性和有效性.     

求解双相和黏弹性介质波传播方程的间断有限元方法及其波场模拟

间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相.     

线性粘弹介质中地震波场数值模拟

本文将以往两种粘弹介质中地震波模拟方法的优点结合起来，以模型理论和积分本构方程为基础，从理论上分析了模型对地震波场的影响；采用交错网格有限差分法对粘弹介质中的地震波进行数值模拟.数值计算结果表明该方法不仅便于计算，同时也便于从力学的角度来分析地震波的传播.数值计算结果与理论分析一致，说明这种方法可以更为有效地模拟粘弹介质中地震波的传播.     

基于旋转交错网格伪谱法和L-S理论的热弹性介质波场模拟与特征分析

基于L-S热弹性理论,采用旋转交错网格伪谱法,实现了均匀各向同性介质的一阶速度-应力-温度微分方程组的数值求解和波场模拟.其中,用时间分裂法解决方程组的刚性问题,用旋转交错伪谱算子计算空间一阶导数,用中心伪谱算子计算空间二阶导数;对于热导率变化比较大的双层介质模型和参考温度随深度按梯度分布的模型,用Crank-Nicolson显式方法取代旋转交错网格伪谱法进行计算;讨论了热耦合波场性质和传播规律,对比了常规伪谱法、交错网格伪谱法和旋转交错网格伪谱法热耦合波场模拟效果.数值模拟结果表明:基于L-S热弹性理论,用时间分裂法结合旋转交错网格伪谱法对均匀各向同性介质的波场计算,能够得到稳定的、高精度的模拟结果,但是在热导率变化剧烈的情况下不能用大时间步长进行求解,而且在参考温度不均匀分布的情况下算法不稳定.本文在网格剖分方式与数值算法的优化组合应用方面进行了探索,为将这些方法推广到孔隙热弹性、热黏弹性和各向异性研究奠定了基础.     

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

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

波场模拟中的数值频散分析与校正策略

波动方程有限差分法正演模拟,对认识地震波传播规律、进行地震属性研究、地震资料地质解释、储层评价等,均具有重要的理论和实际意义.但有限差分法本身固有存在着数值频散问题,数值频散在正演模拟中是一种严重的干扰,会降低波场模拟的精度与分辨率.针对TI介质波场模拟的交错网格有限差分方法,本文从空间网格离散、时间网格离散和算子近似等三个方面对其产生的数值频散进行了分析,并结合其他学者的研究成果给出了TI介质波场模拟中压制数值频散的方法与策略：在已知介质频散关系时,对差分算子可实施算子校正;通过提高差分方程的阶数来提高波场模拟精度;采用流体力学中守恒式方程的通量校正传输方法来压制波场模拟中的数值频散;在实际正演模拟时,采用交错网格高阶有限差分方程,不仅在空间上采用高阶差分,而且在时间上也要采用高阶差分,否则只在单一方向上（空间或时间）提高方程的阶数对压制数值频散也不会取得理想的效果.     

修正辛格式有限元法的地震波场模拟

三角网格有限元法具有网格剖分的灵活性,能有效模拟地震波在复杂介质中的传播.但传统有限元法用于地震波场模拟时计算效率较低,消耗较大计算资源.本文采用改进的核矩阵存储(IKMS)策略以提高有限元法的计算效率,该方法不用组合总体刚度矩阵,且相比于常规有限元法节省成倍的内存.对于时间离散,将有限元离散后的地震波运动方程变换至Hamilton体系,在显式二阶辛Runge-Kutta-Nystr9m(RKN)格式的基础之上加入额外空间离散算子构造修正辛差分格式,通过Taylor展开式得到具有四阶时间精度时间格式,且辛系数全为正数.本文从理论上分析了时空改进方法相比传统辛-有限元方法在频散压制、稳定性提升等方面的优势.数值算例进一步证实本方法具有内存消耗少、稳定性强和数值频散弱等优点.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

改进的高精度拟谱微分方法及其在三维非均匀介质地震传播研究中的应用

改进了用于模拟地震波场的傅里叶拟谱微分方法,改进后的方法精度是常规拟谱方法的4倍,称为改进的傅里叶拟谱方法.在较高数值精度的一阶应力-速度弹性波动方程的基础上,采用该方法和常规拟谱方法对Marmousi模型进行数值求解,结果表明,该方法的数值频散效应明显比常规拟谱方法弱.将该方法与有限元方法在各向异性介质中进行模拟比较,发现该方法的精度接近有限元方法,数值频散效应比有限元方法明显减小,而且可在较大空间网格间距下进行计算,从而提高计算效率.在3-D非均匀介质中的地震波传播数值模拟结果表明,该方法是一种研究复杂非均匀介质中地震波传播问题的高效方法.     

A Monte Carlo study of rainfall forecasting with a stochastic model

A procedure for short-term rainfall forecasting in real-time is developed and a study of the role of sampling on forecast ability is conducted. Ground level rainfall fields are forecasted using a stochastic space-time rainfall model in state-space form. Updating of the rainfall field in real-time is accomplished using a distributed parameter Kalman filter to optimally combine measurement information and forecast model estimates. The influence of sampling density on forecast accuracy is evaluated using a series of a simulated rainfall events generated with the same stochastic rainfall model. Sampling was conducted at five different network spatial densities. The results quantify the influence of sampling network density on real-time rainfall field forecasting. Statistical analyses of the rainfall field residuals illustrate improvement in one hour lead time forecasts at higher measurement densities.     

Earthquake Engineering and Engineering Vibration Aims And Scope

正This journal is established by the Institute of Engineering Mechanics(IEM),China Earthquake Administration,to promote scientific exchange between Chinese and foreign scientists and engineers so as to improve the theory and practice of earthquake hazards mitigation,preparedness,and recovery.To accomplish this purpose,the journal aims to attract a balanced number of papers between Chinese and     

Foreword

Destructive earthquakes have caused great damage in China and the United States and collapsing buildings havecaused many deaths and injuries. The field of earthquake engineering studies earthquake hazards, the occurrence ofearthquakes of various magnitudes, the nature of the ground shaking during an earthquake, the vibration of structuresduring earthquakes, the strengthening of existing structures and the design of new structures to be earthquake resistant,and finally, how to cope with earthquake damage and restore a city to normal functioning. Such efforts are in progressin both countries, but unfortunately, the language barrier interferes with the free flow of information between China andthe Untied States. It would be mutually beneficial if some means could be developed to promote the exchangeof information across the Pacific Ocean. This new journal has been established for this purpose and its success willbe an important step in promoting earthquake engineering in China and the United States.     

WORLD ASSOCIATION FOR SEDIMENTATION AND EROSION RESEARCH(WASER)

正President:Giampaolo Di Silvio,Italy Vice Presidents:Ulrich C.E.Zanke,Germany Zhao-yin Wang,China The World Association for Sedimentation and Erosion Research(WASER),inaugurated on Oct.19,2004,is an independent non-governmental,non-profit organization.The mission of WASER is to promote international co-operation on the study     

Enhancing remote sensing research on global change to improve our understanding on Earth system processes

正Global Change includes climate change and other environmental changes caused by the joint interaction among various layers of Earth. From the positive side, global change provides new opportunities to human and other living forms on Earth. In the meantime, it creates tremendous challenges and negative impact. At present, the negative impacts have reached all primary processes of the global ecosystem and every aspect of human society, especially causing degradation of the ecosystem. For instance, intensive deforestation causes decline of biodiversity; global warming causes sea level rise and increases     

SCIENCE CHINA Earth Sciences SUBJECT INDEX TO VOLUME 57(2014)

正~~