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

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

南沙海域海底浅层第四系地震地层结构与地质意义

科技进步使得地震采集数据量及其精度实现质的飞跃,南沙海域的地震调查,采用国际先进水平的等离子震源,获得地震剖面精度优于3 m局部甚至可在1 m之内,为科学研究、工程、浅层资源、地质灾害预警等研究奠定基础.在南沙第四系识别出了七个地震反射界面,确定500 ms至海底沉积区间内地震层序与三期海平面升降并与冰期、间冰期对应,典型地震剖面展示第四纪至少有三套完整的具有三角洲顶积层、前积层、底积层的地震相证据,与全球的第四纪海平面的变化一致.精细的浅层地震结构表明：南沙的地质现象丰富,浅层断层非常发育,下切河谷、泥石流、滑坡等地质遗迹的形成的地震相清晰,是地质活动频繁的地区.     

地震油气储层的小样本卷积神经网络学习与预测

地震储层预测是油气勘探的重要组成部分,但完成该项工作往往需要经历多个环节,而多工序或长周期的研究分析降低了勘探效率.基于油气藏分布规律及其在地震响应上所具有的特点,本文引入卷积神经网络深度学习方法,用于智能提取、分类并识别地震油气特征.卷积神经网络所具有的强适用性、强泛化能力,使之可以在小样本条件下,对未解释地震数据体进行全局优化提取特征并加以分类,即利用有限的已知含油气井段信息构建卷积核,以地震数据为驱动,借助卷积神经网络提取、识别蕴藏其中的地震油气特征.将本方案应用于模型数据及实际数据的验算,取得了预期效果.通过与实际钻井信息及基于多波地震数据机器学习所预测结果对比,本方案利用实际数据所演算结果与实际情况有较高的吻合度.表明本方案具有一定的可行性,为缩短相关环节的周期提供了一种新的途径.     

一种基于天基的瞬变电磁高分辨探测小尺度太空碎片的方法

目前,太空碎片的天基探测方法一般是激光雷达法.由于激光雷达的辐射频率一般为可见光或红外波段,在分辨厘米量级的太空碎片时具有一定的优势.但激光雷达在探测时也面临一些客观问题：即（1）难以捕捉快速移动的太空碎片;（2）对目标体穿透能力低;（3）外层空间的强干扰环境对激光雷达的不良影响等.这些问题在一定程度上限制了激光雷达在太空碎片探测方面的应用.本文针对太空碎片体积小、预警距离远和强干扰环境等特点,提出使用高性能瞬变电磁辐射源对远距离、小尺度的太空碎片进行探测.通过三维矢量有限元方法,分析了利用高性能瞬变电磁辐射源的优点,对比不同瞬变电磁辐射源的幅频特性,从发射机理上证明高性能瞬变电磁辐射源在辐射能力与频率带宽方面都优于传统瞬变电磁辐射源.并通过调整脉冲宽度,得到最适合探测本文太空碎片模型的辐射脉冲.最后通过电场分布图和多测道图对太空碎片的明显电场分异结果,证实了高性能瞬变电磁辐射源在探测远距离、小尺度太空碎片方面的有效性.     

Piping dynamics in mid‐altitude mountains under a temperate climate: Bieszczady Mountains,eastern Carpathians

Piping has been recognized as an important geomorphic, soil erosion and hydrologic process. It seems that it is far more widespread than it has often been supposed. However, our knowledge about piping dynamics and its quantification currently relies on a limited number of data for mainly loess‐derived areas and marl badlands. Therefore, this research aimed to recognize piping dynamics in mid‐altitude mountains under a temperate climate, where piping occurs in Cambisols, not previously considered as piping‐prone soils. It has been expressed by the estimation of erosion rates due to piping and elongation of pipes in the Bere?nica Wy?na catchment in the Bieszczady Mountains, eastern Carpathians (305 ha, 188 collapsed pipes). The research was based on the monitoring of selected piping systems (1971–1974, 2013–2016). Changes in soil loss vary significantly between different years (up to 27.36 t ha^?1 yr^?1), as well as between the mean short‐term erosion rate (up to 13.10 t ha^?1 yr^?1), and the long‐term (45 years) mean of 1.34 t ha^?1 yr^?1. The elongation of pipes also differs, from no changes to 36 m during one year. The mean total soil loss is 48.8 t ha^?1 in plots, whereas in the whole studied catchment it is 2.0 t ha^?1. Hence, piping is both spatially and temporally dependent. The magnitude of piping in the study area is at least three orders of magnitude higher than surface erosion rates (i.e. sheet and rill erosion) under similar land use (grasslands), and it is comparable to the magnitude of surface soil erosion on arable lands. It means that piping constitutes a significant environmental problem and, wherever it occurs, it is an important, or even the main, sediment source. Copyright © 2017 John Wiley & Sons, Ltd.     

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     

Enhanced identification of a hydrologic model using streamflow and satellite water storage data: A multicriteria sensitivity analysis and optimization approach

