添加正则化项的黏声逆时偏移成像方法研究

在实际地球介质中传播的地震波会产生衰减和频散现象,因此其更接近于黏弹性介质,在地震处理中补偿这些黏性影响是十分必要的。基于波动方程的叠前深度偏移中进行吸收衰减补偿更准确,也更有物理意义,直接求解双程波动方程的逆时偏移(RTM)能够成像大倾角复杂构造,具有诸多优势。然而当考虑吸收衰减补偿时通常会产生不稳定现象,大部分研究都是在逆时偏移的波场延拓中进行波数域的低通滤波来解决这个问题。本文采用广义标准线性固体的黏声波动方程进行吸收衰减补偿的Q--RTM方法,通过添加正则化项的方式来稳定延拓过程。添加正则化项本质上是低通滤波,滤波窗口是指数形式,在时空域有明确的形式,可以阻止发生高频不稳定。与直接滤波相比,正则化参数可以是空变的,因此比较适合剧烈变化的区域,我们还发现震源归一化互相关成像条件更适合Q--RTM方法。      

<Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>9 Tohoku earthquake in 2011 in Japan: precursors uncovered by natural time analysis

This paper reviews the precursory phenomena of the 2011 M _W9 Tohoku earthquake in Japan that emerge solely when we analyze the seismicity data in a new time domain termed natural time. If we do not consider this analysis, important precursory changes cannot be identified and hence are missed. Natural time analysis has the privilege that enables the introduction of an order parameter of seismicity. In this frame, we find that the fluctuations of this parameter exhibit an unprecedented characteristic change, i.e., an evident minimum, approximately two months before Tohoku earthquake, which strikingly is almost simultaneous with unique anomalous geomagnetic field variations recorded mainly on the z component. This is consistent with our finding that such a characteristic change in seismicity appears when a seismic electric signal (SES) activity of the VAN method (from the initials of Varotsos, Alexopoulos, Nomicos) initiates, and provides a direct confirmation of the physical interconnection between SES and seismicity.     

利用FG5绝对重力仪观测数据确定重力垂直梯度

本文建立了利用FG5实测数据求解重力垂直梯度的数据处理模型与算法.通过对多次自由落体实验的下落距离拟合残差叠加求均值,发现下落距离观测量中存在明显的有色噪声.通过对有色噪声的建模,并以剩余残差为依据选取可靠的下落时段,解算测站点的重力垂直梯度.利用本文所提出的数据处理方法分别对FG5-214绝对重力仪在两个测站上的观测数据进行处理,以相对重力仪测量的重力垂直梯度结果为参考值,本文处理得到的重力垂直梯度结果相比于未考虑有色噪声并依据经验选取下落时段的解算方法得到了显著改善.     

基于广义旋转法双九点格式标量波方程数值模拟

前人虽然基于传统旋转法提出了四阶精度最优17点差分格式,用于提高频率域地震波场数值模拟精度,但其仅适用于等网格间距的情形.这大大限制了该格式的使用范围.为了进一步提高17点有限差分格式的数值精度并将其推广到矩形网格,本文基于二阶精度有限差分算子利用广义旋转法提出了双九点格式,由于其差分格点与前人17点格式在分布上一样,所以也可称为二阶精度最优17点格式,但由差分格式构造原理上来讲,称其为双九点格式更妥.前人四阶精度17点格式仅为本文等网格间距情形时的特殊情况,本文方法单位波长网格点数仅需要2.2个即可.本文格式在继承传统旋转法良好几何旋转性质的同时,拥有平均导数方法适用于矩形网格的特点,和平均导数法所得的广义17点格式相比,本文格式数值精度更高,数值频散抑制性能和差分算子对称性更好.同时,本文双九点格式方法和思想对于后人借助传统九点格式的构造方法将其扩展到17点格式求解各类波动方程具有十分重要的意义.     

郯庐断裂带合肥段五河—合肥断裂构造特征

五河—合肥断裂是郯庐断裂带的西边界断裂,该断裂穿过合肥市城区,是1条规模较大、切割较深的隐伏活动断裂.为了研究该断裂的浅部结构特征、空间展布以及断裂活动性,我们利用2015年在合肥盆地完成的深地震反射剖面数据,采用初至波层析成像方法得到了郯庐断裂带合肥段的浅层P波速度结构和构造形态;考虑到仅根据速度结构剖面还难以确定断裂的准确位置、断层上断点埋深、断层的近地表构造组合样式等特征,研究中跨五河—合肥断裂还完成了2条高分辨率的浅层地震反射剖面.研究结果表明：郯庐断裂带合肥段是一个由多条主干断裂构成的复杂构造带,近地表速度结构表现为凹隆相间的构造特征,且沉积盖层厚度明显受到郯庐断裂带分支断裂的影响和控制.五河—合肥断裂在P波速度结构剖面表现为高速和低速区的分界,对断裂两侧的地层沉积具有重要的控制作用,该断裂向下错断了盆地基底,向上错断了埋深21~35 m的中更新统下部地层,其最新活动时代为中更新世早期.研究结果不仅为进一步认识五河—合肥断裂浅部构造形态提供了地震学依据,还可为该区断裂两侧的城镇规划和建设中避让活动断层提供基础资料.     

