中国大陆下扬子及邻区地壳结构基本特征:深地震测深研究综述

系统回顾了20世纪70年代以来在中国大陆下扬子及其邻区开展的深地震测深工程,总结了相关宽角反射/折射地震资料的波组特征及地壳结构的基本特征。下扬子及邻区深地震测深资料普遍具有较清晰的初至震波Pg、上地壳底界面反射波P1、中地壳底界面反射波P3、莫霍界面反射波Pm及上地幔顶部首波Pn,震相连续,易追踪对比。受浅表沉积盖层或基岩等的影响,Pg波常出现局部延迟滞后或走时超前现象。尤其是大别造山带个别炮点的反射能量较弱或同相轴扭曲、波形紊乱等,均与地壳内界面和莫霍面的深度突变或破碎有关。该区域地壳结构大致分为上、中、下3层,但视资料情况中地壳和下地壳又可进一步划分为2个亚层。华北地台和扬子地台地壳厚度30~36km,莫霍界面形态变化较缓,存在局部隆起,下地壳平均速度6.7±0.3km/s。但大别造山带下方地壳厚度32~41km,莫霍界面下凹且出现4~7km垂向错断,下地壳平均速度6.8±0.2km/s。     

Three-dimensional forward modeling and inversion of borehole-to-surface electrical imaging with different power sources

Borehole-to-surface electrical imaging (BSEI) uses a line source and a point source to generate a stable electric field in the ground. In order to study the surface potential of anomalies, three-dimensional forward modeling of point and line sources was conducted by using the finite-difference method and the incomplete Cholesky conjugate gradient (ICCG) method. Then, the damping least square method was used in the 3D inversion of the formation resistivity data. Several geological models were considered in the forward modeling and inversion. The forward modeling results suggest that the potentials generated by the two sources have different surface signatures. The inversion data suggest that the low-resistivity anomaly is outlined better than the high-resistivity anomaly. Moreover, when the point source is under the anomaly, the resistivity anomaly boundaries are better outlined than when using a line source.     

Thermo-hydro-mechanical Modeling of CO2 Sequestration System Around Fault Environment

Geological sequestration of CO₂ (carbon dioxide) shows great potential to reduce Greenhouse gas emissions. However, CO₂ injection into geological formations may give rise to a variety of coupled chemical and physical processes. The thermo-hydro-mechanical (THM) impact of CO₂ injection can induce fault instability, even possibly lead to seismic activities in and around the disposal reservoir. A sequential coupling approach under some assumptions was proposed in the numerical study to investigate the THM behavior of the CO₂ sequestration system concerning the temperature, initial geological stress, injection pressure and CO₂ buoyancy. The fault was treated as a flexible contact model. The effects of CO₂ injection on the mechanical behavior of the faults were investigated. The Drucker-Prager model and the cap model were used to model the constitutive relationship of formations. The numerical results show that injection pressure sensitively affects the relative slip change of the fault. At the initial stage of the sequestration process, the injection pressure plays a key role in affecting the pore pressure of the formations. However, as time continues, the influence of CO₂-induced buoyancy becomes obvious on the pore pressure of the formations. In general, The THM effects of CO₂ geosequestration do not affect the mechanical stability of formations and faults.     

山西裂谷系的深部地球物理探测

综述了山西裂谷系近年来通过重力、工程爆破及测震、大地电磁测深、航磁资料反演、深地震测深等研究该地区地壳与上地幔结构的结果与进展。文中还提出了在本区继续开展地壳与上地幔细结构研究的任务与设想。     

结合高光谱数据反演吉林石头口门水库悬浮物含量和透明度

水中悬浮物含量是评价水环境质量的一个重要参数,它可影响水体透明度、混浊度、水色等光学性质,决定太阳光照在水下的分布和浮游植物对光照的利用,并最终影响水体的初级生产力.本文利用长春市石头口门水库的高光谱实测数据和水质采样分析数据,尝试通过一阶微分法建立悬浮物估测模型,再用估测结果反演透明度信息.结果表明:用590nm处的一阶微分光谱值建立的悬浮物估测模型的决定系数R2和均方根误差(RMSE)分别为0.76、9.09 mg/L,验证模型的R2和RMSE分别为0.79、7.15 mg/L;石头口门水库上游严重的水土流失导致透明度受悬浮物含量影响较大,两者存在明显的指数关系,相关系数r为-0.80,用悬浮物含量建立的透明度估测模型的R2和RMSE分别为0.79和0.12 m,验证模型的R2和RMSE分别为0.82、0.11 m.显著水平均为p＜0.01.研究结果表明,该方法用于石头口门水库悬浮物含量和透明度的定量遥感结果较理想.     

