宁波市地下水位动态与地面沉降预测分析

为了准确预测分析宁波市地下水位动态与地面沉降的发展趋势,建立了宁波市第四纪松散沉积层孔隙地下水流三维数值模拟模型和地面沉降与地下水位多元线性回归模型,预测了2009年底到2020年底的逐月地下水位动态和逐年地面沉降量的变化特征.结果表明,从2013年起,除山区沟谷孔隙潜水地下水位降落漏斗逐渐扩大外,其余孔隙水的地下水流场基本趋于稳定,地下水位年际变化很小,年地面沉降量也逐渐变小,由2012年的5.62 mm/a逐渐下降到2020年的5.54 mm/a,由地下水位下降引起的地面沉降基本得到控制.     

Raman Micro-Spectroscopic Study of Sulfate Ion in the System Na2SO4 – H2O

This work aims to quantify sulfate ion concentrations in the system Na₂SO₄-H₂O using Raman micro-spectroscopy.Raman spectra of sodium sulfate solutions with known concentrations were collected at ambient temperature（293 K） and in the 500 cm^₁-4000 cm^-1 spectral region.The results indicate that the intensity of the SO₄^2- band increases with increasing concentrations of sulfate ion.A linear correlation was found between the concentration of SO₄^2-（c） and parameter I₁,which represents the ratio of the area of the SO₄^2- band to that of the O-H stretching band of water（A_s/A_w）:I₁=-0.00102+0.01538 c.Furthermore,we deconvoluted the O-H stretching band of water(2800 cm^-1-3800 cm^-1) at 3232 and 3430 cm^-1 into two sub-Gaussian bands,and then defined Raman intensity of the two sub-bands as A_Bi(3232 cm^-1) and A_B2(3430 cm^-1),defined the full width of half maximum（FWHM） of the two sub-bands as W_B1(3232 cm^-1) and W_B2(3430 cm^-1).A linear correlation between the concentration of SO₄^2-（c） and parameter I₂,which represents the ratio of Raman intensity of SO₄^2-（A_s）(in 981 cm^-1) to(A_B1+A_B2),was also established:I₂=-0.0111+0.3653 c.However,no correlations were found between concentration of SO₄^2-（c） and FWHM ratios,which includes the ratio of FWHM of SO₄^2-（W_s） to W_B1 W_B2 and W_B1+B2(the sum of W_B1 and W_B2),suggesting that FWHM is not suitable for quantitative studies of sulfate solutions with Raman spectroscopy.A comparison of Raman spectroscopic studies of mixed Na₂SO₄ and NaCI solutions with a constant SO₄^2- concentration and variable CI^- concentrations suggest that the I\ parameter is affected by CI^-,whereas the I₂ parameter was not.Therefore,even if the solution is not purely Na₂SO₄-H₂O,SO₄^2- concentrations can still be calculated from the Raman spectra if the H₂O band is deconvoluted into two sub-bands,making this method potentially applicable to analysis of natural fluid inclusions.     

A new approach of selecting real input ground motions for seismic design: The most unfavourable real seismic design ground motions

This paper presents a new way of selecting real input ground motions for seismic design and analysis of structures based on a comprehensive method for estimating the damage potential of ground motions, which takes into consideration of various ground motion parameters and structural seismic damage criteria in terms of strength, deformation, hysteretic energy and dual damage of Park & Ang damage index. The proposed comprehensive method fully involves the effects of the intensity, frequency content and duration of ground motions and the dynamic characteristics of structures. Then, the concept of the most unfavourable real seismic design ground motion is introduced. Based on the concept, the most unfavourable real seismic design ground motions for rock, stiff soil, medium soil and soft soil site conditions are selected in terms of three typical period ranges of structures. The selected real strong motion records are suitable for seismic analysis of important structures whose failure or collapse will be avoided at a higher level of confidence during the strong earthquake, as they can cause the greatest damage to structures and thereby result in the highest damage potential from an extended real ground motion database for a given site. In addition, this paper also presents the real input design ground motions with medium damage potential, which can be used for the seismic analysis of structures located at the area with low and moderate seismicity. The most unfavourable real seismic design ground motions are verified by analysing the seismic response of structures. It is concluded that the most unfavourable real seismic design ground motion approach can select the real ground motions that can result in the highest damage potential for a given structure and site condition, and the real ground motions can be mainly used for structures whose failure or collapse will be avoided at a higher level of confidence during the strong earthquake. Copyright © 2007 John Wiley & Sons, Ltd.     

Inter- and intra-annual trend analysis of water level and flow in the middle and lower reaches of the Ganjiang River,China

ABSTRACT

The Ganjiang River is the largest tributary of Poyang Lake in China, and its hydrological regime variation greatly affects the utilization of regional water resources and the ecological environment of the lake. In this study, a novel trend analysis method, the Moving Average over Shifting Horizon (MASH), was applied to investigate the inter- and intra-annual trends of flow and water level from 1976 to 2016 at the Xiajiang and the Waizhou hydrological stations in the Ganjiang River. The Significant Change Rate Method (SCRM) was proposed to determine the MASH averaging parameters. The trend analysis results show a statistically significant decrease in water level series throughout the year and the relationship of flow and water level have changed greatly at the Waizhou station. The sediment load reduction, large-scale sand mining and water level decrease of Poyang Lake are identified as the main causes for the water level decrease.     