多通道瞬变电磁法频率域比值算法及视电阻率计算研究

视电阻率在电磁法勘探中起重要作用,可以作为初步解释结果,多通道瞬变电磁法通过大地脉冲响应的峰值时刻定义视电阻率.鉴于传统方法计算峰值时刻视电阻率需要恢复完整大地脉冲响应,本文通过研究大地脉冲响应及其频谱特征,提出一种新的频率域比值法.利用参考谱及其与一般谱的伸缩因子计算一般谱的峰值时刻,进而得到视电阻率.该方法只需大地脉冲响应单个频点幅度谱,从而不用估计整个大地脉冲响应.实际计算时,可以用大地脉冲响应多个频点幅度谱得到平均视电阻率.模拟数据与实测数据结果表明,在一维介质中,频率域比值法可以得到合理的视电阻率.     

使用宽频带地震记录研究2017年6月24日茂县新磨滑坡的动力学过程

北京时间2017年6月24日5时39分左右,四川省茂县叠溪镇新磨村发生大型岩质滑坡.体积约4.3×10⁶ m³的巨型岩体从山顶脱落,顺坡滑行约2.6 km后破碎沉积;碎屑物掩埋了整个新磨村,造成了巨大的人员伤亡和财产损失.本文使用来自滑坡周围的10个地震台站的宽频带观测资料的长周期信号反演了这次滑坡的受力时间函数;同时使用逐步细化的格点搜索方法得到了滑坡的位置,与其真实位置一致;根据反演的受力时间函数计算了滑坡过程中滑体的运动学参数,得到的滑体运动轨迹与实际路径吻合.综合分析地震信号、受力时间函数和运动学参数表明,本次滑坡主运动的持续时间约为79 s;脱落岩体在5∶38∶50.2启动后持续加速,在5∶39∶37.2达到速度峰值,约为52.1 m·s^-1;这段时间内岩体没有明显的破碎;之后,岩体开始铲刮并裹挟古滑坡造成的碎屑沉积物,自身也开始破碎解体,总体开始减速运动,直到5∶40∶9.2主运动停止;此后,小规模的碎屑散落又持续了约10 s的时间.     

Feasibility of signal enhancement with multiple grounded‐wire sources for a frequency‐domain electromagnetic survey

Frequency‐domain electromagnetic methods with a grounded‐wire source are powerful tools in geophysical exploration. However, the signal may be too weak to guarantee the quality of survey data in complex electromagnetic environments, especially when the receiver is located in the air for the newly developed grounded‐source airborne frequency‐domain electromagnetic method. In this paper, a signal enhancement method with multiple sources is proposed to solve this problem. To evaluate the signal enhancement effect, we compared the signals generated by a single source and multiple sources with equal electric moment. The signal differences caused by synchronisation error and separation distance between source elements were analysed, and the methods to achieve maximum signal were introduced. Besides, we discussed the interaction between adjacent source elements to ensure the system safety, including the changes in output current and the safe distance between two sources using a dual‐source model. Lastly, a comprehensive field experiment was designed and conducted to test the multiple‐source method. The data processing results are comparable for single and dual sources, and the signal‐to‐noise ratio of dual source is higher in the field test. The subsurface resistivity structure at the test site is consistent with the previous controlled‐source audio‐frequency magnetotellurics method. These results show that signal enhancement with multiple sources is feasible. This study provides guidance to the application of multiple sources in field surveys when the survey environment is complex and rigorous.     

Development of a distributed hydrological model to facilitate watershed management

To facilitate precise and cost-effective watershed management, a simple yet spatially and temporally distributed hydrological model (DHM-WM) was developed. The DHM-WM is based on the Mishra-Singh version of the curve number method, with several modifications: The spatial distribution of soil moisture was considered in moisture updating; the travel time of surface runoff was calculated on a grid cell basis for routing; a simple tile flow module was included as an option. The DHM-WM was tested on a tile-drained agricultural watershed in Indiana, USA. The model with the tile flow module performed well in the study area, providing a balanced water budget and reasonable flow partitioning. The daily coefficient of determination and Nash-Sutcliffe coefficient were 0.58 and 0.56, for the calibration period, and 0.63 and 0.62 for the validation period. The DHM-WM also provides detailed information about the source areas of flow components, the travel time and pathways of surface runoff.

EDITOR A. Castellarin; ASSOCIATE EDITOR F.-J. Chang