Method for solving dynamic equilibrium equations based on displacement excitation

Mode superposition is a widely used method for solving the dynamic equilibrium equation in structural dynamic analysis. However, the accuracy of this method may be reduced when the dynamic equilibrium equations are set up using displacement excitation. A new method for developing solutions for dynamic equilibrium equations based on displacement excitation is introduced. The dynamic equilibrium equation is decomposed into two parts, namely displacement excitation and velocity excitation, and precise integration and mode superposition methods are combined to solve the equation. Ritz vectors are then used to calculate the static response of the truncated modes of the structure, and a method for determining the number of participating modes is obtained. Using multi-degree-of-freedom systems as two computational examples, the differences in the structural responses obtained from the displacement excitation and acceleration excitation are compared and analyzed. It is shown that the new solution method generates consistent accuracy between the displacement excitation and acceleration excitation.     

3-D multi-parameter type traveltime tomography in a spherical coordinate frame: comparison of double and triple class simultaneous inversions

It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later(primary reflected) arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm(in a spherical coordinate system) with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data(noise free and added noise) indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.     

Numerical simulation of seismic wavefields in TTI media using the rotated staggered-grid compact finite-difference scheme

Numerical simulation in transverse isotropic media with tilted symmetry axis(TTI) using the standard staggered-grid finite-difference scheme(SSG)results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme(RSG) is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme(CSG) is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme(RSGC). The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.     

常用岩土本构模型的选择及对场地动力反应分析结果的影响

合理选择本构模型是土动力学问题数值模拟中的一项重要工作。利用PLAXIS 2D软件的土工实验模拟功能分别对4种常用的岩土本构模型——线弹性模型、摩尔库伦模型、土体硬化模型和小应变土体硬化模型在往复荷载下的理论滞回曲线进行了对比分析,并在此基础上研究了选择不同本构模型对自由场地震反应分析结果的影响以及不同本构模型中各参数的变化对场地动力计算结果的敏感性分析。研究结果为土动力学问题数值模拟中如何选择本构模型和合理判断数值分析结果提供了参考依据。     

Estimation of the Budyko model parameter for small basins in China

At the mean annual scale, water availability of a basin is substantially determined by how much precipitation will be partitioned into evapotranspiration and run-off. The Budyko framework provides a simple but efficient tool to estimate precipitation partitioning at the basin scale. As one form of the Budyko framework, Fu's equation has been widely used to model long-term basin-scale water balance. The major difficulty in applications of Fu's equation is determining how to estimate the curve shape parameter ω efficiently. Previous studies have suggested that the parameter ω is closely related to the long-term vegetation coverage on large river basins globally. However, on small basins, the parameter ω is difficult to estimate due to the diversity of controlling factors. Here, we focused on the estimation of ω for small basins in China. We identified the major factors controlling the basin-specific (calibrated) ω from nine catchment attributes based on a dataset from 206 small basins (≤50,000 km²) across China. Next, we related the calibrated ω to the major factors controlling ω using two statistical models, that is, the multiple linear regression (MLR) model and artificial neural network (ANN) model. We compared and validated the two statistical models using an independent dataset of 80 small basins. The results indicated that in addition to vegetation, other landscape factors (e.g., topography and human activity) need to be considered to capture the variability of ω on small basins better. Contrary to previous findings reached on large basins worldwide, the basin-specific ω and remote sensing-based vegetation greenness index exhibit a significant negative correlation. Compared with the default ω value of 2.6 used in the Budyko curve method, the two statistical models significantly improved the mean annual ET simulations on validation basins by reducing the root mean square error from 114 mm/year to 74.5 mm/year for the MLR model and 70 mm/year for the ANN model. In comparison, the ANN model can provide a better ω estimation than the MLR model.