瞬变电磁法正演计算进展

详细介绍了瞬变电磁法正演计算的方法、现状和发展趋势.瞬变电磁法一维正演计算需要将电磁场从频率域转换至时间域,转换方法有三种,分别是Gaver-Stehfest算法、余弦变换和Guptasarma算法.在这三种方法中,使用较多的是Gaver-Stehfest算法和余弦变换,Gaver-Stehfest算法速度较快,但精度不及余弦变换.瞬变电磁法的数值模拟主要集中于2.5维和三维,使用的数值计算方法有积分方程法、有限差分法、有限单元法和SLDM法.积分方程法主要在三维数值模拟中使用,现已很少使用;有限差分法和有限单元法是目前瞬变电磁法2.5维和三维数值模拟的主要方法;SLDM法主要应用于三维数值模拟.我国瞬变电磁法正演计算成果主要集中在回线源激发的瞬变电磁场一维数值计算和利用有限单元法进行2.5维和三维数值模拟.瞬变电磁法正演计算的发展趋势有:数值算法的改进、提高计算效率和研究地形对瞬变电磁场的影响规律.     

用于波场成像的谱法LU分解

地震波场模拟和偏移成像等有限差分隐格式算法中的重要环节，是实现亥姆霍兹算子表示矩阵H的快速求逆运算. 在螺旋边界条件下，H具有Toeplitz结构的正定厄密矩阵，其快速求逆可由谱法LU分解实现. 本文分析了谱法LU分解对提高计算速度的原理及特点，并着重讨论了在不同类型的介质模型中，采用谱法分解矩阵H时带来的数值误差、误差的分布及其对波场计算的影响. 研究结果表明，对均匀介质而言，矩阵H各列具有相同的非零元素分布，谱法LU分解的误差在吸收边界条件下，不影响波场模拟和成像计算；但对于非均匀介质模型，矩阵H各列具有不同的非零元素分布，谱法LU分解的误差随介质不均匀性程度的增大而增大，势必影响非均匀介质中波场计算. 在波场模拟和成像等有限差分隐格式算法中，采用谱法LU分解完成矩阵求逆时，必须考虑到并尽量减少该方法的误差对波场计算的影响.     

基于波动方程理论的地震波场数值模拟方法综述

地震波场数值模拟的重要性日益凸显,其在地震勘探、各向异性介质属性研究、强地面振动预测、理论地震图合成等方面的应用日渐增多.基于波动方程理论的模拟方法有很多,如:有限差分法、伪谱法、有限元法、谱元法、边界元法等,而现有的综述性文献仅仅是关注某一种或是几种算法,还未有融入最新的文献后系统论述以上五种方法的综述性文献.鉴于此...     

复杂介质地震波传播的褶积微分算子数值模拟h

将佛尔塞（Forsyte）广义正交多项式微分算子地震波场正演模拟算法运用于复杂非均匀介质模型的波场数值模拟中,并比较了该方法在计算效率和计算精度方面与有限差分方法和伪谱法的差异. 数值结果表明,这种广义正交多项式微分算子法计算速度快、精度高,对计算资源需求低,是一种颇具潜力的数值模拟方法.      

2022年1月8日青海门源6.9级强地震动模拟研究

为满足数字化时代应急处置的新要求,基于随机有限断层法的发展,利用青海门源地区地下三维速度结构的模型和v_S30数据,通过逐步迭代射线追踪法、格林函数位移解析相位谱和有限断层法的联合计算,以青海门源6.9级地震为例,获得具有地表土层放大效应的强地震动模拟数据,进而绘制出研究区内PGA和烈度分布模拟图;并与实际台站记录的PGA和实地调查烈度结果相比,其烈度区划范围基本一致,同时也验证此联合计算方法可用于未来地震灾害的快速评估,并为灾后应急救援提供参考。     

用模式匹配法计算电阻率测井响应研究进展

数值模式匹配法是地球物理电磁场计算中的一种半数值、半解析的高效算法,相对于别的数值算法具有明显的优越性,近年来在电法测井数值模拟中得到了成功应用.本文详细地介绍了数值模式匹配法的原理和模式匹配法在电法测井的发展历程和最新应用情况,系统总结了用模式匹配法计算电阻率测井响应的思路和基本步骤,并展望了其在电法测井三维数值模拟中的应用前景.     

基于电场Helmholtz方程的回线源瞬变电磁法三维正演

