利用不完全Cholesky共轭梯度法求解点源三维地电场

点源三维地电场的求解是一大型数值计算问题．本文用有限差分方法求解,最后形成一个线性方程组Ax＝b,这里A是大型稀疏的带状对称矩阵．解大型稀疏方程组的完全Cholesky分解直接算法,一般要求巨大的机器内存来存储系数矩阵A,而且计算速度极慢．因此引入不完全Cholesky共轭梯度（ICCG）算法及按行索引的稀疏存储模式,充分利用A的稀疏性,使得计算速度大大提高,而内存要求则大大减少,因此ICCG算法是地电三维正演的强有力工具．     

THREE-DIMENSIONAL NUMERICAL MODELING OF COHESIVE SEDIMENT TRANSPORT BY TIDAL CURRENT IN PEARL RIVER ESTUARV

IINTRODUCTIONEstUariesareprominentcoastalfeatUres.Estuariesareofgreateconomicssignificancetomankind.Attheseareas,manyharborsandwaterchannelshavetobebuiltforeconomicpurposes.ThedesignandconstrUctionofcoastalstrUctUresinestUariesrequireknowledgeofhydrodynamicsaswellassedimenttransportinsuchregions.ThenatUreofestuariesiscontrolledbyvariouscoastalhydrodynamicprocesses.Undertheactionofhydrodynamics,sedimentdepositionsorerosionswilloccurinestuariesornearcoastalstrUCtures.Tomaintainnavigati…     

基于密集线性台阵的背景噪声成像在明光市城市活断层调查中的应用

近年来,短周期密集台阵被动源探测技术逐渐成为地壳精细结构探测研究中的重要手段.为了探讨其应用于城市活断层调查在勘探尺度的有效性,在横跨安徽省明光市郯庐断裂带东界两条主干断裂位置,沿同期的主动源S波反射剖面,安徽省地震局和中国科学技术大学联合布设了一条由133个短周期地震仪组成的被动源密集测线.我们通过该被动源测线获得了约28天的连续波形记录,并采用一种新的线性台阵背景噪声成像方法——拓距相移法提取了10 Hz~2.0 s宽频的基阶Rayleigh波相速度频散曲线.利用得到的宽频带频散曲线,我们反演了近地表 1.3 km以浅的高分辨率S波速度结构,并与基本重合的主动源S波反射剖面以及地质资料等进行对比,验证了反演结果的可靠性,同时进一步分析了郯庐断裂带东界两条主干断裂的断层特征.结果表明：（1）被动源与主动源反演结果在近地表100 m以浅高度一致,且被动源的反演深度更深;（2）池河—太湖断裂是隐伏于深部的走滑断裂和浅部伸展断裂的综合反映,浅部为倾向北西的正断层,向上穿透第四系地层,可能为活动断裂;（3）拓距相移法能有效拓宽频散曲线的提取频带,保证浅部结构具有高分辨率,同时大大加强对深部结构的约束.该方法能够获得可靠的地下结构反演结果,可以作为城市活断层调查的一种新的有效手段.

     

A local model for analysis of pump and treat systems with vertical barrier walls

We present a mathematical model of local, steady groundwater flow near a vertical barrier wall. Flow features represented in the model include an impermeable arc-shaped barrier wall and multiple wells; distant boundary conditions are not included explicitly, but their effects on the local flow field are modelled by specifying a uniform flow at infinity and a constant areal recharge within a local domain. We develop an explicit closed-form solution to the boundary-value problem using the analytic element method. The solution is an extension of a harmonic solution presented by Anderson and Mesa [Anderson EI, Mesa E. The effects of vertical barrier walls on the hydraulic control of contaminated groundwater. Adv Water Resourc 2006;29(1):89–98] which does not include the effects of recharge. We demonstrate that the general solution with recharge consists of the harmonic solution superposed on a special case of the harmonic solution along with two elementary one-dimensional flow solutions. The results are used to investigate the effects of areal recharge on the capture zone envelopes of the pumping wells and on the reduction in discharge that can be achieved by including a barrier wall in a pump and treat design. We find that the benefits of including an open barrier wall in a design, measured as a reduction in the pumping rate required to contain a plume, increase for higher recharge rates. Dimensionless plots of capture zone envelopes are presented for a practical well and barrier wall configuration.     

利用GRACE卫星Level-1B数据反演陆地水储量变化的方法研究

卫星重力探测技术为监测陆地水储量变化提供了新的技术手段。对利用GRACE卫星Level-1B数据反演地球陆地水储量变化的重力位差法和Mascon方法的数学模型作了详细推导分析,总结两种方法的特点和解算处理步骤。推导过程表明:重力位差法和Mascon方法在反演时只采用卫星飞临研究区域上空时的观测数据,能够提高反演结果的空间分辨率,比传统的球谐系数法更具优势;Mascon方法在解算时还引入了时空约束方程,进一步提高了反演结果的时空分辨率。     