干旱区膜下滴灌棉田表层土壤盐分日内动态特征

利用便携式土壤盐分计测量膜下滴灌棉田不同位置和不同深度的土壤电导率,来研究膜下滴灌条件下表层土壤电导率的日内变化特征.结果显示:干旱区一膜两管膜下滴灌条件下,窄行处土壤电导率值最低,膜间最大,宽行居中,表层土壤电导率值存在显著的日内变化特征,早晚土壤电导率值低,中午电导率值高,上午的上升速率大于下午的下降速率;窄行处相对变化幅度大,膜间绝对变化幅度大,土壤电导率值日内变化幅度随深度呈减小趋势;膜下滴灌情况下,表层土壤电导率值的变化受土壤温度和气温的影响显著,但二者相关性复杂,不同深度、不同时间段的相关性不同,膜间大多呈显著正相关,窄行大多呈非显著正相关,宽行处上午呈非显著正相关,但下午呈显著负相关.     

中国大陆地震地下流体异常特征研究

根据《中国震例》(1966~1999年)在中国大陆发生M_S≥5.0地震188次震例资料 , 将其中与流体异常相关的126个震例中的803项流体异常按异常出现、 异常转折至发震随时间变化和异常数量随震中距距离变化分布的时、 空特征进行统计、 分析。 研究结果表明, 无论是趋势性异常还是短临异常, 在不同的时间段中, 异常数量呈现出明显的阶段性分布特征; 异常数量空间分布会随震中距的变化而不同。 在统计研究中还表明, 两者在时、 空分布方面有较好的一致性。     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

利用奇异谱分析提高极移中长期预报精度

由于空间大地观测数据传输耗时及处理过程复杂, 导致极移测量值的获取存在时延, 无法满足对高精度的极移预报值有重大需求的应用领域. 针对极移复杂的时变特性, 提出一种基于奇异谱分析(singular spectrum analysis, SSA)的预报方法. 首先用SSA分离提取极移时序中的高频组分与低频组分; 其次建立最小二乘(least square, LS)外推与自回归(AutoreGressive, AR)模型对极移高频和低频组分进行组合预报. 结果表明, SSA方法能够准确地分离和提取极移低频和高频组分, 对低频和高频组分组合预报可以显著改善极移的中长期(30--365d)预报精度, 与国际地球自转和参考系服务局(International Earth Rotation and Reference Systems Service, IERS)提供的A公报中的极移预报值相比, SSA方法对极移PMX分量(本初子午线方向)和PMY分量(西90$^\circ$子午线方向)的中长期预报精度改进最高分别可达45.97%和62.44%. 研究结果验证了SSA方法对极移中长期预报改进的有效性.     

大型射电望远镜主反射面调整方法研究

随着毫米波天文学和空间通信的重要性日益提高, 对天线性能提出了越来越高的要求, 而天线性能往往受到其反射器表面精度的限制. 微波全息技术是一种快速有效的检测反射面天线表面轮廓的测量技术. 通过微波全息测量得到天线口径场, 计算天马65m射电望远镜反射面与理想抛物面的偏差. 天马65m射电望远镜的主反射面板是放射状的, 有14圈. 面板的每个角都固定在面板下方促动器的螺栓上进行上下移动, 且相邻面板交点处的拐角共用一个促动器. 采用平面拟合的方法可以计算各块面板拐角处的调整值, 但是同一个促动器会得到4个不同的调整量. 通过平面拟合, 同时以天线照明函数为权重的平差计算方法得到相邻面板拐角的一个平差值, 即天马65m射电望远镜1104个促动器的最佳调整值. 通过多次调整和新算法的应用, 天马65m射电望远镜反射面的面形精度逐渐提高到了0.24mm.     

高层大气密度的午夜极大值特征分析与建模

利用CHAMP (CHAllenging Minisatellite Payload)、GRACE-A (Gravity Recovery and Climate Experiment-A)、SWARM-C (The Earth''s Magnetic Field and Environment Explorers-C)等3颗极轨卫星的资料, 研究360—480km高层大气密度在低纬度区域的午夜极大值(Midnight Density Maximum, MDM)现象. MDM一般出现在23:00- 02:00 LT (Local Time)之间,峰值位置在低纬度15°以内,谷值位置在中纬度35°-45°附近,整体略偏向南半球,振幅约为平均密度的26%.随着高度增大以及太阳辐射水平的增强,MDM振幅呈减小趋势;冬至和夏至日附近的季节效应会减弱MDM振幅,在春秋分日的振幅最大.用3个主流大气模型DTM2000 (Drag Temperature Model 2000), NRLMSISE00 (US Naval Research Laboratory, Mass Spectrom-eter and Incoherent Scatter radar)和JB2008 (Jacchia- Bowman 2008 model)对MDM进行模拟,JB2008没有刻画出MDM现象;另两个模型低估了MDM效应,在360km和480km两个高度DTM2000模型的振幅仅为观测的46%和53%, NRLMSISE00模型仅为观测的33%和26%;模型没有准确刻画出MDM与高度、辐射水平和季节的关系.联合3颗卫星的资料,研究了-种基于地理纬度的6阶勒让德多项式,同时融合地方时和高度因素的经验函数,在振幅和相位上可以较好地刻画MDM特征,相关系数达到0.923,可为大气密度模型的修正提供借鉴,服务于低轨道航天器高精度轨道预报.