Hydrologic model development and calibration have continued in most cases to focus only on accurately reproducing streamflows. However, complex models, for example, the so‐called physically based models, possess large degrees of freedom that, if not constrained properly, may lead to poor model performance when used for prediction. We argue that constraining a model to represent streamflow, which is an integrated resultant of many factors across the watershed, is necessary but by no means sufficient to develop a high‐fidelity model. To address this problem, we develop a framework to utilize the Gravity Recovery and Climate Experiment's (GRACE) total water storage anomaly data as a supplement to streamflows for model calibration, in a multiobjective setting. The VARS method (Variogram Analysis of Response Surfaces) for global sensitivity analysis is used to understand the model behaviour with respect to streamflow and GRACE data, and the BORG multiobjective optimization method is applied for model calibration. Two subbasins of the Saskatchewan River Basin in Western Canada are used as a case study. Results show that the developed framework is superior to the conventional approach of calibration only to streamflows, even when multiple streamflow‐based error functions are simultaneously minimized. It is shown that a range of (possibly false) system trajectories in state variable space can lead to similar (acceptable) model responses. This observation has significant implications for land‐surface and hydrologic model development and, if not addressed properly, may undermine the credibility of the model in prediction. The framework effectively constrains the model behaviour (by constraining posterior parameter space) and results in more credible representation of hydrology across the watershed.     

2016年台湾美浓MW6.4地震震源参数的InSAR和GPS反演

2016年2月6日台湾西南部高雄市美浓区发生了M_W6.4地震.本文结合ALOS2卫星升降轨、Sentinel-1A升轨SAR数据,采用两轨差分干涉技术获取了该区域的同震形变场,形变结果表明震中西北部以抬升为主,最大视线向形变量约为11.2 cm.基于均匀位错模型和多峰值粒子群（MPSO）算法,利用InSAR和GPS形变数据联合反演了美浓地震的断层几何参数,结果表明震源中心位于22.920°N,120.420°E,深度约12 km,发震断层长度约15 km,走向角307°,倾角16.5°,平均滑动角为51.5°,此次地震是以逆冲倾滑兼左旋走滑的破裂模式.利用格网迭代搜索法得到最优倾角为15.7°,GPS和InSAR最优权比为18：1,最优平滑因子为0.06.基于非均匀位错模型,利用非负最小二乘方法进行线性反演,结果显示最大倾滑和走滑量分别为51.7 cm和55.3 cm,对应矩震级为M_W6.38,略小于GCMT （M_W6.4）的结果.通过与已有文献的比较和对该区域断层构造的分析,发现美浓地震的发震断层为单一断层的解释更为合理,我们推测发震断层是位于左镇、后甲里等断层之间的一条东南-西北走向往东北倾斜的盲断层,并初步推测2010年M_W6.3甲仙地震也同该断层有关.     

几种典型近地表有限导体瞬变电磁特征响应研究

近地表有限导体的瞬变电磁响应中包含其形状、结构特征,如何建立这些特征与其响应之间的关系是近地表有限导体探测与识别的核心问题.本文建立了几种典型有限导体响应的三维感应偶极子模型,通过有限导体早期幂函数衰减的幅度与速率、晚期指数衰减的幅度与时间常数共四个模型参数表征有限导体的形状和结构特征,利用最小二乘拟合算法提取这些模型参数,揭示有限导体的形状、结构特征与模型参数之间的内在联系.理论与实验研究结果表明:具有绕球心旋转对称性的球形目标体,其三维感应偶极子响应完全相同,早期响应呈幂函数衰减,晚期呈e指数衰减,信号幅度和时间常数均能很好反映球体半径与材料特性;具有绕轴旋转对称性目标体,其三维感应偶极子响应部分相同,早期幂函数衰减的幅度与速率主要体现目标体的长度特征,晚期e指数衰减的时间常数体现了目标体的外径与壁厚特征;非轴对称性目标体,其三维感应偶极子响应均不相同,各个感应偶极子参数差异较大.本文的研究对提高瞬变电磁系统区分轴对称目标体与金属碎片能力,降低系统探测虚警率具有重要意义.     

Spatial behaviour of Rayleigh waves in layered half‐spaces under active surface sources

In the free state, Rayleigh waves are assumed to travel in the form of planar wavefronts. Under such an assumption, the propagation behaviour of the modes of Rayleigh waves in layered half‐spaces is only frequency dependent. The frequency behaviour, which is often termed as dispersion, is determined by the shear wave velocity profile of layered soils within the depth related to wavelength (or frequency). According to this characteristic, the shear wave velocity profile can be back‐analysed from the dispersion. The technique is widely used in the surface wave testing. However, the wavefronts of Rayleigh waves activated by the surface sources are non‐planar. The geometric discrepancy could result in Rayleigh waves manifesting distance‐dependent behaviour, which is referred to as spatial behaviour in this paper. Conventional analysis ignoring this spatial behaviour could introduce unexpected errors. In order to take the effects of sources on the propagation behaviour into account, a new mathematical model is established for Rayleigh waves in layered elastic media under vertical disc‐like surface sources using the thin‐layer method. The spatial behaviour of the activated modes and the apparent phase velocity, which is the propagation velocity of Rayleigh waves superposed by the multiple modes, are then analysed. Aspects of the spatial behaviour investigated in this paper include the equilibrium path, the particle orbit, and the geometric attenuation of the activated Rayleigh waves. The results presented in this paper can provide some guidelines for developing new inverse mathematical models and algorithms.