Seismic,petrologic, and geodetic analyses of the 1999 dome-forming eruption of Guagua Pichincha volcano,Ecuador

Guagua Pichincha, located 14 km west of Quito, Ecuador, is a stratovolcano bisected by a horseshoe-shaped caldera. In 1999, after some months of phreatic activity, Guagua Pichincha entered into an eruptive period characterized by the extrusion of several dacitic domes, vulcanian eruptions, and pyroclastic flows. We estimated the three-dimensional (3-D) P-wave velocity structure beneath Guagua Pichincha using a tomographic inversion method based on finite-difference calculations of first-arrival times. Hypocenters of volcano-tectonic (VT) earthquakes and long-period (LP) events were relocated using the 3-D P-wave velocity model. A low-velocity anomaly exists beneath the caldera and may represent an active volcanic conduit. Petrologic analysis of eruptive products indicates a magma storage region beneath the caldera, having a vertical extent of 7–8 km with the upper boundary at about sea level. This zone coincides with the source region of deeper VT earthquakes, indicating that a primary magma body exists in this region. LP swarms occurred in a cyclic pattern synchronous with ground deformation during magma extrusions. The correlation between seismicity and ground deformation suggests that both respond to pressure changes caused by the cyclic eruptive behavior of lava domes.     

鲁豫交界地区尾波Q_c值分析研究

该文对濮阳台记录的21次地震进行了地震波衰减研究,讨论了地球介质品质因子随频率变化的关系;利用不同地震采用不同流逝时间的尾波品质因子测定方法,有效去除了由于不同地震或同一地震的地震散射波在纵向介质中穿透深度的不同而导致对品质因子测量结果的影响,从而有利于扑捉到品质因子的异常变化特征与地震孕育过程的关系。利用改进后的模型,对鲁豫交界地区尾波值进行了计算,得到了该地区品质因子随频率的变化关系和区域介质衰减的时间变化特征。     

季节性Kendall检验及其在三门峡水库水质趋势分析中的应用

选择1997—2002年三门峡库区的龙门、潼关和三门峡3个断面的水质监测资料,运用季节性Kendall检验方法分析三门峡水库水质趋势,结果显示：龙门断面,CODMn、NO2-N及Hg的变化趋势不明显,NH3-N呈显著上升趋势,As、Cd和Cu呈显著下降趋势,而Pb则呈高度显著下降趋势;潼关断面,NH3-N、NO2-N和Hg的变化趋势不明显,CODMn呈显著上升趋势,As、Cd和Pb呈显著下降趋势,而Cu则呈高度显著下降趋势;三门峡断面,CODMn、NO2-N、As及Hg的变化趋势不明显,NH3-N污染呈显著上升趋势,Cd呈显著下降趋势,而Cu和Pb则呈高度显著下降趋势。     

Errors and correction of precipitation measurements in China

In order to discover the range of various errors in Chinese precipitation measurements and seek a correction method, 30 precipitation evaluation stations were set up countrywide before 1993. All the stations are reference stations in China. To seek a correction method for wind-induced error, a precipitation correction instrument called the ＂horizontal precipitation gauge＂ was devised beforehand. Field intercomparison observations regarding 29,000 precipitation events have been conducted using one pit gauge, two elevated operational gauges and one horizontal gauge at the above 30 stations. The range of precipitation measurement errors in China is obtained by analysis of intercomparison measurement results. The distribution of random errors and systematic errors in precipitation measurements are studied in this paper. A correction method, especially for wind-induced errors, is developed. The results prove that a correlation of power function exists between the precipitation amount caught by the horizontal gauge and the absolute difference of observations implemented by the operational gauge and pit gauge. The correlation coefficient is 0.99. For operational observations, precipitation correction can be carried out only by parallel observation with a horizontal precipitation gauge. The precipitation accuracy after correction approaches that of the pit gauge. The correction method developed is simple and feasible.     

A CIP-based numerical simulation of free surface flow related to freak waves

An enhanced numerical model for simulating two-dimensional incompressible viscous flow with distorted free surface is reported. The numerical simulation is carried out through the CIP （Constrained Interpolation Profile）-based method, which is described in the paper. A more accurate interface capturing scheme, the VOF/WLIC scheme （VOF：Volume-of-Fluid;WLIC：weighed line interface calculation）, is adopted as the interface capturing method. To assess the developed algorithm and its versatility, a selection of test problems are examined, i.e. the square wave propagation, the Zalesak’s rigid body rotation, dam breaking problem with and without obstacles, wave sloshing in an excited wave tank and interaction between extreme waves and a floating body. Excellent agreements are obtained when numerical results are compared with available analytical, experimental, and other numerical results. These examples demonstrate that the use of the VOF/WLIC scheme in the free surface capturing makes better results and also the proposed CIP-based model is capable of predicting the freak wave-related phenomena.