正演是电磁法勘探野外工作参数选取、室内资料处理与解释的基础,精确、稳定、高效的三维正演算法尤为重要.本文采取先求解拉普拉斯域电场、再由Gaver-Stehfest算法获得时间域磁场的思路,基于电场异常场Helmholtz方程实现了交错网格有限差分法和有限体积法对回线源瞬变电磁法的三维正演.通过对比低阻块状体的积分方程法、时域有限差分法、矢量有限单元法和SLDM法的数值解,验证了交错网格有限差分法和有限体积法的正确性.由于交错网格有限差分法、有限体积法和基于矩形块单元的矢量有限单元法将待求电场均定义在矩形块单元棱边上,因此三种数值算法可采用相同方法进行电场待求量编码、计算背景场和后处理.然而,与矢量有限单元法相比,交错网格有限差分法和有限体积法的系数矩阵更加稀疏,求解效率更高.通过对水平低阻板状体三维模型的数值模拟,我们发现本研究中交错网格有限差分法比有限体积法精度更高;再利用一维解析法求解相应三层层状地电模型的感应电动势,我们还发现两种数值算法和一维解析法计算的感应电动势等值线形状吻合程度高,只是数值范围略有差异.     

BFO-PSO算法下的弹性波数值模拟

地震波正演模拟是地震反演与成像的基础和关键,有限差分算法广泛应用于地震波数值模拟,差分算子的精度直接影响数值模拟的质量和效率.本文提出一种BFO-PSO算法下的有限差分算子优化方法,并应用其进行弹性波数值模拟.首先,将BFO算法中的趋化、复制、驱散三个步骤引入PSO算法,形成具有更好全局搜索能力和更快收敛速度的BFO-PSO混合优化算法;之后构造包含有限差分系数的目标函数,并应用BFO-PSO混合优化算法求取最优解,获得优化的有限差分算子;最后应用此优化的有限差分算子在不同模型上进行弹性波数值模拟.根据频散曲线及数值模拟结果,可以分析得出,BFO-PSO算法优化后的有限差分算子在保证计算效率的同时,具有更高的精度,可以有效压制数值频散,提高数值模拟的精度和效率.     

基于广义正交多项式褶积微分算子的地震波场数值模拟方法

地震波场模拟方法研究对于与波动现象有关的地震学问题的重要性是不言而喻的.就目前现有的各种正演算法来说,精度较高的算法（如有限元法、谱元法、高阶有限差分法等）,其计算速度较慢;计算速度较快的算法（如低阶有限差分法、付氏伪谱法等）计算精度却比较低.为了兼顾地震波场模拟的精度与速度,本文推出了一种快速的、高精度地震波场模拟方法（基于Forsyte广义正交多项式的褶积微分算子法）,该方法是以计算数学中的Forsyte广义正交多项式插值函数为基础,构建一个新的褶积微分算子,并将该算子引入到地震波动方程的一阶速度-应力方程的空间微分运算中去,采用时间交错网格有限差分算子替代普通的差分算子以匹配高精度的褶积微分算子,从而构造一种全新的地震波场数值模拟方法.该方法同时具有广义正交多项式方法的高精度和短算子低阶有限差分算法的高速度.通过对算子长度的调节及算子系数的优化,可同时兼顾波场解的全局信息与局部信息.复杂非均匀介质模型中的波场数值模拟实验证实了该方法的可行性及优越性.     

复杂介质多尺度变网格地震波数值模拟

在应用有限差分法地震波场数值模拟过程中,模拟精度和计算效率是关键问题之一,特别是对于速度变化剧烈的含小尺度地质体的数值模拟尤为重要.为了既能精细刻画介质的局部结构,又能保证模拟的时效性,本文在传统交错网格有限差分算法的基础上,将变网格有限差分方法引入到交错网格高阶差分数值模拟中,对交错网格空间算法进行了改进,并避免了因插值因素影响模拟精度和计算效率.本文采用完全匹配层吸收边界条件实现对可变网格数值模拟的边界处理,根据完全匹配层吸收边界原理和参量的分裂思想,在不同的方向分别加相应的阻尼衰减因子,解决边界的反射问题.此外还分析了数值频散对于模拟结果的影响.文中分别设计了地堑和生物礁两个介质模型对可变交错网格算法的有效性和稳定性进行了验证.数值模拟结果表明,可变交错网格有限差分算法数值模拟能够对地下介质物性的空间变化进行准确刻画,进一步增强数值模拟对复杂介质的适应性,同时为地震数据的波场成像、纵横波联合解释等提供可靠依据.     

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)

正~~