二维速度结构走时反演RAYINVR算法试验

为分析RAYINVR算法对成像结果的影响,分别从初始模型、观测数据、阻尼因子和识别误差4个方面进行数值试验。结果表明,该程序能够进行快速有效的反演,初始模型、观测数据、阻尼因子、数据不确定性等对反演结果和分辨率有一定的影响。     

Finite-difference calculation of traveltimes based on rectangular grid

Introduction The calculation of seismic wave traveltimes is a basic and the most important step in tomo-graphy, seismic wave forward modeling and Kirchhoff prestack depth migration. Limitations withtraditional ray tracing fall into four categories. a) Analytical methods can only realize ray tracingfor simply varying velocity fields, so they have relative small applied-range; b) Shooting methodsof ray tracing can cause shadow zones. When the shadow zones exist the method will invalid; c)…     

1998年夏第二阶段梅雨期乌拉尔山阻塞形势的维持

１９９８年夏季长江流域发生了近５０年来最严重的洪水。洪水形成最直接、最主要的原因是梅雨异常。异常梅雨的形成与东亚夏季风偏弱及热带外环流持续异常有关，其中一个明显特征是乌拉尔山地区长时间维持阻塞形势。本文结合诊断分析和数值试验，从瞬变对基本流的强迫（大气内部强迫）及热带热源强迫（外源强迫）两方面，分析了与第二段梅雨相对应的乌拉尔长时间阻塞的维持机制。利用共轭敏感性分析方法，计算了最有利于乌拉尔阻塞发展和维持的敏感扰动，发现扰动的分布位置，刚好与观测到的，从东大西洋到欧洲区域的异常增强的瞬变活动区相重叠。Ｅ矢量及斜压线性静止波模式的诊断进一步表明，异常期间的增强瞬变活动有利于乌拉尔出现正高度异常。计算了持续异常期间的高空急流及大气加热场，发现北美到大西洋的高空急流及热带加热都出现明显异常。中期天气预报模式 ＩＡＰ Ｔ４２Ｌ９的集合预报试验表明，热带地区的加热异常，尤其是热带中西太平洋和大西洋的加热异常，有利于乌拉尔正高度异常的形成。最后，提出了一种热带异常热源驱动下，瞬变波与定常波双向相互作用的阻塞形成与自维持的可能机制。     

利用Munk公式和海底反射走时反演深海水体速度:理论

为解决深海水体速度建模问题，建立了两个利用Munk公式和海底反射走时反演深海水体速度的反演理论，一个是射线走时反演理论，另一个是波动方程走时反演理论。具体地讲，是建立了地震波走时关于Munk参数的Fréchet导数和海森矩阵公式。为得到射线走时的Fréchet导数和海森矩阵，利用了3种不同的方法，即变分法、微扰法和Taylor级数展开法。这3种方法在理论上完全等价，其最终结果也完全相同。为得到波动方程走时反演理论中的Fréchet导数和海森矩阵，利用复合函数求导的链式法则对文献中的有关结果进行变换，以将其从面向速度本身的Fréchet导数和海森矩阵变换为面向Munk参数的Fréchet导数和海森矩阵。此外，还提出了实现利用Munk公式和海底反射走时构建深海水体声速剖面的5个策略，包括利用解析公式计算海底反射走时的由下至上的射线追踪策略，选用最简单的海底边界条件把数值模拟区域限定为水体的策略，用已知的深海声道轴深度和交替反演策略将四参数反演问题转简化为单参数反演问题的策略，由点到面，利用一维声速剖面公式实现深海水体的维速度建模的策略，以及用最简单的算法实现走时反演的最优化计算策略等。与常规走时反演相比，利用Munk公式进行的海底反射走时反演最多只需要反演4个参数，并可利用交替反演将四参数反演化为单参数反演，大大地减少了走时反演的计算量。     

叠前逆时深度偏移中的激发时间成像条件

与其他偏移方法相比,逆时偏移基于精确的波动方程而不是对其近似,用时间外推来代替深度外推.因此,它具有良好的精度,不受地下构造倾角和介质横向速度变化的限制.激发时间成像条件的求取是叠前逆时偏移的难点之一,本文采用求解程函方程的方法得到地下各点的初至波走时,以此作为叠前逆时偏移的成像条件.基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法精度较高并适用于强纵横向变速的复杂介质.试算结果表明,在复杂介质模型中利用叠前逆时深度偏移收到了很好的成像效果.     

复杂介质地震定位中震源轨迹的计算

在地震定位中常常需要求解震源轨迹,但由于复杂介质中的震源轨迹较为复杂,难以给出其解析解,因此震源轨迹的计算通常仅限于简单介质模型.本文基于最小走时树射线追踪技术,提出了一种计算复杂介质中震源轨迹的方法.为回避发震时间问题,以观测到时差作为震源轨迹的约束条件.首先从模型节点中选出少量理论到时差与观测到时差之绝对差,即双重时差较小的点作为震源轨迹的代表点,然后以其中双重时差最小的点为初始点,在双重时差场中利用最小走时树射线追踪方法计算出初始点到其他震源轨迹代表点的射线路径作为震源轨迹.当选的震源轨迹代表点较多时,得到的震源轨迹较为粗略,此时可去掉射线经过次数较少的代表点的射线路径使震源轨迹更为精细.为减少计算量,对最小走时树射线追踪方法的终止条件做了修正.以一个复杂介质模型中的地震为例,计算了包括速度扰动、到时扰动等不同情况下的震源轨迹,结果表明所提出的震源轨迹计算方法切实可行.     

Establishment of numerical beach-litter hindcast/forecast models: an application to Goto Islands, Japan

This study attempts to establish a system for hindcasting/forecasting the quantity of litter reaching a beach using an ocean circulation model, a two-way particle tracking model (PTM) to find litter sources, and an inverse method to compute litter outflows at each source. Twelve actual beach survey results, and satellite and forecasted wind data were also used. The quantity of beach litter was hindcasted/forecasted using a forward in-time PTM with the surface currents computed in the ocean circulation model driven by satellite-derived/forecasted wind data. Outflows obtained using the inverse method was given for each source in the model. The time series of the hindcasted/forecasted quantity of beach litter were found consistent with the quantity of beach litter determined from sequential webcam images of the actual beach. The accuracy of the model, however, is reduced drastically by intense winds such as typhoons which disturb drifting litter motion.     

虑传播效应的多波保幅AVO正演(英文)

传统的AVO正演只考虑了单一界面的反射系数对地震波波场振幅的影响,忽略了地震波在介质中传播的各种传播效应。通过引入地震波在介质中传播的几何扩散、吸收衰减以及透射损失等传播效应,提出了基于射线理论的水平层状介质多波保幅AVO正演方法。推导了水平层状介质多波几何扩散校正公式,来描述多波在介质中传播的几何扩散效应。通过直接引入复旅行时,而无需借助复速度,建立了复旅行时与品质因子的关系,来描述粘弹介质的吸收衰减。直接求解Zoeppritz方程计算多波的透射系数,用于描述多波在介质中传播时的透射损失。数值计算表明,几何扩散、吸收衰减以及透射损失对多波振幅的影响是随偏移距变化而变化的,多波保幅AVO正演需要考虑波传播效应对反射波振幅的改造。