石家庄地震前兆台网月评价系统

基于VB语言,利用C/S架构,编写石家庄地震前兆台网月评价系统,设置不同功能单元,开展仪器运行监控、台网日常运维监控,通过界面选择相应日期,读取省前兆台网中心发布的《运行监控日报》等Excel文件,使用input box函数等在文本框中输入、统计、计算,并显示地震前兆台网月评价成绩及台网运行情况。     

甘孜-玉树断裂带晚第四纪走滑速率与滑动分解作用

滑动速率是研究断裂运动学特征、地震活动性和区域应变分配的重要参数和依据。前人关于甘孜-玉树断裂带滑动速率的研究结果存在较大差异,因此,其晚第四纪滑动速率有待进一步调查研究。本文基于卫星影像解译和野外实地考察,对甘孜-玉树断裂带西段（玉树断裂）上典型断错地貌点进行测量分析,得到玉树断裂晚第四纪走滑速率为6.6±0.1-7.4±1.2mm/a。通过与前人对甘孜-玉树断裂带东段（甘孜断裂）滑动速率的研究结果进行对比,发现甘孜-玉树断裂带东、西段滑动速率不一致,其原因是甘孜断裂的左旋滑移在向西传递的过程中,一部分应变被分配到了巴塘盆地南缘断裂上。巴塘盆地南缘断裂的存在很好地解释了玉树断裂的走滑速率比甘孜断裂偏低的原因。但是,从区域变形来看,巴塘盆地南缘断裂分配的滑动速率恰好说明了甘孜-玉树断裂带东、西段及鲜水河断裂带的水平构造变形是协调一致的。     

宽频带超导重力仪的研制进展

本文报道一种构建宽频带高分辨率超导重力仪（SG）的方法.与已商业化的GWR超导重力仪类似,新型超导重力仪同样使用磁悬浮超导检验质量来构建弹簧振子结构,不同的是弹簧振子的固有频率更高,可以有效覆盖地脉动频带.同时,利用基于超导量子干涉器件（SQUID）的位移传感技术,大幅提高了位移传感灵敏度,以保证仪器在测量带宽扩展之后仍具有与GWR仪器相当的测量分辨率,从而实现对时变重力信号与地球背景噪声信号的同时高分辨观测.宽频带超导重力仪核心部件垂向超导加速度计（VSA）的设计、组装和初步测试已经完成.结果表明,加速度计的噪声本底为8 μGal/Hz^1/2@0.1 Hz,实现了对地脉动信号的高分辨率测量.而为了将仪器噪声本底降低到GWR仪器水平,加速度计的固有频率需要进一步减小,该部分参数优化工作正在进行.本文将对垂向超导加速度计的工作原理、结构和测试结果进行详细介绍.     

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

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

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

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

Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.     

Effect of compressive strength on the performance of the NEMO-LIM model in Arctic Sea ice simulation

Satellite records show that the extent and thickness of sea ice in the Arctic Ocean have significantly decreased since the early 1970s. The prediction of sea ice is highly important, but accurate simulation of sea ice variations remains highly challenging. For improving model performance, sensitivity experiments were conducted using the coupled ocean and sea ice model (NEMO-LIM), and the simulation results were compared against satellite observations. Moreover, the contribution ratios of dynamic and thermodynamic processes to sea ice variations were analyzed. The results show that the performance of the model in reconstructing the spatial distribution of Arctic sea ice is highly sensitive to ice strength decay constant (C^rhg). By reducing the C^rhg constant, the sea ice compressive strength increases, leading to improved simulated sea ice states. The contribution of thermodynamic processes to sea ice melting was reduced due to less deformation and fracture of sea ice with increased compressive strength. Meanwhile, dynamic processes constrained more sea ice to the central Arctic Ocean and contributed to the increases in ice concentration, reducing the simulation bias in the central Arctic Ocean in summer. The root mean square error (RMSE) between modeled and the CryoSat-2/SMOS satellite observed ice thickness was reduced in the compressive strength-enhanced model solution. The ice thickness, especially of multiyear thick ice, was also reduced and matched with the satellite observation better in the freezing season. These provide an essential foundation on exploring the response of the marine ecosystem and biogeochemical cycling to sea ice changes.

     

Numerical Analysis of Groundwater Ridging Processes Considering Water‐Air Flow in a Hillslope

In this study, a water‐air two‐phase flow model was employed to investigate the formation, extension, and dissipation of groundwater ridging induced by recharge events in a hypothetical hillslope‐riparian zone, considering interactions between the liquid and gas phases in soil voids. The simulation results show that, after a rain begins, the groundwater table near the stream is elevated instantaneously and significantly, thereby generating a pressure gradient driving water toward both the stream (the discharge of groundwater to the stream) and upslope (the extension of groundwater ridging into upslope). Meanwhile, the airflow upslope triggered by the advancing wetting front moves downward gradually. Therefore, the extension of groundwater ridging into upslope and the downward airflow interact within a certain region. After the rain stops, groundwater ridging near the stream declines quickly while the airflow in the lower part of upslope is still moving into the hillslope. Thus, the airflow upslope mitigates the dissipation of groundwater ridging. Additionally, the development of groundwater ridging under different conditions, including rain intensity, intrinsic permeability, capillary fringe height, and initial groundwater table, was analyzed. Changes in intrinsic permeability affect the magnitude of groundwater ridging near the stream, as well as the downward speed of airflow, thereby generating highly complex responses. The capillary fringe is not a controlling factor but an influence factor on the formation of groundwater ridging, which is mainly related to the antecedent moisture. It was demonstrated that groundwater ridging also occurs where an unsaturated zone occurs above the capillary fringe with a subsurface lateral flow.     

强震临震微波动现象初步研究(二)

为了验证和检验强震临震微波动现象及其主要特征,利用甘、青、川、滇、藏等区域台网200余个宽频带数字地震台站的波形资料,建立了实时监控技术系统。通过对2012~2014年间青藏块体发生的24次5级以上地震及部分无震区域台站观测的全程实时跟踪、动态监控,检验和验证了临震微波动现象的重现性和客观性,验证了临震微波动现象的频谱、时间、空间、方向性等特征,进一步得出了临震微波动现象可能与地震强度有关,且震级越大,震前出现持续临震微波动的可能性越大的结论。     

唐山马家沟矿井水位异常分析

地下水动态受水文因素影响较大,对地震和构造活动具有较灵敏的响应。判别并排除各种水文干扰,确认地下水在地震前的异常变化,对提高地震分析预报能力,具有重要作用。马家沟矿井水位动态观测层与地下水开采层为同一含水层,井水位于2010年出现破年变异常,加速持续上升,截至2015年,最大上升幅度约30 m。依据该井水文地质环境特征,根据唐山市区2001-2015年地下水位、降雨量、地下水开采量实测资料,建立合理多元回归模型和三维地下水流动模型,发现地下水开采量减少应为影响马家沟矿井水位动态的可能因素。文中采取的异常识别与分析方法,可为其他类似井孔的地下水动态异常识别及判定提供一定借鉴。